!

! $Id: phys_output_write_mod.F90 4575 2023-06-14 15:06:15Z fairhead $

!

MODULE phys_output_write_mod

  USE phytrac_mod, ONLY : d_tr_cl, d_tr_th, d_tr_cv, d_tr_lessi_impa, &

       d_tr_lessi_nucl, d_tr_insc, d_tr_bcscav, d_tr_evapls, d_tr_ls,  &

       d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav, flux_tr_dry

  ! Author: Abderrahmane IDELKADI (original include file)

  ! Author: Laurent FAIRHEAD (transformation to module/subroutine)

  ! Author: Ulysse GERARD (effective implementation)

CONTAINS

  ! ug Routine pour définir (lors du premier passage) ET sortir les variables

       ok_ade, ok_aie, ok_volcan, ivap, iliq, isol, ibs, ok_sync, &

       ptconv, read_climoz, clevSTD, ptconvth, &

    ! This subroutine does the actual writing of diagnostics that were

    ! defined and initialised in phys_output_mod.F90

    USE dimphy, ONLY: klon, klev, klevp1

    USE infotrac_phy, ONLY: nbtr, nqtot, nqo, type_trac, tracers, niso, ntiso

    USE strings_mod,  ONLY: maxlen

    USE mod_phys_lmdz_para, ONLY: is_north_pole_phy,is_south_pole_phy

    USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat

    USE time_phylmdz_mod, ONLY: day_step_phy, start_time, itau_phy

    USE vertical_layers_mod, ONLY : ap, bp, aps, bps

    USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, &

         o_longitude, o_latitude, &

         o_Ahyb, o_Bhyb,o_Ahyb_bounds, o_Bhyb_bounds, &

         o_Ahyb_mid, o_Bhyb_mid,o_Ahyb_mid_bounds, o_Bhyb_mid_bounds, &

         is_ave, is_sic, o_contfracATM, o_contfracOR, &

         o_aireTER, o_flat, o_slp, o_ptstar, o_pt0, o_tsol, &

         o_t2m, o_t2m_min, o_t2m_max, &

         o_t2m_min_mon, o_t2m_max_mon, &

         o_nt2mout, o_nt2moutfg, &

         o_nq2mout, o_nq2moutfg, &

         o_nu2mout, o_nu2moutfg, &

         o_q2m, o_ustar, o_u10m, o_v10m, &

         o_wind10m, o_wind10max, o_wind100m, o_gusts, o_sicf, &

         o_loadfactor_wind_onshore, o_loadfactor_wind_offshore, &

         o_psol, o_mass, o_qsurf, o_qsol, &

         o_precip, o_rain_fall, o_rain_con, o_ndayrain, o_plul, o_pluc, o_plun, &

         o_snow, o_msnow, o_fsnow, o_evap, o_snowerosion, o_ustart_lic, o_rhosnow_lic, o_bsfall, &

         o_ep,o_epmax_diag, & ! epmax_cape

         o_tops, o_tops0, o_topl, o_topl0, &

         o_SWupTOA, o_SWupTOAclr, o_SWupTOAcleanclr, o_SWdnTOA, o_fdiffSWdnSFC, &

         o_SWdnTOAclr, o_nettop, o_SWup200, &

         o_SWup200clr, o_SWdn200, o_SWdn200clr, &

         o_LWup200, o_LWup200clr, o_LWdn200, &

         o_LWdn200clr, o_sols, o_sols0, &

         o_soll, o_radsol, o_soll0, o_SWupSFC, &

         o_SWupSFCclr, o_SWupSFCcleanclr, o_SWdnSFC, o_SWdnSFCclr, o_SWdnSFCcleanclr, &

         o_LWupSFC, o_LWdnSFC, o_LWupSFCclr, &

         o_LWdnSFCclr, o_LWupTOAcleanclr, o_LWdnSFCcleanclr, o_bils, o_bils_diss, &

         o_bils_ec,o_bils_ech, o_bils_tke, o_bils_kinetic, &

         o_bils_latent, o_bils_enthalp, o_sens, &

         o_fder, o_ffonte, o_fqcalving, o_fqfonte, o_mrroli, o_runofflic, &

         o_taux, o_tauy, o_snowsrf, o_qsnow, &

         o_snowhgt, o_toice, o_sissnow, o_runoff, &

         o_albslw3, o_pourc_srf, o_fract_srf, &

         o_taux_srf, o_tauy_srf, o_tsol_srf, &

         o_evappot_srf, o_ustar_srf, o_u10m_srf, &

         o_v10m_srf, o_t2m_srf, o_evap_srf, &

         o_sens_srf, o_lat_srf, o_flw_srf, &

         o_fsw_srf, o_wbils_srf, o_wbilo_srf, &

         o_wevap_srf, o_wrain_srf, o_wsnow_srf, &

         o_tke_srf, o_tke_max_srf,o_dltpbltke_srf, o_wstar, &

         o_l_mixmin,o_l_mix, &

         o_cdrm, o_cdrh, o_cldl, o_cldm, o_cldh, &

         o_cldt, o_JrNt, o_cldljn, o_cldmjn, &

         o_cldhjn, o_cldtjn, o_cldq, o_lwp, o_iwp, &

         o_ue, o_ve, o_uq, o_vq, o_cape, o_pbase, &

         o_uwat, o_vwat, &

         o_ptop, o_fbase, o_plcl, o_plfc, &

         o_wbeff, o_convoccur, o_cape_max, o_upwd, o_ep,o_epmax_diag, &

         o_Mipsh, o_Ma, &

         o_dnwd, o_dnwd0, o_ftime_deepcv, o_ftime_con, o_mc, &

         o_prw, o_prlw, o_prsw, o_prbsw, o_s_pblh, o_s_pblt, o_s_lcl, &

         o_s_therm, o_uSTDlevs, o_vSTDlevs, &

         o_wSTDlevs, o_zSTDlevs, o_qSTDlevs, &

         o_tSTDlevs, epsfra, o_t_oce_sic, &

         o_ale_bl, o_alp_bl, o_ale_wk, o_alp_wk, &

         o_dtvdf_x    , o_dtvdf_w    , o_dqvdf_x    , o_dqvdf_w    , &

         o_sens_x     , o_sens_w     , o_flat_x     , o_flat_w     , &

         o_delta_tsurf, o_delta_tsurf_srf, &

         o_cdragh_x   , o_cdragh_w   , o_cdragm_x   , o_cdragm_w   , &

         o_kh         , o_kh_x       , o_kh_w       , &

         o_ale, o_alp, o_cin, o_WAPE, o_wake_h, o_cv_gen, o_wake_dens, &

         o_wake_s, o_wake_deltat, o_wake_deltaq, &

         o_wake_omg, o_dtwak, o_dqwak, o_dqwak2d, o_Vprecip, &

         o_qtaa, o_Clwaa, &

         o_ftd, o_fqd, o_wdtrainA, o_wdtrainS, o_wdtrainM, &

         o_n2, o_s2, o_proba_notrig, &

         o_random_notrig, o_ale_bl_stat, &

         o_ale_bl_trig, o_alp_bl_det, &

         o_alp_bl_fluct_m, o_alp_bl_fluct_tke, &

         o_alp_bl_conv, o_alp_bl_stat, &

         o_slab_qflux, o_tslab, o_slab_bils, &

         o_slab_bilg, o_slab_sic, o_slab_tice, &

         o_slab_hdiff, o_slab_ekman, o_slab_gm,  &

         o_weakinv, o_dthmin, o_cldtau, &

         o_cldemi, o_pr_con_l, o_pr_con_i, &

         o_pr_lsc_l, o_pr_lsc_i, o_pr_bs, o_re, o_fl, &

         o_rh2m, &

         o_qsat2m, o_tpot, o_tpote, o_SWnetOR, &

         o_LWdownOR, o_snowl, &

         o_solldown, o_dtsvdfo, o_dtsvdft, &

         o_dtsvdfg, o_dtsvdfi, o_z0m, o_z0h,  o_od443aer, o_od550aer, &

         o_dryod550aer, o_od865aer, o_abs550aer, o_od550lt1aer, &

         o_sconcso4, o_sconcno3, o_sconcoa, o_sconcbc, &

         o_sconcss, o_sconcdust, o_concso4, o_concno3, &

         o_concoa, o_concbc, o_concss, o_concdust, &

         o_loadso4, o_loadoa, o_loadbc, o_loadss, &

         o_loaddust, o_loadno3, o_tausumaero, &

         o_drytausumaero, o_tausumaero_lw, &

         o_topswad, o_topswad0, o_solswad, o_solswad0, &

         o_toplwad, o_toplwad0, o_sollwad, o_sollwad0, &

         o_swtoaas_nat, o_swsrfas_nat, &

         o_swtoacs_nat, o_swtoaas_ant, &

         o_swsrfas_ant, o_swtoacs_ant, &

         o_swsrfcs_ant, o_swtoacf_nat, &

         o_swsrfcf_nat, o_swtoacf_ant, &

         o_swsrfcs_nat, o_swsrfcf_ant, &

         o_swtoacf_zero, o_swsrfcf_zero, &

         o_topswai, o_solswai, o_toplwai, o_sollwai, o_scdnc, &

         o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,&

         o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &

         o_lcc3dstra, o_icc3dcon, o_icc3dstra, &

         o_cldicemxrat, o_cldwatmxrat, o_reffclwtop, o_ec550aer, &

         o_lwcon, o_iwcon, o_temp, o_theta, &

         o_ovapinit, o_ovap, o_oliq, o_ocond, o_geop,o_qbs, &

         o_vitu, o_vitv, o_vitw, o_pres, o_paprs, &

         o_zfull, o_zhalf, o_rneb, o_rnebjn, o_rnebcon, &

         o_rnebls, o_rneblsvol, o_rhum, o_rhl, o_rhi, o_ozone, o_ozone_light, &

         o_pfraclr, o_pfracld, &

         o_duphy, o_dtphy, o_dqphy, o_dqphy2d, o_dqlphy, o_dqlphy2d, &

         o_dqsphy, o_dqsphy2d, o_dqbsphy, o_dqbsphy2d, o_albe_srf, o_z0m_srf, o_z0h_srf, &

         o_ages_srf, o_snow_srf, o_alb1, o_alb2, o_tke, o_tke_dissip, &

         o_tke_max, o_kz, o_kz_max, o_clwcon, &

         o_dtdyn, o_dqdyn, o_dqdyn2d, o_dqldyn, o_dqldyn2d, &

         o_dqsdyn, o_dqsdyn2d, o_dqbsdyn, o_dqbsdyn2d, o_dudyn, o_dvdyn, &

         o_dtcon, o_tntc, o_ducon, o_dvcon, &

         o_dqcon, o_dqcon2d, o_tnhusc, o_tnhusc, o_dtlsc, &

         o_dtlschr, o_dqlsc, o_dqlsc2d, o_beta_prec, &

         o_dtlscth, o_dtlscst, o_dqlscth, o_dqlscth2d, &

         o_dqlscst, o_dqlscst2d, o_plulth, o_plulst, &

         o_ptconvth, o_lmaxth, o_dtvdf, &

         o_dtdis, o_dqvdf, o_dqvdf2d, o_dteva, o_dqeva, o_dqeva2d, &

         o_dqbsvdf, o_dtbs, o_dqbs, o_dqbsbs, &

         o_ptconv, o_ratqs, o_dtthe, &

         o_duthe, o_dvthe, o_ftime_th, &

         o_f_th, o_e_th, o_w_th, o_q_th, &

         o_a_th, o_cloudth_sth, o_cloudth_senv, &

         o_cloudth_sigmath, o_cloudth_sigmaenv, &

         o_d_th, o_f0_th, o_zmax_th, &

         o_dqthe, o_dqthe2d, o_dtajs, o_dqajs, o_dqajs2d, o_dtswr, &

         o_dtsw0, o_dtlwr, o_dtlw0, o_dtec, &

         o_duvdf, o_dvvdf, o_duoro, o_dvoro, &

         o_dtoro, o_dulif, o_dvlif, o_dtlif, &

         o_du_gwd_hines, o_dv_gwd_hines, o_dthin, o_dqch4, o_rsu, &

         o_du_gwd_front, o_dv_gwd_front, &

         o_east_gwstress, o_west_gwstress, &

         o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, o_rsucsaf, o_rsdcsaf, &

         o_rlucs, o_rldcs, o_tnt, o_tntr, &

         o_tntscpbl, o_tnhus, o_tnhusscpbl, &

         o_evu, o_h2o, o_mcd, o_dmc, o_ref_liq, &

         o_ref_ice, o_rsut4co2, o_rlut4co2, &

         o_rsutcs4co2, o_rlutcs4co2, o_rsu4co2, &

         o_rlu4co2, o_rsucs4co2, o_rlucs4co2, &

         o_rsd4co2, o_rld4co2, o_rsdcs4co2, &

         o_rldcs4co2, o_tnondef, o_ta, o_zg, &

         o_hus, o_hur, o_ua, o_va, o_wap, &

         o_psbg, o_tro3, o_tro3_daylight, &

         o_uxv, o_vxq, o_vxT, o_wxq, o_vxphi, &

         o_wxT, o_uxu, o_vxv, o_TxT, o_trac, &

         o_dtr_vdf, o_dtr_the, o_dtr_con, &

         o_dtr_lessi_impa, o_dtr_lessi_nucl, &

         o_dtr_insc, o_dtr_bcscav, o_dtr_evapls, &

         o_dtr_ls, o_dtr_trsp, o_dtr_sscav, o_dtr_dry, &

         o_dtr_sat, o_dtr_uscav, o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &

         o_ustr_gwd_hines,o_vstr_gwd_hines,o_ustr_gwd_rando,o_vstr_gwd_rando, &

         o_ustr_gwd_front,o_vstr_gwd_front, &

         o_sens_prec_liq_oce, o_sens_prec_liq_sic, &

         o_sens_prec_sol_oce, o_sens_prec_sol_sic, &

         o_lat_prec_liq_oce, o_lat_prec_liq_sic, &

         o_lat_prec_sol_oce, o_lat_prec_sol_sic, &

         o_sza, &

