Context Navigation

phys_local_var_mod.F90 @ 3807

Last change on this file since 3807 was 2992, checked in by oboucher, 7 years ago
Adding p, z, and t of tropopause in output files. Stratosphere_mask is streamlined a bit. This includes a treatment of missing values.
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 `Author Date Id Revision`
File size: 44.1 KB

Rev	Line
[1279]	1	!
	2	! $Id: phys_local_var_mod.F90 2992 2017-09-21 09:51:25Z fairhead $
	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)
	16	REAL, SAVE, ALLOCATABLE :: u_seri(:,:), v_seri(:,:)
	17	!$OMP THREADPRIVATE(u_seri, v_seri)
[2561]	18	REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:), l_mix(:,:,:)
	19	!$OMP THREADPRIVATE(l_mixmin, l_mix)
[904]	20
[987]	21	REAL, SAVE, ALLOCATABLE :: tr_seri(:,:,:)
	22	!$OMP THREADPRIVATE(tr_seri)
	23	REAL, SAVE, ALLOCATABLE :: d_t_dyn(:,:), d_q_dyn(:,:)
	24	!$OMP THREADPRIVATE(d_t_dyn, d_q_dyn)
[2496]	25	REAL, SAVE, ALLOCATABLE :: d_ql_dyn(:,:), d_qs_dyn(:,:)
	26	!$OMP THREADPRIVATE(d_ql_dyn, d_qs_dyn)
[2499]	27	REAL, SAVE, ALLOCATABLE :: d_q_dyn2d(:), d_ql_dyn2d(:), d_qs_dyn2d(:)
	28	!$OMP THREADPRIVATE(d_q_dyn2d, d_ql_dyn2d, d_qs_dyn2d)
[1054]	29	REAL, SAVE, ALLOCATABLE :: d_u_dyn(:,:), d_v_dyn(:,:)
	30	!$OMP THREADPRIVATE(d_u_dyn, d_v_dyn)
[1742]	31	REAL, SAVE, ALLOCATABLE :: d_tr_dyn(:,:,:)
	32	!$OMP THREADPRIVATE(d_tr_dyn)
[987]	33	REAL, SAVE, ALLOCATABLE :: d_t_con(:,:),d_q_con(:,:)
	34	!$OMP THREADPRIVATE(d_t_con,d_q_con)
	35	REAL, SAVE, ALLOCATABLE :: d_u_con(:,:),d_v_con(:,:)
	36	!$OMP THREADPRIVATE(d_u_con,d_v_con)
	37	REAL, SAVE, ALLOCATABLE :: d_t_wake(:,:),d_q_wake(:,:)
	38	!$OMP THREADPRIVATE( d_t_wake,d_q_wake)
[2086]	39	REAL, SAVE, ALLOCATABLE :: d_t_lsc(:,:),d_q_lsc(:,:),d_ql_lsc(:,:),d_qi_lsc(:,:)
[2092]	40	!$OMP THREADPRIVATE(d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc)
[2194]	41	REAL, SAVE, ALLOCATABLE :: d_t_lwr(:,:),d_t_lw0(:,:),d_t_swr(:,:),d_t_sw0(:,:)
	42	!$OMP THREADPRIVATE(d_t_lwr,d_t_lw0,d_t_swr,d_t_sw0)
[987]	43	REAL, SAVE, ALLOCATABLE :: d_t_ajsb(:,:), d_q_ajsb(:,:)
	44	!$OMP THREADPRIVATE(d_t_ajsb, d_q_ajsb)
	45	REAL, SAVE, ALLOCATABLE :: d_t_ajs(:,:), d_q_ajs(:,:)
	46	!$OMP THREADPRIVATE(d_t_ajs, d_q_ajs)
	47	REAL, SAVE, ALLOCATABLE :: d_u_ajs(:,:), d_v_ajs(:,:)
	48	!$OMP THREADPRIVATE(d_u_ajs, d_v_ajs)
[2159]	49	!nrlmd<
	50	REAL, SAVE, ALLOCATABLE :: d_t_ajs_w(:,:), d_q_ajs_w(:,:)
	51	!$OMP THREADPRIVATE(d_t_ajs_w, d_q_ajs_w)
	52	REAL, SAVE, ALLOCATABLE :: d_t_ajs_x(:,:), d_q_ajs_x(:,:)
	53	!$OMP THREADPRIVATE(d_t_ajs_x, d_q_ajs_x)
	54	!>nrlmd
[2706]	55	REAL, SAVE, ALLOCATABLE :: d_t_eva(:,:),d_q_eva(:,:),d_ql_eva(:,:),d_qi_eva(:,:)
	56	!$OMP THREADPRIVATE(d_t_eva,d_q_eva,d_ql_eva,d_qi_eva)
[1507]	57	REAL, SAVE, ALLOCATABLE :: d_t_lscst(:,:),d_q_lscst(:,:)
	58	!$OMP THREADPRIVATE(d_t_lscst,d_q_lscst)
	59	REAL, SAVE, ALLOCATABLE :: d_t_lscth(:,:),d_q_lscth(:,:)
	60	!$OMP THREADPRIVATE(d_t_lscth,d_q_lscth)
	61	REAL, SAVE, ALLOCATABLE :: plul_th(:),plul_st(:)
	62	!$OMP THREADPRIVATE(plul_th,plul_st)
[909]	63	!tendances dues a oro et lif
[987]	64	REAL, SAVE, ALLOCATABLE :: d_t_oli(:,:)
	65	!$OMP THREADPRIVATE(d_t_oli)
	66	REAL, SAVE, ALLOCATABLE :: d_u_oli(:,:), d_v_oli(:,:)
	67	!$OMP THREADPRIVATE(d_u_oli, d_v_oli)
[1761]	68	REAL, SAVE, ALLOCATABLE :: d_t_vdf(:,:), d_q_vdf(:,:), d_t_diss(:,:)
	69	!$OMP THREADPRIVATE( d_t_vdf, d_q_vdf,d_t_diss)
[987]	70	REAL, SAVE, ALLOCATABLE :: d_u_vdf(:,:), d_v_vdf(:,:)
	71	!$OMP THREADPRIVATE(d_u_vdf, d_v_vdf)
[2159]	72	!nrlmd+jyg<
	73	REAL, SAVE, ALLOCATABLE :: d_t_vdf_w(:,:), d_q_vdf_w(:,:)
	74	!$OMP THREADPRIVATE( d_t_vdf_w, d_q_vdf_w)
	75	REAL, SAVE, ALLOCATABLE :: d_t_vdf_x(:,:), d_q_vdf_x(:,:)
	76	!$OMP THREADPRIVATE( d_t_vdf_x, d_q_vdf_x)
	77	!>nrlmd+jyg
[987]	78	REAL, SAVE, ALLOCATABLE :: d_t_oro(:,:)
	79	!$OMP THREADPRIVATE(d_t_oro)
	80	REAL, SAVE, ALLOCATABLE :: d_u_oro(:,:), d_v_oro(:,:)
	81	!$OMP THREADPRIVATE(d_u_oro, d_v_oro)
[2897]	82	REAL, SAVE, ALLOCATABLE :: d_t_oro_gw(:,:)
[2898]	83	!$OMP THREADPRIVATE(d_t_oro_gw)
[2897]	84	REAL, SAVE, ALLOCATABLE :: d_u_oro_gw(:,:), d_v_oro_gw(:,:)
	85	!$OMP THREADPRIVATE(d_u_oro_gw, d_v_oro_gw)
[987]	86	REAL, SAVE, ALLOCATABLE :: d_t_lif(:,:)
	87	!$OMP THREADPRIVATE(d_t_lif)
	88	REAL, SAVE, ALLOCATABLE :: d_u_lif(:,:), d_v_lif(:,:)
	89	!$OMP THREADPRIVATE(d_u_lif, d_v_lif)
[1001]	90	! Tendances Ondes de G non oro (runs strato).
[2333]	91	REAL, SAVE, ALLOCATABLE :: du_gwd_hines(:,:)
	92	!$OMP THREADPRIVATE(du_gwd_hines)
	93	REAL, SAVE, ALLOCATABLE :: dv_gwd_hines(:,:)
	94	!$OMP THREADPRIVATE(dv_gwd_hines)
	95	REAL, SAVE, ALLOCATABLE :: dv_gwd_rando(:,:)
	96	!$OMP THREADPRIVATE(dv_gwd_rando)
	97	REAL, SAVE, ALLOCATABLE :: dv_gwd_front(:,:)
	98	!$OMP THREADPRIVATE(dv_gwd_front)
	99	REAL, SAVE, ALLOCATABLE :: east_gwstress(:,:)
	100	!$OMP THREADPRIVATE(east_gwstress)
	101	REAL, SAVE, ALLOCATABLE :: west_gwstress(:,:)
	102	!$OMP THREADPRIVATE(west_gwstress)
[1001]	103	REAL, SAVE, ALLOCATABLE :: d_t_hin(:,:)
	104	!$OMP THREADPRIVATE(d_t_hin)
[2136]	105	! tendance due a l'oxydation du methane
	106	REAL, SAVE, ALLOCATABLE :: d_q_ch4(:,:)
	107	!$OMP THREADPRIVATE(d_q_ch4)
