1 | ! |
---|
2 | ! $Id: phys_state_var_mod.F90 3781 2020-10-26 09:33:06Z musat $ |
---|
3 | ! |
---|
4 | MODULE phys_state_var_mod |
---|
5 | ! Variables sauvegardees pour le startphy.nc |
---|
6 | !====================================================================== |
---|
7 | ! |
---|
8 | ! |
---|
9 | !====================================================================== |
---|
10 | ! Declaration des variables |
---|
11 | USE dimphy |
---|
12 | USE netcdf, only: nf90_fill_real |
---|
13 | INTEGER, PARAMETER :: nlevSTD=17 |
---|
14 | INTEGER, PARAMETER :: nlevSTD8=8 |
---|
15 | INTEGER, PARAMETER :: nlevSTD3=3 |
---|
16 | INTEGER, PARAMETER :: nout=3 |
---|
17 | INTEGER, PARAMETER :: napisccp=1 |
---|
18 | INTEGER, SAVE :: radpas ! radiation is called every "radpas" step |
---|
19 | INTEGER, SAVE :: cvpas ! convection is called every "cvpas" step |
---|
20 | INTEGER, SAVE :: cvpas_0 = 1 ! reference value for cvpas |
---|
21 | INTEGER, SAVE :: wkpas ! wake scheme is called every "wkpas" step |
---|
22 | REAL, PARAMETER :: missing_val_nf90=nf90_fill_real |
---|
23 | !$OMP THREADPRIVATE(radpas) |
---|
24 | !$OMP THREADPRIVATE(cvpas) |
---|
25 | !$OMP THREADPRIVATE(cvpas_0) |
---|
26 | !$OMP THREADPRIVATE(wkpas) |
---|
27 | REAL, SAVE :: phys_tstep=0, solaire_etat0 |
---|
28 | !$OMP THREADPRIVATE(phys_tstep, solaire_etat0) |
---|
29 | |
---|
30 | REAL, ALLOCATABLE, SAVE :: pctsrf(:,:) |
---|
31 | !$OMP THREADPRIVATE(pctsrf) |
---|
32 | REAL, ALLOCATABLE, SAVE :: ftsol(:,:) |
---|
33 | !$OMP THREADPRIVATE(ftsol) |
---|
34 | REAL,ALLOCATABLE,SAVE :: qsol(:),fevap(:,:),z0m(:,:),z0h(:,:),agesno(:,:) |
---|
35 | !$OMP THREADPRIVATE(qsol,fevap,z0m,z0h,agesno) |
---|
36 | !FC drag des arbres |
---|
37 | REAL, ALLOCATABLE, SAVE :: treedrg(:,:,:) |
---|
38 | !$OMP THREADPRIVATE(treedrg) |
---|
39 | |
---|
40 | ! character(len=6), SAVE :: ocean |
---|
41 | !!!!!!$OMP THREADPRIVATE(ocean) |
---|
42 | ! logical, SAVE :: ok_veget |
---|
43 | !!!!!!$OMP THREADPRIVATE(ok_veget) |
---|
44 | REAL, ALLOCATABLE, SAVE :: falb1(:,:), falb2(:,:) |
---|
45 | !$OMP THREADPRIVATE(falb1, falb2) |
---|
46 | |
---|
47 | !albedo SB >>> |
---|
48 | REAL, ALLOCATABLE, SAVE :: falb_dif(:,:,:), falb_dir(:,:,:) |
---|
49 | REAL, ALLOCATABLE, SAVE :: chl_con(:) |
---|
50 | !$OMP THREADPRIVATE(falb_dir,falb_dif,chl_con) |
---|
51 | !albedo SB <<< |
---|
52 | |
---|
53 | |
---|
54 | REAL, ALLOCATABLE, SAVE :: rain_fall(:), snow_fall(:) |
---|
55 | !$OMP THREADPRIVATE( rain_fall, snow_fall) |
---|
56 | REAL, ALLOCATABLE, SAVE :: solsw(:), solswfdiff(:), sollw(:) |
---|
57 | !$OMP THREADPRIVATE(solsw, solswfdiff, sollw) |
---|
58 | REAL, ALLOCATABLE, SAVE :: radsol(:) |
---|
59 | !$OMP THREADPRIVATE(radsol) |
---|
60 | REAL, ALLOCATABLE, SAVE :: swradcorr(:) |
---|
61 | !$OMP THREADPRIVATE(swradcorr) |
---|
62 | |
---|
63 | !clesphy0 param physiq |
---|
64 | ! |
---|
65 | ! Parametres de l'Orographie a l'Echelle Sous-Maille (OESM): |
---|
66 | ! |
---|
67 | REAL, ALLOCATABLE, SAVE :: zmea(:), zstd(:), zsig(:), zgam(:) |
---|
68 | !$OMP THREADPRIVATE(zmea, zstd, zsig, zgam) |
---|
69 | REAL, ALLOCATABLE, SAVE :: zthe(:), zpic(:), zval(:) |
---|
70 | !$OMP THREADPRIVATE(zthe, zpic, zval) |
---|
71 | ! REAL tabcntr0(100) |
---|
72 | REAL, ALLOCATABLE, SAVE :: rugoro(:) |
---|
73 | !$OMP THREADPRIVATE(rugoro) |
---|
74 | REAL, ALLOCATABLE, SAVE :: t_ancien(:,:), q_ancien(:,:) |
---|
75 | !$OMP THREADPRIVATE(t_ancien, q_ancien) |
---|
76 | REAL, ALLOCATABLE, SAVE :: ql_ancien(:,:), qs_ancien(:,:) |
---|
77 | !$OMP THREADPRIVATE(ql_ancien, qs_ancien) |
---|
78 | REAL, ALLOCATABLE, SAVE :: prw_ancien(:), prlw_ancien(:), prsw_ancien(:) |
---|
79 | !$OMP THREADPRIVATE(prw_ancien, prlw_ancien, prsw_ancien) |
---|
80 | REAL, ALLOCATABLE, SAVE :: u_ancien(:,:), v_ancien(:,:) |
---|
81 | !$OMP THREADPRIVATE(u_ancien, v_ancien) |
---|
82 | !!! RomP >>> |
---|
83 | REAL, ALLOCATABLE, SAVE :: tr_ancien(:,:,:) |
---|
84 | !$OMP THREADPRIVATE(tr_ancien) |
---|
85 | !!! RomP <<< |
---|
86 | LOGICAL, SAVE :: ancien_ok |
---|
87 | !$OMP THREADPRIVATE(ancien_ok) |
---|
88 | REAL, ALLOCATABLE, SAVE :: clwcon(:,:),rnebcon(:,:) |
---|
89 | !$OMP THREADPRIVATE(clwcon,rnebcon) |
---|
90 | REAL, ALLOCATABLE, SAVE :: qtc_cv(:,:),sigt_cv(:,:) |
---|
91 | !$OMP THREADPRIVATE(qtc_cv,sigt_cv) |
---|
92 | REAL, ALLOCATABLE, SAVE :: ratqs(:,:) |
---|
93 | !$OMP THREADPRIVATE(ratqs) |
---|
94 | REAL, ALLOCATABLE, SAVE :: pbl_tke(:,:,:) ! turb kinetic energy |
---|
95 | REAL, ALLOCATABLE, SAVE :: coefh(:,:,:) ! Kz enthalpie |
---|
96 | REAL, ALLOCATABLE, SAVE :: coefm(:,:,:) ! Kz momentum |
---|
97 | !$OMP THREADPRIVATE(pbl_tke, coefh,coefm) |
---|
98 | !nrlmd< |
---|
99 | REAL, ALLOCATABLE, SAVE :: delta_tsurf(:,:) ! Surface temperature difference inside-outside cold pool |
---|
100 | !$OMP THREADPRIVATE(delta_tsurf) |
---|
101 | !>nrlmd |
---|
102 | REAL, ALLOCATABLE, SAVE :: zmax0(:), f0(:) ! |
---|
103 | !$OMP THREADPRIVATE(zmax0,f0) |
---|
104 | REAL, ALLOCATABLE, SAVE :: sig1(:,:), w01(:,:) |
---|
105 | !$OMP THREADPRIVATE(sig1,w01) |
---|
106 | REAL, ALLOCATABLE, SAVE :: entr_therm(:,:), fm_therm(:,:) |
---|
107 | !$OMP THREADPRIVATE(entr_therm,fm_therm) |
---|
108 | REAL, ALLOCATABLE, SAVE :: detr_therm(:,:) |
---|
109 | !$OMP THREADPRIVATE(detr_therm) |
---|
110 | !IM 150408 |
---|
111 | ! pour phsystoke avec thermiques |
---|
112 | REAL,ALLOCATABLE,SAVE :: clwcon0th(:,:),rnebcon0th(:,:) |
---|
113 | !$OMP THREADPRIVATE(clwcon0th,rnebcon0th) |
---|
114 | ! radiation outputs |
---|
115 | REAL,ALLOCATABLE,SAVE :: swdnc0(:,:), swdn0(:,:), swdn(:,:) |
---|
116 | !$OMP THREADPRIVATE(swdnc0,swdn0,swdn) |
---|
117 | REAL,ALLOCATABLE,SAVE :: swupc0(:,:), swup0(:,:), swup(:,:) |
---|
118 | !$OMP THREADPRIVATE(swupc0, swup0,swup) |
---|
119 | REAL,ALLOCATABLE,SAVE :: SWdn200clr(:), SWdn200(:) |
---|
120 | !$OMP THREADPRIVATE(SWdn200clr,SWdn200) |
---|
121 | REAL,ALLOCATABLE,SAVE :: SWup200clr(:), SWup200(:) |
---|
122 | !$OMP THREADPRIVATE(SWup200clr,SWup200) |
---|
123 | REAL,ALLOCATABLE,SAVE :: lwdnc0(:,:), lwdn0(:,:), lwdn(:,:) |
---|
124 | !$OMP THREADPRIVATE(lwdnc0,lwdn0,lwdn) |
---|
125 | REAL,ALLOCATABLE,SAVE :: lwupc0(:,:), lwup0(:,:), lwup(:,:) |
---|
126 | !$OMP THREADPRIVATE(lwupc0,lwup0,lwup) |
---|
127 | REAL,ALLOCATABLE,SAVE :: LWdn200clr(:), LWdn200(:) |
---|
128 | !$OMP THREADPRIVATE(LWdn200clr,LWdn200) |
---|
129 | REAL,ALLOCATABLE,SAVE :: LWup200clr(:), LWup200(:) |
---|
130 | !$OMP THREADPRIVATE(LWup200clr,LWup200) |
---|
131 | REAL,ALLOCATABLE,SAVE :: LWdnTOA(:), LWdnTOAclr(:) |
---|
132 | !$OMP THREADPRIVATE(LWdnTOA,LWdnTOAclr) |
---|
133 | ! pressure level |
---|
134 | REAL,ALLOCATABLE,SAVE :: tsumSTD(:,:,:) |
---|
135 | !$OMP THREADPRIVATE(tsumSTD) |
---|
136 | REAL,ALLOCATABLE,SAVE :: usumSTD(:,:,:), vsumSTD(:,:,:) |
---|
137 | !$OMP THREADPRIVATE(usumSTD,vsumSTD) |
---|
138 | REAL,ALLOCATABLE,SAVE :: wsumSTD(:,:,:), phisumSTD(:,:,:) |
---|
139 | !$OMP THREADPRIVATE(wsumSTD,phisumSTD) |
---|
140 | REAL,ALLOCATABLE,SAVE :: qsumSTD(:,:,:), rhsumSTD(:,:,:) |
---|
141 | !$OMP THREADPRIVATE(qsumSTD,rhsumSTD) |
---|
142 | REAL,ALLOCATABLE,SAVE :: tnondef(:,:,:) |
---|
143 | !$OMP THREADPRIVATE(tnondef) |
---|
144 | REAL,ALLOCATABLE,SAVE :: uvsumSTD(:,:,:) |
---|
145 | !$OMP THREADPRIVATE(uvsumSTD) |
---|
146 | REAL,ALLOCATABLE,SAVE :: vqsumSTD(:,:,:) |
---|
147 | !$OMP THREADPRIVATE(vqsumSTD) |
---|
148 | REAL,ALLOCATABLE,SAVE :: vTsumSTD(:,:,:) |
---|
149 | !$OMP THREADPRIVATE(vTsumSTD) |
---|
150 | REAL,ALLOCATABLE,SAVE :: wqsumSTD(:,:,:) |
---|
151 | !$OMP THREADPRIVATE(wqsumSTD) |
---|
152 | REAL,ALLOCATABLE,SAVE :: vphisumSTD(:,:,:) |
---|
153 | !$OMP THREADPRIVATE(vphisumSTD) |
---|
154 | REAL,ALLOCATABLE,SAVE :: wTsumSTD(:,:,:) |
---|
155 | !$OMP THREADPRIVATE(wTsumSTD) |
---|
156 | REAL,ALLOCATABLE,SAVE :: u2sumSTD(:,:,:) |
---|
157 | !$OMP THREADPRIVATE(u2sumSTD) |
---|
158 | REAL,ALLOCATABLE,SAVE :: v2sumSTD(:,:,:) |
---|
159 | !$OMP THREADPRIVATE(v2sumSTD) |
---|
160 | REAL,ALLOCATABLE,SAVE :: T2sumSTD(:,:,:) |
---|
161 | !$OMP THREADPRIVATE(T2sumSTD) |
---|
162 | REAL,ALLOCATABLE,SAVE :: O3sumSTD(:,:,:), O3daysumSTD(:,:,:) |
---|
163 | !$OMP THREADPRIVATE(O3sumSTD,O3daysumSTD) |
---|
164 | !IM begin |
---|
165 | REAL,ALLOCATABLE,SAVE :: wlevSTD(:,:), ulevSTD(:,:), vlevSTD(:,:) |
---|
166 | !$OMP THREADPRIVATE(wlevSTD,ulevSTD,vlevSTD) |
---|
167 | REAL,ALLOCATABLE,SAVE :: tlevSTD(:,:), qlevSTD(:,:), rhlevSTD(:,:) |
---|
168 | !$OMP THREADPRIVATE(tlevSTD,qlevSTD,rhlevSTD) |
---|
169 | REAL,ALLOCATABLE,SAVE :: philevSTD(:,:) |
---|
170 | !$OMP THREADPRIVATE(philevSTD) |
---|
171 | REAL,ALLOCATABLE,SAVE :: uvSTD(:,:) |
---|
172 | !$OMP THREADPRIVATE(uvSTD) |
---|
173 | REAL,ALLOCATABLE,SAVE :: vqSTD(:,:) |
---|
174 | !$OMP THREADPRIVATE(vqSTD) |
---|
175 | REAL,ALLOCATABLE,SAVE :: vTSTD(:,:) |
---|
176 | !$OMP THREADPRIVATE(vTSTD) |
---|
177 | REAL,ALLOCATABLE,SAVE :: wqSTD(:,:) |
---|
178 | !$OMP THREADPRIVATE(wqSTD) |
---|
179 | REAL,ALLOCATABLE,SAVE :: vphiSTD(:,:) |
---|
180 | !$OMP THREADPRIVATE(vphiSTD) |
---|
181 | REAL,ALLOCATABLE,SAVE :: wTSTD(:,:) |
---|
182 | !$OMP THREADPRIVATE(wTSTD) |
---|
183 | REAL,ALLOCATABLE,SAVE :: u2STD(:,:) |
---|
184 | !$OMP THREADPRIVATE(u2STD) |
---|
185 | REAL,ALLOCATABLE,SAVE :: v2STD(:,:) |
---|
186 | !$OMP THREADPRIVATE(v2STD) |
---|
187 | REAL,ALLOCATABLE,SAVE :: T2STD(:,:) |
---|
188 | !$OMP THREADPRIVATE(T2STD) |
---|
189 | REAL,ALLOCATABLE,SAVE :: O3STD(:,:), O3daySTD(:,:) |
---|
190 | !$OMP THREADPRIVATE(O3STD,O3daySTD) |
---|
191 | !IM end |
---|
192 | INTEGER,ALLOCATABLE,SAVE :: seed_old(:,:) |
---|
193 | !$OMP THREADPRIVATE(seed_old) |
---|
194 | REAL,ALLOCATABLE,SAVE :: zuthe(:),zvthe(:) |
---|
195 | !$OMP THREADPRIVATE(zuthe,zvthe) |
---|
196 | REAL,ALLOCATABLE,SAVE :: alb_neig(:) |
---|
197 | !$OMP THREADPRIVATE(alb_neig) |
---|
198 | !cloud base mass flux |
---|
199 | REAL,ALLOCATABLE,SAVE :: ema_cbmf(:) |
---|
200 | !$OMP THREADPRIVATE(ema_cbmf) |
---|
201 | !cloud base pressure & cloud top pressure |
---|
202 | REAL,ALLOCATABLE,SAVE :: ema_pcb(:), ema_pct(:) |
---|
203 | !$OMP THREADPRIVATE(ema_pcb,ema_pct) |
---|
204 | REAL,ALLOCATABLE,SAVE :: Mipsh(:,:) ! mass flux shed from adiab. ascents |
---|
205 | !$OMP THREADPRIVATE(Mipsh) |
---|
206 | REAL,ALLOCATABLE,SAVE :: Ma(:,:) ! undilute upward mass flux |
---|
207 | !$OMP THREADPRIVATE(Ma) |
---|
208 | REAL,ALLOCATABLE,SAVE :: qcondc(:,:) ! in-cld water content from convect |
---|
209 | !$OMP THREADPRIVATE(qcondc) |
---|
210 | REAL,ALLOCATABLE,SAVE :: wd(:) ! sb |
---|
211 | !$OMP THREADPRIVATE(wd) |
---|
212 | REAL,ALLOCATABLE,SAVE :: sigd(:) |
---|
213 | !$OMP THREADPRIVATE(sigd) |
---|
214 | ! |
---|
215 | REAL,ALLOCATABLE,SAVE :: cin(:) |
---|
216 | !$OMP THREADPRIVATE(cin) |
---|
217 | ! ftd : convective heating due to unsaturated downdraughts |
---|
218 | REAL,ALLOCATABLE,SAVE :: ftd(:,:) |
---|
219 | !$OMP THREADPRIVATE(ftd) |
---|
220 | ! fqd : convective moistening due to unsaturated downdraughts |
---|
221 | REAL,ALLOCATABLE,SAVE :: fqd(:,:) |
---|
222 | !$OMP THREADPRIVATE(fqd) |
---|
223 | !34EK |
---|
224 | ! -- Variables de controle de ALE et ALP |
---|
225 | !ALE : Energie disponible pour soulevement : utilisee par la |
---|
226 | ! convection d'Emanuel pour le declenchement et la regulation |
---|
227 | REAL,ALLOCATABLE,SAVE :: ALE(:) |
---|
228 | !$OMP THREADPRIVATE(ALE) |
---|
229 | !ALP : Puissance disponible pour soulevement |
---|
230 | REAL,ALLOCATABLE,SAVE :: ALP(:) |
---|
231 | !$OMP THREADPRIVATE(ALP) |
---|
232 | ! |
---|
233 | ! nouvelles variables pour le couplage convection-couche limite |
---|
234 | REAL,ALLOCATABLE,SAVE :: Ale_bl(:) |
---|
235 | !$OMP THREADPRIVATE(Ale_bl) |
---|
236 | REAL,ALLOCATABLE,SAVE :: Alp_bl(:) |
---|
237 | !$OMP THREADPRIVATE(Alp_bl) |
---|
238 | INTEGER,ALLOCATABLE,SAVE :: lalim_conv(:) |
---|
239 | !$OMP THREADPRIVATE(lalim_conv) |
---|
240 | REAL,ALLOCATABLE,SAVE :: wght_th(:,:) |
---|
241 | !$OMP THREADPRIVATE(wght_th) |
---|
242 | REAL,ALLOCATABLE,SAVE :: ale_wake(:) |
---|
243 | !$OMP THREADPRIVATE(ale_wake) |
---|
244 | REAL,ALLOCATABLE,SAVE :: ale_bl_stat(:) |
---|
245 | !$OMP THREADPRIVATE(ale_bl_stat) |
---|
246 | ! |
---|
247 | ! variables de la wake |
---|
248 | ! wake_deltat : ecart de temperature avec la zone non perturbee |
---|
249 | ! wake_deltaq : ecart d'humidite avec la zone non perturbee |
---|
250 | ! wake_s : fraction surfacique occupee par la poche froide |
---|
251 | ! awake_dens : number of active wakes per unit area |
---|
252 | ! wake_dens : number of wakes per unit area |
---|
253 | ! wake_occ : occurence of wakes (= 1 if wakes occur, =0 otherwise) |
---|
254 | ! wake_Cstar : vitesse d'etalement de la poche |
---|
255 | ! wake_pe : wake potential energy - WAPE |
---|
256 | ! wake_fip : Gust Front Impinging power - ALP |
---|
257 | REAL,ALLOCATABLE,SAVE :: wake_deltat(:,:) |
---|
258 | !$OMP THREADPRIVATE(wake_deltat) |
---|
259 | REAL,ALLOCATABLE,SAVE :: wake_deltaq(:,:) |
---|
260 | !$OMP THREADPRIVATE(wake_deltaq) |
---|
261 | REAL,ALLOCATABLE,SAVE :: wake_s(:) |
---|
262 | !$OMP THREADPRIVATE(wake_s) |
---|
263 | REAL,ALLOCATABLE,SAVE :: awake_dens(:), wake_dens(:) |
---|
264 | !$OMP THREADPRIVATE(awake_dens, wake_dens) |
---|
265 | REAL,ALLOCATABLE,SAVE :: wake_Cstar(:) |
---|
266 | !$OMP THREADPRIVATE(wake_Cstar) |
---|
267 | REAL,ALLOCATABLE,SAVE :: wake_pe(:) |
---|
268 | !$OMP THREADPRIVATE(wake_pe) |
---|
269 | REAL,ALLOCATABLE,SAVE :: wake_fip(:) |
---|
270 | !$OMP THREADPRIVATE(wake_fip) |
---|
271 | ! |
---|
272 | !jyg< |
---|
273 | ! variables related to the spitting of the PBL between wake and |
---|
274 | ! off-wake regions. |
---|
275 | ! wake_delta_pbl_TKE : difference TKE_w - TKE_x |
---|
276 | REAL,ALLOCATABLE,SAVE :: wake_delta_pbl_TKE(:,:,:) |
---|
277 | !$OMP THREADPRIVATE(wake_delta_pbl_TKE) |
---|
278 | !>jyg |
---|
279 | ! |
---|
280 | ! pfrac_impa : Produits des coefs lessivage impaction |
---|
281 | ! pfrac_nucl : Produits des coefs lessivage nucleation |
---|
282 | ! pfrac_1nucl: Produits des coefs lessi nucl (alpha = 1) |
---|
283 | REAL,ALLOCATABLE,SAVE :: pfrac_impa(:,:), pfrac_nucl(:,:) |
---|
284 | !$OMP THREADPRIVATE(pfrac_impa,pfrac_nucl) |
---|
285 | REAL,ALLOCATABLE,SAVE :: pfrac_1nucl(:,:) |
---|
286 | !$OMP THREADPRIVATE(pfrac_1nucl) |
---|
287 | ! |
---|
288 | REAL,ALLOCATABLE,SAVE :: total_rain(:), nday_rain(:) |
---|
289 | !$OMP THREADPRIVATE(total_rain,nday_rain) |
---|
290 | REAL,ALLOCATABLE,SAVE :: paire_ter(:) |
---|
291 | !$OMP THREADPRIVATE(paire_ter) |
---|
292 | ! albsol1: albedo du sol total pour SW visible |
---|
293 | ! albsol2: albedo du sol total pour SW proche IR |
---|
294 | REAL,ALLOCATABLE,SAVE :: albsol1(:), albsol2(:) |
---|
295 | !$OMP THREADPRIVATE(albsol1,albsol2) |
---|
296 | |
---|
297 | !albedo SB >>> |
---|
298 | REAL,ALLOCATABLE,SAVE :: albsol_dif(:,:),albsol_dir(:,:) |
---|
299 | !$OMP THREADPRIVATE(albsol_dif,albsol_dir) |
---|
300 | !albedo SB <<< |
---|
301 | |
---|
302 | |
---|
303 | REAL, ALLOCATABLE, SAVE:: wo(:, :, :) |
---|
304 | ! column-density of ozone in a layer, in kilo-Dobsons |
---|
305 | ! Third dimension has size 1 or 2. |
---|
306 | ! "wo(:, :, 1)" is for the average day-night field, |
---|
307 | ! "wo(:, :, 2)" is for daylight time. |
---|
308 | !$OMP THREADPRIVATE(wo) |
---|
309 | |
---|
310 | ! heat : chauffage solaire |
---|
311 | ! heat0: chauffage solaire ciel clair |
---|
312 | ! cool : refroidissement infrarouge |
---|
313 | ! cool0 : refroidissement infrarouge ciel clair |
---|
314 | ! sollwdown : downward LW flux at surface |
---|
315 | ! sollwdownclr : downward CS LW flux at surface |
---|
316 | ! toplwdown : downward CS LW flux at TOA |
---|
317 | ! toplwdownclr : downward CS LW flux at TOA |
---|
318 | ! heat_volc : chauffage solaire du au volcanisme |
---|
319 | ! cool_volc : refroidissement infrarouge du au volcanisme |
---|
320 | REAL,ALLOCATABLE,SAVE :: clwcon0(:,:),rnebcon0(:,:) |
---|
321 | !$OMP THREADPRIVATE(clwcon0,rnebcon0) |
---|
322 | REAL,ALLOCATABLE,SAVE :: heat(:,:) |
---|
323 | !$OMP THREADPRIVATE(heat) |
---|
324 | REAL,ALLOCATABLE,SAVE :: heat0(:,:) |
---|
325 | !$OMP THREADPRIVATE(heat0) |
---|
326 | REAL,ALLOCATABLE,SAVE :: cool(:,:) |
---|
327 | !$OMP THREADPRIVATE(cool) |
---|
328 | REAL,ALLOCATABLE,SAVE :: cool0(:,:) |
---|
329 | !$OMP THREADPRIVATE(cool0) |
---|
330 | REAL,ALLOCATABLE,SAVE :: heat_volc(:,:) |
---|
331 | !$OMP THREADPRIVATE(heat_volc) |
---|
332 | REAL,ALLOCATABLE,SAVE :: cool_volc(:,:) |
---|
333 | !$OMP THREADPRIVATE(cool_volc) |
---|
334 | REAL,ALLOCATABLE,SAVE :: topsw(:), toplw(:) |
---|
335 | !$OMP THREADPRIVATE(topsw,toplw) |
---|
336 | REAL,ALLOCATABLE,SAVE :: sollwdown(:) |
---|
337 | !$OMP THREADPRIVATE(sollwdown) |
---|
338 | REAL,ALLOCATABLE,SAVE :: gustiness(:) |
---|
339 | !$OMP THREADPRIVATE(gustiness) |
---|
340 | REAL,ALLOCATABLE,SAVE :: sollwdownclr(:) |
---|
341 | !$OMP THREADPRIVATE(sollwdownclr) |
---|
342 | REAL,ALLOCATABLE,SAVE :: toplwdown(:) |
---|
343 | !$OMP THREADPRIVATE(toplwdown) |
---|
344 | REAL,ALLOCATABLE,SAVE :: toplwdownclr(:) |
---|
345 | !$OMP THREADPRIVATE(toplwdownclr) |
---|
346 | REAL,ALLOCATABLE,SAVE :: topsw0(:),toplw0(:),solsw0(:),sollw0(:) |
---|
347 | !$OMP THREADPRIVATE(topsw0,toplw0,solsw0,sollw0) |
---|
348 | REAL,ALLOCATABLE,SAVE :: albpla(:) |
---|
349 | !$OMP THREADPRIVATE(albpla) |
---|
350 | |
---|
351 | !IM ajout variables CFMIP2/CMIP5 |
---|
352 | REAL,ALLOCATABLE,SAVE :: heatp(:,:), coolp(:,:) |
---|
353 | !$OMP THREADPRIVATE(heatp, coolp) |
---|
354 | REAL,ALLOCATABLE,SAVE :: heat0p(:,:), cool0p(:,:) |
---|
355 | !$OMP THREADPRIVATE(heat0p, cool0p) |
---|
356 | REAL,ALLOCATABLE,SAVE :: radsolp(:), topswp(:), toplwp(:) |
---|
357 | !$OMP THREADPRIVATE(radsolp, topswp, toplwp) |
---|
358 | REAL,ALLOCATABLE,SAVE :: albplap(:) |
---|
359 | !$OMP THREADPRIVATE(albplap) |
---|
360 | REAL,ALLOCATABLE,SAVE :: solswp(:), solswfdiffp(:), sollwp(:) |
---|
361 | !$OMP THREADPRIVATE(solswp, solswfdiffp, sollwp) |
---|
362 | REAL,ALLOCATABLE,SAVE :: sollwdownp(:) |
---|
363 | !$OMP THREADPRIVATE(sollwdownp) |
---|
364 | REAL,ALLOCATABLE,SAVE :: topsw0p(:),toplw0p(:) |
---|
365 | REAL,ALLOCATABLE,SAVE :: solsw0p(:),sollw0p(:) |
---|
366 | !$OMP THREADPRIVATE(topsw0p,toplw0p,solsw0p,sollw0p) |
---|
367 | REAL,ALLOCATABLE,SAVE :: lwdnc0p(:,:), lwdn0p(:,:), lwdnp(:,:) |
---|
368 | REAL,ALLOCATABLE,SAVE :: lwupc0p(:,:), lwup0p(:,:), lwupp(:,:) |
---|
369 | !$OMP THREADPRIVATE(lwdnc0p, lwdn0p, lwdnp, lwupc0p, lwup0p, lwupp) |
---|
370 | REAL,ALLOCATABLE,SAVE :: swdnc0p(:,:), swdn0p(:,:), swdnp(:,:) |
---|
371 | REAL,ALLOCATABLE,SAVE :: swupc0p(:,:), swup0p(:,:), swupp(:,:) |
---|
372 | !$OMP THREADPRIVATE(swdnc0p, swdn0p, swdnp, swupc0p, swup0p, swupp) |
---|
373 | |
---|
374 | ! pbase : cloud base pressure |
---|
375 | ! bbase : cloud base buoyancy |
---|
376 | REAL,ALLOCATABLE,SAVE :: cape(:) |
---|
377 | !$OMP THREADPRIVATE(cape) |
---|
378 | REAL,ALLOCATABLE,SAVE :: pbase(:) |
---|
379 | !$OMP THREADPRIVATE(pbase) |
---|
380 | REAL,ALLOCATABLE,SAVE :: bbase(:) |
---|
381 | !$OMP THREADPRIVATE(bbase) |
---|
382 | ! |
---|
383 | REAL,SAVE,ALLOCATABLE :: zqasc(:,:) |
---|
384 | !$OMP THREADPRIVATE( zqasc) |
---|
385 | INTEGER,ALLOCATABLE,SAVE :: ibas_con(:), itop_con(:) |
---|
386 | !$OMP THREADPRIVATE(ibas_con,itop_con) |
---|
387 | REAL,SAVE,ALLOCATABLE :: rain_con(:) |
---|
388 | !$OMP THREADPRIVATE(rain_con) |
---|
389 | REAL,SAVE,ALLOCATABLE :: snow_con(:) |
---|
390 | !$OMP THREADPRIVATE(snow_con) |
---|
391 | ! |
---|
392 | REAL,SAVE,ALLOCATABLE :: rlonPOS(:) |
---|
393 | !$OMP THREADPRIVATE(rlonPOS) |
---|
394 | REAL,SAVE,ALLOCATABLE :: newsst(:) |
---|
395 | !$OMP THREADPRIVATE(newsst) |
---|
396 | REAL,SAVE,ALLOCATABLE :: ustar(:,:),u10m(:,:), v10m(:,:),wstar(:,:) |
---|
397 | !$OMP THREADPRIVATE(ustar,u10m,v10m,wstar) |
---|
398 | ! |
---|
399 | ! ok_ade=T -ADE=topswad-topsw |
---|
400 | ! ok_aie=T -> |
---|
401 | ! ok_ade=T -AIE=topswai-topswad |
---|
402 | ! ok_ade=F -AIE=topswai-topsw |
---|
403 | ! |
---|
404 | !topswad, solswad : Aerosol direct effect |
---|
405 | REAL,SAVE,ALLOCATABLE :: topswad(:), solswad(:) |
---|
406 | !$OMP THREADPRIVATE(topswad,solswad) |
---|
407 | !topswai, solswai : Aerosol indirect effect |
---|
408 | REAL,SAVE,ALLOCATABLE :: topswai(:), solswai(:) |
---|
409 | !$OMP THREADPRIVATE(topswai,solswai) |
---|
410 | |
---|
411 | REAL,SAVE,ALLOCATABLE :: tau_aero(:,:,:,:), piz_aero(:,:,:,:), cg_aero(:,:,:,:) |
---|
412 | !$OMP THREADPRIVATE(tau_aero, piz_aero, cg_aero) |
---|
413 | REAL,SAVE,ALLOCATABLE :: tau_aero_sw_rrtm(:,:,:,:), piz_aero_sw_rrtm(:,:,:,:), cg_aero_sw_rrtm(:,:,:,:) |
---|
414 | !$OMP THREADPRIVATE(tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm) |
---|
415 | REAL,SAVE,ALLOCATABLE :: tau_aero_lw_rrtm(:,:,:,:), piz_aero_lw_rrtm(:,:,:,:), cg_aero_lw_rrtm(:,:,:,:) |
---|
416 | !$OMP THREADPRIVATE(tau_aero_lw_rrtm, piz_aero_lw_rrtm, cg_aero_lw_rrtm) |
---|
417 | REAL,SAVE,ALLOCATABLE :: ccm(:,:,:) |
---|
418 | !$OMP THREADPRIVATE(ccm) |
---|
419 | |
---|
420 | !!! nrlmd le 10/04/2012 |
---|
421 | REAL,SAVE,ALLOCATABLE :: ale_bl_trig(:) |
---|
422 | !$OMP THREADPRIVATE(ale_bl_trig) |
---|
423 | !!! fin nrlmd le 10/04/2012 |
---|
424 | |
---|
425 | REAL, ALLOCATABLE, SAVE:: du_gwd_rando(:, :), du_gwd_front(:, :) |
---|
426 | !$OMP THREADPRIVATE(du_gwd_rando, du_gwd_front) |
---|
427 | ! tendencies on wind due to gravity waves |
---|
428 | |
---|
429 | LOGICAL,SAVE :: is_initialized=.FALSE. |
---|
430 | !$OMP THREADPRIVATE(is_initialized) |
---|
431 | |
---|
432 | CONTAINS |
---|
433 | |
---|
434 | !====================================================================== |
---|
435 | SUBROUTINE phys_state_var_init(read_climoz) |
---|
436 | USE dimphy |
---|
437 | USE aero_mod |
---|
438 | USE infotrac_phy, ONLY : nbtr |
---|
439 | USE indice_sol_mod |
---|
440 | IMPLICIT NONE |
---|
441 | |
---|
442 | integer, intent(in):: read_climoz |
---|
443 | ! read ozone climatology |
---|
444 | ! Allowed values are 0, 1 and 2 |
---|
445 | ! 0: do not read an ozone climatology |
---|
446 | ! 1: read a single ozone climatology that will be used day and night |
---|
447 | ! 2: read two ozone climatologies, the average day and night |
---|
448 | ! climatology and the daylight climatology |
---|
449 | |
---|
450 | include "clesphys.h" |
---|
451 | |
---|
452 | print*, 'is_initialized', is_initialized |
---|
453 | IF (is_initialized) RETURN |
---|
454 | is_initialized=.TRUE. |
---|
455 | ALLOCATE(pctsrf(klon,nbsrf)) |
---|
456 | ALLOCATE(ftsol(klon,nbsrf)) |
---|
457 | ALLOCATE(qsol(klon),fevap(klon,nbsrf)) |
---|
458 | ALLOCATE(z0m(klon,nbsrf+1),z0h(klon,nbsrf+1),agesno(klon,nbsrf)) |
---|
459 | !FC |
---|
460 | ALLOCATE(treedrg(klon,klev,nbsrf)) |
---|
461 | ALLOCATE(falb1(klon,nbsrf)) |
---|
462 | ALLOCATE(falb2(klon,nbsrf)) |
---|
463 | !albedo SB >>> |
---|
464 | print*, 'allocate falb' |
---|
465 | ALLOCATE(falb_dir(klon,nsw,nbsrf),falb_dif(klon,nsw,nbsrf)) |
---|
466 | print*, 'allocate falb good', falb_dir(1,1,1) |
---|
467 | ALLOCATE(chl_con(klon)) |
---|
468 | !albedo SB <<< |
---|
469 | ALLOCATE(rain_fall(klon)) |
---|
470 | ALLOCATE(snow_fall(klon)) |
---|
471 | ALLOCATE(solsw(klon), solswfdiff(klon), sollw(klon)) |
---|
472 | sollw=0.0 |
---|
473 | ALLOCATE(radsol(klon)) |
---|
474 | ALLOCATE(swradcorr(klon)) |
---|
475 | ALLOCATE(zmea(klon), zstd(klon), zsig(klon), zgam(klon)) |
---|
476 | ALLOCATE(zthe(klon), zpic(klon), zval(klon)) |
---|
477 | |
---|
478 | ALLOCATE(rugoro(klon)) |
---|
479 | ALLOCATE(t_ancien(klon,klev), q_ancien(klon,klev)) |
---|
480 | ALLOCATE(ql_ancien(klon,klev), qs_ancien(klon,klev)) |
---|
481 | ALLOCATE(prw_ancien(klon), prlw_ancien(klon), prsw_ancien(klon)) |
---|
482 | ALLOCATE(u_ancien(klon,klev), v_ancien(klon,klev)) |
---|
483 | !!! Rom P >>> |
---|
484 | ALLOCATE(tr_ancien(klon,klev,nbtr)) |
---|
485 | !!! Rom P <<< |
---|
486 | ALLOCATE(clwcon(klon,klev),rnebcon(klon,klev)) |
---|
487 | ALLOCATE(qtc_cv(klon,klev),sigt_cv(klon,klev)) |
---|
488 | ALLOCATE(ratqs(klon,klev)) |
---|
489 | ALLOCATE(pbl_tke(klon,klev+1,nbsrf+1)) |
---|
490 | !nrlmd< |
---|
491 | ALLOCATE(delta_tsurf(klon,nbsrf)) |
---|
492 | !>nrlmd |
---|
493 | ALLOCATE(coefh(klon,klev+1,nbsrf+1)) |
---|
494 | ALLOCATE(coefm(klon,klev+1,nbsrf+1)) |
---|
495 | ALLOCATE(zmax0(klon), f0(klon)) |
---|
496 | ALLOCATE(sig1(klon,klev), w01(klon,klev)) |
---|
497 | ALLOCATE(entr_therm(klon,klev), fm_therm(klon,klev+1)) |
---|
498 | ALLOCATE(detr_therm(klon,klev)) |
---|
499 | ! pour phsystoke avec thermiques |
---|
500 | ALLOCATE(clwcon0th(klon,klev),rnebcon0th(klon,klev)) |
---|
501 | ! radiation outputs |
---|
502 | ALLOCATE(swdnc0(klon,klevp1), swdn0(klon,klevp1), swdn(klon,klevp1)) |
---|
503 | ALLOCATE(swupc0(klon,klevp1), swup0(klon,klevp1), swup(klon,klevp1)) |
---|
504 | ALLOCATE(lwdnc0(klon,klevp1), lwdn0(klon,klevp1), lwdn(klon,klevp1)) |
---|
505 | ALLOCATE(lwupc0(klon,klevp1), lwup0(klon,klevp1), lwup(klon,klevp1)) |
---|
506 | ALLOCATE(SWdn200clr(klon), SWdn200(klon)) |
---|
507 | ALLOCATE(SWup200clr(klon), SWup200(klon)) |
---|
508 | ALLOCATE(LWdn200clr(klon), LWdn200(klon)) |
---|
509 | ALLOCATE(LWup200clr(klon), LWup200(klon)) |
---|
510 | ALLOCATE(LWdnTOA(klon), LWdnTOAclr(klon)) |
---|
511 | ! pressure level |
---|
512 | ALLOCATE(tsumSTD(klon,nlevSTD,nout)) |
---|
513 | ALLOCATE(usumSTD(klon,nlevSTD,nout), vsumSTD(klon,nlevSTD,nout)) |
---|
514 | ALLOCATE(wsumSTD(klon,nlevSTD,nout), phisumSTD(klon,nlevSTD,nout)) |
---|
515 | ALLOCATE(qsumSTD(klon,nlevSTD,nout), rhsumSTD(klon,nlevSTD,nout)) |
---|
516 | ALLOCATE(tnondef(klon,nlevSTD,nout)) |
---|
517 | ALLOCATE(uvsumSTD(klon,nlevSTD,nout)) |
---|
518 | ALLOCATE(vqsumSTD(klon,nlevSTD,nout)) |
---|
519 | ALLOCATE(vTsumSTD(klon,nlevSTD,nout)) |
---|
520 | ALLOCATE(wqsumSTD(klon,nlevSTD,nout)) |
---|
521 | ALLOCATE(vphisumSTD(klon,nlevSTD,nout)) |
---|
522 | ALLOCATE(wTsumSTD(klon,nlevSTD,nout)) |
---|
523 | ALLOCATE(u2sumSTD(klon,nlevSTD,nout)) |
---|
524 | ALLOCATE(v2sumSTD(klon,nlevSTD,nout)) |
---|
525 | ALLOCATE(T2sumSTD(klon,nlevSTD,nout)) |
---|
526 | ALLOCATE(O3sumSTD(klon,nlevSTD,nout)) |
---|
527 | ALLOCATE(O3daysumSTD(klon,nlevSTD,nout)) |
---|
528 | !IM beg |
---|
529 | ALLOCATE(wlevSTD(klon,nlevSTD), ulevSTD(klon,nlevSTD), vlevSTD(klon,nlevSTD)) |
---|
530 | ALLOCATE(tlevSTD(klon,nlevSTD), qlevSTD(klon,nlevSTD), rhlevSTD(klon,nlevSTD)) |
---|
531 | ALLOCATE(philevSTD(klon,nlevSTD)) |
---|
532 | ALLOCATE(uvSTD(klon,nlevSTD),vqSTD(klon,nlevSTD)) |
---|
533 | ALLOCATE(vTSTD(klon,nlevSTD),wqSTD(klon,nlevSTD)) |
---|
534 | ALLOCATE(vphiSTD(klon,nlevSTD),wTSTD(klon,nlevSTD)) |
---|
535 | ALLOCATE(u2STD(klon,nlevSTD),v2STD(klon,nlevSTD)) |
---|
536 | ALLOCATE(T2STD(klon,nlevSTD)) |
---|
537 | ALLOCATE(O3STD(klon,nlevSTD)) |
---|
538 | ALLOCATE(O3daySTD(klon,nlevSTD)) |
---|
539 | !IM end |
---|
540 | ALLOCATE(seed_old(klon,napisccp)) |
---|
541 | ALLOCATE(zuthe(klon),zvthe(klon)) |
---|
542 | ALLOCATE(alb_neig(klon)) |
---|
543 | !cloud base mass flux |
---|
544 | ALLOCATE(ema_cbmf(klon)) |
---|
545 | !cloud base pressure & cloud top pressure |
---|
546 | ALLOCATE(ema_pcb(klon), ema_pct(klon)) |
---|
547 | ! |
---|
548 | ALLOCATE(Mipsh(klon,klev)) |
---|
549 | ALLOCATE(Ma(klon,klev)) |
---|
550 | ALLOCATE(qcondc(klon,klev)) |
---|
551 | ALLOCATE(wd(klon)) |
---|
552 | ALLOCATE(sigd(klon)) |
---|
553 | ALLOCATE(cin(klon), ALE(klon), ALP(klon)) |
---|
554 | ALLOCATE(ftd(klon,klev), fqd(klon,klev)) |
---|
555 | ALLOCATE(Ale_bl(klon)) |
---|
556 | ALLOCATE(ale_wake(klon)) |
---|
557 | ALLOCATE(ale_bl_stat(klon)) |
---|
558 | ale_bl_stat(:)=0 |
---|
559 | ALLOCATE(Alp_bl(klon)) |
---|
560 | ALLOCATE(lalim_conv(klon)) |
---|
561 | ALLOCATE(wght_th(klon,klev)) |
---|
562 | ALLOCATE(wake_deltat(klon,klev), wake_deltaq(klon,klev)) |
---|
563 | ALLOCATE(wake_s(klon), awake_dens(klon), wake_dens(klon)) |
---|
564 | awake_dens = 0. |
---|
565 | ALLOCATE(wake_Cstar(klon)) |
---|
566 | ALLOCATE(wake_pe(klon), wake_fip(klon)) |
---|
567 | !jyg< |
---|
568 | ALLOCATE(wake_delta_pbl_TKE(klon,klev+1,nbsrf+1)) |
---|
569 | !>jyg |
---|
570 | ALLOCATE(pfrac_impa(klon,klev), pfrac_nucl(klon,klev)) |
---|
571 | ALLOCATE(pfrac_1nucl(klon,klev)) |
---|
572 | ALLOCATE(total_rain(klon), nday_rain(klon)) |
---|
573 | ALLOCATE(paire_ter(klon)) |
---|
574 | ALLOCATE(albsol1(klon), albsol2(klon)) |
---|
575 | !albedo SB >>> |
---|
576 | ALLOCATE(albsol_dir(klon,nsw),albsol_dif(klon,nsw)) |
---|
577 | !albedo SB <<< |
---|
578 | |
---|
579 | if (read_climoz <= 1) then |
---|
580 | ALLOCATE(wo(klon,klev, 1)) |
---|
581 | else |
---|
582 | ! read_climoz == 2 |
---|
583 | ALLOCATE(wo(klon,klev, 2)) |
---|
584 | end if |
---|
585 | |
---|
586 | ALLOCATE(clwcon0(klon,klev),rnebcon0(klon,klev)) |
---|
587 | ALLOCATE(heat(klon,klev), heat0(klon,klev)) |
---|
588 | ALLOCATE(cool(klon,klev), cool0(klon,klev)) |
---|
589 | ALLOCATE(heat_volc(klon,klev), cool_volc(klon,klev)) |
---|
590 | ALLOCATE(topsw(klon), toplw(klon)) |
---|
591 | ALLOCATE(sollwdown(klon), sollwdownclr(klon)) |
---|
592 | sollwdown = 0. |
---|
593 | ALLOCATE(toplwdown(klon), toplwdownclr(klon)) |
---|
594 | ALLOCATE(topsw0(klon),toplw0(klon),solsw0(klon),sollw0(klon)) |
---|
595 | sollw0 = 0. |
---|
596 | ALLOCATE(albpla(klon)) |
---|
597 | !IM ajout variables CFMIP2/CMIP5 |
---|
598 | ALLOCATE(heatp(klon,klev), coolp(klon,klev)) |
---|
599 | ALLOCATE(heat0p(klon,klev), cool0p(klon,klev)) |
---|
600 | ALLOCATE(radsolp(klon), topswp(klon), toplwp(klon)) |
---|
601 | ALLOCATE(albplap(klon)) |
---|
602 | ALLOCATE(solswp(klon), solswfdiffp(klon), sollwp(klon)) |
---|
603 | ALLOCATE(gustiness(klon)) |
---|
604 | ALLOCATE(sollwdownp(klon)) |
---|
605 | ALLOCATE(topsw0p(klon),toplw0p(klon)) |
---|
606 | ALLOCATE(solsw0p(klon),sollw0p(klon)) |
---|
607 | ALLOCATE(lwdnc0p(klon,klevp1), lwdn0p(klon,klevp1), lwdnp(klon,klevp1)) |
---|
608 | ALLOCATE(lwupc0p(klon,klevp1), lwup0p(klon,klevp1), lwupp(klon,klevp1)) |
---|
609 | ALLOCATE(swdnc0p(klon,klevp1), swdn0p(klon,klevp1), swdnp(klon,klevp1)) |
---|
610 | ALLOCATE(swupc0p(klon,klevp1), swup0p(klon,klevp1), swupp(klon,klevp1)) |
---|
611 | |
---|
612 | ALLOCATE(cape(klon)) |
---|
613 | ALLOCATE(pbase(klon),bbase(klon)) |
---|
614 | ALLOCATE(zqasc(klon,klev)) |
---|
615 | ALLOCATE(ibas_con(klon), itop_con(klon)) |
---|
616 | ALLOCATE(rain_con(klon), snow_con(klon)) |
---|
617 | ALLOCATE(rlonPOS(klon)) |
---|
618 | ALLOCATE(newsst(klon)) |
---|
619 | ALLOCATE(ustar(klon,nbsrf),u10m(klon,nbsrf), v10m(klon,nbsrf),wstar(klon,nbsrf+1)) |
---|
620 | ALLOCATE(topswad(klon), solswad(klon)) |
---|
621 | ALLOCATE(topswai(klon), solswai(klon)) |
---|
622 | ALLOCATE(tau_aero(klon,klev,naero_grp,nbands),piz_aero(klon,klev,naero_grp,nbands),cg_aero(klon,klev,naero_grp,nbands)) |
---|
623 | ALLOCATE(tau_aero_sw_rrtm(klon,klev,2,nbands_sw_rrtm),piz_aero_sw_rrtm(klon,klev,2,nbands_sw_rrtm)) |
---|
624 | ALLOCATE(cg_aero_sw_rrtm(klon,klev,2,nbands_sw_rrtm)) |
---|
625 | ALLOCATE(tau_aero_lw_rrtm(klon,klev,2,nbands_lw_rrtm),piz_aero_lw_rrtm(klon,klev,2,nbands_lw_rrtm)) |
---|
626 | ALLOCATE(cg_aero_lw_rrtm(klon,klev,2,nbands_lw_rrtm)) |
---|
627 | ALLOCATE(ccm(klon,klev,nbands)) |
---|
628 | |
---|
629 | !!! nrlmd le 10/04/2012 |
---|
630 | ALLOCATE(ale_bl_trig(klon)) |
---|
631 | !!! fin nrlmd le 10/04/2012 |
---|
632 | IF (ok_gwd_rando) THEN |
---|
633 | ALLOCATE(du_gwd_rando(klon, klev)) |
---|
634 | du_gwd_rando(:,:)=0. |
---|
635 | ENDIF |
---|
636 | IF (.not. ok_hines .and. ok_gwd_rando) THEN |
---|
637 | ALLOCATE(du_gwd_front(klon, klev)) |
---|
638 | du_gwd_front(:,:) = 0 !ym missing init |
---|
639 | ENDIF |
---|
640 | END SUBROUTINE phys_state_var_init |
---|
641 | |
---|
642 | !====================================================================== |
---|
643 | SUBROUTINE phys_state_var_end |
---|
644 | !USE dimphy |
---|
645 | USE indice_sol_mod |
---|
646 | IMPLICIT NONE |
---|
647 | include "clesphys.h" |
---|
648 | |
---|
649 | DEALLOCATE(pctsrf, ftsol, falb1, falb2) |
---|
650 | DEALLOCATE(qsol,fevap,z0m,z0h,agesno) |
---|
651 | !FC |
---|
652 | DEALLOCATE(treedrg) |
---|
653 | DEALLOCATE(rain_fall, snow_fall, solsw, solswfdiff, sollw, radsol, swradcorr) |
---|
654 | DEALLOCATE(zmea, zstd, zsig, zgam) |
---|
655 | DEALLOCATE(zthe, zpic, zval) |
---|
656 | DEALLOCATE(rugoro, t_ancien, q_ancien, clwcon, rnebcon) |
---|
657 | DEALLOCATE(qs_ancien, ql_ancien) |
---|
658 | DEALLOCATE(prw_ancien, prlw_ancien, prsw_ancien) |
---|
659 | DEALLOCATE(qtc_cv,sigt_cv) |
---|
660 | DEALLOCATE(u_ancien, v_ancien) |
---|
661 | DEALLOCATE(tr_ancien) !RomP |
---|
662 | DEALLOCATE(ratqs, pbl_tke,coefh,coefm) |
---|
663 | !nrlmd< |
---|
664 | DEALLOCATE(delta_tsurf) |
---|
665 | !>nrlmd |
---|
666 | DEALLOCATE(zmax0, f0) |
---|
667 | DEALLOCATE(sig1, w01) |
---|
668 | DEALLOCATE(entr_therm, fm_therm) |
---|
669 | DEALLOCATE(detr_therm) |
---|
670 | DEALLOCATE(clwcon0th, rnebcon0th) |
---|
671 | ! radiation outputs |
---|
672 | DEALLOCATE(swdnc0, swdn0, swdn) |
---|
673 | DEALLOCATE(swupc0, swup0, swup) |
---|
674 | DEALLOCATE(lwdnc0, lwdn0, lwdn) |
---|
675 | DEALLOCATE(lwupc0, lwup0, lwup) |
---|
676 | DEALLOCATE(SWdn200clr, SWdn200) |
---|
677 | DEALLOCATE(SWup200clr, SWup200) |
---|
678 | DEALLOCATE(LWdn200clr, LWdn200) |
---|
679 | DEALLOCATE(LWup200clr, LWup200) |
---|
680 | DEALLOCATE(LWdnTOA, LWdnTOAclr) |
---|
681 | ! pressure level |
---|
682 | DEALLOCATE(tsumSTD) |
---|
683 | DEALLOCATE(usumSTD, vsumSTD) |
---|
684 | DEALLOCATE(wsumSTD, phisumSTD) |
---|
685 | DEALLOCATE(tnondef) |
---|
686 | DEALLOCATE(qsumSTD, rhsumSTD) |
---|
687 | DEALLOCATE(uvsumSTD) |
---|
688 | DEALLOCATE(vqsumSTD) |
---|
689 | DEALLOCATE(vTsumSTD) |
---|
690 | DEALLOCATE(wqsumSTD) |
---|
691 | DEALLOCATE(vphisumSTD) |
---|
692 | DEALLOCATE(wTsumSTD) |
---|
693 | DEALLOCATE(u2sumSTD) |
---|
694 | DEALLOCATE(v2sumSTD) |
---|
695 | DEALLOCATE(T2sumSTD) |
---|
696 | DEALLOCATE(O3sumSTD) |
---|
697 | DEALLOCATE(O3daysumSTD) |
---|
698 | !IM beg |
---|
699 | DEALLOCATE(wlevSTD,ulevSTD,vlevSTD,tlevSTD,qlevSTD,rhlevSTD,philevSTD) |
---|
700 | DEALLOCATE(uvSTD,vqSTD,vTSTD,wqSTD,vphiSTD,wTSTD,u2STD,v2STD,T2STD,O3STD,O3daySTD) |
---|
701 | !IM end |
---|
702 | DEALLOCATE(seed_old) |
---|
703 | DEALLOCATE(zuthe, zvthe) |
---|
704 | DEALLOCATE(alb_neig) |
---|
705 | DEALLOCATE(ema_cbmf) |
---|
706 | DEALLOCATE(ema_pcb, ema_pct) |
---|
707 | DEALLOCATE(Mipsh, Ma, qcondc) |
---|
708 | DEALLOCATE(wd, sigd) |
---|
709 | DEALLOCATE(cin, ALE, ALP) |
---|
710 | DEALLOCATE(ftd, fqd) |
---|
711 | DEALLOCATE(Ale_bl, Alp_bl) |
---|
712 | DEALLOCATE(ale_wake) |
---|
713 | DEALLOCATE(ale_bl_stat) |
---|
714 | DEALLOCATE(lalim_conv, wght_th) |
---|
715 | DEALLOCATE(wake_deltat, wake_deltaq) |
---|
716 | DEALLOCATE(wake_s, awake_dens, wake_dens) |
---|
717 | DEALLOCATE(wake_Cstar, wake_pe, wake_fip) |
---|
718 | !jyg< |
---|
719 | DEALLOCATE(wake_delta_pbl_TKE) |
---|
720 | !>jyg |
---|
721 | DEALLOCATE(pfrac_impa, pfrac_nucl) |
---|
722 | DEALLOCATE(pfrac_1nucl) |
---|
723 | DEALLOCATE(total_rain, nday_rain) |
---|
724 | DEALLOCATE(paire_ter) |
---|
725 | DEALLOCATE(albsol1, albsol2) |
---|
726 | !albedo SB >>> |
---|
727 | DEALLOCATE(albsol_dir,albsol_dif,falb_dir,falb_dif,chl_con) |
---|
728 | !albedo SB <<< |
---|
729 | DEALLOCATE(wo) |
---|
730 | DEALLOCATE(clwcon0,rnebcon0) |
---|
731 | DEALLOCATE(heat, heat0) |
---|
732 | DEALLOCATE(cool, cool0) |
---|
733 | DEALLOCATE(heat_volc, cool_volc) |
---|
734 | DEALLOCATE(topsw, toplw) |
---|
735 | DEALLOCATE(sollwdown, sollwdownclr) |
---|
736 | DEALLOCATE(gustiness) |
---|
737 | DEALLOCATE(toplwdown, toplwdownclr) |
---|
738 | DEALLOCATE(topsw0,toplw0,solsw0,sollw0) |
---|
739 | DEALLOCATE(albpla) |
---|
740 | !IM ajout variables CFMIP2/CMIP5 |
---|
741 | DEALLOCATE(heatp, coolp) |
---|
742 | DEALLOCATE(heat0p, cool0p) |
---|
743 | DEALLOCATE(radsolp, topswp, toplwp) |
---|
744 | DEALLOCATE(albplap) |
---|
745 | DEALLOCATE(solswp, solswfdiffp, sollwp) |
---|
746 | DEALLOCATE(sollwdownp) |
---|
747 | DEALLOCATE(topsw0p,toplw0p) |
---|
748 | DEALLOCATE(solsw0p,sollw0p) |
---|
749 | DEALLOCATE(lwdnc0p, lwdn0p, lwdnp) |
---|
750 | DEALLOCATE(lwupc0p, lwup0p, lwupp) |
---|
751 | DEALLOCATE(swdnc0p, swdn0p, swdnp) |
---|
752 | DEALLOCATE(swupc0p, swup0p, swupp) |
---|
753 | DEALLOCATE(cape) |
---|
754 | DEALLOCATE(pbase,bbase) |
---|
755 | DEALLOCATE(zqasc) |
---|
756 | DEALLOCATE(ibas_con, itop_con) |
---|
757 | DEALLOCATE(rain_con, snow_con) |
---|
758 | DEALLOCATE(rlonPOS) |
---|
759 | DEALLOCATE(newsst) |
---|
760 | DEALLOCATE(ustar,u10m, v10m,wstar) |
---|
761 | DEALLOCATE(topswad, solswad) |
---|
762 | DEALLOCATE(topswai, solswai) |
---|
763 | DEALLOCATE(tau_aero,piz_aero,cg_aero) |
---|
764 | DEALLOCATE(tau_aero_sw_rrtm,piz_aero_sw_rrtm,cg_aero_sw_rrtm) |
---|
765 | DEALLOCATE(tau_aero_lw_rrtm,piz_aero_lw_rrtm,cg_aero_lw_rrtm) |
---|
766 | DEALLOCATE(ccm) |
---|
767 | if (ok_gwd_rando) DEALLOCATE(du_gwd_rando) |
---|
768 | if (.not. ok_hines .and. ok_gwd_rando) DEALLOCATE(du_gwd_front) |
---|
769 | |
---|
770 | !!! nrlmd le 10/04/2012 |
---|
771 | DEALLOCATE(ale_bl_trig) |
---|
772 | !!! fin nrlmd le 10/04/2012 |
---|
773 | is_initialized=.FALSE. |
---|
774 | END SUBROUTINE phys_state_var_end |
---|
775 | |
---|
776 | END MODULE phys_state_var_mod |
---|