! Marine

         o_map_prop_hc, o_map_prop_hist, o_map_emis_hc, o_map_iwp_hc, &

         o_map_deltaz_hc, o_map_pcld_hc, o_map_tcld_hc, &

         o_map_emis_hist, o_map_iwp_hist, o_map_deltaz_hist, &

         o_map_rad_hist, &

         o_map_emis_Cb, o_map_pcld_Cb, o_map_tcld_Cb, &

         o_map_emis_ThCi, o_map_pcld_ThCi, o_map_tcld_ThCi, &

         o_map_emis_Anv, o_map_pcld_Anv, o_map_tcld_Anv, &

         o_map_ntot, o_map_hc,o_map_hist,o_map_Cb,o_map_ThCi,o_map_Anv, &

#ifdef ISO

! Isotopes

         o_xtprecip,o_xtplul,o_xtpluc,o_xtovap,o_xtoliq,o_xtcond, &

         o_xtevap,o_dxtdyn,o_dxtldyn,o_dxtcon,o_dxtlsc,o_dxteva, &

         o_dxtajs,o_dxtvdf,o_dxtthe, o_dxtch4, &

         o_dxtprod_nucl,o_dxtcosmo,o_dxtdecroiss, &

         o_xtevap_srf, &

#endif

! Tropopause

         o_alt_tropo, &

         o_p_tropopause, o_z_tropopause, o_t_tropopause,  &

         o_col_O3_strato, o_col_O3_tropo,                 &

