Context Navigation

phys_local_var_mod.F90 @ 4799

Last change on this file since 4799 was 4753, checked in by lfalletti, 14 months ago
Add d_q_emiss declaration for Reprobus, which was missing.
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: 54.0 KB

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

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

Download in other formats:

Original Format