[1001]	108
[909]	109	! tendance du a la conersion Ec -> E thermique
[987]	110	REAL, SAVE, ALLOCATABLE :: d_t_ec(:,:)
	111	!$OMP THREADPRIVATE(d_t_ec)
	112	REAL, SAVE, ALLOCATABLE :: d_ts(:,:), d_tr(:,:,:)
	113	!$OMP THREADPRIVATE(d_ts, d_tr)
[1279]	114
	115	! diagnostique pour le rayonnement
	116	REAL, SAVE, ALLOCATABLE :: topswad_aero(:), solswad_aero(:) ! diag
	117	!$OMP THREADPRIVATE(topswad_aero,solswad_aero)
	118	REAL, SAVE, ALLOCATABLE :: topswai_aero(:), solswai_aero(:) ! diag
	119	!$OMP THREADPRIVATE(topswai_aero,solswai_aero)
	120	REAL, SAVE, ALLOCATABLE :: topswad0_aero(:), solswad0_aero(:) ! diag
	121	!$OMP THREADPRIVATE(topswad0_aero,solswad0_aero)
	122	REAL, SAVE, ALLOCATABLE :: topsw_aero(:,:), solsw_aero(:,:) ! diag
	123	!$OMP THREADPRIVATE(topsw_aero,solsw_aero)
	124	REAL, SAVE, ALLOCATABLE :: topsw0_aero(:,:), solsw0_aero(:,:) ! diag
	125	!$OMP THREADPRIVATE(topsw0_aero,solsw0_aero)
	126	REAL, SAVE, ALLOCATABLE :: topswcf_aero(:,:), solswcf_aero(:,:) ! diag
	127	!$OMP THREADPRIVATE(topswcf_aero,solswcf_aero)
[2146]	128	! LW radiation diagnostics CK
	129	REAL, SAVE, ALLOCATABLE :: toplwad_aero(:), sollwad_aero(:) ! diag
	130	!$OMP THREADPRIVATE(toplwad_aero,sollwad_aero)
	131	REAL, SAVE, ALLOCATABLE :: toplwai_aero(:), sollwai_aero(:) ! diag
	132	!$OMP THREADPRIVATE(toplwai_aero,sollwai_aero)
	133	REAL, SAVE, ALLOCATABLE :: toplwad0_aero(:), sollwad0_aero(:) ! diag
	134	!$OMP THREADPRIVATE(toplwad0_aero,sollwad0_aero)
[1989]	135	! Special RRTM
	136	REAL, SAVE, ALLOCATABLE :: ZLWFT0_i(:,:), ZSWFT0_i(:,:) ! diag
	137	!$OMP THREADPRIVATE(ZLWFT0_i,ZSWFT0_i)
	138	REAL, SAVE, ALLOCATABLE :: ZFLDN0(:,:), ZFLUP0(:,:) ! diag
	139	!$OMP THREADPRIVATE(ZFLDN0,ZFLUP0)
	140	REAL, SAVE, ALLOCATABLE :: ZFSDN0(:,:), ZFSUP0(:,:) ! diag
	141	!$OMP THREADPRIVATE(ZFSDN0,ZFSUP0)
	142	!
[2536]	143	REAL, SAVE, ALLOCATABLE :: stratomask(:,:)
	144	!$OMP THREADPRIVATE(stratomask)
[1279]	145	REAL, SAVE, ALLOCATABLE :: tausum_aero(:,:,:)
	146	!$OMP THREADPRIVATE(tausum_aero)
[2854]	147	REAL, SAVE, ALLOCATABLE :: drytausum_aero(:,:)
	148	!$OMP THREADPRIVATE(drytausum_aero)
[1279]	149	REAL, SAVE, ALLOCATABLE :: tau3d_aero(:,:,:,:)
	150	!$OMP THREADPRIVATE(tau3d_aero)
[1337]	151	REAL, SAVE, ALLOCATABLE :: scdnc(:,:)
	152	!$OMP THREADPRIVATE(scdnc)
	153	REAL, SAVE, ALLOCATABLE :: cldncl(:)
	154	!$OMP THREADPRIVATE(cldncl)
	155	REAL, SAVE, ALLOCATABLE :: reffclwtop(:)
	156	!$OMP THREADPRIVATE(reffclwtop)
	157	REAL, SAVE, ALLOCATABLE :: lcc(:)
	158	!$OMP THREADPRIVATE(lcc)
	159	REAL, SAVE, ALLOCATABLE :: reffclws(:,:)
	160	!$OMP THREADPRIVATE(reffclws)
	161	REAL, SAVE, ALLOCATABLE :: reffclwc(:,:)
	162	!$OMP THREADPRIVATE(reffclwc)
	163	REAL, SAVE, ALLOCATABLE :: cldnvi(:)
	164	!$OMP THREADPRIVATE(cldnvi)
	165	REAL, SAVE, ALLOCATABLE :: lcc3d(:,:)
	166	!$OMP THREADPRIVATE(lcc3d)
	167	REAL, SAVE, ALLOCATABLE :: lcc3dcon(:,:)
	168	!$OMP THREADPRIVATE(lcc3dcon)
	169	REAL, SAVE, ALLOCATABLE :: lcc3dstra(:,:)
	170	!$OMP THREADPRIVATE(lcc3dstra)
[2806]	171	REAL, SAVE, ALLOCATABLE :: od443aer(:)
	172	!$OMP THREADPRIVATE(od443aer)
[1337]	173	REAL, SAVE, ALLOCATABLE :: od550aer(:)
	174	!$OMP THREADPRIVATE(od550aer)
[2854]	175	REAL, SAVE, ALLOCATABLE :: dryod550aer(:)
	176	!$OMP THREADPRIVATE(dryod550aer)
[2842]	177	REAL, SAVE, ALLOCATABLE :: abs550aer(:)
	178	!$OMP THREADPRIVATE(abs550aer)
[1337]	179	REAL, SAVE, ALLOCATABLE :: od865aer(:)
	180	!$OMP THREADPRIVATE(od865aer)
	181	REAL, SAVE, ALLOCATABLE :: ec550aer(:,:)
	182	!$OMP THREADPRIVATE(ec550aer)
	183	REAL, SAVE, ALLOCATABLE :: od550lt1aer(:)
	184	!$OMP THREADPRIVATE(od550lt1aer)
	185	REAL, SAVE, ALLOCATABLE :: sconcso4(:)
	186	!$OMP THREADPRIVATE(sconcso4)
[2146]	187	REAL, SAVE, ALLOCATABLE :: sconcno3(:)
	188	!$OMP THREADPRIVATE(sconcno3)
[1337]	189	REAL, SAVE, ALLOCATABLE :: sconcoa(:)
	190	!$OMP THREADPRIVATE(sconcoa)
	191	REAL, SAVE, ALLOCATABLE :: sconcbc(:)
	192	!$OMP THREADPRIVATE(sconcbc)
	193	REAL, SAVE, ALLOCATABLE :: sconcss(:)
	194	!$OMP THREADPRIVATE(sconcss)
	195	REAL, SAVE, ALLOCATABLE :: sconcdust(:)
	196	!$OMP THREADPRIVATE(sconcdust)
	197	REAL, SAVE, ALLOCATABLE :: concso4(:,:)
	198	!$OMP THREADPRIVATE(concso4)
[2146]	199	REAL, SAVE, ALLOCATABLE :: concno3(:,:)
	200	!$OMP THREADPRIVATE(concno3)
[1337]	201	REAL, SAVE, ALLOCATABLE :: concoa(:,:)
	202	!$OMP THREADPRIVATE(concoa)
	203	REAL, SAVE, ALLOCATABLE :: concbc(:,:)
	204	!$OMP THREADPRIVATE(concbc)
	205	REAL, SAVE, ALLOCATABLE :: concss(:,:)
	206	!$OMP THREADPRIVATE(concss)
	207	REAL, SAVE, ALLOCATABLE :: concdust(:,:)
	208	!$OMP THREADPRIVATE(concdust)
	209	REAL, SAVE, ALLOCATABLE :: loadso4(:)
	210	!$OMP THREADPRIVATE(loadso4)
	211	REAL, SAVE, ALLOCATABLE :: loadoa(:)
	212	!$OMP THREADPRIVATE(loadoa)
	213	REAL, SAVE, ALLOCATABLE :: loadbc(:)
	214	!$OMP THREADPRIVATE(loadbc)
	215	REAL, SAVE, ALLOCATABLE :: loadss(:)
	216	!$OMP THREADPRIVATE(loadss)
	217	REAL, SAVE, ALLOCATABLE :: loaddust(:)
	218	!$OMP THREADPRIVATE(loaddust)
[2823]	219	REAL, SAVE, ALLOCATABLE :: loadno3(:)
	220	!$OMP THREADPRIVATE(loadno3)
