source: codes/icosagcm/trunk/param_sets/make_rundefs.sh @ 217

Last change on this file since 217 was 217, checked in by dubos, 10 years ago

Removed "automatic" physics + bugfix

  • Property svn:executable set to *
File size: 8.1 KB
Line 
1#!/bin/bash
2
3function print()
4{
5    # prints one comment line followed by any number of "key = value" lines
6    echo "# $1 : $2 (default=$3)"
7    shift; shift; shift;
8    while (($# > 0)) ; do
9        echo "$1 = $2"
10        shift; shift;
11    done
12    echo
13}
14
15################################### Planet #######################################
16 
17function radius { print 'Planetary radius' real 6371220 radius $1 
18}
19function g      { print Gravity real 9.80616 g $1 
20}
21function omega { print 'Planetary rotation rate' real 7.292E-5 omega $1 
22}
23function kappa { print 'kappa=Rd/cpp' real 0.2857143 kappa $1 
24}
25function cpp     { print 'Heat capacity at constant pressure' real 1004.70885 cpp $1
26}
27function preff { print 'Reference pressure' real 101325 preff $1 
28}
29function scale_factor { print 'Scale factor for small planet experiment' real 1. scale_factor $1 
30}
31
32##################################### Mesh #####################################
33
34function nbp   { print 'Number of subdivisions on a main triangle' integer 40 nbp $1 
35}
36function split { print 'Sub splitting of main rhombus' integer 1 nsplit_i $1 nsplit_j $2 
37}
38function llm   { print 'Number of vertical layers' integer 19 llm $1 
39}
40function disvert { print 'Vertical grid' '[std|ncar|ncarl30]' std disvert $1 
41}
42function ncar_dz { print 'Model layer thickness in m, used if disvert=ncar' real 300 ncar_dz $1
43}
44function optim_it { print 'Mesh optimisation : number of iterations' integer 0 optim_it $1 
45}
46
47################################# Numerics ###################################
48
49function scheme { print 'Time-stepping scheme' '[euler|leapfrog_matsuno|runge_kutta]' runge_kutta scheme $1 
50}
51function time_step { print 'Time step in s' real 480 dt $1 
52}
53function matsuno { print 'Matsuno period' integer 5 matsuno_period $1 
54}
55function itau_adv { print 'Advection called every itau_adv time steps' integer 2 itau_adv $1 
56}
57function nqtot { print 'Number of tracers' integer 1 nqtot $1
58}
59function caldyn { print 'Equations solved' '[gcm|adv]' gcm caldyn $1 
60}
61function caldyn_conserv { print 'Discretization of Coriolis term' '[energy|enstophy]' energy caldyn_conserv $1 
62}
63function caldyn_eta { print 'Vertical coordinate' '[eta_mass|eta_lag]' mass caldyn_eta $1 
64}
65function boussinesq { print 'Boussinesq approximation, requires eta_lag' '[.FALSE.|.TRUE.]' .FALSE. boussinesq $1
66}
67############################### Time and output ###############################
68
69function time_style() { print 'Time style' '[none|lmd]' lmd time_style $1 
70}
71function ndays() { print 'Number of days' integer ?? ndays $1 
72}
73function run_length() { print 'Run length in s' real ?? run_length $1 
74}
75function day_step { print 'Number of time steps per day' integer ?? day_step $1 
76}
77function write_period { print 'Interval in s between two outputs' integer ?? write_period $1 
78}
79function out_pressure_level { print 'Pressure level to which interpolate T' real 0 out_pression_level $1
80}
81
82################################ Physics ######################################
83
84function etat0 { 
85    OPTS=$(printf '\n#   [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear|dcmip3]')
86    print 'Initial state' "$OPTS" jablonowsky06 etat0 $1
87}
88function dissip
89{
90    print 'Dissipation time for grad(div)' real 5000 tau_graddiv $1
91    print 'Exponent of grad(div) disspation' integer 1 nitergdiv $2
92    print 'Dissipation time for curl(curl)' real 5000 tau_gradrot $3
93    print 'Exponent of curl(curl) disspation' integer 1 nitergrot $4
94    print 'Dissipation time for div(grad)' real 5000 tau_divgrad $5
95    print 'Exponent of div(grad) disspation' integer 1 niterdivgrad $6
96}
97function guided_type { print 'Type of guiding terms' '[none|dcmip1]' none guided_type $1 
98}
99
100function physics { print 'Physics package' '[none|held_suarez|dcmip]' none physics $1
101}
102
103function itau_physics { print 'Call physics every itau_physics time steps' 'integer' 1 itau_physics $1
104}
105
106#################################### DCMIP #################################
107
108function dcmip1_def()
109{
110    OPTS=$(printf '\n#    [const|slotted_cyl|cos_bell|dbl_cos_bell_q1|dbl_cos_bell_q2|complement|hadley|dcmip11]')
111    print 'DCMIP 1.1 advection test, initial tracer field' "$OPTS" cos_bell dcmip1_shape $1
112    print 'DCMIP 1.1 advection test, wind field' '[solid|deform|hadley]' deform dcmip1_shape $2
113}
114function dcmip_refs()
115{
116    print 'Reference temperature T0 for DCMIP test cases' real 300 ncar_T0 $1
117    print 'Reference pressure p0 for DCMIP test cases, used by disvert if disver=ncar' real 1e5 ncar_p0 $2
118    print 'Exponent for B(eta), used by disvert if disvert=ncar' integer 1 ncar_disvert_c $3
119}
120function dcmip_physics()
121{
122    physics dcmip
123    print 'DCMIP physics' '[0|1]' 1 dcmip_physics $1
124}
125function dcmip4_testcase { print 'DCMIP 4 Test case' '[1|2]' 1 dcmip4_testcase $1 
126}
127
128####################### Functions generating parameter files ####################
129
130function params {
131    # executes any number of "command arg" commands
132    while (($# > 0)) ; do
133        $1 $2
134        shift; shift;
135    done
136}
137
138function group { 
139    # prints one separation line then executes any number of "command arg" commands
140    echo
141    echo "#---------------- $1 ----------------"
142    echo
143    shift
144    params $*
145}
146
147function williamson91()
148{
149    group 'Mesh' nbp 40 llm 1 optim_it 100 ; split 2 2
150    group 'Numerics' time_step 480 caldyn_conserv enstrophy
151    group 'Time and output' out_pressure_level 85000 write_period 21600 run_length 864000
152    group 'Physical parameters' radius 6.37122e6 g 9.80616 caldyn_eta eta_lag boussinesq .TRUE. etat0 williamson91.6
153    dissip 5000 2 5000 2 5000 2   
154}
155
156function held_suarez()
157{
158    group 'Mesh' nbp 40 llm 19 disvert std optim_it 1000 ; split 2 2
159    group 'Numerics' itau_adv 1 day_step 180
160    group 'Time' time_style lmd ndays 1200 write_period 86400
161    group 'Planet' radius 6.371e6 g 9.8 omega 7.292e-5 kappa 0.2857143 cpp 1004 preff 1e5
162    group 'Physical parameters' etat0 held_suarez physics held_suarez
163    dissip 50000 2 50000 2 50000 2
164}
165
166function const_dcmip
167{
168    group 'Planet' radius 6.37122e6 g 9.80616 kappa 0.2857143 cpp 1004.5 preff 1e5 omega $1 
169    group 'Constants' ; dcmip_refs 300 1e5 1
170}
171
172function dcmip11
173{
174    echo 'INCLUDEDEF = const.def'
175    group 'Mesh' nbp 80 llm 60 disvert ncar ncar_dz 200 optim_it 100 ; split 2 2
176    group 'Numerics' caldyn adv time_step 180
177    group 'Time' run_length 1036800 write_period 86400
178    group 'Physical parameters' guided_type dcmip1 ; dcmip1_def dcmip11 deform
179}
180
181function dcmip21
182{
183    echo 'INCLUDEDEF = const.def'
184    group 'Mesh' nbp 64 llm 60 disvert ncar ncar_dz 500 optim_it 100 ; split 2 2
185    group 'Numerics' time_step 300
186    dissip 1800 1 1800 2 1800 2
187    group 'Time' run_length 3.6e6 write_period 50000
188    group 'Physical parameters' scale_factor 500 guided_type dcmip1 etat0 dcmip2_schaer_noshear
189}
190
191function dcmip3
192{
193    const_dcmip 0
194    group 'Mesh' nbp 80 llm 10 disvert ncar ncar_dz 1000 optim_it 100 ; split 2 2
195    group 'Numerics' time_step 240
196    dissip 1800 1 1800 2 1800 2
197    group 'Time' run_length 4.5e5 write_period 12500
198    group 'Physical parameters' omega 0 scale_factor 125 etat0 dcmip3
199}
200
201function dcmip4x
202{
203    echo 'INCLUDEDEF = const.def'
204    group 'Mesh' nbp 80 llm 30 disvert ncarl30 optim_it 100 ; split 2 2
205    group 'Numerics' time_step 240
206    dissip 1800 1 1800 2 1800 2
207    group 'Time' run_length 3.6e6 write_period 21600
208    group 'Physical parameters' nqtot 2 etat0 dcmip4 dcmip4_testcase $1
209}
210
211function dcmip42
212{
213    dcmip4x 2
214    dcmip_physics 1
215}
216
217function dcmip51
218{
219    echo 'INCLUDEDEF = const.def'
220    group 'Mesh' nbp 64 llm 30 disvert ncarl30 optim_it 1000 ; split 2 2
221    group 'Numerics' time_step 300
222    dissip 10800 1 10800 2 10800 2
223    group 'Time' run_length 864000 write_period 7200
224    group 'Physical parameters' nqtot 1 etat0 dcmip5 dcmip_physics 0 itau_physics 6
225}
226   
227########## Main program : generate all *.def files in param_sets ##############
228
229function mkcd { cd $ROOT ; mkdir -p $1 ; cd $1
230}
231
232set -o errexit
233ROOT=$PWD
234
235mkcd climate/Held_Suarez
236held_suarez > run.def
237
238mkcd shallow_water/williamson91
239williamson91 > run6.def
240
241mkcd dcmip2012
242const_dcmip 7.292e-5 > const.def
243dcmip11 > run11.def
244dcmip21 > run21.def
245dcmip3 > run3.def
246dcmip4x 1 > run41.def
247dcmip42 > run42.def
248dcmip51 > run51.def
Note: See TracBrowser for help on using the repository browser.