User: knarf Date: 2011/02/10 09:23 AM

Modified: /trunk/src/ GRHydro_Prim2Con.F90

Log: fix declaration of xtemp, which has to be an array of 1 according to the function declaration of prim2con_hot

File Changes:

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

File [modified]: GRHydro_Prim2Con.F90 Delta lines: +3 -3 =================================================================== --- trunk/src/GRHydro_Prim2Con.F90 2011-02-09 19:23:22 UTC (rev 216) +++ trunk/src/GRHydro_Prim2Con.F90 2011-02-10 15:23:57 UTC (rev 217) @@ -46,7 +46,7 @@ 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 @@ ! 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 @@ ! 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,&