[2953]	221	REAL, SAVE, ALLOCATABLE :: load_tmp1(:), load_tmp2(:)
	222	!$OMP THREADPRIVATE(load_tmp1, load_tmp2)
	223	REAL, SAVE, ALLOCATABLE :: load_tmp3(:), load_tmp4(:)
	224	!$OMP THREADPRIVATE(load_tmp3, load_tmp4)
	225	REAL, SAVE, ALLOCATABLE :: load_tmp5(:), load_tmp6(:)
	226	!$OMP THREADPRIVATE(load_tmp5, load_tmp6)
	227	REAL, SAVE, ALLOCATABLE :: load_tmp7(:), load_tmp8(:)
	228	!$OMP THREADPRIVATE(load_tmp7, load_tmp8)
	229	REAL, SAVE, ALLOCATABLE :: load_tmp9(:), load_tmp10(:)
	230	!$OMP THREADPRIVATE(load_tmp9, load_tmp10)
[1279]	231
[1539]	232	!IM ajout variables CFMIP2/CMIP5
	233	REAL,ALLOCATABLE,SAVE :: topswad_aerop(:), solswad_aerop(:)
	234	!$OMP THREADPRIVATE(topswad_aerop, solswad_aerop)
	235	REAL,ALLOCATABLE,SAVE :: topswai_aerop(:), solswai_aerop(:)
	236	!$OMP THREADPRIVATE(topswai_aerop, solswai_aerop)
	237	REAL,ALLOCATABLE,SAVE :: topswad0_aerop(:), solswad0_aerop(:)
	238	!$OMP THREADPRIVATE(topswad0_aerop, solswad0_aerop)
	239	REAL,ALLOCATABLE,SAVE :: topsw_aerop(:,:), topsw0_aerop(:,:)
	240	!$OMP THREADPRIVATE(topsw_aerop, topsw0_aerop)
	241	REAL,ALLOCATABLE,SAVE :: solsw_aerop(:,:), solsw0_aerop(:,:)
	242	!$OMP THREADPRIVATE(solsw_aerop, solsw0_aerop)
	243	REAL,ALLOCATABLE,SAVE :: topswcf_aerop(:,:), solswcf_aerop(:,:)
	244	!$OMP THREADPRIVATE(topswcf_aerop, solswcf_aerop)
	245
[2146]	246	! additional LW variables CK
	247	REAL,ALLOCATABLE,SAVE :: toplwad_aerop(:), sollwad_aerop(:)
	248	!$OMP THREADPRIVATE(toplwad_aerop, sollwad_aerop)
	249	REAL,ALLOCATABLE,SAVE :: toplwai_aerop(:), sollwai_aerop(:)
	250	!$OMP THREADPRIVATE(toplwai_aerop, sollwai_aerop)
	251	REAL,ALLOCATABLE,SAVE :: toplwad0_aerop(:), sollwad0_aerop(:)
	252	!$OMP THREADPRIVATE(toplwad0_aerop, sollwad0_aerop)
[1797]	253
[2897]	254	!Ajout de celles n??cessaires au phys_output_write_mod
[2385]	255	REAL, SAVE, ALLOCATABLE :: tal1(:), pal1(:), pab1(:), pab2(:)
	256	!$OMP THREADPRIVATE(tal1, pal1, pab1, pab2)
	257	REAL, SAVE, ALLOCATABLE :: ptstar(:), pt0(:), slp(:)
	258	!$OMP THREADPRIVATE(ptstar, pt0, slp)
[1797]	259	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: sens, flwp, fiwp
	260	!$OMP THREADPRIVATE(sens, flwp, fiwp)
[2635]	261	!!
	262	!! Wake variables
	263	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: ale_wake, alp_wake
	264	!$OMP THREADPRIVATE(ale_wake, alp_wake)
[2849]	265	!!jyg! REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wake_h,wake_k
	266	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wake_h
	267	INTEGER,ALLOCATABLE,SAVE,DIMENSION(:) :: wake_k
[2732]	268	!$OMP THREADPRIVATE(wake_h,wake_k)
[2635]	269	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wake_omg
	270	!$OMP THREADPRIVATE(wake_omg)
	271	REAL, SAVE, ALLOCATABLE,DIMENSION(:,:) :: d_deltat_wk, d_deltaq_wk
	272	!$OMP THREADPRIVATE(d_deltat_wk, d_deltaq_wk)
	273	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: d_s_wk, d_dens_wk
	274	!$OMP THREADPRIVATE(d_s_wk, d_dens_wk)
	275	REAL, SAVE, ALLOCATABLE,DIMENSION(:,:) :: d_deltat_wk_gw, d_deltaq_wk_gw
	276	!$OMP THREADPRIVATE(d_deltat_wk_gw, d_deltaq_wk_gw)
	277	REAL, SAVE, ALLOCATABLE,DIMENSION(:,:) :: d_deltat_vdf, d_deltaq_vdf
	278	!$OMP THREADPRIVATE(d_deltat_vdf, d_deltaq_vdf)
	279	!!! REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: d_s_vdf, d_dens_vdf
	280	!!!OMP THREADPRIVATE(d_s_vdf, d_dens_vdf)
	281	REAL, SAVE, ALLOCATABLE,DIMENSION(:,:) :: d_deltat_the, d_deltaq_the
	282	!$OMP THREADPRIVATE(d_deltat_the, d_deltaq_the)
	283	!!! REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: d_s_the, d_dens_the
	284	!!!OMP THREADPRIVATE(d_s_the, d_dens_the)
	285	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: d_deltat_ajs_cv, d_deltaq_ajs_cv
	286	!$OMP THREADPRIVATE(d_deltat_ajs_cv, d_deltaq_ajs_cv)
	287	!! End of Wake variables
	288	!!
	289	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: bils
	290	!$OMP THREADPRIVATE(bils)
[1797]	291	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cdragm, cdragh
	292	!$OMP THREADPRIVATE(cdragm, cdragh)
[2243]	293	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cldh, cldl, cldm, cldq, cldt, qsat2m
	294	!$OMP THREADPRIVATE(cldh, cldl, cldm, cldq, cldt, qsat2m )
[1924]	295	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cldhjn, cldljn, cldmjn,cldtjn
	296	!$OMP THREADPRIVATE(cldhjn, cldljn, cldmjn, cldtjn)
	297	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: JrNt
	298	!$OMP THREADPRIVATE(JrNt)
[2496]	299	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: dthmin, evap, fder, plcl, plfc, prw, prlw, prsw
	300	!$OMP THREADPRIVATE(dthmin, evap, fder, plcl, plfc, prw, prlw, prsw)
[1797]	301	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zustar, zu10m, zv10m, rh2m
	302	!$OMP THREADPRIVATE(zustar, zu10m, zv10m, rh2m)
	303	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: s_lcl, s_pblh, s_pblt, s_therm
	304	!$OMP THREADPRIVATE(s_lcl, s_pblh, s_pblt, s_therm)
[2159]	305	!
	306	!nrlmd+jyg<
	307	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: s_pblh_x, s_pblh_w
	308	!$OMP THREADPRIVATE(s_pblh_x, s_pblh_w)
	309	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: s_lcl_x, s_lcl_w
	310	!$OMP THREADPRIVATE(s_lcl_x, s_lcl_w)
	311	!>nrlmd+jyg
	312	!
[1797]	313	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: slab_wfbils
	314	!$OMP THREADPRIVATE(slab_wfbils)
	315	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: tpot, tpote, ue, uq, ve, vq, zxffonte
	316	!$OMP THREADPRIVATE(tpot, tpote, ue, uq, ve, vq, zxffonte)
	317	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfqcalving
	318	!$OMP THREADPRIVATE(zxfqcalving)
[2243]	319	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfluxlat, zxtsol, snow_lsc, zxfqfonte
	320	!$OMP THREADPRIVATE(zxfluxlat, zxtsol, snow_lsc, zxfqfonte)
[2518]	321	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxrunofflic
	322	!$OMP THREADPRIVATE(zxrunofflic)
[2516]	323	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxqsurf, rain_lsc, rain_num
	324	!$OMP THREADPRIVATE(zxqsurf, rain_lsc, rain_num)
[2159]	325	!
	326	!jyg+nrlmd<
	327	!!!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
	328	! c
	329	! Declarations liees a la couche limite differentiee w-x c
	330	! c
	331	!!!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
	332	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: sens_x, sens_w
	333	!$OMP THREADPRIVATE(sens_x, sens_w)
	334	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfluxlat_x, zxfluxlat_w
	335	!$OMP THREADPRIVATE(zxfluxlat_x, zxfluxlat_w)
[2181]	336	!jyg<
[2897]	337	!!! Entr\E9es suppl\E9mentaires couche-limite
[2181]	338	!! REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: t_x, t_w
	339	!!!$OMP THREADPRIVATE(t_x, t_w)
	340	!! REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: q_x, q_w
	341	!!!$OMP THREADPRIVATE(q_x, q_w)
	342	!>jyg
	343	!!! Sorties ferret
[2159]	344	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dtvdf_x, dtvdf_w
	345	!$OMP THREADPRIVATE(dtvdf_x, dtvdf_w)
	346	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dqvdf_x, dqvdf_w
	347	!$OMP THREADPRIVATE(dqvdf_x, dqvdf_w)
[2897]	348	! Variables suppl\E9mentaires dans physiq.F relative au splitting de la surface
[2159]	349	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: pbl_tke_input
	350	!$OMP THREADPRIVATE(pbl_tke_input)
	351	! Entree supplementaire Thermiques :
