source: LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90 @ 4582

Last change on this file since 4582 was 4576, checked in by Laurent Fairhead, 17 months ago

Forgot modifications in phylmdiso
LF
File size: 62.5 KB
RevLine 
[3927]1!
[3940]2! $Id: phys_local_var_mod.F90 3888 2021-05-05 10:50:37Z jyg $
[3927]3!
4      MODULE phys_local_var_mod
5! Variables locales pour effectuer les appels en serie
6!======================================================================
7!
8!
9!======================================================================
10! Declaration des variables
11
12      REAL, SAVE, ALLOCATABLE :: t_seri(:,:), q_seri(:,:)
13      !$OMP THREADPRIVATE(t_seri, q_seri)
14      REAL, SAVE, ALLOCATABLE :: ql_seri(:,:),qs_seri(:,:)
15      !$OMP THREADPRIVATE(ql_seri,qs_seri)
[4523]16      REAL, SAVE, ALLOCATABLE :: qbs_seri(:,:)
17      !$OMP THREADPRIVATE(qbs_seri)
[3927]18      REAL, SAVE, ALLOCATABLE :: u_seri(:,:), v_seri(:,:)
19      !$OMP THREADPRIVATE(u_seri, v_seri)
[4065]20      REAL, SAVE, ALLOCATABLE :: rneb_seri(:,:)
21      !$OMP THREADPRIVATE(rneb_seri)
22      REAL, SAVE, ALLOCATABLE :: d_rneb_dyn(:,:)
23      !$OMP THREADPRIVATE(d_rneb_dyn)
24      REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:),l_mix(:,:,:),tke_dissip(:,:,:),wprime(:,:,:)
[4040]25      !$OMP THREADPRIVATE(l_mixmin, l_mix, tke_dissip,wprime)
[3927]26      REAL, SAVE, ALLOCATABLE :: tr_seri(:,:,:)
27      !$OMP THREADPRIVATE(tr_seri)
28      REAL, SAVE, ALLOCATABLE :: d_t_dyn(:,:), d_q_dyn(:,:)
29      !$OMP THREADPRIVATE(d_t_dyn, d_q_dyn)
[4523]30      REAL, SAVE, ALLOCATABLE :: d_ql_dyn(:,:), d_qs_dyn(:,:), d_qbs_dyn(:,:)
31      !$OMP THREADPRIVATE(d_ql_dyn, d_qs_dyn, d_qbs_dyn)
32      REAL, SAVE, ALLOCATABLE :: d_q_dyn2d(:), d_ql_dyn2d(:), d_qs_dyn2d(:), d_qbs_dyn2d(:)
33      !$OMP THREADPRIVATE(d_q_dyn2d, d_ql_dyn2d, d_qs_dyn2d, d_qbs_dyn2d)
[3927]34      REAL, SAVE, ALLOCATABLE :: d_u_dyn(:,:), d_v_dyn(:,:)
35      !$OMP THREADPRIVATE(d_u_dyn, d_v_dyn)
36      REAL, SAVE, ALLOCATABLE :: d_tr_dyn(:,:,:)
37      !$OMP THREADPRIVATE(d_tr_dyn)
38      REAL, SAVE, ALLOCATABLE :: d_t_con(:,:),d_q_con(:,:)
39      !$OMP THREADPRIVATE(d_t_con,d_q_con)
40      REAL, SAVE, ALLOCATABLE :: d_u_con(:,:),d_v_con(:,:)
41      !$OMP THREADPRIVATE(d_u_con,d_v_con)
42      REAL, SAVE, ALLOCATABLE :: d_t_wake(:,:),d_q_wake(:,:)
43      !$OMP THREADPRIVATE( d_t_wake,d_q_wake)
44      REAL, SAVE, ALLOCATABLE :: d_t_lsc(:,:),d_q_lsc(:,:),d_ql_lsc(:,:),d_qi_lsc(:,:)
45      !$OMP THREADPRIVATE(d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc)
46      REAL, SAVE, ALLOCATABLE :: d_t_lwr(:,:),d_t_lw0(:,:),d_t_swr(:,:),d_t_sw0(:,:)
47      !$OMP THREADPRIVATE(d_t_lwr,d_t_lw0,d_t_swr,d_t_sw0)
48      REAL, SAVE, ALLOCATABLE :: d_t_ajsb(:,:), d_q_ajsb(:,:)
49      !$OMP THREADPRIVATE(d_t_ajsb, d_q_ajsb)
50      REAL, SAVE, ALLOCATABLE :: d_t_ajs(:,:), d_q_ajs(:,:)
51      !$OMP THREADPRIVATE(d_t_ajs, d_q_ajs)
52      REAL, SAVE, ALLOCATABLE :: d_u_ajs(:,:), d_v_ajs(:,:)
53      !$OMP THREADPRIVATE(d_u_ajs, d_v_ajs)
54!nrlmd<
55      REAL, SAVE, ALLOCATABLE :: d_t_ajs_w(:,:), d_q_ajs_w(:,:)
56      !$OMP THREADPRIVATE(d_t_ajs_w, d_q_ajs_w)
57      REAL, SAVE, ALLOCATABLE :: d_t_ajs_x(:,:), d_q_ajs_x(:,:)
58      !$OMP THREADPRIVATE(d_t_ajs_x, d_q_ajs_x)
59!>nrlmd
60      REAL, SAVE, ALLOCATABLE :: d_t_eva(:,:),d_q_eva(:,:),d_ql_eva(:,:),d_qi_eva(:,:)
61      !$OMP THREADPRIVATE(d_t_eva,d_q_eva,d_ql_eva,d_qi_eva)
62      REAL, SAVE, ALLOCATABLE :: d_t_lscst(:,:),d_q_lscst(:,:)
63      !$OMP THREADPRIVATE(d_t_lscst,d_q_lscst)
64      REAL, SAVE, ALLOCATABLE :: d_t_lscth(:,:),d_q_lscth(:,:)
65      !$OMP THREADPRIVATE(d_t_lscth,d_q_lscth)
66      REAL, SAVE, ALLOCATABLE :: plul_th(:),plul_st(:)
67      !$OMP THREADPRIVATE(plul_th,plul_st)
68!tendances dues a oro et lif
69      REAL, SAVE, ALLOCATABLE :: d_t_oli(:,:)
70      !$OMP THREADPRIVATE(d_t_oli)
71      REAL, SAVE, ALLOCATABLE :: d_u_oli(:,:), d_v_oli(:,:)
72      !$OMP THREADPRIVATE(d_u_oli, d_v_oli)
[4523]73      REAL, SAVE, ALLOCATABLE :: d_t_vdf(:,:), d_q_vdf(:,:), d_qbs_vdf(:,:), d_t_diss(:,:)
74      !$OMP THREADPRIVATE( d_t_vdf, d_q_vdf, d_qbs_vdf, d_t_diss)
[3927]75      REAL, SAVE, ALLOCATABLE :: d_u_vdf(:,:), d_v_vdf(:,:)
76      !$OMP THREADPRIVATE(d_u_vdf, d_v_vdf)
77!nrlmd+jyg<
78      REAL, SAVE, ALLOCATABLE :: d_t_vdf_w(:,:), d_q_vdf_w(:,:)
79      !$OMP THREADPRIVATE( d_t_vdf_w, d_q_vdf_w)
80      REAL, SAVE, ALLOCATABLE :: d_t_vdf_x(:,:), d_q_vdf_x(:,:)
81      !$OMP THREADPRIVATE( d_t_vdf_x, d_q_vdf_x)
[4523]82      REAL, SAVE, ALLOCATABLE :: d_t_bs(:,:), d_q_bs(:,:), d_qbs_bs(:,:)
83      !$OMP THREADPRIVATE( d_t_bs,d_q_bs, d_qbs_bs)
[3927]84!>nrlmd+jyg
85      REAL, SAVE, ALLOCATABLE :: d_t_oro(:,:)
86      !$OMP THREADPRIVATE(d_t_oro)
87      REAL, SAVE, ALLOCATABLE :: d_u_oro(:,:), d_v_oro(:,:)
88      !$OMP THREADPRIVATE(d_u_oro, d_v_oro)
89      REAL, SAVE, ALLOCATABLE :: d_t_oro_gw(:,:)
90      !$OMP THREADPRIVATE(d_t_oro_gw)
91      REAL, SAVE, ALLOCATABLE :: d_u_oro_gw(:,:), d_v_oro_gw(:,:)
92      !$OMP THREADPRIVATE(d_u_oro_gw, d_v_oro_gw)
93      REAL, SAVE, ALLOCATABLE :: d_t_lif(:,:)
94      !$OMP THREADPRIVATE(d_t_lif)
95      REAL, SAVE, ALLOCATABLE :: d_u_lif(:,:), d_v_lif(:,:)
96      !$OMP THREADPRIVATE(d_u_lif, d_v_lif)
97! Tendances Ondes de G non oro (runs strato).
98      REAL, SAVE, ALLOCATABLE :: du_gwd_hines(:,:)
99      !$OMP THREADPRIVATE(du_gwd_hines)
100      REAL, SAVE, ALLOCATABLE :: dv_gwd_hines(:,:)
101      !$OMP THREADPRIVATE(dv_gwd_hines)
102      REAL, SAVE, ALLOCATABLE :: dv_gwd_rando(:,:)
103      !$OMP THREADPRIVATE(dv_gwd_rando)
104      REAL, SAVE, ALLOCATABLE :: dv_gwd_front(:,:)
105      !$OMP THREADPRIVATE(dv_gwd_front)
106      REAL, SAVE, ALLOCATABLE :: east_gwstress(:,:)
107      !$OMP THREADPRIVATE(east_gwstress)
108      REAL, SAVE, ALLOCATABLE :: west_gwstress(:,:)
109      !$OMP THREADPRIVATE(west_gwstress)
110      REAL, SAVE, ALLOCATABLE :: d_t_hin(:,:)
111      !$OMP THREADPRIVATE(d_t_hin)
112! tendance due a l'oxydation du methane
113      REAL, SAVE, ALLOCATABLE :: d_q_ch4(:,:)
114      !$OMP THREADPRIVATE(d_q_ch4)
115#ifdef ISO
116      REAL, SAVE, ALLOCATABLE :: xt_seri(:,:,:)
117      !$OMP THREADPRIVATE( xt_seri)
118      REAL, SAVE, ALLOCATABLE :: xtl_seri(:,:,:)
119      !$OMP THREADPRIVATE( xtl_seri)
120      REAL, SAVE, ALLOCATABLE :: xts_seri(:,:,:)
121      !$OMP THREADPRIVATE( xts_seri)
122      REAL, SAVE, ALLOCATABLE :: d_xt_eva(:,:,:)
123      !$OMP THREADPRIVATE( d_xt_eva)
124      REAL, SAVE, ALLOCATABLE :: d_xtl_eva(:,:,:)
125      !$OMP THREADPRIVATE( d_xtl_eva)
126      REAL, SAVE, ALLOCATABLE :: d_xti_eva(:,:,:)
127      !$OMP THREADPRIVATE( d_xti_eva)
128      REAL, SAVE, ALLOCATABLE :: d_xt_vdf(:,:,:)
129      !$OMP THREADPRIVATE( d_xt_vdf)
130      REAL, SAVE, ALLOCATABLE :: d_xt_dyn(:,:,:)
131      !$OMP THREADPRIVATE( d_xt_dyn)
132      REAL, SAVE, ALLOCATABLE :: d_xtl_dyn(:,:,:), d_xts_dyn(:,:,:)
133      !$OMP THREADPRIVATE(d_xtl_dyn, d_xts_dyn)
134      REAL, SAVE, ALLOCATABLE :: d_xt_con(:,:,:)
135      !$OMP THREADPRIVATE( d_xt_con)
136      REAL, SAVE, ALLOCATABLE :: d_xt_wake(:,:,:)
137      !$OMP THREADPRIVATE( d_xt_wake)
138      REAL, SAVE, ALLOCATABLE :: d_xt_lsc(:,:,:),d_xtl_lsc(:,:,:),d_xti_lsc(:,:,:)
139      !$OMP THREADPRIVATE( d_xt_lsc,d_xtl_lsc,d_xti_lsc)
140      REAL, SAVE, ALLOCATABLE :: d_xt_ajsb(:,:,:)
141      !$OMP THREADPRIVATE( d_xt_ajsb)
142      REAL, SAVE, ALLOCATABLE :: d_xt_ajs(:,:,:)
143      !$OMP THREADPRIVATE( d_xt_ajs)
144      REAL, SAVE, ALLOCATABLE :: d_xt_ajs_w(:,:,:), d_xt_ajs_x(:,:,:)
145      !$OMP THREADPRIVATE(d_xt_ajs_w, d_xt_ajs_x)
146      REAL, SAVE, ALLOCATABLE :: d_xt_vdf_w(:,:,:), d_xt_vdf_x(:,:,:)
147      !$OMP THREADPRIVATE(d_xt_vdf_w, d_xt_vdf_x)
148      REAL, SAVE, ALLOCATABLE :: d_xt_ch4(:,:,:)
149      !$OMP THREADPRIVATE( d_xt_ch4)
150      REAL, SAVE, ALLOCATABLE :: d_xt_prod_nucl(:,:,:)
151      !$OMP THREADPRIVATE( d_xt_prod_nucl)
152      REAL, SAVE, ALLOCATABLE :: d_xt_cosmo(:,:,:)
153      !$OMP THREADPRIVATE( d_xt_cosmo)
154      REAL, SAVE, ALLOCATABLE :: d_xt_decroiss(:,:,:)
155      !$OMP THREADPRIVATE( d_xt_decroiss)
156#endif
157
158! tendance du a la conersion Ec -> E thermique
159      REAL, SAVE, ALLOCATABLE :: d_t_ec(:,:)
160      !$OMP THREADPRIVATE(d_t_ec)
161      REAL, SAVE, ALLOCATABLE :: d_ts(:,:), d_tr(:,:,:)
162      !$OMP THREADPRIVATE(d_ts, d_tr)
163
164! diagnostique pour le rayonnement
165      REAL, SAVE, ALLOCATABLE :: topswad_aero(:),  solswad_aero(:)      ! diag
166      !$OMP THREADPRIVATE(topswad_aero,solswad_aero)
167      REAL, SAVE, ALLOCATABLE :: topswai_aero(:),  solswai_aero(:)      ! diag
168      !$OMP THREADPRIVATE(topswai_aero,solswai_aero)
169      REAL, SAVE, ALLOCATABLE :: topswad0_aero(:), solswad0_aero(:)     ! diag
170      !$OMP THREADPRIVATE(topswad0_aero,solswad0_aero)
171      REAL, SAVE, ALLOCATABLE :: topsw_aero(:,:),  solsw_aero(:,:)      ! diag
172      !$OMP THREADPRIVATE(topsw_aero,solsw_aero)
173      REAL, SAVE, ALLOCATABLE :: topsw0_aero(:,:), solsw0_aero(:,:)     ! diag
174      !$OMP THREADPRIVATE(topsw0_aero,solsw0_aero)
175      REAL, SAVE, ALLOCATABLE :: topswcf_aero(:,:),  solswcf_aero(:,:)  ! diag
176      !$OMP THREADPRIVATE(topswcf_aero,solswcf_aero)
177! LW radiation diagnostics CK
178      REAL, SAVE, ALLOCATABLE :: toplwad_aero(:),  sollwad_aero(:)      ! diag
179      !$OMP THREADPRIVATE(toplwad_aero,sollwad_aero)
180      REAL, SAVE, ALLOCATABLE :: toplwai_aero(:),  sollwai_aero(:)      ! diag
181      !$OMP THREADPRIVATE(toplwai_aero,sollwai_aero)
182      REAL, SAVE, ALLOCATABLE :: toplwad0_aero(:), sollwad0_aero(:)     ! diag
183      !$OMP THREADPRIVATE(toplwad0_aero,sollwad0_aero)
184! Special RRTM
185      REAL, SAVE, ALLOCATABLE :: ZLWFT0_i(:,:),  ZSWFT0_i(:,:)      ! diag
186      !$OMP THREADPRIVATE(ZLWFT0_i,ZSWFT0_i)
187      REAL, SAVE, ALLOCATABLE :: ZFLDN0(:,:),  ZFLUP0(:,:)      ! diag
188      !$OMP THREADPRIVATE(ZFLDN0,ZFLUP0)
189      REAL, SAVE, ALLOCATABLE :: ZFSDN0(:,:),  ZFSUP0(:,:)      ! diag
190      !$OMP THREADPRIVATE(ZFSDN0,ZFSUP0)
191!
192      REAL, SAVE, ALLOCATABLE :: stratomask(:,:)
193      !$OMP THREADPRIVATE(stratomask)
194      REAL, SAVE, ALLOCATABLE :: tausum_aero(:,:,:)
195      !$OMP THREADPRIVATE(tausum_aero)
196      REAL, SAVE, ALLOCATABLE :: drytausum_aero(:,:)
197      !$OMP THREADPRIVATE(drytausum_aero)
198      REAL, SAVE, ALLOCATABLE :: tau3d_aero(:,:,:,:)
199      !$OMP THREADPRIVATE(tau3d_aero)
200      REAL, SAVE, ALLOCATABLE :: scdnc(:,:)
201      !$OMP THREADPRIVATE(scdnc)
202      REAL, SAVE, ALLOCATABLE :: dNovrN(:)
203      !$OMP THREADPRIVATE(dNovrN)
204      REAL, SAVE, ALLOCATABLE :: cldncl(:)
205      !$OMP THREADPRIVATE(cldncl)
206      REAL, SAVE, ALLOCATABLE :: reffclwtop(:)
207      !$OMP THREADPRIVATE(reffclwtop)
208      REAL, SAVE, ALLOCATABLE :: lcc(:)
209      !$OMP THREADPRIVATE(lcc)
210      REAL, SAVE, ALLOCATABLE :: reffclws(:,:)
211      !$OMP THREADPRIVATE(reffclws)
212      REAL, SAVE, ALLOCATABLE :: reffclwc(:,:)
213      !$OMP THREADPRIVATE(reffclwc)
214      REAL, SAVE, ALLOCATABLE :: cldnvi(:)
215      !$OMP THREADPRIVATE(cldnvi)
216      REAL, SAVE, ALLOCATABLE :: lcc3d(:,:)
217      !$OMP THREADPRIVATE(lcc3d)
218      REAL, SAVE, ALLOCATABLE :: lcc3dcon(:,:)
219      !$OMP THREADPRIVATE(lcc3dcon)
220      REAL, SAVE, ALLOCATABLE :: lcc3dstra(:,:)
221      !$OMP THREADPRIVATE(lcc3dstra)
222      REAL, SAVE, ALLOCATABLE :: icc3dcon(:,:)
223      !$OMP THREADPRIVATE(icc3dcon)
224      REAL, SAVE, ALLOCATABLE :: icc3dstra(:,:)
225      !$OMP THREADPRIVATE(icc3dstra)
226      REAL, SAVE, ALLOCATABLE :: zfice(:,:)
227      !$OMP THREADPRIVATE(zfice)
228      REAL, SAVE, ALLOCATABLE :: od443aer(:)
229      !$OMP THREADPRIVATE(od443aer)
230      REAL, SAVE, ALLOCATABLE :: od550aer(:)
231      !$OMP THREADPRIVATE(od550aer)
232      REAL, SAVE, ALLOCATABLE :: dryod550aer(:)
233      !$OMP THREADPRIVATE(dryod550aer)
234      REAL, SAVE, ALLOCATABLE :: abs550aer(:)
235      !$OMP THREADPRIVATE(abs550aer)
236      REAL, SAVE, ALLOCATABLE :: od865aer(:)
237      !$OMP THREADPRIVATE(od865aer)
238      REAL, SAVE, ALLOCATABLE :: ec550aer(:,:)
239      !$OMP THREADPRIVATE(ec550aer)
240      REAL, SAVE, ALLOCATABLE :: od550lt1aer(:)
241      !$OMP THREADPRIVATE(od550lt1aer)
242      REAL, SAVE, ALLOCATABLE :: sconcso4(:)
243      !$OMP THREADPRIVATE(sconcso4)
244      REAL, SAVE, ALLOCATABLE :: sconcno3(:)
245      !$OMP THREADPRIVATE(sconcno3)
246      REAL, SAVE, ALLOCATABLE :: sconcoa(:)
247      !$OMP THREADPRIVATE(sconcoa)
248      REAL, SAVE, ALLOCATABLE :: sconcbc(:)
249      !$OMP THREADPRIVATE(sconcbc)
250      REAL, SAVE, ALLOCATABLE :: sconcss(:)
251      !$OMP THREADPRIVATE(sconcss)
252      REAL, SAVE, ALLOCATABLE :: sconcdust(:)
253      !$OMP THREADPRIVATE(sconcdust)
254      REAL, SAVE, ALLOCATABLE :: concso4(:,:)
255      !$OMP THREADPRIVATE(concso4)
256      REAL, SAVE, ALLOCATABLE :: concno3(:,:)
257      !$OMP THREADPRIVATE(concno3)
258      REAL, SAVE, ALLOCATABLE :: concoa(:,:)
259      !$OMP THREADPRIVATE(concoa)
260      REAL, SAVE, ALLOCATABLE :: concbc(:,:)
261      !$OMP THREADPRIVATE(concbc)
262      REAL, SAVE, ALLOCATABLE :: concss(:,:)
263      !$OMP THREADPRIVATE(concss)
264      REAL, SAVE, ALLOCATABLE :: concdust(:,:)
265      !$OMP THREADPRIVATE(concdust)
266      REAL, SAVE, ALLOCATABLE :: loadso4(:)
267      !$OMP THREADPRIVATE(loadso4)
268      REAL, SAVE, ALLOCATABLE :: loadoa(:)
269      !$OMP THREADPRIVATE(loadoa)
270      REAL, SAVE, ALLOCATABLE :: loadbc(:)
271      !$OMP THREADPRIVATE(loadbc)
272      REAL, SAVE, ALLOCATABLE :: loadss(:)
273      !$OMP THREADPRIVATE(loadss)
274      REAL, SAVE, ALLOCATABLE :: loaddust(:)
275      !$OMP THREADPRIVATE(loaddust)
276      REAL, SAVE, ALLOCATABLE :: loadno3(:)
277      !$OMP THREADPRIVATE(loadno3)
278      REAL, SAVE, ALLOCATABLE :: load_tmp1(:), load_tmp2(:)
279      !$OMP THREADPRIVATE(load_tmp1, load_tmp2)
280      REAL, SAVE, ALLOCATABLE :: load_tmp3(:), load_tmp4(:)
281      !$OMP THREADPRIVATE(load_tmp3, load_tmp4)
282      REAL, SAVE, ALLOCATABLE :: load_tmp5(:), load_tmp6(:)
283      !$OMP THREADPRIVATE(load_tmp5, load_tmp6)
284      REAL, SAVE, ALLOCATABLE :: load_tmp7(:), load_tmp8(:)
285      !$OMP THREADPRIVATE(load_tmp7, load_tmp8)
286      REAL, SAVE, ALLOCATABLE :: load_tmp9(:), load_tmp10(:)
287      !$OMP THREADPRIVATE(load_tmp9, load_tmp10)
288
289!IM ajout variables CFMIP2/CMIP5
290      REAL,ALLOCATABLE,SAVE :: topswad_aerop(:), solswad_aerop(:)
291!$OMP THREADPRIVATE(topswad_aerop, solswad_aerop)
292      REAL,ALLOCATABLE,SAVE :: topswai_aerop(:), solswai_aerop(:)
293!$OMP THREADPRIVATE(topswai_aerop, solswai_aerop)
294      REAL,ALLOCATABLE,SAVE :: topswad0_aerop(:), solswad0_aerop(:)
295!$OMP THREADPRIVATE(topswad0_aerop, solswad0_aerop)
296      REAL,ALLOCATABLE,SAVE :: topsw_aerop(:,:), topsw0_aerop(:,:) 
297!$OMP THREADPRIVATE(topsw_aerop, topsw0_aerop)
298      REAL,ALLOCATABLE,SAVE :: solsw_aerop(:,:), solsw0_aerop(:,:)
299!$OMP THREADPRIVATE(solsw_aerop, solsw0_aerop)
300      REAL,ALLOCATABLE,SAVE :: topswcf_aerop(:,:), solswcf_aerop(:,:)
301!$OMP THREADPRIVATE(topswcf_aerop, solswcf_aerop)
302
303! additional LW variables CK
304      REAL,ALLOCATABLE,SAVE :: toplwad_aerop(:), sollwad_aerop(:)
305!$OMP THREADPRIVATE(toplwad_aerop, sollwad_aerop)
306      REAL,ALLOCATABLE,SAVE :: toplwai_aerop(:), sollwai_aerop(:)
307!$OMP THREADPRIVATE(toplwai_aerop, sollwai_aerop)
308      REAL,ALLOCATABLE,SAVE :: toplwad0_aerop(:), sollwad0_aerop(:)
309!$OMP THREADPRIVATE(toplwad0_aerop, sollwad0_aerop)
310
311!Ajout de celles n??cessaires au phys_output_write_mod
312      REAL, SAVE, ALLOCATABLE :: tal1(:), pal1(:), pab1(:), pab2(:)
313!$OMP THREADPRIVATE(tal1, pal1, pab1, pab2)
314      REAL, SAVE, ALLOCATABLE :: ptstar(:), pt0(:), slp(:)
315!$OMP THREADPRIVATE(ptstar, pt0, slp)
316      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: sens, flwp, fiwp
317!$OMP THREADPRIVATE(sens, flwp, fiwp)
318!!
[4576]319!FC
320      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zxfluxt, zxfluxq
321!$OMP THREADPRIVATE(zxfluxt, zxfluxq)
322!FC
[3927]323!!         Wake variables
324      REAL,ALLOCATABLE,SAVE,DIMENSION(:)            :: alp_wake
325!$OMP THREADPRIVATE(alp_wake)           
326!!jyg!      REAL,ALLOCATABLE,SAVE,DIMENSION(:)            :: wake_h,wake_k
327      REAL,ALLOCATABLE,SAVE,DIMENSION(:)            :: wake_h
328      INTEGER,ALLOCATABLE,SAVE,DIMENSION(:)         :: wake_k
329!$OMP THREADPRIVATE(wake_h,wake_k)                       
330    REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)            :: wake_omg
331!$OMP THREADPRIVATE(wake_omg)                     
332    REAL, SAVE, ALLOCATABLE,DIMENSION(:,:)          :: d_deltat_wk, d_deltaq_wk
333!$OMP THREADPRIVATE(d_deltat_wk, d_deltaq_wk)
334      REAL,ALLOCATABLE,SAVE,DIMENSION(:)            :: d_s_wk, d_dens_a_wk, d_dens_wk
335!$OMP THREADPRIVATE(d_s_wk, d_dens_a_wk, d_dens_wk)
336    REAL, SAVE, ALLOCATABLE,DIMENSION(:,:)          :: d_deltat_wk_gw, d_deltaq_wk_gw
337!$OMP THREADPRIVATE(d_deltat_wk_gw, d_deltaq_wk_gw)
338    REAL, SAVE, ALLOCATABLE,DIMENSION(:,:)          :: d_deltat_vdf, d_deltaq_vdf
339!$OMP THREADPRIVATE(d_deltat_vdf, d_deltaq_vdf)
340!!!      REAL,ALLOCATABLE,SAVE,DIMENSION(:)          :: d_s_vdf, d_dens_vdf
341!!!OMP THREADPRIVATE(d_s_vdf, d_dens_vdf)
342    REAL, SAVE, ALLOCATABLE,DIMENSION(:,:)          :: d_deltat_the, d_deltaq_the
343!$OMP THREADPRIVATE(d_deltat_the, d_deltaq_the)
344!!!      REAL,ALLOCATABLE,SAVE,DIMENSION(:)          :: d_s_the, d_dens_the
345!!!OMP THREADPRIVATE(d_s_the, d_dens_the)
346      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)           :: d_deltat_ajs_cv, d_deltaq_ajs_cv
[3940]347!$OMP THREADPRIVATE(d_deltat_ajs_cv, d_deltaq_ajs_cv)                       
[3927]348#ifdef ISO
349    REAL, SAVE, ALLOCATABLE,DIMENSION(:,:,:)          :: d_deltaxt_wk
350!$OMP THREADPRIVATE(d_deltaxt_wk)
351    REAL, SAVE, ALLOCATABLE,DIMENSION(:,:,:)          :: d_deltaxt_wk_gw
352!$OMP THREADPRIVATE(d_deltaxt_wk_gw)
353    REAL, SAVE, ALLOCATABLE,DIMENSION(:,:,:)          ::  d_deltaxt_the
354!$OMP THREADPRIVATE(d_deltaxt_the)
355    REAL, SAVE, ALLOCATABLE,DIMENSION(:,:,:)          ::  d_deltaxt_vdf
356!$OMP THREADPRIVATE(d_deltaxt_vdf)
357      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)           ::  d_deltaxt_ajs_cv
358!$OMP THREADPRIVATE(d_deltaxt_ajs_cv)
359#endif                       
360!!         End of Wake variables
361!!
362      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: bils
363!$OMP THREADPRIVATE(bils)
364      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cdragm, cdragh
365!$OMP THREADPRIVATE(cdragm, cdragh)
366      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cldh, cldl, cldm, cldq, cldt, qsat2m
367!$OMP THREADPRIVATE(cldh, cldl, cldm, cldq, cldt, qsat2m )
[3940]368!AS: cldhjn, cldljn, cldmjn,cldtjn pas utilisés en tant que variables, juste noms de diagnostics
[3927]369      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: JrNt
370!$OMP THREADPRIVATE(JrNt)
[4523]371      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: dthmin, evap, snowerosion, fder, plcl, plfc, prw, prlw, prsw, prbsw
372!$OMP THREADPRIVATE(dthmin, evap, snowerosion, fder, plcl, plfc, prw, prlw, prsw, prbsw)
[3927]373      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zustar, zu10m, zv10m, rh2m
374!$OMP THREADPRIVATE(zustar, zu10m, zv10m, rh2m)
375      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: s_lcl, s_pblh, s_pblt, s_therm
376!$OMP THREADPRIVATE(s_lcl, s_pblh, s_pblt, s_therm)
377!
378!nrlmd+jyg<
379      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: s_pblh_x, s_pblh_w
380!$OMP THREADPRIVATE(s_pblh_x, s_pblh_w)
381      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: s_lcl_x, s_lcl_w
382!$OMP THREADPRIVATE(s_lcl_x, s_lcl_w)
383!>nrlmd+jyg
384!
385      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: slab_wfbils
386!$OMP THREADPRIVATE(slab_wfbils)
387      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: tpot, tpote, ue, uq, uwat, ve, vq, vwat, zxffonte
388!$OMP THREADPRIVATE(tpot, tpote, ue, uq, uwat, ve, vq, vwat, zxffonte)
[4523]389      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxustartlic, zxrhoslic
390!$OMP THREADPRIVATE(zxustartlic, zxrhoslic)
[3927]391      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfqcalving
392!$OMP THREADPRIVATE(zxfqcalving)
393      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfluxlat, zxtsol, snow_lsc, zxfqfonte
394!$OMP THREADPRIVATE(zxfluxlat, zxtsol, snow_lsc, zxfqfonte)
395      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxrunofflic
396!$OMP THREADPRIVATE(zxrunofflic)
397      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxqsurf, rain_lsc, rain_num
398!$OMP THREADPRIVATE(zxqsurf, rain_lsc, rain_num)
399#ifdef ISO
400      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: xtevap,xtprw
401!$OMP THREADPRIVATE(xtevap,xtprw)
402      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: h1_diag,runoff_diag
403      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: xtrunoff_diag
404!$OMP THREADPRIVATE(h1_diag,runoff_diag,xtrunoff_diag)
405      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zxfxtcalving
406!$OMP THREADPRIVATE(zxfxtcalving)
407      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: xtsnow_lsc, zxfxtfonte
408!$OMP THREADPRIVATE(xtsnow_lsc, zxfxtfonte)
409      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zxxtrunofflic
410!$OMP THREADPRIVATE(zxxtrunofflic)
411      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: xtrain_lsc
412!$OMP THREADPRIVATE(xtrain_lsc)
413#endif
414!
415!jyg+nrlmd<
416!!!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
417!                                                                          c
418!       Declarations liees a la couche limite differentiee w-x             c
419!                                                                          c
420!!!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
421      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: sens_x, sens_w
422!$OMP THREADPRIVATE(sens_x, sens_w)
423      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfluxlat_x, zxfluxlat_w
424!$OMP THREADPRIVATE(zxfluxlat_x, zxfluxlat_w)
[3940]425      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: delta_qsurf
426!$OMP THREADPRIVATE(delta_qsurf)
[3927]427!jyg<
428!!! Entrees supplementaires couche-limite
429!!      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: t_x, t_w
430!!!$OMP THREADPRIVATE(t_x, t_w)
431!!      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: q_x, q_w
432!!!$OMP THREADPRIVATE(q_x, q_w)
433!>jyg
434!!! Sorties ferret
435      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dtvdf_x, dtvdf_w
436!$OMP THREADPRIVATE(dtvdf_x, dtvdf_w)
437      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dqvdf_x, dqvdf_w
438!$OMP THREADPRIVATE(dqvdf_x, dqvdf_w)
439! Variables supplementaires dans physiq.F relative au splitting de la surface
440      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: pbl_tke_input
441!$OMP THREADPRIVATE(pbl_tke_input)
442! Entree supplementaire Thermiques :
443      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: t_therm, q_therm, u_therm, v_therm
444!$OMP THREADPRIVATE(t_therm, q_therm, u_therm, v_therm)
445      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cdragh_x, cdragh_w
446!$OMP THREADPRIVATE(cdragh_x, cdragh_w)
447      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cdragm_x, cdragm_w
448!$OMP THREADPRIVATE(cdragm_x, cdragm_w)
449      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: kh, kh_x, kh_w
450!$OMP THREADPRIVATE(kh, kh_x, kh_w)
451#ifdef ISO
452      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: dxtvdf_x, dxtvdf_w
453!$OMP THREADPRIVATE(dxtvdf_x, dxtvdf_w)
454      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: xt_therm
455!$OMP THREADPRIVATE(xt_therm)
456#endif
457!!!
458!!!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
459      LOGICAL, SAVE, ALLOCATABLE :: ptconv(:,:)
460      !$OMP THREADPRIVATE(ptconv)
461      REAL, SAVE, ALLOCATABLE :: ratqsc(:,:)
462      !$OMP THREADPRIVATE(ratqsc)
463!>jyg+nrlmd
464  !
465      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wbeff, convoccur, zmax_th, zq2m, zt2m
466!$OMP THREADPRIVATE(wbeff, convoccur, zmax_th, zq2m, zt2m)
467      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zt2m_min_mon, zt2m_max_mon
468!$OMP THREADPRIVATE(zt2m_min_mon, zt2m_max_mon)
469      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: t2m_min_mon, t2m_max_mon
470!$OMP THREADPRIVATE(t2m_min_mon, t2m_max_mon)
471      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: weak_inversion
472!$OMP THREADPRIVATE(weak_inversion)
473      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_conv,alp_bl_det
474!$OMP THREADPRIVATE(alp_bl_conv,alp_bl_det)
475      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_fluct_m,alp_bl_fluct_tke
476!$OMP THREADPRIVATE(alp_bl_fluct_m,alp_bl_fluct_tke)
477      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_stat, n2, s2
478!$OMP THREADPRIVATE(alp_bl_stat, n2, s2)
479      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: proba_notrig, random_notrig
480!$OMP THREADPRIVATE(proba_notrig, random_notrig)
481      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: fsolsw, wfbils, wfbilo
482!$OMP THREADPRIVATE(fsolsw, wfbils, wfbilo)
483      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wfevap, wfrain, wfsnow
484!$OMP THREADPRIVATE(wfevap,wfrain,wfsnow)
485      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: t2m, fluxlat, fsollw,evap_pot
486!$OMP THREADPRIVATE(t2m, fluxlat, fsollw,evap_pot)
487      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dnwd0, omega
488!$OMP THREADPRIVATE(dnwd0, omega)
489      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: epmax_diag ! epmax_cape
490!$OMP THREADPRIVATE(epmax_diag)
491!
492!  Deep convective variables used in phytrac
493      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: ep ! epmax_cape
494!$OMP THREADPRIVATE(ep)
495      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: dnwd, upwd
496!$OMP THREADPRIVATE(dnwd, upwd)
497      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: pmflxr, pmflxs
498!$OMP THREADPRIVATE(pmflxr, pmflxs)
499      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: wdtrainA, wdtrainS, wdtrainM
500!$OMP THREADPRIVATE(wdtrainA, wdtrainS, wdtrainM)
501      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: da, mp
502!$OMP THREADPRIVATE(da, mp)
503      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: wght_cvfd
504!$OMP THREADPRIVATE(wght_cvfd)
505      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:):: phi, phi2, elij, epmlmMm
506!$OMP THREADPRIVATE(phi, phi2, elij, epmlmMm)
507      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: d1a, dam
508!$OMP THREADPRIVATE(d1a, dam)
509      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: ev
510!$OMP THREADPRIVATE(ev)
511      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: qtaa
512!$OMP THREADPRIVATE(qtaa)
513      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: clw
514!$OMP THREADPRIVATE(clw)
515      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: eplaMm
516!$OMP THREADPRIVATE(eplaMm)
517      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:):: sij
518!$OMP THREADPRIVATE(sij)
519#ifdef ISO
520      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: xtwdtrainA
521!$OMP THREADPRIVATE(xtwdtrainA)
522      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: xtev
523!$OMP THREADPRIVATE(xtev)
524      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: xttaa
525!$OMP THREADPRIVATE(xttaa)
526      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: xtclw
527!$OMP THREADPRIVATE(xtclw)
528#ifdef DIAGISO
529      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: qlp
530!$OMP THREADPRIVATE(qlp)
531      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: qvp
532!$OMP THREADPRIVATE(qvp)
533      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: fq_detrainement
534!$OMP THREADPRIVATE(fq_detrainement)
535      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: fq_ddft
536!$OMP THREADPRIVATE(fq_ddft)
537      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: fq_fluxmasse
538!$OMP THREADPRIVATE(fq_fluxmasse)
539      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: fq_evapprecip
540!$OMP THREADPRIVATE(fq_evapprecip)
541      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: f_detrainement
542!$OMP THREADPRIVATE(f_detrainement)
543      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: q_detrainement
544!$OMP THREADPRIVATE(q_detrainement)
545      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: xt_detrainement
546!$OMP THREADPRIVATE(xt_detrainement)
547      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: xtlp
548!$OMP THREADPRIVATE(xtlp)
549      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: xtvp
550!$OMP THREADPRIVATE(xtvp)
551      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: q_the
552!$OMP THREADPRIVATE(q_the)
553      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: xt_the
554!$OMP THREADPRIVATE(xt_the)
555      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: fxt_detrainement
556!$OMP THREADPRIVATE(fxt_detrainement)
557      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: fxt_ddft
558!$OMP THREADPRIVATE(fxt_ddft)
559      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: fxt_fluxmasse
560!$OMP THREADPRIVATE(fxt_fluxmasse)
561      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:)  :: fxt_evapprecip
562!$OMP THREADPRIVATE(fxt_evapprecip)
563#endif
564#endif
565!
566!      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: coefh, coefm, lambda_th
567      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: lambda_th
568!$OMP THREADPRIVATE(lambda_th)
569      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: cldemi, cldfra, cldtau, fiwc, fl, re, flwc
570!$OMP THREADPRIVATE(cldemi, cldfra, cldtau, fiwc, fl, re, flwc)
[4118]571      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: qlth, qith, qsith, wiceth
572!$OMP THREADPRIVATE(qlth, qith, qsith, wiceth)
[3927]573      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq, ref_ice, theta, zphi
574!$OMP THREADPRIVATE(ref_liq, ref_ice, theta, zphi)
575      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq_pi, ref_ice_pi
576!$OMP THREADPRIVATE(ref_liq_pi, ref_ice_pi)
[3940]577      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zx_rh, zx_rhl, zx_rhi
578!$OMP THREADPRIVATE(zx_rh, zx_rhl, zx_rhi)
[4523]579      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: prfl, psfl, fraca, bsfl
580!$OMP THREADPRIVATE(prfl, psfl, fraca, bsfl)
[3927]581      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: Vprecip, zw2
582!$OMP THREADPRIVATE(Vprecip, zw2)
583      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: pmfd, pmfu
584!$OMP THREADPRIVATE(pmfd, pmfu)
585      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: fluxt, fluxu, fluxv
586!$OMP THREADPRIVATE(fluxt, fluxu, fluxv)
587      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: uwriteSTD, vwriteSTD, wwriteSTD
588!$OMP THREADPRIVATE(uwriteSTD, vwriteSTD, wwriteSTD)
589      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: phiwriteSTD, qwriteSTD, twriteSTD, rhwriteSTD
590!$OMP THREADPRIVATE(phiwriteSTD, qwriteSTD, twriteSTD, rhwriteSTD)
591
592! ug et d'autres encore:
593      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: beta_prec
594!$OMP THREADPRIVATE(beta_prec)
595      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rneb,rnebjn,rneblsvol
596!$OMP THREADPRIVATE(rneb,rnebjn,rneblsvol)
[4530]597      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: pfraclr,pfracld
598!$OMP THREADPRIVATE(pfraclr,pfracld)
[4562]599      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: distcltop
600!$OMP THREADPRIVATE(distcltop)
[3927]601
602! variables de sorties MM
603      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: zxsnow,snowhgt,qsnow,to_ice
604      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: sissnow,runoff,albsol3_lic
605!$OMP THREADPRIVATE(zxsnow,snowhgt,qsnow,to_ice)
606!$OMP THREADPRIVATE(sissnow,runoff,albsol3_lic)
607#ifdef ISO
608      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: zxxtsnow
609!$OMP THREADPRIVATE(zxxtsnow)     
610      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: xtVprecip,xtVprecipi
611!$OMP THREADPRIVATE(xtVprecip,xtVprecipi)
612      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: pxtrfl, pxtsfl
613!$OMP THREADPRIVATE(pxtrfl, pxtsfl)
614#endif
615
616      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: p_tropopause, z_tropopause, t_tropopause
617!$OMP THREADPRIVATE(p_tropopause, z_tropopause, t_tropopause)
618
[3940]619      INTEGER,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zn2mout
620!$OMP THREADPRIVATE(zn2mout)
621
[4065]622      REAL, SAVE, ALLOCATABLE :: qclr(:,:)
623      !$OMP THREADPRIVATE(qclr)
624      REAL, SAVE, ALLOCATABLE :: qcld(:,:)
625      !$OMP THREADPRIVATE(qcld)
626      REAL, SAVE, ALLOCATABLE :: qss(:,:)
627      !$OMP THREADPRIVATE(qss)
628      REAL, SAVE, ALLOCATABLE :: qvc(:,:)
629      !$OMP THREADPRIVATE(qvc)
630      REAL, SAVE, ALLOCATABLE :: rnebclr(:,:)
631      !$OMP THREADPRIVATE(rnebclr)
632      REAL, SAVE, ALLOCATABLE :: rnebss(:,:)
633      !$OMP THREADPRIVATE(rnebss)
634      REAL, SAVE, ALLOCATABLE :: gamma_ss(:,:)
635      !$OMP THREADPRIVATE(gamma_ss)
636      REAL, SAVE, ALLOCATABLE :: N1_ss(:,:)
637      !$OMP THREADPRIVATE(N1_ss)
638      REAL, SAVE, ALLOCATABLE :: N2_ss(:,:)
639      !$OMP THREADPRIVATE(N2_ss)
640      REAL, SAVE, ALLOCATABLE :: drneb_sub(:,:)
641      !$OMP THREADPRIVATE(drneb_sub)
642      REAL, SAVE, ALLOCATABLE :: drneb_con(:,:)
643      !$OMP THREADPRIVATE(drneb_con)
644      REAL, SAVE, ALLOCATABLE :: drneb_tur(:,:)
645      !$OMP THREADPRIVATE(drneb_tur)
646      REAL, SAVE, ALLOCATABLE :: drneb_avi(:,:)
647      !$OMP THREADPRIVATE(drneb_avi)
648      REAL, SAVE, ALLOCATABLE :: zqsatl(:,:)
649      !$OMP THREADPRIVATE(zqsatl)
650      REAL, SAVE, ALLOCATABLE :: zqsats(:,:)
651      !$OMP THREADPRIVATE(zqsats)
652      REAL, SAVE, ALLOCATABLE :: Tcontr(:,:)
653      !$OMP THREADPRIVATE(Tcontr)
654      REAL, SAVE, ALLOCATABLE :: qcontr(:,:)
655      !$OMP THREADPRIVATE(qcontr)
656      REAL, SAVE, ALLOCATABLE :: qcontr2(:,:)
657      !$OMP THREADPRIVATE(qcontr2)
658      REAL, SAVE, ALLOCATABLE :: fcontrN(:,:)
659      !$OMP THREADPRIVATE(fcontrN)
660      REAL, SAVE, ALLOCATABLE :: fcontrP(:,:)
661      !$OMP THREADPRIVATE(fcontrP)
662
[3927]663#ifdef CPP_StratAer
664!
665! variables for stratospheric aerosol
666      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: R2SO4
667!$OMP THREADPRIVATE(R2SO4)
668      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: DENSO4
669!$OMP THREADPRIVATE(DENSO4)
670      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: f_r_wet
671!$OMP THREADPRIVATE(f_r_wet)
672      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: decfluxaer
673!$OMP THREADPRIVATE(decfluxaer)
674      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: mdw
675!$OMP THREADPRIVATE(mdw)
676      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: OCS_lifetime
677!$OMP THREADPRIVATE(OCS_lifetime)
678      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: SO2_lifetime
679!$OMP THREADPRIVATE(SO2_lifetime)
680      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: alpha_bin
681!$OMP THREADPRIVATE(alpha_bin)
682      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: piz_bin
683!$OMP THREADPRIVATE(piz_bin)
684      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: cg_bin
685!$OMP THREADPRIVATE(cg_bin)
686      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_550
687!$OMP THREADPRIVATE(tau_strat_550)
688      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_1020
689!$OMP THREADPRIVATE(tau_strat_1020)
690      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tausum_strat
691!$OMP THREADPRIVATE(tausum_strat)
692      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: surf_PM25_sulf
693!$OMP THREADPRIVATE(surf_PM25_sulf)
694      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: vsed_aer
695!$OMP THREADPRIVATE(vsed_aer)
696!
697!---3D budget variables
698      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_nucl
699!$OMP THREADPRIVATE(budg_3D_nucl)
700      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_cond_evap
701!$OMP THREADPRIVATE(budg_3D_cond_evap)
702      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_ocs_to_so2
703!$OMP THREADPRIVATE(budg_3D_ocs_to_so2)
704      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_so2_to_h2so4
705!$OMP THREADPRIVATE(budg_3D_so2_to_h2so4)
706      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_backgr_ocs
707!$OMP THREADPRIVATE(budg_3D_backgr_ocs)
708      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_backgr_so2
709!$OMP THREADPRIVATE(budg_3D_backgr_so2)
710!
711!---2D budget variables
712      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_ocs
713      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_so2
714      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_h2so4
715      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_part
716!$OMP THREADPRIVATE(budg_dep_dry_ocs,budg_dep_dry_so2,budg_dep_dry_h2so4,budg_dep_dry_part)
717      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_ocs
718      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_so2
719      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_h2so4
720      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_part
721!$OMP THREADPRIVATE(budg_dep_wet_ocs,budg_dep_wet_so2,budg_dep_wet_h2so4,budg_dep_wet_part)
722      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_ocs
723      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_so2
724      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_h2so4
725      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_part
726!$OMP THREADPRIVATE(budg_emi_ocs,budg_emi_so2,budg_emi_h2so4,budg_emi_part)
727      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_ocs_to_so2
728      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_so2_to_h2so4
729      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_h2so4_to_part
730!$OMP THREADPRIVATE(budg_ocs_to_so2,budg_so2_to_h2so4,budg_h2so4_to_part)
731      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_sed_part
732!$OMP THREADPRIVATE(budg_sed_part)
733#endif
734
735CONTAINS
736
737!======================================================================
738SUBROUTINE phys_local_var_init
739USE dimphy
740USE infotrac_phy, ONLY : nbtr
741#ifdef ISO
[4143]742USE infotrac_phy, ONLY : ntraciso=>ntiso,niso
[3927]743#endif
744USE aero_mod
745USE indice_sol_mod
746USE phys_output_var_mod
747USE phys_state_var_mod
748
749IMPLICIT NONE
[4523]750      ALLOCATE(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev), qbs_seri(klon,klev))
[3927]751      ALLOCATE(u_seri(klon,klev),v_seri(klon,klev))
[4040]752      ALLOCATE(l_mixmin(klon,klev+1,nbsrf),l_mix(klon,klev+1,nbsrf),tke_dissip(klon,klev+1,nbsrf),wprime(klon,klev+1,nbsrf))
753      l_mix(:,:,:)=0.;l_mixmin(:,:,:)=0.;tke_dissip(:,:,:)=0.;wprime(:,:,:)=0. ! doit etre initialse car pas toujours remplis
[3927]754
755      ALLOCATE(tr_seri(klon,klev,nbtr))
756      ALLOCATE(d_t_dyn(klon,klev),d_q_dyn(klon,klev))
[4523]757      ALLOCATE(d_ql_dyn(klon,klev),d_qs_dyn(klon,klev), d_qbs_dyn(klon,klev))
758      ALLOCATE(d_q_dyn2d(klon),d_ql_dyn2d(klon),d_qs_dyn2d(klon), d_qbs_dyn2d(klon))
[3927]759      ALLOCATE(d_u_dyn(klon,klev),d_v_dyn(klon,klev))
760      ALLOCATE(d_tr_dyn(klon,klev,nbtr))                   !RomP
761      ALLOCATE(d_t_con(klon,klev),d_q_con(klon,klev))
762      ALLOCATE(d_u_con(klon,klev),d_v_con(klon,klev))
763      ALLOCATE(d_t_wake(klon,klev),d_q_wake(klon,klev))
764      ALLOCATE(d_t_lsc(klon,klev),d_q_lsc(klon,klev))
765      ALLOCATE(d_t_lwr(klon,klev),d_t_lw0(klon,klev))
766      ALLOCATE(d_t_swr(klon,klev),d_t_sw0(klon,klev))
767      ALLOCATE(d_ql_lsc(klon,klev),d_qi_lsc(klon,klev))
768      ALLOCATE(d_t_ajsb(klon,klev),d_q_ajsb(klon,klev))
769      ALLOCATE(d_t_ajs(klon,klev),d_q_ajs(klon,klev))
770!nrlmd<
771      ALLOCATE(d_t_ajs_w(klon,klev),d_q_ajs_w(klon,klev))
772      ALLOCATE(d_t_ajs_x(klon,klev),d_q_ajs_x(klon,klev))
773!>nrlmd
774      ALLOCATE(d_u_ajs(klon,klev),d_v_ajs(klon,klev))
775      ALLOCATE(d_t_eva(klon,klev),d_q_eva(klon,klev))
776      ALLOCATE(d_ql_eva(klon,klev),d_qi_eva(klon,klev))
777      ALLOCATE(d_t_lscst(klon,klev),d_q_lscst(klon,klev))
778      ALLOCATE(d_t_lscth(klon,klev),d_q_lscth(klon,klev))
779      ALLOCATE(plul_st(klon),plul_th(klon))
780      ALLOCATE(d_t_vdf(klon,klev),d_q_vdf(klon,klev),d_t_diss(klon,klev))
[4523]781      ALLOCATE (d_qbs_vdf(klon,klev))
782      ALLOCATE(d_t_bs(klon,klev),d_q_bs(klon,klev),d_qbs_bs(klon,klev))
[3927]783      ALLOCATE(d_t_vdf_w(klon,klev),d_q_vdf_w(klon,klev))
784      ALLOCATE(d_t_vdf_x(klon,klev),d_q_vdf_x(klon,klev))
785#ifdef ISO
786      allocate(xt_seri(ntraciso,klon,klev))
787      allocate(xtl_seri(ntraciso,klon,klev))
788      allocate(xts_seri(ntraciso,klon,klev))
789      allocate(d_xt_dyn(ntraciso,klon,klev))
790      allocate(d_xtl_dyn(ntraciso,klon,klev))
791      allocate(d_xts_dyn(ntraciso,klon,klev))
792      allocate(d_xt_con(ntraciso,klon,klev))
793      allocate(d_xt_wake(ntraciso,klon,klev))
794      allocate(d_xt_lsc(ntraciso,klon,klev))
795      allocate(d_xtl_lsc(ntraciso,klon,klev))
796      allocate(d_xti_lsc(ntraciso,klon,klev))
797      allocate(d_xt_ajsb(ntraciso,klon,klev))
798      allocate(d_xt_ajs(ntraciso,klon,klev))
799      allocate(d_xt_ajs_w(ntraciso,klon,klev))
800      allocate(d_xt_ajs_x(ntraciso,klon,klev))
801      allocate(d_xt_eva(ntraciso,klon,klev))
802      allocate(d_xtl_eva(ntraciso,klon,klev))
803      allocate(d_xti_eva(ntraciso,klon,klev))
804      allocate(d_xt_vdf(ntraciso,klon,klev)) 
805      allocate(d_xt_vdf_w(ntraciso,klon,klev))
806      allocate(d_xt_vdf_x(ntraciso,klon,klev))
807      allocate(d_xt_ch4(ntraciso,klon,klev))
808      allocate(d_xt_prod_nucl(ntraciso,klon,klev))
809      allocate(d_xt_cosmo(ntraciso,klon,klev))
810      allocate(d_xt_decroiss(ntraciso,klon,klev))
811#endif
812
813      ALLOCATE(d_u_vdf(klon,klev),d_v_vdf(klon,klev))
814      ALLOCATE(d_t_oli(klon,klev),d_t_oro(klon,klev))
815      ALLOCATE(d_u_oli(klon,klev),d_v_oli(klon,klev))
816      ALLOCATE(d_u_oro(klon,klev),d_v_oro(klon,klev))
817      ALLOCATE(d_u_oro_gw(klon,klev),d_v_oro_gw(klon,klev))
818      ALLOCATE(d_t_oro_gw(klon,klev))
819      ALLOCATE(d_t_lif(klon,klev),d_t_ec(klon,klev))
820      ALLOCATE(d_u_lif(klon,klev),d_v_lif(klon,klev))
821      ALLOCATE(d_ts(klon,nbsrf), d_tr(klon,klev,nbtr))
822! Special RRTM
823      ALLOCATE(ZLWFT0_i(klon,klev+1),ZSWFT0_i(klon,klev+1),ZFLDN0(klon,klev+1))
824      ZFLDN0= 0.
825      ALLOCATE(ZFLUP0(klon,klev+1),ZFSDN0(klon,klev+1),ZFSUP0(klon,klev+1))
826!
827      ALLOCATE(topswad_aero(klon), solswad_aero(klon))
828      ALLOCATE(topswai_aero(klon), solswai_aero(klon))
829      ALLOCATE(topswad0_aero(klon), solswad0_aero(klon))
830      ALLOCATE(toplwad_aero(klon), sollwad_aero(klon))
831      ALLOCATE(toplwai_aero(klon), sollwai_aero(klon))
832      ALLOCATE(toplwad0_aero(klon), sollwad0_aero(klon))
833      ALLOCATE(topsw_aero(klon,naero_grp), solsw_aero(klon,naero_grp))
834      ALLOCATE(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp))
835      ALLOCATE(topswcf_aero(klon,3), solswcf_aero(klon,3))
836      ALLOCATE(du_gwd_hines(klon,klev),dv_gwd_hines(klon,klev))
837      ALLOCATE(dv_gwd_rando(klon,klev),dv_gwd_front(klon,klev))
838      ALLOCATE(east_gwstress(klon,klev),west_gwstress(klon,klev))
839      east_gwstress(:,:)=0 !ym missing init
840      west_gwstress(:,:)=0 !ym missing init
841      ALLOCATE(d_t_hin(klon,klev))
842      ALLOCATE(d_q_ch4(klon,klev))
843      ALLOCATE(stratomask(klon,klev))
844      ALLOCATE(tausum_aero(klon,nwave,naero_tot))
845      ALLOCATE(drytausum_aero(klon,naero_tot))
846      ALLOCATE(tau3d_aero(klon,klev,nwave,naero_tot))
847      ALLOCATE(scdnc(klon, klev))
848      ALLOCATE(dNovrN(klon))
849      ALLOCATE(cldncl(klon))
850      ALLOCATE(reffclwtop(klon))
851      ALLOCATE(lcc(klon))
852      ALLOCATE(reffclws(klon, klev))
853      ALLOCATE(reffclwc(klon, klev))
854      ALLOCATE(cldnvi(klon))
855      ALLOCATE(lcc3d(klon, klev))
856      ALLOCATE(lcc3dcon(klon, klev))
857      ALLOCATE(lcc3dstra(klon, klev))
858      ALLOCATE(icc3dcon(klon, klev))
859      ALLOCATE(icc3dstra(klon, klev))
860      ALLOCATE(zfice(klon, klev))
861      ALLOCATE(od443aer(klon))
862      ALLOCATE(od550aer(klon))
863      ALLOCATE(od865aer(klon))
864      ALLOCATE(dryod550aer(klon))
865      dryod550aer(:) = 0.
866      ALLOCATE(abs550aer(klon))
867      abs550aer(:) = 0.
868      ALLOCATE(ec550aer(klon,klev))
869      ALLOCATE(od550lt1aer(klon))
870      ALLOCATE(sconcso4(klon))
871      ALLOCATE(sconcno3(klon))
872      ALLOCATE(sconcoa(klon))
873      ALLOCATE(sconcbc(klon))
874      ALLOCATE(sconcss(klon))
875      ALLOCATE(sconcdust(klon))
876      ALLOCATE(concso4(klon,klev))
877      ALLOCATE(concno3(klon,klev))
878      ALLOCATE(concoa(klon,klev))
879      ALLOCATE(concbc(klon,klev))
880      ALLOCATE(concss(klon,klev))
881      ALLOCATE(concdust(klon,klev))
882      ALLOCATE(loadso4(klon))
883      ALLOCATE(loadoa(klon))
884      ALLOCATE(loadbc(klon))
885      ALLOCATE(loadss(klon))
886      ALLOCATE(loaddust(klon))
887      ALLOCATE(loadno3(klon))
888      ALLOCATE(load_tmp1(klon))
889      ALLOCATE(load_tmp2(klon))
890      ALLOCATE(load_tmp3(klon))
891      ALLOCATE(load_tmp4(klon))
892      ALLOCATE(load_tmp5(klon))
893      ALLOCATE(load_tmp6(klon))
894      ALLOCATE(load_tmp7(klon))
895      ALLOCATE(load_tmp8(klon))
896      ALLOCATE(load_tmp9(klon))
897      ALLOCATE(load_tmp10(klon))
898
899!IM ajout variables CFMIP2/CMIP5
900      ALLOCATE(topswad_aerop(klon), solswad_aerop(klon))
901      ALLOCATE(topswai_aerop(klon), solswai_aerop(klon))
902      ALLOCATE(topswad0_aerop(klon), solswad0_aerop(klon))
903      ALLOCATE(topsw_aerop(klon,naero_grp), topsw0_aerop(klon,naero_grp))
904      ALLOCATE(solsw_aerop(klon,naero_grp), solsw0_aerop(klon,naero_grp))
905      ALLOCATE(topswcf_aerop(klon,naero_grp), solswcf_aerop(klon,naero_grp))
906
907! additional LW variables CK
908      ALLOCATE(toplwad_aerop(klon), sollwad_aerop(klon))
909      ALLOCATE(toplwai_aerop(klon), sollwai_aerop(klon))
910      ALLOCATE(toplwad0_aerop(klon), sollwad0_aerop(klon))
911
912! FH Ajout de celles necessaires au phys_output_write_mod
913
914      ALLOCATE(tal1(klon), pal1(klon), pab1(klon), pab2(klon))
915      ALLOCATE(ptstar(klon),pt0(klon),slp(klon))
916!!
917!!          Wake variables
918      ALLOCATE(alp_wake(klon))
919      ALLOCATE(wake_h(klon),wake_k(klon))
920      ALLOCATE(wake_omg(klon, klev))
921      ALLOCATE(d_deltat_wk(klon, klev), d_deltaq_wk(klon, klev))
922      ALLOCATE(d_s_wk(klon), d_dens_a_wk(klon), d_dens_wk(klon))
923      ALLOCATE(d_deltat_wk_gw(klon, klev), d_deltaq_wk_gw(klon, klev))
924      ALLOCATE(d_deltat_vdf(klon, klev), d_deltaq_vdf(klon, klev))
925!!      ALLOCATE( d_s_vdf(klon), d_dens_vdf(klon))
926      ALLOCATE(d_deltat_the(klon, klev), d_deltaq_the(klon, klev))
927!!      ALLOCATE( d_s_the(klon), d_dens_the(klon))
928      ALLOCATE(d_deltat_ajs_cv(klon, klev), d_deltaq_ajs_cv(klon, klev))
929#ifdef ISO
930      ALLOCATE(d_deltaxt_wk(ntraciso,klon, klev))
931      ALLOCATE(d_deltaxt_wk_gw(ntraciso,klon, klev))
932      ALLOCATE(d_deltaxt_the(ntraciso,klon, klev))
933      ALLOCATE(d_deltaxt_vdf(ntraciso,klon, klev))
934      ALLOCATE(d_deltaxt_ajs_cv(ntraciso,klon, klev))
935#endif
936!!         End of wake variables
937!!
938      ALLOCATE(bils(klon))
939      ALLOCATE(cdragm(klon), cdragh(klon), cldh(klon), cldl(klon))
940      ALLOCATE(cldm(klon), cldq(klon), cldt(klon), qsat2m(klon))
941      ALLOCATE(JrNt(klon))
[4523]942      ALLOCATE(dthmin(klon), evap(klon), snowerosion(klon), fder(klon), plcl(klon), plfc(klon))
943      ALLOCATE(prw(klon), prlw(klon), prsw(klon), prbsw(klon), zustar(klon), zu10m(klon), zv10m(klon), rh2m(klon))
[3927]944      ALLOCATE(s_lcl(klon))
945      ALLOCATE(s_pblh(klon), s_pblt(klon), s_therm(klon))
946!
947!nrlmd+jyg<
948      ALLOCATE(s_pblh_x(klon), s_pblh_w(klon))
949      ALLOCATE(s_lcl_x(klon), s_lcl_w(klon))
950!>nrlmd+jyg
951!
952      ALLOCATE(slab_wfbils(klon), tpot(klon), tpote(klon), ue(klon))
953      ALLOCATE(uq(klon), ve(klon), vq(klon), zxffonte(klon))
954      ALLOCATE(uwat(klon), vwat(klon))
955      ALLOCATE(zxfqcalving(klon), zxfluxlat(klon))
956      ALLOCATE(zxtsol(klon), snow_lsc(klon), zxfqfonte(klon), zxqsurf(klon))
957      ALLOCATE(zxrunofflic(klon))
958      ALLOCATE(rain_lsc(klon))
959      ALLOCATE(rain_num(klon))
[4118]960      ALLOCATE(qlth(klon,klev), qith(klon,klev), qsith(klon,klev), wiceth(klon,klev))
[4003]961      !
[3927]962#ifdef ISO
963      ALLOCATE(xtevap(ntraciso,klon))
964      ALLOCATE(xtprw(ntraciso,klon))
965      ALLOCATE(zxfxtcalving(niso,klon))
966      ALLOCATE(xtsnow_lsc(ntraciso,klon), zxfxtfonte(niso,klon))
967      ALLOCATE(zxxtrunofflic(niso,klon))
968      ALLOCATE(xtrain_lsc(ntraciso,klon))
969      ALLOCATE(xtrunoff_diag(niso,klon))
970      ALLOCATE(h1_diag(klon),runoff_diag(klon))
971#endif
972!
973      ALLOCATE(sens_x(klon), sens_w(klon))
974      ALLOCATE(zxfluxlat_x(klon), zxfluxlat_w(klon))
[3940]975      ALLOCATE(delta_qsurf(klon))
[3927]976!jyg<
977!!      ALLOCATE(t_x(klon,klev), t_w(klon,klev))
978!!      ALLOCATE(q_x(klon,klev), q_w(klon,klev))
979!>jyg
980      ALLOCATE(dtvdf_x(klon,klev), dtvdf_w(klon,klev))
981      dtvdf_x = 0 ; dtvdf_w=0 ;   !ym missing init
982      ALLOCATE(dqvdf_x(klon,klev), dqvdf_w(klon,klev))
983      dqvdf_x = 0 ; dqvdf_w=0 ;   !ym missing init
984      ALLOCATE(pbl_tke_input(klon,klev+1,nbsrf))
985      ALLOCATE(t_therm(klon,klev), q_therm(klon,klev),u_therm(klon,klev), v_therm(klon,klev))
986      ALLOCATE(cdragh_x(klon), cdragh_w(klon))
987      ALLOCATE(cdragm_x(klon), cdragm_w(klon))
988      ALLOCATE(kh(klon), kh_x(klon), kh_w(klon))
989#ifdef ISO
990      ALLOCATE(dxtvdf_x(ntraciso,klon,klev), dxtvdf_w(ntraciso,klon,klev))
991      ALLOCATE(xt_therm(ntraciso,klon,klev))
992#endif
993!
994      ALLOCATE(ptconv(klon,klev))
995      ALLOCATE(ratqsc(klon,klev))
996!
997      ALLOCATE(wbeff(klon), convoccur(klon), zmax_th(klon))
998      ALLOCATE(zq2m(klon), zt2m(klon), weak_inversion(klon))
999      ALLOCATE(zt2m_min_mon(klon), zt2m_max_mon(klon))
1000      ALLOCATE(t2m_min_mon(klon), t2m_max_mon(klon))
1001      ALLOCATE(sens(klon), flwp(klon), fiwp(klon))
1002      ALLOCATE(alp_bl_conv(klon), alp_bl_det(klon))
1003      alp_bl_conv(:)=0 ;  alp_bl_det(:)=0
1004      ALLOCATE(alp_bl_fluct_m(klon), alp_bl_fluct_tke(klon))
1005      alp_bl_fluct_m(:)=0 ; alp_bl_fluct_tke(:)= 0.
1006      ALLOCATE(alp_bl_stat(klon), n2(klon), s2(klon))
1007      alp_bl_stat(:)=0
1008      ALLOCATE(proba_notrig(klon), random_notrig(klon))
1009
1010      ALLOCATE(dnwd0(klon, klev))
1011!      ALLOCATE(upwd(klon, klev), omega(klon, klev), coefh(klon, klev))
1012      ALLOCATE(omega(klon, klev))
1013      ALLOCATE(epmax_diag(klon)) ! epmax_cape
1014!      ALLOCATE(coefm(klon, klev), lambda_th(klon, klev), cldemi(klon, klev))
1015      ALLOCATE(lambda_th(klon, klev), cldemi(klon, klev))
1016      ALLOCATE(cldfra(klon, klev), cldtau(klon, klev), fiwc(klon, klev))
1017      ALLOCATE(fl(klon, klev), re(klon, klev), flwc(klon, klev))
1018      ALLOCATE(ref_liq(klon, klev), ref_ice(klon, klev), theta(klon, klev))
1019      ALLOCATE(ref_liq_pi(klon, klev), ref_ice_pi(klon, klev))
[3940]1020      ALLOCATE(zphi(klon, klev), zx_rh(klon, klev), zx_rhl(klon,klev), zx_rhi(klon,klev))
[4072]1021      zx_rhl(:,:)=0.; zx_rhi(:,:)=0. ! because not always defined
[3927]1022      ALLOCATE(pmfd(klon, klev), pmfu(klon, klev))
1023
1024      ALLOCATE(t2m(klon, nbsrf), fluxlat(klon, nbsrf))
1025      ALLOCATE(fsollw(klon, nbsrf))
1026      ALLOCATE(fsolsw(klon, nbsrf), wfbils(klon, nbsrf), wfbilo(klon, nbsrf))
1027      ALLOCATE(wfevap(klon, nbsrf), wfrain(klon,nbsrf), wfsnow(klon, nbsrf))
1028      ALLOCATE(evap_pot(klon, nbsrf))
[4576]1029! FC
1030      ALLOCATE(zxfluxq(klon,klev),zxfluxt(klon,klev))
[3927]1031!
[4576]1032!
[3927]1033!  Deep convective variables used in phytrac
1034      ALLOCATE(pmflxr(klon, klev+1), pmflxs(klon, klev+1))
1035      ALLOCATE(wdtrainA(klon,klev),wdtrainS(klon,klev),wdtrainM(klon,klev))
1036      ALLOCATE(dnwd(klon, klev), upwd(klon, klev) )
1037      ALLOCATE(ep(klon,klev))                          ! epmax_cape
1038      ALLOCATE(da(klon,klev), mp(klon,klev) )
1039      ALLOCATE(phi(klon,klev,klev) )
1040      ALLOCATE(wght_cvfd(klon,klev) )
1041      ALLOCATE(phi2(klon,klev,klev) )
1042      ALLOCATE(d1a(klon,klev), dam(klon,klev))
1043      ALLOCATE(ev(klon,klev) )
1044      ALLOCATE(elij(klon,klev,klev) )
1045      ALLOCATE(qtaa(klon,klev) )
1046      ALLOCATE(clw(klon,klev) )
1047      ALLOCATE(epmlmMm(klon,klev,klev), eplaMm(klon,klev) )
1048      ALLOCATE(sij(klon,klev,klev) )
1049#ifdef ISO
1050      ALLOCATE(xtwdtrainA(ntraciso,klon,klev))
1051      ALLOCATE(xtev(ntraciso,klon,klev) )
1052      ALLOCATE(xttaa(ntraciso,klon,klev) )
1053      ALLOCATE(xtclw(ntraciso,klon,klev) )
1054#ifdef DIAGISO
1055      ALLOCATE(qlp(klon,klev))
1056      ALLOCATE(qvp(klon,klev))
1057      ALLOCATE(fq_detrainement(klon,klev))
1058      ALLOCATE(fq_ddft(klon,klev))
1059      ALLOCATE(fq_fluxmasse(klon,klev))
1060      ALLOCATE(fq_evapprecip(klon,klev))
1061      ALLOCATE(f_detrainement(klon,klev), q_detrainement(klon,klev))
1062      ALLOCATE(xtlp(ntraciso,klon,klev))
1063      ALLOCATE(xtvp(ntraciso,klon,klev))
1064      ALLOCATE(q_the(klon,klev), xt_the(ntraciso,klon,klev))
1065      ALLOCATE(fxt_detrainement(ntraciso,klon,klev))
1066      ALLOCATE(fxt_ddft(ntraciso,klon,klev))
1067      ALLOCATE(fxt_fluxmasse(ntraciso,klon,klev))
1068      ALLOCATE(fxt_evapprecip(ntraciso,klon,klev))
1069      ALLOCATE(xt_detrainement(ntraciso,klon,klev))
1070#endif
1071#endif
1072
1073      ALLOCATE(prfl(klon, klev+1))
1074      ALLOCATE(psfl(klon, klev+1), fraca(klon, klev+1), Vprecip(klon, klev+1))
[4523]1075      ALLOCATE(bsfl(klon,klev+1))
[3927]1076      ALLOCATE(zw2(klon, klev+1))
1077
1078      ALLOCATE(fluxu(klon, klev, nbsrf), fluxv(klon, klev, nbsrf))
1079      ALLOCATE(fluxt(klon, klev, nbsrf))
1080
1081      ALLOCATE(uwriteSTD(klon,nlevSTD,nfiles), vwriteSTD(klon,nlevSTD,nfiles))
1082      ALLOCATE(wwriteSTD(klon,nlevSTD,nfiles), phiwriteSTD(klon,nlevSTD,nfiles))
1083      ALLOCATE(qwriteSTD(klon,nlevSTD,nfiles), twriteSTD(klon,nlevSTD,nfiles))
1084      ALLOCATE(rhwriteSTD(klon,nlevSTD,nfiles))
1085
1086! ug et d'autres encore:
1087      ALLOCATE(beta_prec(klon,klev))
1088      ALLOCATE(rneb(klon,klev),rnebjn(klon,klev),rneblsvol(klon,klev))
[4530]1089      ALLOCATE(pfraclr(klon,klev),pfracld(klon,klev))
1090      pfraclr(:,:)=0. ; pfracld(:,:)=0. ! because not always defined
[4562]1091      ALLOCATE(distcltop(klon,klev))
[3927]1092
1093      ALLOCATE (zxsnow(klon),snowhgt(klon),qsnow(klon),to_ice(klon))
1094      ALLOCATE (sissnow(klon),runoff(klon),albsol3_lic(klon))
1095#ifdef ISO
1096      ALLOCATE (zxxtsnow(niso,klon))
1097      ALLOCATE(xtVprecip(ntraciso,klon, klev+1),xtVprecipi(ntraciso,klon, klev+1))
1098      ALLOCATE(pxtsfl(ntraciso,klon, klev+1),pxtrfl(ntraciso,klon, klev+1))
1099#endif
1100
1101      ALLOCATE (p_tropopause(klon))
1102      ALLOCATE (z_tropopause(klon))
1103      ALLOCATE (t_tropopause(klon))
1104
[3940]1105      ALLOCATE(zn2mout(klon,6))
1106
[4065]1107! Supersaturation
1108      ALLOCATE(rneb_seri(klon,klev))
1109      ALLOCATE(d_rneb_dyn(klon,klev))
1110      ALLOCATE(qclr(klon,klev), qcld(klon,klev), qss(klon,klev), qvc(klon,klev))
1111      ALLOCATE(rnebclr(klon,klev), rnebss(klon,klev), gamma_ss(klon,klev))
1112      ALLOCATE(N1_ss(klon,klev), N2_ss(klon,klev))
1113      ALLOCATE(drneb_sub(klon,klev), drneb_con(klon,klev), drneb_tur(klon,klev), drneb_avi(klon,klev))
1114      ALLOCATE(zqsatl(klon,klev), zqsats(klon,klev))
1115      ALLOCATE(Tcontr(klon,klev), qcontr(klon,klev), qcontr2(klon,klev), fcontrN(klon,klev), fcontrP(klon,klev))
1116
[3927]1117#ifdef CPP_StratAer
1118      ALLOCATE (R2SO4(klon,klev))
1119      ALLOCATE (DENSO4(klon,klev))
1120      ALLOCATE (f_r_wet(klon,klev))
1121      ALLOCATE (decfluxaer(klon,nbtr))
1122      ALLOCATE (mdw(nbtr))
1123      ALLOCATE (budg_3D_nucl(klon,klev))
1124      ALLOCATE (budg_3D_cond_evap(klon,klev))
1125      ALLOCATE (budg_3D_ocs_to_so2(klon,klev))
1126      ALLOCATE (budg_3D_so2_to_h2so4(klon,klev))
1127      ALLOCATE (budg_3D_backgr_ocs(klon,klev))
1128      ALLOCATE (budg_3D_backgr_so2(klon,klev))
1129      ALLOCATE (OCS_lifetime(klon,klev))
1130      ALLOCATE (SO2_lifetime(klon,klev))
1131      ALLOCATE (alpha_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave,nbtr))
1132      ALLOCATE (piz_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave,nbtr))
1133      ALLOCATE (cg_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave,nbtr))
1134      ALLOCATE (tau_strat_550(klon,klev))
1135      ALLOCATE (tau_strat_1020(klon,klev))
1136      ALLOCATE (tausum_strat(klon,3))
1137      ALLOCATE (budg_dep_dry_ocs(klon))
1138      ALLOCATE (budg_dep_wet_ocs(klon))
1139      ALLOCATE (budg_dep_dry_so2(klon))
1140      ALLOCATE (budg_dep_wet_so2(klon))
1141      ALLOCATE (budg_dep_dry_h2so4(klon))
1142      ALLOCATE (budg_dep_wet_h2so4(klon))
1143      ALLOCATE (budg_dep_dry_part(klon))
1144      ALLOCATE (budg_dep_wet_part(klon))
1145      ALLOCATE (budg_emi_ocs(klon))
1146      ALLOCATE (budg_emi_so2(klon))
1147      ALLOCATE (budg_emi_h2so4(klon))
1148      ALLOCATE (budg_emi_part(klon))
1149      ALLOCATE (budg_ocs_to_so2(klon))
1150      ALLOCATE (budg_so2_to_h2so4(klon))
1151      ALLOCATE (budg_h2so4_to_part(klon))
1152      ALLOCATE (budg_sed_part(klon))
1153      ALLOCATE (surf_PM25_sulf(klon))
1154      ALLOCATE (vsed_aer(klon,klev))
1155#endif
1156
1157END SUBROUTINE phys_local_var_init
1158
1159!======================================================================
1160SUBROUTINE phys_local_var_end
1161USE dimphy
1162USE indice_sol_mod
1163IMPLICIT NONE
[4523]1164      DEALLOCATE(t_seri,q_seri,ql_seri,qs_seri, qbs_seri)
[3927]1165      DEALLOCATE(u_seri,v_seri)
[4040]1166      DEALLOCATE(l_mixmin,l_mix, tke_dissip,wprime)
[3927]1167
1168      DEALLOCATE(tr_seri)
1169      DEALLOCATE(d_t_dyn,d_q_dyn)
[4523]1170      DEALLOCATE(d_ql_dyn,d_qs_dyn, d_qbs_dyn)
1171      DEALLOCATE(d_q_dyn2d,d_ql_dyn2d,d_qs_dyn2d, d_qbs_dyn2d)
[3927]1172      DEALLOCATE(d_u_dyn,d_v_dyn)
1173      DEALLOCATE(d_tr_dyn)                      !RomP
1174      DEALLOCATE(d_t_con,d_q_con)
1175      DEALLOCATE(d_u_con,d_v_con)
1176      DEALLOCATE(d_t_wake,d_q_wake)
1177      DEALLOCATE(d_t_lsc,d_q_lsc)
1178      DEALLOCATE(d_t_lwr,d_t_lw0)
1179      DEALLOCATE(d_t_swr,d_t_sw0)
1180      DEALLOCATE(d_ql_lsc,d_qi_lsc)
1181      DEALLOCATE(d_t_ajsb,d_q_ajsb)
1182      DEALLOCATE(d_t_ajs,d_q_ajs)
1183!nrlmd<
1184      DEALLOCATE(d_t_ajs_w,d_q_ajs_w)
1185      DEALLOCATE(d_t_ajs_x,d_q_ajs_x)
1186!>nrlmd
1187      DEALLOCATE(d_u_ajs,d_v_ajs)
1188      DEALLOCATE(d_t_eva,d_q_eva)
1189      DEALLOCATE(d_ql_eva,d_qi_eva)
1190      DEALLOCATE(d_t_lscst,d_q_lscst)
1191      DEALLOCATE(d_t_lscth,d_q_lscth)
1192      DEALLOCATE(plul_st,plul_th)
1193      DEALLOCATE(d_t_vdf,d_q_vdf,d_t_diss)
[4523]1194      DEALLOCATE(d_qbs_vdf)
1195      DEALLOCATE(d_t_bs,d_q_bs,d_qbs_bs)
[3927]1196#ifdef ISO
1197      deallocate(xt_seri,xtl_seri,xts_seri)
1198      DEALLOCATE(d_xtl_eva,d_xti_eva)
1199      deallocate(d_xt_dyn,d_xtl_dyn,d_xts_dyn)
1200      deallocate(d_xt_con)
1201      deallocate(d_xt_wake)
1202      deallocate(d_xt_lsc)
1203      deallocate(d_xtl_lsc,d_xti_lsc)
1204      deallocate(d_xt_ajsb)
1205      deallocate(d_xt_ajs)
1206      deallocate(d_xt_ajs_w,d_xt_ajs_x)
1207      deallocate(d_xt_eva)
1208      deallocate(d_xtl_eva)
1209      deallocate(d_xti_eva)
1210      deallocate(d_xt_vdf)
1211      deallocate(d_xt_vdf_w,d_xt_vdf_x)
1212      deallocate(d_xt_ch4)
1213      deallocate(d_xt_prod_nucl)
1214      deallocate(d_xt_cosmo)
1215      deallocate(d_xt_decroiss)
1216#endif
1217
1218      DEALLOCATE(d_u_vdf,d_v_vdf)
1219      DEALLOCATE(d_t_oli,d_t_oro)
1220      DEALLOCATE(d_u_oli,d_v_oli)
1221      DEALLOCATE(d_u_oro,d_v_oro)
1222      DEALLOCATE(d_t_oro_gw)
1223      DEALLOCATE(d_u_oro_gw,d_v_oro_gw)
1224      DEALLOCATE(d_t_lif,d_t_ec)
1225      DEALLOCATE(d_u_lif,d_v_lif)
1226      DEALLOCATE(d_ts, d_tr)
1227      DEALLOCATE(topswad_aero,solswad_aero)
1228      DEALLOCATE(topswai_aero,solswai_aero)
1229      DEALLOCATE(topswad0_aero,solswad0_aero)
1230      DEALLOCATE(toplwad_aero,sollwad_aero)
1231      DEALLOCATE(toplwai_aero,sollwai_aero)
1232      DEALLOCATE(toplwad0_aero,sollwad0_aero)
1233      DEALLOCATE(topsw_aero,solsw_aero)
1234      DEALLOCATE(topsw0_aero,solsw0_aero)
1235      DEALLOCATE(topswcf_aero,solswcf_aero)
1236      DEALLOCATE(stratomask)
1237      DEALLOCATE(tausum_aero)
1238      DEALLOCATE(drytausum_aero)
1239      DEALLOCATE(tau3d_aero)
1240      DEALLOCATE(scdnc)
1241      DEALLOCATE(dNovrN)
1242      DEALLOCATE(cldncl)
1243      DEALLOCATE(reffclwtop)
1244      DEALLOCATE(lcc)
1245      DEALLOCATE(reffclws)
1246      DEALLOCATE(reffclwc)
1247      DEALLOCATE(cldnvi)
1248      DEALLOCATE(lcc3d)
1249      DEALLOCATE(lcc3dcon)
1250      DEALLOCATE(lcc3dstra)
1251      DEALLOCATE(icc3dcon)
1252      DEALLOCATE(icc3dstra)
1253      DEALLOCATE(zfice)
1254      DEALLOCATE(od443aer)
1255      DEALLOCATE(od550aer)
1256      DEALLOCATE(od865aer)
1257      DEALLOCATE(dryod550aer)
1258      DEALLOCATE(abs550aer)
1259      DEALLOCATE(ec550aer)
1260      DEALLOCATE(od550lt1aer)
1261      DEALLOCATE(sconcso4)
1262      DEALLOCATE(sconcno3)
1263      DEALLOCATE(sconcoa)
1264      DEALLOCATE(sconcbc)
1265      DEALLOCATE(sconcss)
1266      DEALLOCATE(sconcdust)
1267      DEALLOCATE(concso4)
1268      DEALLOCATE(concno3)
1269      DEALLOCATE(concoa)
1270      DEALLOCATE(concbc)
1271      DEALLOCATE(concss)
1272      DEALLOCATE(concdust)
1273      DEALLOCATE(loadso4)
1274      DEALLOCATE(loadoa)
1275      DEALLOCATE(loadbc)
1276      DEALLOCATE(loadss)
1277      DEALLOCATE(loaddust)
1278      DEALLOCATE(loadno3)
1279      DEALLOCATE(load_tmp1)
1280      DEALLOCATE(load_tmp2)
1281      DEALLOCATE(load_tmp3)
1282      DEALLOCATE(load_tmp4)
1283      DEALLOCATE(load_tmp5)
1284      DEALLOCATE(load_tmp6)
1285      DEALLOCATE(load_tmp7)
1286      DEALLOCATE(load_tmp8)
1287      DEALLOCATE(load_tmp9)
1288      DEALLOCATE(load_tmp10)
1289      DEALLOCATE(du_gwd_hines,dv_gwd_hines,d_t_hin)
1290      DEALLOCATE(d_q_ch4)
1291      DEALLOCATE(dv_gwd_rando,dv_gwd_front)
1292      DEALLOCATE(east_gwstress,west_gwstress)
1293
1294!IM ajout variables CFMIP2/CMIP5
1295      DEALLOCATE(topswad_aerop, solswad_aerop)
1296      DEALLOCATE(topswai_aerop, solswai_aerop)
1297      DEALLOCATE(topswad0_aerop, solswad0_aerop)
1298      DEALLOCATE(topsw_aerop, topsw0_aerop)
1299      DEALLOCATE(solsw_aerop, solsw0_aerop)
1300      DEALLOCATE(topswcf_aerop, solswcf_aerop)
1301
1302!CK LW diagnostics
1303      DEALLOCATE(toplwad_aerop, sollwad_aerop)
1304      DEALLOCATE(toplwai_aerop, sollwai_aerop)
1305      DEALLOCATE(toplwad0_aerop, sollwad0_aerop)
1306
1307! FH Ajout de celles necessaires au phys_output_write_mod
1308      DEALLOCATE(tal1, pal1, pab1, pab2)
1309      DEALLOCATE(ptstar, pt0, slp)
1310!
1311      DEALLOCATE(alp_wake)
1312      DEALLOCATE(wake_h,wake_k)
1313      DEALLOCATE(wake_omg)
1314      DEALLOCATE(d_deltat_wk, d_deltaq_wk)
1315      DEALLOCATE(d_s_wk, d_dens_a_wk, d_dens_wk)
1316      DEALLOCATE(d_deltat_wk_gw, d_deltaq_wk_gw)
1317      DEALLOCATE(d_deltat_vdf, d_deltaq_vdf)
1318!!      DEALLOCATE( d_s_vdf, d_dens_vdf)
1319      DEALLOCATE(d_deltat_the, d_deltaq_the)
1320!!      DEALLOCATE( d_s_the, d_dens_the)
1321      DEALLOCATE(d_deltat_ajs_cv, d_deltaq_ajs_cv)
1322#ifdef ISO
1323      DEALLOCATE(d_deltaxt_wk)
1324      DEALLOCATE(d_deltaxt_wk_gw)
1325      DEALLOCATE(d_deltaxt_ajs_cv)
1326      DEALLOCATE(d_deltaxt_vdf)
1327#endif
1328!
1329      DEALLOCATE(bils)
1330      DEALLOCATE(cdragm, cdragh, cldh, cldl)
1331      DEALLOCATE(cldm, cldq, cldt, qsat2m)
[3940]1332      DEALLOCATE(JrNt)
[4523]1333      DEALLOCATE(dthmin, evap, snowerosion, fder, plcl, plfc)
1334      DEALLOCATE(prw, prlw, prsw, prbsw, zustar, zu10m, zv10m, rh2m, s_lcl)
[3927]1335      DEALLOCATE(s_pblh, s_pblt, s_therm)
1336!
1337!nrlmd+jyg<
1338      DEALLOCATE(s_pblh_x, s_pblh_w)
1339      DEALLOCATE(s_lcl_x, s_lcl_w)
1340!>nrlmd+jyg
1341!
1342      DEALLOCATE(slab_wfbils, tpot, tpote, ue)
1343      DEALLOCATE(uq, ve, vq, zxffonte)
1344      DEALLOCATE(uwat, vwat)
1345      DEALLOCATE(zxfqcalving, zxfluxlat)
1346      DEALLOCATE(zxrunofflic)
[4523]1347      DEALLOCATE(zxustartlic, zxrhoslic)
[3927]1348      DEALLOCATE(zxtsol, snow_lsc, zxfqfonte, zxqsurf)
1349      DEALLOCATE(rain_lsc)
1350      DEALLOCATE(rain_num)
[4118]1351      DEALLOCATE(qlth, qith, qsith, wiceth)
[3927]1352!
1353      DEALLOCATE(sens_x, sens_w)
1354      DEALLOCATE(zxfluxlat_x, zxfluxlat_w)
[3940]1355      DEALLOCATE(delta_qsurf)
[3927]1356!jyg<
1357!!      DEALLOCATE(t_x, t_w)
1358!!      DEALLOCATE(q_x, q_w)
1359!>jyg
1360      DEALLOCATE(d_t_vdf_x, d_t_vdf_w)
1361      DEALLOCATE(d_q_vdf_x, d_q_vdf_w)
1362      DEALLOCATE(pbl_tke_input)
1363      DEALLOCATE(t_therm, q_therm, u_therm, v_therm)
1364      DEALLOCATE(cdragh_x, cdragh_w)
1365      DEALLOCATE(cdragm_x, cdragm_w)
1366      DEALLOCATE(kh, kh_x, kh_w)
1367#ifdef ISO
1368      DEALLOCATE(xtevap,xtprw)
1369      DEALLOCATE(zxfxtcalving)
1370      DEALLOCATE(zxxtrunofflic)
1371      DEALLOCATE(xtsnow_lsc, zxfxtfonte)
1372      DEALLOCATE(xtrain_lsc)
1373      DEALLOCATE(dxtvdf_x, dxtvdf_w)
1374      DEALLOCATE(xt_therm)
1375      DEALLOCATE(h1_diag,runoff_diag,xtrunoff_diag)
1376#endif
1377!
1378      DEALLOCATE(ptconv)
1379      DEALLOCATE(ratqsc)
1380!
1381      DEALLOCATE(wbeff, convoccur, zmax_th)
1382      DEALLOCATE(zq2m, zt2m, weak_inversion)
1383      DEALLOCATE(zt2m_min_mon, zt2m_max_mon)
1384      DEALLOCATE(t2m_min_mon, t2m_max_mon)
1385      DEALLOCATE(sens, flwp, fiwp)
1386      DEALLOCATE(alp_bl_conv,alp_bl_det)
1387      DEALLOCATE(alp_bl_fluct_m,alp_bl_fluct_tke)
1388      DEALLOCATE(alp_bl_stat, n2, s2)
1389      DEALLOCATE(proba_notrig, random_notrig)
[4576]1390!FC
1391      DEALLOCATE(zxfluxq,zxfluxt)
[3927]1392
1393      DEALLOCATE(dnwd0)
1394!      DEALLOCATE(upwd, omega, coefh)
1395      DEALLOCATE(omega)
1396      DEALLOCATE(epmax_diag)
1397!      DEALLOCATE(coefm, lambda_th, cldemi)
1398      DEALLOCATE(lambda_th, cldemi)
1399      DEALLOCATE(cldfra, cldtau, fiwc)
1400      DEALLOCATE(fl, re, flwc)
1401      DEALLOCATE(ref_liq, ref_ice, theta)
1402      DEALLOCATE(ref_liq_pi, ref_ice_pi)
[3940]1403      DEALLOCATE(zphi, zx_rh, zx_rhl, zx_rhi)
[3927]1404      DEALLOCATE(pmfd, pmfu)
1405
1406      DEALLOCATE(t2m, fluxlat)
1407      DEALLOCATE(fsollw, evap_pot)
1408      DEALLOCATE(fsolsw, wfbils, wfbilo)
1409      DEALLOCATE(wfevap,wfrain,wfsnow)
1410
1411      DEALLOCATE(pmflxr, pmflxs)
1412      DEALLOCATE(wdtrainA, wdtrainS, wdtrainM)
1413      DEALLOCATE(upwd, dnwd)
1414      DEALLOCATE(ep)
1415      DEALLOCATE(da, mp )
1416      DEALLOCATE(phi )
1417      DEALLOCATE(wght_cvfd )
1418      DEALLOCATE(phi2 )
1419      DEALLOCATE(d1a, dam)
1420      DEALLOCATE(ev )
1421      DEALLOCATE(elij )
1422      DEALLOCATE(qtaa )
1423      DEALLOCATE(clw )
1424      DEALLOCATE(epmlmMm, eplaMm )
1425      DEALLOCATE(sij )
1426#ifdef ISO
1427      DEALLOCATE(xtwdtrainA)
1428      DEALLOCATE(xttaa )
1429      DEALLOCATE(xtclw )
1430      DEALLOCATE(xtev )
1431#ifdef DIAGISO
1432      DEALLOCATE(qlp)
1433      DEALLOCATE(qvp)
1434      DEALLOCATE(fq_detrainement)
1435      DEALLOCATE(fq_ddft)
1436      DEALLOCATE(fq_fluxmasse)
1437      DEALLOCATE(fq_evapprecip)
1438      DEALLOCATE(f_detrainement,q_detrainement)
1439      DEALLOCATE(xtlp)
1440      DEALLOCATE(xtvp)
1441      DEALLOCATE(q_the,xt_the)
1442      DEALLOCATE(fxt_detrainement)
1443      DEALLOCATE(fxt_ddft)
1444      DEALLOCATE(fxt_fluxmasse)
1445      DEALLOCATE(fxt_evapprecip)
1446      DEALLOCATE(xt_detrainement)
1447#endif
1448#endif
1449
1450
[4523]1451      DEALLOCATE(prfl, psfl, bsfl, fraca, Vprecip)
[3927]1452      DEALLOCATE(zw2)
1453
1454      DEALLOCATE(fluxu, fluxv)
1455      DEALLOCATE(fluxt)
1456
1457      DEALLOCATE(uwriteSTD, vwriteSTD)
1458      DEALLOCATE(wwriteSTD, phiwriteSTD)
1459      DEALLOCATE(qwriteSTD, twriteSTD, rhwriteSTD)
1460
1461! ug et d'autres encore:
1462      DEALLOCATE(beta_prec)
1463      DEALLOCATE(rneb)
[4530]1464      DEALLOCATE(pfraclr,pfracld)
[4562]1465      DEALLOCATE(distcltop)
[3927]1466      DEALLOCATE (zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic)
1467#ifdef ISO
1468      DEALLOCATE (zxxtsnow,xtVprecip,xtVprecipi,pxtrfl,pxtsfl)
1469#endif
1470
1471      DEALLOCATE (p_tropopause)
1472      DEALLOCATE (z_tropopause)
1473      DEALLOCATE (t_tropopause)
[3940]1474      DEALLOCATE(zn2mout)
[3927]1475
[4065]1476! Supersaturation
1477      DEALLOCATE(rneb_seri)
1478      DEALLOCATE(d_rneb_dyn)
1479      DEALLOCATE(qclr, qcld, qss, qvc)
1480      DEALLOCATE(rnebclr, rnebss, gamma_ss)
1481      DEALLOCATE(N1_ss, N2_ss)
1482      DEALLOCATE(drneb_sub, drneb_con, drneb_tur, drneb_avi)
1483      DEALLOCATE(zqsatl, zqsats)
1484      DEALLOCATE(Tcontr, qcontr, qcontr2, fcontrN, fcontrP)
1485
[3927]1486#ifdef CPP_StratAer
1487! variables for strat. aerosol CK
1488      DEALLOCATE (R2SO4)
1489      DEALLOCATE (DENSO4)
1490      DEALLOCATE (f_r_wet)
1491      DEALLOCATE (decfluxaer)
1492      DEALLOCATE (mdw)
1493      DEALLOCATE (SO2_lifetime)
1494      DEALLOCATE (OCS_lifetime)
1495      DEALLOCATE (alpha_bin)
1496      DEALLOCATE (piz_bin)
1497      DEALLOCATE (cg_bin)
1498      DEALLOCATE (tau_strat_550)
1499      DEALLOCATE (tau_strat_1020)
1500      DEALLOCATE (tausum_strat)
1501      DEALLOCATE (surf_PM25_sulf)
1502      DEALLOCATE (vsed_aer)
1503      DEALLOCATE (budg_3D_ocs_to_so2)
1504      DEALLOCATE (budg_3D_so2_to_h2so4)
1505      DEALLOCATE (budg_3D_backgr_ocs)
1506      DEALLOCATE (budg_3D_backgr_so2)
1507      DEALLOCATE (budg_3D_nucl)
1508      DEALLOCATE (budg_3D_cond_evap)
1509      DEALLOCATE (budg_dep_dry_ocs)
1510      DEALLOCATE (budg_dep_wet_ocs)
1511      DEALLOCATE (budg_dep_dry_so2)
1512      DEALLOCATE (budg_dep_wet_so2)
1513      DEALLOCATE (budg_dep_dry_h2so4)
1514      DEALLOCATE (budg_dep_wet_h2so4)
1515      DEALLOCATE (budg_dep_dry_part)
1516      DEALLOCATE (budg_dep_wet_part)
1517      DEALLOCATE (budg_emi_ocs)
1518      DEALLOCATE (budg_emi_so2)
1519      DEALLOCATE (budg_emi_h2so4)
1520      DEALLOCATE (budg_emi_part)
1521      DEALLOCATE (budg_ocs_to_so2)
1522      DEALLOCATE (budg_so2_to_h2so4)
1523      DEALLOCATE (budg_h2so4_to_part)
1524      DEALLOCATE (budg_sed_part)
1525#endif
1526
1527END SUBROUTINE phys_local_var_end
1528
1529END MODULE phys_local_var_mod
Note: See TracBrowser for help on using the repository browser.