diff options
Diffstat (limited to 'src/GRHydro_Marquina.F90')
-rw-r--r-- | src/GRHydro_Marquina.F90 | 20 |
1 files changed, 8 insertions, 12 deletions
diff --git a/src/GRHydro_Marquina.F90 b/src/GRHydro_Marquina.F90 index 58dd742..ac54b0a 100644 --- a/src/GRHydro_Marquina.F90 +++ b/src/GRHydro_Marquina.F90 @@ -111,19 +111,15 @@ subroutine GRHydro_Marquina(CCTK_ARGUMENTS) !!$ Set metric terms at interface - if (shift_state .ne. 0) then - if (flux_direction == 1) then - avg_beta = 0.5d0 * (betax(i+xoffset,j+yoffset,k+zoffset) + & - betax(i,j,k)) - else if (flux_direction == 2) then - avg_beta = 0.5d0 * (betay(i+xoffset,j+yoffset,k+zoffset) + & - betay(i,j,k)) - else - avg_beta = 0.5d0 * (betaz(i+xoffset,j+yoffset,k+zoffset) + & - betaz(i,j,k)) - end if + if (flux_direction == 1) then + avg_beta = 0.5d0 * (betax(i+xoffset,j+yoffset,k+zoffset) + & + betax(i,j,k)) + else if (flux_direction == 2) then + avg_beta = 0.5d0 * (betay(i+xoffset,j+yoffset,k+zoffset) + & + betay(i,j,k)) else - avg_beta = 0.d0 + avg_beta = 0.5d0 * (betaz(i+xoffset,j+yoffset,k+zoffset) + & + betaz(i,j,k)) end if avg_alp = 0.5 * (alp(i,j,k) + alp(i+xoffset,j+yoffset,k+zoffset)) |