[2607]	352	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: t_therm, q_therm, u_therm, v_therm
	353	!$OMP THREADPRIVATE(t_therm, q_therm, u_therm, v_therm)
[2159]	354	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cdragh_x, cdragh_w
	355	!$OMP THREADPRIVATE(cdragh_x, cdragh_w)
	356	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cdragm_x, cdragm_w
	357	!$OMP THREADPRIVATE(cdragm_x, cdragm_w)
	358	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: kh, kh_x, kh_w
	359	!$OMP THREADPRIVATE(kh, kh_x, kh_w)
	360	!!!
	361	!!!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
[2877]	362	LOGICAL, SAVE, ALLOCATABLE :: ptconv(:,:)
	363	!$OMP THREADPRIVATE(ptconv)
[2159]	364	!>jyg+nrlmd
	365	!
[2824]	366	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wbeff, convoccur, zmax_th, zq2m, zt2m
	367	!$OMP THREADPRIVATE(wbeff, convoccur, zmax_th, zq2m, zt2m)
[2103]	368	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zt2m_min_mon, zt2m_max_mon
	369	!$OMP THREADPRIVATE(zt2m_min_mon, zt2m_max_mon)
	370	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: t2m_min_mon, t2m_max_mon
	371	!$OMP THREADPRIVATE(t2m_min_mon, t2m_max_mon)
[1797]	372	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: weak_inversion
	373	!$OMP THREADPRIVATE(weak_inversion)
	374	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: ale_bl_stat,alp_bl_conv,alp_bl_det
	375	!$OMP THREADPRIVATE(ale_bl_stat,alp_bl_conv,alp_bl_det)
	376	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_fluct_m,alp_bl_fluct_tke
	377	!$OMP THREADPRIVATE(alp_bl_fluct_m,alp_bl_fluct_tke)
	378	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_stat, n2, s2
	379	!$OMP THREADPRIVATE(alp_bl_stat, n2, s2)
	380	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: proba_notrig, random_notrig
	381	!$OMP THREADPRIVATE(proba_notrig, random_notrig)
	382	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: fsolsw, wfbils, wfbilo
	383	!$OMP THREADPRIVATE(fsolsw, wfbils, wfbilo)
[2670]	384	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wfevap, wfrain, wfsnow
	385	!$OMP THREADPRIVATE(wfevap,wfrain,wfsnow)
[2243]	386	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: t2m, fluxlat, fsollw,evap_pot
	387	!$OMP THREADPRIVATE(t2m, fluxlat, fsollw,evap_pot)
[1797]	388	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dnwd, dnwd0, upwd, omega
	389	!$OMP THREADPRIVATE(dnwd, dnwd0, upwd, omega)
[2481]	390	REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: epmax_diag ! epmax_cape
	391	!$OMP THREADPRIVATE(epmax_diag)
	392	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ep ! epmax_cape
	393	!$OMP THREADPRIVATE(ep)
[1797]	394	! REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: coefh, coefm, lambda_th
	395	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: lambda_th
	396	!$OMP THREADPRIVATE(lambda_th)
	397	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: cldemi, cldfra, cldtau, fiwc, fl, re, flwc
	398	!$OMP THREADPRIVATE(cldemi, cldfra, cldtau, fiwc, fl, re, flwc)
	399	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq, ref_ice, theta, zphi
	400	!$OMP THREADPRIVATE(ref_liq, ref_ice, theta, zphi)
[1989]	401	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq_pi, ref_ice_pi
	402	!$OMP THREADPRIVATE(ref_liq_pi, ref_ice_pi)
[2635]	403	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zx_rh
	404	!$OMP THREADPRIVATE(zx_rh)
[1797]	405	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: pmflxr, pmflxs, prfl, psfl, fraca
	406	!$OMP THREADPRIVATE(pmflxr, pmflxs, prfl, psfl, fraca)
	407	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: Vprecip, zw2
	408	!$OMP THREADPRIVATE(Vprecip, zw2)
	409	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: pmfd, pmfu
	410	!$OMP THREADPRIVATE(pmfd, pmfu)
	411	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: fluxt, fluxu, fluxv
	412	!$OMP THREADPRIVATE(fluxt, fluxu, fluxv)
	413	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: uwriteSTD, vwriteSTD, wwriteSTD
	414	!$OMP THREADPRIVATE(uwriteSTD, vwriteSTD, wwriteSTD)
[1828]	415	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: phiwriteSTD, qwriteSTD, twriteSTD, rhwriteSTD
	416	!$OMP THREADPRIVATE(phiwriteSTD, qwriteSTD, twriteSTD, rhwriteSTD)
[1797]	417
	418	! ug et d'autres encore:
	419	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wdtrainA, wdtrainM
	420	!$OMP THREADPRIVATE(wdtrainA, wdtrainM)
	421	REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: beta_prec
	422	!$OMP THREADPRIVATE(beta_prec)
[2945]	423	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rneb,rnebjn,rneblsvol
	424	!$OMP THREADPRIVATE(rneb,rnebjn,rneblsvol)
[1797]	425
[1865]	426	! variables de sorties MM
[1872]	427	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: zxsnow,snowhgt,qsnow,to_ice
	428	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: sissnow,runoff,albsol3_lic
	429	!$OMP THREADPRIVATE(zxsnow,snowhgt,qsnow,to_ice)
	430	!$OMP THREADPRIVATE(sissnow,runoff,albsol3_lic)
[1865]	431
[2992]	432	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: p_tropopause, z_tropopause, t_tropopause
	433	!$OMP THREADPRIVATE(p_tropopause, z_tropopause, t_tropopause)
	434
[2690]	435	#ifdef CPP_StratAer
[2752]	436	!
	437	! variables for stratospheric aerosol
[2690]	438	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: R2SO4
	439	!$OMP THREADPRIVATE(R2SO4)
	440	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: DENSO4
	441	!$OMP THREADPRIVATE(DENSO4)
	442	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: f_r_wet
	443	!$OMP THREADPRIVATE(f_r_wet)
	444	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: decfluxaer
	445	!$OMP THREADPRIVATE(decfluxaer)
	446	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: mdw
	447	!$OMP THREADPRIVATE(mdw)
	448	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: OCS_lifetime
	449	!$OMP THREADPRIVATE(OCS_lifetime)
	450	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: SO2_lifetime
	451	!$OMP THREADPRIVATE(SO2_lifetime)
	452	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: alpha_bin
	453	!$OMP THREADPRIVATE(alpha_bin)
	454	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: piz_bin
	455	!$OMP THREADPRIVATE(piz_bin)
	456	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: cg_bin
	457	!$OMP THREADPRIVATE(cg_bin)
	458	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_550
	459	!$OMP THREADPRIVATE(tau_strat_550)
	460	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_1020
	461	!$OMP THREADPRIVATE(tau_strat_1020)
	462	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tausum_strat
	463	!$OMP THREADPRIVATE(tausum_strat)
	464	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: surf_PM25_sulf
	465	!$OMP THREADPRIVATE(surf_PM25_sulf)
	466	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: vsed_aer
	467	!$OMP THREADPRIVATE(vsed_aer)
[2752]	468	!
	469	!---3D budget variables
	470	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_nucl
	471	!$OMP THREADPRIVATE(budg_3D_nucl)
	472	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_cond_evap
	473	!$OMP THREADPRIVATE(budg_3D_cond_evap)
	474	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_ocs_to_so2
	475	!$OMP THREADPRIVATE(budg_3D_ocs_to_so2)
	476	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_so2_to_h2so4
	477	!$OMP THREADPRIVATE(budg_3D_so2_to_h2so4)
	478	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_backgr_ocs
	479	!$OMP THREADPRIVATE(budg_3D_backgr_ocs)
	480	REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_backgr_so2
	481	!$OMP THREADPRIVATE(budg_3D_backgr_so2)
	482	!
	483	!---2D budget variables
	484	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_h2so4
	485	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_part
	486	!$OMP THREADPRIVATE(budg_dep_dry_h2so4,budg_dep_dry_part)
	487	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_h2so4
	488	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_part
	489	!$OMP THREADPRIVATE(budg_dep_wet_h2so4,budg_dep_wet_part)
	490	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_ocs
	491	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_so2
	492	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_h2so4
	493	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_part
	494	!$OMP THREADPRIVATE(budg_emi_ocs,budg_emi_so2,budg_emi_h2so4,budg_emi_part)
	495	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_ocs_to_so2
	496	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_so2_to_h2so4
	497	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_h2so4_to_part
	498	!$OMP THREADPRIVATE(budg_ocs_to_so2,budg_so2_to_h2so4,budg_h2so4_to_part)
	499	REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_sed_part
	500	!$OMP THREADPRIVATE(budg_sed_part)
[2690]	501	#endif
	502
[904]	503	CONTAINS
	504
	505	!======================================================================
	506	SUBROUTINE phys_local_var_init
[1785]	507	USE dimphy
[2320]	508	USE infotrac_phy, ONLY : nbtr
[1279]	509	USE aero_mod
[1785]	510	USE indice_sol_mod
[1797]	511	USE phys_output_var_mod
	512	USE phys_state_var_mod
