source: codes/icosagcm/branches/SATURN_DYNAMICO/TEST/run.def @ 220

Last change on this file since 220 was 220, checked in by ymipsl, 10 years ago

Creating temporary dynamico/lmdz/saturn branche

YM
File size: 8.1 KB
Line 
1###########################################################################
2### RUN.DEF ### RUN.DEF ### RUN.DEF ### RUN.DEF ### RUN.DEF ### RUN.DEF ###
3###########################################################################
4### A MERGED VERSION TO BE USED WITH OLD OR RECENT LMDZ CORE
5### i.e. compiles gcm.def+run.def in LMDZ5 with run.def in old planeto LMDZ
6### (*) indicates same parameters in old and recent LMDZ core
7###########################################################################
8
9###########################################################################
10### NUMERICAL INTEGRATIONS
11###########################################################################
12### nday --> number of days of integration (*)
13### iperiod --> period for the Matsuno step (in dynamical steps) (*)
14### day_step --> number of steps per day (multiple of iperiod) (*)
15### purmats --> temporal integration: Matsuno or Matsuno-leapfrog? (*)
16### grireg --> regular grid or not (*)
17### adjust --> flag to activate load optimization
18### use_filtre_fft --> flag to activate FFT filtering
19###########################################################################
20nday = 100
21iperiod = 5
22day_step = 280
23purmats = .false.
24grireg = .true.
25adjust = n
26use_filtre_fft = n
27
28###########################################################################
29### WHAT RELATES TO PHYSICS
30###########################################################################
31### iphysiq --> period for calling physics (in dynamical steps) (*)
32### ecritphy --> frequency for outputs in diagfi (in dynamical steps) (*)
33### planet_type --> the planet
34### disvert_type --> type of vertical coordinates (2 for planeto)
35### hybrid --> use hybrid coordinates or not? (*)
36### iflag_phys=0 --> no physics (e.g. shallow water)
37### iflag_phys=1 --> physics (e.g. phylmd)
38### iflag_phys=2 --> simple newtonian relaxation
39### physic --> flag with or without physics (*)
40###########################################################################
41iphysiq = 140
42ecritphy = 2800
43planet_type = generic
44disvert_type = 2
45hybrid = .false.
46iflag_phys = 1
47physic = .true.
48
49###########################################################################
50### INITIALIZATION and NUDGING
51###########################################################################
52### read_start --> with or without start files ?
53###                (without: initialization by iniacademic)
54### raz_date=0 --> use the date in start file
55### raz_date=1 --> erase the date and start from scratch
56### ok_guide --> flag for nudging towards winds (needs relax.nc)
57### facwind --> time constant for relaxation
58###             5.e8 medium 1.e7-5.e6 strong 1.e6 very fast
59###########################################################################
60read_start = y
61raz_date = 0
62ok_guide = n
63facwind = 1.e6
64
65###########################################################################
66### HORIZONTAL DISSIPATION
67###########################################################################
68### dissip_period --> period call to dissipation (in dynamical steps)
69### - must be a multiple of iperiod
70### - 0 means automatic computation
71### - use idissip instead if old LMDZ core (idissip ignored in LMDZ5)
72### dissip_fac_mid --> enhancement factor in middle atmosphere
73### - set to 1 to prevent dissipation enhancement
74### lstardis --> star or not star dissipation operator (*)
75### nitergdiv --> number of iterations for gradiv operator (*)
76### nitergrot --> number of iterations for nxgradrot operator (*)
77### niterh --> number of iterations for divgrad operator (*)
78### tetagdiv --> dissipation timescale of smallest wvl: u,v (gradiv) (*)
79### tetagrot --> dissipation timescale of smallest wvl: u,v (nxgradrot) (*)
80### tetatemp --> dissipation timescale of smallest wvl: h (divgrad) (*)
81### coefdis --> coefficient for gamdissip (*)
82###########################################################################
83dissip_period = 5
84dissip_fac_mid = 2.
85lstardis = .true.
86nitergdiv = 1
87nitergrot = 2
88niterh = 2
89#tetagdiv = 5400. 8000. 160000. 1000. 2500. 4000. 160000.
90tetagdiv = 1000
91#tetagrot = 7200. 8000. 160000. 1000. 2500. 4000. 160000.
92tetagrot = 2000
93#tetatemp = 7200. 8000. 160000. 1000. 2500. 4000. 160000.
94tetatemp = 2000
95
96###########################################################################
97### ENHANCEMENT OF HORIZONTAL DISSIPATION IN UPPERMOST LEVELS
98###########################################################################
99### - not activated if ok_strato=n
100### - dissip_fac_up --> enhancement factor
101### - dissip_deltaz --> altitude range of transition (km)
102### - dissip_hdelta --> scale height at transition (km)
103### - dissip_pupstart --> pressure at bottom of transition region (Pa)
104###########################################################################
105ok_strato = y
106dissip_fac_up = 2.
107dissip_deltaz = 150.
108dissip_hdelta = 50.
109dissip_pupstart = 50.
110
111###########################################################################
112### DAMPING ("SPONGE") LAYER AT MODEL TOP
113###########################################################################
114### iflag_top_bound=0 for no sponge
115### iflag_top_bound=1 for sponge over 4 topmost layers
116### iflag_top_bound=2 for sponge from top to ~1% of top layer pressure
117### mode_top_bound=0 for no sponge
118### mode_top_bound=1 for u,v --> 0
119### mode_top_bound=2 for u,v --> zonal mean
120### mode_top_bound=3 for u,v,h --> zonal mean
121### tau_top_bound --> damping frequency at last layer
122### e-5 medium / e-4 strong yet reasonable / e-3 very strong
123###########################################################################
124iflag_top_bound = 2
125mode_top_bound = 1
126tau_top_bound = 1.e-4
127
128###########################################################################
129### OUTPUT FROM DYNAMICS
130###########################################################################
131### prt_level --> level of control prints (1-10)
132### ok_dyn_ins --> flag instantaneous outputs (dyn_hist.nc)
133### iecri --> frequency of instantaneous outputs (in dynamical steps)
134### ok_dyn_ave --> flag averaged outputs (dyn_hist_ave.nc)
135### periodav --> period of output of averages in dynamics (in days)
136### ok_dynzon --> flag outputs dynzon
137### NB: iecri and periodav used to be in old LMDZ core although useless
138### iconser --> frequency of outputs of control variables (in dyn steps)(*)
139###########################################################################
140prt_level = 1
141ok_dyn_ins = n
142iecri = 480
143ok_dyn_ave = n
144periodav = 1
145ok_dynzon = n
146iconser = 120
147
148###########################################################################
149### ZOOMED GRID (same parameters in old and recent LMDZ core)
150###########################################################################
151### clon --> longitude of center of zoom (deg) (*)
152### clat --> latitude of center of zoom (deg) (*)
153### grossismx --> zoom factor in longitude (*)
154### grossismy --> zoom factor in latitude (*)
155### dzoomx --> longitude extent of zoom zone (fraction of total area) (*)
156### dzoomy --> latitude extent of zoom zone (fraction of total area) (*)
157### taux --> "stifness" of zoom in longitude (*)
158### tauy --> "stifness" of zoom in latitude (*)
159### fxyhypb --> f(y) hyperbolic (T) or sinusoidal (F) (*)
160### ysinus --> f(y) = sin(lat) (T) or lat (F) (*)
161###########################################################################
162clon = 0.
163clat = 0.
164grossismx = 1.
165grossismy = 1.
166dzoomx = 0.
167dzoomy = 0.
168taux = 2.
169tauy = 2.
170fxyhypb = .false.
171ysinus = .false.
172
173###########################################################################
174### LEFTOVERS
175### calend=earth_360d ## specific calendar
176### resetvarc=n   ## reinit des variables de controle
177### iflag_trac=0  ## avec ou sans traceurs
178### dayref=10120  ## initial date
179### anneeref=2000 ## year of initial state
180### --- old version of spong layer
181## Avec sponge layer
182#  callsponge  = .true.
183## Sponge:  mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
184#  mode_sponge= 2
185## Sponge:  hauteur de sponge (km)
186#  hsponge= 50
187## Sponge:  tetasponge (secondes)
188#  tetasponge = 50000
189
190###########################################################################
191### INCLUDE OTHER DEF FILES (physics, specific settings, etc...)
192###########################################################################
193INCLUDEDEF=callphys.def
194#INCLUDEDEF=modif.def
Note: See TracBrowser for help on using the repository browser.