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