Context Navigation

phys_output_write_mod.F90 @ 5447

Last change on this file since 5447 was 5447, checked in by jyg, 10 hours ago
output the convective wet deposit of tracers
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:keywords set to `Id`
File size: 121.9 KB

Rev	Line
[1798]	1	!
[2429]	2	! $Id: phys_output_write_mod.F90 5447 2024-12-22 16:48:23Z jyg $
[1798]	3	!
	4	MODULE phys_output_write_mod
	5
[1898]	6	USE phytrac_mod, ONLY : d_tr_cl, d_tr_th, d_tr_cv, d_tr_lessi_impa, &
	7	d_tr_lessi_nucl, d_tr_insc, d_tr_bcscav, d_tr_evapls, d_tr_ls, &
[5447]	8	d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav, flux_tr_wet, flux_tr_dry
[1813]	9
[1898]	10	! Author: Abderrahmane IDELKADI (original include file)
	11	! Author: Laurent FAIRHEAD (transformation to module/subroutine)
	12	! Author: Ulysse GERARD (effective implementation)
[1813]	13
[1898]	14	CONTAINS
[1798]	15
[3121]	16	! ug Routine pour définir (lors du premier passage) ET sortir les variables
[1898]	17	SUBROUTINE phys_output_write(itap, pdtphys, paprs, pphis, &
	18	pplay, lmax_th, aerosol_couple, &
[4523]	19	ok_ade, ok_aie, ok_volcan, ivap, iliq, isol, ibs, ok_sync, &
[1915]	20	ptconv, read_climoz, clevSTD, ptconvth, &
[4608]	21	d_u, d_t, qx, d_qx, zmasse, flag_aerosol, flag_aerosol_strat, ok_cdnc, t, u1, v1)
[1798]	22
[1898]	23	! This subroutine does the actual writing of diagnostics that were
	24	! defined and initialised in phys_output_mod.F90
[1798]	25
[2704]	26	USE dimphy, ONLY: klon, klev, klevp1
[4389]	27	USE infotrac_phy, ONLY: nbtr, nqtot, nqo, type_trac, tracers, niso, ntiso
[4046]	28	USE strings_mod, ONLY: maxlen
[2429]	29	USE mod_phys_lmdz_para, ONLY: is_north_pole_phy,is_south_pole_phy
[2344]	30	USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
[2704]	31	USE time_phylmdz_mod, ONLY: day_step_phy, start_time, itau_phy
[3003]	32	USE vertical_layers_mod, ONLY : ap, bp, aps, bps
[4608]	33	USE phystokenc_mod, ONLY: offline
[2704]	34	USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, &
[3059]	35	o_longitude, o_latitude, &
[3474]	36	o_Ahyb, o_Bhyb,o_Ahyb_bounds, o_Bhyb_bounds, &
	37	o_Ahyb_mid, o_Bhyb_mid,o_Ahyb_mid_bounds, o_Bhyb_mid_bounds, &
[1898]	38	is_ave, is_sic, o_contfracATM, o_contfracOR, &
[2385]	39	o_aireTER, o_flat, o_slp, o_ptstar, o_pt0, o_tsol, &
[1898]	40	o_t2m, o_t2m_min, o_t2m_max, &
[2103]	41	o_t2m_min_mon, o_t2m_max_mon, &
[3817]	42	o_nt2mout, o_nt2moutfg, &
	43	o_nq2mout, o_nq2moutfg, &
	44	o_nu2mout, o_nu2moutfg, &
[1898]	45	o_q2m, o_ustar, o_u10m, o_v10m, &
[3622]	46	o_wind10m, o_wind10max, o_wind100m, o_gusts, o_sicf, &
[3691]	47	o_loadfactor_wind_onshore, o_loadfactor_wind_offshore, &
[1898]	48	o_psol, o_mass, o_qsurf, o_qsol, &
[3219]	49	o_precip, o_rain_fall, o_rain_con, o_ndayrain, o_plul, o_pluc, o_plun, &
[4835]	50	o_snow, o_msnow, o_fsnow, o_evap, o_snowerosion, o_ustart_lic, o_qsalt_lic, o_rhosnow_lic, o_bsfall, &
[4523]	51	o_ep,o_epmax_diag, & ! epmax_cape
[1898]	52	o_tops, o_tops0, o_topl, o_topl0, &
[3756]	53	o_SWupTOA, o_SWupTOAclr, o_SWupTOAcleanclr, o_SWdnTOA, o_fdiffSWdnSFC, &
[1898]	54	o_SWdnTOAclr, o_nettop, o_SWup200, &
	55	o_SWup200clr, o_SWdn200, o_SWdn200clr, &
[5204]	56	o_LWupTOA, o_LWupTOAclr, &
[1898]	57	o_LWup200, o_LWup200clr, o_LWdn200, &
	58	o_LWdn200clr, o_sols, o_sols0, &
	59	o_soll, o_radsol, o_soll0, o_SWupSFC, &
[3082]	60	o_SWupSFCclr, o_SWupSFCcleanclr, o_SWdnSFC, o_SWdnSFCclr, o_SWdnSFCcleanclr, &
[1898]	61	o_LWupSFC, o_LWdnSFC, o_LWupSFCclr, &
[3106]	62	o_LWdnSFCclr, o_LWupTOAcleanclr, o_LWdnSFCcleanclr, o_bils, o_bils_diss, &
[2042]	63	o_bils_ec,o_bils_ech, o_bils_tke, o_bils_kinetic, &
[1898]	64	o_bils_latent, o_bils_enthalp, o_sens, &
[3095]	65	o_fder, o_ffonte, o_fqcalving, o_fqfonte, o_mrroli, o_runofflic, &
[1898]	66	o_taux, o_tauy, o_snowsrf, o_qsnow, &
[5109]	67	! SN runoff_diag
	68	o_snowhgt, o_toice, o_sissnow, o_runoff, o_runoff_diag, &
[1898]	69	o_albslw3, o_pourc_srf, o_fract_srf, &
	70	o_taux_srf, o_tauy_srf, o_tsol_srf, &
	71	o_evappot_srf, o_ustar_srf, o_u10m_srf, &
	72	o_v10m_srf, o_t2m_srf, o_evap_srf, &
	73	o_sens_srf, o_lat_srf, o_flw_srf, &
	74	o_fsw_srf, o_wbils_srf, o_wbilo_srf, &
[2670]	75	o_wevap_srf, o_wrain_srf, o_wsnow_srf, &
[2159]	76	o_tke_srf, o_tke_max_srf,o_dltpbltke_srf, o_wstar, &
[2561]	77	o_l_mixmin,o_l_mix, &
[1898]	78	o_cdrm, o_cdrh, o_cldl, o_cldm, o_cldh, &
[1924]	79	o_cldt, o_JrNt, o_cldljn, o_cldmjn, &
	80	o_cldhjn, o_cldtjn, o_cldq, o_lwp, o_iwp, &
	81	o_ue, o_ve, o_uq, o_vq, o_cape, o_pbase, &
[3257]	82	o_uwat, o_vwat, &
[1898]	83	o_ptop, o_fbase, o_plcl, o_plfc, &
[3496]	84	o_wbeff, o_convoccur, o_cape_max, o_upwd, o_ep,o_epmax_diag, &
	85	o_Mipsh, o_Ma, &
[3171]	86	o_dnwd, o_dnwd0, o_ftime_deepcv, o_ftime_con, o_mc, &
[4887]	87	o_prw, o_prlw, o_prsw, o_prbsw, o_water_budget, o_s_pblh, o_s_pblt, o_s_lcl, &
[1898]	88	o_s_therm, o_uSTDlevs, o_vSTDlevs, &
	89	o_wSTDlevs, o_zSTDlevs, o_qSTDlevs, &
	90	o_tSTDlevs, epsfra, o_t_oce_sic, &
	91	o_ale_bl, o_alp_bl, o_ale_wk, o_alp_wk, &
[2393]	92	o_dtvdf_x , o_dtvdf_w , o_dqvdf_x , o_dqvdf_w , &
	93	o_sens_x , o_sens_w , o_flat_x , o_flat_w , &
[3888]	94	o_delta_tsurf, o_delta_tsurf_srf, &
[2393]	95	o_cdragh_x , o_cdragh_w , o_cdragm_x , o_cdragm_w , &
	96	o_kh , o_kh_x , o_kh_w , &
[3208]	97	o_ale, o_alp, o_cin, o_WAPE, o_wake_h, o_cv_gen, o_wake_dens, &
[1898]	98	o_wake_s, o_wake_deltat, o_wake_deltaq, &
[2486]	99	o_wake_omg, o_dtwak, o_dqwak, o_dqwak2d, o_Vprecip, &
[3496]	100	o_qtaa, o_Clwaa, &
	101	o_ftd, o_fqd, o_wdtrainA, o_wdtrainS, o_wdtrainM, &
[4843]	102	o_n2, o_s2, o_strig, o_zcong, o_zlcl_th, o_proba_notrig, &
[1898]	103	o_random_notrig, o_ale_bl_stat, &
	104	o_ale_bl_trig, o_alp_bl_det, &
	105	o_alp_bl_fluct_m, o_alp_bl_fluct_tke, &
[2057]	106	o_alp_bl_conv, o_alp_bl_stat, &
	107	o_slab_qflux, o_tslab, o_slab_bils, &
[2209]	108	o_slab_bilg, o_slab_sic, o_slab_tice, &
[3002]	109	o_slab_hdiff, o_slab_ekman, o_slab_gm, &
[1898]	110	o_weakinv, o_dthmin, o_cldtau, &
	111	o_cldemi, o_pr_con_l, o_pr_con_i, &
[4523]	112	o_pr_lsc_l, o_pr_lsc_i, o_pr_bs, o_re, o_fl, &
[3165]	113	o_rh2m, &
[1898]	114	o_qsat2m, o_tpot, o_tpote, o_SWnetOR, &
[3332]	115	o_LWdownOR, o_snowl, &
[1898]	116	o_solldown, o_dtsvdfo, o_dtsvdft, &
[2952]	117	o_dtsvdfg, o_dtsvdfi, o_z0m, o_z0h, o_od443aer, o_od550aer, &
[2854]	118	o_dryod550aer, o_od865aer, o_abs550aer, o_od550lt1aer, &
[2146]	119	o_sconcso4, o_sconcno3, o_sconcoa, o_sconcbc, &
	120	o_sconcss, o_sconcdust, o_concso4, o_concno3, &
[1898]	121	o_concoa, o_concbc, o_concss, o_concdust, &
	122	o_loadso4, o_loadoa, o_loadbc, o_loadss, &
[2854]	123	o_loaddust, o_loadno3, o_tausumaero, &
	124	o_drytausumaero, o_tausumaero_lw, &
[2146]	125	o_topswad, o_topswad0, o_solswad, o_solswad0, &
	126	o_toplwad, o_toplwad0, o_sollwad, o_sollwad0, &
[1898]	127	o_swtoaas_nat, o_swsrfas_nat, &
	128	o_swtoacs_nat, o_swtoaas_ant, &
	129	o_swsrfas_ant, o_swtoacs_ant, &
	130	o_swsrfcs_ant, o_swtoacf_nat, &
	131	o_swsrfcf_nat, o_swtoacf_ant, &
	132	o_swsrfcs_nat, o_swsrfcf_ant, &
	133	o_swtoacf_zero, o_swsrfcf_zero, &
[2780]	134	o_topswai, o_solswai, o_toplwai, o_sollwai, o_scdnc, &
[3117]	135	o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,&
[1898]	136	o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
[3121]	137	o_lcc3dstra, o_icc3dcon, o_icc3dstra, &
[3124]	138	o_cldicemxrat, o_cldwatmxrat, o_reffclwtop, o_ec550aer, &
[1898]	139	o_lwcon, o_iwcon, o_temp, o_theta, &
[5383]	140	o_ovapinit, o_ovap, o_oliq, o_ocond, o_oice, o_geop,o_qbs, &
[1898]	141	o_vitu, o_vitv, o_vitw, o_pres, o_paprs, &
[1924]	142	o_zfull, o_zhalf, o_rneb, o_rnebjn, o_rnebcon, &
[3780]	143	o_rnebls, o_rneblsvol, o_rhum, o_rhl, o_rhi, o_ozone, o_ozone_light, &
[5007]	144	o_pfraclr, o_pfracld, o_cldfraliq, o_sigma2_icefracturb, o_mean_icefracturb, &
[4819]	145	o_qrainlsc, o_qsnowlsc, o_dqreva, o_dqrauto, o_dqrcol, o_dqrmelt, o_dqrfreez, &
[4803]	146	o_dqssub, o_dqsauto, o_dqsagg, o_dqsrim, o_dqsmelt, o_dqsfreez, &
[2665]	147	o_duphy, o_dtphy, o_dqphy, o_dqphy2d, o_dqlphy, o_dqlphy2d, &
[4523]	148	o_dqsphy, o_dqsphy2d, o_dqbsphy, o_dqbsphy2d, o_albe_srf, o_z0m_srf, o_z0h_srf, &
[3780]	149	o_ages_srf, o_snow_srf, o_alb1, o_alb2, o_tke, o_tke_dissip, &
[5039]	150	o_tke_max, o_kz, o_kz_max, o_clwcon, o_tke_shear, o_tke_buoy, o_tke_trans, &
[2496]	151	o_dtdyn, o_dqdyn, o_dqdyn2d, o_dqldyn, o_dqldyn2d, &
[4523]	152	o_dqsdyn, o_dqsdyn2d, o_dqbsdyn, o_dqbsdyn2d, o_dudyn, o_dvdyn, &
[1898]	153	o_dtcon, o_tntc, o_ducon, o_dvcon, &
[2486]	154	o_dqcon, o_dqcon2d, o_tnhusc, o_tnhusc, o_dtlsc, &
	155	o_dtlschr, o_dqlsc, o_dqlsc2d, o_beta_prec, &
	156	o_dtlscth, o_dtlscst, o_dqlscth, o_dqlscth2d, &
	157	o_dqlscst, o_dqlscst2d, o_plulth, o_plulst, &
[1898]	158	o_ptconvth, o_lmaxth, o_dtvdf, &
[2486]	159	o_dtdis, o_dqvdf, o_dqvdf2d, o_dteva, o_dqeva, o_dqeva2d, &
[4523]	160	o_dqbsvdf, o_dtbs, o_dqbs, o_dqbsbs, &
[1903]	161	o_ptconv, o_ratqs, o_dtthe, &
	162	o_duthe, o_dvthe, o_ftime_th, &
[1898]	163	o_f_th, o_e_th, o_w_th, o_q_th, &
[2958]	164	o_a_th, o_cloudth_sth, o_cloudth_senv, &
	165	o_cloudth_sigmath, o_cloudth_sigmaenv, &
	166	o_d_th, o_f0_th, o_zmax_th, &
[2486]	167	o_dqthe, o_dqthe2d, o_dtajs, o_dqajs, o_dqajs2d, o_dtswr, &
[1898]	168	o_dtsw0, o_dtlwr, o_dtlw0, o_dtec, &
	169	o_duvdf, o_dvvdf, o_duoro, o_dvoro, &
	170	o_dtoro, o_dulif, o_dvlif, o_dtlif, &
[2333]	171	o_du_gwd_hines, o_dv_gwd_hines, o_dthin, o_dqch4, o_rsu, &
	172	o_du_gwd_front, o_dv_gwd_front, &
	173	o_east_gwstress, o_west_gwstress, &
[3082]	174	o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, o_rsucsaf, o_rsdcsaf, &
[1898]	175	o_rlucs, o_rldcs, o_tnt, o_tntr, &
	176	o_tntscpbl, o_tnhus, o_tnhusscpbl, &
	177	o_evu, o_h2o, o_mcd, o_dmc, o_ref_liq, &
	178	o_ref_ice, o_rsut4co2, o_rlut4co2, &
	179	o_rsutcs4co2, o_rlutcs4co2, o_rsu4co2, &
	180	o_rlu4co2, o_rsucs4co2, o_rlucs4co2, &
	181	o_rsd4co2, o_rld4co2, o_rsdcs4co2, &
	182	o_rldcs4co2, o_tnondef, o_ta, o_zg, &
	183	o_hus, o_hur, o_ua, o_va, o_wap, &
	184	o_psbg, o_tro3, o_tro3_daylight, &
	185	o_uxv, o_vxq, o_vxT, o_wxq, o_vxphi, &
	186	o_wxT, o_uxu, o_vxv, o_TxT, o_trac, &
	187	o_dtr_vdf, o_dtr_the, o_dtr_con, &
	188	o_dtr_lessi_impa, o_dtr_lessi_nucl, &
	189	o_dtr_insc, o_dtr_bcscav, o_dtr_evapls, &
[2704]	190	o_dtr_ls, o_dtr_trsp, o_dtr_sscav, o_dtr_dry, &
[5447]	191	o_dtr_sat, o_dtr_uscav, o_dtr_wet_con, &
	192	o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
[2333]	193	o_ustr_gwd_hines,o_vstr_gwd_hines,o_ustr_gwd_rando,o_vstr_gwd_rando, &
[2538]	194	o_ustr_gwd_front,o_vstr_gwd_front, &
	195	o_sens_prec_liq_oce, o_sens_prec_liq_sic, &
	196	o_sens_prec_sol_oce, o_sens_prec_sol_sic, &
	197	o_lat_prec_liq_oce, o_lat_prec_liq_sic, &
[2580]	198	o_lat_prec_sol_oce, o_lat_prec_sol_sic, &
[3110]	199	o_sza, &
[2580]	200	! Marine
	201	o_map_prop_hc, o_map_prop_hist, o_map_emis_hc, o_map_iwp_hc, &
	202	o_map_deltaz_hc, o_map_pcld_hc, o_map_tcld_hc, &
	203	o_map_emis_hist, o_map_iwp_hist, o_map_deltaz_hist, &
	204	o_map_rad_hist, &
	205	o_map_emis_Cb, o_map_pcld_Cb, o_map_tcld_Cb, &
	206	o_map_emis_ThCi, o_map_pcld_ThCi, o_map_tcld_ThCi, &
	207	o_map_emis_Anv, o_map_pcld_Anv, o_map_tcld_Anv, &
	208	o_map_ntot, o_map_hc,o_map_hist,o_map_Cb,o_map_ThCi,o_map_Anv, &
[4056]	209	#ifdef ISO
	210	! Isotopes
	211	o_xtprecip,o_xtplul,o_xtpluc,o_xtovap,o_xtoliq,o_xtcond, &
[5109]	212	o_xtrunoff_diag, &
[4056]	213	o_xtevap,o_dxtdyn,o_dxtldyn,o_dxtcon,o_dxtlsc,o_dxteva, &
	214	o_dxtajs,o_dxtvdf,o_dxtthe, o_dxtch4, &
	215	o_dxtprod_nucl,o_dxtcosmo,o_dxtdecroiss, &
[4525]	216	o_xtevap_srf, &
[4056]	217	#endif
	218	! Tropopause
[2992]	219	o_alt_tropo, &
[2993]	220	o_p_tropopause, o_z_tropopause, o_t_tropopause, &
[3453]	221	o_col_O3_strato, o_col_O3_tropo, &
[5204]	222	!-- LSCP - condensation and ice supersaturation variables
	223	o_cfseri, o_dcfdyn, o_rvcseri, o_drvcdyn, &
	224	o_qsub, o_qissr, o_qcld, o_subfra, o_issrfra, o_gammacond, &
	225	o_dcfsub, o_dcfcon, o_dcfmix, o_dqiadj, o_dqisub, o_dqicon, o_dqimix, &
	226	o_dqvcadj, o_dqvcsub, o_dqvccon, o_dqvcmix, o_qsatl, o_qsati, &
	227	!-- LSCP - aviation variables
[4059]	228	o_Tcontr, o_qcontr, o_qcontr2, o_fcontrN, o_fcontrP, &
[5204]	229	o_dcfavi, o_dqiavi, o_dqvcavi, o_flight_dist, o_flight_h2o, &
[3453]	230	!--interactive CO2
[3857]	231	o_flx_co2_ocean, o_flx_co2_ocean_cor, &
	232	o_flx_co2_land, o_flx_co2_land_cor, &
	233	o_flx_co2_ff, o_flx_co2_bb, &
[3815]	234	o_delta_sst, o_delta_sal, o_ds_ns, o_dt_ns, o_dter, o_dser, o_tkt, &
[4575]	235	o_tks, o_taur, o_sss, &
	236	!FC
	237	o_zxfluxt,o_zxfluxq
[1798]	238
[4677]	239	#ifdef CPP_ECRAD
	240	USE phys_output_ctrlout_mod, ONLY: &
	241	o_soll0_s2,o_soll_s2,o_sols0_s2,o_sols_s2, &
[4773]	242	o_topl0_s2,o_topl_s2,o_tops0_s2,o_tops_s2, &
	243	o_SWupTOA_s2,o_SWupTOAclr_s2,o_cloud_cover_sw, &
	244	o_cloud_cover_sw_s2,o_SWdnTOA_s2,o_SWdnTOAclr_s2, &
	245	o_LWupSFCclr_s2, o_LWdnSFCclr_s2, o_SWupSFC_s2, &
	246	o_SWupSFCclr_s2, o_SWdnSFC_s2, o_SWdnSFCclr_s2, &
[4810]	247	o_LWupSFC_s2, o_LWdnSFC_s2, o_rlu_s2, o_rld_s2, &
[4808]	248	o_rlucs_s2, o_rldcs_s2, o_rsu_s2, o_rsd_s2, &
	249	o_rsucs_s2, o_rsdcs_s2
[4677]	250	#endif
	251
[4998]	252	USE infotrac_phy, ONLY: nbtr_bin
[5252]	253	USE phys_output_ctrlout_mod, ONLY: &
[2752]	254	o_budg_3D_nucl, o_budg_3D_cond_evap, o_budg_3D_ocs_to_so2, o_budg_3D_so2_to_h2so4, &
	255	o_budg_sed_part, o_R2SO4, o_OCS_lifetime, o_SO2_lifetime, &
[5252]	256	o_budg_3D_backgr_ocs, o_budg_3D_backgr_so2, &
[3100]	257	o_budg_dep_dry_ocs, o_budg_dep_wet_ocs, &
	258	o_budg_dep_dry_so2, o_budg_dep_wet_so2, &
[2752]	259	o_budg_dep_dry_h2so4, o_budg_dep_wet_h2so4, &
[5252]	260	o_budg_dep_dry_part, o_budg_dep_wet_part, &
	261	o_budg_emi_ocs, o_budg_emi_so2, o_budg_emi_h2so4, o_budg_emi_part, &
[2752]	262	o_budg_ocs_to_so2, o_budg_so2_to_h2so4, o_budg_h2so4_to_part, &
[2690]	263	o_surf_PM25_sulf, o_ext_strat_550, o_tau_strat_550, &
[4998]	264	o_vsed_aer, o_tau_strat_1020, o_ext_strat_1020, o_f_r_wet, &
[5150]	265	o_SAD_sulfate, o_reff_sulfate, o_sulfmmr, o_nd_mode, o_sulfmmr_mode
[4059]	266
[4803]	267	USE lmdz_lscp_ini, ONLY: ok_poprecip
	268
[3479]	269	USE phys_output_ctrlout_mod, ONLY: o_heat_volc, o_cool_volc !NL
	270	USE phys_state_var_mod, ONLY: heat_volc, cool_volc !NL