!--aviation & supersaturation

         o_oclr, o_ocld, o_oss, o_ovc, o_rnebss, o_rnebclr, o_rnebseri, o_gammass, &

         o_N1_ss, o_N2_ss, o_qsatl, o_qsats, &

         o_drnebsub, o_drnebcon, o_drnebtur, o_drnebavi, o_flight_m, o_flight_h2o, &

         o_Tcontr, o_qcontr, o_qcontr2, o_fcontrN, o_fcontrP, &

!--interactive CO2

         o_flx_co2_ocean, o_flx_co2_ocean_cor, &

         o_flx_co2_land, o_flx_co2_land_cor, &

         o_flx_co2_ff, o_flx_co2_bb, &

         o_delta_sst, o_delta_sal, o_ds_ns, o_dt_ns, o_dter, o_dser, o_tkt, &

         o_tks, o_taur, o_sss, &

!FC

         o_zxfluxt,o_zxfluxq

#ifdef CPP_StratAer

    USE phys_output_ctrlout_mod, ONLY:  &

         o_budg_3D_nucl, o_budg_3D_cond_evap, o_budg_3D_ocs_to_so2, o_budg_3D_so2_to_h2so4, &

         o_budg_sed_part, o_R2SO4, o_OCS_lifetime, o_SO2_lifetime, &

         o_budg_3D_backgr_ocs, o_budg_3D_backgr_so2, &

         o_budg_dep_dry_ocs, o_budg_dep_wet_ocs, &

         o_budg_dep_dry_so2, o_budg_dep_wet_so2, &

         o_budg_dep_dry_h2so4, o_budg_dep_wet_h2so4, &

         o_budg_dep_dry_part, o_budg_dep_wet_part, &

         o_budg_emi_ocs, o_budg_emi_so2, o_budg_emi_h2so4, o_budg_emi_part, &

         o_budg_ocs_to_so2, o_budg_so2_to_h2so4, o_budg_h2so4_to_part, &

         o_surf_PM25_sulf, o_ext_strat_550, o_tau_strat_550, &

         o_vsed_aer, o_tau_strat_1020, o_ext_strat_1020, o_f_r_wet

