diff options
Diffstat (limited to 'src/GRHydro_HLLE.F90')
-rw-r--r-- | src/GRHydro_HLLE.F90 | 20 |
1 files changed, 10 insertions, 10 deletions
diff --git a/src/GRHydro_HLLE.F90 b/src/GRHydro_HLLE.F90 index bfe3cd0..deb37e3 100644 --- a/src/GRHydro_HLLE.F90 +++ b/src/GRHydro_HLLE.F90 @@ -216,7 +216,7 @@ subroutine GRHydro_HLLE(CCTK_ARGUMENTS) call SpaceMask_GetStateBits(trivial, "Hydro_RiemannProblemZ", & &"trivial") else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if if(evolve_temper.eq.1.and.reconstruct_temper.eq.1) then @@ -272,7 +272,7 @@ subroutine GRHydro_HLLE(CCTK_ARGUMENTS) avg_beta = 0.5d0 * (beta3(i+xoffset,j+yoffset,k+zoffset) + & beta3(i,j,k)) else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if avg_alp = 0.5 * (alp(i,j,k) + alp(i+xoffset,j+yoffset,k+zoffset)) @@ -311,7 +311,7 @@ subroutine GRHydro_HLLE(CCTK_ARGUMENTS) velzplus(i,j,k),pressplus(i,j,k),& avg_det,avg_alp,avg_beta) else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if else !!! The end of this branch is right at the bottom of the routine @@ -327,7 +327,7 @@ subroutine GRHydro_HLLE(CCTK_ARGUMENTS) else if (flux_direction == 3) then usendh = uzzh else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if !!$ Calculate the jumps in the conserved variables @@ -498,7 +498,7 @@ subroutine GRHydro_HLLE(CCTK_ARGUMENTS) pressminus(i+xoffset,j+yoffset,k+zoffset),& avg_det,avg_alp,avg_beta) else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if !!$ Compute H viscosity if requested @@ -536,7 +536,7 @@ subroutine GRHydro_HLLE(CCTK_ARGUMENTS) etaY(i,j-1,k, vup, eos_c),& etaY(i,j-1,k+1, vup, eos_c)) else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if ! modify eigenvalues of Roe's matrix by computed H viscosity @@ -678,7 +678,7 @@ subroutine GRHydro_HLLE_Tracer(CCTK_ARGUMENTS) call SpaceMask_GetStateBits(trivial, "Hydro_RiemannProblemZ", & &"trivial") else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil @@ -714,7 +714,7 @@ subroutine GRHydro_HLLE_Tracer(CCTK_ARGUMENTS) avg_beta = 0.5d0 * (beta3(i+xoffset,j+yoffset,k+zoffset) + & beta3(i,j,k)) else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if avg_alp = 0.5 * (alp(i,j,k) + alp(i+xoffset,j+yoffset,k+zoffset)) @@ -740,7 +740,7 @@ subroutine GRHydro_HLLE_Tracer(CCTK_ARGUMENTS) else if (flux_direction == 3) then usendh = uzzh else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if !!$ Calculate the jumps in the conserved variables @@ -813,7 +813,7 @@ subroutine GRHydro_HLLE_Tracer(CCTK_ARGUMENTS) fminus(:) = (velzminus(i+xoffset,j+yoffset,k+zoffset) - avg_beta / avg_alp) * & cons_tracerminus(i+xoffset,j+yoffset,k+zoffset,:) else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if !!$ Find minimum and maximum wavespeeds |