[2580]	271
[4523]	272	USE phys_state_var_mod, ONLY: pctsrf, rain_fall, snow_fall, bs_fall,&
[2243]	273	qsol, z0m, z0h, fevap, agesno, &
[4556]	274	nday_rain, ndayrain_mth, rain_con, snow_con, &
[3756]	275	topsw, toplw, toplw0, swup, swdn, solswfdiff, &
[3082]	276	topsw0, swupc0, swdnc0, swup0, swdn0, SWup200, SWup200clr, &
[1898]	277	SWdn200, SWdn200clr, LWup200, LWup200clr, &
	278	LWdn200, LWdn200clr, solsw, solsw0, sollw, &
[2366]	279	radsol, swradcorr, sollw0, sollwdown, sollw, gustiness, &
[3106]	280	sollwdownclr, lwdnc0, lwdn0, ftsol, ustar, u10m, &
[2159]	281	v10m, pbl_tke, wake_delta_pbl_TKE, &
[2393]	282	delta_tsurf, &
[2159]	283	wstar, cape, ema_pcb, ema_pct, &
[3496]	284	ema_cbmf, Mipsh, Ma, fm_therm, ale_bl, alp_bl, ale, &
[3956]	285	alp, cin, wake_pe, wake_dens, cv_gen, wake_s, wake_deltat, &
[1898]	286	wake_deltaq, ftd, fqd, ale_bl_trig, albsol1, &
[3080]	287	ale_wake, ale_bl_stat, &
[1898]	288	rnebcon, wo, falb1, albsol2, coefh, clwcon0, &
[2194]	289	ratqs, entr_therm, zqasc, detr_therm, f0, &
[3106]	290	lwup, lwdn, lwupc0, lwup0, coefm, &
	291	swupp, lwupp, swupc0p, swup0p, lwupc0p, lwup0p, swdnp, lwdnp, &
	292	swdnc0p, swdn0p, lwdnc0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
[1898]	293	vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &
	294	vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &
[2333]	295	T2sumSTD, nlevSTD, du_gwd_rando, du_gwd_front, &
[2271]	296	ulevSTD, vlevSTD, wlevSTD, philevSTD, qlevSTD, tlevSTD, &
[3815]	297	rhlevSTD, O3STD, O3daySTD, uvSTD, vqSTD, vTSTD, wqSTD, vphiSTD, &
	298	wTSTD, u2STD, v2STD, T2STD, missing_val_nf90, delta_sal, ds_ns, &
[4056]	299	#ifdef ISO
[4737]	300	xtrain_con, xtsnow_con, xtrain_fall, xtsnow_fall, fxtevap, &
[4056]	301	#endif
[4370]	302	dt_ns, delta_sst, dter, dser
[4677]	303
	304	! AI 08 2023 pour ECRAD 3Deffect
	305	#ifdef CPP_ECRAD
	306	USE phys_state_var_mod, ONLY: &
	307	sollw0_s2,sollw_s2,solsw0_s2,solsw_s2, &
[4773]	308	toplw0_s2,toplw_s2,topsw0_s2,topsw_s2, &
	309	toplw0_s2,toplw_s2,topsw0_s2,topsw_s2, &
	310	swup0_s2,swup_s2,swdn_s2,swdn0_s2,sollwdownclr_s2, &
[4808]	311	sollwdown_s2,lwdn0_s2,lwup_s2,lwdn_s2,lwup0_s2,lwdn0_s2
[4773]	312	USE phys_output_var_mod, ONLY: cloud_cover_sw, &
	313	cloud_cover_sw_s2
[4677]	314	#endif
[1898]	315
[2704]	316	USE phys_local_var_mod, ONLY: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, &
[3817]	317	zn2mout, t2m_min_mon, t2m_max_mon, evap, &
[4835]	318	snowerosion, zxustartlic, zxrhoslic, zxqsaltlic, &
[5039]	319	l_mixmin,l_mix, pbl_eps, tke_shear, tke_buoy, tke_trans, &
[2243]	320	zu10m, zv10m, zq2m, zustar, zxqsurf, &
[2516]	321	rain_lsc, rain_num, snow_lsc, bils, sens, fder, &
[2517]	322	zxffonte, zxfqcalving, zxfqfonte, zxrunofflic, fluxu, &
[1898]	323	fluxv, zxsnow, qsnow, snowhgt, to_ice, &
[5109]	324	! SN runoff_diag
	325	sissnow, runoff, runoff_diag, albsol3_lic, evap_pot, &
[2243]	326	t2m, fluxt, fluxlat, fsollw, fsolsw, &
[4737]	327	wfbils, wfevap, &
[2670]	328	cdragm, cdragh, cldl, cldm, &
[3779]	329	cldh, cldt, JrNt, & ! only output names: cldljn,cldmjn,cldhjn,cldtjn
	330	cldq, flwp, fiwp, ue, ve, uq, vq, &
[3257]	331	uwat, vwat, &
[4887]	332	plcl, plfc, wbeff, convoccur, upwd, dnwd, dnwd0, prw, prlw, prsw, prbsw, water_budget, &
[1898]	333	s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
	334	vwriteSTD, wwriteSTD, phiwriteSTD, qwriteSTD, &
[3080]	335	twriteSTD, alp_wake, &
[3179]	336	!! dtvdf_x ,dtvdf_w ,dqvdf_x ,dqvdf_w , &
	337	d_t_vdf_x ,d_t_vdf_w ,d_q_vdf_x ,d_q_vdf_w , &
[2393]	338	sens_x ,sens_w ,zxfluxlat_x,zxfluxlat_w, &
	339	cdragh_x ,cdragh_w ,cdragm_x ,cdragm_w , &
	340	kh ,kh_x ,kh_w , &
[3956]	341	wake_h, &
[3496]	342	wake_omg, d_t_wake, d_q_wake, Vprecip, qtaa, Clw, &
[4843]	343	wdtrainA, wdtrainS, wdtrainM, n2, s2, strig, zcong, zlcl_th, proba_notrig, &
[3080]	344	random_notrig, &
[5204]	345	cf_seri, d_cf_dyn, rvc_seri, d_rvc_dyn, &
	346	qsub, qissr, qcld, subfra, issrfra, gamma_cond, &
	347	dcf_sub, dcf_con, dcf_mix, &
	348	dqi_adj, dqi_sub, dqi_con, dqi_mix, &
	349	dqvc_adj, dqvc_sub, dqvc_con, dqvc_mix, &
	350	qsatliq, qsatice, &
[4059]	351	Tcontr, qcontr, qcontr2, fcontrN, fcontrP, &
[5204]	352	dcf_avi, dqi_avi, dqvc_avi, flight_dist, flight_h2o, &
[1898]	353	alp_bl_det, alp_bl_fluct_m, alp_bl_conv, &
	354	alp_bl_stat, alp_bl_fluct_tke, slab_wfbils, &
	355	weak_inversion, dthmin, cldtau, cldemi, &
[4523]	356	pmflxr, pmflxs, prfl, psfl,bsfl, re, fl, rh2m, &
[2854]	357	qsat2m, tpote, tpot, d_ts, od443aer, od550aer, dryod550aer, &
[2842]	358	od865aer, abs550aer, od550lt1aer, sconcso4, sconcno3, &
[2146]	359	sconcoa, sconcbc, sconcss, sconcdust, concso4, concno3, &
[1898]	360	concoa, concbc, concss, concdust, loadso4, &
[2854]	361	loadoa, loadbc, loadss, loaddust, loadno3, tausum_aero, drytausum_aero, &
[1898]	362	topswad_aero, topswad0_aero, solswad_aero, &
	363	solswad0_aero, topsw_aero, solsw_aero, &
	364	topsw0_aero, solsw0_aero, topswcf_aero, &
	365	solswcf_aero, topswai_aero, solswai_aero, &
[2146]	366	toplwad_aero, toplwad0_aero, sollwad_aero, &
	367	sollwad0_aero, toplwai_aero, sollwai_aero, &
[4703]	368	stratomask,&
	369	zfice, &
[1898]	370	ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
[4523]	371	ql_seri, qs_seri, qbs_seri, tr_seri, qbs_seri,&
[2284]	372	zphi, u_seri, v_seri, omega, cldfra, &
[5007]	373	rneb, rnebjn, rneblsvol, &
	374	zx_rh, zx_rhl, zx_rhi, &
	375	pfraclr, pfracld, cldfraliq, sigma2_icefracturb, mean_icefracturb, &
[4830]	376	qraindiag, qsnowdiag, dqreva, dqssub, &
[4803]	377	dqrauto,dqrcol,dqrmelt,dqrfreez, &
	378	dqsauto,dqsagg,dqsrim,dqsmelt,dqsfreez, &
	379	d_t_dyn, &
[4523]	380	d_q_dyn, d_ql_dyn, d_qs_dyn, d_qbs_dyn, &
	381	d_q_dyn2d, d_ql_dyn2d, d_qs_dyn2d, d_qbs_dyn2d, &
[1898]	382	d_u_dyn, d_v_dyn, d_t_con, d_t_ajsb, d_t_ajs, &
[1943]	383	d_u_ajs, d_v_ajs, &
[1898]	384	d_u_con, d_v_con, d_q_con, d_q_ajs, d_t_lsc, &
[2194]	385	d_t_lwr,d_t_lw0,d_t_swr,d_t_sw0, &
[1898]	386	d_t_eva, d_q_lsc, beta_prec, d_t_lscth, &
	387	d_t_lscst, d_q_lscth, d_q_lscst, plul_th, &
	388	plul_st, d_t_vdf, d_t_diss, d_q_vdf, d_q_eva, &
[5050]	389	d_t_bsss, d_q_bsss, d_qbs_bsss, d_qbs_vdf, &
[1898]	390	zw2, fraca, zmax_th, d_q_ajsb, d_t_ec, d_u_vdf, &
	391	d_v_vdf, d_u_oro, d_v_oro, d_t_oro, d_u_lif, &
[2333]	392	d_v_lif, d_t_lif, du_gwd_hines, dv_gwd_hines, d_t_hin, &
	393	dv_gwd_rando, dv_gwd_front, &
	394	east_gwstress, west_gwstress, &
[2481]	395	d_q_ch4, pmfd, pmfu, ref_liq, ref_ice, rhwriteSTD, &
[4056]	396	#ifdef ISO
	397	xtrain_lsc, xtsnow_lsc, xt_seri, xtl_seri,xts_seri,xtevap, &
	398	d_xt_dyn,d_xtl_dyn,d_xt_con,d_xt_vdf,d_xt_ajsb, &
	399	d_xt_lsc,d_xt_eva,d_xt_ch4, &
	400	d_xt_ajs, d_xt_ajsb, &
	401	d_xt_prod_nucl,d_xt_cosmo,d_xt_decroiss, &
[5109]	402	xtrunoff_diag, &
[4056]	403	#endif
[2992]	404	ep, epmax_diag, & ! epmax_cape
[4575]	405	p_tropopause, t_tropopause, z_tropopause, &
[4608]	406	zxfluxt,zxfluxq, &
	407	! offline
	408	da, mp, phi, wght_cvfd
[4703]	409	USE phys_output_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, &
	410	reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra
	411
[2704]	412	USE phys_local_var_mod, ONLY: &
[2752]	413	budg_3D_nucl, budg_3D_cond_evap, budg_3D_ocs_to_so2, budg_3D_so2_to_h2so4, &
	414	budg_sed_part, R2SO4, OCS_lifetime, SO2_lifetime, &
	415	budg_3D_backgr_ocs, budg_3D_backgr_so2, &
[3100]	416	budg_dep_dry_ocs, budg_dep_wet_ocs, &
	417	budg_dep_dry_so2, budg_dep_wet_so2, &
[2752]	418	budg_dep_dry_h2so4, budg_dep_wet_h2so4, &
	419	budg_dep_dry_part, budg_dep_wet_part, &
	420	budg_emi_ocs, budg_emi_so2, budg_emi_h2so4, budg_emi_part, &
	421	budg_ocs_to_so2, budg_so2_to_h2so4, budg_h2so4_to_part, &
[2992]	422	surf_PM25_sulf, tau_strat_550, tausum_strat, &
[4998]	423	vsed_aer, tau_strat_1020, f_r_wet, &
[5150]	424	SAD_sulfate, reff_sulfate, sulfmmr, nd_mode, sulfmmr_mode
[2690]	425
[3453]	426	USE carbon_cycle_mod, ONLY: fco2_ff, fco2_bb, fco2_land, fco2_ocean
[3857]	427	USE carbon_cycle_mod, ONLY: fco2_ocean_cor, fco2_land_cor
[3453]	428
[2704]	429	USE phys_output_var_mod, ONLY: vars_defined, snow_o, zfra_o, bils_diss, &
[2042]	430	bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
[2854]	431	itau_con, nfiles, clef_files, nid_files, dryaod_diag, &
[2333]	432	zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
[3110]	433	zustr_gwd_front, zvstr_gwd_front, sza_o, &
[2580]	434	sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
[2958]	435	cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
[2580]	436	! Marine
	437	map_prop_hc, map_prop_hist, &
	438	map_emis_hc,map_iwp_hc,map_deltaz_hc,&
	439	map_pcld_hc,map_tcld_hc,&
	440	map_emis_hist,map_iwp_hist,map_deltaz_hist,&
	441	map_rad_hist,&
	442	map_ntot,map_hc,map_hist,&
	443	map_Cb,map_ThCi,map_Anv,&
	444	map_emis_Cb,map_pcld_Cb,map_tcld_Cb,&
	445	map_emis_ThCi,map_pcld_ThCi,map_tcld_ThCi,&
	446	map_emis_Anv,map_pcld_Anv,map_tcld_Anv, &
[2989]	447	alt_tropo, &
	448	!Ionela
[4370]	449	ok_4xCO2atm, tkt, tks, taur, sss
[2580]	450
[3002]	451	USE ocean_slab_mod, ONLY: nslay, tslab, slab_bilg, tice, seaice, &
	452	slab_ekman,slab_hdiff,slab_gm,dt_ekman, dt_hdiff, dt_gm, dt_qflux
[2704]	453	USE pbl_surface_mod, ONLY: snow
	454	USE indice_sol_mod, ONLY: nbsrf
[4056]	455	#ifdef ISO
[5109]	456	USE isotopes_mod, ONLY: iso_HTO, isoName
[4056]	457	#endif
[3059]	458	USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg
[3792]	459	USE surface_data, ONLY: type_ocean, version_ocean, ok_veget, landice_opt
[2704]	460	USE aero_mod, ONLY: naero_tot, id_STRAT_phy
	461	USE ioipsl, ONLY: histend, histsync
	462	USE iophy, ONLY: set_itau_iophy, histwrite_phy
[5084]	463	USE netcdf, ONLY: nf90_fill_real
[2319]	464	USE print_control_mod, ONLY: prt_level,lunout
[1825]	465	! ug Pour les sorties XIOS
[4619]	466	USE lmdz_xios
[5310]	467	use wxios_mod, ONLY: wxios_closedef, missing_val_xios=>missing_val, wxios_set_context
[2704]	468	USE phys_cal_mod, ONLY : mth_len
[1825]	469
[2645]	470	#ifdef CPP_RRTM
	471	USE YOESW, ONLY : RSUN
	472	#endif
[5296]	473	USE compbl_mod_h
	474	USE alpale_mod
[5284]	475	USE clesphys_mod_h
[2966]	476	USE tracinca_mod, ONLY: config_inca
[4737]	477	USE config_ocean_skin_m, ONLY: activate_ocean_skin
[2114]	478
[3622]	479	USE vertical_layers_mod, ONLY: presnivs
[5252]	480	USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_STRATAER
[3622]	481
[5285]	482	USE yomcst_mod_h
[5296]	483	USE phys_constants_mod, ONLY: dobson_u
[5274]	484	IMPLICIT NONE
[1798]	485
	486
[5274]	487
[1898]	488	! Input
[4523]	489	INTEGER :: itap, ivap, iliq, isol, ibs, read_climoz
[1798]	490	INTEGER, DIMENSION(klon) :: lmax_th
	491	LOGICAL :: aerosol_couple, ok_sync
[3630]	492	LOGICAL :: ok_ade, ok_aie, ok_volcan
[1798]	493	LOGICAL, DIMENSION(klon, klev) :: ptconv, ptconvth
	494	REAL :: pdtphys
	495	CHARACTER (LEN=4), DIMENSION(nlevSTD) :: clevSTD
[1971]	496	REAL, DIMENSION(klon,nlevSTD) :: zx_tmp_fi3d_STD
[1798]	497	REAL, DIMENSION(klon) :: pphis
[2665]	498	REAL, DIMENSION(klon, klev) :: pplay, d_u, d_t
[1798]	499	REAL, DIMENSION(klon, klev+1) :: paprs
	500	REAL, DIMENSION(klon,klev,nqtot) :: qx, d_qx
[2344]	501	REAL, DIMENSION(klon, klev) :: zmasse
[2530]	502	INTEGER :: flag_aerosol_strat
[2003]	503	INTEGER :: flag_aerosol
	504	LOGICAL :: ok_cdnc
[4608]	505	REAL, DIMENSION(klon,klev) :: t ! output for phystoken - offline flux
	506	REAL, DIMENSION(klon) :: u1, v1 ! output for phystoken - offline flux
	507
[1828]	508	REAL, DIMENSION(3) :: freq_moyNMC
[1798]	509
[1898]	510	! Local
[1798]	511	INTEGER :: itau_w
[4056]	512	INTEGER :: i, iinit, iinitend=1, iff, iq, nsrf, k, ll, naero
[3691]	513	REAL, DIMENSION (klon) :: zx_tmp_fi2d, zpt_conv2d, wind100m
[1798]	514	REAL, DIMENSION (klon,klev) :: zx_tmp_fi3d, zpt_conv
	515	REAL, DIMENSION (klon,klev+1) :: zx_tmp_fi3d1
[3048]	516	REAL, DIMENSION (klon,NSW) :: zx_tmp_fi3dsp
[1798]	517	CHARACTER (LEN=4) :: bb2
[2344]	518	INTEGER, DIMENSION(nbp_lon*nbp_lat) :: ndex2d
	519	INTEGER, DIMENSION(nbp_lonnbp_latklev) :: ndex3d
[2271]	520	! REAL, PARAMETER :: missing_val=nf90_fill_real
[3474]	521	REAL, DIMENSION(klev+1,2) :: Ahyb_bounds, Bhyb_bounds
	522	REAL, DIMENSION(klev,2) :: Ahyb_mid_bounds, Bhyb_mid_bounds
[3003]	523	INTEGER :: ilev
[3622]	524	INTEGER, SAVE :: kmax_100m
	525	!$OMP THREADPRIVATE(kmax_100m)
[3691]	526	REAL :: x
[2271]	527	REAL :: missing_val
[1915]	528	REAL, PARAMETER :: un_jour=86400.
[4608]	529	CHARACTER(len=12) :: nvar
[4056]	530	INTEGER :: ISW, itr, ixt, it
[2645]	531	CHARACTER*1 ch1
[4046]	532	CHARACTER(LEN=maxlen) :: varname, dn
[4608]	533	REAL, DIMENSION(klon,klev) :: coefh_stok
	534
[2752]	535	LOGICAL, PARAMETER :: debug_strataer=.FALSE.
[4293]	536	CHARACTER(LEN=maxlen) :: unt
[5109]	537
	538	#ifdef ISO
	539	CHARACTER(LEN=maxlen) :: outiso
	540	#endif
	541
[3622]	542	REAL,DIMENSION(klon,klev) :: z, dz
	543	REAL,DIMENSION(klon) :: zrho, zt
[2704]	544
[1898]	545	! On calcul le nouveau tau:
[2551]	546	itau_w = itau_phy + itap
[1898]	547	! On le donne à iophy pour que les histwrite y aient accès:
	548	CALL set_itau_iophy(itau_w)
	549
[3435]	550	! IF (.NOT.vars_defined) THEN
[1898]	551	iinitend = 1
[3435]	552	! ELSE
	553	! iinitend = 1
	554	! ENDIF
[1898]	555
[4619]	556	IF (using_xios) CALL wxios_set_context
[3435]	557
[4619]	558	IF (using_xios) THEN
	559	missing_val=missing_val_xios
	560	ELSE
	561	missing_val=missing_val_nf90
	562	ENDIF
[3748]	563
[3622]	564	IF (.NOT.vars_defined) THEN
	565	kmax_100m=1
	566	DO k=1, klev-1
[5084]	567	IF (presnivs(k).GT.0.97*101325.) kmax_100m = k !--finding out max level for 100 m with a good margin
[3622]	568	ENDDO
	569	ENDIF
	570
[3474]	571	Ahyb_bounds(1,1) = 0.
	572	Ahyb_bounds(1,2) = aps(1)
	573	Bhyb_bounds(1,1) = 1.
	574	Bhyb_bounds(1,2) = bps(1)
[4737]	575
[3474]	576	DO ilev=2,klev
	577	Ahyb_bounds(ilev,1) = aps(ilev-1)
	578	Ahyb_bounds(ilev,2) = aps(ilev)
	579	Bhyb_bounds(ilev,1) = bps(ilev-1)
	580	Bhyb_bounds(ilev,2) = bps(ilev)
	581	ENDDO
	582
[4737]	583	Ahyb_bounds(klev+1,1) = aps(klev)
	584	Ahyb_bounds(klev+1,2) = 0.
	585	Bhyb_bounds(klev+1,1) = bps(klev)
	586	Bhyb_bounds(klev+1,2) = 0.
	587
[3474]	588	DO ilev=1, klev
	589	Ahyb_mid_bounds(ilev,1) = ap(ilev)
	590	Ahyb_mid_bounds(ilev,2) = ap(ilev+1)
	591	Bhyb_mid_bounds(ilev,1) = bp(ilev)
	592	Bhyb_mid_bounds(ilev,2) = bp(ilev+1)
[4737]	593	ENDDO
[3003]	594
[4046]	595
[1898]	596	! ug la boucle qui suit ne sert qu'une fois, pour l'initialisation, sinon il n'y a toujours qu'un seul passage:
	597	DO iinit=1, iinitend
[2780]	598	! print *,'IFF iinit=', iinit, iinitend
[4619]	599	IF (using_xios) THEN
	600	!$OMP MASTER
	601	IF (vars_defined) THEN
	602	IF (prt_level >= 10) then
	603	write(lunout,*)"phys_output_write: call xios_update_calendar, itau_w=",itau_w
	604	ENDIF
	605	! CALL xios_update_calendar(itau_w)
	606	CALL xios_update_calendar(itap)
	607	ENDIF
	608	!$OMP END MASTER
	609	!$OMP BARRIER
[2690]	610	ENDIF
[4619]	611
[1898]	612	! On procède à l'écriture ou à la définition des nombreuses variables:
[1813]	613	!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1898]	614	CALL histwrite_phy(o_phis, pphis)
[1798]	615
[2429]	616	zx_tmp_fi2d = cell_area
[2690]	617	IF (is_north_pole_phy) then
[2429]	618	zx_tmp_fi2d(1) = cell_area(1)/nbp_lon
[2690]	619	ENDIF
	620	IF (is_south_pole_phy) then
[2429]	621	zx_tmp_fi2d(klon) = cell_area(klon)/nbp_lon
[2690]	622	ENDIf
[2429]	623	CALL histwrite_phy(o_aire, zx_tmp_fi2d)
[4056]	624
[1898]	625	IF (vars_defined) THEN
	626	DO i=1, klon
	627	zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
	628	ENDDO
	629	ENDIF
[1813]	630
[1898]	631	CALL histwrite_phy(o_contfracATM, zx_tmp_fi2d)
	632	CALL histwrite_phy(o_contfracOR, pctsrf(:,is_ter))