[1279]	513
[904]	514	IMPLICIT NONE
[2854]	515	ALLOCATE(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev))
	516	ALLOCATE(u_seri(klon,klev),v_seri(klon,klev))
	517	ALLOCATE(l_mixmin(klon,klev,nbsrf), l_mix(klon,klev,nbsrf))
[2574]	518	l_mix(:,:,:)=0. ; l_mixmin(:,:,:)=0. ! doit etre initialse car pas toujours remplis
[909]	519
[2854]	520	ALLOCATE(tr_seri(klon,klev,nbtr))
	521	ALLOCATE(d_t_dyn(klon,klev),d_q_dyn(klon,klev))
	522	ALLOCATE(d_ql_dyn(klon,klev),d_qs_dyn(klon,klev))
	523	ALLOCATE(d_q_dyn2d(klon),d_ql_dyn2d(klon),d_qs_dyn2d(klon))
	524	ALLOCATE(d_u_dyn(klon,klev),d_v_dyn(klon,klev))
	525	ALLOCATE(d_tr_dyn(klon,klev,nbtr)) !RomP
	526	ALLOCATE(d_t_con(klon,klev),d_q_con(klon,klev))
	527	ALLOCATE(d_u_con(klon,klev),d_v_con(klon,klev))
	528	ALLOCATE(d_t_wake(klon,klev),d_q_wake(klon,klev))
	529	ALLOCATE(d_t_lsc(klon,klev),d_q_lsc(klon,klev))
	530	ALLOCATE(d_t_lwr(klon,klev),d_t_lw0(klon,klev))
	531	ALLOCATE(d_t_swr(klon,klev),d_t_sw0(klon,klev))
	532	ALLOCATE(d_ql_lsc(klon,klev),d_qi_lsc(klon,klev))
	533	ALLOCATE(d_t_ajsb(klon,klev),d_q_ajsb(klon,klev))
	534	ALLOCATE(d_t_ajs(klon,klev),d_q_ajs(klon,klev))
[2159]	535	!nrlmd<
[2854]	536	ALLOCATE(d_t_ajs_w(klon,klev),d_q_ajs_w(klon,klev))
	537	ALLOCATE(d_t_ajs_x(klon,klev),d_q_ajs_x(klon,klev))
[2159]	538	!>nrlmd
[2854]	539	ALLOCATE(d_u_ajs(klon,klev),d_v_ajs(klon,klev))
	540	ALLOCATE(d_t_eva(klon,klev),d_q_eva(klon,klev))
	541	ALLOCATE(d_ql_eva(klon,klev),d_qi_eva(klon,klev))
	542	ALLOCATE(d_t_lscst(klon,klev),d_q_lscst(klon,klev))
	543	ALLOCATE(d_t_lscth(klon,klev),d_q_lscth(klon,klev))
	544	ALLOCATE(plul_st(klon),plul_th(klon))
	545	ALLOCATE(d_t_vdf(klon,klev),d_q_vdf(klon,klev),d_t_diss(klon,klev))
[2159]	546	!nrlmd+jyg<
[2854]	547	ALLOCATE(d_t_vdf_w(klon,klev),d_q_vdf_w(klon,klev))
	548	ALLOCATE(d_t_vdf_x(klon,klev),d_q_vdf_x(klon,klev))
[2159]	549	!>nrlmd+jyg
[2854]	550	ALLOCATE(d_u_vdf(klon,klev),d_v_vdf(klon,klev))
	551	ALLOCATE(d_t_oli(klon,klev),d_t_oro(klon,klev))
	552	ALLOCATE(d_u_oli(klon,klev),d_v_oli(klon,klev))
	553	ALLOCATE(d_u_oro(klon,klev),d_v_oro(klon,klev))
[2897]	554	ALLOCATE(d_u_oro_gw(klon,klev),d_v_oro_gw(klon,klev))
	555	ALLOCATE(d_t_oro_gw(klon,klev))
[2854]	556	ALLOCATE(d_t_lif(klon,klev),d_t_ec(klon,klev))
	557	ALLOCATE(d_u_lif(klon,klev),d_v_lif(klon,klev))
	558	ALLOCATE(d_ts(klon,nbsrf), d_tr(klon,klev,nbtr))
[1989]	559	! Special RRTM
[2854]	560	ALLOCATE(ZLWFT0_i(klon,klev+1),ZSWFT0_i(klon,klev+1),ZFLDN0(klon,klev+1))
	561	ALLOCATE(ZFLUP0(klon,klev+1),ZFSDN0(klon,klev+1),ZFSUP0(klon,klev+1))
[1989]	562	!
[2854]	563	ALLOCATE(topswad_aero(klon), solswad_aero(klon))
	564	ALLOCATE(topswai_aero(klon), solswai_aero(klon))
	565	ALLOCATE(topswad0_aero(klon), solswad0_aero(klon))
	566	ALLOCATE(toplwad_aero(klon), sollwad_aero(klon))
	567	ALLOCATE(toplwai_aero(klon), sollwai_aero(klon))
	568	ALLOCATE(toplwad0_aero(klon), sollwad0_aero(klon))
	569	ALLOCATE(topsw_aero(klon,naero_grp), solsw_aero(klon,naero_grp))
	570	ALLOCATE(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp))
	571	ALLOCATE(topswcf_aero(klon,3), solswcf_aero(klon,3))
	572	ALLOCATE(du_gwd_hines(klon,klev),dv_gwd_hines(klon,klev))
	573	ALLOCATE(dv_gwd_rando(klon,klev),dv_gwd_front(klon,klev))
	574	ALLOCATE(east_gwstress(klon,klev),west_gwstress(klon,klev))
	575	ALLOCATE(d_t_hin(klon,klev))
	576	ALLOCATE(d_q_ch4(klon,klev))
	577	ALLOCATE(stratomask(klon,klev))
	578	ALLOCATE(tausum_aero(klon,nwave,naero_tot))
	579	ALLOCATE(drytausum_aero(klon,naero_tot))
	580	ALLOCATE(tau3d_aero(klon,klev,nwave,naero_tot))
	581	ALLOCATE(scdnc(klon, klev))
	582	ALLOCATE(cldncl(klon))
	583	ALLOCATE(reffclwtop(klon))
	584	ALLOCATE(lcc(klon))
	585	ALLOCATE(reffclws(klon, klev))
	586	ALLOCATE(reffclwc(klon, klev))
	587	ALLOCATE(cldnvi(klon))
	588	ALLOCATE(lcc3d(klon, klev))
	589	ALLOCATE(lcc3dcon(klon, klev))
	590	ALLOCATE(lcc3dstra(klon, klev))
	591	ALLOCATE(od443aer(klon))
	592	ALLOCATE(od550aer(klon))
	593	ALLOCATE(od865aer(klon))
	594	ALLOCATE(dryod550aer(klon))
	595	ALLOCATE(abs550aer(klon))
	596	ALLOCATE(ec550aer(klon,klev))
	597	ALLOCATE(od550lt1aer(klon))
	598	ALLOCATE(sconcso4(klon))
	599	ALLOCATE(sconcno3(klon))
	600	ALLOCATE(sconcoa(klon))
	601	ALLOCATE(sconcbc(klon))
	602	ALLOCATE(sconcss(klon))
	603	ALLOCATE(sconcdust(klon))
	604	ALLOCATE(concso4(klon,klev))
	605	ALLOCATE(concno3(klon,klev))
	606	ALLOCATE(concoa(klon,klev))
	607	ALLOCATE(concbc(klon,klev))
	608	ALLOCATE(concss(klon,klev))
	609	ALLOCATE(concdust(klon,klev))
	610	ALLOCATE(loadso4(klon))
	611	ALLOCATE(loadoa(klon))
	612	ALLOCATE(loadbc(klon))
	613	ALLOCATE(loadss(klon))
	614	ALLOCATE(loaddust(klon))
	615	ALLOCATE(loadno3(klon))
	616	ALLOCATE(load_tmp1(klon))
	617	ALLOCATE(load_tmp2(klon))
	618	ALLOCATE(load_tmp3(klon))
[2953]	619	ALLOCATE(load_tmp4(klon))
	620	ALLOCATE(load_tmp5(klon))
	621	ALLOCATE(load_tmp6(klon))
	622	ALLOCATE(load_tmp7(klon))
	623	ALLOCATE(load_tmp8(klon))
	624	ALLOCATE(load_tmp9(klon))
	625	ALLOCATE(load_tmp10(klon))
[1279]	626
[1539]	627	!IM ajout variables CFMIP2/CMIP5
	628	ALLOCATE(topswad_aerop(klon), solswad_aerop(klon))
	629	ALLOCATE(topswai_aerop(klon), solswai_aerop(klon))
	630	ALLOCATE(topswad0_aerop(klon), solswad0_aerop(klon))
	631	ALLOCATE(topsw_aerop(klon,naero_grp), topsw0_aerop(klon,naero_grp))
	632	ALLOCATE(solsw_aerop(klon,naero_grp), solsw0_aerop(klon,naero_grp))
	633	ALLOCATE(topswcf_aerop(klon,naero_grp), solswcf_aerop(klon,naero_grp))
	634
