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@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