[2645]	633	!
[4619]	634	IF (using_xios) THEN
	635
	636	CALL histwrite_phy("R_ecc",R_ecc)
	637	CALL histwrite_phy("R_peri",R_peri)
	638	CALL histwrite_phy("R_incl",R_incl)
	639	CALL histwrite_phy("solaire",solaire)
	640	CALL histwrite_phy(o_Ahyb, ap)
	641	CALL histwrite_phy(o_Bhyb, bp)
	642	CALL histwrite_phy(o_Ahyb_bounds, Ahyb_bounds)
	643	CALL histwrite_phy(o_Bhyb_bounds, Bhyb_bounds)
	644	CALL histwrite_phy(o_Ahyb_mid, aps)
	645	CALL histwrite_phy(o_Bhyb_mid, bps)
	646	CALL histwrite_phy(o_Ahyb_mid_bounds, Ahyb_mid_bounds)
	647	CALL histwrite_phy(o_Bhyb_mid_bounds, Bhyb_mid_bounds)
	648	CALL histwrite_phy(o_longitude, longitude_deg)
	649	CALL histwrite_phy(o_latitude, latitude_deg)
[2645]	650	!
	651	#ifdef CPP_RRTM
[4737]	652	IF (iflag_rrtm.EQ.1) THEN
	653	DO ISW=1, NSW
	654	WRITE(ch1,'(i1)') ISW
	655	! zx_tmp_0d=RSUN(ISW)
	656	! CALL histwrite_phy("rsun"//ch1,zx_tmp_0d)
	657	CALL histwrite_phy("rsun"//ch1,RSUN(ISW))
	658	ENDDO
	659	ENDIF
[2645]	660	#endif
	661	!
[4737]	662	CALL histwrite_phy("co2_ppm",co2_ppm)
	663	CALL histwrite_phy("CH4_ppb",CH4_ppb)
	664	CALL histwrite_phy("N2O_ppb",N2O_ppb)
	665	CALL histwrite_phy("CFC11_ppt",CFC11_ppt)
	666	CALL histwrite_phy("CFC12_ppt",CFC12_ppt)
[2645]	667	!
[4737]	668	ENDIF !using_xios
[2645]	669
[1813]	670	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2580]	671	! Simulateur AIRS
[4737]	672	IF (ok_airs) then
	673	CALL histwrite_phy(o_alt_tropo,alt_tropo)
[2580]	674
[4737]	675	CALL histwrite_phy(o_map_prop_hc,map_prop_hc)
	676	CALL histwrite_phy(o_map_prop_hist,map_prop_hist)
[2580]	677
[4737]	678	CALL histwrite_phy(o_map_emis_hc,map_emis_hc)
	679	CALL histwrite_phy(o_map_iwp_hc,map_iwp_hc)
	680	CALL histwrite_phy(o_map_deltaz_hc,map_deltaz_hc)
	681	CALL histwrite_phy(o_map_pcld_hc,map_pcld_hc)
	682	CALL histwrite_phy(o_map_tcld_hc,map_tcld_hc)
[2580]	683
[4737]	684	CALL histwrite_phy(o_map_emis_hist,map_emis_hist)
	685	CALL histwrite_phy(o_map_iwp_hist,map_iwp_hist)
	686	CALL histwrite_phy(o_map_deltaz_hist,map_deltaz_hist)
[2580]	687
[4737]	688	CALL histwrite_phy(o_map_ntot,map_ntot)
	689	CALL histwrite_phy(o_map_hc,map_hc)
	690	CALL histwrite_phy(o_map_hist,map_hist)
[2580]	691
[4737]	692	CALL histwrite_phy(o_map_Cb,map_Cb)
	693	CALL histwrite_phy(o_map_ThCi,map_ThCi)
	694	CALL histwrite_phy(o_map_Anv,map_Anv)
[2580]	695
[4737]	696	CALL histwrite_phy(o_map_emis_Cb,map_emis_Cb)
	697	CALL histwrite_phy(o_map_pcld_Cb,map_pcld_Cb)
	698	CALL histwrite_phy(o_map_tcld_Cb,map_tcld_Cb)
[2580]	699
[4737]	700	CALL histwrite_phy(o_map_emis_ThCi,map_emis_ThCi)
	701	CALL histwrite_phy(o_map_pcld_ThCi,map_pcld_ThCi)
	702	CALL histwrite_phy(o_map_tcld_ThCi,map_tcld_ThCi)
[2580]	703
[4737]	704	CALL histwrite_phy(o_map_emis_Anv,map_emis_Anv)
	705	CALL histwrite_phy(o_map_pcld_Anv,map_pcld_Anv)
	706	CALL histwrite_phy(o_map_tcld_Anv,map_tcld_Anv)
	707	ENDIF
[2580]	708
[3110]	709	CALL histwrite_phy(o_sza, sza_o)
[1898]	710	CALL histwrite_phy(o_flat, zxfluxlat)
[2385]	711	CALL histwrite_phy(o_ptstar, ptstar)
	712	CALL histwrite_phy(o_pt0, pt0)
[1898]	713	CALL histwrite_phy(o_slp, slp)
	714	CALL histwrite_phy(o_tsol, zxtsol)
[3817]	715	CALL histwrite_phy(o_t2m, zt2m)
	716	CALL histwrite_phy(o_t2m_min, zt2m)
	717	CALL histwrite_phy(o_t2m_max, zt2m)
[2146]	718	CALL histwrite_phy(o_t2m_max_mon, t2m_max_mon)
	719	CALL histwrite_phy(o_t2m_min_mon, t2m_min_mon)
[1813]	720
[1898]	721	IF (vars_defined) THEN
	722	DO i=1, klon
[4737]	723	zx_tmp_fi2d(i)=REAL(zn2mout(i,1))
[1898]	724	ENDDO
	725	ENDIF
[3817]	726	CALL histwrite_phy(o_nt2mout, zx_tmp_fi2d)
	727
	728	IF (vars_defined) THEN
	729	DO i=1, klon
[4737]	730	zx_tmp_fi2d(i)=REAL(zn2mout(i,2))
[3817]	731	ENDDO
	732	ENDIF
	733	CALL histwrite_phy(o_nt2moutfg, zx_tmp_fi2d)
	734
	735	IF (vars_defined) THEN
	736	DO i=1, klon
[4737]	737	zx_tmp_fi2d(i)=REAL(zn2mout(i,3))
[3817]	738	ENDDO
	739	ENDIF
	740	CALL histwrite_phy(o_nq2mout, zx_tmp_fi2d)
	741
	742	IF (vars_defined) THEN
	743	DO i=1, klon
[4737]	744	zx_tmp_fi2d(i)=REAL(zn2mout(i,4))
[3817]	745	ENDDO
	746	ENDIF
	747	CALL histwrite_phy(o_nq2moutfg, zx_tmp_fi2d)
	748
	749	IF (vars_defined) THEN
	750	DO i=1, klon
[4737]	751	zx_tmp_fi2d(i)=REAL(zn2mout(i,5))
[3817]	752	ENDDO
	753	ENDIF
	754	CALL histwrite_phy(o_nu2mout, zx_tmp_fi2d)
	755
	756	IF (vars_defined) THEN
	757	DO i=1, klon
[4737]	758	zx_tmp_fi2d(i)=REAL(zn2mout(i,6))
[3817]	759	ENDDO
	760	ENDIF
	761	CALL histwrite_phy(o_nu2moutfg, zx_tmp_fi2d)
	762
	763	IF (vars_defined) THEN
	764	DO i=1, klon
	765	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
	766	ENDDO
	767	ENDIF
[1898]	768	CALL histwrite_phy(o_wind10m, zx_tmp_fi2d)
[1813]	769
[1898]	770	IF (vars_defined) THEN
	771	DO i=1, klon
[3817]	772	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
[1898]	773	ENDDO
	774	ENDIF
	775	CALL histwrite_phy(o_wind10max, zx_tmp_fi2d)
[1813]	776
[2240]	777	CALL histwrite_phy(o_gusts, gustiness)
	778
[1898]	779	IF (vars_defined) THEN
[3622]	780	DO k = 1, kmax_100m !--we could stop much lower
	781	zrho(:) = pplay(:,k)/t_seri(:,k)/RD ! air density in kg/m3
	782	dz(:,k) = (paprs(:,k)-paprs(:,k+1))/zrho(:)/RG ! layer thickness in m
	783	IF (k==1) THEN
	784	z(:,1) = (paprs(:,1)-pplay(:,1))/zrho(:)/RG ! altitude middle of first layer in m
	785	zt(:) = dz(:,1) ! altitude top of first layer in m
	786	ELSE
	787	z(:,k) = zt(:) + (paprs(:,k)-pplay(:,k))/zrho(:)/RG ! altitude middle of layer k in m
	788	zt(:) = zt(:) + dz(:,k) ! altitude top of layer k in m
	789	ENDIF
	790	ENDDO
[3692]	791	wind100m(:)=missing_val
[3622]	792	DO k=1, kmax_100m-1 !--we could stop much lower
	793	DO i=1,klon
[5084]	794	IF (z(i,k).LT.100..AND.z(i,k+1).GE.100.) THEN
[3691]	795	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 + &
	796	(v_seri(i,k)+(100.-z(i,k))/(z(i,k+1)-z(i,k))(v_seri(i,k+1)-v_seri(i,k)))*2.0 )
[3622]	797	ENDIF
	798	ENDDO
	799	ENDDO
	800	ENDIF
[3691]	801	CALL histwrite_phy(o_wind100m, wind100m)
[3622]	802
	803	IF (vars_defined) THEN
[3691]	804	!--polynomial fit for 14,Vestas,1074,V136/3450 kW windmill - Olivier
	805	DO i=1,klon
[5084]	806	IF (pctsrf(i,is_ter).GT.0.05 .AND. wind100m(i).NE.missing_val) THEN
[3691]	807	x=wind100m(i)
[5084]	808	IF (x.LE.3.0 .OR. x.GE.22.5) THEN
[3691]	809	zx_tmp_fi2d(i)=0.0
[5084]	810	ELSE IF (x.GE.10.0) THEN
[3691]	811	zx_tmp_fi2d(i)=1.0
	812	ELSE
[3702]	813	zx_tmp_fi2d(i)= 10.73 + x(-14.69 + x(8.339 + x(-2.59 + x(0.4893 + x(-0.05898 + x(0.004627 + &
	814	x(-0.0002352 + x(7.478e-06 + x(-1.351e-07 + x(1.059e-09))))))))))
[3691]	815	zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
	816	ENDIF
	817	ELSE
[3692]	818	zx_tmp_fi2d(i)=missing_val
[3691]	819	ENDIF
	820	ENDDO
	821	ENDIF
	822	CALL histwrite_phy(o_loadfactor_wind_onshore, zx_tmp_fi2d)
	823
	824	IF (vars_defined) THEN
	825	!--polynomial fit for 14,Vestas,867,V164/8000 kW - Olivier
	826	DO i=1,klon
[5084]	827	IF (pctsrf(i,is_oce).GT.0.05 .AND. wind100m(i).NE.missing_val) THEN
[3691]	828	x=wind100m(i)
[5084]	829	IF (x.LE.3.0 .OR. x.GE.25.5) THEN
[3691]	830	zx_tmp_fi2d(i)=0.0
[5084]	831	ELSE IF (x.GE.12.5) THEN
[3691]	832	zx_tmp_fi2d(i)=1.0
	833	ELSE
[3702]	834	zx_tmp_fi2d(i)= 20.59 + x(-22.39 + x(10.25 + x(-2.601 + x(0.4065 + x(-0.04099 + x(0.002716 + &
	835	x(-0.0001175 + x(3.195e-06 + x(-4.959e-08 + x(3.352e-10))))))))))
[3691]	836	zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
	837	ENDIF
	838	ELSE
[3692]	839	zx_tmp_fi2d(i)=missing_val
[3691]	840	ENDIF
	841	ENDDO
	842	ENDIF
	843	CALL histwrite_phy(o_loadfactor_wind_offshore, zx_tmp_fi2d)
	844
	845	IF (vars_defined) THEN
[1898]	846	DO i = 1, klon
	847	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
	848	ENDDO
	849	ENDIF
	850	CALL histwrite_phy(o_sicf, zx_tmp_fi2d)
[3817]	851	CALL histwrite_phy(o_q2m, zq2m)
[5363]	852	!IF (vars_defined) zx_tmp_fi2d = zustar
	853	!set ustar output variable as directly related to actual surface stress
	854	IF (vars_defined) THEN
	855	zx_tmp_fi2d=0.
	856	DO nsrf=1,nbsrf
	857	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)((t_seri(:,1)RD/pplay(:,1))*0.5)(fluxu(:,1,nsrf)2+fluxv(:,1,nsrf)2)**0.25
	858	ENDDO
	859	ENDIF
[3817]	860	CALL histwrite_phy(o_ustar, zx_tmp_fi2d)
	861	CALL histwrite_phy(o_u10m, zu10m)
	862	CALL histwrite_phy(o_v10m, zv10m)
[1813]	863
[1898]	864	IF (vars_defined) THEN
	865	DO i = 1, klon
	866	zx_tmp_fi2d(i) = paprs(i,1)
	867	ENDDO
	868	ENDIF
	869	CALL histwrite_phy(o_psol, zx_tmp_fi2d)
	870	CALL histwrite_phy(o_mass, zmasse)
	871	CALL histwrite_phy(o_qsurf, zxqsurf)
[1813]	872
[1898]	873	IF (.NOT. ok_veget) THEN
	874	CALL histwrite_phy(o_qsol, qsol)
	875	ENDIF
[1813]	876
[1898]	877	IF (vars_defined) THEN
[4523]	878	IF (ok_bs) THEN
	879	DO i = 1, klon
	880	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) + bs_fall(i)
	881	ENDDO
	882	ELSE
	883	DO i = 1, klon
[1898]	884	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
[4523]	885	ENDDO
	886	ENDIF
[1898]	887	ENDIF
[1813]	888
[1898]	889	CALL histwrite_phy(o_precip, zx_tmp_fi2d)
[3095]	890	CALL histwrite_phy(o_rain_fall, rain_fall)
[4556]	891	CALL histwrite_phy(o_ndayrain, ndayrain_mth)
[1813]	892
[2481]	893	! epmax_cape:
[2483]	894	! CALL histwrite_phy(o_epmax_diag, epmax_diag)
[2481]	895	CALL histwrite_phy(o_ep, ep)
	896
[1898]	897	IF (vars_defined) THEN
	898	DO i = 1, klon
	899	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
	900	ENDDO
	901	ENDIF
	902	CALL histwrite_phy(o_plul, zx_tmp_fi2d)
[2516]	903	CALL histwrite_phy(o_plun, rain_num)
[1813]	904
[1898]	905	IF (vars_defined) THEN
	906	DO i = 1, klon
	907	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
	908	ENDDO
	909	ENDIF
	910	CALL histwrite_phy(o_pluc, zx_tmp_fi2d)
[3216]	911	CALL histwrite_phy(o_rain_con, rain_con)
[1898]	912	CALL histwrite_phy(o_snow, snow_fall)
[2209]	913	CALL histwrite_phy(o_msnow, zxsnow)
[1898]	914	CALL histwrite_phy(o_fsnow, zfra_o)
	915	CALL histwrite_phy(o_evap, evap)
[3435]	916
[4523]	917	IF (ok_bs) THEN
	918	CALL histwrite_phy(o_bsfall, bs_fall)
	919	CALL histwrite_phy(o_snowerosion, snowerosion)
	920	CALL histwrite_phy(o_ustart_lic, zxustartlic)
	921	CALL histwrite_phy(o_rhosnow_lic, zxrhoslic)
[4835]	922	CALL histwrite_phy(o_qsalt_lic, zxqsaltlic)
[4523]	923	ENDIF
	924
[3435]	925	IF (vars_defined) THEN
	926	zx_tmp_fi2d = topsw*swradcorr
	927	ENDIF
	928	CALL histwrite_phy(o_tops, zx_tmp_fi2d)
	929
	930	IF (vars_defined) THEN
	931	zx_tmp_fi2d = topsw0*swradcorr
	932	ENDIF
	933	CALL histwrite_phy(o_tops0, zx_tmp_fi2d)
	934
[1898]	935	CALL histwrite_phy(o_topl, toplw)
	936	CALL histwrite_phy(o_topl0, toplw0)
[1813]	937
[4608]	938	! offline
[4619]	939	IF (using_xios) THEN
	940	IF (offline) THEN
[4608]	941
[4619]	942	coefh_stok(:,1) = cdragh(:)
	943	coefh_stok(:,2:klev) = coefh(:,2:klev, is_ave)
[4608]	944
[4619]	945	CALL histwrite_phy('upwd_stok', upwd)
	946	CALL histwrite_phy('t_stok', t)
	947	CALL histwrite_phy('fm_th_stok', fm_therm(:,1:klev))
	948	CALL histwrite_phy('en_th_stok', entr_therm)
	949	CALL histwrite_phy('da_stok',da )
	950	CALL histwrite_phy('mp_stok',mp )
	951	CALL histwrite_phy('dnwd_stok', dnwd)
	952	CALL histwrite_phy('wght_stok', wght_cvfd)
	953	CALL histwrite_phy('coefh_stok', coefh_stok)
	954	CALL histwrite_phy('yu1_stok', u1)
	955	CALL histwrite_phy('yv1_stok', v1)
[4608]	956
[4619]	957	DO k=1,klev
	958	IF (k<10) THEN
	959	WRITE(nvar,'(i1)') k
	960	ELSE IF (k<100) THEN
	961	WRITE(nvar,'(i2)') k
	962	ELSE
	963	WRITE(nvar,'(i3)') k
[4737]	964	ENDIF
[4619]	965	nvar='phi_lev'//trim(nvar)
	966	CALL histwrite_phy(nvar,phi(:,:,k))
[4737]	967	ENDDO
[4608]	968
[4619]	969	ENDIF
[4608]	970	ENDIF
	971
[1898]	972	IF (vars_defined) THEN
[2366]	973	zx_tmp_fi2d(:) = swup(:,klevp1)*swradcorr(:)
[1898]	974	ENDIF
	975	CALL histwrite_phy(o_SWupTOA, zx_tmp_fi2d)
[1813]	976
[1898]	977	IF (vars_defined) THEN
[2366]	978	zx_tmp_fi2d(:) = swup0(:,klevp1)*swradcorr(:)
[1898]	979	ENDIF
	980	CALL histwrite_phy(o_SWupTOAclr, zx_tmp_fi2d)
[1813]	981
[4389]	982	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	983	IF (vars_defined) THEN
	984	zx_tmp_fi2d(:) = swupc0(:,klevp1)*swradcorr(:)
	985	ENDIF
	986	CALL histwrite_phy(o_SWupTOAcleanclr, zx_tmp_fi2d)
[3082]	987	ENDIF
	988
	989	IF (vars_defined) THEN
[2366]	990	zx_tmp_fi2d(:) = swdn(:,klevp1)*swradcorr(:)
[1898]	991	ENDIF
	992	CALL histwrite_phy(o_SWdnTOA, zx_tmp_fi2d)
[1813]	993
[1898]	994	IF (vars_defined) THEN
[2366]	995	zx_tmp_fi2d(:) = swdn0(:,klevp1)*swradcorr(:)
[1898]	996	ENDIF
	997	CALL histwrite_phy(o_SWdnTOAclr, zx_tmp_fi2d)
[1813]	998
[1898]	999	IF (vars_defined) THEN
[2366]	1000	zx_tmp_fi2d(:) = topsw(:)*swradcorr(:)-toplw(:)
[1898]	1001	ENDIF
	1002	CALL histwrite_phy(o_nettop, zx_tmp_fi2d)
[3435]	1003
	1004	IF (vars_defined) THEN
	1005	zx_tmp_fi2d = SWup200*swradcorr
	1006	ENDIF
	1007	CALL histwrite_phy(o_SWup200, zx_tmp_fi2d)
	1008
	1009	IF (vars_defined) THEN
	1010	zx_tmp_fi2d = SWup200clr*swradcorr
	1011	ENDIF
	1012	CALL histwrite_phy(o_SWup200clr, zx_tmp_fi2d)
	1013
	1014	IF (vars_defined) THEN
	1015	zx_tmp_fi2d = SWdn200*swradcorr
	1016	ENDIF
	1017	CALL histwrite_phy(o_SWdn200, zx_tmp_fi2d)
	1018
	1019
	1020	IF (vars_defined) THEN
	1021	zx_tmp_fi2d = SWdn200clr*swradcorr
	1022	ENDIF
	1023	CALL histwrite_phy(o_SWdn200clr, zx_tmp_fi2d)
	1024
[1898]	1025	CALL histwrite_phy(o_LWup200, LWup200)
	1026	CALL histwrite_phy(o_LWup200clr, LWup200clr)
	1027	CALL histwrite_phy(o_LWdn200, LWdn200)
	1028	CALL histwrite_phy(o_LWdn200clr, LWdn200clr)
[3435]	1029
	1030	IF (vars_defined) THEN
	1031	zx_tmp_fi2d = solsw*swradcorr
	1032	ENDIF
	1033	CALL histwrite_phy(o_sols, zx_tmp_fi2d)
	1034
	1035	IF (vars_defined) THEN
	1036	zx_tmp_fi2d = solsw0*swradcorr
	1037	ENDIF
	1038	CALL histwrite_phy(o_sols0, zx_tmp_fi2d)
[1898]	1039	CALL histwrite_phy(o_soll, sollw)
[2366]	1040	CALL histwrite_phy(o_soll0, sollw0)
[1898]	1041	CALL histwrite_phy(o_radsol, radsol)
[1813]	1042
[1898]	1043	IF (vars_defined) THEN
[2366]	1044	zx_tmp_fi2d(:) = swup(:,1)*swradcorr(:)
[1898]	1045	ENDIF
	1046	CALL histwrite_phy(o_SWupSFC, zx_tmp_fi2d)
[1813]	1047
[1898]	1048	IF (vars_defined) THEN
[2366]	1049	zx_tmp_fi2d(:) = swup0(:,1)*swradcorr(:)
[1898]	1050	ENDIF
	1051	CALL histwrite_phy(o_SWupSFCclr, zx_tmp_fi2d)
[1813]	1052
[4389]	1053	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	1054	IF (vars_defined) THEN
	1055	zx_tmp_fi2d(:) = swupc0(:,1)*swradcorr(:)
	1056	ENDIF
	1057	CALL histwrite_phy(o_SWupSFCcleanclr, zx_tmp_fi2d)
[3082]	1058	ENDIF
	1059
	1060	IF (vars_defined) THEN
[2366]	1061	zx_tmp_fi2d(:) = swdn(:,1)*swradcorr(:)
[1898]	1062	ENDIF
	1063	CALL histwrite_phy(o_SWdnSFC, zx_tmp_fi2d)
[1813]	1064
[1898]	1065	IF (vars_defined) THEN
[2366]	1066	zx_tmp_fi2d(:) = swdn0(:,1)*swradcorr(:)
[1898]	1067	ENDIF
	1068	CALL histwrite_phy(o_SWdnSFCclr, zx_tmp_fi2d)
[1813]	1069
[4389]	1070	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	1071	IF (vars_defined) THEN
	1072	zx_tmp_fi2d(:) = swdnc0(:,1)*swradcorr(:)
	1073	ENDIF
	1074	CALL histwrite_phy(o_SWdnSFCcleanclr, zx_tmp_fi2d)
[3082]	1075	ENDIF
	1076
[3756]	1077	CALL histwrite_phy(o_fdiffSWdnSFC, solswfdiff)
	1078
[3082]	1079	IF (vars_defined) THEN
[2366]	1080	zx_tmp_fi2d(:)=sollwdown(:)-sollw(:)
[1898]	1081	ENDIF
	1082	CALL histwrite_phy(o_LWupSFC, zx_tmp_fi2d)
	1083	CALL histwrite_phy(o_LWdnSFC, sollwdown)
