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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/07/06 01:11 PM

Modified:
 /trunk/src/
  GRHydro_Prim2Con.F90, GRHydro_Prim2ConM.F90

Log:
 GRHydro: pass correct temperature o prim2con call on coarse levels
 
 Fix a bug in prim2con_hot that was causing temp to go negative near the
 shock in 3D CCSN simulations with WENO reconstruction.
 
 From: abdik <abdik at tapir.caltech.edu>

File Changes:

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

File [modified]: GRHydro_Prim2Con.F90
Delta lines: +1 -0
===================================================================
--- trunk/src/GRHydro_Prim2Con.F90	2013-07-06 18:11:00 UTC (rev 553)
+++ trunk/src/GRHydro_Prim2Con.F90	2013-07-06 18:11:04 UTC (rev 554)
@@ -378,6 +378,7 @@
            ! away), we use the average temperature between cells and call
            ! the EOS with keytemp=1
            lkeytemp=1
+           temp = temp0
            call EOS_Omni_press(handle,lkeytemp,GRHydro_eos_rf_prec,n,&
                 drho,deps,temp,ye,dpress,keyerr,anyerr)  
            lkeytemp=0

File [modified]: GRHydro_Prim2ConM.F90
Delta lines: +1 -0
===================================================================
--- trunk/src/GRHydro_Prim2ConM.F90	2013-07-06 18:11:00 UTC (rev 553)
+++ trunk/src/GRHydro_Prim2ConM.F90	2013-07-06 18:11:04 UTC (rev 554)
@@ -294,6 +294,7 @@
         ! away), we use the average temperature between cells and call
         ! the EOS with keytemp=1
         keytemp=1
+        temp = temp0
         call EOS_Omni_press(handle,keytemp,GRHydro_eos_rf_prec,n,&
              drho,deps,temp,ye,dpress,keyerr,anyerr)  
         keytemp=0