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 | |
---|