[1813]	1084
[1898]	1085	IF (vars_defined) THEN
	1086	sollwdownclr(1:klon) = -1.*lwdn0(1:klon,1)
	1087	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
	1088	ENDIF
	1089	CALL histwrite_phy(o_LWupSFCclr, zx_tmp_fi2d)
	1090	CALL histwrite_phy(o_LWdnSFCclr, sollwdownclr)
[5204]	1091
	1092	IF (vars_defined) THEN
	1093	zx_tmp_fi2d(:) = lwup(:,klevp1)
	1094	ENDIF
	1095	CALL histwrite_phy(o_LWupTOA, zx_tmp_fi2d)
	1096
	1097	IF (vars_defined) THEN
	1098	zx_tmp_fi2d(:) = lwup0(:,klevp1)
	1099	ENDIF
	1100	CALL histwrite_phy(o_LWupTOAclr, zx_tmp_fi2d)
[3106]	1101
[4389]	1102	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	1103	IF (vars_defined) THEN
	1104	zx_tmp_fi2d(:) = lwupc0(:,klevp1)
	1105	ENDIF
	1106	CALL histwrite_phy(o_LWupTOAcleanclr, zx_tmp_fi2d)
[3106]	1107	ENDIF
[4106]	1108
[4389]	1109	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	1110	IF (vars_defined) THEN
	1111	zx_tmp_fi2d(:) = -1.*lwdnc0(:,1)
	1112	ENDIF
	1113	CALL histwrite_phy(o_LWdnSFCcleanclr, zx_tmp_fi2d)
[3106]	1114	ENDIF
	1115
[4677]	1116	!AI 08 2023 Ecrad 3Deffect
	1117	#ifdef CPP_ECRAD
[4773]	1118	CALL histwrite_phy(o_cloud_cover_sw, cloud_cover_sw)
[4737]	1119	IF (ok_3Deffect) THEN
[4773]	1120	CALL histwrite_phy(o_cloud_cover_sw_s2, cloud_cover_sw_s2)
	1121	IF (vars_defined) THEN
[4677]	1122	zx_tmp_fi2d = solsw_s2*swradcorr
	1123	ENDIF
	1124	CALL histwrite_phy(o_sols_s2, zx_tmp_fi2d)
	1125	IF (vars_defined) THEN
	1126	zx_tmp_fi2d = solsw0_s2*swradcorr
	1127	ENDIF
	1128	CALL histwrite_phy(o_sols0_s2, zx_tmp_fi2d)
	1129	CALL histwrite_phy(o_soll_s2, sollw_s2)
	1130	CALL histwrite_phy(o_soll0_s2, sollw0_s2)
	1131	IF (vars_defined) THEN
	1132	zx_tmp_fi2d = topsw_s2*swradcorr
	1133	ENDIF
	1134	CALL histwrite_phy(o_tops_s2, zx_tmp_fi2d)
	1135
	1136	IF (vars_defined) THEN
	1137	zx_tmp_fi2d = topsw0_s2*swradcorr
	1138	ENDIF
	1139	CALL histwrite_phy(o_tops0_s2, zx_tmp_fi2d)
	1140
	1141	CALL histwrite_phy(o_topl_s2, toplw_s2)
	1142	CALL histwrite_phy(o_topl0_s2, toplw0_s2)
[4773]	1143
	1144	IF (vars_defined) THEN
	1145	zx_tmp_fi2d(:) = swup_s2(:,klevp1)*swradcorr(:)
	1146	ENDIF
	1147	CALL histwrite_phy(o_SWupTOA_s2, zx_tmp_fi2d)
	1148
	1149	IF (vars_defined) THEN
	1150	zx_tmp_fi2d(:) = swup0_s2(:,klevp1)*swradcorr(:)
	1151	ENDIF
	1152	CALL histwrite_phy(o_SWupTOAclr_s2, zx_tmp_fi2d)
	1153
	1154	IF (vars_defined) THEN
	1155	zx_tmp_fi2d(:) = swdn_s2(:,klevp1)*swradcorr(:)
	1156	ENDIF
	1157	CALL histwrite_phy(o_SWdnTOA_s2, zx_tmp_fi2d)
	1158
	1159	IF (vars_defined) THEN
	1160	zx_tmp_fi2d(:) = swdn0_s2(:,klevp1)*swradcorr(:)
	1161	ENDIF
	1162	CALL histwrite_phy(o_SWdnTOAclr_s2, zx_tmp_fi2d)
	1163
	1164	IF (vars_defined) THEN
	1165	zx_tmp_fi2d(:)=sollwdown_s2(:)-sollw_s2(:)
	1166	ENDIF
	1167	CALL histwrite_phy(o_LWupSFC_s2, zx_tmp_fi2d)
	1168	CALL histwrite_phy(o_LWdnSFC_s2, sollwdown_s2)
	1169
	1170	IF (vars_defined) THEN
	1171	sollwdownclr_s2(1:klon) = -1.*lwdn0_s2(1:klon,1)
	1172	zx_tmp_fi2d(1:klon)=sollwdownclr_s2(1:klon)-sollw0_s2(1:klon)
	1173	ENDIF
	1174	CALL histwrite_phy(o_LWupSFCclr_s2, zx_tmp_fi2d)
	1175	CALL histwrite_phy(o_LWdnSFCclr_s2, sollwdownclr_s2)
	1176
	1177	IF (vars_defined) THEN
	1178	zx_tmp_fi2d(:) = swup_s2(:,1)*swradcorr(:)
	1179	ENDIF
	1180	CALL histwrite_phy(o_SWupSFC_s2, zx_tmp_fi2d)
	1181
	1182	IF (vars_defined) THEN
	1183	zx_tmp_fi2d(:) = swup0_s2(:,1)*swradcorr(:)
	1184	ENDIF
	1185	CALL histwrite_phy(o_SWupSFCclr_s2, zx_tmp_fi2d)
	1186
	1187	IF (vars_defined) THEN
	1188	zx_tmp_fi2d(:) = swdn_s2(:,1)*swradcorr(:)
	1189	ENDIF
	1190	CALL histwrite_phy(o_SWdnSFC_s2, zx_tmp_fi2d)
	1191
	1192	IF (vars_defined) THEN
	1193	zx_tmp_fi2d(:) = swdn0_s2(:,1)*swradcorr(:)
	1194	ENDIF
	1195	CALL histwrite_phy(o_SWdnSFCclr_s2, zx_tmp_fi2d)
[4808]	1196
	1197	IF (vars_defined) THEN
	1198	DO k=1, klevp1
	1199	zx_tmp_fi3d1(:,k)=swup_s2(:,k)*swradcorr(:)
	1200	ENDDO
	1201	ENDIF
	1202	CALL histwrite_phy(o_rsu_s2, zx_tmp_fi3d1)
	1203
	1204	IF (vars_defined) THEN
	1205	DO k=1, klevp1
	1206	zx_tmp_fi3d1(:,k)=swdn_s2(:,k)*swradcorr(:)
	1207	ENDDO
	1208	ENDIF
	1209	CALL histwrite_phy(o_rsd_s2, zx_tmp_fi3d1)
	1210
	1211	IF (vars_defined) THEN
	1212	DO k=1, klevp1
	1213	zx_tmp_fi3d1(:,k)=swup0_s2(:,k)*swradcorr(:)
	1214	ENDDO
	1215	ENDIF
	1216	CALL histwrite_phy(o_rsucs_s2, zx_tmp_fi3d1)
	1217
	1218	IF (vars_defined) THEN
	1219	DO k=1, klevp1
	1220	zx_tmp_fi3d1(:,k)=swdn0_s2(:,k)*swradcorr(:)
	1221	ENDDO
	1222	ENDIF
	1223	CALL histwrite_phy(o_rsdcs_s2, zx_tmp_fi3d1)
	1224
	1225	CALL histwrite_phy(o_rlu_s2, lwup_s2)
	1226	CALL histwrite_phy(o_rld_s2, lwdn_s2)
	1227	CALL histwrite_phy(o_rlucs_s2, lwup0_s2)
	1228	CALL histwrite_phy(o_rldcs_s2, lwdn0_s2)
[4773]	1229	ENDIF !ok_3Deffect
[4677]	1230	#endif
	1231
[1898]	1232	CALL histwrite_phy(o_bils, bils)
	1233	CALL histwrite_phy(o_bils_diss, bils_diss)
	1234	CALL histwrite_phy(o_bils_ec, bils_ec)
[2816]	1235	CALL histwrite_phy(o_bils_ech, bils_ech)
[1898]	1236	CALL histwrite_phy(o_bils_tke, bils_tke)
	1237	CALL histwrite_phy(o_bils_kinetic, bils_kinetic)
	1238	CALL histwrite_phy(o_bils_latent, bils_latent)
	1239	CALL histwrite_phy(o_bils_enthalp, bils_enthalp)
[1813]	1240
[1898]	1241	IF (vars_defined) THEN
	1242	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
	1243	ENDIF
	1244	CALL histwrite_phy(o_sens, zx_tmp_fi2d)
	1245	CALL histwrite_phy(o_fder, fder)
	1246	CALL histwrite_phy(o_ffonte, zxffonte)
	1247	CALL histwrite_phy(o_fqcalving, zxfqcalving)
	1248	CALL histwrite_phy(o_fqfonte, zxfqfonte)
[3095]	1249	IF (vars_defined) THEN
	1250	zx_tmp_fi2d(1:klon)=(zxfqfonte(1:klon)+rain_fall(1:klon))*pctsrf(1:klon,is_lic)
	1251	ENDIF
	1252	CALL histwrite_phy(o_mrroli, zx_tmp_fi2d)
[2517]	1253	CALL histwrite_phy(o_runofflic, zxrunofflic)
[1898]	1254	IF (vars_defined) THEN
	1255	zx_tmp_fi2d=0.
	1256	DO nsrf=1,nbsrf
	1257	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf)
	1258	ENDDO
	1259	ENDIF
	1260	CALL histwrite_phy(o_taux, zx_tmp_fi2d)
[1813]	1261
[1898]	1262	IF (vars_defined) THEN
	1263	zx_tmp_fi2d=0.
	1264	DO nsrf=1,nbsrf
	1265	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf)
	1266	ENDDO
	1267	ENDIF
	1268	CALL histwrite_phy(o_tauy, zx_tmp_fi2d)
[1813]	1269
[1898]	1270	DO nsrf = 1, nbsrf
[2670]	1271
[4737]	1272	IF (vars_defined) zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
[1898]	1273	CALL histwrite_phy(o_pourc_srf(nsrf), zx_tmp_fi2d)
[4737]	1274	IF (vars_defined) zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
[1898]	1275	CALL histwrite_phy(o_fract_srf(nsrf), zx_tmp_fi2d)
[4737]	1276	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
[1898]	1277	CALL histwrite_phy(o_taux_srf(nsrf), zx_tmp_fi2d)
[4737]	1278	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
[1898]	1279	CALL histwrite_phy(o_tauy_srf(nsrf), zx_tmp_fi2d)
[4737]	1280	IF (vars_defined) zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
[1898]	1281	CALL histwrite_phy(o_tsol_srf(nsrf), zx_tmp_fi2d)
[4737]	1282	IF (vars_defined) zx_tmp_fi2d(1 : klon) = evap_pot( 1 : klon, nsrf)
[1898]	1283	CALL histwrite_phy(o_evappot_srf(nsrf), zx_tmp_fi2d)
	1284	IF (vars_defined) zx_tmp_fi2d(1 : klon) = ustar(1 : klon, nsrf)
	1285	CALL histwrite_phy(o_ustar_srf(nsrf), zx_tmp_fi2d)
	1286	IF (vars_defined) zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf)
	1287	CALL histwrite_phy(o_u10m_srf(nsrf), zx_tmp_fi2d)
	1288	IF (vars_defined) zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf)
	1289	CALL histwrite_phy(o_v10m_srf(nsrf), zx_tmp_fi2d)
	1290	IF (vars_defined) zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf)
	1291	CALL histwrite_phy(o_t2m_srf(nsrf), zx_tmp_fi2d)
	1292	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf)
	1293	CALL histwrite_phy(o_evap_srf(nsrf), zx_tmp_fi2d)
[4737]	1294	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
[1898]	1295	CALL histwrite_phy(o_sens_srf(nsrf), zx_tmp_fi2d)
[4737]	1296	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
[1898]	1297	CALL histwrite_phy(o_lat_srf(nsrf), zx_tmp_fi2d)
[4737]	1298	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
[1898]	1299	CALL histwrite_phy(o_flw_srf(nsrf), zx_tmp_fi2d)
[4737]	1300	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
[1898]	1301	CALL histwrite_phy(o_fsw_srf(nsrf), zx_tmp_fi2d)
[4737]	1302	IF (vars_defined) zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
[1898]	1303	CALL histwrite_phy(o_wbils_srf(nsrf), zx_tmp_fi2d)
[4737]	1304	IF (vars_defined) zx_tmp_fi2d(1 : klon) = (fevap(1:klon,nsrf)-(rain_fall(1:klon)+snow_fall(1:klon)))*pctsrf(1:klon,nsrf)
[1898]	1305	CALL histwrite_phy(o_wbilo_srf(nsrf), zx_tmp_fi2d)
[4737]	1306	IF (vars_defined) zx_tmp_fi2d(1 : klon) = wfevap( 1 : klon, nsrf)
[2670]	1307	CALL histwrite_phy(o_wevap_srf(nsrf), zx_tmp_fi2d)
[4737]	1308	IF (vars_defined) zx_tmp_fi2d(1 : klon) = rain_fall(1:klon)*pctsrf(1:klon,nsrf)
[2670]	1309	CALL histwrite_phy(o_wrain_srf(nsrf), zx_tmp_fi2d)
[4737]	1310	IF (vars_defined) zx_tmp_fi2d(1 : klon) = snow_fall(1:klon)*pctsrf(1:klon,nsrf)
[2670]	1311	CALL histwrite_phy(o_wsnow_srf(nsrf), zx_tmp_fi2d)
[1865]	1312
[1898]	1313	IF (iflag_pbl > 1) THEN
	1314	CALL histwrite_phy(o_tke_srf(nsrf), pbl_tke(:,1:klev,nsrf))
[3575]	1315	!CALL histwrite_phy(o_l_mix(nsrf), l_mix(:,1:klev,nsrf))
[2561]	1316	CALL histwrite_phy(o_l_mixmin(nsrf), l_mixmin(:,1:klev,nsrf))
[1898]	1317	CALL histwrite_phy(o_tke_max_srf(nsrf), pbl_tke(:,1:klev,nsrf))
	1318	ENDIF
[2159]	1319	!jyg<
[2255]	1320	IF (iflag_pbl > 1 .AND. iflag_wake>=1 .AND. iflag_pbl_split >=1) THEN
[2159]	1321	CALL histwrite_phy(o_dltpbltke_srf(nsrf), wake_delta_pbl_TKE(:,1:klev,nsrf))
	1322	ENDIF
	1323	!>jyg
[2952]	1324	! IF (iflag_pbl > 1 .AND. ifl_pbltree >=1 ) THEN
	1325	! CALL histwrite_phy(o_treedrg_srf(nsrf), treedrg(:,1:klev,nsrf))
	1326	! ENDIF
[1813]	1327
[1898]	1328	ENDDO
[4984]	1329
	1330
[3780]	1331	IF (iflag_pbl > 1) THEN
	1332	zx_tmp_fi3d=0.
	1333	IF (vars_defined) THEN
	1334	DO nsrf=1,nbsrf
	1335	DO k=1,klev
	1336	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) &
[4881]	1337	+pctsrf(:,nsrf)*pbl_eps(:,k,nsrf)
[3780]	1338	ENDDO
	1339	ENDDO
	1340	ENDIF
	1341
[5039]	1342	CALL histwrite_phy(o_tke_dissip, zx_tmp_fi3d)
	1343
	1344	zx_tmp_fi3d=0.
	1345	IF (vars_defined) THEN
	1346	DO nsrf=1,nbsrf
	1347	DO k=1,klev
	1348	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) &
	1349	+pctsrf(:,nsrf)*tke_shear(:,k,nsrf)
	1350	ENDDO
	1351	ENDDO
	1352	ENDIF
	1353
	1354	CALL histwrite_phy(o_tke_shear, zx_tmp_fi3d)
	1355
	1356	zx_tmp_fi3d=0.
	1357	IF (vars_defined) THEN
	1358	DO nsrf=1,nbsrf
	1359	DO k=1,klev
	1360	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) &
	1361	+pctsrf(:,nsrf)*tke_buoy(:,k,nsrf)
	1362	ENDDO
	1363	ENDDO
	1364	ENDIF
	1365
	1366	CALL histwrite_phy(o_tke_buoy, zx_tmp_fi3d)
	1367
	1368
	1369	zx_tmp_fi3d=0.
	1370	IF (vars_defined) THEN
	1371	DO nsrf=1,nbsrf
	1372	DO k=1,klev
	1373	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) &
	1374	+pctsrf(:,nsrf)*tke_trans(:,k,nsrf)
	1375	ENDDO
	1376	ENDDO
	1377	ENDIF
	1378
	1379	CALL histwrite_phy(o_tke_trans, zx_tmp_fi3d)
	1380
[3780]	1381	ENDIF
[2538]	1382
	1383	IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_liq_o(1 : klon, 1)
	1384	CALL histwrite_phy(o_sens_prec_liq_oce, zx_tmp_fi2d)
	1385	IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_liq_o(1 : klon, 2)
	1386	CALL histwrite_phy(o_sens_prec_liq_sic, zx_tmp_fi2d)
	1387	IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_sol_o(1 : klon, 1)
	1388	CALL histwrite_phy(o_sens_prec_sol_oce, zx_tmp_fi2d)
	1389	IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_sol_o(1 : klon, 2)
	1390	CALL histwrite_phy(o_sens_prec_sol_sic, zx_tmp_fi2d)
	1391
	1392	IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_liq_o(1 : klon, 1)
	1393	CALL histwrite_phy(o_lat_prec_liq_oce, zx_tmp_fi2d)
	1394	IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_liq_o(1 : klon, 2)
	1395	CALL histwrite_phy(o_lat_prec_liq_sic, zx_tmp_fi2d)
	1396	IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_sol_o(1 : klon, 1)
	1397	CALL histwrite_phy(o_lat_prec_sol_oce, zx_tmp_fi2d)
	1398	IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_sol_o(1 : klon, 2)
	1399	CALL histwrite_phy(o_lat_prec_sol_sic, zx_tmp_fi2d)
	1400
[1898]	1401	DO nsrf=1,nbsrf+1
	1402	CALL histwrite_phy(o_wstar(nsrf), wstar(1 : klon, nsrf))
	1403	ENDDO
[1813]	1404
[1898]	1405	CALL histwrite_phy(o_cdrm, cdragm)
	1406	CALL histwrite_phy(o_cdrh, cdragh)
	1407	CALL histwrite_phy(o_cldl, cldl)
	1408	CALL histwrite_phy(o_cldm, cldm)
	1409	CALL histwrite_phy(o_cldh, cldh)
	1410	CALL histwrite_phy(o_cldt, cldt)
[1924]	1411	CALL histwrite_phy(o_JrNt, JrNt)
[3435]	1412
	1413	IF (vars_defined) zx_tmp_fi2d=cldl*JrNt
	1414	CALL histwrite_phy(o_cldljn, zx_tmp_fi2d)
	1415
	1416	IF (vars_defined) zx_tmp_fi2d=cldm*JrNt
	1417	CALL histwrite_phy(o_cldmjn, zx_tmp_fi2d)
	1418
	1419	IF (vars_defined) zx_tmp_fi2d=cldh*JrNt
	1420	CALL histwrite_phy(o_cldhjn, zx_tmp_fi2d)
	1421
	1422	IF (vars_defined) zx_tmp_fi2d=cldt*JrNt
	1423	CALL histwrite_phy(o_cldtjn, zx_tmp_fi2d)
	1424
[1898]	1425	CALL histwrite_phy(o_cldq, cldq)
	1426	IF (vars_defined) zx_tmp_fi2d(1:klon) = flwp(1:klon)
	1427	CALL histwrite_phy(o_lwp, zx_tmp_fi2d)
	1428	IF (vars_defined) zx_tmp_fi2d(1:klon) = fiwp(1:klon)
	1429	CALL histwrite_phy(o_iwp, zx_tmp_fi2d)
	1430	CALL histwrite_phy(o_ue, ue)
	1431	CALL histwrite_phy(o_ve, ve)
	1432	CALL histwrite_phy(o_uq, uq)
	1433	CALL histwrite_phy(o_vq, vq)
[3257]	1434	CALL histwrite_phy(o_uwat, uwat)
	1435	CALL histwrite_phy(o_vwat, vwat)
[5084]	1436	IF (iflag_con.GE.3) THEN ! sb
[1898]	1437	CALL histwrite_phy(o_cape, cape)
	1438	CALL histwrite_phy(o_pbase, ema_pcb)
	1439	CALL histwrite_phy(o_ptop, ema_pct)
	1440	CALL histwrite_phy(o_fbase, ema_cbmf)
[2690]	1441	IF (iflag_con /= 30) THEN
[1898]	1442	CALL histwrite_phy(o_plcl, plcl)
	1443	CALL histwrite_phy(o_plfc, plfc)
	1444	CALL histwrite_phy(o_wbeff, wbeff)
[2824]	1445	CALL histwrite_phy(o_convoccur, convoccur)
[2690]	1446	ENDIF
[1816]	1447
[1898]	1448	CALL histwrite_phy(o_cape_max, cape)
[1813]	1449
[1898]	1450	CALL histwrite_phy(o_upwd, upwd)
	1451	CALL histwrite_phy(o_Ma, Ma)
	1452	CALL histwrite_phy(o_dnwd, dnwd)
	1453	CALL histwrite_phy(o_dnwd0, dnwd0)
[3171]	1454	!! The part relative to the frequency of occurence of convection
	1455	!! is now grouped with the part relative to thermals and shallow
	1456	!! convection (output of the 3 fields: ftime_deepcv, ftime_th and
	1457	!! ftime_con).
[3167]	1458	IF (vars_defined) THEN
[2690]	1459	IF (iflag_thermals>=1)THEN
[4109]	1460	zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm(:,1:klev)
[1898]	1461	ELSE
[4109]	1462	zx_tmp_fi3d=dnwd+dnwd0+upwd
[1898]	1463	ENDIF
	1464	ENDIF
	1465	CALL histwrite_phy(o_mc, zx_tmp_fi3d)
	1466	ENDIF !iflag_con .GE. 3
	1467	CALL histwrite_phy(o_prw, prw)
[2496]	1468	CALL histwrite_phy(o_prlw, prlw)
	1469	CALL histwrite_phy(o_prsw, prsw)
[4523]	1470	IF (ok_bs) THEN
	1471	CALL histwrite_phy(o_prbsw, prbsw)
	1472	ENDIF
[4887]	1473	CALL histwrite_phy(o_water_budget, water_budget)
[1898]	1474	CALL histwrite_phy(o_s_pblh, s_pblh)
	1475	CALL histwrite_phy(o_s_pblt, s_pblt)
	1476	CALL histwrite_phy(o_s_lcl, s_lcl)
	1477	CALL histwrite_phy(o_s_therm, s_therm)
	1478	!IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
	1479	! IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN
	1480	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1481	! $o_s_capCL%name,itau_w,s_capCL)
	1482	! ENDIF
	1483	! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
	1484	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1485	! $o_s_oliqCL%name,itau_w,s_oliqCL)
	1486	! ENDIF
	1487	! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
	1488	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1489	! $o_s_cteiCL%name,itau_w,s_cteiCL)
	1490	! ENDIF
	1491	! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
	1492	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1493	! $o_s_trmb1%name,itau_w,s_trmb1)
	1494	! ENDIF
	1495	! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
	1496	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1497	! $o_s_trmb2%name,itau_w,s_trmb2)
	1498	! ENDIF
	1499	! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
	1500	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1501	! $o_s_trmb3%name,itau_w,s_trmb3)
	1502	! ENDIF
