-
Notifications
You must be signed in to change notification settings - Fork 4
/
Copy pathPARAM.in.test.restart.OHPT
252 lines (202 loc) · 3.99 KB
/
PARAM.in.test.restart.OHPT
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
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
#COMPONENTMAP
OH 0 -1 1 CompMap runs on all PE-s
PT 1 2 1 CompMap runs on 1 PE
#DESCRIPTION
Param/PARAM.in.test.restart.OHPT: Restart single-ion OH and PT StringDescription
#LOOKUPTABLE
PT StringCompTable
ChargeExchange NameTable
load NameCommand
OH/Param/ExchangeRate.dat NameFile
ascii TypeFile
#TIMEACCURATE
F IsTimeAccurate
#INCLUDE
RESTART.in
#BEGIN_COMP PT ---------------------------------------------------------------
#RESTART
T isRestart
#END_COMP PT -----------------------------------------------------------------
#BEGIN_COMP OH ---------------------------------------------------------------
#INCLUDE
OH/restartIN/restart.H
#GRAVITY
F UseGravity (rest of parameters read if true)
0 iDirGravity(0 - central, 1 - X, 2 - Y, 3 - Z direction)
MINIMUMDENSITY
1e-7 RhoMin
MINIMUMPRESSURE
1e-20 pMin
#BODY
T UseBody
30.0 rBody (AU)
0.01 BodyNDim [/cc] number density for ions
1.09E5 BodyTDim [K] temperature for ions
#PLOTFILENAME
year NameMaxTimeUnit
#GRIDBLOCKALL
3500 MaxBlock
#GRIDLEVEL
2 nLevel
all StringShape
#GRIDLEVEL
5 nLevel (usually 11)
initial StringShape
#GRIDLEVEL
11 nLevel
shell0 StringShape
30.0 RadiusInner
30.2 Radius
#GRIDRESOLUTION
35.937 Resolution
box StringShape
0.0 xMinBox
-600. yMinBox
-600. zMinBox
-600. xMaxBox
600. yMaxBox
600. zMaxBox
#GRIDRESOLUTION
17.969 Resolution
box StringShape
-153.0 xMinBox
-340.0 yMinBox
-340.0 zMinBox
-600.0 xMaxBox
340.0 yMaxBox
340.0 zMaxBox
#GRIDRESOLUTION
8.984 Resolution
box StringShape
-109.0 xMinBox
-200.0 yMinBox
-200.0 zMinBox
-400.0 xMaxBox
200.0 yMaxBox
200.0 zMaxBox
#GRIDRESOLUTION
4.492 Resolution
box StringShape
-59.64 xMinBox
-136.0 yMinBox
-136.0 zMinBox
-272.0 xMaxBox
136.0 yMaxBox
136.0 zMaxBox
#GRIDRESOLUTION
3.0 Resolution
box StringShape
-71.31 xMinBox
-85.0 yMinBox
-85.0 zMinBox
-170.0 xMaxBox
85.0 yMaxBox
85.0 zMaxBox
#GRIDRESOLUTION
1.5 Resolution
box StringShape
-20.00 xMinBox
-42.0 yMinBox
-42.0 zMinBox
-68.0 xMaxBox
42.0 yMaxBox
42.0 zMaxBox
#GRIDRESOLUTION
0.7 Resolution
box StringShape
-5.0 xMinBox
-40.0 yMinBox
-40.0 zMinBox
-40.0 xMaxBox
40.0 yMaxBox
40.0 zMaxBox
#INNERBOUNDARY
user TypeBcBody
#OUTERBOUNDARY
user TypeBc1
outflow TypeBc2
outflow TypeBc3
outflow TypeBc4
outflow TypeBc5
outflow TypeBc6
NEUTRALFLUID
F DoConserveNeutrals
Linde TypeFluxNeutral
#SCHEME
1 nOrder
Linde TypeFlux
#TIMESTEPPING
1 nStage (1 or 2)
0.80 CflExpl
#NONCONSERVATIVE
T UseNonConservative
#CONSERVATIVECRITERIA
0 nConservCrit
#SAVERESTART
T DoSaveRestart
1000 DnSaveRestart
-1.0 DtSaveRestart
#SAVELOGFILE
T DoSaveLogfile
RAW StringLog
1 DnSaveLogfile
-1. DtSaveLogfile
#SAVEINITIAL
F DoSaveInitial
#SAVEPLOT
2 nPlotFile
y=0 MHD idl StringPlot
10 DnSavePlot
-1.0 DtSavePlot
-1.0 DxSavePlot
z=0 MHD idl StringPlot
10 DnSavePlot
-1.0 DtSavePlot
-1.0 DxSavePlot
#USERSWITCH
+source_expl +init +ic StringSwitch
#USERINPUTBEGIN -----------------------------------
#SOLARWINDH
0.00874 SWH_N_dim [/cc]
1.0868E5 SWH_T_dim [K]
417.07 SWH_Ux_dim [km/s]
0.0 SWH_Uy_dim [km/s]
0.0 SWH_Uz_dim [km/s]
7.17E-03 SWH_Bx_dim [nT]
0.0 SWH_By_dim [nT]
0.0 SWH_Bz_dim [nT]
#VLISW
0.06 VLISW_N_dim [n/cc]
6519.0 VLISW_T_dim
26.3 VLISW_Ux_dim
0.0 VLISW_Uy_dim
-2.3 VLISW_Uz_dim
0.3093 VLISW_Bx_dim
-0.080 VLISW_By_dim
-0.2988 VLISW_Bz_dim
#USERINPUTEND -------------------------------------
#END_COMP OH -----------------------------------------------------------------
#COUPLE2
OH NameComp1
PT NameComp2
1 DnCouple
-1.0 DtCouple
#STOP
20 MaxIter
-1.0 TimeMax
#RUN ##########################################################################
#TIMEACCURATE
T IsTimeAccurate
#SAVERESTART
T DoSaveRestart
-1 DnSaveRestart
5.0 year DtSaveRestart
! coupling with 1.0 year does not work
#COUPLE2
OH NameComp1
PT NameComp2
-1 DnCouple
0.1 year DtCouple
#STOP
-1 MaxIter
5.0 year TimeMax