From 292d02d83d2070dfb28939ac3b2afd3818c3ed40 Mon Sep 17 00:00:00 2001 From: knarf Date: Thu, 10 Feb 2011 15:23:57 +0000 Subject: fix declaration of xtemp, which has to be an array of 1 according to the function declaration of prim2con_hot git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@217 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Prim2Con.F90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'src/GRHydro_Prim2Con.F90') 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,& -- cgit v1.2.3