[1813]	1503
[4619]	1504	IF (.NOT. using_xios) THEN
	1505	IF (.NOT.ok_all_xml) THEN
	1506	! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX:
	1507	! Champs interpolles sur des niveaux de pression
	1508	DO iff=1, nfiles
	1509	ll=0
	1510	DO k=1, nlevSTD
[5267]	1511	bb2=clevSTD(k)
[4619]	1512	IF (bb2.EQ."850".OR.bb2.EQ."700".OR. &
	1513	bb2.EQ."500".OR.bb2.EQ."200".OR. &
	1514	bb2.EQ."100".OR. &
	1515	bb2.EQ."50".OR.bb2.EQ."10") THEN
	1516
	1517	! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1518	ll=ll+1
	1519	CALL histwrite_phy(o_uSTDlevs(ll),uwriteSTD(:,k,iff), iff)
	1520	CALL histwrite_phy(o_vSTDlevs(ll),vwriteSTD(:,k,iff), iff)
	1521	CALL histwrite_phy(o_wSTDlevs(ll),wwriteSTD(:,k,iff), iff)
	1522	CALL histwrite_phy(o_zSTDlevs(ll),phiwriteSTD(:,k,iff), iff)
	1523	CALL histwrite_phy(o_qSTDlevs(ll),qwriteSTD(:,k,iff), iff)
	1524	CALL histwrite_phy(o_tSTDlevs(ll),twriteSTD(:,k,iff), iff)
	1525
	1526	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
	1527	ENDDO
	1528	ENDDO
	1529	ENDIF
	1530	ENDIF
	1531
[5267]	1532
[4619]	1533	IF (using_xios) THEN
	1534	IF (ok_all_xml) THEN
	1535	!XIOS CALL xios_get_field_attr("u850",default_value=missing_val)
	1536	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1537	ll=0
	1538	DO k=1, nlevSTD
	1539	bb2=clevSTD(k)
[5084]	1540	IF (bb2.EQ."850".OR.bb2.EQ."700".OR. &
	1541	bb2.EQ."500".OR.bb2.EQ."200".OR. &
	1542	bb2.EQ."100".OR. &
	1543	bb2.EQ."50".OR.bb2.EQ."10") THEN
[4619]	1544	ll=ll+1
	1545	CALL histwrite_phy(o_uSTDlevs(ll),ulevSTD(:,k))
	1546	CALL histwrite_phy(o_vSTDlevs(ll),vlevSTD(:,k))
	1547	CALL histwrite_phy(o_wSTDlevs(ll),wlevSTD(:,k))
	1548	CALL histwrite_phy(o_zSTDlevs(ll),philevSTD(:,k))
	1549	CALL histwrite_phy(o_qSTDlevs(ll),qlevSTD(:,k))
	1550	CALL histwrite_phy(o_tSTDlevs(ll),tlevSTD(:,k))
	1551	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
	1552	ENDDO
	1553	ENDIF
	1554	ENDIF
[1813]	1555
[1898]	1556	IF (vars_defined) THEN
	1557	DO i=1, klon
[5084]	1558	IF (pctsrf(i,is_oce).GT.epsfra.OR. &
	1559	pctsrf(i,is_sic).GT.epsfra) THEN
[1898]	1560	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ &
	1561	ftsol(i, is_sic) * pctsrf(i,is_sic))/ &
	1562	(pctsrf(i,is_oce)+pctsrf(i,is_sic))
	1563	ELSE
	1564	zx_tmp_fi2d(i) = 273.15
	1565	ENDIF
	1566	ENDDO
[1813]	1567	ENDIF
[1898]	1568	CALL histwrite_phy(o_t_oce_sic, zx_tmp_fi2d)
[1813]	1569
[1898]	1570	! Couplage convection-couche limite
[5084]	1571	IF (iflag_con.GE.3) THEN
[1898]	1572	IF (iflag_coupl>=1) THEN
	1573	CALL histwrite_phy(o_ale_bl, ale_bl)
	1574	CALL histwrite_phy(o_alp_bl, alp_bl)
	1575	ENDIF !iflag_coupl>=1
	1576	ENDIF !(iflag_con.GE.3)
	1577	! Wakes
[5084]	1578	IF (iflag_con.EQ.3) THEN
[3523]	1579	CALL histwrite_phy(o_Mipsh, Mipsh)
[1898]	1580	IF (iflag_wake>=1) THEN
	1581	CALL histwrite_phy(o_ale_wk, ale_wake)
	1582	CALL histwrite_phy(o_alp_wk, alp_wake)
[2393]	1583	IF (iflag_pbl_split>=1) THEN
[3179]	1584	!! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_x(1:klon,1:klev)/pdtphys
	1585	!! CALL histwrite_phy(o_dtvdf_x ,zx_tmp_fi3d)
	1586	!! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_w(1:klon,1:klev)/pdtphys
	1587	!! CALL histwrite_phy(o_dtvdf_w ,zx_tmp_fi3d)
	1588	!! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_x(1:klon,1:klev)/pdtphys
	1589	!! CALL histwrite_phy(o_dqvdf_x ,zx_tmp_fi3d)
	1590	!! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_w(1:klon,1:klev)/pdtphys
	1591	!
	1592	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf_x(1:klon,1:klev)/pdtphys
[2393]	1593	CALL histwrite_phy(o_dtvdf_x ,zx_tmp_fi3d)
[3179]	1594	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf_w(1:klon,1:klev)/pdtphys
[2393]	1595	CALL histwrite_phy(o_dtvdf_w ,zx_tmp_fi3d)
[3179]	1596	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf_x(1:klon,1:klev)/pdtphys
[2393]	1597	CALL histwrite_phy(o_dqvdf_x ,zx_tmp_fi3d)
[3179]	1598	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf_w(1:klon,1:klev)/pdtphys
	1599	!
[2393]	1600	CALL histwrite_phy(o_dqvdf_w ,zx_tmp_fi3d)
[3888]	1601	IF (vars_defined) zx_tmp_fi2d(1:klon)=-1*sens_x(1:klon)
	1602	CALL histwrite_phy(o_sens_x ,zx_tmp_fi2d)
	1603	IF (vars_defined) zx_tmp_fi2d(1:klon)=-1*sens_w(1:klon)
	1604	CALL histwrite_phy(o_sens_w ,zx_tmp_fi2d)
[2393]	1605	CALL histwrite_phy(o_flat_x ,zxfluxlat_x)
	1606	CALL histwrite_phy(o_flat_w ,zxfluxlat_w)
[3888]	1607	zx_tmp_fi2d=0.
	1608	IF (vars_defined) THEN
	1609	DO nsrf=1,nbsrf
	1610	zx_tmp_fi2d(:)=zx_tmp_fi2d(:) &
	1611	+pctsrf(:,nsrf)*delta_tsurf(:,nsrf)
	1612	ENDDO
	1613	ENDIF
	1614	CALL histwrite_phy(o_delta_tsurf,zx_tmp_fi2d)
[2393]	1615	CALL histwrite_phy(o_cdragh_x ,cdragh_x )
	1616	CALL histwrite_phy(o_cdragh_w ,cdragh_w )
	1617	CALL histwrite_phy(o_cdragm_x ,cdragm_x )
	1618	CALL histwrite_phy(o_cdragm_w ,cdragm_w )
	1619	CALL histwrite_phy(o_kh ,kh )
	1620	CALL histwrite_phy(o_kh_x ,kh_x )
[4843]	1621	CALL histwrite_phy(o_strig, strig)
	1622	CALL histwrite_phy(o_zcong, zcong)
	1623	CALL histwrite_phy(o_zlcl_th, zlcl_th)
[2393]	1624	CALL histwrite_phy(o_kh_w ,kh_w )
	1625	ENDIF ! (iflag_pbl_split>=1)
[1898]	1626	CALL histwrite_phy(o_ale, ale)
	1627	CALL histwrite_phy(o_alp, alp)
	1628	CALL histwrite_phy(o_cin, cin)
	1629	CALL histwrite_phy(o_WAPE, wake_pe)
[3208]	1630	CALL histwrite_phy(o_cv_gen, cv_gen)
[1898]	1631	CALL histwrite_phy(o_wake_h, wake_h)
[3208]	1632	CALL histwrite_phy(o_wake_dens, wake_dens)
[1898]	1633	CALL histwrite_phy(o_wake_s, wake_s)
	1634	CALL histwrite_phy(o_wake_deltat, wake_deltat)
	1635	CALL histwrite_phy(o_wake_deltaq, wake_deltaq)
	1636	CALL histwrite_phy(o_wake_omg, wake_omg)
	1637	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev) &
	1638	/pdtphys
	1639	CALL histwrite_phy(o_dtwak, zx_tmp_fi3d)
	1640	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
	1641	CALL histwrite_phy(o_dqwak, zx_tmp_fi3d)
[3435]	1642	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
[2486]	1643	CALL histwrite_phy(o_dqwak2d, zx_tmp_fi2d)
[1898]	1644	ENDIF ! iflag_wake>=1
	1645	CALL histwrite_phy(o_ftd, ftd)
	1646	CALL histwrite_phy(o_fqd, fqd)
[2253]	1647	ENDIF !(iflag_con.EQ.3)
[5084]	1648	IF (iflag_con.EQ.3.OR.iflag_con.EQ.30) THEN
[1898]	1649	! sortie RomP convection descente insaturee iflag_con=30
[2253]	1650	! etendue a iflag_con=3 (jyg)
[1898]	1651	CALL histwrite_phy(o_Vprecip, Vprecip)
[3496]	1652	CALL histwrite_phy(o_qtaa, qtaa)
	1653	CALL histwrite_phy(o_clwaa, clw)
[1898]	1654	CALL histwrite_phy(o_wdtrainA, wdtrainA)
[3496]	1655	CALL histwrite_phy(o_wdtrainS, wdtrainS)
[1898]	1656	CALL histwrite_phy(o_wdtrainM, wdtrainM)
	1657	ENDIF !(iflag_con.EQ.3.or.iflag_con.EQ.30)
[1813]	1658	!!! nrlmd le 10/04/2012
[1898]	1659	IF (iflag_trig_bl>=1) THEN
	1660	CALL histwrite_phy(o_n2, n2)
	1661	CALL histwrite_phy(o_s2, s2)
	1662	CALL histwrite_phy(o_proba_notrig, proba_notrig)
	1663	CALL histwrite_phy(o_random_notrig, random_notrig)
	1664	CALL histwrite_phy(o_ale_bl_stat, ale_bl_stat)
	1665	CALL histwrite_phy(o_ale_bl_trig, ale_bl_trig)
[1813]	1666	ENDIF !(iflag_trig_bl>=1)
[1898]	1667	IF (iflag_clos_bl>=1) THEN
	1668	CALL histwrite_phy(o_alp_bl_det, alp_bl_det)
	1669	CALL histwrite_phy(o_alp_bl_fluct_m, alp_bl_fluct_m)
	1670	CALL histwrite_phy(o_alp_bl_fluct_tke, &
	1671	alp_bl_fluct_tke)
	1672	CALL histwrite_phy(o_alp_bl_conv, alp_bl_conv)
	1673	CALL histwrite_phy(o_alp_bl_stat, alp_bl_stat)
[1813]	1674	ENDIF !(iflag_clos_bl>=1)
	1675	!!! fin nrlmd le 10/04/2012
[2057]	1676	! Output of slab ocean variables
[1898]	1677	IF (type_ocean=='slab ') THEN
[3002]	1678	CALL histwrite_phy(o_slab_bils, slab_wfbils)
[5084]	1679	IF (nslay.EQ.1) THEN
[3950]	1680	IF (vars_defined) zx_tmp_fi2d(:)=tslab(:,1)
[2057]	1681	CALL histwrite_phy(o_tslab, zx_tmp_fi2d)
[3950]	1682	IF (vars_defined) zx_tmp_fi2d(:)=dt_qflux(:,1)
[3002]	1683	CALL histwrite_phy(o_slab_qflux, zx_tmp_fi2d)
[2057]	1684	ELSE
[2656]	1685	CALL histwrite_phy(o_tslab, tslab(:,1:nslay))
[3002]	1686	CALL histwrite_phy(o_slab_qflux, dt_qflux(:,1:nslay))
[2690]	1687	ENDIF
[2209]	1688	IF (version_ocean=='sicINT') THEN
	1689	CALL histwrite_phy(o_slab_bilg, slab_bilg)
	1690	CALL histwrite_phy(o_slab_tice, tice)
	1691	CALL histwrite_phy(o_slab_sic, seaice)
[2690]	1692	ENDIF
[3002]	1693	IF (slab_gm) THEN
	1694	CALL histwrite_phy(o_slab_gm, dt_gm(:,1:nslay))
[3453]	1695	ENDIF
[2656]	1696	IF (slab_hdiff) THEN
[5084]	1697	IF (nslay.EQ.1) THEN
[3950]	1698	IF (vars_defined) zx_tmp_fi2d(:)=dt_hdiff(:,1)
[2656]	1699	CALL histwrite_phy(o_slab_hdiff, zx_tmp_fi2d)
	1700	ELSE
	1701	CALL histwrite_phy(o_slab_hdiff, dt_hdiff(:,1:nslay))
[2690]	1702	ENDIF
	1703	ENDIF
[5084]	1704	IF (slab_ekman.GT.0) THEN
	1705	IF (nslay.EQ.1) THEN
[3950]	1706	IF (vars_defined) zx_tmp_fi2d(:)=dt_ekman(:,1)
[2656]	1707	CALL histwrite_phy(o_slab_ekman, zx_tmp_fi2d)
	1708	ELSE
	1709	CALL histwrite_phy(o_slab_ekman, dt_ekman(:,1:nslay))
[2690]	1710	ENDIF
	1711	ENDIF
[1898]	1712	ENDIF !type_ocean == force/slab
	1713	CALL histwrite_phy(o_weakinv, weak_inversion)
	1714	CALL histwrite_phy(o_dthmin, dthmin)
	1715	CALL histwrite_phy(o_cldtau, cldtau)
	1716	CALL histwrite_phy(o_cldemi, cldemi)
	1717	CALL histwrite_phy(o_pr_con_l, pmflxr(:,1:klev))
	1718	CALL histwrite_phy(o_pr_con_i, pmflxs(:,1:klev))
	1719	CALL histwrite_phy(o_pr_lsc_l, prfl(:,1:klev))
	1720	CALL histwrite_phy(o_pr_lsc_i, psfl(:,1:klev))
	1721	CALL histwrite_phy(o_re, re)
	1722	CALL histwrite_phy(o_fl, fl)
[4523]	1723
	1724	IF (ok_bs) THEN
	1725	CALL histwrite_phy(o_pr_bs, bsfl(:,1:klev))
	1726	ENDIF
	1727
[1898]	1728	IF (vars_defined) THEN
	1729	DO i=1, klon
[5084]	1730	IF (zt2m(i).LE.273.15) then
[3929]	1731	zx_tmp_fi2d(i)=MAX(0.,rh2m(i)*100.)
	1732	ELSE
	1733	zx_tmp_fi2d(i)=MAX(0.,MIN(100.,rh2m(i)*100.))
	1734	ENDIF
[1898]	1735	ENDDO
	1736	ENDIF
	1737	CALL histwrite_phy(o_rh2m, zx_tmp_fi2d)
[1813]	1738
[3165]	1739	! IF (vars_defined) THEN
	1740	! DO i=1, klon
	1741	! zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
	1742	! ENDDO
	1743	! ENDIF
	1744	! CALL histwrite_phy(o_rh2m_min, zx_tmp_fi2d)
[1813]	1745
[3165]	1746	! IF (vars_defined) THEN
	1747	! DO i=1, klon
	1748	! zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
	1749	! ENDDO
	1750	! ENDIF
	1751	! CALL histwrite_phy(o_rh2m_max, zx_tmp_fi2d)
[1813]	1752
[3817]	1753	CALL histwrite_phy(o_qsat2m, qsat2m)
[1898]	1754	CALL histwrite_phy(o_tpot, tpot)
	1755	CALL histwrite_phy(o_tpote, tpote)
	1756	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
	1757	CALL histwrite_phy(o_SWnetOR, zx_tmp_fi2d)
	1758	CALL histwrite_phy(o_LWdownOR, sollwdown)
	1759	CALL histwrite_phy(o_snowl, snow_lsc)
	1760	CALL histwrite_phy(o_solldown, sollwdown)
	1761	CALL histwrite_phy(o_dtsvdfo, d_ts(:,is_oce))
	1762	CALL histwrite_phy(o_dtsvdft, d_ts(:,is_ter))
	1763	CALL histwrite_phy(o_dtsvdfg, d_ts(:,is_lic))
	1764	CALL histwrite_phy(o_dtsvdfi, d_ts(:,is_sic))
[2243]	1765	CALL histwrite_phy(o_z0m, z0m(:,nbsrf+1))
	1766	CALL histwrite_phy(o_z0h, z0h(:,nbsrf+1))
[2952]	1767
[2854]	1768	! od550 per species
[2146]	1769	!--OLIVIER
	1770	!This is warranted by treating INCA aerosols as offline aerosols
[4489]	1771	#ifndef CPP_ECRAD
[5084]	1772	IF (flag_aerosol.GT.0) THEN
[4389]	1773	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	1774
	1775	CALL histwrite_phy(o_od443aer, od443aer)
	1776	CALL histwrite_phy(o_od550aer, od550aer)
	1777	CALL histwrite_phy(o_od865aer, od865aer)
	1778	CALL histwrite_phy(o_abs550aer, abs550aer)
	1779	CALL histwrite_phy(o_od550lt1aer, od550lt1aer)
	1780	CALL histwrite_phy(o_sconcso4, sconcso4)
	1781	CALL histwrite_phy(o_sconcno3, sconcno3)
	1782	CALL histwrite_phy(o_sconcoa, sconcoa)
	1783	CALL histwrite_phy(o_sconcbc, sconcbc)
	1784	CALL histwrite_phy(o_sconcss, sconcss)
	1785	CALL histwrite_phy(o_sconcdust, sconcdust)
	1786	CALL histwrite_phy(o_concso4, concso4)
	1787	CALL histwrite_phy(o_concno3, concno3)
	1788	CALL histwrite_phy(o_concoa, concoa)
	1789	CALL histwrite_phy(o_concbc, concbc)
	1790	CALL histwrite_phy(o_concss, concss)
	1791	CALL histwrite_phy(o_concdust, concdust)
	1792	CALL histwrite_phy(o_loadso4, loadso4)
	1793	CALL histwrite_phy(o_loadoa, loadoa)
	1794	CALL histwrite_phy(o_loadbc, loadbc)
	1795	CALL histwrite_phy(o_loadss, loadss)
	1796	CALL histwrite_phy(o_loaddust, loaddust)
	1797	CALL histwrite_phy(o_loadno3, loadno3)
	1798	CALL histwrite_phy(o_dryod550aer, dryod550aer)
	1799	DO naero = 1, naero_tot-1
	1800	CALL histwrite_phy(o_drytausumaero(naero),drytausum_aero(:,naero))
[4737]	1801	ENDDO
[4106]	1802	ENDIF
[1898]	1803	ENDIF
[3630]	1804	!--STRAT AER
[5084]	1805	IF (flag_aerosol.GT.0.OR.flag_aerosol_strat.GT.0) THEN
[3630]	1806	DO naero = 1, naero_tot
	1807	CALL histwrite_phy(o_tausumaero(naero),tausum_aero(:,2,naero))
[4737]	1808	ENDDO
[3630]	1809	ENDIF
[5084]	1810	IF (flag_aerosol_strat.GT.0) THEN
[3630]	1811	CALL histwrite_phy(o_tausumaero_lw,tausum_aero(:,6,id_STRAT_phy))
	1812	ENDIF
[2992]	1813
	1814	CALL histwrite_phy(o_p_tropopause, p_tropopause)
	1815	CALL histwrite_phy(o_t_tropopause, t_tropopause)
	1816	CALL histwrite_phy(o_z_tropopause, z_tropopause)
	1817
[2993]	1818	! ThL -- In the following, we assume read_climoz == 1
[3435]	1819	IF (vars_defined) THEN
	1820	zx_tmp_fi2d = 0.0 ! Computation for strato, added ThL
	1821	DO k=1, klev
	1822	zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * stratomask(:,k) * 1.e3
[4737]	1823	ENDDO
[3435]	1824	ENDIF
[2993]	1825	CALL histwrite_phy(o_col_O3_strato, zx_tmp_fi2d) ! Added ThL
[3435]	1826
	1827	IF (vars_defined) THEN
	1828	zx_tmp_fi2d = 0.0 ! Computation for tropo, added ThL
	1829	DO k=1, klev
	1830	zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * (1.0-stratomask(:,k)) * 1.e3
[4737]	1831	ENDDO
[3435]	1832	ENDIF
[2993]	1833	CALL histwrite_phy(o_col_O3_tropo, zx_tmp_fi2d) ! Added ThL
	1834	! end add ThL
	1835
[5252]	1836	IF (CPPKEY_STRATAER) THEN
[4389]	1837	IF (type_trac=='coag') THEN
[2690]	1838	CALL histwrite_phy(o_R2SO4, R2SO4)
	1839	CALL histwrite_phy(o_OCS_lifetime, OCS_lifetime)
	1840	CALL histwrite_phy(o_SO2_lifetime, SO2_lifetime)
[2752]	1841	CALL histwrite_phy(o_budg_3D_backgr_ocs, budg_3D_backgr_ocs)
	1842	CALL histwrite_phy(o_budg_3D_backgr_so2, budg_3D_backgr_so2)
	1843	CALL histwrite_phy(o_budg_3D_ocs_to_so2, budg_3D_ocs_to_so2)
	1844	CALL histwrite_phy(o_budg_3D_so2_to_h2so4, budg_3D_so2_to_h2so4)
	1845	CALL histwrite_phy(o_budg_3D_nucl, budg_3D_nucl)
	1846	CALL histwrite_phy(o_budg_3D_cond_evap, budg_3D_cond_evap)
[3100]	1847	CALL histwrite_phy(o_budg_dep_dry_ocs, budg_dep_dry_ocs)
	1848	CALL histwrite_phy(o_budg_dep_wet_ocs, budg_dep_wet_ocs)
	1849	CALL histwrite_phy(o_budg_dep_dry_so2, budg_dep_dry_so2)
	1850	CALL histwrite_phy(o_budg_dep_wet_so2, budg_dep_wet_so2)
[2752]	1851	CALL histwrite_phy(o_budg_dep_dry_h2so4, budg_dep_dry_h2so4)
	1852	CALL histwrite_phy(o_budg_dep_wet_h2so4, budg_dep_wet_h2so4)
	1853	CALL histwrite_phy(o_budg_dep_dry_part, budg_dep_dry_part)
	1854	CALL histwrite_phy(o_budg_dep_wet_part, budg_dep_wet_part)
	1855	CALL histwrite_phy(o_budg_emi_ocs, budg_emi_ocs)
	1856	CALL histwrite_phy(o_budg_emi_so2, budg_emi_so2)
	1857	CALL histwrite_phy(o_budg_emi_h2so4, budg_emi_h2so4)
	1858	CALL histwrite_phy(o_budg_emi_part, budg_emi_part)
	1859	CALL histwrite_phy(o_budg_ocs_to_so2, budg_ocs_to_so2)
	1860	CALL histwrite_phy(o_budg_so2_to_h2so4, budg_so2_to_h2so4)
	1861	CALL histwrite_phy(o_budg_h2so4_to_part, budg_h2so4_to_part)
	1862	CALL histwrite_phy(o_budg_sed_part, budg_sed_part)
