Project

General

Profile

cosmo_NORWAY22_job.sh

Runscript (generating INPUT_* files) for 25km run - Andreas Dobler, 08/24/2016 11:18 AM

 
1
#!/bin/bash
2
#PBS -A mifa01mk
3
#PBS -m n
4
#PBS -N cosmo_no22
5
#PBS -l select=32:ncpus=32:mpiprocs=16
6
#PBS -o /work/dobler/COSMO/cclm_011/cclm_out.out
7
#PBS -j oe
8
#PBS -l walltime=90:00:00
9

    
10
cd /work/dobler/COSMO/cclm_011
11
rm -f YU* M_*
12

    
13
#################################################
14
# global settings
15
# DINI = Initialisation date (cold start)
16
# DSTA = Starting date of this run (restart, cont. run)
17
# DSTP = Stopping date of this run
18
# HOUTINC = Hours increment to write output files
19
#################################################
20
DINI=1968010100
21
DSTA=1969051706
22
DSTP=2001010100
23

    
24
HSTART=`/home/metno/dobler/bin/cfu get_hours $DINI $DSTA`
25
HSTOP=`/home/metno/dobler/bin/cfu get_hours $DINI $DSTP`
26
HOUTINC=6
27
HINC=`expr $HSTOP - $HSTART`
28
HRESINC=`expr $HINC / 23`
29
#################################################
30
# cat together the INPUT*-files
31
#################################################
32

    
33
cat > INPUT_ORG << end_input_org
34
 &LMGRID
35
  startlat_tot  = -7.78, startlon_tot  = -44.78,
36
  pollat=22.0,        pollon=-120.0,
37
  dlat=0.22,          dlon=0.22,
38
  ie_tot=248,         je_tot=123,        ke_tot=40,
39
 /
40
 &RUNCTL
41
  hstart = $HSTART.0,  hstop  = $HSTOP.0,   dt      = 120.0,    ydate_ini='$DINI',
42
  nprocx = 32,  nprocy = 16,   nprocio = 0,
43
  lphys      = .TRUE.,   luse_rttov    = .FALSE.,   luseobs = .FALSE.,  leps = .FALSE.,
44
  lreorder   = .FALSE.,  lreproduce    = .TRUE.,   itype_timing = 4,
45
  ldatatypes = .FALSE.,  ltime_barrier = .FALSE.,  ncomm_type=3,
46
  nboundlines=3, idbg_level = 2, lartif_data=.FALSE,
47
  ldfi=.FALSE., ldebug_io=.FALSE., lprintdeb_all=.FALSE.,
48
  hincmxt=24, hincmxu=24,
49
  ldump_ascii=.FALSE.,
50
  itype_calendar=0,
51
 /
52
 &TUNING
53
  c_soil   =   1.0,
54
  clc_diag =   0.5,
55
  crsmin   = 150.0,
56
  qc0      =   0.0,
57
  q_crit   =   4.0,
58
  qi0      =   0.0,
59
  rat_can  =   1.0,
60
  rat_lam  =   1.0,
61
  tur_len  = 500.0,
62
  v0snow   =  20.0,
63
  wichfakt =   0.0,
64
  tkmmin   =   1.0,
65
  tkhmin   =  0.35,
66
  rlam_heat=0.5249,
67
  uc1      =0.0626,
68
  radfac   =  0.50,
69
  soilhyd  =  1.62,
70
  fac_rootdp2= 0.9,
71
  entr_sc  =1.86E-04,
72
 /
73
end_input_org
74

    
75
cat > INPUT_INI  << end_input_ini
76
 &INICTL
77
  ndfi  = 2, 
78
  nfilt = 1,
79
  tspan = 600.0,
80
  dtbak = 15.0,
81
  dtfwd = 15.0,
82
  taus  = 600.0,
83
 /
84
end_input_ini
85

    
86
cat > INPUT_IO  << end_input_io
87
 &IOCTL
88
  lasync_io=.FALSE., ngribout=3,
89
  yform_read='ncdf',
90
  lbdclim=.TRUE.
91
  ymode_write='w  ',
92
  nhour_restart=$HSTART,$HSTOP,$HRESINC,
93
  ydir_restart='/work/dobler/COSMO/cclm_011/output_restart',
94
  ytunit_restart='d',
95
  yncglob_institution="Meteorologisk institutt, Oslo, Norway",
96
  yncglob_title="NORWAY22 with CCLM_5.0 dr. by MPI-ESM-LR_historical_r1i1p1",
97
  yncglob_source="cosmo_131108_5.00_clm6, int2lm_131101_2.00_clm4",
98
  yncglob_project_id="MET-DOWNSCALE",
99
  yncglob_experiment_id="MET011",
100
  ncglob_realization=1,
101
  yncglob_contact="andreasd@met.no",
102
  yncglob_references="http://www.clm-community.eu/",
103
 /
104
 &DATABASE
105
 /
106
 &GRIBIN
107
  ydirini='/work/dobler/COSMO/cclm_011/input',
108
  ydirhhl='/work/dobler/COSMO/cclm_011/input',
109
  ydirbd='/work/dobler/COSMO/cclm_011/input',
110
  lana_qi=.TRUE.,
111
  llb_qi=.TRUE.,
112
  hincbound=6,