#endif

    USE ice_sursat_mod, ONLY: flight_m, flight_h2o

    USE phys_output_ctrlout_mod, ONLY: o_heat_volc, o_cool_volc !NL

    USE phys_state_var_mod, ONLY: heat_volc, cool_volc !NL

    USE phys_state_var_mod, ONLY: pctsrf, rain_fall, snow_fall, bs_fall,&

         qsol, z0m, z0h, fevap, agesno, &

         nday_rain, ndayrain_mth, rain_con, snow_con, &

         topsw, toplw, toplw0, swup, swdn, solswfdiff, &

         topsw0, swupc0, swdnc0, swup0, swdn0, SWup200, SWup200clr, &

         SWdn200, SWdn200clr, LWup200, LWup200clr, &

         LWdn200, LWdn200clr, solsw, solsw0, sollw, &

         radsol, swradcorr, sollw0, sollwdown, sollw, gustiness, &

         sollwdownclr, lwdnc0, lwdn0, ftsol, ustar, u10m, &

         v10m, pbl_tke, wake_delta_pbl_TKE, &

         delta_tsurf, &

         wstar, cape, ema_pcb, ema_pct, &

         ema_cbmf, Mipsh, Ma, fm_therm, ale_bl, alp_bl, ale, &

         alp, cin, wake_pe, wake_dens, cv_gen, wake_s, wake_deltat, &

         wake_deltaq, ftd, fqd, ale_bl_trig, albsol1, &

         ale_wake, ale_bl_stat, &

         rnebcon, wo, falb1, albsol2, coefh, clwcon0, &

         ratqs, entr_therm, zqasc, detr_therm, f0, &

         lwup, lwdn, lwupc0, lwup0, coefm, &

         swupp, lwupp, swupc0p, swup0p, lwupc0p, lwup0p, swdnp, lwdnp, &

         swdnc0p, swdn0p, lwdnc0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &

         vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &

         vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &

         T2sumSTD, nlevSTD, du_gwd_rando, du_gwd_front, &

         ulevSTD, vlevSTD, wlevSTD, philevSTD, qlevSTD, tlevSTD, &

         rhlevSTD, O3STD, O3daySTD, uvSTD, vqSTD, vTSTD, wqSTD, vphiSTD, &

         wTSTD, u2STD, v2STD, T2STD, missing_val_nf90, delta_sal, ds_ns, &

