GRHD: impose atmosphere when momentum admissibility condition is violated...

GRHD: impose atmosphere when momentum admissibility condition is violated (!224)

The condition is from arXiv:1306.4953 where they suggest to rescale the momemtum to fix this.
I think that rescaling is as ambiguous as resetting to atmosphere.
I decided to to the latter now as this is how I wrote it up in the thesis.

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)