1 | |
---|
2 | ! Declarations specifiques au cas Toga |
---|
3 | character*80 :: fich_toga |
---|
4 | ! integer nlev_prof |
---|
5 | ! parameter (nlev_prof = 41) |
---|
6 | integer nlev_toga, nt_toga |
---|
7 | parameter (nlev_toga=41, nt_toga=480) |
---|
8 | integer year_ini_toga, day_ini_toga, mth_ini_toga |
---|
9 | real day_ju_ini_toga ! Julian day of toga coare first day |
---|
10 | parameter (year_ini_toga=1992) |
---|
11 | parameter (mth_ini_toga=11) |
---|
12 | parameter (day_ini_toga=1) ! 1erNov1992 |
---|
13 | real dt_toga |
---|
14 | parameter (dt_toga=6.*3600.) |
---|
15 | !! |
---|
16 | integer year_print, month_print, day_print |
---|
17 | real sec_print |
---|
18 | !! |
---|
19 | real ts_toga(nt_toga) |
---|
20 | real plev_toga(nlev_toga,nt_toga),w_toga(nlev_toga,nt_toga) |
---|
21 | real t_toga(nlev_toga,nt_toga),q_toga(nlev_toga,nt_toga) |
---|
22 | real u_toga(nlev_toga,nt_toga),v_toga(nlev_toga,nt_toga) |
---|
23 | real ht_toga(nlev_toga,nt_toga),vt_toga(nlev_toga,nt_toga) |
---|
24 | real hq_toga(nlev_toga,nt_toga),vq_toga(nlev_toga,nt_toga) |
---|
25 | |
---|
26 | real ts_prof |
---|
27 | real plev_prof(nlev_toga),w_prof(nlev_toga) |
---|
28 | real t_prof(nlev_toga),q_prof(nlev_toga) |
---|
29 | real u_prof(nlev_toga),v_prof(nlev_toga) |
---|
30 | real ht_prof(nlev_toga),vt_prof(nlev_toga) |
---|
31 | real hq_prof(nlev_toga),vq_prof(nlev_toga) |
---|
32 | |
---|
33 | real w_mod(llm), t_mod(llm),q_mod(llm) |
---|
34 | real u_mod(llm),v_mod(llm), ht_mod(llm),vt_mod(llm),ug_mod(llm),vg_mod(llm) |
---|
35 | real hq_mod(llm),vq_mod(llm),qv_mod(llm),ql_mod(llm),qt_mod(llm) |
---|
36 | real th_mod(llm) |
---|
37 | |
---|
38 | !real ts_cur |
---|
39 | !common /sst_forcing/ts_cur ! also in read_tsurf1d.F |
---|
40 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
41 | ! Declarations specifiques au cas RICO |
---|
42 | character*80 :: fich_rico |
---|
43 | integer nlev_rico |
---|
44 | |
---|
45 | parameter (nlev_rico=81) |
---|
46 | real ts_rico,ps_rico |
---|
47 | real w_rico(llm) |
---|
48 | real t_rico(llm),q_rico(llm) |
---|
49 | real u_rico(llm),v_rico(llm) |
---|
50 | real dth_rico(llm) |
---|
51 | real dqh_rico(llm) |
---|
52 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
53 | ! Declarations specifiques au cas TWPice |
---|
54 | character*80 :: fich_twpice |
---|
55 | integer nlev_twpi, nt_twpi |
---|
56 | parameter (nlev_twpi=40, nt_twpi=215) |
---|
57 | integer year_ini_twpi, day_ini_twpi, mth_ini_twpi |
---|
58 | real heure_ini_twpi |
---|
59 | real day_ju_ini_twpi ! Julian day of twpice first day |
---|
60 | parameter (year_ini_twpi=2006) |
---|
61 | parameter (mth_ini_twpi=1) |
---|
62 | parameter (day_ini_twpi=17) ! 17 = 17Jan2006 |
---|
63 | parameter (heure_ini_twpi=10800.) !3h en secondes |
---|
64 | real dt_twpi |
---|
65 | parameter (dt_twpi=3.*3600.) |
---|
66 | |
---|
67 | real ts_twpi(nt_twpi) |
---|
68 | real plev_twpi(nlev_twpi,nt_twpi),w_twpi(nlev_twpi,nt_twpi) |
---|
69 | real t_twpi(nlev_twpi,nt_twpi),q_twpi(nlev_twpi,nt_twpi) |
---|
70 | real u_twpi(nlev_twpi,nt_twpi),v_twpi(nlev_twpi,nt_twpi) |
---|
71 | real ht_twpi(nlev_twpi,nt_twpi),vt_twpi(nlev_twpi,nt_twpi) |
---|
72 | real hq_twpi(nlev_twpi,nt_twpi),vq_twpi(nlev_twpi,nt_twpi) |
---|
73 | |
---|
74 | real ts_proftwp |
---|
75 | real plev_proftwp(nlev_twpi),w_proftwp(nlev_twpi) |
---|
76 | real t_proftwp(nlev_twpi),q_proftwp(nlev_twpi) |
---|
77 | real u_proftwp(nlev_twpi),v_proftwp(nlev_twpi) |
---|
78 | real ht_proftwp(nlev_twpi),vt_proftwp(nlev_twpi) |
---|
79 | real hq_proftwp(nlev_twpi),vq_proftwp(nlev_twpi) |
---|
80 | |
---|
81 | |
---|
82 | |
---|
83 | !Declarations specifiques au cas FIRE |
---|
84 | character*80 :: fich_fire |
---|
85 | integer nlev_fire, nt_fire |
---|
86 | parameter (nlev_fire=120, nt_fire=1) |
---|
87 | integer year_ini_fire, day_ini_fire, mth_ini_fire |
---|
88 | real heure_ini_fire |
---|
89 | parameter (year_ini_fire=1987) |
---|
90 | parameter (mth_ini_fire=7) |
---|
91 | parameter (day_ini_fire=14) ! 14 = 14Juil1987 |
---|
92 | parameter (heure_ini_fire=0.) !0h en secondes |
---|
93 | |
---|
94 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
95 | !Declarations specifiques au cas GABLS4 (MPL 20141023) |
---|
96 | character*80 :: fich_gabls4 |
---|
97 | integer nlev_gabls4, nt_gabls4, nsol_gabls4 |
---|
98 | parameter (nlev_gabls4=90, nt_gabls4=37, nsol_gabls4=19) |
---|
99 | integer year_ini_gabls4, day_ini_gabls4, mth_ini_gabls4 |
---|
100 | real heure_ini_gabls4 |
---|
101 | real day_ju_ini_gabls4 ! Julian day of gabls4 first day |
---|
102 | parameter (year_ini_gabls4=2009) |
---|
103 | parameter (mth_ini_gabls4=12) |
---|
104 | parameter (day_ini_gabls4=11) ! 11 = 11 decembre 2009 |
---|
105 | parameter (heure_ini_gabls4=0.) !0UTC en secondes |
---|
106 | real dt_gabls4 |
---|
107 | parameter (dt_gabls4=3600.) ! 1 forcage ttes les heures |
---|
108 | |
---|
109 | !profils initiaux: |
---|
110 | real plev_gabls4(nlev_gabls4) |
---|
111 | real zz_gabls4(nlev_gabls4) |
---|
112 | real th_gabls4(nlev_gabls4),t_gabls4(nlev_gabls4),qv_gabls4(nlev_gabls4) |
---|
113 | real u_gabls4(nlev_gabls4), v_gabls4(nlev_gabls4) |
---|
114 | real depth_sn_gabls4(nsol_gabls4),tsnow_gabls4(nsol_gabls4),snow_dens_gabls4(nsol_gabls4) |
---|
115 | real t_gabi(nlev_gabls4),qv_gabi(nlev_gabls4) |
---|
116 | real u_gabi(nlev_gabls4), v_gabi(nlev_gabls4),ug_gabi(nlev_gabls4), vg_gabi(nlev_gabls4) |
---|
117 | real ht_gabi(nlev_gabls4),hq_gabi(nlev_gabls4),poub(nlev_gabls4) |
---|
118 | |
---|
119 | !forcings |
---|
120 | real ht_gabls4(nlev_gabls4,nt_gabls4),hq_gabls4(nlev_gabls4,nt_gabls4) |
---|
121 | real ug_gabls4(nlev_gabls4,nt_gabls4),vg_gabls4(nlev_gabls4,nt_gabls4) |
---|
122 | real tg_gabls4(nt_gabls4) |
---|
123 | real ht_profg(nlev_gabls4),hq_profg(nlev_gabls4) |
---|
124 | real ug_profg(nlev_gabls4),vg_profg(nlev_gabls4) |
---|
125 | real tg_profg |
---|
126 | |
---|
127 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
128 | |
---|
129 | !Declarations specifiques au cas DICE (MPL 02072013) |
---|
130 | character*80 :: fich_dice |
---|
131 | integer nlev_dice, nt_dice |
---|
132 | parameter (nlev_dice=70, nt_dice=145) |
---|
133 | integer year_ini_dice, day_ini_dice, mth_ini_dice |
---|
134 | real heure_ini_dice |
---|
135 | real day_ju_ini_dice ! Julian day of dice first day |
---|
136 | parameter (year_ini_dice=1999) |
---|
137 | parameter (mth_ini_dice=10) |
---|
138 | parameter (day_ini_dice=23) ! 23 = 23 october 1999 |
---|
139 | parameter (heure_ini_dice=68400.) !19UTC en secondes |
---|
140 | real dt_dice |
---|
141 | parameter (dt_dice=0.5*3600.) ! 1 forcage ttes les demi-heures |
---|
142 | |
---|
143 | !profils initiaux: |
---|
144 | real plev_dice(nlev_dice) |
---|
145 | |
---|
146 | real zz_dice(nlev_dice) |
---|
147 | real t_dice(nlev_dice),qv_dice(nlev_dice) |
---|
148 | real u_dice(nlev_dice), v_dice(nlev_dice),o3_dice(nlev_dice) |
---|
149 | real ht_dice(nlev_dice,nt_dice) |
---|
150 | real hq_dice(nlev_dice,nt_dice), hu_dice(nlev_dice,nt_dice) |
---|
151 | real hv_dice(nlev_dice,nt_dice) |
---|
152 | real w_dice(nlev_dice,nt_dice),omega_dice(nlev_dice,nt_dice) |
---|
153 | real o3_mod(llm),hu_mod(llm),hv_mod(llm) |
---|
154 | real t_dicei(nlev_dice),qv_dicei(nlev_dice) |
---|
155 | real u_dicei(nlev_dice), v_dicei(nlev_dice),o3_dicei(nlev_dice) |
---|
156 | real ht_dicei(nlev_dice) |
---|
157 | real hq_dicei(nlev_dice), hu_dicei(nlev_dice) |
---|
158 | real hv_dicei(nlev_dice) |
---|
159 | real w_dicei(nlev_dice),omega_dicei(nlev_dice) |
---|
160 | |
---|
161 | |
---|
162 | !forcings |
---|
163 | real shf_dice(nt_dice),lhf_dice(nt_dice) |
---|
164 | real lwup_dice(nt_dice),swup_dice(nt_dice) |
---|
165 | real tg_dice(nt_dice),ustar_dice(nt_dice),psurf_dice(nt_dice) |
---|
166 | real ug_dice(nt_dice),vg_dice(nt_dice) |
---|
167 | |
---|
168 | real shf_prof,lhf_prof,lwup_prof,swup_prof,tg_prof |
---|
169 | real ustar_prof,psurf_prof,cdrag |
---|
170 | real ht_profd(nlev_dice),hq_profd(nlev_dice),hu_profd(nlev_dice) |
---|
171 | real hv_profd(nlev_dice),w_profd(nlev_dice) |
---|
172 | real omega_profd(nlev_dice),ug_profd,vg_profd |
---|
173 | |
---|
174 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
175 | ! Declarations specifiques au cas GCSSold |
---|
176 | character*80 :: fich_gcssold_ctl |
---|
177 | character*80 :: fich_gcssold_dat |
---|
178 | real ht_gcssold(llm),hq_gcssold(llm),hw_gcssold(llm) |
---|
179 | real hu_gcssold(llm) |
---|
180 | real hv_gcssold(llm) |
---|
181 | real hthturb_gcssold(llm) |
---|
182 | real hqturb_gcssold(llm) |
---|
183 | real Ts_gcssold |
---|
184 | real dtime_frcg |
---|
185 | logical :: Turb_fcg_gcssold |
---|
186 | |
---|
187 | common /turb_forcing/ & |
---|
188 | & dtime_frcg,hthturb_gcssold, hqturb_gcssold,Turb_fcg_gcssold |
---|
189 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
190 | ! Declarations specifiques au cas Arm_cu |
---|
191 | character*80 :: fich_armcu |
---|
192 | |
---|
193 | |
---|
194 | integer nlev_armcu, nt_armcu |
---|
195 | parameter (nlev_armcu=40, nt_armcu=31) |
---|
196 | integer year_ini_armcu, day_ini_armcu, mth_ini_armcu |
---|
197 | real heure_ini_armcu |
---|
198 | real day_ju_ini_armcu ! Julian day of armcu case first day |
---|
199 | parameter (year_ini_armcu=1997) |
---|
200 | parameter (mth_ini_armcu=6) |
---|
201 | parameter (day_ini_armcu=21) ! 172 = 21 juin 1997 |
---|
202 | parameter (heure_ini_armcu=41400) ! 11:30 en secondes |
---|
203 | real dt_armcu |
---|
204 | parameter (dt_armcu=1.*1800.) ! forcages donnes ttes les demi-heures par ifa_armcu.txt |
---|
205 | real sens_armcu(nt_armcu),flat_armcu(nt_armcu) |
---|
206 | real adv_theta_armcu(nt_armcu),rad_theta_armcu(nt_armcu) |
---|
207 | real adv_qt_armcu(nt_armcu) |
---|
208 | real theta_mod(llm),rv_mod(llm),play_mod(llm) |
---|
209 | ! profc comme "profil armcu" |
---|
210 | |
---|
211 | ! forcages interpoles dans le temps |
---|
212 | real adv_theta_prof,rad_theta_prof,adv_qt_prof |
---|
213 | real sens_prof,flat_prof,fact |
---|
214 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
215 | ! declarations specifiques au cas Sandu |
---|
216 | character*80 :: fich_sandu |
---|
217 | ! integer nlev_prof |
---|
218 | ! parameter (nlev_prof = 41) |
---|
219 | integer nlev_sandu, nt_sandu |
---|
220 | parameter (nlev_sandu=87, nt_sandu=13) |
---|
221 | integer year_ini_sandu, day_ini_sandu, mth_ini_sandu |
---|
222 | real day_ju_ini_sandu ! Julian day of sandu case first day |
---|
223 | parameter (year_ini_sandu=2006) |
---|
224 | parameter (mth_ini_sandu=7) |
---|
225 | parameter (day_ini_sandu=15) ! 196 = 15 juillet 2006 |
---|
226 | real dt_sandu, tau_sandu |
---|
227 | logical :: trouve_700=.true. |
---|
228 | parameter (dt_sandu=6.*3600.) ! forcages donnes ttes les 6 heures par ifa_sandu.txt |
---|
229 | parameter (tau_sandu=30000*3600.) ! temps de relaxation u,v,thetal,qt vers profil init et au dessus 700hPa |
---|
230 | !! |
---|
231 | real ts_sandu(nt_sandu) |
---|
232 | ! profs comme "profil sandu" |
---|
233 | real plev_profs(nlev_sandu) |
---|
234 | real t_profs(nlev_sandu),thl_profs(nlev_sandu) |
---|
235 | real q_profs(nlev_sandu) |
---|
236 | real u_profs(nlev_sandu),v_profs(nlev_sandu),w_profs(nlev_sandu) |
---|
237 | real omega_profs(nlev_sandu),o3mmr_profs(nlev_sandu) |
---|
238 | |
---|
239 | real, dimension(llm) :: relax_u,relax_v,relax_thl |
---|
240 | real, dimension(llm,2) :: relax_q |
---|
241 | |
---|
242 | real thl_mod(llm),omega_mod(llm),o3mmr_mod(llm),tke_mod(llm) |
---|
243 | !vertical advection computation |
---|
244 | real d_t_z(llm),d_th_z(llm), d_q_z(llm) |
---|
245 | real d_t_dyn_z(llm),d_th_dyn_z(llm), d_q_dyn_z(llm) |
---|
246 | real d_u_z(llm),d_v_z(llm) |
---|
247 | real d_u_dyn(llm),d_v_dyn(llm) |
---|
248 | real d_u_dyn_z(llm),d_v_dyn_z(llm) |
---|
249 | real d_u_adv(llm),d_v_adv(llm) |
---|
250 | real zz(llm) |
---|
251 | real zfact |
---|
252 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
253 | ! Declarations specifiques au cas Astex |
---|
254 | character*80 :: fich_astex |
---|
255 | integer nlev_astex, nt_astex |
---|
256 | parameter (nlev_astex=34, nt_astex=49) |
---|
257 | integer year_ini_astex, day_ini_astex, mth_ini_astex |
---|
258 | real day_ju_ini_astex ! Julian day of astex case first day |
---|
259 | parameter (year_ini_astex=1992) |
---|
260 | parameter (mth_ini_astex=6) |
---|
261 | parameter (day_ini_astex=13) ! 165 = 13 juin 1992 |
---|
262 | real dt_astex |
---|
263 | parameter (dt_astex=3600.) ! forcages donnes ttes les heures par ifa_astex.txt |
---|
264 | real ts_astex(nt_astex),div_astex(nt_astex),ug_astex(nt_astex) |
---|
265 | real vg_astex(nt_astex),ufa_astex(nt_astex),vfa_astex(nt_astex) |
---|
266 | real div_prof,ug_prof,vg_prof,ufa_prof,vfa_prof |
---|
267 | ! profa comme "profil astex" |
---|
268 | real plev_profa(nlev_astex) |
---|
269 | real t_profa(nlev_astex),thl_profa(nlev_astex) |
---|
270 | real qv_profa(nlev_astex),ql_profa(nlev_astex) |
---|
271 | real qt_profa(nlev_astex),o3mmr_profa(nlev_astex) |
---|
272 | real u_profa(nlev_astex),v_profa(nlev_astex),w_profa(nlev_astex) |
---|
273 | real tke_profa(nlev_astex) |
---|
274 | |
---|
275 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
276 | !Declarations specifiques au cas standard |
---|
277 | |
---|
278 | real w_mod_cas(llm), t_mod_cas(llm),q_mod_cas(llm) |
---|
279 | real theta_mod_cas(llm),thl_mod_cas(llm),thv_mod_cas(llm) |
---|
280 | real qv_mod_cas(llm),ql_mod_cas(llm),qi_mod_cas(llm) |
---|
281 | real ug_mod_cas(llm),vg_mod_cas(llm) |
---|
282 | real u_mod_cas(llm),v_mod_cas(llm) |
---|
283 | real omega_mod_cas(llm) |
---|
284 | real ht_mod_cas(llm),vt_mod_cas(llm),dt_mod_cas(llm),dtrad_mod_cas(llm) |
---|
285 | real hth_mod_cas(llm),vth_mod_cas(llm),dth_mod_cas(llm) |
---|
286 | real hq_mod_cas(llm),vq_mod_cas(llm),dq_mod_cas(llm) |
---|
287 | real hu_mod_cas(llm),vu_mod_cas(llm),du_mod_cas(llm) |
---|
288 | real hv_mod_cas(llm),vv_mod_cas(llm),dv_mod_cas(llm) |
---|
289 | integer day_ini_cas |
---|
290 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
291 | |
---|
292 | |
---|