LMDZ
1D_decl_cases.h
Go to the documentation of this file.
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)
8  parameter (nlev_toga=41, nt_toga=480)
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 !!
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 
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)
36  real hq_mod(llm),vq_mod(llm),qv_mod(llm),ql_mod(llm),qt_mod(llm)
37  real th_mod(llm)
38 
39  real ts_cur
40  common /sst_forcing/ts_cur ! also in read_tsurf1d.F
41 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
42 ! Declarations specifiques au cas RICO
43  character*80 :: fich_rico
44  integer nlev_rico
45 
46  parameter (nlev_rico=81)
47  real ts_rico,ps_rico
48  real w_rico(llm)
49  real t_rico(llm),q_rico(llm)
50  real u_rico(llm),v_rico(llm)
51  real dth_rico(llm)
52  real dqh_rico(llm)
53 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
54 ! Declarations specifiques au cas TWPice
55  character*80 :: fich_twpice
56  integer nlev_twpi, nt_twpi
57  parameter (nlev_twpi=40, nt_twpi=215)
58  integer year_ini_twpi, day_ini_twpi, mth_ini_twpi
59  real heure_ini_twpi
60  real day_ju_ini_twpi ! Julian day of twpice first day
61  parameter (year_ini_twpi=2006)
62  parameter (mth_ini_twpi=1)
63  parameter (day_ini_twpi=17) ! 17 = 17Jan2006
64  parameter (heure_ini_twpi=10800.) !3h en secondes
65  real dt_twpi
66  parameter (dt_twpi=3.*3600.)
67 
68  real ts_twpi(nt_twpi)
69  real plev_twpi(nlev_twpi,nt_twpi),w_twpi(nlev_twpi,nt_twpi)
70  real t_twpi(nlev_twpi,nt_twpi),q_twpi(nlev_twpi,nt_twpi)
71  real u_twpi(nlev_twpi,nt_twpi),v_twpi(nlev_twpi,nt_twpi)
72  real ht_twpi(nlev_twpi,nt_twpi),vt_twpi(nlev_twpi,nt_twpi)
73  real hq_twpi(nlev_twpi,nt_twpi),vq_twpi(nlev_twpi,nt_twpi)
74 
75  real ts_proftwp
76  real plev_proftwp(nlev_twpi),w_proftwp(nlev_twpi)
77  real t_proftwp(nlev_twpi),q_proftwp(nlev_twpi)
78  real u_proftwp(nlev_twpi),v_proftwp(nlev_twpi)
79  real ht_proftwp(nlev_twpi),vt_proftwp(nlev_twpi)
80  real hq_proftwp(nlev_twpi),vq_proftwp(nlev_twpi)
81 
82 
83 
84 !Declarations specifiques au cas FIRE
85  character*80 :: fich_fire
86  integer nlev_fire, nt_fire
87  parameter (nlev_fire=120, nt_fire=1)
88  integer year_ini_fire, day_ini_fire, mth_ini_fire
89  real heure_ini_fire
90  parameter (year_ini_fire=1987)
91  parameter (mth_ini_fire=7)
92  parameter (day_ini_fire=14) ! 14 = 14Juil1987
93  parameter (heure_ini_fire=0.) !0h en secondes
94 
95 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
96 !Declarations specifiques au cas DICE (MPL 02072013)
97  character*80 :: fich_dice
99  parameter (nlev_dice=70, nt_dice=145)
100  integer year_ini_dice, day_ini_dice, mth_ini_dice
101  real heure_ini_dice
102  real day_ju_ini_dice ! Julian day of dice first day
103  parameter (year_ini_dice=1999)
104  parameter (mth_ini_dice=10)
105  parameter (day_ini_dice=23) ! 23 = 23 october 1999
106  parameter (heure_ini_dice=68400.) !19UTC en secondes
107  real dt_dice
108  parameter (dt_dice=0.5*3600.) ! 1 forcage ttes les demi-heures
109 
110 !profils initiaux:
111  real plev_dice(nlev_dice)
112 
113  real zz_dice(nlev_dice)
114  real th_dice(nlev_dice),qv_dice(nlev_dice)
115  real u_dice(nlev_dice), v_dice(nlev_dice),o3_dice(nlev_dice)
116  real ht_dice(nlev_dice,nt_dice)
117  real hq_dice(nlev_dice,nt_dice), hu_dice(nlev_dice,nt_dice)
118  real hv_dice(nlev_dice,nt_dice)
119  real w_dice(nlev_dice,nt_dice),omega_dice(nlev_dice,nt_dice)
120  real o3_mod(llm),hu_mod(llm),hv_mod(llm)
121  real th_dicei(nlev_dice),qv_dicei(nlev_dice)
122  real u_dicei(nlev_dice), v_dicei(nlev_dice),o3_dicei(nlev_dice)
123  real ht_dicei(nlev_dice)
124  real hq_dicei(nlev_dice), hu_dicei(nlev_dice)
125  real hv_dicei(nlev_dice)
126  real w_dicei(nlev_dice),omega_dicei(nlev_dice)
127 
128 
129 !forcings
130  real shf_dice(nt_dice),lhf_dice(nt_dice)
131  real lwup_dice(nt_dice),swup_dice(nt_dice)
132  real tg_dice(nt_dice),ustar_dice(nt_dice),psurf_dice(nt_dice)
133  real ug_dice(nt_dice),vg_dice(nt_dice)
134 
137  real ht_profd(nlev_dice),hq_profd(nlev_dice),hu_profd(nlev_dice)
138  real hv_profd(nlev_dice),w_profd(nlev_dice)
139  real omega_profd(nlev_dice),ug_profd,vg_profd
140 
141 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
142 ! Declarations specifiques au cas GCSSold
143  character*80 :: fich_gcssold_ctl
144  character*80 :: fich_gcssold_dat
145  real ht_gcssold(llm),hq_gcssold(llm),hw_gcssold(llm)
146  real hu_gcssold(llm)
147  real hv_gcssold(llm)
148  real hthturb_gcssold(llm)
149  real hqturb_gcssold(llm)
151  real dtime_frcg
152  logical :: Turb_fcg_gcssold
153 
154  common /turb_forcing/ &
155  & dtime_frcg,hthturb_gcssold, hqturb_gcssold,Turb_fcg_gcssold
156 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
157 ! Declarations specifiques au cas Arm_cu
158  character*80 :: fich_armcu
159 
160 
161  integer nlev_armcu, nt_armcu
162  parameter (nlev_armcu=40, nt_armcu=31)
163  integer year_ini_armcu, day_ini_armcu, mth_ini_armcu
164  real heure_ini_armcu
165  real day_ju_ini_armcu ! Julian day of armcu case first day
166  parameter (year_ini_armcu=1997)
167  parameter (mth_ini_armcu=6)
168  parameter (day_ini_armcu=21) ! 172 = 21 juin 1997
169  parameter (heure_ini_armcu=41400) ! 11:30 en secondes
170  real dt_armcu
171  parameter (dt_armcu=1.*1800.) ! forcages donnes ttes les demi-heures par ifa_armcu.txt
172  real sens_armcu(nt_armcu),flat_armcu(nt_armcu)
173  real adv_theta_armcu(nt_armcu),rad_theta_armcu(nt_armcu)
174  real adv_qt_armcu(nt_armcu)
175  real theta_mod(llm),rv_mod(llm),play_mod(llm)
176 ! profc comme "profil armcu"
177 
178 ! forcages interpoles dans le temps
179  real adv_theta_prof,rad_theta_prof,adv_qt_prof
180  real sens_prof,flat_prof,fact
181 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
182 ! declarations specifiques au cas Sandu
183  character*80 :: fich_sandu
184 ! integer nlev_prof
185 ! parameter (nlev_prof = 41)
186  integer nlev_sandu, nt_sandu
187  parameter (nlev_sandu=87, nt_sandu=13)
188  integer year_ini_sandu, day_ini_sandu, mth_ini_sandu
189  real day_ju_ini_sandu ! Julian day of sandu case first day
190  parameter (year_ini_sandu=2006)
191  parameter (mth_ini_sandu=7)
192  parameter (day_ini_sandu=15) ! 196 = 15 juillet 2006
193  real dt_sandu, tau_sandu
194  logical :: trouve_700=.true.
195  parameter (dt_sandu=6.*3600.) ! forcages donnes ttes les 6 heures par ifa_sandu.txt
196  parameter (tau_sandu=3600.) ! temps de relaxation u,v,thetal,qt vers profil init et au dessus 700hPa
197 !!
198  real ts_sandu(nt_sandu)
199 ! profs comme "profil sandu"
200  real plev_profs(nlev_sandu)
201  real t_profs(nlev_sandu),thl_profs(nlev_sandu)
202  real q_profs(nlev_sandu)
203  real u_profs(nlev_sandu),v_profs(nlev_sandu),w_profs(nlev_sandu)
204  real omega_profs(nlev_sandu),o3mmr_profs(nlev_sandu)
205 
206  real, dimension(llm) :: relax_u,relax_v,relax_thl
207  real, dimension(llm,2) :: relax_q
208 
209  real thl_mod(llm),omega_mod(llm),o3mmr_mod(llm),tke_mod(llm)
210 !vertical advection computation
211  real d_t_z(llm), d_q_z(llm)
212  real d_t_dyn_z(llm), d_q_dyn_z(llm)
213  real d_u_z(llm),d_v_z(llm)
214  real d_u_dyn(llm),d_v_dyn(llm)
215  real d_u_dyn_z(llm),d_v_dyn_z(llm)
216  real d_u_adv(llm),d_v_adv(llm)
217  real zz(llm)
218  real zfact
219 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
220 ! Declarations specifiques au cas Astex
221  character*80 :: fich_astex
222  integer nlev_astex, nt_astex
223  parameter (nlev_astex=34, nt_astex=49)
224  integer year_ini_astex, day_ini_astex, mth_ini_astex
225  real day_ju_ini_astex ! Julian day of astex case first day
226  parameter (year_ini_astex=1992)
227  parameter (mth_ini_astex=6)
228  parameter (day_ini_astex=13) ! 165 = 13 juin 1992
229  real dt_astex
230  parameter (dt_astex=3600.) ! forcages donnes ttes les heures par ifa_astex.txt
231  real ts_astex(nt_astex),div_astex(nt_astex),ug_astex(nt_astex)
232  real vg_astex(nt_astex),ufa_astex(nt_astex),vfa_astex(nt_astex)
233  real div_prof,ug_prof,vg_prof,ufa_prof,vfa_prof
234 ! profa comme "profil astex"
235  real plev_profa(nlev_astex)
236  real t_profa(nlev_astex),thl_profa(nlev_astex)
237  real qv_profa(nlev_astex),ql_profa(nlev_astex)
238  real qt_profa(nlev_astex),o3mmr_profa(nlev_astex)
239  real u_profa(nlev_astex),v_profa(nlev_astex),w_profa(nlev_astex)
240  real tke_profa(nlev_astex)
241 
242 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
243 !Declarations specifiques au cas standard
244 
245  real w_mod_cas(llm), t_mod_cas(llm),q_mod_cas(llm)
246  real ug_mod_cas(llm),vg_mod_cas(llm)
247  real u_mod_cas(llm),v_mod_cas(llm)
248  real ht_mod_cas(llm),vt_mod_cas(llm),dt_mod_cas(llm),dtrad_mod_cas(llm)
249  real hq_mod_cas(llm),vq_mod_cas(llm),dq_mod_cas(llm)
250  real hu_mod_cas(llm),vu_mod_cas(llm),du_mod_cas(llm)
251  real hv_mod_cas(llm),vv_mod_cas(llm),dv_mod_cas(llm)
252  integer day_ini_cas
253 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
254 
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice w_dice
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice lhf_prof
!$Id Ts_gcssold
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd hv_profd
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd ht_profd
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice ustar_dice
!$Id hv_gcssold
subroutine read_tsurf1d(knon, sst_out)
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga u_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice ug_dice
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga plev_prof
do llm!au dessus on relaxe vers profil init!on fait l hypothese que dans ce il n y a plus d eau liq au dessus!donc la relaxation en thetal et qt devient relaxation en tempe et qv l dq1 relax relax_q(l, 1)
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga w_toga
real, parameter demi
Definition: VARphy.F90:12
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vt_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga v_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice ht_dice
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd w_profd omega_profd!do llm!print llm l llm
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga q_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof ht_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd w_profd omega_profd!do llm!print llm l llm!print llm l omega_mod
do llm!au dessus on relaxe vers profil init!on fait l hypothese que dans ce il n y a plus d eau liq au dessus!donc la relaxation en thetal et qt devient relaxation en tempe et qv l dq1 relax l dq2 relax l dt relax relax_thl(l)
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga t_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice lwup_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice hu_dice
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd hq_profd
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga t_toga
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice lhf_dice
do llm!au dessus on relaxe vers profil init!on fait l hypothese que dans ce il n y a plus d eau liq au dessus!donc la relaxation en thetal et qt devient relaxation en tempe et qv tau_sandu
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof ug_profd
!$Id && ht_gcssold
!$Id && hu_gcssold
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref day_ju_ini_toga
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof ustar_prof
!Declarations specifiques au cas Toga character *::fich_toga!integer nlev_prof nt_toga day_ini_toga
Definition: 1D_decl_cases.h:9
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice swup_prof
!$Id && hthturb_gcssold
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm &&& day
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof psurf_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref nt_toga
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof hq_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref year_ini_toga
!$Id hq_gcssold
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice psurf_dice
!$Id fich_gcssold_dat
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd w_profd omega_profd!do llm!print llm l omega_profd
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm u(l)
!$Id fich_gcssold_ctl
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga ts_prof
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL ulevSTD CALL &zphi philevSTD CALL &zx_rh rhlevSTD!DO klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon du jour ou toutes les read_climoz CALL true
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga ts_toga
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice shf_dice
!$Id nlev_fire
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vt_toga
subroutine cdrag(knon, nsrf, speed, t1, q1, zgeop1, psol, tsurf, qsurf, z0m, z0h, pcfm, pcfh, zri, pref)
Definition: cdrag.F90:8
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice shf_prof
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice tg_dice
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref year_ini_dice
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga ht_toga
nsplit_thermals!nrlmd le iflag_clos_bl tau_trig_deep real::s_trig!fin nrlmd le fact_thermals_ed_dz integer
Definition: thermcell.h:9
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga plev_toga
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga q_toga
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga nlev_toga
!Declarations specifiques au cas Toga character *::fich_toga!integer nlev_prof nt_toga parameter(nlev_toga=41, nt_toga=480) integer year_ini_toga
do llm!au dessus on relaxe vers profil init!on fait l hypothese que dans ce il n y a plus d eau liq au dessus!donc la relaxation en thetal et qt devient relaxation en tempe et qv relax_u(l)
do llm!au dessus on relaxe vers profil init!on fait l hypothese que dans ce il n y a plus d eau liq au dessus!donc la relaxation en thetal et qt devient relaxation en tempe et qv u_mod(l)!if(l.ge.llm700) then relax_q(l
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice hv_dice
do llm!au dessus de
!$Id hqturb_gcssold
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga hq_toga
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice nlev_dice
real, dimension(:,:), allocatable, save d_u_dyn
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice hq_dice
real, dimension(:,:), allocatable, save d_v_dyn
!IM for NMC files!real nfiles!real nfiles!real nfiles!real nlevSTD3 real nlevSTD3 real nlevSTD3 real nlevSTD3 real nlevSTD8 real nlevSTD8 real nlevSTD8 real nlevSTD8 real
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice lwup_dice
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga v_toga
do llm!au dessus on relaxe vers profil init!on fait l hypothese que dans ce cas
!$Id hw_gcssold
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref nt_dice
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref day_ju_ini_dice
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd hu_profd