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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/05/29 11:35 AM

Modified:
 /trunk/src/
  GRHydro_ENOReconstruct_drv.F90, GRHydro_PPMMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv.F90, GRHydro_TVDReconstruct_drv.F90, GRHydro_WENOReconstruct_drv.F90

Log:
 GRHydro: initialy set y_e_plus/minus to cell center value in reconstruct
 
 This ensures that all EOS calls have some sensible y_e to work with.
 
 From: Christian Ott <cott at bethe.tapir.caltech.edu>

File Changes:

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

File [modified]: GRHydro_ENOReconstruct_drv.F90
Delta lines: +5 -2
===================================================================
--- trunk/src/GRHydro_ENOReconstruct_drv.F90	2013-05-29 16:35:08 UTC (rev 530)
+++ trunk/src/GRHydro_ENOReconstruct_drv.F90	2013-05-29 16:35:12 UTC (rev 531)
@@ -157,8 +157,11 @@
            endif
 
            if (evolve_Y_e .ne. 0) then
-              Y_e_plus(i,j,k) = 0.0d0
-              Y_e_minus(i,j,k) = 0.0d0
+              ! set this to the cell center values
+              ! to make sure we have good Y_e even in
+              ! the boundary region (for full GF EOS calls)       
+              Y_e_plus(i,j,k) = Y_e(i,j,k)
+              Y_e_minus(i,j,k) = Y_e(i,j,k)
            endif
 
            if (evolve_temper .ne. 0) then

File [modified]: GRHydro_PPMMReconstruct_drv.F90
Delta lines: +8 -2
===================================================================
--- trunk/src/GRHydro_PPMMReconstruct_drv.F90	2013-05-29 16:35:08 UTC (rev 530)
+++ trunk/src/GRHydro_PPMMReconstruct_drv.F90	2013-05-29 16:35:12 UTC (rev 531)
@@ -175,10 +175,16 @@
            endif
 
            if (evolve_Y_e .ne. 0) then
-              Y_e_plus(i,j,k) = 0.0d0
-              Y_e_minus(i,j,k) = 0.0d0
+              Y_e_plus(i,j,k) = Y_e(i,j,k)
+              Y_e_minus(i,j,k) = Y_e(i,j,k)
            endif
 
+           if(evolve_temper .ne. 0) then
+              tempplus(i,j,k) = temperature(i,j,k)
+              tempminus(i,j,k) = temperature(i,j,k)
+           endif
+
+
         enddo
      enddo
   enddo

File [modified]: GRHydro_PPMReconstruct_drv.F90
Delta lines: +5 -2
===================================================================
--- trunk/src/GRHydro_PPMReconstruct_drv.F90	2013-05-29 16:35:08 UTC (rev 530)
+++ trunk/src/GRHydro_PPMReconstruct_drv.F90	2013-05-29 16:35:12 UTC (rev 531)
@@ -142,8 +142,11 @@
            endif
 
            if (evolve_Y_e .ne. 0) then
-              Y_e_plus(i,j,k) = 0.0d0
-              Y_e_minus(i,j,k) = 0.0d0
+              ! set this to the cell center values
+              ! to make sure we have good Y_e even in
+              ! the boundary region (for full GF EOS calls)
+              Y_e_plus(i,j,k) = Y_e(i,j,k)
+              Y_e_minus(i,j,k) = Y_e(i,j,k)
            endif
 
            if (evolve_temper .ne. 0) then

File [modified]: GRHydro_TVDReconstruct_drv.F90
Delta lines: +5 -2
===================================================================
--- trunk/src/GRHydro_TVDReconstruct_drv.F90	2013-05-29 16:35:08 UTC (rev 530)
+++ trunk/src/GRHydro_TVDReconstruct_drv.F90	2013-05-29 16:35:12 UTC (rev 531)
@@ -189,8 +189,11 @@
            endif
 
            if (evolve_Y_e .ne. 0) then
-              Y_e_plus(i,j,k) = 0.0d0
-              Y_e_minus(i,j,k) = 0.0d0
+              ! set this to the cell center values
+              ! to make sure we have good Y_e even in
+              ! the boundary region (for full GF EOS calls)
+              Y_e_plus(i,j,k) = Y_e(i,j,k)
+              Y_e_minus(i,j,k) = Y_e(i,j,k)
            endif
            
            if(evolve_temper .ne. 0) then

File [modified]: GRHydro_WENOReconstruct_drv.F90
Delta lines: +5 -2
===================================================================
--- trunk/src/GRHydro_WENOReconstruct_drv.F90	2013-05-29 16:35:08 UTC (rev 530)
+++ trunk/src/GRHydro_WENOReconstruct_drv.F90	2013-05-29 16:35:12 UTC (rev 531)
@@ -234,8 +234,11 @@
            endif
 
            if (evolve_Y_e .ne. 0) then
-              Y_e_plus(i,j,k) = 0.0d0
-              Y_e_minus(i,j,k) = 0.0d0
+              ! set this to the cell center values
+              ! to make sure we have good Y_e even in 
+              ! the boundary region (for full GF EOS calls)
+              Y_e_plus(i,j,k) = Y_e(i,j,k)
+              Y_e_minus(i,j,k) = Y_e(i,j,k)
            endif
 
            if (evolve_temper .ne. 0) then