! ! $Id: phys_local_var_mod.F90 2877 2017-05-12 10:16:53Z evignon $ ! MODULE phys_local_var_mod ! Variables locales pour effectuer les appels en serie !====================================================================== ! ! !====================================================================== ! Declaration des variables REAL, SAVE, ALLOCATABLE :: t_seri(:,:), q_seri(:,:) !$OMP THREADPRIVATE(t_seri, q_seri) REAL, SAVE, ALLOCATABLE :: ql_seri(:,:),qs_seri(:,:) !$OMP THREADPRIVATE(ql_seri,qs_seri) REAL, SAVE, ALLOCATABLE :: u_seri(:,:), v_seri(:,:) !$OMP THREADPRIVATE(u_seri, v_seri) REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:), l_mix(:,:,:) !$OMP THREADPRIVATE(l_mixmin, l_mix) REAL, SAVE, ALLOCATABLE :: tr_seri(:,:,:) !$OMP THREADPRIVATE(tr_seri) REAL, SAVE, ALLOCATABLE :: d_t_dyn(:,:), d_q_dyn(:,:) !$OMP THREADPRIVATE(d_t_dyn, d_q_dyn) REAL, SAVE, ALLOCATABLE :: d_ql_dyn(:,:), d_qs_dyn(:,:) !$OMP THREADPRIVATE(d_ql_dyn, d_qs_dyn) REAL, SAVE, ALLOCATABLE :: d_q_dyn2d(:), d_ql_dyn2d(:), d_qs_dyn2d(:) !$OMP THREADPRIVATE(d_q_dyn2d, d_ql_dyn2d, d_qs_dyn2d) REAL, SAVE, ALLOCATABLE :: d_u_dyn(:,:), d_v_dyn(:,:) !$OMP THREADPRIVATE(d_u_dyn, d_v_dyn) REAL, SAVE, ALLOCATABLE :: d_tr_dyn(:,:,:) !$OMP THREADPRIVATE(d_tr_dyn) REAL, SAVE, ALLOCATABLE :: d_t_con(:,:),d_q_con(:,:) !$OMP THREADPRIVATE(d_t_con,d_q_con) REAL, SAVE, ALLOCATABLE :: d_u_con(:,:),d_v_con(:,:) !$OMP THREADPRIVATE(d_u_con,d_v_con) REAL, SAVE, ALLOCATABLE :: d_t_wake(:,:),d_q_wake(:,:) !$OMP THREADPRIVATE( d_t_wake,d_q_wake) REAL, SAVE, ALLOCATABLE :: d_t_lsc(:,:),d_q_lsc(:,:),d_ql_lsc(:,:),d_qi_lsc(:,:) !$OMP THREADPRIVATE(d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc) REAL, SAVE, ALLOCATABLE :: d_t_lwr(:,:),d_t_lw0(:,:),d_t_swr(:,:),d_t_sw0(:,:) !$OMP THREADPRIVATE(d_t_lwr,d_t_lw0,d_t_swr,d_t_sw0) REAL, SAVE, ALLOCATABLE :: d_t_ajsb(:,:), d_q_ajsb(:,:) !$OMP THREADPRIVATE(d_t_ajsb, d_q_ajsb) REAL, SAVE, ALLOCATABLE :: d_t_ajs(:,:), d_q_ajs(:,:) !$OMP THREADPRIVATE(d_t_ajs, d_q_ajs) REAL, SAVE, ALLOCATABLE :: d_u_ajs(:,:), d_v_ajs(:,:) !$OMP THREADPRIVATE(d_u_ajs, d_v_ajs) !nrlmd< REAL, SAVE, ALLOCATABLE :: d_t_ajs_w(:,:), d_q_ajs_w(:,:) !$OMP THREADPRIVATE(d_t_ajs_w, d_q_ajs_w) REAL, SAVE, ALLOCATABLE :: d_t_ajs_x(:,:), d_q_ajs_x(:,:) !$OMP THREADPRIVATE(d_t_ajs_x, d_q_ajs_x) !>nrlmd REAL, SAVE, ALLOCATABLE :: d_t_eva(:,:),d_q_eva(:,:),d_ql_eva(:,:),d_qi_eva(:,:) !$OMP THREADPRIVATE(d_t_eva,d_q_eva,d_ql_eva,d_qi_eva) REAL, SAVE, ALLOCATABLE :: d_t_lscst(:,:),d_q_lscst(:,:) !$OMP THREADPRIVATE(d_t_lscst,d_q_lscst) REAL, SAVE, ALLOCATABLE :: d_t_lscth(:,:),d_q_lscth(:,:) !$OMP THREADPRIVATE(d_t_lscth,d_q_lscth) REAL, SAVE, ALLOCATABLE :: plul_th(:),plul_st(:) !$OMP THREADPRIVATE(plul_th,plul_st) !tendances dues a oro et lif REAL, SAVE, ALLOCATABLE :: d_t_oli(:,:) !$OMP THREADPRIVATE(d_t_oli) REAL, SAVE, ALLOCATABLE :: d_u_oli(:,:), d_v_oli(:,:) !$OMP THREADPRIVATE(d_u_oli, d_v_oli) REAL, SAVE, ALLOCATABLE :: d_t_vdf(:,:), d_q_vdf(:,:), d_t_diss(:,:) !$OMP THREADPRIVATE( d_t_vdf, d_q_vdf,d_t_diss) REAL, SAVE, ALLOCATABLE :: d_u_vdf(:,:), d_v_vdf(:,:) !$OMP THREADPRIVATE(d_u_vdf, d_v_vdf) !nrlmd+jyg< REAL, SAVE, ALLOCATABLE :: d_t_vdf_w(:,:), d_q_vdf_w(:,:) !$OMP THREADPRIVATE( d_t_vdf_w, d_q_vdf_w) REAL, SAVE, ALLOCATABLE :: d_t_vdf_x(:,:), d_q_vdf_x(:,:) !$OMP THREADPRIVATE( d_t_vdf_x, d_q_vdf_x) !>nrlmd+jyg REAL, SAVE, ALLOCATABLE :: d_t_oro(:,:) !$OMP THREADPRIVATE(d_t_oro) REAL, SAVE, ALLOCATABLE :: d_u_oro(:,:), d_v_oro(:,:) !$OMP THREADPRIVATE(d_u_oro, d_v_oro) REAL, SAVE, ALLOCATABLE :: d_t_lif(:,:) !$OMP THREADPRIVATE(d_t_lif) REAL, SAVE, ALLOCATABLE :: d_u_lif(:,:), d_v_lif(:,:) !$OMP THREADPRIVATE(d_u_lif, d_v_lif) ! Tendances Ondes de G non oro (runs strato). REAL, SAVE, ALLOCATABLE :: du_gwd_hines(:,:) !$OMP THREADPRIVATE(du_gwd_hines) REAL, SAVE, ALLOCATABLE :: dv_gwd_hines(:,:) !$OMP THREADPRIVATE(dv_gwd_hines) REAL, SAVE, ALLOCATABLE :: dv_gwd_rando(:,:) !$OMP THREADPRIVATE(dv_gwd_rando) REAL, SAVE, ALLOCATABLE :: dv_gwd_front(:,:) !$OMP THREADPRIVATE(dv_gwd_front) REAL, SAVE, ALLOCATABLE :: east_gwstress(:,:) !$OMP THREADPRIVATE(east_gwstress) REAL, SAVE, ALLOCATABLE :: west_gwstress(:,:) !$OMP THREADPRIVATE(west_gwstress) REAL, SAVE, ALLOCATABLE :: d_t_hin(:,:) !$OMP THREADPRIVATE(d_t_hin) ! tendance due a l'oxydation du methane REAL, SAVE, ALLOCATABLE :: d_q_ch4(:,:) !$OMP THREADPRIVATE(d_q_ch4) ! tendance du a la conersion Ec -> E thermique REAL, SAVE, ALLOCATABLE :: d_t_ec(:,:) !$OMP THREADPRIVATE(d_t_ec) REAL, SAVE, ALLOCATABLE :: d_ts(:,:), d_tr(:,:,:) !$OMP THREADPRIVATE(d_ts, d_tr) ! diagnostique pour le rayonnement REAL, SAVE, ALLOCATABLE :: topswad_aero(:), solswad_aero(:) ! diag !$OMP THREADPRIVATE(topswad_aero,solswad_aero) REAL, SAVE, ALLOCATABLE :: topswai_aero(:), solswai_aero(:) ! diag !$OMP THREADPRIVATE(topswai_aero,solswai_aero) REAL, SAVE, ALLOCATABLE :: topswad0_aero(:), solswad0_aero(:) ! diag !$OMP THREADPRIVATE(topswad0_aero,solswad0_aero) REAL, SAVE, ALLOCATABLE :: topsw_aero(:,:), solsw_aero(:,:) ! diag !$OMP THREADPRIVATE(topsw_aero,solsw_aero) REAL, SAVE, ALLOCATABLE :: topsw0_aero(:,:), solsw0_aero(:,:) ! diag !$OMP THREADPRIVATE(topsw0_aero,solsw0_aero) REAL, SAVE, ALLOCATABLE :: topswcf_aero(:,:), solswcf_aero(:,:) ! diag !$OMP THREADPRIVATE(topswcf_aero,solswcf_aero) ! LW radiation diagnostics CK REAL, SAVE, ALLOCATABLE :: toplwad_aero(:), sollwad_aero(:) ! diag !$OMP THREADPRIVATE(toplwad_aero,sollwad_aero) REAL, SAVE, ALLOCATABLE :: toplwai_aero(:), sollwai_aero(:) ! diag !$OMP THREADPRIVATE(toplwai_aero,sollwai_aero) REAL, SAVE, ALLOCATABLE :: toplwad0_aero(:), sollwad0_aero(:) ! diag !$OMP THREADPRIVATE(toplwad0_aero,sollwad0_aero) ! Special RRTM REAL, SAVE, ALLOCATABLE :: ZLWFT0_i(:,:), ZSWFT0_i(:,:) ! diag !$OMP THREADPRIVATE(ZLWFT0_i,ZSWFT0_i) REAL, SAVE, ALLOCATABLE :: ZFLDN0(:,:), ZFLUP0(:,:) ! diag !$OMP THREADPRIVATE(ZFLDN0,ZFLUP0) REAL, SAVE, ALLOCATABLE :: ZFSDN0(:,:), ZFSUP0(:,:) ! diag !$OMP THREADPRIVATE(ZFSDN0,ZFSUP0) ! REAL, SAVE, ALLOCATABLE :: stratomask(:,:) !$OMP THREADPRIVATE(stratomask) REAL, SAVE, ALLOCATABLE :: tausum_aero(:,:,:) !$OMP THREADPRIVATE(tausum_aero) REAL, SAVE, ALLOCATABLE :: drytausum_aero(:,:) !$OMP THREADPRIVATE(drytausum_aero) REAL, SAVE, ALLOCATABLE :: tau3d_aero(:,:,:,:) !$OMP THREADPRIVATE(tau3d_aero) REAL, SAVE, ALLOCATABLE :: scdnc(:,:) !$OMP THREADPRIVATE(scdnc) REAL, SAVE, ALLOCATABLE :: cldncl(:) !$OMP THREADPRIVATE(cldncl) REAL, SAVE, ALLOCATABLE :: reffclwtop(:) !$OMP THREADPRIVATE(reffclwtop) REAL, SAVE, ALLOCATABLE :: lcc(:) !$OMP THREADPRIVATE(lcc) REAL, SAVE, ALLOCATABLE :: reffclws(:,:) !$OMP THREADPRIVATE(reffclws) REAL, SAVE, ALLOCATABLE :: reffclwc(:,:) !$OMP THREADPRIVATE(reffclwc) REAL, SAVE, ALLOCATABLE :: cldnvi(:) !$OMP THREADPRIVATE(cldnvi) REAL, SAVE, ALLOCATABLE :: lcc3d(:,:) !$OMP THREADPRIVATE(lcc3d) REAL, SAVE, ALLOCATABLE :: lcc3dcon(:,:) !$OMP THREADPRIVATE(lcc3dcon) REAL, SAVE, ALLOCATABLE :: lcc3dstra(:,:) !$OMP THREADPRIVATE(lcc3dstra) REAL, SAVE, ALLOCATABLE :: od443aer(:) !$OMP THREADPRIVATE(od443aer) REAL, SAVE, ALLOCATABLE :: od550aer(:) !$OMP THREADPRIVATE(od550aer) REAL, SAVE, ALLOCATABLE :: dryod550aer(:) !$OMP THREADPRIVATE(dryod550aer) REAL, SAVE, ALLOCATABLE :: abs550aer(:) !$OMP THREADPRIVATE(abs550aer) REAL, SAVE, ALLOCATABLE :: od865aer(:) !$OMP THREADPRIVATE(od865aer) REAL, SAVE, ALLOCATABLE :: ec550aer(:,:) !$OMP THREADPRIVATE(ec550aer) REAL, SAVE, ALLOCATABLE :: od550lt1aer(:) !$OMP THREADPRIVATE(od550lt1aer) REAL, SAVE, ALLOCATABLE :: sconcso4(:) !$OMP THREADPRIVATE(sconcso4) REAL, SAVE, ALLOCATABLE :: sconcno3(:) !$OMP THREADPRIVATE(sconcno3) REAL, SAVE, ALLOCATABLE :: sconcoa(:) !$OMP THREADPRIVATE(sconcoa) REAL, SAVE, ALLOCATABLE :: sconcbc(:) !$OMP THREADPRIVATE(sconcbc) REAL, SAVE, ALLOCATABLE :: sconcss(:) !$OMP THREADPRIVATE(sconcss) REAL, SAVE, ALLOCATABLE :: sconcdust(:) !$OMP THREADPRIVATE(sconcdust) REAL, SAVE, ALLOCATABLE :: concso4(:,:) !$OMP THREADPRIVATE(concso4) REAL, SAVE, ALLOCATABLE :: concno3(:,:) !$OMP THREADPRIVATE(concno3) REAL, SAVE, ALLOCATABLE :: concoa(:,:) !$OMP THREADPRIVATE(concoa) REAL, SAVE, ALLOCATABLE :: concbc(:,:) !$OMP THREADPRIVATE(concbc) REAL, SAVE, ALLOCATABLE :: concss(:,:) !$OMP THREADPRIVATE(concss) REAL, SAVE, ALLOCATABLE :: concdust(:,:) !$OMP THREADPRIVATE(concdust) REAL, SAVE, ALLOCATABLE :: loadso4(:) !$OMP THREADPRIVATE(loadso4) REAL, SAVE, ALLOCATABLE :: loadoa(:) !$OMP THREADPRIVATE(loadoa) REAL, SAVE, ALLOCATABLE :: loadbc(:) !$OMP THREADPRIVATE(loadbc) REAL, SAVE, ALLOCATABLE :: loadss(:) !$OMP THREADPRIVATE(loadss) REAL, SAVE, ALLOCATABLE :: loaddust(:) !$OMP THREADPRIVATE(loaddust) REAL, SAVE, ALLOCATABLE :: loadno3(:) !$OMP THREADPRIVATE(loadno3) REAL, SAVE, ALLOCATABLE :: load_tmp1(:) !$OMP THREADPRIVATE(load_tmp1) REAL, SAVE, ALLOCATABLE :: load_tmp2(:) !$OMP THREADPRIVATE(load_tmp2) REAL, SAVE, ALLOCATABLE :: load_tmp3(:) !$OMP THREADPRIVATE(load_tmp3) !IM ajout variables CFMIP2/CMIP5 REAL,ALLOCATABLE,SAVE :: topswad_aerop(:), solswad_aerop(:) !$OMP THREADPRIVATE(topswad_aerop, solswad_aerop) REAL,ALLOCATABLE,SAVE :: topswai_aerop(:), solswai_aerop(:) !$OMP THREADPRIVATE(topswai_aerop, solswai_aerop) REAL,ALLOCATABLE,SAVE :: topswad0_aerop(:), solswad0_aerop(:) !$OMP THREADPRIVATE(topswad0_aerop, solswad0_aerop) REAL,ALLOCATABLE,SAVE :: topsw_aerop(:,:), topsw0_aerop(:,:) !$OMP THREADPRIVATE(topsw_aerop, topsw0_aerop) REAL,ALLOCATABLE,SAVE :: solsw_aerop(:,:), solsw0_aerop(:,:) !$OMP THREADPRIVATE(solsw_aerop, solsw0_aerop) REAL,ALLOCATABLE,SAVE :: topswcf_aerop(:,:), solswcf_aerop(:,:) !$OMP THREADPRIVATE(topswcf_aerop, solswcf_aerop) ! additional LW variables CK REAL,ALLOCATABLE,SAVE :: toplwad_aerop(:), sollwad_aerop(:) !$OMP THREADPRIVATE(toplwad_aerop, sollwad_aerop) REAL,ALLOCATABLE,SAVE :: toplwai_aerop(:), sollwai_aerop(:) !$OMP THREADPRIVATE(toplwai_aerop, sollwai_aerop) REAL,ALLOCATABLE,SAVE :: toplwad0_aerop(:), sollwad0_aerop(:) !$OMP THREADPRIVATE(toplwad0_aerop, sollwad0_aerop) !Ajout de celles nécessaires au phys_output_write_mod REAL, SAVE, ALLOCATABLE :: tal1(:), pal1(:), pab1(:), pab2(:) !$OMP THREADPRIVATE(tal1, pal1, pab1, pab2) REAL, SAVE, ALLOCATABLE :: ptstar(:), pt0(:), slp(:) !$OMP THREADPRIVATE(ptstar, pt0, slp) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: sens, flwp, fiwp !$OMP THREADPRIVATE(sens, flwp, fiwp) !! !! Wake variables REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: ale_wake, alp_wake !$OMP THREADPRIVATE(ale_wake, alp_wake) !!jyg! REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wake_h,wake_k REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wake_h INTEGER,ALLOCATABLE,SAVE,DIMENSION(:) :: wake_k !$OMP THREADPRIVATE(wake_h,wake_k) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wake_omg !$OMP THREADPRIVATE(wake_omg) REAL, SAVE, ALLOCATABLE,DIMENSION(:,:) :: d_deltat_wk, d_deltaq_wk !$OMP THREADPRIVATE(d_deltat_wk, d_deltaq_wk) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: d_s_wk, d_dens_wk !$OMP THREADPRIVATE(d_s_wk, d_dens_wk) REAL, SAVE, ALLOCATABLE,DIMENSION(:,:) :: d_deltat_wk_gw, d_deltaq_wk_gw !$OMP THREADPRIVATE(d_deltat_wk_gw, d_deltaq_wk_gw) REAL, SAVE, ALLOCATABLE,DIMENSION(:,:) :: d_deltat_vdf, d_deltaq_vdf !$OMP THREADPRIVATE(d_deltat_vdf, d_deltaq_vdf) !!! REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: d_s_vdf, d_dens_vdf !!!OMP THREADPRIVATE(d_s_vdf, d_dens_vdf) REAL, SAVE, ALLOCATABLE,DIMENSION(:,:) :: d_deltat_the, d_deltaq_the !$OMP THREADPRIVATE(d_deltat_the, d_deltaq_the) !!! REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: d_s_the, d_dens_the !!!OMP THREADPRIVATE(d_s_the, d_dens_the) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: d_deltat_ajs_cv, d_deltaq_ajs_cv !$OMP THREADPRIVATE(d_deltat_ajs_cv, d_deltaq_ajs_cv) !! End of Wake variables !! REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: bils !$OMP THREADPRIVATE(bils) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cdragm, cdragh !$OMP THREADPRIVATE(cdragm, cdragh) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cldh, cldl, cldm, cldq, cldt, qsat2m !$OMP THREADPRIVATE(cldh, cldl, cldm, cldq, cldt, qsat2m ) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cldhjn, cldljn, cldmjn,cldtjn !$OMP THREADPRIVATE(cldhjn, cldljn, cldmjn, cldtjn) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: JrNt !$OMP THREADPRIVATE(JrNt) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: dthmin, evap, fder, plcl, plfc, prw, prlw, prsw !$OMP THREADPRIVATE(dthmin, evap, fder, plcl, plfc, prw, prlw, prsw) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zustar, zu10m, zv10m, rh2m !$OMP THREADPRIVATE(zustar, zu10m, zv10m, rh2m) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: s_lcl, s_pblh, s_pblt, s_therm !$OMP THREADPRIVATE(s_lcl, s_pblh, s_pblt, s_therm) ! !nrlmd+jyg< REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: s_pblh_x, s_pblh_w !$OMP THREADPRIVATE(s_pblh_x, s_pblh_w) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: s_lcl_x, s_lcl_w !$OMP THREADPRIVATE(s_lcl_x, s_lcl_w) !>nrlmd+jyg ! REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: slab_wfbils !$OMP THREADPRIVATE(slab_wfbils) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: tpot, tpote, ue, uq, ve, vq, zxffonte !$OMP THREADPRIVATE(tpot, tpote, ue, uq, ve, vq, zxffonte) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfqcalving !$OMP THREADPRIVATE(zxfqcalving) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfluxlat, zxtsol, snow_lsc, zxfqfonte !$OMP THREADPRIVATE(zxfluxlat, zxtsol, snow_lsc, zxfqfonte) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxrunofflic !$OMP THREADPRIVATE(zxrunofflic) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxqsurf, rain_lsc, rain_num !$OMP THREADPRIVATE(zxqsurf, rain_lsc, rain_num) ! !jyg+nrlmd< !!!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc ! c ! Declarations liees a la couche limite differentiee w-x c ! c !!!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: sens_x, sens_w !$OMP THREADPRIVATE(sens_x, sens_w) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfluxlat_x, zxfluxlat_w !$OMP THREADPRIVATE(zxfluxlat_x, zxfluxlat_w) !jyg< !!! Entrées supplémentaires couche-limite !! REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: t_x, t_w !!!$OMP THREADPRIVATE(t_x, t_w) !! REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: q_x, q_w !!!$OMP THREADPRIVATE(q_x, q_w) !>jyg !!! Sorties ferret REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dtvdf_x, dtvdf_w !$OMP THREADPRIVATE(dtvdf_x, dtvdf_w) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dqvdf_x, dqvdf_w !$OMP THREADPRIVATE(dqvdf_x, dqvdf_w) ! Variables supplémentaires dans physiq.F relative au splitting de la surface REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: pbl_tke_input !$OMP THREADPRIVATE(pbl_tke_input) ! Entree supplementaire Thermiques : REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: t_therm, q_therm, u_therm, v_therm !$OMP THREADPRIVATE(t_therm, q_therm, u_therm, v_therm) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cdragh_x, cdragh_w !$OMP THREADPRIVATE(cdragh_x, cdragh_w) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cdragm_x, cdragm_w !$OMP THREADPRIVATE(cdragm_x, cdragm_w) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: kh, kh_x, kh_w !$OMP THREADPRIVATE(kh, kh_x, kh_w) !!! !!!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc LOGICAL, SAVE, ALLOCATABLE :: ptconv(:,:) !$OMP THREADPRIVATE(ptconv) !>jyg+nrlmd ! REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wbeff, convoccur, zmax_th, zq2m, zt2m !$OMP THREADPRIVATE(wbeff, convoccur, zmax_th, zq2m, zt2m) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zt2m_min_mon, zt2m_max_mon !$OMP THREADPRIVATE(zt2m_min_mon, zt2m_max_mon) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: t2m_min_mon, t2m_max_mon !$OMP THREADPRIVATE(t2m_min_mon, t2m_max_mon) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: weak_inversion !$OMP THREADPRIVATE(weak_inversion) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: ale_bl_stat,alp_bl_conv,alp_bl_det !$OMP THREADPRIVATE(ale_bl_stat,alp_bl_conv,alp_bl_det) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_fluct_m,alp_bl_fluct_tke !$OMP THREADPRIVATE(alp_bl_fluct_m,alp_bl_fluct_tke) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_stat, n2, s2 !$OMP THREADPRIVATE(alp_bl_stat, n2, s2) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: proba_notrig, random_notrig !$OMP THREADPRIVATE(proba_notrig, random_notrig) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: fsolsw, wfbils, wfbilo !$OMP THREADPRIVATE(fsolsw, wfbils, wfbilo) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wfevap, wfrain, wfsnow !$OMP THREADPRIVATE(wfevap,wfrain,wfsnow) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: t2m, fluxlat, fsollw,evap_pot !$OMP THREADPRIVATE(t2m, fluxlat, fsollw,evap_pot) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dnwd, dnwd0, upwd, omega !$OMP THREADPRIVATE(dnwd, dnwd0, upwd, omega) REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: epmax_diag ! epmax_cape !$OMP THREADPRIVATE(epmax_diag) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ep ! epmax_cape !$OMP THREADPRIVATE(ep) ! REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: coefh, coefm, lambda_th REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: lambda_th !$OMP THREADPRIVATE(lambda_th) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: cldemi, cldfra, cldtau, fiwc, fl, re, flwc !$OMP THREADPRIVATE(cldemi, cldfra, cldtau, fiwc, fl, re, flwc) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq, ref_ice, theta, zphi !$OMP THREADPRIVATE(ref_liq, ref_ice, theta, zphi) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq_pi, ref_ice_pi !$OMP THREADPRIVATE(ref_liq_pi, ref_ice_pi) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zx_rh !$OMP THREADPRIVATE(zx_rh) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: pmflxr, pmflxs, prfl, psfl, fraca !$OMP THREADPRIVATE(pmflxr, pmflxs, prfl, psfl, fraca) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: Vprecip, zw2 !$OMP THREADPRIVATE(Vprecip, zw2) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: pmfd, pmfu !$OMP THREADPRIVATE(pmfd, pmfu) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: fluxt, fluxu, fluxv !$OMP THREADPRIVATE(fluxt, fluxu, fluxv) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: uwriteSTD, vwriteSTD, wwriteSTD !$OMP THREADPRIVATE(uwriteSTD, vwriteSTD, wwriteSTD) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: phiwriteSTD, qwriteSTD, twriteSTD, rhwriteSTD !$OMP THREADPRIVATE(phiwriteSTD, qwriteSTD, twriteSTD, rhwriteSTD) ! ug et d'autres encore: REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wdtrainA, wdtrainM !$OMP THREADPRIVATE(wdtrainA, wdtrainM) REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: beta_prec !$OMP THREADPRIVATE(beta_prec) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rneb,rnebjn !$OMP THREADPRIVATE(rneb,rnebjn) ! variables de sorties MM REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: zxsnow,snowhgt,qsnow,to_ice REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: sissnow,runoff,albsol3_lic !$OMP THREADPRIVATE(zxsnow,snowhgt,qsnow,to_ice) !$OMP THREADPRIVATE(sissnow,runoff,albsol3_lic) #ifdef CPP_StratAer ! ! variables for stratospheric aerosol REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: R2SO4 !$OMP THREADPRIVATE(R2SO4) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: DENSO4 !$OMP THREADPRIVATE(DENSO4) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: f_r_wet !$OMP THREADPRIVATE(f_r_wet) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: decfluxaer !$OMP THREADPRIVATE(decfluxaer) REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: mdw !$OMP THREADPRIVATE(mdw) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: OCS_lifetime !$OMP THREADPRIVATE(OCS_lifetime) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: SO2_lifetime !$OMP THREADPRIVATE(SO2_lifetime) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: alpha_bin !$OMP THREADPRIVATE(alpha_bin) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: piz_bin !$OMP THREADPRIVATE(piz_bin) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: cg_bin !$OMP THREADPRIVATE(cg_bin) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_550 !$OMP THREADPRIVATE(tau_strat_550) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_1020 !$OMP THREADPRIVATE(tau_strat_1020) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tausum_strat !$OMP THREADPRIVATE(tausum_strat) REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: surf_PM25_sulf !$OMP THREADPRIVATE(surf_PM25_sulf) REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: p_tropopause !$OMP THREADPRIVATE(p_tropopause) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: vsed_aer !$OMP THREADPRIVATE(vsed_aer) ! !---3D budget variables REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_nucl !$OMP THREADPRIVATE(budg_3D_nucl) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_cond_evap !$OMP THREADPRIVATE(budg_3D_cond_evap) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_ocs_to_so2 !$OMP THREADPRIVATE(budg_3D_ocs_to_so2) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_so2_to_h2so4 !$OMP THREADPRIVATE(budg_3D_so2_to_h2so4) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_backgr_ocs !$OMP THREADPRIVATE(budg_3D_backgr_ocs) REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: budg_3D_backgr_so2 !$OMP THREADPRIVATE(budg_3D_backgr_so2) ! !---2D budget variables REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_h2so4 REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_part !$OMP THREADPRIVATE(budg_dep_dry_h2so4,budg_dep_dry_part) REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_h2so4 REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_part !$OMP THREADPRIVATE(budg_dep_wet_h2so4,budg_dep_wet_part) REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_ocs REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_so2 REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_h2so4 REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_part !$OMP THREADPRIVATE(budg_emi_ocs,budg_emi_so2,budg_emi_h2so4,budg_emi_part) REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_ocs_to_so2 REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_so2_to_h2so4 REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_h2so4_to_part !$OMP THREADPRIVATE(budg_ocs_to_so2,budg_so2_to_h2so4,budg_h2so4_to_part) REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_sed_part !$OMP THREADPRIVATE(budg_sed_part) #endif REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: pr_tropopause !$OMP THREADPRIVATE(pr_tropopause) CONTAINS !====================================================================== SUBROUTINE phys_local_var_init USE dimphy USE infotrac_phy, ONLY : nbtr USE aero_mod USE indice_sol_mod USE phys_output_var_mod USE phys_state_var_mod IMPLICIT NONE ALLOCATE(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev)) ALLOCATE(u_seri(klon,klev),v_seri(klon,klev)) ALLOCATE(l_mixmin(klon,klev,nbsrf), l_mix(klon,klev,nbsrf)) l_mix(:,:,:)=0. ; l_mixmin(:,:,:)=0. ! doit etre initialse car pas toujours remplis ALLOCATE(tr_seri(klon,klev,nbtr)) ALLOCATE(d_t_dyn(klon,klev),d_q_dyn(klon,klev)) ALLOCATE(d_ql_dyn(klon,klev),d_qs_dyn(klon,klev)) ALLOCATE(d_q_dyn2d(klon),d_ql_dyn2d(klon),d_qs_dyn2d(klon)) ALLOCATE(d_u_dyn(klon,klev),d_v_dyn(klon,klev)) ALLOCATE(d_tr_dyn(klon,klev,nbtr)) !RomP ALLOCATE(d_t_con(klon,klev),d_q_con(klon,klev)) ALLOCATE(d_u_con(klon,klev),d_v_con(klon,klev)) ALLOCATE(d_t_wake(klon,klev),d_q_wake(klon,klev)) ALLOCATE(d_t_lsc(klon,klev),d_q_lsc(klon,klev)) ALLOCATE(d_t_lwr(klon,klev),d_t_lw0(klon,klev)) ALLOCATE(d_t_swr(klon,klev),d_t_sw0(klon,klev)) ALLOCATE(d_ql_lsc(klon,klev),d_qi_lsc(klon,klev)) ALLOCATE(d_t_ajsb(klon,klev),d_q_ajsb(klon,klev)) ALLOCATE(d_t_ajs(klon,klev),d_q_ajs(klon,klev)) !nrlmd< ALLOCATE(d_t_ajs_w(klon,klev),d_q_ajs_w(klon,klev)) ALLOCATE(d_t_ajs_x(klon,klev),d_q_ajs_x(klon,klev)) !>nrlmd ALLOCATE(d_u_ajs(klon,klev),d_v_ajs(klon,klev)) ALLOCATE(d_t_eva(klon,klev),d_q_eva(klon,klev)) ALLOCATE(d_ql_eva(klon,klev),d_qi_eva(klon,klev)) ALLOCATE(d_t_lscst(klon,klev),d_q_lscst(klon,klev)) ALLOCATE(d_t_lscth(klon,klev),d_q_lscth(klon,klev)) ALLOCATE(plul_st(klon),plul_th(klon)) ALLOCATE(d_t_vdf(klon,klev),d_q_vdf(klon,klev),d_t_diss(klon,klev)) !nrlmd+jyg< ALLOCATE(d_t_vdf_w(klon,klev),d_q_vdf_w(klon,klev)) ALLOCATE(d_t_vdf_x(klon,klev),d_q_vdf_x(klon,klev)) !>nrlmd+jyg ALLOCATE(d_u_vdf(klon,klev),d_v_vdf(klon,klev)) ALLOCATE(d_t_oli(klon,klev),d_t_oro(klon,klev)) ALLOCATE(d_u_oli(klon,klev),d_v_oli(klon,klev)) ALLOCATE(d_u_oro(klon,klev),d_v_oro(klon,klev)) ALLOCATE(d_t_lif(klon,klev),d_t_ec(klon,klev)) ALLOCATE(d_u_lif(klon,klev),d_v_lif(klon,klev)) ALLOCATE(d_ts(klon,nbsrf), d_tr(klon,klev,nbtr)) ! Special RRTM ALLOCATE(ZLWFT0_i(klon,klev+1),ZSWFT0_i(klon,klev+1),ZFLDN0(klon,klev+1)) ALLOCATE(ZFLUP0(klon,klev+1),ZFSDN0(klon,klev+1),ZFSUP0(klon,klev+1)) ! ALLOCATE(topswad_aero(klon), solswad_aero(klon)) ALLOCATE(topswai_aero(klon), solswai_aero(klon)) ALLOCATE(topswad0_aero(klon), solswad0_aero(klon)) ALLOCATE(toplwad_aero(klon), sollwad_aero(klon)) ALLOCATE(toplwai_aero(klon), sollwai_aero(klon)) ALLOCATE(toplwad0_aero(klon), sollwad0_aero(klon)) ALLOCATE(topsw_aero(klon,naero_grp), solsw_aero(klon,naero_grp)) ALLOCATE(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp)) ALLOCATE(topswcf_aero(klon,3), solswcf_aero(klon,3)) ALLOCATE(du_gwd_hines(klon,klev),dv_gwd_hines(klon,klev)) ALLOCATE(dv_gwd_rando(klon,klev),dv_gwd_front(klon,klev)) ALLOCATE(east_gwstress(klon,klev),west_gwstress(klon,klev)) ALLOCATE(d_t_hin(klon,klev)) ALLOCATE(d_q_ch4(klon,klev)) ALLOCATE(stratomask(klon,klev)) ALLOCATE(tausum_aero(klon,nwave,naero_tot)) ALLOCATE(drytausum_aero(klon,naero_tot)) ALLOCATE(tau3d_aero(klon,klev,nwave,naero_tot)) ALLOCATE(scdnc(klon, klev)) ALLOCATE(cldncl(klon)) ALLOCATE(reffclwtop(klon)) ALLOCATE(lcc(klon)) ALLOCATE(reffclws(klon, klev)) ALLOCATE(reffclwc(klon, klev)) ALLOCATE(cldnvi(klon)) ALLOCATE(lcc3d(klon, klev)) ALLOCATE(lcc3dcon(klon, klev)) ALLOCATE(lcc3dstra(klon, klev)) ALLOCATE(od443aer(klon)) ALLOCATE(od550aer(klon)) ALLOCATE(od865aer(klon)) ALLOCATE(dryod550aer(klon)) ALLOCATE(abs550aer(klon)) ALLOCATE(ec550aer(klon,klev)) ALLOCATE(od550lt1aer(klon)) ALLOCATE(sconcso4(klon)) ALLOCATE(sconcno3(klon)) ALLOCATE(sconcoa(klon)) ALLOCATE(sconcbc(klon)) ALLOCATE(sconcss(klon)) ALLOCATE(sconcdust(klon)) ALLOCATE(concso4(klon,klev)) ALLOCATE(concno3(klon,klev)) ALLOCATE(concoa(klon,klev)) ALLOCATE(concbc(klon,klev)) ALLOCATE(concss(klon,klev)) ALLOCATE(concdust(klon,klev)) ALLOCATE(loadso4(klon)) ALLOCATE(loadoa(klon)) ALLOCATE(loadbc(klon)) ALLOCATE(loadss(klon)) ALLOCATE(loaddust(klon)) ALLOCATE(loadno3(klon)) ALLOCATE(load_tmp1(klon)) ALLOCATE(load_tmp2(klon)) ALLOCATE(load_tmp3(klon)) !IM ajout variables CFMIP2/CMIP5 ALLOCATE(topswad_aerop(klon), solswad_aerop(klon)) ALLOCATE(topswai_aerop(klon), solswai_aerop(klon)) ALLOCATE(topswad0_aerop(klon), solswad0_aerop(klon)) ALLOCATE(topsw_aerop(klon,naero_grp), topsw0_aerop(klon,naero_grp)) ALLOCATE(solsw_aerop(klon,naero_grp), solsw0_aerop(klon,naero_grp)) ALLOCATE(topswcf_aerop(klon,naero_grp), solswcf_aerop(klon,naero_grp)) ! additional LW variables CK ALLOCATE(toplwad_aerop(klon), sollwad_aerop(klon)) ALLOCATE(toplwai_aerop(klon), sollwai_aerop(klon)) ALLOCATE(toplwad0_aerop(klon), sollwad0_aerop(klon)) ! FH Ajout de celles nécessaires au phys_output_write_mod ALLOCATE(tal1(klon), pal1(klon), pab1(klon), pab2(klon)) ALLOCATE(ptstar(klon),pt0(klon),slp(klon)) !! !! Wake variables ALLOCATE(ale_wake(klon), alp_wake(klon)) ale_wake(:)=0. ALLOCATE(wake_h(klon),wake_k(klon)) ALLOCATE(wake_omg(klon, klev)) ALLOCATE(d_deltat_wk(klon, klev), d_deltaq_wk(klon, klev)) ALLOCATE(d_s_wk(klon), d_dens_wk(klon)) ALLOCATE(d_deltat_wk_gw(klon, klev), d_deltaq_wk_gw(klon, klev)) ALLOCATE(d_deltat_vdf(klon, klev), d_deltaq_vdf(klon, klev)) !! ALLOCATE( d_s_vdf(klon), d_dens_vdf(klon)) ALLOCATE(d_deltat_the(klon, klev), d_deltaq_the(klon, klev)) !! ALLOCATE( d_s_the(klon), d_dens_the(klon)) ALLOCATE(d_deltat_ajs_cv(klon, klev), d_deltaq_ajs_cv(klon, klev)) !! End of wake variables !! ALLOCATE(bils(klon)) ALLOCATE(cdragm(klon), cdragh(klon), cldh(klon), cldl(klon)) ALLOCATE(cldm(klon), cldq(klon), cldt(klon), qsat2m(klon)) ALLOCATE(cldhjn(klon), cldljn(klon), cldmjn(klon), cldtjn(klon)) ALLOCATE(JrNt(klon)) ALLOCATE(dthmin(klon), evap(klon), fder(klon), plcl(klon), plfc(klon)) ALLOCATE(prw(klon), prlw(klon), prsw(klon), zustar(klon), zu10m(klon), zv10m(klon), rh2m(klon)) ALLOCATE(s_lcl(klon)) ALLOCATE(s_pblh(klon), s_pblt(klon), s_therm(klon)) ! !nrlmd+jyg< ALLOCATE(s_pblh_x(klon), s_pblh_w(klon)) ALLOCATE(s_lcl_x(klon), s_lcl_w(klon)) !>nrlmd+jyg ! ALLOCATE(slab_wfbils(klon), tpot(klon), tpote(klon), ue(klon)) ALLOCATE(uq(klon), ve(klon), vq(klon), zxffonte(klon)) ALLOCATE(zxfqcalving(klon), zxfluxlat(klon)) ALLOCATE(zxtsol(klon), snow_lsc(klon), zxfqfonte(klon), zxqsurf(klon)) ALLOCATE(zxrunofflic(klon)) ALLOCATE(rain_lsc(klon)) ALLOCATE(rain_num(klon)) ! ALLOCATE(sens_x(klon), sens_w(klon)) ALLOCATE(zxfluxlat_x(klon), zxfluxlat_w(klon)) !jyg< !! ALLOCATE(t_x(klon,klev), t_w(klon,klev)) !! ALLOCATE(q_x(klon,klev), q_w(klon,klev)) !>jyg ALLOCATE(dtvdf_x(klon,klev), dtvdf_w(klon,klev)) ALLOCATE(dqvdf_x(klon,klev), dqvdf_w(klon,klev)) ALLOCATE(pbl_tke_input(klon,klev+1,nbsrf)) ALLOCATE(t_therm(klon,klev), q_therm(klon,klev),u_therm(klon,klev), v_therm(klon,klev)) ALLOCATE(cdragh_x(klon), cdragh_w(klon)) ALLOCATE(cdragm_x(klon), cdragm_w(klon)) ALLOCATE(kh(klon), kh_x(klon), kh_w(klon)) ! ALLOCATE(ptconv(klon,klev)) ! ALLOCATE(wbeff(klon), convoccur(klon), zmax_th(klon)) ALLOCATE(zq2m(klon), zt2m(klon), weak_inversion(klon)) ALLOCATE(zt2m_min_mon(klon), zt2m_max_mon(klon)) ALLOCATE(t2m_min_mon(klon), t2m_max_mon(klon)) ALLOCATE(sens(klon), flwp(klon), fiwp(klon)) ALLOCATE(ale_bl_stat(klon), alp_bl_conv(klon), alp_bl_det(klon)) ALLOCATE(alp_bl_fluct_m(klon), alp_bl_fluct_tke(klon)) ALLOCATE(alp_bl_stat(klon), n2(klon), s2(klon)) ALLOCATE(proba_notrig(klon), random_notrig(klon)) ALLOCATE(dnwd(klon, klev), dnwd0(klon, klev)) ! ALLOCATE(upwd(klon, klev), omega(klon, klev), coefh(klon, klev)) ALLOCATE(upwd(klon, klev), omega(klon, klev)) ALLOCATE(epmax_diag(klon)) ! epmax_cape ALLOCATE(ep(klon,klev)) ! epmax_cape ! ALLOCATE(coefm(klon, klev), lambda_th(klon, klev), cldemi(klon, klev)) ALLOCATE(lambda_th(klon, klev), cldemi(klon, klev)) ALLOCATE(cldfra(klon, klev), cldtau(klon, klev), fiwc(klon, klev)) ALLOCATE(fl(klon, klev), re(klon, klev), flwc(klon, klev)) ALLOCATE(ref_liq(klon, klev), ref_ice(klon, klev), theta(klon, klev)) ALLOCATE(ref_liq_pi(klon, klev), ref_ice_pi(klon, klev)) ALLOCATE(zphi(klon, klev), zx_rh(klon, klev)) ALLOCATE(pmfd(klon, klev), pmfu(klon, klev)) ALLOCATE(t2m(klon, nbsrf), fluxlat(klon, nbsrf)) ALLOCATE(fsollw(klon, nbsrf)) ALLOCATE(fsolsw(klon, nbsrf), wfbils(klon, nbsrf), wfbilo(klon, nbsrf)) ALLOCATE(wfevap(klon, nbsrf), wfrain(klon,nbsrf), wfsnow(klon, nbsrf)) ALLOCATE(evap_pot(klon, nbsrf)) ALLOCATE(pmflxr(klon, klev+1), pmflxs(klon, klev+1), prfl(klon, klev+1)) ALLOCATE(psfl(klon, klev+1), fraca(klon, klev+1), Vprecip(klon, klev+1)) ALLOCATE(zw2(klon, klev+1)) ALLOCATE(fluxu(klon, klev, nbsrf), fluxv(klon, klev, nbsrf)) ALLOCATE(fluxt(klon, klev, nbsrf)) ALLOCATE(uwriteSTD(klon,nlevSTD,nfiles), vwriteSTD(klon,nlevSTD,nfiles)) ALLOCATE(wwriteSTD(klon,nlevSTD,nfiles), phiwriteSTD(klon,nlevSTD,nfiles)) ALLOCATE(qwriteSTD(klon,nlevSTD,nfiles), twriteSTD(klon,nlevSTD,nfiles)) ALLOCATE(rhwriteSTD(klon,nlevSTD,nfiles)) ! ug et d'autres encore: ALLOCATE(wdtrainA(klon,klev),wdtrainM(klon,klev)) ALLOCATE(beta_prec(klon,klev)) ALLOCATE(rneb(klon,klev),rnebjn(klon,klev)) ALLOCATE (zxsnow(klon),snowhgt(klon),qsnow(klon),to_ice(klon)) ALLOCATE (sissnow(klon),runoff(klon),albsol3_lic(klon)) #ifdef CPP_StratAer ALLOCATE (R2SO4(klon,klev)) ALLOCATE (DENSO4(klon,klev)) ALLOCATE (f_r_wet(klon,klev)) ALLOCATE (decfluxaer(klon,nbtr)) ALLOCATE (mdw(nbtr)) ALLOCATE (budg_3D_nucl(klon,klev)) ALLOCATE (budg_3D_cond_evap(klon,klev)) ALLOCATE (budg_3D_ocs_to_so2(klon,klev)) ALLOCATE (budg_3D_so2_to_h2so4(klon,klev)) ALLOCATE (budg_3D_backgr_ocs(klon,klev)) ALLOCATE (budg_3D_backgr_so2(klon,klev)) ALLOCATE (OCS_lifetime(klon,klev)) ALLOCATE (SO2_lifetime(klon,klev)) ALLOCATE (alpha_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave+nwave_lw,nbtr)) ALLOCATE (piz_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave+nwave_lw,nbtr)) ALLOCATE (cg_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave+nwave_lw,nbtr)) ALLOCATE (tau_strat_550(klon,klev)) ALLOCATE (tau_strat_1020(klon,klev)) ALLOCATE (tausum_strat(klon,3)) ALLOCATE (budg_dep_dry_h2so4(klon)) ALLOCATE (budg_dep_wet_h2so4(klon)) ALLOCATE (budg_dep_dry_part(klon)) ALLOCATE (budg_dep_wet_part(klon)) ALLOCATE (budg_emi_ocs(klon)) ALLOCATE (budg_emi_so2(klon)) ALLOCATE (budg_emi_h2so4(klon)) ALLOCATE (budg_emi_part(klon)) ALLOCATE (budg_ocs_to_so2(klon)) ALLOCATE (budg_so2_to_h2so4(klon)) ALLOCATE (budg_h2so4_to_part(klon)) ALLOCATE (budg_sed_part(klon)) ALLOCATE (surf_PM25_sulf(klon)) ALLOCATE (p_tropopause(klon)) ALLOCATE (vsed_aer(klon,klev)) #endif ALLOCATE (pr_tropopause(klon)) END SUBROUTINE phys_local_var_init !====================================================================== SUBROUTINE phys_local_var_end USE dimphy USE indice_sol_mod IMPLICIT NONE DEALLOCATE(t_seri,q_seri,ql_seri,qs_seri) DEALLOCATE(u_seri,v_seri) DEALLOCATE(l_mixmin,l_mix) DEALLOCATE(tr_seri) DEALLOCATE(d_t_dyn,d_q_dyn) DEALLOCATE(d_ql_dyn,d_qs_dyn) DEALLOCATE(d_q_dyn2d,d_ql_dyn2d,d_qs_dyn2d) DEALLOCATE(d_u_dyn,d_v_dyn) DEALLOCATE(d_tr_dyn) !RomP DEALLOCATE(d_t_con,d_q_con) DEALLOCATE(d_u_con,d_v_con) DEALLOCATE(d_t_wake,d_q_wake) DEALLOCATE(d_t_lsc,d_q_lsc) DEALLOCATE(d_t_lwr,d_t_lw0) DEALLOCATE(d_t_swr,d_t_sw0) DEALLOCATE(d_ql_lsc,d_qi_lsc) DEALLOCATE(d_t_ajsb,d_q_ajsb) DEALLOCATE(d_t_ajs,d_q_ajs) !nrlmd< DEALLOCATE(d_t_ajs_w,d_q_ajs_w) DEALLOCATE(d_t_ajs_x,d_q_ajs_x) !>nrlmd DEALLOCATE(d_u_ajs,d_v_ajs) DEALLOCATE(d_t_eva,d_q_eva) DEALLOCATE(d_ql_eva,d_qi_eva) DEALLOCATE(d_t_lscst,d_q_lscst) DEALLOCATE(d_t_lscth,d_q_lscth) DEALLOCATE(plul_st,plul_th) DEALLOCATE(d_t_vdf,d_q_vdf,d_t_diss) !nrlmd+jyg< DEALLOCATE(d_t_vdf_w,d_q_vdf_w) DEALLOCATE(d_t_vdf_x,d_q_vdf_x) !>nrlmd+jyg DEALLOCATE(d_u_vdf,d_v_vdf) DEALLOCATE(d_t_oli,d_t_oro) DEALLOCATE(d_u_oli,d_v_oli) DEALLOCATE(d_u_oro,d_v_oro) DEALLOCATE(d_t_lif,d_t_ec) DEALLOCATE(d_u_lif,d_v_lif) DEALLOCATE(d_ts, d_tr) DEALLOCATE(topswad_aero,solswad_aero) DEALLOCATE(topswai_aero,solswai_aero) DEALLOCATE(topswad0_aero,solswad0_aero) DEALLOCATE(toplwad_aero,sollwad_aero) DEALLOCATE(toplwai_aero,sollwai_aero) DEALLOCATE(toplwad0_aero,sollwad0_aero) DEALLOCATE(topsw_aero,solsw_aero) DEALLOCATE(topsw0_aero,solsw0_aero) DEALLOCATE(topswcf_aero,solswcf_aero) DEALLOCATE(stratomask) DEALLOCATE(tausum_aero) DEALLOCATE(drytausum_aero) DEALLOCATE(tau3d_aero) DEALLOCATE(scdnc) DEALLOCATE(cldncl) DEALLOCATE(reffclwtop) DEALLOCATE(lcc) DEALLOCATE(reffclws) DEALLOCATE(reffclwc) DEALLOCATE(cldnvi) DEALLOCATE(lcc3d) DEALLOCATE(lcc3dcon) DEALLOCATE(lcc3dstra) DEALLOCATE(od443aer) DEALLOCATE(od550aer) DEALLOCATE(od865aer) DEALLOCATE(dryod550aer) DEALLOCATE(abs550aer) DEALLOCATE(ec550aer) DEALLOCATE(od550lt1aer) DEALLOCATE(sconcso4) DEALLOCATE(sconcno3) DEALLOCATE(sconcoa) DEALLOCATE(sconcbc) DEALLOCATE(sconcss) DEALLOCATE(sconcdust) DEALLOCATE(concso4) DEALLOCATE(concno3) DEALLOCATE(concoa) DEALLOCATE(concbc) DEALLOCATE(concss) DEALLOCATE(concdust) DEALLOCATE(loadso4) DEALLOCATE(loadoa) DEALLOCATE(loadbc) DEALLOCATE(loadss) DEALLOCATE(loaddust) DEALLOCATE(loadno3) DEALLOCATE(load_tmp1) DEALLOCATE(load_tmp2) DEALLOCATE(load_tmp3) DEALLOCATE(du_gwd_hines,dv_gwd_hines,d_t_hin) DEALLOCATE(d_q_ch4) DEALLOCATE(dv_gwd_rando,dv_gwd_front) DEALLOCATE(east_gwstress,west_gwstress) !IM ajout variables CFMIP2/CMIP5 DEALLOCATE(topswad_aerop, solswad_aerop) DEALLOCATE(topswai_aerop, solswai_aerop) DEALLOCATE(topswad0_aerop, solswad0_aerop) DEALLOCATE(topsw_aerop, topsw0_aerop) DEALLOCATE(solsw_aerop, solsw0_aerop) DEALLOCATE(topswcf_aerop, solswcf_aerop) !CK LW diagnostics DEALLOCATE(toplwad_aerop, sollwad_aerop) DEALLOCATE(toplwai_aerop, sollwai_aerop) DEALLOCATE(toplwad0_aerop, sollwad0_aerop) ! FH Ajout de celles nécessaires au phys_output_write_mod DEALLOCATE(tal1, pal1, pab1, pab2) DEALLOCATE(ptstar, pt0, slp) ! DEALLOCATE(ale_wake, alp_wake) DEALLOCATE(wake_h,wake_k) DEALLOCATE(wake_omg) DEALLOCATE(d_deltat_wk, d_deltaq_wk) DEALLOCATE(d_s_wk, d_dens_wk) DEALLOCATE(d_deltat_wk_gw, d_deltaq_wk_gw) DEALLOCATE(d_deltat_vdf, d_deltaq_vdf) !! DEALLOCATE( d_s_vdf, d_dens_vdf) DEALLOCATE(d_deltat_the, d_deltaq_the) !! DEALLOCATE( d_s_the, d_dens_the) DEALLOCATE(d_deltat_ajs_cv, d_deltaq_ajs_cv) ! DEALLOCATE(bils) DEALLOCATE(cdragm, cdragh, cldh, cldl) DEALLOCATE(cldm, cldq, cldt, qsat2m) DEALLOCATE(cldljn, cldmjn, cldhjn, cldtjn, JrNt) DEALLOCATE(dthmin, evap, fder, plcl, plfc) DEALLOCATE(prw, prlw, prsw, zustar, zu10m, zv10m, rh2m, s_lcl) DEALLOCATE(s_pblh, s_pblt, s_therm) ! !nrlmd+jyg< DEALLOCATE(s_pblh_x, s_pblh_w) DEALLOCATE(s_lcl_x, s_lcl_w) !>nrlmd+jyg ! DEALLOCATE(slab_wfbils, tpot, tpote, ue) DEALLOCATE(uq, ve, vq, zxffonte) DEALLOCATE(zxfqcalving, zxfluxlat) DEALLOCATE(zxrunofflic) DEALLOCATE(zxtsol, snow_lsc, zxfqfonte, zxqsurf) DEALLOCATE(rain_lsc) DEALLOCATE(rain_num) ! DEALLOCATE(sens_x, sens_w) DEALLOCATE(zxfluxlat_x, zxfluxlat_w) !jyg< !! DEALLOCATE(t_x, t_w) !! DEALLOCATE(q_x, q_w) !>jyg DEALLOCATE(dtvdf_x, dtvdf_w) DEALLOCATE(dqvdf_x, dqvdf_w) DEALLOCATE(pbl_tke_input) DEALLOCATE(t_therm, q_therm, u_therm, v_therm) DEALLOCATE(cdragh_x, cdragh_w) DEALLOCATE(cdragm_x, cdragm_w) DEALLOCATE(kh, kh_x, kh_w) ! DEALLOCATE(ptconv) ! DEALLOCATE(wbeff, convoccur, zmax_th) DEALLOCATE(zq2m, zt2m, weak_inversion) DEALLOCATE(zt2m_min_mon, zt2m_max_mon) DEALLOCATE(t2m_min_mon, t2m_max_mon) DEALLOCATE(sens, flwp, fiwp) DEALLOCATE(ale_bl_stat,alp_bl_conv,alp_bl_det) DEALLOCATE(alp_bl_fluct_m,alp_bl_fluct_tke) DEALLOCATE(alp_bl_stat, n2, s2) DEALLOCATE(proba_notrig, random_notrig) DEALLOCATE(dnwd, dnwd0) ! DEALLOCATE(upwd, omega, coefh) DEALLOCATE(upwd, omega) DEALLOCATE(epmax_diag) DEALLOCATE(ep) ! DEALLOCATE(coefm, lambda_th, cldemi) DEALLOCATE(lambda_th, cldemi) DEALLOCATE(cldfra, cldtau, fiwc) DEALLOCATE(fl, re, flwc) DEALLOCATE(ref_liq, ref_ice, theta) DEALLOCATE(ref_liq_pi, ref_ice_pi) DEALLOCATE(zphi, zx_rh) DEALLOCATE(pmfd, pmfu) DEALLOCATE(t2m, fluxlat) DEALLOCATE(fsollw, evap_pot) DEALLOCATE(fsolsw, wfbils, wfbilo) DEALLOCATE(wfevap,wfrain,wfsnow) DEALLOCATE(pmflxr, pmflxs, prfl) DEALLOCATE(psfl, fraca, Vprecip) DEALLOCATE(zw2) DEALLOCATE(fluxu, fluxv) DEALLOCATE(fluxt) DEALLOCATE(uwriteSTD, vwriteSTD) DEALLOCATE(wwriteSTD, phiwriteSTD) DEALLOCATE(qwriteSTD, twriteSTD, rhwriteSTD) ! ug et d'autres encore: DEALLOCATE(wdtrainA, wdtrainM) DEALLOCATE(beta_prec) DEALLOCATE(rneb) DEALLOCATE (zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic) #ifdef CPP_StratAer ! variables for strat. aerosol CK DEALLOCATE (R2SO4) DEALLOCATE (DENSO4) DEALLOCATE (f_r_wet) DEALLOCATE (decfluxaer) DEALLOCATE (mdw) DEALLOCATE (SO2_lifetime) DEALLOCATE (OCS_lifetime) DEALLOCATE (alpha_bin) DEALLOCATE (piz_bin) DEALLOCATE (cg_bin) DEALLOCATE (tau_strat_550) DEALLOCATE (tau_strat_1020) DEALLOCATE (tausum_strat) DEALLOCATE (surf_PM25_sulf) DEALLOCATE (p_tropopause) DEALLOCATE (vsed_aer) DEALLOCATE (budg_3D_ocs_to_so2) DEALLOCATE (budg_3D_so2_to_h2so4) DEALLOCATE (budg_3D_backgr_ocs) DEALLOCATE (budg_3D_backgr_so2) DEALLOCATE (budg_3D_nucl) DEALLOCATE (budg_3D_cond_evap) DEALLOCATE (budg_dep_dry_h2so4) DEALLOCATE (budg_dep_wet_h2so4) DEALLOCATE (budg_dep_dry_part) DEALLOCATE (budg_dep_wet_part) DEALLOCATE (budg_emi_ocs) DEALLOCATE (budg_emi_so2) DEALLOCATE (budg_emi_h2so4) DEALLOCATE (budg_emi_part) DEALLOCATE (budg_ocs_to_so2) DEALLOCATE (budg_so2_to_h2so4) DEALLOCATE (budg_h2so4_to_part) DEALLOCATE (budg_sed_part) #endif DEALLOCATE (pr_tropopause) END SUBROUTINE phys_local_var_end END MODULE phys_local_var_mod