Context Navigation

phys_local_var_mod.F90 @ 4601

Last change on this file since 4601 was 4576, checked in by Laurent Fairhead, 17 months ago
Forgot modifications in phylmdiso LF
File size: 62.5 KB

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

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format