source: LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90 @ 5649

Last change on this file since 5649 was 5649, checked in by lebasn, 7 weeks ago

StratAer?: bugfix wrong dimension allocate for extinction outputs.

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