[2690]	1863	CALL histwrite_phy(o_surf_PM25_sulf, surf_PM25_sulf)
	1864	CALL histwrite_phy(o_vsed_aer, vsed_aer)
	1865	CALL histwrite_phy(o_f_r_wet, f_r_wet)
	1866	CALL histwrite_phy(o_ext_strat_550, tau_strat_550)
	1867	CALL histwrite_phy(o_ext_strat_1020, tau_strat_1020)
	1868	CALL histwrite_phy(o_tau_strat_550, tausum_strat(:,1))
	1869	CALL histwrite_phy(o_tau_strat_1020, tausum_strat(:,2))
[4998]	1870	CALL histwrite_phy(o_SAD_sulfate, SAD_sulfate)
[5150]	1871	CALL histwrite_phy(o_reff_sulfate, reff_sulfate)
[4998]	1872	CALL histwrite_phy(o_sulfmmr, sulfmmr)
	1873	! All BINs fields
	1874	DO itr = 1, nbtr_bin
	1875	CALL histwrite_phy(o_nd_mode(itr), nd_mode(:,:,itr))
	1876	CALL histwrite_phy(o_sulfmmr_mode(itr), sulfmmr_mode(:,:,itr))
	1877	ENDDO !--itr
[2690]	1878	ENDIF
[5252]	1879	END IF
[3479]	1880	!NL
	1881	IF (ok_volcan .AND. ok_ade) THEN
	1882	DO k=1, klev
	1883	IF (vars_defined) zx_tmp_fi3d(:,k)=heat_volc(:,k)*swradcorr(:)
	1884	ENDDO
	1885	CALL histwrite_phy(o_heat_volc, zx_tmp_fi3d)
	1886	DO k=1, klev
[3480]	1887	IF (vars_defined) zx_tmp_fi3d(:,k)=cool_volc(:,k)
[3479]	1888	ENDDO
	1889	CALL histwrite_phy(o_cool_volc, zx_tmp_fi3d)
	1890	ENDIF
[1813]	1891	IF (ok_ade) THEN
[3435]	1892	IF (vars_defined) zx_tmp_fi2d(:)=topswad_aero*swradcorr
	1893	CALL histwrite_phy(o_topswad, zx_tmp_fi2d)
	1894
	1895	IF (vars_defined) zx_tmp_fi2d(:)=topswad0_aero*swradcorr
	1896	CALL histwrite_phy(o_topswad0, zx_tmp_fi2d)
	1897
	1898	IF (vars_defined) zx_tmp_fi2d(:)=solswad_aero*swradcorr
	1899	CALL histwrite_phy(o_solswad, zx_tmp_fi2d)
	1900
	1901	IF (vars_defined) zx_tmp_fi2d(:)=solswad0_aero*swradcorr
	1902	CALL histwrite_phy(o_solswad0, zx_tmp_fi2d)
	1903
[4389]	1904	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	1905
	1906	CALL histwrite_phy(o_toplwad, toplwad_aero)
	1907	CALL histwrite_phy(o_toplwad0, toplwad0_aero)
	1908	CALL histwrite_phy(o_sollwad, sollwad_aero)
	1909	CALL histwrite_phy(o_sollwad0, sollwad0_aero)
	1910	ENDIF
[1898]	1911	!====MS forcing diagnostics
[3435]	1912	!ym warning : topsw_aero, solsw_aero, topsw0_aero, solsw0_aero are not defined by model
	1913	!ym => init to 0 in radlwsw_m.F90 ztopsw_aero, zsolsw_aero, ztopsw0_aero, zsolsw0_aero
	1914
[3630]	1915	IF (vars_defined) zx_tmp_fi2d(:)=topsw_aero(:,1)*swradcorr(:)
	1916	CALL histwrite_phy(o_swtoaas_nat,zx_tmp_fi2d)
	1917	IF (vars_defined) zx_tmp_fi2d(:)=solsw_aero(:,1)*swradcorr(:)
	1918	CALL histwrite_phy(o_swsrfas_nat,zx_tmp_fi2d)
	1919	IF (vars_defined) zx_tmp_fi2d(:)=topsw0_aero(:,1)*swradcorr(:)
	1920	CALL histwrite_phy(o_swtoacs_nat,zx_tmp_fi2d)
	1921	IF (vars_defined) zx_tmp_fi2d(:)=solsw0_aero(:,1)*swradcorr(:)
	1922	CALL histwrite_phy(o_swsrfcs_nat,zx_tmp_fi2d)
	1923	!ant
	1924	IF (vars_defined) zx_tmp_fi2d(:)=topsw_aero(:,2)*swradcorr(:)
	1925	CALL histwrite_phy(o_swtoaas_ant,zx_tmp_fi2d)
	1926	IF (vars_defined) zx_tmp_fi2d(:)=solsw_aero(:,2)*swradcorr(:)
	1927	CALL histwrite_phy(o_swsrfas_ant,zx_tmp_fi2d)
	1928	IF (vars_defined) zx_tmp_fi2d(:)=topsw0_aero(:,2)*swradcorr(:)
	1929	CALL histwrite_phy(o_swtoacs_ant,zx_tmp_fi2d)
	1930	IF (vars_defined) zx_tmp_fi2d(:)=solsw0_aero(:,2)*swradcorr(:)
	1931	CALL histwrite_phy(o_swsrfcs_ant,zx_tmp_fi2d)
	1932	!cf
	1933	IF (.not. aerosol_couple) THEN
	1934	IF (vars_defined) zx_tmp_fi2d(:)=topswcf_aero(:,1)*swradcorr(:)
	1935	CALL histwrite_phy(o_swtoacf_nat,zx_tmp_fi2d)
	1936	IF (vars_defined) zx_tmp_fi2d(:)=solswcf_aero(:,1)*swradcorr(:)
	1937	CALL histwrite_phy(o_swsrfcf_nat,zx_tmp_fi2d)
	1938	IF (vars_defined) zx_tmp_fi2d(:)=topswcf_aero(:,2)*swradcorr(:)
	1939	CALL histwrite_phy(o_swtoacf_ant,zx_tmp_fi2d)
	1940	IF (vars_defined) zx_tmp_fi2d(:)=solswcf_aero(:,2)*swradcorr(:)
	1941	CALL histwrite_phy(o_swsrfcf_ant,zx_tmp_fi2d)
	1942	IF (vars_defined) zx_tmp_fi2d(:)=topswcf_aero(:,3)*swradcorr(:)
	1943	CALL histwrite_phy(o_swtoacf_zero,zx_tmp_fi2d)
	1944	IF (vars_defined) zx_tmp_fi2d(:)=solswcf_aero(:,3)*swradcorr(:)
	1945	CALL histwrite_phy(o_swsrfcf_zero,zx_tmp_fi2d)
	1946	ENDIF
[1898]	1947	!====MS forcing diagnostics
[1813]	1948	ENDIF
	1949	IF (ok_aie) THEN
[3435]	1950	IF (vars_defined) zx_tmp_fi2d(:)= topswai_aero*swradcorr
	1951	CALL histwrite_phy(o_topswai, zx_tmp_fi2d)
	1952
	1953	IF (vars_defined) zx_tmp_fi2d(:)=toplwai_aero*swradcorr
	1954	CALL histwrite_phy(o_toplwai, zx_tmp_fi2d)
	1955
	1956	IF (vars_defined) zx_tmp_fi2d(:)=solswai_aero*swradcorr
	1957	CALL histwrite_phy(o_solswai, zx_tmp_fi2d)
	1958
	1959	IF (vars_defined) zx_tmp_fi2d(:)=sollwai_aero*swradcorr
	1960	CALL histwrite_phy(o_sollwai, zx_tmp_fi2d)
[2003]	1961	ENDIF
[5084]	1962	IF (flag_aerosol.GT.0.AND.ok_cdnc) THEN
[1898]	1963	CALL histwrite_phy(o_scdnc, scdnc)
	1964	CALL histwrite_phy(o_cldncl, cldncl)
	1965	CALL histwrite_phy(o_reffclws, reffclws)
	1966	CALL histwrite_phy(o_reffclwc, reffclwc)
	1967	CALL histwrite_phy(o_cldnvi, cldnvi)
	1968	CALL histwrite_phy(o_lcc, lcc)
	1969	CALL histwrite_phy(o_lcc3d, lcc3d)
	1970	CALL histwrite_phy(o_lcc3dcon, lcc3dcon)
	1971	CALL histwrite_phy(o_lcc3dstra, lcc3dstra)
[3121]	1972	CALL histwrite_phy(o_icc3dcon, icc3dcon)
	1973	CALL histwrite_phy(o_icc3dstra, icc3dstra)
[3124]	1974	CALL histwrite_phy(o_cldicemxrat, zfice)
[3590]	1975	IF (vars_defined) zx_tmp_fi3d(:,:)=1-zfice(:,:)
	1976	CALL histwrite_phy(o_cldwatmxrat, zx_tmp_fi3d)
[1898]	1977	CALL histwrite_phy(o_reffclwtop, reffclwtop)
[1813]	1978	ENDIF
[1898]	1979	! Champs 3D:
[1813]	1980	IF (ok_ade .OR. ok_aie) then
[4389]	1981	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	1982	CALL histwrite_phy(o_ec550aer, ec550aer)
	1983	ENDIF
[1813]	1984	ENDIF
[4489]	1985
[1898]	1986	CALL histwrite_phy(o_lwcon, flwc)
	1987	CALL histwrite_phy(o_iwcon, fiwc)
	1988	CALL histwrite_phy(o_temp, t_seri)
	1989	CALL histwrite_phy(o_theta, theta)
	1990	CALL histwrite_phy(o_ovapinit, qx(:,:,ivap))
	1991	CALL histwrite_phy(o_ovap, q_seri)
	1992	CALL histwrite_phy(o_oliq, ql_seri)
[4575]	1993	!FC
	1994	CALL histwrite_phy(o_zxfluxt, zxfluxt)
	1995	CALL histwrite_phy(o_zxfluxq, zx_tmp_fi3d)
	1996	!FC
[3435]	1997
	1998	IF (vars_defined) zx_tmp_fi3d = ql_seri+qs_seri
	1999	CALL histwrite_phy(o_ocond, zx_tmp_fi3d)
[5383]	2000
	2001	IF (vars_defined) zx_tmp_fi3d = qs_seri
	2002	CALL histwrite_phy(o_oice, zx_tmp_fi3d)
[3435]	2003
[1898]	2004	CALL histwrite_phy(o_geop, zphi)
	2005	CALL histwrite_phy(o_vitu, u_seri)
	2006	CALL histwrite_phy(o_vitv, v_seri)
	2007	CALL histwrite_phy(o_vitw, omega)
	2008	CALL histwrite_phy(o_pres, pplay)
	2009	CALL histwrite_phy(o_paprs, paprs(:,1:klev))
[3435]	2010
	2011	IF (vars_defined) zx_tmp_fi3d = zphi/RG
	2012	CALL histwrite_phy(o_zfull,zx_tmp_fi3d)
[1813]	2013
[4523]	2014	IF (ok_bs) THEN
	2015	CALL histwrite_phy(o_qbs, qbs_seri)
	2016	ENDIF
	2017
[4619]	2018	IF (using_xios) THEN
[3048]	2019	!solbnd begin
	2020	#ifdef CPP_RRTM
[4619]	2021	IF (iflag_rrtm.EQ.1) THEN
	2022	IF (vars_defined) THEN
	2023	DO ISW=1, NSW
	2024	zx_tmp_fi3dsp(:,ISW) = swdn(:,klevp1)swradcorr(:)RSUN(ISW)
	2025	ENDDO
	2026	CALL histwrite_phy(o_solbnd, zx_tmp_fi3dsp)
	2027	ENDIF
	2028	ENDIF
[3048]	2029	#endif
	2030	!solbnd end
[4619]	2031	ENDIF
[3066]	2032	#endif
[3048]	2033
[5084]	2034	IF (flag_aerosol_strat.EQ.2) THEN
[2541]	2035	CALL histwrite_phy(o_stratomask, stratomask)
	2036	ENDIF
	2037
[1898]	2038	IF (vars_defined) THEN
[2377]	2039	zx_tmp_fi3d(:,1)= pphis(:)/RG
	2040	DO k = 2, klev
	2041	DO i = 1, klon
[3234]	2042	zx_tmp_fi3d(i,k) = zphi(i,k-1)/RG + &
[2377]	2043	(zphi(i,k)-zphi(i,k-1))/RG * &
[3234]	2044	(paprs(i,k)-pplay(i,k-1))/(pplay(i,k)-pplay(i,k-1))
[2377]	2045	ENDDO
	2046	ENDDO
[1898]	2047	ENDIF
	2048	CALL histwrite_phy(o_zhalf, zx_tmp_fi3d)
	2049	CALL histwrite_phy(o_rneb, cldfra)
	2050	CALL histwrite_phy(o_rnebcon, rnebcon)
	2051	CALL histwrite_phy(o_rnebls, rneb)
[2945]	2052	CALL histwrite_phy(o_rneblsvol, rneblsvol)
[1938]	2053	IF (vars_defined) THEN
	2054	DO k=1, klev
	2055	DO i=1, klon
	2056	zx_tmp_fi3d(i,k)=cldfra(i,k)*JrNt(i)
	2057	ENDDO
	2058	ENDDO
	2059	ENDIF
[1924]	2060	CALL histwrite_phy(o_rnebjn, zx_tmp_fi3d)
[1898]	2061	CALL histwrite_phy(o_rhum, zx_rh)
[5084]	2062	IF (iflag_ice_thermo .GT. 0) THEN
[3799]	2063	IF (vars_defined) zx_tmp_fi3d = zx_rhl * 100.
	2064	CALL histwrite_phy(o_rhl, zx_tmp_fi3d)
	2065	IF (vars_defined) zx_tmp_fi3d = zx_rhi * 100.
	2066	CALL histwrite_phy(o_rhi, zx_tmp_fi3d)
[3787]	2067	ENDIF
[4530]	2068
	2069	IF (ok_new_lscp) THEN
	2070	CALL histwrite_phy(o_pfraclr, pfraclr)
	2071	CALL histwrite_phy(o_pfracld, pfracld)
[5007]	2072	CALL histwrite_phy(o_cldfraliq, cldfraliq)
	2073	CALL histwrite_phy(o_sigma2_icefracturb, sigma2_icefracturb)
	2074	CALL histwrite_phy(o_mean_icefracturb, mean_icefracturb)
[4803]	2075	IF (ok_poprecip) THEN
[4830]	2076	CALL histwrite_phy(o_qrainlsc, qraindiag)
	2077	CALL histwrite_phy(o_qsnowlsc, qsnowdiag)
[4803]	2078	CALL histwrite_phy(o_dqreva, dqreva)
	2079	CALL histwrite_phy(o_dqrauto, dqrauto)
	2080	CALL histwrite_phy(o_dqrcol, dqrcol)
	2081	CALL histwrite_phy(o_dqrmelt, dqrmelt)
	2082	CALL histwrite_phy(o_dqrfreez, dqrfreez)
	2083	CALL histwrite_phy(o_dqssub, dqssub)
	2084	CALL histwrite_phy(o_dqsauto, dqsauto)
	2085	CALL histwrite_phy(o_dqsagg, dqsagg)
	2086	CALL histwrite_phy(o_dqsmelt, dqsmelt)
	2087	CALL histwrite_phy(o_dqsfreez, dqsfreez)
	2088	CALL histwrite_phy(o_dqsrim, dqsrim)
[5384]	2089	ELSE
	2090	CALL histwrite_phy(o_dqreva, dqreva)
	2091	CALL histwrite_phy(o_dqssub, dqssub)
	2092	CALL histwrite_phy(o_dqrauto, dqrauto)
	2093	CALL histwrite_phy(o_dqsauto, dqsauto)
[4803]	2094	ENDIF
[4530]	2095	ENDIF
[3780]	2096
[5204]	2097	!-- LSCP - condensation and supersaturation variables
	2098	IF (ok_ice_supersat) THEN
	2099	CALL histwrite_phy(o_cfseri, cf_seri)
	2100	CALL histwrite_phy(o_dcfdyn, d_cf_dyn)
	2101	CALL histwrite_phy(o_rvcseri, rvc_seri)
	2102	CALL histwrite_phy(o_drvcdyn, d_rvc_dyn)
	2103	CALL histwrite_phy(o_qsub, qsub)
	2104	CALL histwrite_phy(o_qissr, qissr)
	2105	CALL histwrite_phy(o_qcld, qcld)
	2106	CALL histwrite_phy(o_subfra, subfra)
	2107	CALL histwrite_phy(o_issrfra, issrfra)
	2108	CALL histwrite_phy(o_gammacond, gamma_cond)
	2109	CALL histwrite_phy(o_dcfsub, dcf_sub)
	2110	CALL histwrite_phy(o_dcfcon, dcf_con)
	2111	CALL histwrite_phy(o_dcfmix, dcf_mix)
	2112	CALL histwrite_phy(o_dqiadj, dqi_adj)
	2113	CALL histwrite_phy(o_dqisub, dqi_sub)
	2114	CALL histwrite_phy(o_dqicon, dqi_con)
	2115	CALL histwrite_phy(o_dqimix, dqi_mix)
	2116	CALL histwrite_phy(o_dqvcadj, dqvc_adj)
	2117	CALL histwrite_phy(o_dqvcsub, dqvc_sub)
	2118	CALL histwrite_phy(o_dqvccon, dqvc_con)
	2119	CALL histwrite_phy(o_dqvcmix, dqvc_mix)
	2120	CALL histwrite_phy(o_qsatl, qsatliq)
	2121	CALL histwrite_phy(o_qsati, qsatice)
	2122	ENDIF
	2123	!-- LSCP - aviation variables
	2124	IF (ok_plane_contrail) THEN
[4070]	2125	CALL histwrite_phy(o_Tcontr, Tcontr)
	2126	CALL histwrite_phy(o_qcontr, qcontr)
	2127	CALL histwrite_phy(o_qcontr2, qcontr2)
	2128	CALL histwrite_phy(o_fcontrN, fcontrN)
	2129	CALL histwrite_phy(o_fcontrP, fcontrP)
[5204]	2130	CALL histwrite_phy(o_dcfavi, dcf_avi)
	2131	CALL histwrite_phy(o_dqiavi, dqi_avi)
	2132	CALL histwrite_phy(o_dqvcavi, dqvc_avi)
	2133	CALL histwrite_phy(o_flight_dist, flight_dist)
[4070]	2134	ENDIF
[4062]	2135	IF (ok_plane_h2o) THEN
	2136	CALL histwrite_phy(o_flight_h2o, flight_h2o)
	2137	ENDIF
[3435]	2138
	2139	IF (vars_defined) zx_tmp_fi3d = wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd
	2140	CALL histwrite_phy(o_ozone, zx_tmp_fi3d)
[1813]	2141
[1898]	2142	IF (read_climoz == 2) THEN
[3435]	2143	IF (vars_defined) zx_tmp_fi3d = wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd
	2144	CALL histwrite_phy(o_ozone_light, zx_tmp_fi3d)
[1898]	2145	ENDIF
[1813]	2146
[2665]	2147	CALL histwrite_phy(o_duphy, d_u)
	2148
[1898]	2149	CALL histwrite_phy(o_dtphy, d_t)
[2496]	2150
[1898]	2151	CALL histwrite_phy(o_dqphy, d_qx(:,:,ivap))
[3435]	2152	IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,ivap),zmasse,zx_tmp_fi2d)
[2486]	2153	CALL histwrite_phy(o_dqphy2d, zx_tmp_fi2d)
[2496]	2154
	2155	CALL histwrite_phy(o_dqlphy, d_qx(:,:,iliq))
[3435]	2156	IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,iliq),zmasse,zx_tmp_fi2d)
[2496]	2157	CALL histwrite_phy(o_dqlphy2d, zx_tmp_fi2d)
	2158
[5084]	2159	IF (nqo.EQ.3) THEN
[2496]	2160	CALL histwrite_phy(o_dqsphy, d_qx(:,:,isol))
[3435]	2161	IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,isol),zmasse,zx_tmp_fi2d)
[2496]	2162	CALL histwrite_phy(o_dqsphy2d, zx_tmp_fi2d)
	2163	ELSE
	2164	zx_tmp_fi3d=0.0
	2165	CALL histwrite_phy(o_dqsphy, zx_tmp_fi3d)
	2166	zx_tmp_fi2d=0.0
	2167	CALL histwrite_phy(o_dqsphy2d, zx_tmp_fi2d)
	2168	ENDIF
	2169
[4523]	2170
	2171	IF (ok_bs) THEN
	2172	CALL histwrite_phy(o_dqbsphy, d_qx(:,:,ibs))
	2173	IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,ibs),zmasse,zx_tmp_fi2d)
	2174	CALL histwrite_phy(o_dqbsphy2d, zx_tmp_fi2d)
	2175	ELSE
	2176	zx_tmp_fi3d=0.0
	2177	CALL histwrite_phy(o_dqbsphy, zx_tmp_fi3d)
	2178	zx_tmp_fi2d=0.0
	2179	CALL histwrite_phy(o_dqbsphy2d, zx_tmp_fi2d)
	2180	ENDIF
	2181
[1898]	2182	DO nsrf=1, nbsrf
	2183	IF (vars_defined) zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
	2184	CALL histwrite_phy(o_albe_srf(nsrf), zx_tmp_fi2d)
[2243]	2185	IF (vars_defined) zx_tmp_fi2d(1 : klon) = z0m( 1 : klon, nsrf)
	2186	CALL histwrite_phy(o_z0m_srf(nsrf), zx_tmp_fi2d)
	2187	IF (vars_defined) zx_tmp_fi2d(1 : klon) = z0h( 1 : klon, nsrf)
	2188	CALL histwrite_phy(o_z0h_srf(nsrf), zx_tmp_fi2d)
[1898]	2189	IF (vars_defined) zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
	2190	CALL histwrite_phy(o_ages_srf(nsrf), zx_tmp_fi2d)
[2209]	2191	IF (vars_defined) zx_tmp_fi2d(1 : klon) = snow( 1 : klon, nsrf)
	2192	CALL histwrite_phy(o_snow_srf(nsrf), zx_tmp_fi2d)
[1898]	2193	ENDDO !nsrf=1, nbsrf
	2194	CALL histwrite_phy(o_alb1, albsol1)
	2195	CALL histwrite_phy(o_alb2, albsol2)
	2196	!FH Sorties pour la couche limite
[2690]	2197	IF (iflag_pbl>1) THEN
[1898]	2198	zx_tmp_fi3d=0.
	2199	IF (vars_defined) THEN
[2690]	2200	DO nsrf=1,nbsrf
	2201	DO k=1,klev
[1898]	2202	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) &
	2203	+pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
[2690]	2204	ENDDO
	2205	ENDDO
[1898]	2206	ENDIF
	2207	CALL histwrite_phy(o_tke, zx_tmp_fi3d)
[3780]	2208	CALL histwrite_phy(o_tke_max, zx_tmp_fi3d)
[1898]	2209	ENDIF
[1813]	2210
[1898]	2211	CALL histwrite_phy(o_kz, coefh(:,:,is_ave))
[1813]	2212
[1898]	2213	CALL histwrite_phy(o_kz_max, coefh(:,:,is_ave))
[1813]	2214
[1898]	2215	CALL histwrite_phy(o_clwcon, clwcon0)
[4737]	2216
[1898]	2217	CALL histwrite_phy(o_dtdyn, d_t_dyn)
[2496]	2218
[1898]	2219	CALL histwrite_phy(o_dqdyn, d_q_dyn)
[2496]	2220
[2499]	2221	CALL histwrite_phy(o_dqdyn2d,d_q_dyn2d)
	2222