[2146]	635	! additional LW variables CK
	636	ALLOCATE(toplwad_aerop(klon), sollwad_aerop(klon))
	637	ALLOCATE(toplwai_aerop(klon), sollwai_aerop(klon))
	638	ALLOCATE(toplwad0_aerop(klon), sollwad0_aerop(klon))
	639
[2897]	640	! FH Ajout de celles n??cessaires au phys_output_write_mod
[1797]	641
[2385]	642	ALLOCATE(tal1(klon), pal1(klon), pab1(klon), pab2(klon))
	643	ALLOCATE(ptstar(klon),pt0(klon),slp(klon))
[2635]	644	!!
	645	!! Wake variables
	646	ALLOCATE(ale_wake(klon), alp_wake(klon))
[2827]	647	ale_wake(:)=0.
[2732]	648	ALLOCATE(wake_h(klon),wake_k(klon))
[2635]	649	ALLOCATE(wake_omg(klon, klev))
	650	ALLOCATE(d_deltat_wk(klon, klev), d_deltaq_wk(klon, klev))
	651	ALLOCATE(d_s_wk(klon), d_dens_wk(klon))
	652	ALLOCATE(d_deltat_wk_gw(klon, klev), d_deltaq_wk_gw(klon, klev))
	653	ALLOCATE(d_deltat_vdf(klon, klev), d_deltaq_vdf(klon, klev))
	654	!! ALLOCATE( d_s_vdf(klon), d_dens_vdf(klon))
	655	ALLOCATE(d_deltat_the(klon, klev), d_deltaq_the(klon, klev))
	656	!! ALLOCATE( d_s_the(klon), d_dens_the(klon))
	657	ALLOCATE(d_deltat_ajs_cv(klon, klev), d_deltaq_ajs_cv(klon, klev))
	658	!! End of wake variables
	659	!!
	660	ALLOCATE(bils(klon))
[1797]	661	ALLOCATE(cdragm(klon), cdragh(klon), cldh(klon), cldl(klon))
[2243]	662	ALLOCATE(cldm(klon), cldq(klon), cldt(klon), qsat2m(klon))
[1924]	663	ALLOCATE(cldhjn(klon), cldljn(klon), cldmjn(klon), cldtjn(klon))
	664	ALLOCATE(JrNt(klon))
[1797]	665	ALLOCATE(dthmin(klon), evap(klon), fder(klon), plcl(klon), plfc(klon))
[2496]	666	ALLOCATE(prw(klon), prlw(klon), prsw(klon), zustar(klon), zu10m(klon), zv10m(klon), rh2m(klon))
[1797]	667	ALLOCATE(s_lcl(klon))
	668	ALLOCATE(s_pblh(klon), s_pblt(klon), s_therm(klon))
[2159]	669	!
	670	!nrlmd+jyg<
	671	ALLOCATE(s_pblh_x(klon), s_pblh_w(klon))
	672	ALLOCATE(s_lcl_x(klon), s_lcl_w(klon))
	673	!>nrlmd+jyg
	674	!
[1797]	675	ALLOCATE(slab_wfbils(klon), tpot(klon), tpote(klon), ue(klon))
	676	ALLOCATE(uq(klon), ve(klon), vq(klon), zxffonte(klon))
[2243]	677	ALLOCATE(zxfqcalving(klon), zxfluxlat(klon))
[1797]	678	ALLOCATE(zxtsol(klon), snow_lsc(klon), zxfqfonte(klon), zxqsurf(klon))
[2518]	679	ALLOCATE(zxrunofflic(klon))
[2159]	680	ALLOCATE(rain_lsc(klon))
[2516]	681	ALLOCATE(rain_num(klon))
[2159]	682	!
	683	ALLOCATE(sens_x(klon), sens_w(klon))
	684	ALLOCATE(zxfluxlat_x(klon), zxfluxlat_w(klon))
[2181]	685	!jyg<
	686	!! ALLOCATE(t_x(klon,klev), t_w(klon,klev))
	687	!! ALLOCATE(q_x(klon,klev), q_w(klon,klev))
	688	!>jyg
[2159]	689	ALLOCATE(dtvdf_x(klon,klev), dtvdf_w(klon,klev))
	690	ALLOCATE(dqvdf_x(klon,klev), dqvdf_w(klon,klev))
	691	ALLOCATE(pbl_tke_input(klon,klev+1,nbsrf))
[2607]	692	ALLOCATE(t_therm(klon,klev), q_therm(klon,klev),u_therm(klon,klev), v_therm(klon,klev))
[2159]	693	ALLOCATE(cdragh_x(klon), cdragh_w(klon))
	694	ALLOCATE(cdragm_x(klon), cdragm_w(klon))
	695	ALLOCATE(kh(klon), kh_x(klon), kh_w(klon))
	696	!
[2877]	697	ALLOCATE(ptconv(klon,klev))
	698	!
[2824]	699	ALLOCATE(wbeff(klon), convoccur(klon), zmax_th(klon))
[1797]	700	ALLOCATE(zq2m(klon), zt2m(klon), weak_inversion(klon))
[2103]	701	ALLOCATE(zt2m_min_mon(klon), zt2m_max_mon(klon))
	702	ALLOCATE(t2m_min_mon(klon), t2m_max_mon(klon))
[1797]	703	ALLOCATE(sens(klon), flwp(klon), fiwp(klon))
	704	ALLOCATE(ale_bl_stat(klon), alp_bl_conv(klon), alp_bl_det(klon))
	705	ALLOCATE(alp_bl_fluct_m(klon), alp_bl_fluct_tke(klon))
	706	ALLOCATE(alp_bl_stat(klon), n2(klon), s2(klon))
	707	ALLOCATE(proba_notrig(klon), random_notrig(klon))
	708
	709	ALLOCATE(dnwd(klon, klev), dnwd0(klon, klev))
	710	! ALLOCATE(upwd(klon, klev), omega(klon, klev), coefh(klon, klev))
	711	ALLOCATE(upwd(klon, klev), omega(klon, klev))
[2481]	712	ALLOCATE(epmax_diag(klon)) ! epmax_cape
	713	ALLOCATE(ep(klon,klev)) ! epmax_cape
[1797]	714	! ALLOCATE(coefm(klon, klev), lambda_th(klon, klev), cldemi(klon, klev))
	715	ALLOCATE(lambda_th(klon, klev), cldemi(klon, klev))
	716	ALLOCATE(cldfra(klon, klev), cldtau(klon, klev), fiwc(klon, klev))
	717	ALLOCATE(fl(klon, klev), re(klon, klev), flwc(klon, klev))
	718	ALLOCATE(ref_liq(klon, klev), ref_ice(klon, klev), theta(klon, klev))
[1989]	719	ALLOCATE(ref_liq_pi(klon, klev), ref_ice_pi(klon, klev))
[2635]	720	ALLOCATE(zphi(klon, klev), zx_rh(klon, klev))
[1797]	721	ALLOCATE(pmfd(klon, klev), pmfu(klon, klev))
	722
[2243]	723	ALLOCATE(t2m(klon, nbsrf), fluxlat(klon, nbsrf))
	724	ALLOCATE(fsollw(klon, nbsrf))
[1797]	725	ALLOCATE(fsolsw(klon, nbsrf), wfbils(klon, nbsrf), wfbilo(klon, nbsrf))
[2670]	726	ALLOCATE(wfevap(klon, nbsrf), wfrain(klon,nbsrf), wfsnow(klon, nbsrf))
[1797]	727	ALLOCATE(evap_pot(klon, nbsrf))
	728
	729	ALLOCATE(pmflxr(klon, klev+1), pmflxs(klon, klev+1), prfl(klon, klev+1))
	730	ALLOCATE(psfl(klon, klev+1), fraca(klon, klev+1), Vprecip(klon, klev+1))
	731	ALLOCATE(zw2(klon, klev+1))
	732
	733	ALLOCATE(fluxu(klon, klev, nbsrf), fluxv(klon, klev, nbsrf))
	734	ALLOCATE(fluxt(klon, klev, nbsrf))
	735
	736	ALLOCATE(uwriteSTD(klon,nlevSTD,nfiles), vwriteSTD(klon,nlevSTD,nfiles))
	737	ALLOCATE(wwriteSTD(klon,nlevSTD,nfiles), phiwriteSTD(klon,nlevSTD,nfiles))
	738	ALLOCATE(qwriteSTD(klon,nlevSTD,nfiles), twriteSTD(klon,nlevSTD,nfiles))
[1828]	739	ALLOCATE(rhwriteSTD(klon,nlevSTD,nfiles))
[1797]	740
	741	! ug et d'autres encore:
	742	ALLOCATE(wdtrainA(klon,klev),wdtrainM(klon,klev))
	743	ALLOCATE(beta_prec(klon,klev))
[2945]	744	ALLOCATE(rneb(klon,klev),rnebjn(klon,klev),rneblsvol(klon,klev))
[1797]	745
	746
[1872]	747	ALLOCATE (zxsnow(klon),snowhgt(klon),qsnow(klon),to_ice(klon))
	748	ALLOCATE (sissnow(klon),runoff(klon),albsol3_lic(klon))
[1865]	749
[2992]	750	ALLOCATE (p_tropopause(klon))
	751	ALLOCATE (z_tropopause(klon))
	752	ALLOCATE (t_tropopause(klon))
	753
