| 1 | ! Declarations specifiques au cas Toga |
|---|
| 2 | character*80 :: fich_toga |
|---|
| 3 | ! integer nlev_prof |
|---|
| 4 | ! parameter (nlev_prof = 41) |
|---|
| 5 | integer nlev_toga, nt_toga |
|---|
| 6 | parameter (nlev_toga=41, nt_toga=480) |
|---|
| 7 | integer year_ini_toga, day_ini_toga, mth_ini_toga |
|---|
| 8 | real day_ju_ini_toga ! Julian day of toga coare first day |
|---|
| 9 | parameter (year_ini_toga=1992) |
|---|
| 10 | parameter (mth_ini_toga=11) |
|---|
| 11 | parameter (day_ini_toga=1) ! 1erNov1992 |
|---|
| 12 | real dt_toga |
|---|
| 13 | parameter (dt_toga=6.*3600.) |
|---|
| 14 | !! |
|---|
| 15 | integer year_print, month_print, day_print |
|---|
| 16 | real sec_print |
|---|
| 17 | !! |
|---|
| 18 | integer nn |
|---|
| 19 | integer it_toga1, it_toga2 |
|---|
| 20 | real time_toga1,time_toga2 |
|---|
| 21 | |
|---|
| 22 | real ts_toga(nt_toga) |
|---|
| 23 | real plev_toga(nlev_toga,nt_toga),w_toga(nlev_toga,nt_toga) |
|---|
| 24 | real t_toga(nlev_toga,nt_toga),q_toga(nlev_toga,nt_toga) |
|---|
| 25 | real u_toga(nlev_toga,nt_toga),v_toga(nlev_toga,nt_toga) |
|---|
| 26 | real ht_toga(nlev_toga,nt_toga),vt_toga(nlev_toga,nt_toga) |
|---|
| 27 | real hq_toga(nlev_toga,nt_toga),vq_toga(nlev_toga,nt_toga) |
|---|
| 28 | |
|---|
| 29 | real ts_prof |
|---|
| 30 | real plev_prof(nlev_toga),w_prof(nlev_toga) |
|---|
| 31 | real t_prof(nlev_toga),q_prof(nlev_toga) |
|---|
| 32 | real u_prof(nlev_toga),v_prof(nlev_toga) |
|---|
| 33 | real ht_prof(nlev_toga),vt_prof(nlev_toga) |
|---|
| 34 | real hq_prof(nlev_toga),vq_prof(nlev_toga) |
|---|
| 35 | |
|---|
| 36 | real plev_mod(llm),w_mod(llm), t_mod(llm),q_mod(llm) |
|---|
| 37 | real u_mod(llm),v_mod(llm), ht_mod(llm),vt_mod(llm) |
|---|
| 38 | real hq_mod(llm),vq_mod(llm),qv_mod(llm),ql_mod(llm),qt_mod(llm) |
|---|
| 39 | |
|---|
| 40 | real ts_cur |
|---|
| 41 | common /sst_forcing/ts_cur ! also in read_tsurf1d.F |
|---|
| 42 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
|---|
| 43 | ! Declarations specifiques au cas RICO |
|---|
| 44 | character*80 :: fich_rico |
|---|
| 45 | integer nlev_rico |
|---|
| 46 | |
|---|
| 47 | parameter (nlev_rico=81) |
|---|
| 48 | real ts_rico,ps_rico |
|---|
| 49 | real w_rico(llm) |
|---|
| 50 | real t_rico(llm),q_rico(llm) |
|---|
| 51 | real u_rico(llm),v_rico(llm) |
|---|
| 52 | real dth_rico(llm) |
|---|
| 53 | real dqh_rico(llm) |
|---|
| 54 | real du_age(llm),dv_age(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 | ! Declarations specifiques au cas GCSSold |
|---|
| 85 | character*80 :: fich_gcssold_ctl |
|---|
| 86 | character*80 :: fich_gcssold_dat |
|---|
| 87 | real ht_gcssold(llm),hq_gcssold(llm),hw_gcssold(llm) |
|---|
| 88 | real hu_gcssold(llm) |
|---|
| 89 | real hv_gcssold(llm) |
|---|
| 90 | real hthturb_gcssold(llm) |
|---|
| 91 | real hqturb_gcssold(llm) |
|---|
| 92 | real Ts_gcssold |
|---|
| 93 | real dtime_frcg |
|---|
| 94 | cAl1 logical :: imp_fcg_gcssold |
|---|
| 95 | c logical :: ts_fcg_gcssold |
|---|
| 96 | c logical :: Tp_fcg_gcssold |
|---|
| 97 | logical :: Turb_fcg_gcssold |
|---|
| 98 | common /turb_forcing/ dtime_frcg, |
|---|
| 99 | $ Turb_fcg_gcssold, hthturb_gcssold, hqturb_gcssold |
|---|
| 100 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
|---|
| 101 | ! Declarations specifiques au cas Arm_cu |
|---|
| 102 | character*80 :: fich_armcu |
|---|
| 103 | |
|---|
| 104 | |
|---|
| 105 | integer nlev_armcu, nt_armcu |
|---|
| 106 | parameter (nlev_armcu=40, nt_armcu=31) |
|---|
| 107 | integer year_ini_armcu, day_ini_armcu, mth_ini_armcu |
|---|
| 108 | real heure_ini_armcu |
|---|
| 109 | real day_ju_ini_armcu ! Julian day of armcu case first day |
|---|
| 110 | parameter (year_ini_armcu=1997) |
|---|
| 111 | parameter (mth_ini_armcu=6) |
|---|
| 112 | parameter (day_ini_armcu=21) ! 172 = 21 juin 1997 |
|---|
| 113 | parameter (heure_ini_armcu=41400) ! 11:30 en secondes |
|---|
| 114 | real dt_armcu |
|---|
| 115 | parameter (dt_armcu=1.*1800.) ! forcages donnes ttes les demi-heures par ifa_armcu.txt |
|---|
| 116 | real sens_armcu(nt_armcu),flat_armcu(nt_armcu) |
|---|
| 117 | real adv_theta_armcu(nt_armcu),rad_theta_armcu(nt_armcu) |
|---|
| 118 | real adv_qt_armcu(nt_armcu) |
|---|
| 119 | real theta_mod(llm),rv_mod(llm),play_mod(llm) |
|---|
| 120 | real d_t_dyn_ls(llm),d_q_dyn_ls(llm) |
|---|
| 121 | ! profc comme "profil armcu" |
|---|
| 122 | real h_profc,play_profc,t_profc,th_profc,plev_profc |
|---|
| 123 | real u_profc,v_profc,qv_profc,rv_profc |
|---|
| 124 | |
|---|
| 125 | ! forcages interpoles dans le temps |
|---|
| 126 | real adv_theta_prof,rad_theta_prof,adv_qt_prof |
|---|
| 127 | real sens_prof,flat_prof,fact |
|---|
| 128 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
|---|
| 129 | |
|---|