[2496]	2223	CALL histwrite_phy(o_dqldyn, d_ql_dyn)
	2224
[2499]	2225	CALL histwrite_phy(o_dqldyn2d, d_ql_dyn2d)
	2226
[2496]	2227	CALL histwrite_phy(o_dqsdyn, d_qs_dyn)
	2228
[2499]	2229	CALL histwrite_phy(o_dqsdyn2d, d_qs_dyn2d)
	2230
[4523]	2231	IF (ok_bs) THEN
[4737]	2232	CALL histwrite_phy(o_dqbsdyn, d_qbs_dyn)
	2233	CALL histwrite_phy(o_dqbsdyn2d, d_qbs_dyn2d)
[4523]	2234	ENDIF
	2235
[1898]	2236	CALL histwrite_phy(o_dudyn, d_u_dyn)
	2237	CALL histwrite_phy(o_dvdyn, d_v_dyn)
[1813]	2238
[1898]	2239	IF (vars_defined) THEN
	2240	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
	2241	ENDIF
	2242	CALL histwrite_phy(o_dtcon, zx_tmp_fi3d)
[5084]	2243	IF (iflag_thermals.EQ.0) THEN
[1898]	2244	IF (vars_defined) THEN
	2245	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
	2246	d_t_ajsb(1:klon,1:klev)/pdtphys
	2247	ENDIF
	2248	CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
[5084]	2249	ELSE IF(iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
[1898]	2250	IF (vars_defined) THEN
	2251	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
	2252	d_t_ajs(1:klon,1:klev)/pdtphys + &
	2253	d_t_wake(1:klon,1:klev)/pdtphys
	2254	ENDIF
	2255	CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
[3453]	2256	ENDIF
[1898]	2257	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
	2258	CALL histwrite_phy(o_ducon, zx_tmp_fi3d)
	2259	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
	2260	CALL histwrite_phy(o_dvcon, zx_tmp_fi3d)
	2261	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
	2262	CALL histwrite_phy(o_dqcon, zx_tmp_fi3d)
[3435]	2263	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
[2486]	2264	CALL histwrite_phy(o_dqcon2d, zx_tmp_fi2d)
[1813]	2265
[5084]	2266	IF (iflag_thermals.EQ.0) THEN
[1898]	2267	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
	2268	CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d)
[5084]	2269	ELSE IF (iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
[1898]	2270	IF (vars_defined) THEN
	2271	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys + &
	2272	d_q_ajs(1:klon,1:klev)/pdtphys + &
	2273	d_q_wake(1:klon,1:klev)/pdtphys
	2274	ENDIF
	2275	CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d)
	2276	ENDIF
[1813]	2277
[1898]	2278	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
	2279	CALL histwrite_phy(o_dtlsc, zx_tmp_fi3d)
[4737]	2280	IF (vars_defined) zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+d_t_eva(1:klon,1:klev))/pdtphys
[1898]	2281	CALL histwrite_phy(o_dtlschr, zx_tmp_fi3d)
	2282	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
	2283	CALL histwrite_phy(o_dqlsc, zx_tmp_fi3d)
[3435]	2284	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
[2486]	2285	CALL histwrite_phy(o_dqlsc2d, zx_tmp_fi2d)
[1898]	2286	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=beta_prec(1:klon,1:klev)
	2287	CALL histwrite_phy(o_beta_prec, zx_tmp_fi3d)
[1813]	2288	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1898]	2289	! Sorties specifiques a la separation thermiques/non thermiques
[2690]	2290	IF (iflag_thermals>=1) THEN
	2291	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys
[1898]	2292	CALL histwrite_phy(o_dtlscth, zx_tmp_fi3d)
[2690]	2293	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys
[1898]	2294	CALL histwrite_phy(o_dtlscst, zx_tmp_fi3d)
[2690]	2295	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys
[1898]	2296	CALL histwrite_phy(o_dqlscth, zx_tmp_fi3d)
[3435]	2297	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
[2486]	2298	CALL histwrite_phy(o_dqlscth2d, zx_tmp_fi2d)
[2690]	2299	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys
[1898]	2300	CALL histwrite_phy(o_dqlscst, zx_tmp_fi3d)
[3435]	2301	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
[2486]	2302	CALL histwrite_phy(o_dqlscst2d, zx_tmp_fi2d)
[1898]	2303	CALL histwrite_phy(o_plulth, plul_th)
	2304	CALL histwrite_phy(o_plulst, plul_st)
	2305	IF (vars_defined) THEN
[3453]	2306	DO i=1,klon
[3171]	2307	zx_tmp_fi2d(1:klon)=lmax_th(:)
[3453]	2308	ENDDO
[3171]	2309	ENDIF
	2310	CALL histwrite_phy(o_lmaxth, zx_tmp_fi2d)
	2311	IF (vars_defined) THEN
[2690]	2312	DO k=1,klev
	2313	DO i=1,klon
	2314	IF (ptconvth(i,k)) THEN
[1898]	2315	zx_tmp_fi3d(i,k)=1.
[2690]	2316	ELSE
[1898]	2317	zx_tmp_fi3d(i,k)=0.
[2690]	2318	ENDIF
	2319	ENDDO
	2320	ENDDO
[1898]	2321	ENDIF
	2322	CALL histwrite_phy(o_ptconvth, zx_tmp_fi3d)
[2690]	2323	ENDIF ! iflag_thermals>=1
[1813]	2324	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[3171]	2325	zpt_conv = 0.
	2326	WHERE (ptconv) zpt_conv = 1.
	2327	CALL histwrite_phy(o_ptconv, zpt_conv)
	2328	!! IF (vars_defined) zx_tmp_fi2d=float(itau_con)/float(itap)
	2329	!! CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
	2330	IF (vars_defined) THEN
	2331	zpt_conv2d(:) = 0.
	2332	DO k=1,klev
	2333	WHERE (ptconv(:,k)) zpt_conv2d(:) = 1.
	2334	ENDDO
	2335	ENDIF
	2336	CALL histwrite_phy(o_ftime_deepcv, zpt_conv2d)
	2337	IF (vars_defined) THEN
	2338	zx_tmp_fi2d(:) = 0.
	2339	DO k=1,klev
	2340	WHERE (ptconvth(:,k)) zx_tmp_fi2d(:) = 1.
	2341	ENDDO
	2342	ENDIF
	2343	CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
	2344	IF (vars_defined) THEN
	2345	zx_tmp_fi2d(:) = max(zx_tmp_fi2d(:),zpt_conv2d(:))
	2346	ENDIF
	2347	CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
	2348	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1898]	2349	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
	2350	CALL histwrite_phy(o_dtvdf, zx_tmp_fi3d)
[2690]	2351	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_diss(1:klon,1:klev)/pdtphys
[1898]	2352	CALL histwrite_phy(o_dtdis, zx_tmp_fi3d)
	2353	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
	2354	CALL histwrite_phy(o_dqvdf, zx_tmp_fi3d)
[3435]	2355	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
[2486]	2356	CALL histwrite_phy(o_dqvdf2d, zx_tmp_fi2d)
[1898]	2357	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
	2358	CALL histwrite_phy(o_dteva, zx_tmp_fi3d)
	2359	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
	2360	CALL histwrite_phy(o_dqeva, zx_tmp_fi3d)
[3435]	2361	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
[2486]	2362	CALL histwrite_phy(o_dqeva2d, zx_tmp_fi2d)
[1898]	2363	CALL histwrite_phy(o_ratqs, ratqs)
	2364	IF (vars_defined) THEN
	2365	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys - &
	2366	d_t_ajsb(1:klon,1:klev)/pdtphys
	2367	ENDIF
	2368	CALL histwrite_phy(o_dtthe, zx_tmp_fi3d)
[1903]	2369	IF (vars_defined) THEN
[1943]	2370	zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys
[1903]	2371	ENDIF
	2372	CALL histwrite_phy(o_duthe, zx_tmp_fi3d)
	2373	IF (vars_defined) THEN
[1943]	2374	zx_tmp_fi3d(1:klon,1:klev)=d_v_ajs(1:klon,1:klev)/pdtphys
[1903]	2375	ENDIF
	2376	CALL histwrite_phy(o_dvthe, zx_tmp_fi3d)
	2377
[4523]	2378	IF (ok_bs) THEN
	2379	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_qbs_vdf(1:klon,1:klev)/pdtphys
	2380	CALL histwrite_phy(o_dqbsvdf, zx_tmp_fi3d)
[5050]	2381	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_qbs_bsss(1:klon,1:klev)/pdtphys
[4523]	2382	CALL histwrite_phy(o_dqbsbs, zx_tmp_fi3d)
[5050]	2383	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_bsss(1:klon,1:klev)/pdtphys
[4523]	2384	CALL histwrite_phy(o_dqbs, zx_tmp_fi3d)
[5050]	2385	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_bsss(1:klon,1:klev)/pdtphys
[4523]	2386	CALL histwrite_phy(o_dtbs, zx_tmp_fi3d)
	2387	ENDIF
	2388
[1813]	2389	IF (iflag_thermals>=1) THEN
[1898]	2390	! Pour l instant 0 a y reflichir pour les thermiques
[3171]	2391	! regroupe avec ftime_deepcv et ftime_con
	2392	!!zx_tmp_fi2d=0.
	2393	!!CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
[1898]	2394	CALL histwrite_phy(o_f_th, fm_therm)
	2395	CALL histwrite_phy(o_e_th, entr_therm)
	2396	CALL histwrite_phy(o_w_th, zw2)
	2397	CALL histwrite_phy(o_q_th, zqasc)
	2398	CALL histwrite_phy(o_a_th, fraca)
[2958]	2399	CALL histwrite_phy(o_cloudth_sth, cloudth_sth)
	2400	CALL histwrite_phy(o_cloudth_senv, cloudth_senv)
	2401	CALL histwrite_phy(o_cloudth_sigmath, cloudth_sigmath)
	2402	CALL histwrite_phy(o_cloudth_sigmaenv, cloudth_sigmaenv)
[1898]	2403	CALL histwrite_phy(o_d_th, detr_therm)
	2404	CALL histwrite_phy(o_f0_th, f0)
	2405	CALL histwrite_phy(o_zmax_th, zmax_th)
	2406	IF (vars_defined) THEN
[4737]	2407	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys-d_q_ajsb(1:klon,1:klev)/pdtphys
[1898]	2408	ENDIF
	2409	CALL histwrite_phy(o_dqthe, zx_tmp_fi3d)
[3435]	2410	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
[2486]	2411	CALL histwrite_phy(o_dqthe2d, zx_tmp_fi2d)
[1898]	2412	ENDIF !iflag_thermals
[2690]	2413	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
[1898]	2414	CALL histwrite_phy(o_dtajs, zx_tmp_fi3d)
[2690]	2415	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
[1898]	2416	CALL histwrite_phy(o_dqajs, zx_tmp_fi3d)
[3435]	2417	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
[2486]	2418	CALL histwrite_phy(o_dqajs2d, zx_tmp_fi2d)
[2690]	2419	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys
[1898]	2420	CALL histwrite_phy(o_dtswr, zx_tmp_fi3d)
[2690]	2421	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_sw0(1:klon,1:klev)/pdtphys
[1898]	2422	CALL histwrite_phy(o_dtsw0, zx_tmp_fi3d)
[2690]	2423	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lwr(1:klon,1:klev)/pdtphys
[1898]	2424	CALL histwrite_phy(o_dtlwr, zx_tmp_fi3d)
[2690]	2425	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lw0(1:klon,1:klev)/pdtphys
[1898]	2426	CALL histwrite_phy(o_dtlw0, zx_tmp_fi3d)
[2690]	2427	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)/pdtphys
[1898]	2428	CALL histwrite_phy(o_dtec, zx_tmp_fi3d)
[2690]	2429	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
[1898]	2430	CALL histwrite_phy(o_duvdf, zx_tmp_fi3d)
[2690]	2431	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
[1898]	2432	CALL histwrite_phy(o_dvvdf, zx_tmp_fi3d)
[1813]	2433	IF (ok_orodr) THEN
[2690]	2434	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
[1898]	2435	CALL histwrite_phy(o_duoro, zx_tmp_fi3d)
[2690]	2436	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
[1898]	2437	CALL histwrite_phy(o_dvoro, zx_tmp_fi3d)
[2690]	2438	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys
[1898]	2439	CALL histwrite_phy(o_dtoro, zx_tmp_fi3d)
[1813]	2440	ENDIF
[1898]	2441	IF (ok_orolf) THEN
[2690]	2442	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
[1898]	2443	CALL histwrite_phy(o_dulif, zx_tmp_fi3d)
[1903]	2444
[2690]	2445	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
[1903]	2446	CALL histwrite_phy(o_dvlif, zx_tmp_fi3d)
	2447
[2690]	2448	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys
[1903]	2449	CALL histwrite_phy(o_dtlif, zx_tmp_fi3d)
[1813]	2450	ENDIF
	2451
	2452	IF (ok_hines) THEN
[3435]	2453	IF (vars_defined) zx_tmp_fi3d=du_gwd_hines/pdtphys
	2454	CALL histwrite_phy(o_du_gwd_hines, zx_tmp_fi3d)
	2455
	2456	IF (vars_defined) zx_tmp_fi3d= dv_gwd_hines/pdtphys
	2457	CALL histwrite_phy(o_dv_gwd_hines, zx_tmp_fi3d)
	2458
[2690]	2459	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys
[1898]	2460	CALL histwrite_phy(o_dthin, zx_tmp_fi3d)
[2333]	2461	CALL histwrite_phy(o_ustr_gwd_hines, zustr_gwd_hines)
	2462	CALL histwrite_phy(o_vstr_gwd_hines, zvstr_gwd_hines)
[2690]	2463	ENDIF
[2333]	2464
[2690]	2465	IF (.not. ok_hines .and. ok_gwd_rando) THEN
[3435]	2466	IF (vars_defined) zx_tmp_fi3d=du_gwd_front / pdtphys
	2467	CALL histwrite_phy(o_du_gwd_front, zx_tmp_fi3d)
	2468
	2469	IF (vars_defined) zx_tmp_fi3d=dv_gwd_front / pdtphys
	2470	CALL histwrite_phy(o_dv_gwd_front, zx_tmp_fi3d)
	2471
[2333]	2472	CALL histwrite_phy(o_ustr_gwd_front, zustr_gwd_front)
	2473	CALL histwrite_phy(o_vstr_gwd_front, zvstr_gwd_front)
[1898]	2474	ENDIF
[1938]	2475
[2690]	2476	IF (ok_gwd_rando) THEN
[3435]	2477	IF (vars_defined) zx_tmp_fi3d=du_gwd_rando / pdtphys
	2478	CALL histwrite_phy(o_du_gwd_rando, zx_tmp_fi3d)
	2479
	2480	IF (vars_defined) zx_tmp_fi3d=dv_gwd_rando / pdtphys
	2481	CALL histwrite_phy(o_dv_gwd_rando, zx_tmp_fi3d)
[2333]	2482	CALL histwrite_phy(o_ustr_gwd_rando, zustr_gwd_rando)
[1938]	2483	CALL histwrite_phy(o_vstr_gwd_rando, zvstr_gwd_rando)
[2333]	2484	CALL histwrite_phy(o_east_gwstress, east_gwstress )
	2485	CALL histwrite_phy(o_west_gwstress, west_gwstress )
[2690]	2486	ENDIF
[1938]	2487
[2690]	2488	IF (ok_qch4) THEN
[4755]	2489	IF (vars_defined) zx_tmp_fi3d=d_q_ch4
[3435]	2490	CALL histwrite_phy(o_dqch4, zx_tmp_fi3d)
[2136]	2491	ENDIF
[3435]	2492
	2493	IF (vars_defined) THEN
	2494	DO k=1, klevp1
	2495	zx_tmp_fi3d1(:,k)=swup(:,k)*swradcorr(:)
	2496	ENDDO
	2497	ENDIF
	2498
	2499	CALL histwrite_phy(o_rsu, zx_tmp_fi3d1)
[2136]	2500
[3435]	2501	IF (vars_defined) THEN
	2502	DO k=1, klevp1
	2503	zx_tmp_fi3d1(:,k)=swdn(:,k)*swradcorr(:)
	2504	ENDDO
	2505	ENDIF
	2506
[2366]	2507	CALL histwrite_phy(o_rsd, zx_tmp_fi3d1)
[3435]	2508
	2509	IF (vars_defined) THEN
	2510	DO k=1, klevp1
	2511	zx_tmp_fi3d1(:,k)=swup0(:,k)*swradcorr(:)
	2512	ENDDO
	2513	ENDIF
	2514
[2366]	2515	CALL histwrite_phy(o_rsucs, zx_tmp_fi3d1)
[3435]	2516
[4389]	2517	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	2518	IF (vars_defined) THEN
	2519	DO k=1, klevp1
	2520	zx_tmp_fi3d1(:,k)=swupc0(:,k)*swradcorr(:)
	2521	ENDDO
	2522	ENDIF
	2523	CALL histwrite_phy(o_rsucsaf, zx_tmp_fi3d1)
[3435]	2524	ENDIF
	2525
	2526	IF (vars_defined) THEN
	2527	DO k=1, klevp1
	2528	zx_tmp_fi3d1(:,k)=swdn0(:,k)*swradcorr(:)
	2529	ENDDO
	2530	ENDIF
[2366]	2531	CALL histwrite_phy(o_rsdcs, zx_tmp_fi3d1)
[3435]	2532
[4389]	2533	IF (type_trac/='inca' .OR. config_inca=='aeNP') THEN
[4106]	2534	IF (vars_defined) THEN
	2535	DO k=1, klevp1
	2536	zx_tmp_fi3d1(:,k)=swdnc0(:,k)*swradcorr(:)
	2537	ENDDO
	2538	ENDIF
	2539	CALL histwrite_phy(o_rsdcsaf, zx_tmp_fi3d1)
[3435]	2540	ENDIF
[2366]	2541
[1898]	2542	CALL histwrite_phy(o_rlu, lwup)
	2543	CALL histwrite_phy(o_rld, lwdn)
	2544	CALL histwrite_phy(o_rlucs, lwup0)
	2545	CALL histwrite_phy(o_rldcs, lwdn0)
[2366]	2546
[2690]	2547	IF (vars_defined) THEN
[1898]	2548	zx_tmp_fi3d(1:klon,1:klev)=d_t(1:klon,1:klev)+ &
	2549	d_t_dyn(1:klon,1:klev)
	2550	ENDIF
	2551	CALL histwrite_phy(o_tnt, zx_tmp_fi3d)
[2366]	2552
[2690]	2553	IF (vars_defined) THEN
[2194]	2554	zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys + &
	2555	d_t_lwr(1:klon,1:klev)/pdtphys
[1898]	2556	ENDIF
	2557	CALL histwrite_phy(o_tntr, zx_tmp_fi3d)
[2690]	2558	IF (vars_defined) THEN
[1898]	2559	zx_tmp_fi3d(1:klon,1:klev)= (d_t_lsc(1:klon,1:klev)+ &
	2560	d_t_eva(1:klon,1:klev)+ &
	2561	d_t_vdf(1:klon,1:klev))/pdtphys
	2562	ENDIF
	2563	CALL histwrite_phy(o_tntscpbl, zx_tmp_fi3d)
[2690]	2564	IF (vars_defined) THEN
[1898]	2565	zx_tmp_fi3d(1:klon,1:klev)=d_qx(1:klon,1:klev,ivap)+ &
	2566	d_q_dyn(1:klon,1:klev)
	2567	ENDIF
	2568	CALL histwrite_phy(o_tnhus, zx_tmp_fi3d)
[2690]	2569	IF (vars_defined) THEN
[1898]	2570	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys+ &
	2571	d_q_eva(1:klon,1:klev)/pdtphys
	2572	ENDIF
	2573	CALL histwrite_phy(o_tnhusscpbl, zx_tmp_fi3d)
	2574	CALL histwrite_phy(o_evu, coefm(:,:,is_ave))
[2690]	2575	IF (vars_defined) THEN
[1898]	2576	zx_tmp_fi3d(1:klon,1:klev)=q_seri(1:klon,1:klev)+ &
	2577	ql_seri(1:klon,1:klev)
	2578	ENDIF
	2579	CALL histwrite_phy(o_h2o, zx_tmp_fi3d)
[2690]	2580	IF (iflag_con >= 3) THEN
	2581	IF (vars_defined) THEN
[1813]	2582	zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+ &
	2583	dnwd0(1:klon,1:klev))
[1898]	2584	ENDIF
	2585	CALL histwrite_phy(o_mcd, zx_tmp_fi3d)
[2690]	2586	IF (vars_defined) THEN
[1813]	2587	zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) + &
	2588	dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev)
[1898]	2589	ENDIF
	2590	CALL histwrite_phy(o_dmc, zx_tmp_fi3d)
[2690]	2591	ELSE IF (iflag_con == 2) THEN
[1898]	2592	CALL histwrite_phy(o_mcd, pmfd)
[3435]	2593	IF (vars_defined) zx_tmp_fi3d = pmfu + pmfd
	2594	CALL histwrite_phy(o_dmc, zx_tmp_fi3d)
[2690]	2595	ENDIF
[1898]	2596	CALL histwrite_phy(o_ref_liq, ref_liq)
	2597	CALL histwrite_phy(o_ref_ice, ref_ice)
[2989]	2598	!
	2599	IF (ok_4xCO2atm) THEN
[2690]	2600	IF (vars_defined) zx_tmp_fi2d(:) = swupp(:,klevp1)*swradcorr(:)
[1898]	2601	CALL histwrite_phy(o_rsut4co2, zx_tmp_fi2d)
[2690]	2602	IF (vars_defined) zx_tmp_fi2d(:) = lwupp(:,klevp1)
[1898]	2603	CALL histwrite_phy(o_rlut4co2, zx_tmp_fi2d)
[2690]	2604	IF (vars_defined) zx_tmp_fi2d(:) = swup0p(:,klevp1)*swradcorr(:)
[1898]	2605	CALL histwrite_phy(o_rsutcs4co2, zx_tmp_fi2d)
[2690]	2606	IF (vars_defined) zx_tmp_fi2d(:) = lwup0p(:,klevp1)
[1898]	2607	CALL histwrite_phy(o_rlutcs4co2, zx_tmp_fi2d)
[3435]	2608	IF (vars_defined) THEN
	2609	DO k=1, klevp1
	2610	zx_tmp_fi3d1(:,k)=swupp(:,k)*swradcorr(:)
	2611	ENDDO
	2612	ENDIF
[2366]	2613	CALL histwrite_phy(o_rsu4co2, zx_tmp_fi3d1)
[3435]	2614	IF (vars_defined) THEN
	2615	DO k=1, klevp1
	2616	zx_tmp_fi3d1(:,k)=swup0p(:,k)*swradcorr(:)
	2617	ENDDO
	2618	ENDIF
[2366]	2619	CALL histwrite_phy(o_rsucs4co2, zx_tmp_fi3d1)
[3435]	2620	IF (vars_defined) THEN
	2621	DO k=1, klevp1
	2622	zx_tmp_fi3d1(:,k)=swdnp(:,k)*swradcorr(:)
	2623	ENDDO
	2624	ENDIF
[2366]	2625	CALL histwrite_phy(o_rsd4co2, zx_tmp_fi3d1)
[3435]	2626	IF (vars_defined) THEN
	2627	DO k=1, klevp1
	2628	zx_tmp_fi3d1(:,k)=swdn0p(:,k)*swradcorr(:)
	2629	ENDDO
	2630	ENDIF
