source: LMDZ6/branches/blowing_snow/libf/phylmd/phys_output_write_mod.F90 @ 5327

Last change on this file since 5327 was 4485, checked in by evignon, 20 months ago

premier commit pour l'ajout de la neige soufflee sur la nouvelle branche

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