[2690]	754	#ifdef CPP_StratAer
	755	ALLOCATE (R2SO4(klon,klev))
	756	ALLOCATE (DENSO4(klon,klev))
	757	ALLOCATE (f_r_wet(klon,klev))
	758	ALLOCATE (decfluxaer(klon,nbtr))
	759	ALLOCATE (mdw(nbtr))
[2752]	760	ALLOCATE (budg_3D_nucl(klon,klev))
	761	ALLOCATE (budg_3D_cond_evap(klon,klev))
	762	ALLOCATE (budg_3D_ocs_to_so2(klon,klev))
	763	ALLOCATE (budg_3D_so2_to_h2so4(klon,klev))
	764	ALLOCATE (budg_3D_backgr_ocs(klon,klev))
	765	ALLOCATE (budg_3D_backgr_so2(klon,klev))
[2690]	766	ALLOCATE (OCS_lifetime(klon,klev))
	767	ALLOCATE (SO2_lifetime(klon,klev))
	768	ALLOCATE (alpha_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave+nwave_lw,nbtr))
	769	ALLOCATE (piz_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave+nwave_lw,nbtr))
	770	ALLOCATE (cg_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave+nwave_lw,nbtr))
	771	ALLOCATE (tau_strat_550(klon,klev))
	772	ALLOCATE (tau_strat_1020(klon,klev))
	773	ALLOCATE (tausum_strat(klon,3))
[2752]	774	ALLOCATE (budg_dep_dry_h2so4(klon))
	775	ALLOCATE (budg_dep_wet_h2so4(klon))
	776	ALLOCATE (budg_dep_dry_part(klon))
	777	ALLOCATE (budg_dep_wet_part(klon))
	778	ALLOCATE (budg_emi_ocs(klon))
	779	ALLOCATE (budg_emi_so2(klon))
	780	ALLOCATE (budg_emi_h2so4(klon))
	781	ALLOCATE (budg_emi_part(klon))
	782	ALLOCATE (budg_ocs_to_so2(klon))
	783	ALLOCATE (budg_so2_to_h2so4(klon))
	784	ALLOCATE (budg_h2so4_to_part(klon))
	785	ALLOCATE (budg_sed_part(klon))
[2690]	786	ALLOCATE (surf_PM25_sulf(klon))
	787	ALLOCATE (vsed_aer(klon,klev))
	788	#endif
[1865]	789
[904]	790	END SUBROUTINE phys_local_var_init
	791
	792	!======================================================================
	793	SUBROUTINE phys_local_var_end
[1785]	794	USE dimphy
	795	USE indice_sol_mod
[904]	796	IMPLICIT NONE
[2854]	797	DEALLOCATE(t_seri,q_seri,ql_seri,qs_seri)
	798	DEALLOCATE(u_seri,v_seri)
	799	DEALLOCATE(l_mixmin,l_mix)
[909]	800
[2854]	801	DEALLOCATE(tr_seri)
	802	DEALLOCATE(d_t_dyn,d_q_dyn)
	803	DEALLOCATE(d_ql_dyn,d_qs_dyn)
	804	DEALLOCATE(d_q_dyn2d,d_ql_dyn2d,d_qs_dyn2d)
	805	DEALLOCATE(d_u_dyn,d_v_dyn)
	806	DEALLOCATE(d_tr_dyn) !RomP
	807	DEALLOCATE(d_t_con,d_q_con)
	808	DEALLOCATE(d_u_con,d_v_con)
	809	DEALLOCATE(d_t_wake,d_q_wake)
	810	DEALLOCATE(d_t_lsc,d_q_lsc)
	811	DEALLOCATE(d_t_lwr,d_t_lw0)
	812	DEALLOCATE(d_t_swr,d_t_sw0)
	813	DEALLOCATE(d_ql_lsc,d_qi_lsc)
	814	DEALLOCATE(d_t_ajsb,d_q_ajsb)
	815	DEALLOCATE(d_t_ajs,d_q_ajs)
[2159]	816	!nrlmd<
[2854]	817	DEALLOCATE(d_t_ajs_w,d_q_ajs_w)
	818	DEALLOCATE(d_t_ajs_x,d_q_ajs_x)
[2159]	819	!>nrlmd
[2854]	820	DEALLOCATE(d_u_ajs,d_v_ajs)
	821	DEALLOCATE(d_t_eva,d_q_eva)
	822	DEALLOCATE(d_ql_eva,d_qi_eva)
	823	DEALLOCATE(d_t_lscst,d_q_lscst)
	824	DEALLOCATE(d_t_lscth,d_q_lscth)
	825	DEALLOCATE(plul_st,plul_th)
	826	DEALLOCATE(d_t_vdf,d_q_vdf,d_t_diss)
[2159]	827	!nrlmd+jyg<
[2854]	828	DEALLOCATE(d_t_vdf_w,d_q_vdf_w)
	829	DEALLOCATE(d_t_vdf_x,d_q_vdf_x)
[2159]	830	!>nrlmd+jyg
[2854]	831	DEALLOCATE(d_u_vdf,d_v_vdf)
	832	DEALLOCATE(d_t_oli,d_t_oro)
	833	DEALLOCATE(d_u_oli,d_v_oli)
	834	DEALLOCATE(d_u_oro,d_v_oro)
[2897]	835	DEALLOCATE(d_t_oro_gw)
	836	DEALLOCATE(d_u_oro_gw,d_v_oro_gw)
[2854]	837	DEALLOCATE(d_t_lif,d_t_ec)
	838	DEALLOCATE(d_u_lif,d_v_lif)
	839	DEALLOCATE(d_ts, d_tr)
	840	DEALLOCATE(topswad_aero,solswad_aero)
	841	DEALLOCATE(topswai_aero,solswai_aero)
	842	DEALLOCATE(topswad0_aero,solswad0_aero)
	843	DEALLOCATE(toplwad_aero,sollwad_aero)
	844	DEALLOCATE(toplwai_aero,sollwai_aero)
	845	DEALLOCATE(toplwad0_aero,sollwad0_aero)
	846	DEALLOCATE(topsw_aero,solsw_aero)
	847	DEALLOCATE(topsw0_aero,solsw0_aero)
	848	DEALLOCATE(topswcf_aero,solswcf_aero)
	849	DEALLOCATE(stratomask)
	850	DEALLOCATE(tausum_aero)
	851	DEALLOCATE(drytausum_aero)
	852	DEALLOCATE(tau3d_aero)
	853	DEALLOCATE(scdnc)
	854	DEALLOCATE(cldncl)
	855	DEALLOCATE(reffclwtop)
	856	DEALLOCATE(lcc)
	857	DEALLOCATE(reffclws)
	858	DEALLOCATE(reffclwc)
	859	DEALLOCATE(cldnvi)
	860	DEALLOCATE(lcc3d)
	861	DEALLOCATE(lcc3dcon)
	862	DEALLOCATE(lcc3dstra)
	863	DEALLOCATE(od443aer)
	864	DEALLOCATE(od550aer)
	865	DEALLOCATE(od865aer)
	866	DEALLOCATE(dryod550aer)
	867	DEALLOCATE(abs550aer)
	868	DEALLOCATE(ec550aer)
	869	DEALLOCATE(od550lt1aer)
	870	DEALLOCATE(sconcso4)
	871	DEALLOCATE(sconcno3)
	872	DEALLOCATE(sconcoa)
	873	DEALLOCATE(sconcbc)
	874	DEALLOCATE(sconcss)
	875	DEALLOCATE(sconcdust)
	876	DEALLOCATE(concso4)
	877	DEALLOCATE(concno3)
	878	DEALLOCATE(concoa)
	879	DEALLOCATE(concbc)
	880	DEALLOCATE(concss)
	881	DEALLOCATE(concdust)
	882	DEALLOCATE(loadso4)
	883	DEALLOCATE(loadoa)
	884	DEALLOCATE(loadbc)
	885	DEALLOCATE(loadss)
	886	DEALLOCATE(loaddust)
	887	DEALLOCATE(loadno3)
	888	DEALLOCATE(load_tmp1)
	889	DEALLOCATE(load_tmp2)
	890	DEALLOCATE(load_tmp3)
[2953]	891	DEALLOCATE(load_tmp4)
	892	DEALLOCATE(load_tmp5)
	893	DEALLOCATE(load_tmp6)
	894	DEALLOCATE(load_tmp7)
	895	DEALLOCATE(load_tmp8)
	896	DEALLOCATE(load_tmp9)
	897	DEALLOCATE(load_tmp10)
[2854]	898	DEALLOCATE(du_gwd_hines,dv_gwd_hines,d_t_hin)
	899	DEALLOCATE(d_q_ch4)
	900	DEALLOCATE(dv_gwd_rando,dv_gwd_front)
	901	DEALLOCATE(east_gwstress,west_gwstress)
[1279]	902
[1539]	903	!IM ajout variables CFMIP2/CMIP5
[2854]	904	DEALLOCATE(topswad_aerop, solswad_aerop)
	905	DEALLOCATE(topswai_aerop, solswai_aerop)
	906	DEALLOCATE(topswad0_aerop, solswad0_aerop)
	907	DEALLOCATE(topsw_aerop, topsw0_aerop)
	908	DEALLOCATE(solsw_aerop, solsw0_aerop)
	909	DEALLOCATE(topswcf_aerop, solswcf_aerop)
