1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
/*@@
@file GRHydro_Minima.F90
@date Mon Feb 25 11:43:36 2002
@author
@desc
Sets up the scalars used for the atmosphere, before initial data.
@enddesc
@@*/
#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"
/*@@
@routine GRHydro_Minima_Setup
@date Mon Feb 25 11:25:27 2002
@author Ian Hawke
@desc
Before initial data, set up the scalar GRHydro_rho_min used for the atmosphere.
This is computed only from parameters.
@enddesc
@calls
@calledby
@history
Modified 30 Aug 2006 by Luca Baiotti
@endhistory
@@*/
subroutine GRHydro_Rho_Minima_Setup(CCTK_ARGUMENTS)
implicit none
DECLARE_CCTK_ARGUMENTS
DECLARE_CCTK_PARAMETERS
if (initial_rho_abs_min > 0.0) then
GRHydro_rho_min = initial_rho_abs_min
else if (initial_rho_rel_min > 0.0) then
GRHydro_rho_min = GRHydro_rho_central * initial_rho_rel_min
else if (rho_abs_min > 0.d0) then
GRHydro_rho_min = rho_abs_min
else
GRHydro_rho_min = GRHydro_rho_central * rho_rel_min
end if
if (initial_atmosphere_factor > 0.0) GRHydro_rho_min = GRHydro_rho_min * initial_atmosphere_factor
return
end subroutine GRHydro_Rho_Minima_Setup
/*@@
@routine GRHydro_Check_Rho_Minimum
@date Mon Jul 7 16:35:45 2008
@author Luca Baiotti
@desc
Check whether at some point rho < GRHydro_rho_min and print a warning in case.
@enddesc
@calls
@calledby
@history
@endhistory
@@*/
subroutine GRHydro_Check_Rho_Minimum(CCTK_ARGUMENTS)
implicit none
DECLARE_CCTK_ARGUMENTS
DECLARE_CCTK_PARAMETERS
CCTK_INT i,j,k
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
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
call CCTK_WARN(2,warnline)
write(warnline,'(a25,g15.6)') 'rho: ',rho(i,j,k)
call CCTK_WARN(2,warnline)
write(warnline,'(a25,4g15.6)') 'coordinates: x,y,z,r:',x(i,j,k),y(i,j,k),z(i,j,k),r(i,j,k)
call CCTK_WARN(2,warnline)
end if
end do
end do
end do
return
end subroutine GRHydro_Check_Rho_Minimum
/*@@
@routine GRHydro__Change_Rho_Minimum_At_Recovery
@date Thu Aug 14 17:11:32 2008
@author Luca Baiotti
@desc
Change, via a parameter, the value of GRHydro_rho_min at recovery.
@enddesc
@calls
@calledby
@history
@endhistory
@@*/
subroutine GRHydro_Change_Rho_Minimum_At_Recovery(CCTK_ARGUMENTS)
implicit none
DECLARE_CCTK_ARGUMENTS
DECLARE_CCTK_PARAMETERS
GRHydro_rho_min = rho_abs_min_after_recovery
return
end subroutine GRHydro_Change_Rho_Minimum_At_Recovery
|