[2366]	2631	CALL histwrite_phy(o_rsdcs4co2, zx_tmp_fi3d1)
[1898]	2632	CALL histwrite_phy(o_rlu4co2, lwupp)
	2633	CALL histwrite_phy(o_rlucs4co2, lwup0p)
	2634	CALL histwrite_phy(o_rld4co2, lwdnp)
	2635	CALL histwrite_phy(o_rldcs4co2, lwdn0p)
[2989]	2636	ENDIF !ok_4xCO2atm
[1828]	2637	!!!!!!!!!!!! Sorties niveaux de pression NMC !!!!!!!!!!!!!!!!!!!!
[4619]	2638	IF (.NOT. using_xios) THEN
[5267]	2639	IF (.NOT.ok_all_xml) THEN
[4619]	2640	! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX:
	2641	! Champs interpolles sur des niveaux de pression
	2642	DO iff=7, nfiles-1 !--OB: here we deal with files 7,8,9
[1828]	2643
[4619]	2644	CALL histwrite_phy(o_tnondef,tnondef(:,:,iff-6),iff)
	2645	CALL histwrite_phy(o_ta,twriteSTD(:,:,iff-6),iff)
	2646	CALL histwrite_phy(o_zg,phiwriteSTD(:,:,iff-6),iff)
	2647	CALL histwrite_phy(o_hus,qwriteSTD(:,:,iff-6),iff)
	2648	CALL histwrite_phy(o_hur,rhwriteSTD(:,:,iff-6),iff)
	2649	CALL histwrite_phy(o_ua,uwriteSTD(:,:,iff-6),iff)
	2650	CALL histwrite_phy(o_va,vwriteSTD(:,:,iff-6),iff)
	2651	CALL histwrite_phy(o_wap,wwriteSTD(:,:,iff-6),iff)
[2690]	2652	IF (vars_defined) THEN
[4619]	2653	DO k=1, nlevSTD
	2654	DO i=1, klon
	2655	IF (tnondef(i,k,iff-6).NE.missing_val) THEN
	2656	IF (freq_outNMC(iff-6).LT.0) THEN
	2657	freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6)
	2658	ELSE
	2659	freq_moyNMC(iff-6)=freq_outNMC(iff-6)/freq_calNMC(iff-6)
	2660	ENDIF
	2661	zx_tmp_fi3d_STD(i,k) = (100.*tnondef(i,k,iff-6))/freq_moyNMC(iff-6)
	2662	ELSE
	2663	zx_tmp_fi3d_STD(i,k) = missing_val
	2664	ENDIF
	2665	ENDDO
	2666	ENDDO
	2667	ENDIF
	2668	CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD,iff)
	2669	IF (vars_defined) THEN
	2670	DO k=1, nlevSTD
	2671	DO i=1, klon
	2672	IF (O3sumSTD(i,k,iff-6).NE.missing_val) THEN
	2673	zx_tmp_fi3d_STD(i,k) = O3sumSTD(i,k,iff-6) * 1.e+9
	2674	ELSE
	2675	zx_tmp_fi3d_STD(i,k) = missing_val
	2676	ENDIF
	2677	ENDDO
	2678	ENDDO !k=1, nlevSTD
	2679	ENDIF
	2680	CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD,iff)
	2681	IF (read_climoz == 2) THEN
	2682	IF (vars_defined) THEN
	2683	DO k=1, nlevSTD
[1898]	2684	DO i=1, klon
[2690]	2685	IF (O3daysumSTD(i,k,iff-6).NE.missing_val) THEN
[1971]	2686	zx_tmp_fi3d_STD(i,k) = O3daysumSTD(i,k,iff-6) * 1.e+9
[1898]	2687	ELSE
[1971]	2688	zx_tmp_fi3d_STD(i,k) = missing_val
[1898]	2689	ENDIF
	2690	ENDDO
[4619]	2691	ENDDO !k=1, nlevSTD
	2692	ENDIF
	2693	CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d_STD,iff)
[1898]	2694	ENDIF
[4619]	2695	CALL histwrite_phy(o_uxv,uvsumSTD(:,:,iff-6),iff)
	2696	CALL histwrite_phy(o_vxq,vqsumSTD(:,:,iff-6),iff)
	2697	CALL histwrite_phy(o_vxT,vTsumSTD(:,:,iff-6),iff)
	2698	CALL histwrite_phy(o_wxq,wqsumSTD(:,:,iff-6),iff)
	2699	CALL histwrite_phy(o_vxphi,vphisumSTD(:,:,iff-6),iff)
	2700	CALL histwrite_phy(o_wxT,wTsumSTD(:,:,iff-6),iff)
	2701	CALL histwrite_phy(o_uxu,u2sumSTD(:,:,iff-6),iff)
	2702	CALL histwrite_phy(o_vxv,v2sumSTD(:,:,iff-6),iff)
	2703	CALL histwrite_phy(o_TxT,T2sumSTD(:,:,iff-6),iff)
	2704	ENDDO !nfiles
	2705	ENDIF
	2706	ENDIF !.NOT. using_xios
[2282]	2707
[5267]	2708
[4619]	2709	IF (using_xios) THEN
	2710	IF (ok_all_xml) THEN
	2711	! DO iff=7, nfiles
	2712
[2282]	2713	! CALL histwrite_phy(o_tnondef,tnondef(:,:,3))
	2714	CALL histwrite_phy(o_ta,tlevSTD(:,:))
	2715	CALL histwrite_phy(o_zg,philevSTD(:,:))
	2716	CALL histwrite_phy(o_hus,qlevSTD(:,:))
	2717	CALL histwrite_phy(o_hur,rhlevSTD(:,:))
	2718	CALL histwrite_phy(o_ua,ulevSTD(:,:))
	2719	CALL histwrite_phy(o_va,vlevSTD(:,:))
	2720	CALL histwrite_phy(o_wap,wlevSTD(:,:))
[2690]	2721	! IF (vars_defined) THEN
[2282]	2722	! DO k=1, nlevSTD
	2723	! DO i=1, klon
[2690]	2724	! IF (tnondef(i,k,3).NE.missing_val) THEN
	2725	! IF (freq_outNMC(iff-6).LT.0) THEN
[2282]	2726	! freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6)
	2727	! ELSE
	2728	! freq_moyNMC(iff-6)=freq_outNMC(iff-6)/freq_calNMC(iff-6)
	2729	! ENDIF
	2730	! zx_tmp_fi3d_STD(i,k) = (100.*tnondef(i,k,3))/freq_moyNMC(iff-6)
	2731	! ELSE
	2732	! zx_tmp_fi3d_STD(i,k) = missing_val
	2733	! ENDIF
	2734	! ENDDO
	2735	! ENDDO
	2736	! ENDIF
	2737	! CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD)
[2690]	2738	IF (vars_defined) THEN
[2282]	2739	DO k=1, nlevSTD
	2740	DO i=1, klon
[5084]	2741	IF (O3STD(i,k).NE.missing_val) THEN
[2282]	2742	zx_tmp_fi3d_STD(i,k) = O3STD(i,k) * 1.e+9
	2743	ELSE
	2744	zx_tmp_fi3d_STD(i,k) = missing_val
	2745	ENDIF
	2746	ENDDO
	2747	ENDDO !k=1, nlevSTD
	2748	ENDIF
	2749	CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD)
[2690]	2750	IF (read_climoz == 2) THEN
	2751	IF (vars_defined) THEN
[2282]	2752	DO k=1, nlevSTD
	2753	DO i=1, klon
[5084]	2754	IF (O3daySTD(i,k).NE.missing_val) THEN
[2282]	2755	zx_tmp_fi3d_STD(i,k) = O3daySTD(i,k) * 1.e+9
	2756	ELSE
	2757	zx_tmp_fi3d_STD(i,k) = missing_val
	2758	ENDIF
	2759	ENDDO
	2760	ENDDO !k=1, nlevSTD
	2761	ENDIF
	2762	CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d_STD)
[2690]	2763	ENDIF
[2282]	2764	CALL histwrite_phy(o_uxv,uvSTD(:,:))
	2765	CALL histwrite_phy(o_vxq,vqSTD(:,:))
	2766	CALL histwrite_phy(o_vxT,vTSTD(:,:))
	2767	CALL histwrite_phy(o_wxq,wqSTD(:,:))
	2768	CALL histwrite_phy(o_vxphi,vphiSTD(:,:))
	2769	CALL histwrite_phy(o_wxT,wTSTD(:,:))
	2770	CALL histwrite_phy(o_uxu,u2STD(:,:))
	2771	CALL histwrite_phy(o_vxv,v2STD(:,:))
	2772	CALL histwrite_phy(o_TxT,T2STD(:,:))
	2773	! ENDDO !nfiles
[4619]	2774	ENDIF
	2775	ENDIF !using_xios
[1828]	2776	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2973]	2777	IF (iflag_phytrac == 1 ) then
[4170]	2778	!
[4389]	2779	IF (type_trac == 'co2i') THEN
[4056]	2780	itr = 0
	2781	DO iq = 1, nqtot
[4071]	2782	IF(.NOT.tracers(iq)%isInPhysics) CYCLE
[4056]	2783	itr = itr + 1
[4170]	2784	! write(,) 'phys_output_write_mod 2370: itr=',itr
[2704]	2785	!--3D fields
[4056]	2786	CALL histwrite_phy(o_trac(itr), tr_seri(:,:,itr))
	2787	CALL histwrite_phy(o_dtr_vdf(itr),d_tr_cl(:,:,itr))
	2788	CALL histwrite_phy(o_dtr_the(itr),d_tr_th(:,:,itr))
	2789	CALL histwrite_phy(o_dtr_con(itr),d_tr_cv(:,:,itr))
[4170]	2790	!--2D fields
	2791	!--CO2 burden
[1898]	2792	zx_tmp_fi2d=0.
[2690]	2793	IF (vars_defined) THEN
[1898]	2794	DO k=1,klev
[4056]	2795	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,itr)
[1898]	2796	ENDDO
	2797	ENDIF
[4056]	2798	CALL histwrite_phy(o_trac_cum(itr), zx_tmp_fi2d)
[3550]	2799	ENDDO !--iq
[4170]	2800	!--CO2 net fluxes
	2801	CALL histwrite_phy(o_flx_co2_land, fco2_land)
	2802	CALL histwrite_phy(o_flx_co2_ocean, fco2_ocean)
	2803	CALL histwrite_phy(o_flx_co2_ocean_cor, fco2_ocean_cor)
	2804	CALL histwrite_phy(o_flx_co2_land_cor, fco2_land_cor)
	2805	CALL histwrite_phy(o_flx_co2_ff, fco2_ff)
	2806	CALL histwrite_phy(o_flx_co2_bb, fco2_bb)
	2807
[4389]	2808	ELSE IF (type_trac == 'inco') THEN
[4056]	2809	itr = 0
	2810	DO iq = 1, nqtot
[4071]	2811	IF(.NOT.tracers(iq)%isInPhysics) CYCLE
[4170]	2812	itr = itr+1
	2813	IF(tracers(iq)%component /= 'co2i') CYCLE
[3453]	2814	!--3D fields
[4170]	2815	CALL histwrite_phy(o_trac (itr),tr_seri(:,:,itr))
[4056]	2816	CALL histwrite_phy(o_dtr_vdf(itr),d_tr_cl(:,:,itr))
	2817	CALL histwrite_phy(o_dtr_the(itr),d_tr_th(:,:,itr))
	2818	CALL histwrite_phy(o_dtr_con(itr),d_tr_cv(:,:,itr))
[3453]	2819	!--2D fields
	2820	!--CO2 burden
	2821	zx_tmp_fi2d=0.
	2822	IF (vars_defined) THEN
	2823	DO k=1,klev
[4056]	2824	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,itr)
[3453]	2825	ENDDO
	2826	ENDIF
[4056]	2827	CALL histwrite_phy(o_trac_cum(itr), zx_tmp_fi2d)
[3550]	2828	ENDDO !--iq
	2829	!--CO2 net fluxes
	2830	CALL histwrite_phy(o_flx_co2_land, fco2_land)
	2831	CALL histwrite_phy(o_flx_co2_ocean, fco2_ocean)
[3857]	2832	CALL histwrite_phy(o_flx_co2_ocean_cor, fco2_ocean_cor)
	2833	CALL histwrite_phy(o_flx_co2_land_cor, fco2_land_cor)
[3550]	2834	CALL histwrite_phy(o_flx_co2_ff, fco2_ff)
	2835	CALL histwrite_phy(o_flx_co2_bb, fco2_bb)
[3125]	2836
[4389]	2837	ELSE IF (ANY(type_trac==['lmdz','coag'])) THEN
[4071]	2838	itr = 0
[4056]	2839	DO iq = 1, nqtot
[4071]	2840	IF(.NOT.tracers(iq)%isInPhysics) CYCLE
[4170]	2841	itr = itr + 1
	2842	! write(,) 'phys_output_write_mod 2337: itr=',itr
[3865]	2843	!--3D fields
[4170]	2844	CALL histwrite_phy(o_trac(itr), tr_seri(:,:,itr))
[4056]	2845	CALL histwrite_phy(o_dtr_vdf(itr),d_tr_cl(:,:,itr))
	2846	CALL histwrite_phy(o_dtr_the(itr),d_tr_th(:,:,itr))
	2847	CALL histwrite_phy(o_dtr_con(itr),d_tr_cv(:,:,itr))
[4170]	2848	CALL histwrite_phy(o_dtr_lessi_impa(itr),d_tr_lessi_impa(:,:,itr))
	2849	CALL histwrite_phy(o_dtr_lessi_nucl(itr),d_tr_lessi_nucl(:,:,itr))
	2850	CALL histwrite_phy(o_dtr_insc(itr),d_tr_insc(:,:,itr))
	2851	CALL histwrite_phy(o_dtr_bcscav(itr),d_tr_bcscav(:,:,itr))
	2852	CALL histwrite_phy(o_dtr_evapls(itr),d_tr_evapls(:,:,itr))
	2853	CALL histwrite_phy(o_dtr_ls(itr),d_tr_ls(:,:,itr))
	2854	CALL histwrite_phy(o_dtr_trsp(itr),d_tr_trsp(:,:,itr))
	2855	CALL histwrite_phy(o_dtr_sscav(itr),d_tr_sscav(:,:,itr))
	2856	CALL histwrite_phy(o_dtr_sat(itr),d_tr_sat(:,:,itr))
	2857	CALL histwrite_phy(o_dtr_uscav(itr),d_tr_uscav(:,:,itr))
	2858	!--2D fields
[5447]	2859	CALL histwrite_phy(o_dtr_wet_con(itr), flux_tr_wet(:,itr))
[4170]	2860	CALL histwrite_phy(o_dtr_dry(itr), flux_tr_dry(:,itr))
[3865]	2861	zx_tmp_fi2d=0.
	2862	IF (vars_defined) THEN
	2863	DO k=1,klev
[4056]	2864	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,itr)
[3865]	2865	ENDDO
	2866	ENDIF
[4056]	2867	CALL histwrite_phy(o_trac_cum(itr), zx_tmp_fi2d)
[3865]	2868	ENDDO !--iq
[4170]	2869	ENDIF !--type_trac
[2973]	2870	ENDIF !(iflag_phytrac==1)
[1813]	2871
[3815]	2872	if (activate_ocean_skin >= 1) then
	2873	CALL histwrite_phy(o_delta_sst, delta_sst)
	2874	CALL histwrite_phy(o_delta_sal, delta_sal)
	2875	CALL histwrite_phy(o_ds_ns, ds_ns)
	2876	CALL histwrite_phy(o_dt_ns, dt_ns)
	2877	CALL histwrite_phy(o_dter, dter)
	2878	CALL histwrite_phy(o_dser, dser)
	2879	CALL histwrite_phy(o_tkt, tkt)
	2880	CALL histwrite_phy(o_tks, tks)
	2881	CALL histwrite_phy(o_taur, taur)
	2882	CALL histwrite_phy(o_sss, sss)
	2883	end if
	2884
[5109]	2885	!! runoff land bucket - ajout S. Nguyen 23 07 2024
	2886	CALL histwrite_phy(o_runoff_diag, runoff_diag)
	2887
[4056]	2888	#ifdef ISO
[4984]	2889	!write(,) 'tmp phys_output_write: ntiso=',ntiso
[5022]	2890
[5109]	2891	DO ixt = 1, ntiso
[4984]	2892	!write(,) 'ixt,o_xtovap(ixt)=',ixt,o_xtovap(ixt)
[4056]	2893	IF (vars_defined) zx_tmp_fi2d(:) = xtrain_fall(ixt,:) + xtsnow_fall(ixt,:)
	2894	CALL histwrite_phy(o_xtprecip(ixt), zx_tmp_fi2d)
	2895
	2896	IF (vars_defined) zx_tmp_fi2d(:) = xtrain_lsc(ixt,:) + xtsnow_lsc(ixt,:)
	2897	CALL histwrite_phy(o_xtplul(ixt), zx_tmp_fi2d)
	2898
	2899	IF (vars_defined) zx_tmp_fi2d(:) = xtrain_con(ixt,:) + xtsnow_con(ixt,:)
	2900	CALL histwrite_phy(o_xtpluc(ixt), zx_tmp_fi2d)
	2901	CALL histwrite_phy(o_xtevap(ixt), xtevap(ixt,:))
	2902	CALL histwrite_phy(o_xtovap(ixt), xt_seri(ixt,:,:))
	2903	CALL histwrite_phy(o_xtoliq(ixt), xtl_seri(ixt,:,:))
	2904
[5022]	2905	!! runoff land bucket - ajout S. Nguyen 25 avril 2024
	2906	CALL histwrite_phy(o_xtrunoff_diag(ixt), xtrunoff_diag(ixt,:))
	2907
	2908
[4525]	2909	DO nsrf = 1, nbsrf ! ajout Camille 8 mai 2023
	2910	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fxtevap(ixt,:, nsrf)
	2911	CALL histwrite_phy(o_xtevap_srf(ixt,nsrf), zx_tmp_fi2d)
	2912	ENDDO
	2913
[4056]	2914	IF (vars_defined) zx_tmp_fi3d(:,:)=xtl_seri(ixt,:,:)+xts_seri(ixt,:,:)
	2915	CALL histwrite_phy(o_xtcond(ixt), zx_tmp_fi3d)
	2916	CALL histwrite_phy(o_dxtdyn(ixt), d_xt_dyn(ixt,:,:))
	2917	CALL histwrite_phy(o_dxtldyn(ixt), d_xtl_dyn(ixt,:,:))
	2918
	2919	IF (vars_defined) zx_tmp_fi3d(:,:)=d_xt_con(ixt,:,:)/pdtphys
	2920	CALL histwrite_phy(o_dxtcon(ixt), zx_tmp_fi3d)
	2921
	2922	IF (vars_defined) zx_tmp_fi3d(:,:)=d_xt_lsc(ixt,:,:)/pdtphys
	2923	CALL histwrite_phy(o_dxtlsc(ixt), zx_tmp_fi3d)
	2924
	2925	IF (vars_defined) zx_tmp_fi3d(:,:)=d_xt_eva(ixt,:,:)/pdtphys
	2926	CALL histwrite_phy(o_dxteva(ixt), zx_tmp_fi3d)
	2927
	2928	IF (vars_defined) zx_tmp_fi3d(:,:)=d_xt_vdf(ixt,:,:)/pdtphys
	2929	CALL histwrite_phy(o_dxtvdf(ixt), zx_tmp_fi3d)
	2930
	2931	IF (vars_defined) zx_tmp_fi3d(:,:)=d_xt_ajsb(ixt,:,:)/pdtphys
	2932	CALL histwrite_phy(o_dxtajs(ixt), zx_tmp_fi3d)
	2933
	2934	IF (vars_defined) zx_tmp_fi3d(:,:)=(d_xt_ajs(ixt,:,:)-d_xt_ajsb(ixt,:,:))/pdtphys
	2935	CALL histwrite_phy(o_dxtthe(ixt), zx_tmp_fi3d)
	2936
	2937	IF (ok_qch4) THEN
	2938	IF (vars_defined) zx_tmp_fi3d(:,:)=d_xt_ch4(ixt,:,:)/pdtphys
	2939	CALL histwrite_phy(o_dxtch4(ixt), zx_tmp_fi3d)
[4737]	2940	ENDIF
[4056]	2941
	2942	IF (ixt == iso_HTO) THEN
	2943	IF (vars_defined) zx_tmp_fi3d(:,:)=d_xt_prod_nucl(ixt,:,:)/pdtphys
	2944	CALL histwrite_phy(o_dxtprod_nucl(ixt), zx_tmp_fi3d)
	2945
	2946	IF (vars_defined) zx_tmp_fi3d(:,:)=d_xt_cosmo(ixt,:,:)/pdtphys
	2947	CALL histwrite_phy(o_dxtcosmo(ixt), zx_tmp_fi3d)
	2948
	2949	IF (vars_defined) zx_tmp_fi3d(:,:)=d_xt_decroiss(ixt,:,:)/pdtphys
	2950	CALL histwrite_phy(o_dxtdecroiss(ixt), zx_tmp_fi3d)
[4737]	2951	ENDIF
[4056]	2952
	2953	!write(,) 'phys_output_write_mod 2531'
[4737]	2954	ENDDO
[4056]	2955	#endif
	2956
[2690]	2957	IF (.NOT.vars_defined) THEN
[1898]	2958	!$OMP MASTER
[2097]	2959	#ifndef CPP_IOIPSL_NO_OUTPUT
[1898]	2960	DO iff=1,nfiles
	2961	IF (clef_files(iff)) THEN
	2962	CALL histend(nid_files(iff))
	2963	ndex2d = 0
	2964	ndex3d = 0
	2965	ENDIF ! clef_files
	2966	ENDDO ! iff
[1852]	2967	#endif
[5109]	2968
	2969	!SN activate water isotopes present in tracer.def
	2970	#ifdef ISO
	2971	DO ixt = 1, ntiso
	2972	outiso = TRIM(isoName(ixt))
	2973	i = INDEX(outiso, '_', .TRUE.)
	2974	outiso = outiso(1:i-1)//outiso(i+1:LEN_TRIM(outiso))
	2975
[5138]	2976	CALL xios_set_fieldgroup_attr("iso2D_"//TRIM(outiso), enabled=.TRUE.)
	2977	CALL xios_set_fieldgroup_attr("iso3D_"//TRIM(outiso), enabled=.TRUE.)
	2978
[5109]	2979	ENDDO
	2980	#endif
[1898]	2981	!On finalise l'initialisation:
[4619]	2982	IF (using_xios) CALL wxios_closedef()
	2983
[1898]	2984	!$OMP END MASTER
	2985	!$OMP BARRIER
	2986	vars_defined = .TRUE.
[1821]	2987
[2704]	2988	ENDIF !--.NOT.vars_defined
	2989
[2690]	2990	ENDDO
[1821]	2991
[2690]	2992	IF (vars_defined) THEN
[1898]	2993	! On synchronise les fichiers pour IOIPSL
[2097]	2994	#ifndef CPP_IOIPSL_NO_OUTPUT
[1898]	2995	!$OMP MASTER
	2996	DO iff=1,nfiles
[1798]	2997	IF (ok_sync .AND. clef_files(iff)) THEN
[1898]	2998	CALL histsync(nid_files(iff))
[1798]	2999	ENDIF
[4737]	3000	ENDDO
[1898]	3001	!$OMP END MASTER
[1852]	3002	#endif
[1798]	3003	ENDIF
	3004
[1898]	3005	END SUBROUTINE phys_output_write
[1798]	3006
[1898]	3007	END MODULE phys_output_write_mod

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

Download in other formats:

Original Format