Context Navigation

phys_output_write_mod.F90 @ 3780

Last change on this file since 3780 was 3780, checked in by evignon, 4 years ago
Premiere comission Etienne: changements pour le 1D (forcage en Ts au dessus des continents) et inclusion drag arbres dans yamada4_num=6
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: 100.6 KB

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

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

Download in other formats:

Original Format