diff --git a/src/GRHD/cons2prim.jl b/src/GRHD/cons2prim.jl
index ad547ea134a2aa81af54a4350a624982dd9e38c2..11b6bfa6378bbcb958229af517811b978479aa7a 100644
--- a/src/GRHD/cons2prim.jl
+++ b/src/GRHD/cons2prim.jl
@@ -191,14 +191,15 @@ end
     k = r / (1.0+q)
     if !(0.0 ≤ k ≤ 1.0) # admissibility constraint from arXiv:1306.4953, Appendix
       c2p_init_admissible = 1.0
-      # rescale Si such that r = 1+q
-      rescale = 0.99*(1+q)/r
-      Sx *= rescale
-      Sz *= rescale
-      rx = Sx / D
-      rz = Sz / D
-      r2 = γuuxx*rx*rx+γuuzz*rz*rz
-      r  = sqrt(r2)
+      @goto impose_atmosphere
+      # # rescale Si such that r = 1+q
+      # rescale = 0.99*(1+q)/r
+      # Sx *= rescale
+      # Sz *= rescale
+      # rx = Sx / D
+      # rz = Sz / D
+      # r2 = γuuxx*rx*rx+γuuzz*rz*rz
+      # r  = sqrt(r2)
     end
     z0 = r/h0
     v0 = z0 / sqrt(1.0 + z0^2)
diff --git a/test/IntegrationTests/refs/grhd_bondi_infall_avmda/output1d.h5 b/test/IntegrationTests/refs/grhd_bondi_infall_avmda/output1d.h5
index 220135c6a0abbf3b14330cca8d8708b86cf1097e..b81e3c59a825686dabb1aa484efd470af2a683b8 100644
Binary files a/test/IntegrationTests/refs/grhd_bondi_infall_avmda/output1d.h5 and b/test/IntegrationTests/refs/grhd_bondi_infall_avmda/output1d.h5 differ