-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathnh3d.dat
263 lines (262 loc) · 12.2 KB
/
nh3d.dat
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
253
254
255
256
257
258
259
260
261
262
263
nh3d
# 1.Nh3d data file, P Miranda August 2006, [email protected]
# 2. Corrected by V. Stepanenko, June, 2007, [email protected]
# this is an example of nh3d.dat
# note that all lines belonging to optional parameters may be
# omitted and then default values are taken
# in this file after the value of a given parameter its default
# value is shown between (), when available.
# order of lines is irrelevant
#
verbose 0 (0) level of screen output: -1(none) 0(a bit) 1(a lot) 2(all) 3 (debug)
####################################################################
#following parameters MUST be given:
#grid (nx,ny,ns)
nx 241 ! nx-1 must be fft'able (2**n * 3**k * 5**m)
ny 161 ! ny-1 must be fft'able
ns 47
dx 5000.
#6200.
dy 5000.
#5800.
dt 2.
ntime 86400 !number of time steps of the experiment
ptop 40000. !pressure at the top of the model
#####################################################################
sigmagrid
0.
0.041232639
0.081597222
0.12109375
0.159722222
0.197482639
0.234375
0.270399306
0.305555556
0.33984375
0.373263889
0.405815972
0.4375
0.468315972
0.498263889
0.52734375
0.555555556
0.582899306
0.609375
0.634982639
0.659722222
0.68359375
0.706597222
0.728732639
0.75
0.770399306
0.789930556
0.80859375
0.826388889
0.843315972
0.859375
0.874565972
0.888888889
0.90234375
0.914930556
0.926649306
0.9375
0.947482639
0.956597222
0.96484375
0.972222222
0.978732639
0.984375
0.989149306
0.993055556
0.996527778
1
#simpleatmos 273. 1.0e-4 10 0 0 1.0e5 !simplest possible profile: pot.temp_surface,N**2,U,V,RH,p_surface (all constant in vertical)
timeini 20001101060000 ! timeini must be set before profile
profile
260 1 1 0 0
12 12 12 12
0 260 0 0.0 0 -10. 0 1.1
1000 265 1000 0.0 1000 -10. 1000 0.9
2000 270 2000 0.0 2000 -10. 2000 0.7
3000 275 3000 0.0 3000 -10. 3000 0.3
4000 280 4000 0.0 4000 -10. 4000 0
5000 285 5000 0.0 5000 -10. 5000 0
6000 290 6000 0.0 6000 -10. 6000 0
7000 295 7000 0.0 7000 -10. 7000 0
8000 300 8000 0.0 8000 -10. 8000 0
9000 305 9000 0.0 9000 -10. 9000 0
10000 310 10000 0.0 10000 -10. 10000 0
11000 315 11000 0.0 11000 -10. 11000 0
12000 320 12000 0.0 12000 -10. 12000 0
###################################################################
fngrd 'results/svalbard_5/test' ! output grid files (grd format)
fnmap 'map/svalbard_5' ! input grid files for surface parameters (soil model)
#arctic_60km_12x_18ice_18water
#fnmap 'map/lake_map' ! input grid files for surface parameters (soil model)
fnout 'results/svalbard_5/test' ! output data (unformatted) plus other output fields (e.g. profiles)
fnrad 'series.rad'
#####################################################################
#optional parameters (all have default values)
divrefstate .false. (.false.) !Choose false (true is not yet operational)
ifluxcor 0 (1) !what to do about total mass (0) let it go; (1) preserve mass - generaly a good idea; (2) balance mass (not yet operational, need divrefstate=.true.)
adjust .false. (.false.) !advanced, if .true. uses Room et al (2001) adjustment scheme
tfct .true.
centered .false.
dst3 .false.
xlatit 0. (55.) !latitude in degrees fcor is computed accordingly
fcor -999. (-999.)
fcori 1.432e-4
ug0 0.0
vg0 -10.
ifsoil 0 (0) !switch for soil model (0-off, 1-on, 2-simple soil without heat flux)
tfix 1 switch for simple surface with fixed temperature (if ifsoil=0, if radpar=0 and if iflake=0)
z0hz0 .1 (1.) z0h/z0 only used if ifsoil=1
isqv 1 (0) switch for water vapour (0-off, 1-on)
isqc 0 (0) switch for cloud water (0-off, 1-on)
isqr 0 (0) switch for rain water (0-off, 1-on)
isqi 0 (0) switch for cloud ice and snow in clouds (0-off,1-on)
# diffusion parameters
iodif 13 (0) diffusion switch (0-no, 1- constant, 2- def depend., 12-deformation dep but only vertical)
difcof 0. (0.) constant diffusion (to add if iodif=2)
rikey 1. (1.) if 1. diffusion depends also on the richardson number
uvdif 1. (1.) apply diffusion to the reference state wind (us,vs)
tsdif 1. (1.) apply diffusion to the reference state theta profile (pts)
cdcoef 0. (0.) drag coefficient for a bulk surface model
iophis 3
xice 25
ifperturb 0
ifdirichlet 0
##########
uvpert .false. ! switcher for perturbation of wind field
pnudg 1800 ! duration of nudging (in timesteps)
tau_nudg 0.001 ! time scale of nudging
nperturb 2 ! time step for the start of nudging
f_chord 30 ! length scale of perturbation (in kilometers)
xm_pert 7.5 ! amplitude of perturbation (in ms-1)
#########
iobuux 1
iobuuy 0
iobvvx 0
iobvvy 1
iobptx 0
iobpty 0
iobppx 1
iobppy 1
#########
#####################################################################
#topography iomtyp,hmount,xmount,ymount,hxwid,hywid,hexp
# if iomtyp=4 (default) the grid file [fnmap].top contains the orography
# and the rest of the parameters are not needed
# iomtyp=1:
# h=hmount*(1+((x-xmount)/hxwid)**2)**(-hexp)
# iomtyp=2:
# h=hmount*(1+((x-xmount)/hxwid)**2+((y-ymount)/hywid)**2)**(-hexp)
#topography 4
#MOUNTAIN
#1 1000 32000. 8000. 5000. 5000. 1.
#
#####################################################################
#smoothing parameters
numsmo 1 (0) interval beteween smoothing, 0- no smoothing
perdam .true. (.true.) smooth perturtabtions only
dohsmo .true. (.true.) horizontal smoothing on (if numsmo>0)
dovsmo .false. (.true.) vertical smoothing on (if numsmo>0)
dovsmt .false. (.true.) vertical smoothing for pt on (if numsmo>0)
hdampm 0.015 (0.015) smoothing parameter for hsmo (horizontal filter) value below the absorption layer
hdampt 0.0625 (0.0625) smoothing parameter for hsmo maximum value at the topo of the absorption layer
vdamp 0.001 (0.001) smoothing parameter for vsmo (vertical filter) value below the absorption layer
vdampt 0.001 (0.001) smoothing parameter for vsmt maximum value at the topo of the absorption layer
#absorption layer
taum 1.0e8 (1.0e8) relaxation time at the bottom of the sponge layer
ntaut 100 (100) relaxation time (in timesteps) at the top of the sponge layer
zm 8000. (8000.) bottom of the sponge lay im meters
zt 15000. (18000.) top of the sponge lay
#lateral sponge (this is needed for changing b.c. or to guarantee its maintenance in time, if -1 it is not used)
nxsponge1 -1 (-1)
nxsponge2 -1 (-1)
nysponge1 -1 (-1)
nysponge2 -1 (-1)
ntauspon 100 ! relaxation time (in timesteps) at the end of the sponge layer
#
#####################################################################
#output switchs
forout .true. (.true.) produce formatted output on file [fngrd].out
grdout .true. (.false.) formatted output is in grd format [fngrd]*.#
iotfor 3600 (999999) interval between formatted outputs (last one always on), if forout=.true.
iotunf 7200 (999999) interval between unfformatted outputs (last one always on), if unfout=.true.
unfout .false. (.false.) produce unfformatted output on files [fnout].u[*]
drgprt .false. (.false.) produce drag output on file [fngrd].dxy
iotdrg 4800 (1) interval between drag outputs
iotscrout 10 (1) interval between screen outputs
domom .true. (.false.) produce gridded output at z=const on file [fnout].mom
mxspan .true. (.false.) !if .true. the domain for the drag computation is centred at x=y=0 (that is at xmount, ymount for the analytical topography) this is needed for drag calculation with analytical mountains.
iotduv 7200 (20) interval between momentum flux outputs
nlev 300 (ns) number of levels of mom flux and grid output
dzlev 20. (500.) distance between levels of mom flux and grid output
proout .true. (.true.) produce profile output on file [fnout].pro
iniout .true. (.true.) produce formatted output at first time step
#chose variable group for formatted output:
norout .true. (.true.) " u,v,w,pt,theta,qv,cond,evap,qc,qr (water variables if they are switched on). u,v,w come out in their grid and also in the thermodynamic grid (um,vm,wm) interpolated.
surout .true. (.false) " (all surface variables)
morout .true. (.false) " phi,wsig(vertical vel in sigma coord),rh(relative humidity)
refout .true. (.false.) " us,vs,pts,phis (reference state, maybe interesting in ref state is varying in time in samp_for.ref, and ifuprefs=1)
difout .true. (.false) " Diffusion variables
driout .true. (.false) " Richardson number
inrout .false. (.false) " (debug)
fluout .false. (.true) " (debug)
ppbout .false. (.false) " (debug)
duvout .false. (.false) " (debug)
uubout .false. (.false) " (debug)
ptbout .false. (.false) " (debug)
ddiout .false. (.false) " (debug)
defout .false. (.false) " (debug)
desout .false. (.false) " (debug)
elfout .false. (.false) " (debug)
elgout .false. (.false) " (debug)
eneout .false. (.false) " (debug)
echeck .false. (.false) " (debug)
masout .false. (.false.) " (debug)
# output selection (planes)
iomod 0 (iomod=0 means that the following selection is on, iomod=1 all sigma planes for debug, =2 all x planes, =3 all y planes)
ixp -1,-1,-1,-1,-1,-1,-1,-1,-1,-1 select x-planes for formatted output
iyp 20,40,60,-1,-1,-1,-1,-1,-1,-1 select y-planes for formatted output
isp 46,-1,-1,-1,-1,-1,-1,-1,-1,-1 select sigma-planes for formatted output
zplan -1,-1,-1,-1.,-1.,-1.,-1.,-1.,-1.,-1. select z-planes for formatted output
zconst .true. (.true.) vertical cross-sections are interpolated for a xz or yz grid
outpoint 3 30 46 !location of point output goes to file [fngrd]xnnnynnnsnnn.loc, several points may be selected each with one command
#outprof 32 43 !location of profile output goes to file [fngrd]pfxnnnynnn.loc, several points may be selected each with one command
#outprof 3 19 !location of profile output goes to file [fngrd]pfxnnnynnn.loc, several points may be selected each with one command
#outprof 6 20 !location of profile output goes to file [fngrd]pfxnnnynnn.loc, several points may be selected each with one command
#outprof 6 21 !location of profile output goes to file [fngrd]pfxnnnynnn.loc, several points may be selected each with one command
#####################################################################
#other options (in general not changed)
itepsuf .true. (.true.) ! only relevant with rotation
filt 0.1 (0.1) time filter (aselin)
ppfilt 0.1 (0.1) time filter (aselin) for psuf
nupsur 1 ! interval between updates of ref state at the surface
nuppts 1 ! interval between updates of ref state; this is not relevant if ifuprefs=1: in this case the update is always made
nupdif 1 ! interval between updates of diffusion variables
#####################################################################
# NEW SWITCHES & PARAMETERS (V. Stepanenko, May 2007)
#####################################################################
iflake 0 (0) ! switch for lake model (0 - ISBA lake, 1 - LAKE model for lake)
nlakecall 60 (1) ! the period for updating lake parameters, timesteps (if iflake is set to 1)
#####################################################################
# radiation
# radpar 1 - prescribed fluxes of radiation at the surface, no heating in the atmosphere
# radpar 2 - parameterized fluxes of radiation at the surface, no heating in the atmosphere
# radpar 3 - Clirad parameterization (full parameterization of radiation processes in atmosphere)
# nradcall - number of time steps that separate two subsequent calls of Clirad (if radpar set to 3)
nradcall 60 ! the period for updating radiation fluxes, timsteps
radpar 0
longwave 0 ! relevant if radpar = 3, 1 - longwave parameterization is on, 0 - off
shortwave 0 ! relevant if radpar = 3, 1 - shortwave parameterization is on, 0 - off
rat_par 1 ! relevant if radpar = 2
tau_rg 2. ! relevant if radpar = 2
#####################################################################
# SWITCHES FOR PASSIVE TRACER (D. Mikushin & V. Stepanenko, May 2007)
#####################################################################
impur 0
#####################################################################
end