#ifdef ISO

        xtrain_con, xtsnow_con, xtrain_fall, xtsnow_fall, fxtevap, &

#endif

         dt_ns, delta_sst, dter, dser

    USE phys_local_var_mod, ONLY: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, &

         zn2mout, t2m_min_mon, t2m_max_mon, evap, &

         snowerosion, zxustartlic, zxrhoslic, &

         l_mixmin,l_mix, tke_dissip, &

         zu10m, zv10m, zq2m, zustar, zxqsurf, &

         rain_lsc, rain_num, snow_lsc, bils, sens, fder, &

         zxffonte, zxfqcalving, zxfqfonte, zxrunofflic, fluxu, &

         fluxv, zxsnow, qsnow, snowhgt, to_ice, &

         sissnow, runoff, albsol3_lic, evap_pot, &

         t2m, fluxt, fluxlat, fsollw, fsolsw, &

         wfbils, wfbilo, wfevap, wfrain, wfsnow, &

         cdragm, cdragh, cldl, cldm, &

         cldh, cldt, JrNt,   & ! only output names: cldljn,cldmjn,cldhjn,cldtjn

         cldq, flwp, fiwp, ue, ve, uq, vq, &

         uwat, vwat, &

         rneb_seri, d_rneb_dyn, &

         plcl, plfc, wbeff, convoccur, upwd, dnwd, dnwd0, prw, prlw, prsw,prbsw, &

         s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &

         vwriteSTD, wwriteSTD, phiwriteSTD, qwriteSTD, &

         twriteSTD, alp_wake, &

!!         dtvdf_x    ,dtvdf_w    ,dqvdf_x    ,dqvdf_w    , &

         d_t_vdf_x    ,d_t_vdf_w    ,d_q_vdf_x    ,d_q_vdf_w    , &

         sens_x     ,sens_w     ,zxfluxlat_x,zxfluxlat_w, &

         cdragh_x   ,cdragh_w   ,cdragm_x   ,cdragm_w   , &

         kh         ,kh_x       ,kh_w       , &

         wake_h, &

         wake_omg, d_t_wake, d_q_wake, Vprecip, qtaa, Clw, &

         wdtrainA, wdtrainS, wdtrainM, n2, s2, proba_notrig, &

         random_notrig, &

         qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &

         N1_ss, N2_ss, zqsatl, zqsats, &

         Tcontr, qcontr, qcontr2, fcontrN, fcontrP, &

         drneb_sub, drneb_con, drneb_tur, drneb_avi, &

         alp_bl_det, alp_bl_fluct_m, alp_bl_conv, &

         alp_bl_stat, alp_bl_fluct_tke, slab_wfbils, &

         weak_inversion, dthmin, cldtau, cldemi, &

         pmflxr, pmflxs, prfl, psfl,bsfl, re, fl, rh2m, &

         qsat2m, tpote, tpot, d_ts, od443aer, od550aer, dryod550aer, &

         od865aer, abs550aer, od550lt1aer, sconcso4, sconcno3, &

         sconcoa, sconcbc, sconcss, sconcdust, concso4, concno3, &

         concoa, concbc, concss, concdust, loadso4, &

         loadoa, loadbc, loadss, loaddust, loadno3, tausum_aero, drytausum_aero, &

         topswad_aero, topswad0_aero, solswad_aero, &

         solswad0_aero, topsw_aero, solsw_aero, &

         topsw0_aero, solsw0_aero, topswcf_aero, &

         solswcf_aero, topswai_aero, solswai_aero, &

         toplwad_aero, toplwad0_aero, sollwad_aero, &

         sollwad0_aero, toplwai_aero, sollwai_aero, &

         scdnc, cldncl, reffclws, reffclwc, cldnvi, stratomask,&

         lcc, lcc3d, lcc3dcon, lcc3dstra, &

         icc3dcon, icc3dstra, zfice, reffclwtop, &

         ec550aer, flwc, fiwc, t_seri, theta, q_seri, &

         ql_seri, qs_seri, qbs_seri, tr_seri, qbs_seri,&

         zphi, u_seri, v_seri, omega, cldfra, &

         rneb, rnebjn, rneblsvol, zx_rh, zx_rhl, zx_rhi, &

         pfraclr, pfracld, d_t_dyn,  &

         d_q_dyn,  d_ql_dyn, d_qs_dyn, d_qbs_dyn,  &

         d_q_dyn2d,  d_ql_dyn2d, d_qs_dyn2d, d_qbs_dyn2d, &

         d_u_dyn, d_v_dyn, d_t_con, d_t_ajsb, d_t_ajs, &

         d_u_ajs, d_v_ajs, &

         d_u_con, d_v_con, d_q_con, d_q_ajs, d_t_lsc, &

         d_t_lwr,d_t_lw0,d_t_swr,d_t_sw0, &

         d_t_eva, d_q_lsc, beta_prec, d_t_lscth, &

         d_t_lscst, d_q_lscth, d_q_lscst, plul_th, &

         plul_st, d_t_vdf, d_t_diss, d_q_vdf, d_q_eva, &

         d_t_bs, d_q_bs, d_qbs_bs, d_qbs_vdf, &

         zw2, fraca, zmax_th, d_q_ajsb, d_t_ec, d_u_vdf, &

         d_v_vdf, d_u_oro, d_v_oro, d_t_oro, d_u_lif, &

         d_v_lif, d_t_lif, du_gwd_hines, dv_gwd_hines, d_t_hin, &

         dv_gwd_rando, dv_gwd_front, &

         east_gwstress, west_gwstress, &

         d_q_ch4, pmfd, pmfu, ref_liq, ref_ice, rhwriteSTD, &

