[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 449)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/01/14 08:23 AM

Modified:
 /trunk/src/
  GRHydro_Con2PrimM.F90

Log:
 GRHydro: call pointwise con2primM with polytrope eos handle deep inside AH upon inversion failure of gamma law eos.
 
 From: Bruno Coutinho Mundim <bcmsma at astro.rit.edu>

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_Con2PrimM.F90
Delta lines: +28 -0
===================================================================
--- trunk/src/GRHydro_Con2PrimM.F90	2013-01-14 14:23:16 UTC (rev 448)
+++ trunk/src/GRHydro_Con2PrimM.F90	2013-01-14 14:23:19 UTC (rev 449)
@@ -344,6 +344,34 @@
                    uxx,uxy,uxz,uyy,uyz,uzz,det, &
                    epsnegative,GRHydro_C2P_failed(i,j,k))
 
+              if(sdet.ge.sqrtdet_thr) then
+                if(GRHydro_C2P_failed(i,j,k).ne.0) then
+                  GRHydro_C2P_failed(i,j,k) = 0
+
+                  rho_tmp = rho(i,j,k)
+                  press_tmp = press(i,j,k)
+                  eps_tmp = eps(i,j,k)
+                  velx_tmp = vup(i,j,k,1)
+                  vely_tmp = vup(i,j,k,2)
+                  velz_tmp = vup(i,j,k,3)
+                  w_lorentz_tmp = w_lorentz(i,j,k)
+                  Bvecx_tmp = Bprim(i,j,k,1)
+                  Bvecy_tmp = Bprim(i,j,k,2)
+                  Bvecz_tmp = Bprim(i,j,k,3)
+
+                  call GRHydro_Con2PrimM_pt(GRHydro_polytrope_handle, keytemp, &
+                       GRHydro_eos_rf_prec, local_gam(1), dens(i,j,k), &
+                       scon(i,j,k,1),scon(i,j,k,2),scon(i,j,k,3), tau(i,j,k), &
+                       Bcons(i,j,k,1),Bcons(i,j,k,2),Bcons(i,j,k,3),xye(1), &
+                       xtemp(1),rho_tmp,velx_tmp,vely_tmp,velz_tmp,&
+                       eps_tmp,press_tmp,Bvecx_tmp,Bvecy_tmp,Bvecz_tmp,b2,&
+                       w_lorentz_tmp,g11(i,j,k),g12(i,j,k),g13(i,j,k),&
+                       g22(i,j,k),g23(i,j,k),g33(i,j,k), &
+                       uxx,uxy,uxz,uyy,uyz,uzz,det, &
+                       epsnegative,GRHydro_C2P_failed(i,j,k))
+                end if
+              end if
+
            else    ! if(evolve_temper.eq.0) then
 
               rho_tmp = rho(i,j,k)