diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:27 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:27 +0000 |
| commit | c141ed811d4a5b84ecc077999d1fec0a3c598b95 (patch) | |
| tree | 92f2d30c9ae9fbed6e307a802b841f75eadab0b9 /src/GRHydro_Minima.F90 | |
| parent | d5c4a808aeefa7fe95cc25cb7648fd7938561b8d (diff) | |
GRHydro: Introduce variable atmopshere level / tolerance as function of radius. Off by default.
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@453 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Minima.F90')
| -rw-r--r-- | src/GRHydro_Minima.F90 | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/src/GRHydro_Minima.F90 b/src/GRHydro_Minima.F90 index 5f870cd..905a65c 100644 --- a/src/GRHydro_Minima.F90 +++ b/src/GRHydro_Minima.F90 @@ -11,6 +11,8 @@ #include "cctk_Parameters.h" #include "cctk_Arguments.h" +#include "GRHydro_Macros.h" + /*@@ @routine GRHydro_Minima_Setup @date Mon Feb 25 11:25:27 2002 @@ -74,17 +76,19 @@ subroutine GRHydro_Check_Rho_Minimum(CCTK_ARGUMENTS) DECLARE_CCTK_PARAMETERS CCTK_INT i,j,k + CCTK_REAL dummy1 character(len=100) warnline do i=1,cctk_lsh(1) do j=1,cctk_lsh(2) do k=1,cctk_lsh(3) - if (rho(i,j,k) < GRHydro_rho_min) then + SET_ATMO_MIN(dummy1, GRHydro_rho_min, r(i,j,k)) + if (rho(i,j,k) < dummy1) then call CCTK_WARN(2,"rho<GRHydro_rho_min!!!") write(warnline,'(a28,i2)') 'on carpet reflevel: ',GRHydro_reflevel call CCTK_WARN(2,warnline) - write(warnline,'(a25,g15.6)') 'GRHydro_rho_min: ',GRHydro_rho_min + write(warnline,'(a25,g15.6)') 'GRHydro_rho_min: ', dummy1 call CCTK_WARN(2,warnline) write(warnline,'(a25,g15.6)') 'rho: ',rho(i,j,k) call CCTK_WARN(2,warnline) |