diff options
| -rw-r--r-- | src/GRHydro_Prim2Con.F90 | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/GRHydro_Prim2Con.F90 b/src/GRHydro_Prim2Con.F90 index 0533e27..24d70ac 100644 --- a/src/GRHydro_Prim2Con.F90 +++ b/src/GRHydro_Prim2Con.F90 @@ -46,7 +46,7 @@ subroutine primitive2conservative(CCTK_ARGUMENTS) integer :: i, j, k CCTK_REAL :: gxxl,gxyl,gxzl,gyyl,gyzl,gzzl,avg_detl,& gxxr,gxyr,gxzr,gyyr,gyzr,gzzr,avg_detr - CCTK_REAL :: xtemp + CCTK_REAL :: xtemp(1) if(evolve_temper.ne.1) then !$OMP PARALLEL DO PRIVATE(i, j, avg_detl, avg_detr,& @@ -121,7 +121,7 @@ subroutine primitive2conservative(CCTK_ARGUMENTS) ! variable 'temperature' in prim2con at the interfaces ! We will instead use an average temperature as an initial ! guess. - xtemp = 0.5d0*(temperature(i,j,k) + & + xtemp(1) = 0.5d0*(temperature(i,j,k) + & temperature(i-xoffset,j-yoffset,k-zoffset)) call prim2con_hot(GRHydro_eos_handle, GRHydro_reflevel,& i,j,k,x(i,j,k),y(i,j,k),z(i,j,k),gxxl,gxyl,gxzl,gyyl,& @@ -138,7 +138,7 @@ subroutine primitive2conservative(CCTK_ARGUMENTS) ! variable 'temperature' in prim2con at the interfaces ! We will instead use an average temperature as an initial ! guess. - xtemp = 0.5d0*(temperature(i,j,k) + & + xtemp(1) = 0.5d0*(temperature(i,j,k) + & temperature(i+xoffset,j+yoffset,k+zoffset)) call prim2con_hot(GRHydro_eos_handle, GRHydro_reflevel, & i,j,k,x(i,j,k),y(i,j,k),z(i,j,k),gxxr,gxyr,gxzr,& |