Context Navigation

phys_output_write_mod.F90 @ 5209

Last change on this file since 5209 was 5204, checked in by Laurent Fairhead, 7 weeks ago
Integrating A.Borella's work on cirrus in the trunk
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.2 KB

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

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

Download in other formats:

Original Format