-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathps1d.dat
81 lines (79 loc) · 2.47 KB
/
ps1d.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
ps1d
###############
nz 600 # number of vertical grid points
ztop 2400 # height of the domain top (meters)
ntime 3456000 # total time steps
dt 0.1 # time step (sec)
##### sea ice / snow model parameters ##################
nl 25 # number of levels in ice layer
ice_h 2. # ice layer thickness (meters)
nls 15 # number of levels in snow layer
snow_h 0.3 # snow layer thickness (meters) (set to 0 if no snow)
frac 1 # sea ice fraction
#### vertical turbulent diffusion ###############
# options for iodif:
# 3 - local (METRAS)
# 4 - nonlocal Troen and Mahrt (1986)
# 5 - nonlocal Lupkes and Schlunzen (1996)
# 6 - nonlocal + entrainment Noh et al. (2003)
# 7 - nonlocal + entrainment Lock (2000)
# 8 - local (INM climate model)
# 9 - linear fluxes (mixed-layer model)
#################################################
iodif 3 # turbulent diffusion (see options above)
implicit .true. # implicit scheme for diffusion, otherwise expl scheme is used
rad_par 2 # radiation (0 - off; 1 - prescribed(DYCOMS); 2 - Goddard
vadv .true. # vertical advection
tfix 0 # constant surface temperature
ifhle 0 # constant surface fluxes (1/0)
seaice 1 # sea ice model coupled to atm (1/0)
qif 0 # moisture on/off (1/0)
ifwr 0 # cloud water and rain on/off (1/0)
ifmf 0 # cloud snow on/off (1/0)
icetime 10 # specific for cold air outbreaks
###############
pa 100000. # pressure at the surface (Pa)
ug 15. # u-geostrophic
vg 0. # v-geostrophic
phi 80. # latitude
ts 241.15 # surface temperature
###############
#zgrid # uncomment this to read the vertical grid from this file
20
50
90
140
200
270
350
440
540
##########
inifile 1 # initial profiles: 1 - from file (see below), 0 - set in the code
##########
profile # theta u v q
256 1 1 0 0
10 10 10 10 # number of levels -1 in initial profs
0 256.8 15. 0 0
200 257.5 15. 0 0
400 263.2 15. 0 0
600 266.5 15. 0 0
800 268.3 15. 0 0
1000 270. 15. 0 0
2000 276.5 15. 0 0
3000 283 15. 0 0
4000 289.5 15. 0 0
14000 359.5 15. 0 0
####################################################################
iftf .false.
############
tforcing # K/hour
6
0 -0.18
100 -0.1875
200 -0.1875
400 -0.1333
600 -0.0417
800 0.0625
1000 0.0625
end