#ifdef ISO

        xtrain_lsc, xtsnow_lsc, xt_seri, xtl_seri,xts_seri,xtevap, &

        d_xt_dyn,d_xtl_dyn,d_xt_con,d_xt_vdf,d_xt_ajsb, &

        d_xt_lsc,d_xt_eva,d_xt_ch4, &

        d_xt_ajs, d_xt_ajsb, &

        d_xt_prod_nucl,d_xt_cosmo,d_xt_decroiss, &

#endif

         ep, epmax_diag, &  ! epmax_cape

         p_tropopause, t_tropopause, z_tropopause, &

         zxfluxt,zxfluxq

!FC

#ifdef CPP_StratAer

    USE phys_local_var_mod, ONLY:  &

         budg_3D_nucl, budg_3D_cond_evap, budg_3D_ocs_to_so2, budg_3D_so2_to_h2so4, &

         budg_sed_part, R2SO4, OCS_lifetime, SO2_lifetime, &

         budg_3D_backgr_ocs, budg_3D_backgr_so2, &

         budg_dep_dry_ocs, budg_dep_wet_ocs, &

         budg_dep_dry_so2, budg_dep_wet_so2, &

         budg_dep_dry_h2so4, budg_dep_wet_h2so4, &

         budg_dep_dry_part, budg_dep_wet_part, &

         budg_emi_ocs, budg_emi_so2, budg_emi_h2so4, budg_emi_part, &

         budg_ocs_to_so2, budg_so2_to_h2so4, budg_h2so4_to_part, &

         surf_PM25_sulf, tau_strat_550, tausum_strat, &

         vsed_aer, tau_strat_1020, f_r_wet

#endif

    USE carbon_cycle_mod, ONLY: fco2_ff, fco2_bb, fco2_land, fco2_ocean

    USE carbon_cycle_mod, ONLY: fco2_ocean_cor, fco2_land_cor

    USE phys_output_var_mod, ONLY: vars_defined, snow_o, zfra_o, bils_diss, &

         bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &

         itau_con, nfiles, clef_files, nid_files, dryaod_diag, &

         zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &

         zustr_gwd_front, zvstr_gwd_front, sza_o,    &

         sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &

         cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &

! Marine

         map_prop_hc, map_prop_hist, &

         map_emis_hc,map_iwp_hc,map_deltaz_hc,&

         map_pcld_hc,map_tcld_hc,&

         map_emis_hist,map_iwp_hist,map_deltaz_hist,&

         map_rad_hist,&

         map_ntot,map_hc,map_hist,&

         map_Cb,map_ThCi,map_Anv,&

         map_emis_Cb,map_pcld_Cb,map_tcld_Cb,&

         map_emis_ThCi,map_pcld_ThCi,map_tcld_ThCi,&

         map_emis_Anv,map_pcld_Anv,map_tcld_Anv, &

         alt_tropo, &

!Ionela

         ok_4xCO2atm, tkt, tks, taur, sss

    USE ocean_slab_mod, ONLY: nslay, tslab, slab_bilg, tice, seaice, &

        slab_ekman,slab_hdiff,slab_gm,dt_ekman, dt_hdiff, dt_gm, dt_qflux

    USE pbl_surface_mod, ONLY: snow

    USE indice_sol_mod, ONLY: nbsrf

#ifdef ISO

    USE isotopes_mod, ONLY: iso_HTO

#endif

    USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg

    USE surface_data, ONLY: type_ocean, version_ocean, ok_veget, landice_opt

    USE aero_mod, ONLY: naero_tot, id_STRAT_phy

    USE ioipsl, ONLY: histend, histsync

    USE iophy, ONLY: set_itau_iophy, histwrite_phy

    USE netcdf, ONLY: nf90_fill_real

    USE print_control_mod, ONLY: prt_level,lunout