[1539]	910
[2146]	911	!CK LW diagnostics
[2854]	912	DEALLOCATE(toplwad_aerop, sollwad_aerop)
	913	DEALLOCATE(toplwai_aerop, sollwai_aerop)
	914	DEALLOCATE(toplwad0_aerop, sollwad0_aerop)
[1797]	915
[2897]	916	! FH Ajout de celles n??cessaires au phys_output_write_mod
[2385]	917	DEALLOCATE(tal1, pal1, pab1, pab2)
	918	DEALLOCATE(ptstar, pt0, slp)
[2635]	919	!
	920	DEALLOCATE(ale_wake, alp_wake)
[2732]	921	DEALLOCATE(wake_h,wake_k)
[2635]	922	DEALLOCATE(wake_omg)
	923	DEALLOCATE(d_deltat_wk, d_deltaq_wk)
	924	DEALLOCATE(d_s_wk, d_dens_wk)
	925	DEALLOCATE(d_deltat_wk_gw, d_deltaq_wk_gw)
	926	DEALLOCATE(d_deltat_vdf, d_deltaq_vdf)
	927	!! DEALLOCATE( d_s_vdf, d_dens_vdf)
	928	DEALLOCATE(d_deltat_the, d_deltaq_the)
	929	!! DEALLOCATE( d_s_the, d_dens_the)
	930	DEALLOCATE(d_deltat_ajs_cv, d_deltaq_ajs_cv)
	931	!
	932	DEALLOCATE(bils)
[1797]	933	DEALLOCATE(cdragm, cdragh, cldh, cldl)
[2243]	934	DEALLOCATE(cldm, cldq, cldt, qsat2m)
[1924]	935	DEALLOCATE(cldljn, cldmjn, cldhjn, cldtjn, JrNt)
[1797]	936	DEALLOCATE(dthmin, evap, fder, plcl, plfc)
[2496]	937	DEALLOCATE(prw, prlw, prsw, zustar, zu10m, zv10m, rh2m, s_lcl)
[1797]	938	DEALLOCATE(s_pblh, s_pblt, s_therm)
[2159]	939	!
	940	!nrlmd+jyg<
	941	DEALLOCATE(s_pblh_x, s_pblh_w)
	942	DEALLOCATE(s_lcl_x, s_lcl_w)
	943	!>nrlmd+jyg
	944	!
[1797]	945	DEALLOCATE(slab_wfbils, tpot, tpote, ue)
	946	DEALLOCATE(uq, ve, vq, zxffonte)
[2243]	947	DEALLOCATE(zxfqcalving, zxfluxlat)
[2518]	948	DEALLOCATE(zxrunofflic)
[1797]	949	DEALLOCATE(zxtsol, snow_lsc, zxfqfonte, zxqsurf)
[2159]	950	DEALLOCATE(rain_lsc)
[2516]	951	DEALLOCATE(rain_num)
[2159]	952	!
	953	DEALLOCATE(sens_x, sens_w)
	954	DEALLOCATE(zxfluxlat_x, zxfluxlat_w)
[2181]	955	!jyg<
	956	!! DEALLOCATE(t_x, t_w)
	957	!! DEALLOCATE(q_x, q_w)
	958	!>jyg
[2159]	959	DEALLOCATE(dtvdf_x, dtvdf_w)
	960	DEALLOCATE(dqvdf_x, dqvdf_w)
	961	DEALLOCATE(pbl_tke_input)
[2607]	962	DEALLOCATE(t_therm, q_therm, u_therm, v_therm)
[2159]	963	DEALLOCATE(cdragh_x, cdragh_w)
	964	DEALLOCATE(cdragm_x, cdragm_w)
	965	DEALLOCATE(kh, kh_x, kh_w)
	966	!
[2877]	967	DEALLOCATE(ptconv)
	968	!
[2824]	969	DEALLOCATE(wbeff, convoccur, zmax_th)
[1797]	970	DEALLOCATE(zq2m, zt2m, weak_inversion)
[2103]	971	DEALLOCATE(zt2m_min_mon, zt2m_max_mon)
	972	DEALLOCATE(t2m_min_mon, t2m_max_mon)
[1797]	973	DEALLOCATE(sens, flwp, fiwp)
	974	DEALLOCATE(ale_bl_stat,alp_bl_conv,alp_bl_det)
	975	DEALLOCATE(alp_bl_fluct_m,alp_bl_fluct_tke)
	976	DEALLOCATE(alp_bl_stat, n2, s2)
	977	DEALLOCATE(proba_notrig, random_notrig)
	978
	979	DEALLOCATE(dnwd, dnwd0)
	980	! DEALLOCATE(upwd, omega, coefh)
	981	DEALLOCATE(upwd, omega)
[2481]	982	DEALLOCATE(epmax_diag)
	983	DEALLOCATE(ep)
[1797]	984	! DEALLOCATE(coefm, lambda_th, cldemi)
	985	DEALLOCATE(lambda_th, cldemi)
	986	DEALLOCATE(cldfra, cldtau, fiwc)
	987	DEALLOCATE(fl, re, flwc)
	988	DEALLOCATE(ref_liq, ref_ice, theta)
[1989]	989	DEALLOCATE(ref_liq_pi, ref_ice_pi)
[2635]	990	DEALLOCATE(zphi, zx_rh)
[1797]	991	DEALLOCATE(pmfd, pmfu)
	992
[2243]	993	DEALLOCATE(t2m, fluxlat)
	994	DEALLOCATE(fsollw, evap_pot)
[1797]	995	DEALLOCATE(fsolsw, wfbils, wfbilo)
[2670]	996	DEALLOCATE(wfevap,wfrain,wfsnow)
[1797]	997
	998	DEALLOCATE(pmflxr, pmflxs, prfl)
	999	DEALLOCATE(psfl, fraca, Vprecip)
	1000	DEALLOCATE(zw2)
	1001
	1002	DEALLOCATE(fluxu, fluxv)
	1003	DEALLOCATE(fluxt)
	1004
	1005	DEALLOCATE(uwriteSTD, vwriteSTD)
	1006	DEALLOCATE(wwriteSTD, phiwriteSTD)
[1828]	1007	DEALLOCATE(qwriteSTD, twriteSTD, rhwriteSTD)
[1797]	1008
	1009	! ug et d'autres encore:
	1010	DEALLOCATE(wdtrainA, wdtrainM)
	1011	DEALLOCATE(beta_prec)
	1012	DEALLOCATE(rneb)
[1865]	1013	DEALLOCATE (zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic)
[1797]	1014
[2992]	1015	DEALLOCATE (p_tropopause)
	1016	DEALLOCATE (z_tropopause)
	1017	DEALLOCATE (t_tropopause)
	1018
[2690]	1019	#ifdef CPP_StratAer
	1020	! variables for strat. aerosol CK
	1021	DEALLOCATE (R2SO4)
	1022	DEALLOCATE (DENSO4)
	1023	DEALLOCATE (f_r_wet)
	1024	DEALLOCATE (decfluxaer)
	1025	DEALLOCATE (mdw)
	1026	DEALLOCATE (SO2_lifetime)
	1027	DEALLOCATE (OCS_lifetime)
	1028	DEALLOCATE (alpha_bin)
	1029	DEALLOCATE (piz_bin)
	1030	DEALLOCATE (cg_bin)
	1031	DEALLOCATE (tau_strat_550)
	1032	DEALLOCATE (tau_strat_1020)
	1033	DEALLOCATE (tausum_strat)
	1034	DEALLOCATE (surf_PM25_sulf)
	1035	DEALLOCATE (vsed_aer)
[2752]	1036	DEALLOCATE (budg_3D_ocs_to_so2)
	1037	DEALLOCATE (budg_3D_so2_to_h2so4)
	1038	DEALLOCATE (budg_3D_backgr_ocs)
	1039	DEALLOCATE (budg_3D_backgr_so2)
	1040	DEALLOCATE (budg_3D_nucl)
	1041	DEALLOCATE (budg_3D_cond_evap)
	1042	DEALLOCATE (budg_dep_dry_h2so4)
	1043	DEALLOCATE (budg_dep_wet_h2so4)
	1044	DEALLOCATE (budg_dep_dry_part)
	1045	DEALLOCATE (budg_dep_wet_part)
	1046	DEALLOCATE (budg_emi_ocs)
	1047	DEALLOCATE (budg_emi_so2)
	1048	DEALLOCATE (budg_emi_h2so4)
	1049	DEALLOCATE (budg_emi_part)
	1050	DEALLOCATE (budg_ocs_to_so2)
	1051	DEALLOCATE (budg_so2_to_h2so4)
	1052	DEALLOCATE (budg_h2so4_to_part)
	1053	DEALLOCATE (budg_sed_part)
[2690]	1054	#endif
	1055
[904]	1056	END SUBROUTINE phys_local_var_end
	1057
[2898]	1058	END MODULE phys_local_var_mod

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

Download in other formats:

Original Format