113
  lchkini=.TRUE.,
114
  lan_t_so0=.TRUE.,
115
  lan_t_snow=.TRUE.,
116
  lan_t_cl=.TRUE.,
117
  lan_w_snow=.TRUE.,
118
  lan_w_i=.TRUE.,
119
  lan_w_cl=.TRUE.,
120
  lan_vio3=.TRUE.,
121
  lan_hmo3=.TRUE.,
122
  lan_plcov=.TRUE.,
123
  lan_lai=.TRUE.,
124
  lan_rootdp=.TRUE.,
125
  ytunitbd='d',
126
 /
127
 &GRIBOUT
128
  hcomb=$HSTART,$HSTOP,$HOUTINC,
129
  yvarml='ALHFL_S', 'ALWD_S', 'ALWU_S', 'ASHFL_S', 'ASOD_T', 'ASOB_T', 'ASOB_S', 'ATHB_T', 'ATHB_S', 'ASWDIFD_S', 'ASWDIFU_S', 'ASWDIR_S', 'CLCT', 'CLCL', 'CLCM', 'CLCH', 'CLC_CON', 'DURSUN', 'PMSL', 'PS', 'QV_2M', 'T_2M', 'TMAX_2M', 'TMIN_2M', 'T_G', 'U_10M', 'V_10M', 'RELHUM_2M', 'ALB_RAD', 'SNOW_CON', 'SNOW_GSP', 'RAIN_CON', 'RAIN_GSP', 'TOT_PREC', 'H_SNOW', 'RUNOFF_G', 'RUNOFF_S', 'VABSMX_10M', 'VMAX_10M', 
130
  yvarpl=' ',
131
  yvarzl=' ',
132
  luvmasspoint=.TRUE.,
133
  lcheck=.FALSE.,
134
  lwrite_const=.FALSE.,
135
  ytunit='d',
136
  yform_write='ncdf',
137
  ydir='/work/dobler/COSMO/cclm_011/output_restart',
138
 /
139
 &GRIBOUT
140
  hcomb=$HSTART,$HSTOP,$HOUTINC,
141
  yvarml=' ',
142
  yvarpl='QV','T','U','V',
143
  plev=200.,500.,800.,850.,
144
  yvarzl=' ',
145
  luvmasspoint=.TRUE.,
146
  lcheck=.FALSE.,
147
  lwrite_const=.FALSE.,
148
  ytunit='d',
149
  yform_write='ncdf',
150
  ydir='/work/dobler/COSMO/cclm_011/output_restart_plev',
151
 /
152
 &GRIBOUT
153
  hcomb=$HSTART,$HSTOP,3,
154
  yvarml='FRESHSNW','PP','QC','QI','QR','QS','QV','QV_S','T','T_S','T_SNOW','T_SO','U','V','W','W_I',
155
         'W_SNOW','W_SO',
156
  yvarpl=' ', 
157
  yvarzl=' ',
158
  luvmasspoint=.FALSE.,
159
  lcheck = .FALSE.,
160
  lwrite_const=.TRUE.,
161
  ytunit='d',
162
  yform_write='ncdf',
163
  ydir='/work/dobler/COSMO/cclm_011/output_nesting',
164
 /
165

    
166
end_input_io
167

    
168
cat > INPUT_DYN << end_input_dyn
169
 &DYNCTL
170
  l2tls=.true.,
171
  alphaass=0.53,
172
  rlwidth=250000.0,
173
  l_diff_Smag=.FALSE., 
174
  y_scalar_advect='Bott2', 
175
  ldyn_bbc=.TRUE.,
176
  nrdtau=6,
177
  itype_fast_waves=1,
178
  itype_bbc_w=1,
179
  iadv_order=3,
180
 /
181
end_input_dyn
182

    
183
cat > INPUT_PHY << end_input_phy
184
 &PHYCTL
185
  ico2_rad=8,
186
  ke_soil=9,
187
  czml_soil=0.005,0.025,0.07,0.16,0.34,0.7,1.42,2.86,5.74,11.5,
188
  czbot_w_so=4.0,
189
  nincconv=2,
190
  limpltkediff=.FALSE.,
191
  lseaice=.FALSE.,
192
  llake=.FALSE.,
193
  lsso=.TRUE.,
194
  itype_heatcond=2,
195
  itype_evsl=3,
196
  itype_albedo=2,
197
  itype_root=2,
198
 /
199
end_input_phy
200

    
201
cat > INPUT_DIA << end_input_dia
202
 &DIACTL
203
 /
204
end_input_dia
205

    
206
cat > INPUT_EPS << end_input_eps
207
 &EPSCTL
208
 /
209
end_input_eps
210

    
211
cat > INPUT_ASS << end_input_ass
212
 &NUDGING
213
  lnudge  =.false.,
214
end_input_ass
215

    
216
#################################################
217
# run the program
218
#################################################
219

    
220
#################################
221
module load netcdf mpt
222

    
223
mpiexec_mpt cclm.exe
224

    
225
#################################################
226
# cleanup
227
#################################################
228

    
229
#rm -f INPUT_ORG INPUT_IO INPUT_DYN INPUT_DIA INPUT_PHY INPUT_INI
230
#rm -f INPUT_ASS INPUT_SAT INPUT_EPS
231