#ifdef CPP_XIOS

    ! ug Pour les sorties XIOS

    USE xios

    USE wxios, ONLY: wxios_closedef, missing_val, wxios_set_context

#endif

    USE phys_cal_mod, ONLY : mth_len

#ifdef CPP_RRTM

    USE YOESW, ONLY : RSUN

#endif

    USE tracinca_mod, ONLY: config_inca

    use config_ocean_skin_m, only: activate_ocean_skin

    USE vertical_layers_mod, ONLY: presnivs

    IMPLICIT NONE

    INCLUDE "clesphys.h"

    INCLUDE "alpale.h"

    INCLUDE "compbl.h"

    INCLUDE "YOMCST.h"

    ! Input

    INTEGER :: itap, ivap, iliq, isol, ibs, read_climoz

    INTEGER, DIMENSION(klon) :: lmax_th

    LOGICAL :: aerosol_couple, ok_sync

    LOGICAL :: ok_ade, ok_aie, ok_volcan

    LOGICAL, DIMENSION(klon, klev) :: ptconv, ptconvth

    REAL :: pdtphys

    CHARACTER (LEN=4), DIMENSION(nlevSTD) :: clevSTD

    REAL, DIMENSION(klon) :: pphis

    REAL, DIMENSION(klon, klev) :: pplay, d_u, d_t

    REAL, DIMENSION(klon, klev+1) :: paprs

    REAL, DIMENSION(klon,klev,nqtot) :: qx, d_qx

    REAL, DIMENSION(klon, klev) :: zmasse

    INTEGER :: flag_aerosol_strat

    INTEGER :: flag_aerosol

    LOGICAL :: ok_cdnc

    REAL, DIMENSION(3) :: freq_moyNMC

    ! Local

    INTEGER :: itau_w

    INTEGER :: i, iinit, iinitend=1, iff, iq, nsrf, k, ll, naero

    REAL, DIMENSION (klon,NSW) :: zx_tmp_fi3dsp

    CHARACTER (LEN=4)              :: bb2

    REAL, PARAMETER :: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2

!   REAL, PARAMETER :: missing_val=nf90_fill_real

    INTEGER :: ilev

    INTEGER, SAVE :: kmax_100m

!$OMP THREADPRIVATE(kmax_100m)

    REAL :: x

#ifndef CPP_XIOS

    REAL :: missing_val

#endif

    REAL, PARAMETER :: un_jour=86400.

    INTEGER :: ISW, itr, ixt, it

    CHARACTER*1 ch1

    CHARACTER(LEN=maxlen) :: varname, dn

#ifdef CPP_XIOS

    TYPE(xios_fieldgroup) :: group_handle

    TYPE(xios_field) :: child

#endif

#ifdef CPP_StratAer

    LOGICAL, PARAMETER :: debug_strataer=.FALSE.

    CHARACTER(LEN=maxlen) :: unt

#endif

    ! On calcul le nouveau tau:

    ! On le donne à iophy pour que les histwrite y aient accès:

 !   IF (.NOT.vars_defined) THEN

 !   ELSE

 !      iinitend = 1

 !   ENDIF

#ifdef CPP_XIOS

    CALL wxios_set_context

#endif

#ifndef CPP_XIOS

    missing_val=missing_val_nf90

#endif

      ENDDO

    ENDIF

    ENDDO

    END DO

#ifdef CPP_XIOS

#ifdef CPP_StratAer

!$OMP MASTER

   IF (.NOT.vars_defined) THEN

          !On ajoute les variables 3D traceurs par l interface fortran

          CALL xios_get_handle("fields_strataer_trac_3D", group_handle)

          ! On boucle sur les traceurs pour les ajouter au groupe puis fixer les attributs

          DO iq = 1, nqtot

            IF(.NOT.(tracers(iq)%isAdvected .AND. tracers(iq)%isInPhysics)) CYCLE

            dn = 'd'//TRIM(tracers(iq)%name)//'_'

            WRITE (lunout,*) 'XIOS var=', nqo, iq, nqtot, tracers(iq)%name

            unt = "kg kg-1"

            varname=trim(tracers(iq)%name)

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            unt = "kg kg-1 s-1"

            varname=TRIM(dn)//'vdf'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'the'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'con'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'lessi_impa'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'lessi_nucl'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'insc'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'bcscav'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'evapls'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'ls'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'trsp'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'sscav'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'sat'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            varname=TRIM(dn)//'uscav'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

          END DO

          !On ajoute les variables 2D traceurs par l interface fortran

          CALL xios_get_handle("fields_strataer_trac_2D", group_handle)

          ! On boucle sur les traceurs pour les ajouter au groupe puis fixer les attributs

          DO iq = 1, nqtot

            IF(.NOT.(tracers(iq)%isAdvected .AND. tracers(iq)%isInPhysics)) CYCLE

            unt = "kg m-2"

            varname='cum'//trim(tracers(iq)%name)

            WRITE (lunout,*) 'XIOS var=', iq, nqtot, varname

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

            unt = "kg m-2 s-1"

            varname='cumd'//trim(tracers(iq)%name)//'_dry'

            CALL xios_add_child(group_handle, child, varname)

            CALL xios_set_attr(child, name=varname, unit=unt)

          ENDDO

    ENDIF

