source: LMDZ6/branches/contrails/libf/phylmd/phys_local_var_mod.F90 @ 5579

Last change on this file since 5579 was 5579, checked in by aborella, 3 months ago

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