!$OMP END MASTER

#endif

#endif

    ! ug la boucle qui suit ne sert qu'une fois, pour l'initialisation, sinon il n'y a toujours qu'un seul passage:

!      print *,'IFF iinit=', iinit, iinitend

#ifdef CPP_XIOS

       !$OMP MASTER

       IF (vars_defined) THEN

          IF (prt_level >= 10) then

             write(lunout,*)"phys_output_write: call xios_update_calendar, itau_w=",itau_w

          ENDIF

!          CALL xios_update_calendar(itau_w)

          CALL xios_update_calendar(itap)

       ENDIF

       !$OMP END MASTER

       !$OMP BARRIER

#endif

       ! On procède à l'écriture ou à la définition des nombreuses variables:

!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

       ENDIF

       ENDIf

          ENDDO

       ENDIF

!

#ifdef CPP_XIOS

       CALL histwrite_phy("R_ecc",R_ecc)

       CALL histwrite_phy("R_peri",R_peri)

       CALL histwrite_phy("R_incl",R_incl)

       CALL histwrite_phy("solaire",solaire)

       CALL histwrite_phy(o_Ahyb, ap)

       CALL histwrite_phy(o_Bhyb, bp)

       CALL histwrite_phy(o_Ahyb_bounds, Ahyb_bounds)

       CALL histwrite_phy(o_Bhyb_bounds, Bhyb_bounds)

       CALL histwrite_phy(o_Ahyb_mid, aps)

       CALL histwrite_phy(o_Bhyb_mid, bps)

       CALL histwrite_phy(o_Ahyb_mid_bounds, Ahyb_mid_bounds)

       CALL histwrite_phy(o_Bhyb_mid_bounds, Bhyb_mid_bounds)

       CALL histwrite_phy(o_longitude, longitude_deg)

       CALL histwrite_phy(o_latitude, latitude_deg)

!

#ifdef CPP_RRTM

      IF (iflag_rrtm.EQ.1) THEN

        DO ISW=1, NSW

          WRITE(ch1,'(i1)') ISW

!         zx_tmp_0d=RSUN(ISW)

!         CALL histwrite_phy("rsun"//ch1,zx_tmp_0d)

          CALL histwrite_phy("rsun"//ch1,RSUN(ISW))

        ENDDO

      ENDIF

#endif

!

       CALL histwrite_phy("co2_ppm",co2_ppm)

       CALL histwrite_phy("CH4_ppb",CH4_ppb)

       CALL histwrite_phy("N2O_ppb",N2O_ppb)

       CALL histwrite_phy("CFC11_ppt",CFC11_ppt)

       CALL histwrite_phy("CFC12_ppt",CFC12_ppt)

!

#endif

!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Simulateur AIRS

     ENDIF

          ENDDO

       ENDIF

          ENDDO

       ENDIF

          ENDDO

       ENDIF

          ENDDO

       ENDIF

          ENDDO

       ENDIF

          ENDDO

       ENDIF

          ENDDO

       ENDIF

          ENDDO

       ENDIF

            ELSE

            ENDIF

          ENDDO

                wind100m(i)=SQRT( (u_seri(i,k)+(100.-z(i,k))/(z(i,k+1)-z(i,k))*(u_seri(i,k+1)-u_seri(i,k)))**2.0 + &

              ENDIF

            ENDDO

          ENDDO

       ENDIF

         !--polynomial fit for 14,Vestas,1074,V136/3450 kW windmill - Olivier

             x=wind100m(i)

             ELSE

               zx_tmp_fi2d(i)= 10.73 + x*(-14.69 + x*(8.339 + x*(-2.59 + x*(0.4893 + x*(-0.05898 + x*(0.004627 + &

             ENDIF

           ELSE

           ENDIF

         ENDDO

       ENDIF

         !--polynomial fit for 14,Vestas,867,V164/8000 kW - Olivier

             x=wind100m(i)

             ELSE

               zx_tmp_fi2d(i)= 20.59 + x*(-22.39 + x*(10.25 + x*(-2.601 + x*(0.4065 + x*(-0.04099 + x*(0.002716 + &

             ENDIF

           ELSE

           ENDIF

         ENDDO

       ENDIF

          ENDDO

       ENDIF

          ENDDO

       ENDIF

       ENDIF

             ENDDO

          ELSE

             ENDDO

          ENDIF

       ENDIF

       ! epmax_cape:

!       CALL histwrite_phy(o_epmax_diag, epmax_diag)

          ENDDO

       ENDIF

          ENDDO

       ENDIF

       ENDIF

       ENDIF

       ENDIF

       ENDIF

       ENDIF

          ENDIF