Context Navigation

phys_output_write_mod.F90 @ 5428

Last change on this file since 5428 was 3799, checked in by lguez, 4 years ago

Bug fix: no operation in arg of histwrite_phy

An operation like zx_rhl*100. in the actual argument of
histwrite_phy crashes the code if it is compiled in debug mode: at
the first call of phys_output_write_mod, zx_rhl is not defined.

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.8 KB

Line
1	!
2	! $Id: phys_output_write_mod.F90 3799 2021-01-11 23:03:09Z fhourdin $
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, landice_opt
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 (landice_opt .GE. 1) 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	IF (iflag_ice_thermo .GT. 0) THEN
1714	IF (vars_defined) zx_tmp_fi3d = zx_rhl * 100.
1715	CALL histwrite_phy(o_rhl, zx_tmp_fi3d)
1716	IF (vars_defined) zx_tmp_fi3d = zx_rhi * 100.
1717	CALL histwrite_phy(o_rhi, zx_tmp_fi3d)
1718	ENDIF
1719
1720
1721	IF (vars_defined) zx_tmp_fi3d = wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd
1722	CALL histwrite_phy(o_ozone, zx_tmp_fi3d)
1723
1724	IF (read_climoz == 2) THEN
1725	IF (vars_defined) zx_tmp_fi3d = wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd
1726	CALL histwrite_phy(o_ozone_light, zx_tmp_fi3d)
1727	ENDIF
1728
1729	CALL histwrite_phy(o_duphy, d_u)
1730
1731	CALL histwrite_phy(o_dtphy, d_t)
1732
1733	CALL histwrite_phy(o_dqphy, d_qx(:,:,ivap))
1734	IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,ivap),zmasse,zx_tmp_fi2d)
1735	CALL histwrite_phy(o_dqphy2d, zx_tmp_fi2d)
1736
1737	CALL histwrite_phy(o_dqlphy, d_qx(:,:,iliq))
1738	IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,iliq),zmasse,zx_tmp_fi2d)
1739	CALL histwrite_phy(o_dqlphy2d, zx_tmp_fi2d)
1740
1741	IF (nqo.EQ.3) THEN
1742	CALL histwrite_phy(o_dqsphy, d_qx(:,:,isol))
1743	IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,isol),zmasse,zx_tmp_fi2d)
1744	CALL histwrite_phy(o_dqsphy2d, zx_tmp_fi2d)
1745	ELSE
1746	zx_tmp_fi3d=0.0
1747	CALL histwrite_phy(o_dqsphy, zx_tmp_fi3d)
1748	zx_tmp_fi2d=0.0
1749	CALL histwrite_phy(o_dqsphy2d, zx_tmp_fi2d)
1750	ENDIF
1751
1752	DO nsrf=1, nbsrf
1753	IF (vars_defined) zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
1754	CALL histwrite_phy(o_albe_srf(nsrf), zx_tmp_fi2d)
1755	IF (vars_defined) zx_tmp_fi2d(1 : klon) = z0m( 1 : klon, nsrf)
1756	CALL histwrite_phy(o_z0m_srf(nsrf), zx_tmp_fi2d)
1757	IF (vars_defined) zx_tmp_fi2d(1 : klon) = z0h( 1 : klon, nsrf)
1758	CALL histwrite_phy(o_z0h_srf(nsrf), zx_tmp_fi2d)
1759	IF (vars_defined) zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1760	CALL histwrite_phy(o_ages_srf(nsrf), zx_tmp_fi2d)
1761	IF (vars_defined) zx_tmp_fi2d(1 : klon) = snow( 1 : klon, nsrf)
1762	CALL histwrite_phy(o_snow_srf(nsrf), zx_tmp_fi2d)
1763	ENDDO !nsrf=1, nbsrf
1764	CALL histwrite_phy(o_alb1, albsol1)
1765	CALL histwrite_phy(o_alb2, albsol2)
1766	!FH Sorties pour la couche limite
1767	IF (iflag_pbl>1) THEN
1768	zx_tmp_fi3d=0.
1769	IF (vars_defined) THEN
1770	DO nsrf=1,nbsrf
1771	DO k=1,klev
1772	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) &
1773	+pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
1774	ENDDO
1775	ENDDO
1776	ENDIF
1777	CALL histwrite_phy(o_tke, zx_tmp_fi3d)
1778	CALL histwrite_phy(o_tke_max, zx_tmp_fi3d)
1779
1780	ENDIF
1781
1782	CALL histwrite_phy(o_kz, coefh(:,:,is_ave))
1783
1784	CALL histwrite_phy(o_kz_max, coefh(:,:,is_ave))
1785
1786	CALL histwrite_phy(o_clwcon, clwcon0)
1787	CALL histwrite_phy(o_dtdyn, d_t_dyn)
1788
1789	CALL histwrite_phy(o_dqdyn, d_q_dyn)
1790
1791	CALL histwrite_phy(o_dqdyn2d,d_q_dyn2d)
1792
1793	CALL histwrite_phy(o_dqldyn, d_ql_dyn)
1794
1795	CALL histwrite_phy(o_dqldyn2d, d_ql_dyn2d)
1796
1797	CALL histwrite_phy(o_dqsdyn, d_qs_dyn)
1798
1799	CALL histwrite_phy(o_dqsdyn2d, d_qs_dyn2d)
1800
1801	CALL histwrite_phy(o_dudyn, d_u_dyn)
1802	CALL histwrite_phy(o_dvdyn, d_v_dyn)
1803
1804	IF (vars_defined) THEN
1805	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
1806	ENDIF
1807	CALL histwrite_phy(o_dtcon, zx_tmp_fi3d)
1808	IF (iflag_thermals.EQ.0) THEN
1809	IF (vars_defined) THEN
1810	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
1811	d_t_ajsb(1:klon,1:klev)/pdtphys
1812	ENDIF
1813	CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
1814	ELSE IF(iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
1815	IF (vars_defined) THEN
1816	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
1817	d_t_ajs(1:klon,1:klev)/pdtphys + &
1818	d_t_wake(1:klon,1:klev)/pdtphys
1819	ENDIF
1820	CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
1821	ENDIF
1822	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1823	CALL histwrite_phy(o_ducon, zx_tmp_fi3d)
1824	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
1825	CALL histwrite_phy(o_dvcon, zx_tmp_fi3d)
1826	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1827	CALL histwrite_phy(o_dqcon, zx_tmp_fi3d)
1828	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1829	CALL histwrite_phy(o_dqcon2d, zx_tmp_fi2d)
1830
1831	IF (iflag_thermals.EQ.0) THEN
1832	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1833	CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d)
1834	ELSE IF (iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
1835	IF (vars_defined) THEN
1836	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys + &
1837	d_q_ajs(1:klon,1:klev)/pdtphys + &
1838	d_q_wake(1:klon,1:klev)/pdtphys
1839	ENDIF
1840	CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d)
1841	ENDIF
1842
1843	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
1844	CALL histwrite_phy(o_dtlsc, zx_tmp_fi3d)
1845	IF (vars_defined) zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+ &
1846	d_t_eva(1:klon,1:klev))/pdtphys
1847	CALL histwrite_phy(o_dtlschr, zx_tmp_fi3d)
1848	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1849	CALL histwrite_phy(o_dqlsc, zx_tmp_fi3d)
1850	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1851	CALL histwrite_phy(o_dqlsc2d, zx_tmp_fi2d)
1852	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=beta_prec(1:klon,1:klev)
1853	CALL histwrite_phy(o_beta_prec, zx_tmp_fi3d)
1854	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1855	! Sorties specifiques a la separation thermiques/non thermiques
1856	IF (iflag_thermals>=1) THEN
1857	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys
1858	CALL histwrite_phy(o_dtlscth, zx_tmp_fi3d)
1859	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys
1860	CALL histwrite_phy(o_dtlscst, zx_tmp_fi3d)
1861	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys
1862	CALL histwrite_phy(o_dqlscth, zx_tmp_fi3d)
1863	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1864	CALL histwrite_phy(o_dqlscth2d, zx_tmp_fi2d)
1865	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys
1866	CALL histwrite_phy(o_dqlscst, zx_tmp_fi3d)
1867	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1868	CALL histwrite_phy(o_dqlscst2d, zx_tmp_fi2d)
1869	CALL histwrite_phy(o_plulth, plul_th)
1870	CALL histwrite_phy(o_plulst, plul_st)
1871	IF (vars_defined) THEN
1872	DO i=1,klon
1873	zx_tmp_fi2d(1:klon)=lmax_th(:)
1874	ENDDO
1875	ENDIF
1876	CALL histwrite_phy(o_lmaxth, zx_tmp_fi2d)
1877	IF (vars_defined) THEN
1878	DO k=1,klev
1879	DO i=1,klon
1880	IF (ptconvth(i,k)) THEN
1881	zx_tmp_fi3d(i,k)=1.
1882	ELSE
1883	zx_tmp_fi3d(i,k)=0.
1884	ENDIF
1885	ENDDO
1886	ENDDO
1887	ENDIF
1888	CALL histwrite_phy(o_ptconvth, zx_tmp_fi3d)
1889	ENDIF ! iflag_thermals>=1
1890	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1891	zpt_conv = 0.
1892	WHERE (ptconv) zpt_conv = 1.
1893	CALL histwrite_phy(o_ptconv, zpt_conv)
1894	!! IF (vars_defined) zx_tmp_fi2d=float(itau_con)/float(itap)
1895	!! CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
1896	IF (vars_defined) THEN
1897	zpt_conv2d(:) = 0.
1898	DO k=1,klev
1899	WHERE (ptconv(:,k)) zpt_conv2d(:) = 1.
1900	ENDDO
1901	ENDIF
1902	CALL histwrite_phy(o_ftime_deepcv, zpt_conv2d)
1903	IF (vars_defined) THEN
1904	zx_tmp_fi2d(:) = 0.
1905	DO k=1,klev
1906	WHERE (ptconvth(:,k)) zx_tmp_fi2d(:) = 1.
1907	ENDDO
1908	ENDIF
1909	CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
1910	IF (vars_defined) THEN
1911	zx_tmp_fi2d(:) = max(zx_tmp_fi2d(:),zpt_conv2d(:))
1912	ENDIF
1913	CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
1914	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1915	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1916	CALL histwrite_phy(o_dtvdf, zx_tmp_fi3d)
1917	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_diss(1:klon,1:klev)/pdtphys
1918	CALL histwrite_phy(o_dtdis, zx_tmp_fi3d)
1919	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1920	CALL histwrite_phy(o_dqvdf, zx_tmp_fi3d)
1921	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1922	CALL histwrite_phy(o_dqvdf2d, zx_tmp_fi2d)
1923	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1924	CALL histwrite_phy(o_dteva, zx_tmp_fi3d)
1925	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1926	CALL histwrite_phy(o_dqeva, zx_tmp_fi3d)
1927	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1928	CALL histwrite_phy(o_dqeva2d, zx_tmp_fi2d)
1929	CALL histwrite_phy(o_ratqs, ratqs)
1930	IF (vars_defined) THEN
1931	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys - &
1932	d_t_ajsb(1:klon,1:klev)/pdtphys
1933	ENDIF
1934	CALL histwrite_phy(o_dtthe, zx_tmp_fi3d)
1935	IF (vars_defined) THEN
1936	zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys
1937	ENDIF
1938	CALL histwrite_phy(o_duthe, zx_tmp_fi3d)
1939	IF (vars_defined) THEN
1940	zx_tmp_fi3d(1:klon,1:klev)=d_v_ajs(1:klon,1:klev)/pdtphys
1941	ENDIF
1942	CALL histwrite_phy(o_dvthe, zx_tmp_fi3d)
1943
1944	IF (iflag_thermals>=1) THEN
1945	! Pour l instant 0 a y reflichir pour les thermiques
1946	! regroupe avec ftime_deepcv et ftime_con
1947	!!zx_tmp_fi2d=0.
1948	!!CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
1949	CALL histwrite_phy(o_f_th, fm_therm)
1950	CALL histwrite_phy(o_e_th, entr_therm)
1951	CALL histwrite_phy(o_w_th, zw2)
1952	CALL histwrite_phy(o_q_th, zqasc)
1953	CALL histwrite_phy(o_a_th, fraca)
1954	CALL histwrite_phy(o_cloudth_sth, cloudth_sth)
1955	CALL histwrite_phy(o_cloudth_senv, cloudth_senv)
1956	CALL histwrite_phy(o_cloudth_sigmath, cloudth_sigmath)
1957	CALL histwrite_phy(o_cloudth_sigmaenv, cloudth_sigmaenv)
1958	CALL histwrite_phy(o_d_th, detr_therm)
1959	CALL histwrite_phy(o_f0_th, f0)
1960	CALL histwrite_phy(o_zmax_th, zmax_th)
1961	IF (vars_defined) THEN
1962	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys - &
1963	d_q_ajsb(1:klon,1:klev)/pdtphys
1964	ENDIF
1965	CALL histwrite_phy(o_dqthe, zx_tmp_fi3d)
1966	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1967	CALL histwrite_phy(o_dqthe2d, zx_tmp_fi2d)
1968	ENDIF !iflag_thermals
1969	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
1970	CALL histwrite_phy(o_dtajs, zx_tmp_fi3d)
1971	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
1972	CALL histwrite_phy(o_dqajs, zx_tmp_fi3d)
1973	IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1974	CALL histwrite_phy(o_dqajs2d, zx_tmp_fi2d)
1975	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys
1976	CALL histwrite_phy(o_dtswr, zx_tmp_fi3d)
1977	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_sw0(1:klon,1:klev)/pdtphys
1978	CALL histwrite_phy(o_dtsw0, zx_tmp_fi3d)
1979	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lwr(1:klon,1:klev)/pdtphys
1980	CALL histwrite_phy(o_dtlwr, zx_tmp_fi3d)
1981	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lw0(1:klon,1:klev)/pdtphys
1982	CALL histwrite_phy(o_dtlw0, zx_tmp_fi3d)
1983	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)/pdtphys
1984	CALL histwrite_phy(o_dtec, zx_tmp_fi3d)
1985	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1986	CALL histwrite_phy(o_duvdf, zx_tmp_fi3d)
1987	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1988	CALL histwrite_phy(o_dvvdf, zx_tmp_fi3d)
1989	IF (ok_orodr) THEN
1990	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1991	CALL histwrite_phy(o_duoro, zx_tmp_fi3d)
1992	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1993	CALL histwrite_phy(o_dvoro, zx_tmp_fi3d)
1994	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys
1995	CALL histwrite_phy(o_dtoro, zx_tmp_fi3d)
1996	ENDIF
1997	IF (ok_orolf) THEN
1998	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
1999	CALL histwrite_phy(o_dulif, zx_tmp_fi3d)
2000
2001	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
2002	CALL histwrite_phy(o_dvlif, zx_tmp_fi3d)
2003
2004	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys
2005	CALL histwrite_phy(o_dtlif, zx_tmp_fi3d)
2006	ENDIF
2007
2008	IF (ok_hines) THEN
2009	IF (vars_defined) zx_tmp_fi3d=du_gwd_hines/pdtphys
2010	CALL histwrite_phy(o_du_gwd_hines, zx_tmp_fi3d)
2011
2012	IF (vars_defined) zx_tmp_fi3d= dv_gwd_hines/pdtphys
2013	CALL histwrite_phy(o_dv_gwd_hines, zx_tmp_fi3d)
2014
2015	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys
2016	CALL histwrite_phy(o_dthin, zx_tmp_fi3d)
2017	CALL histwrite_phy(o_ustr_gwd_hines, zustr_gwd_hines)
2018	CALL histwrite_phy(o_vstr_gwd_hines, zvstr_gwd_hines)
2019	ENDIF
2020
2021	IF (.not. ok_hines .and. ok_gwd_rando) THEN
2022	IF (vars_defined) zx_tmp_fi3d=du_gwd_front / pdtphys
2023	CALL histwrite_phy(o_du_gwd_front, zx_tmp_fi3d)
2024
2025	IF (vars_defined) zx_tmp_fi3d=dv_gwd_front / pdtphys
2026	CALL histwrite_phy(o_dv_gwd_front, zx_tmp_fi3d)
2027
2028	CALL histwrite_phy(o_ustr_gwd_front, zustr_gwd_front)
2029	CALL histwrite_phy(o_vstr_gwd_front, zvstr_gwd_front)
2030	ENDIF
2031
2032	IF (ok_gwd_rando) THEN
2033	IF (vars_defined) zx_tmp_fi3d=du_gwd_rando / pdtphys
2034	CALL histwrite_phy(o_du_gwd_rando, zx_tmp_fi3d)
2035
2036	IF (vars_defined) zx_tmp_fi3d=dv_gwd_rando / pdtphys
2037	CALL histwrite_phy(o_dv_gwd_rando, zx_tmp_fi3d)
2038	CALL histwrite_phy(o_ustr_gwd_rando, zustr_gwd_rando)
2039	CALL histwrite_phy(o_vstr_gwd_rando, zvstr_gwd_rando)
2040	CALL histwrite_phy(o_east_gwstress, east_gwstress )
2041	CALL histwrite_phy(o_west_gwstress, west_gwstress )
2042	ENDIF
2043
2044	IF (ok_qch4) THEN
2045	IF (vars_defined) zx_tmp_fi3d=d_q_ch4 / pdtphys
2046	CALL histwrite_phy(o_dqch4, zx_tmp_fi3d)
2047	ENDIF
2048
2049	IF (vars_defined) THEN
2050	DO k=1, klevp1
2051	zx_tmp_fi3d1(:,k)=swup(:,k)*swradcorr(:)
2052	ENDDO
2053	ENDIF
2054
2055	CALL histwrite_phy(o_rsu, zx_tmp_fi3d1)
2056
2057
2058	IF (vars_defined) THEN
2059	DO k=1, klevp1
2060	zx_tmp_fi3d1(:,k)=swdn(:,k)*swradcorr(:)
2061	ENDDO
2062	ENDIF
2063
2064	CALL histwrite_phy(o_rsd, zx_tmp_fi3d1)
2065
2066	IF (vars_defined) THEN
2067	DO k=1, klevp1
2068	zx_tmp_fi3d1(:,k)=swup0(:,k)*swradcorr(:)
2069	ENDDO
2070	ENDIF
2071
2072	CALL histwrite_phy(o_rsucs, zx_tmp_fi3d1)
2073
2074	IF (vars_defined) THEN
2075	DO k=1, klevp1
2076	zx_tmp_fi3d1(:,k)=swupc0(:,k)*swradcorr(:)
2077	ENDDO
2078	ENDIF
2079	CALL histwrite_phy(o_rsucsaf, zx_tmp_fi3d1)
2080
2081	IF (vars_defined) THEN
2082	DO k=1, klevp1
2083	zx_tmp_fi3d1(:,k)=swdn0(:,k)*swradcorr(:)
2084	ENDDO
2085	ENDIF
2086	CALL histwrite_phy(o_rsdcs, zx_tmp_fi3d1)
2087
2088
2089	IF (vars_defined) THEN
2090	DO k=1, klevp1
2091	zx_tmp_fi3d1(:,k)=swdnc0(:,k)*swradcorr(:)
2092	ENDDO
2093	ENDIF
2094	CALL histwrite_phy(o_rsdcsaf, zx_tmp_fi3d1)
2095
2096	CALL histwrite_phy(o_rlu, lwup)
2097	CALL histwrite_phy(o_rld, lwdn)
2098	CALL histwrite_phy(o_rlucs, lwup0)
2099	CALL histwrite_phy(o_rldcs, lwdn0)
2100
2101	IF (vars_defined) THEN
2102	zx_tmp_fi3d(1:klon,1:klev)=d_t(1:klon,1:klev)+ &
2103	d_t_dyn(1:klon,1:klev)
2104	ENDIF
2105	CALL histwrite_phy(o_tnt, zx_tmp_fi3d)
2106
2107	IF (vars_defined) THEN
2108	zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys + &
2109	d_t_lwr(1:klon,1:klev)/pdtphys
2110	ENDIF
2111	CALL histwrite_phy(o_tntr, zx_tmp_fi3d)
2112	IF (vars_defined) THEN
2113	zx_tmp_fi3d(1:klon,1:klev)= (d_t_lsc(1:klon,1:klev)+ &
2114	d_t_eva(1:klon,1:klev)+ &
2115	d_t_vdf(1:klon,1:klev))/pdtphys
2116	ENDIF
2117	CALL histwrite_phy(o_tntscpbl, zx_tmp_fi3d)
2118	IF (vars_defined) THEN
2119	zx_tmp_fi3d(1:klon,1:klev)=d_qx(1:klon,1:klev,ivap)+ &
2120	d_q_dyn(1:klon,1:klev)
2121	ENDIF
2122	CALL histwrite_phy(o_tnhus, zx_tmp_fi3d)
2123	IF (vars_defined) THEN
2124	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys+ &
2125	d_q_eva(1:klon,1:klev)/pdtphys
2126	ENDIF
2127	CALL histwrite_phy(o_tnhusscpbl, zx_tmp_fi3d)
2128	CALL histwrite_phy(o_evu, coefm(:,:,is_ave))
2129	IF (vars_defined) THEN
2130	zx_tmp_fi3d(1:klon,1:klev)=q_seri(1:klon,1:klev)+ &
2131	ql_seri(1:klon,1:klev)
2132	ENDIF
2133	CALL histwrite_phy(o_h2o, zx_tmp_fi3d)
2134	IF (iflag_con >= 3) THEN
2135	IF (vars_defined) THEN
2136	zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+ &
2137	dnwd0(1:klon,1:klev))
2138	ENDIF
2139	CALL histwrite_phy(o_mcd, zx_tmp_fi3d)
2140	IF (vars_defined) THEN
2141	zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) + &
2142	dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev)
2143	ENDIF
2144	CALL histwrite_phy(o_dmc, zx_tmp_fi3d)
2145	ELSE IF (iflag_con == 2) THEN
2146	CALL histwrite_phy(o_mcd, pmfd)
2147	IF (vars_defined) zx_tmp_fi3d = pmfu + pmfd
2148	CALL histwrite_phy(o_dmc, zx_tmp_fi3d)
2149	ENDIF
2150	CALL histwrite_phy(o_ref_liq, ref_liq)
2151	CALL histwrite_phy(o_ref_ice, ref_ice)
2152	!
2153	IF (ok_4xCO2atm) THEN
2154	IF (vars_defined) zx_tmp_fi2d(:) = swupp(:,klevp1)*swradcorr(:)
2155	CALL histwrite_phy(o_rsut4co2, zx_tmp_fi2d)
2156	IF (vars_defined) zx_tmp_fi2d(:) = lwupp(:,klevp1)
2157	CALL histwrite_phy(o_rlut4co2, zx_tmp_fi2d)
2158	IF (vars_defined) zx_tmp_fi2d(:) = swup0p(:,klevp1)*swradcorr(:)
2159	CALL histwrite_phy(o_rsutcs4co2, zx_tmp_fi2d)
2160	IF (vars_defined) zx_tmp_fi2d(:) = lwup0p(:,klevp1)
2161	CALL histwrite_phy(o_rlutcs4co2, zx_tmp_fi2d)
2162	IF (vars_defined) THEN
2163	DO k=1, klevp1
2164	zx_tmp_fi3d1(:,k)=swupp(:,k)*swradcorr(:)
2165	ENDDO
2166	ENDIF
2167	CALL histwrite_phy(o_rsu4co2, zx_tmp_fi3d1)
2168	IF (vars_defined) THEN
2169	DO k=1, klevp1
2170	zx_tmp_fi3d1(:,k)=swup0p(:,k)*swradcorr(:)
2171	ENDDO
2172	ENDIF
2173	CALL histwrite_phy(o_rsucs4co2, zx_tmp_fi3d1)
2174	IF (vars_defined) THEN
2175	DO k=1, klevp1
2176	zx_tmp_fi3d1(:,k)=swdnp(:,k)*swradcorr(:)
2177	ENDDO
2178	ENDIF
2179	CALL histwrite_phy(o_rsd4co2, zx_tmp_fi3d1)
2180	IF (vars_defined) THEN
2181	DO k=1, klevp1
2182	zx_tmp_fi3d1(:,k)=swdn0p(:,k)*swradcorr(:)
2183	ENDDO
2184	ENDIF
2185	CALL histwrite_phy(o_rsdcs4co2, zx_tmp_fi3d1)
2186	CALL histwrite_phy(o_rlu4co2, lwupp)
2187	CALL histwrite_phy(o_rlucs4co2, lwup0p)
2188	CALL histwrite_phy(o_rld4co2, lwdnp)
2189	CALL histwrite_phy(o_rldcs4co2, lwdn0p)
2190	ENDIF !ok_4xCO2atm
2191	!!!!!!!!!!!! Sorties niveaux de pression NMC !!!!!!!!!!!!!!!!!!!!
2192	#ifdef CPP_IOIPSL
2193	#ifndef CPP_XIOS
2194	IF (.NOT.ok_all_xml) THEN
2195	! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX:
2196	! Champs interpolles sur des niveaux de pression
2197	DO iff=7, nfiles-1 !--OB: here we deal with files 7,8,9
2198
2199	CALL histwrite_phy(o_tnondef,tnondef(:,:,iff-6),iff)
2200	CALL histwrite_phy(o_ta,twriteSTD(:,:,iff-6),iff)
2201	CALL histwrite_phy(o_zg,phiwriteSTD(:,:,iff-6),iff)
2202	CALL histwrite_phy(o_hus,qwriteSTD(:,:,iff-6),iff)
2203	CALL histwrite_phy(o_hur,rhwriteSTD(:,:,iff-6),iff)
2204	CALL histwrite_phy(o_ua,uwriteSTD(:,:,iff-6),iff)
2205	CALL histwrite_phy(o_va,vwriteSTD(:,:,iff-6),iff)
2206	CALL histwrite_phy(o_wap,wwriteSTD(:,:,iff-6),iff)
2207	IF (vars_defined) THEN
2208	DO k=1, nlevSTD
2209	DO i=1, klon
2210	IF (tnondef(i,k,iff-6).NE.missing_val) THEN
2211	IF (freq_outNMC(iff-6).LT.0) THEN
2212	freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6)
2213	ELSE
2214	freq_moyNMC(iff-6)=freq_outNMC(iff-6)/freq_calNMC(iff-6)
2215	ENDIF
2216	zx_tmp_fi3d_STD(i,k) = (100.*tnondef(i,k,iff-6))/freq_moyNMC(iff-6)
2217	ELSE
2218	zx_tmp_fi3d_STD(i,k) = missing_val
2219	ENDIF
2220	ENDDO
2221	ENDDO
2222	ENDIF
2223	CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD,iff)
2224	IF (vars_defined) THEN
2225	DO k=1, nlevSTD
2226	DO i=1, klon
2227	IF (O3sumSTD(i,k,iff-6).NE.missing_val) THEN
2228	zx_tmp_fi3d_STD(i,k) = O3sumSTD(i,k,iff-6) * 1.e+9
2229	ELSE
2230	zx_tmp_fi3d_STD(i,k) = missing_val
2231	ENDIF
2232	ENDDO
2233	ENDDO !k=1, nlevSTD
2234	ENDIF
2235	CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD,iff)
2236	IF (read_climoz == 2) THEN
2237	IF (vars_defined) THEN
2238	DO k=1, nlevSTD
2239	DO i=1, klon
2240	IF (O3daysumSTD(i,k,iff-6).NE.missing_val) THEN
2241	zx_tmp_fi3d_STD(i,k) = O3daysumSTD(i,k,iff-6) * 1.e+9
2242	ELSE
2243	zx_tmp_fi3d_STD(i,k) = missing_val
2244	ENDIF
2245	ENDDO
2246	ENDDO !k=1, nlevSTD
2247	ENDIF
2248	CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d_STD,iff)
2249	endif
2250	CALL histwrite_phy(o_uxv,uvsumSTD(:,:,iff-6),iff)
2251	CALL histwrite_phy(o_vxq,vqsumSTD(:,:,iff-6),iff)
2252	CALL histwrite_phy(o_vxT,vTsumSTD(:,:,iff-6),iff)
2253	CALL histwrite_phy(o_wxq,wqsumSTD(:,:,iff-6),iff)
2254	CALL histwrite_phy(o_vxphi,vphisumSTD(:,:,iff-6),iff)
2255	CALL histwrite_phy(o_wxT,wTsumSTD(:,:,iff-6),iff)
2256	CALL histwrite_phy(o_uxu,u2sumSTD(:,:,iff-6),iff)
2257	CALL histwrite_phy(o_vxv,v2sumSTD(:,:,iff-6),iff)
2258	CALL histwrite_phy(o_TxT,T2sumSTD(:,:,iff-6),iff)
2259	ENDDO !nfiles
2260	ENDIF
2261	#endif
2262	#endif
2263	#ifdef CPP_XIOS
2264	IF (ok_all_xml) THEN
2265	! DO iff=7, nfiles
2266
2267	! CALL histwrite_phy(o_tnondef,tnondef(:,:,3))
2268	CALL histwrite_phy(o_ta,tlevSTD(:,:))
2269	CALL histwrite_phy(o_zg,philevSTD(:,:))
2270	CALL histwrite_phy(o_hus,qlevSTD(:,:))
2271	CALL histwrite_phy(o_hur,rhlevSTD(:,:))
2272	CALL histwrite_phy(o_ua,ulevSTD(:,:))
2273	CALL histwrite_phy(o_va,vlevSTD(:,:))
2274	CALL histwrite_phy(o_wap,wlevSTD(:,:))
2275	! IF (vars_defined) THEN
2276	! DO k=1, nlevSTD
2277	! DO i=1, klon
2278	! IF (tnondef(i,k,3).NE.missing_val) THEN
2279	! IF (freq_outNMC(iff-6).LT.0) THEN
2280	! freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6)
2281	! ELSE
2282	! freq_moyNMC(iff-6)=freq_outNMC(iff-6)/freq_calNMC(iff-6)
2283	! ENDIF
2284	! zx_tmp_fi3d_STD(i,k) = (100.*tnondef(i,k,3))/freq_moyNMC(iff-6)
2285	! ELSE
2286	! zx_tmp_fi3d_STD(i,k) = missing_val
2287	! ENDIF
2288	! ENDDO
2289	! ENDDO
2290	! ENDIF
2291	! CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD)
2292	IF (vars_defined) THEN
2293	DO k=1, nlevSTD
2294	DO i=1, klon
2295	IF (O3STD(i,k).NE.missing_val) THEN
2296	zx_tmp_fi3d_STD(i,k) = O3STD(i,k) * 1.e+9
2297	ELSE
2298	zx_tmp_fi3d_STD(i,k) = missing_val
2299	ENDIF
2300	ENDDO
2301	ENDDO !k=1, nlevSTD
2302	ENDIF
2303	CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD)
2304	IF (read_climoz == 2) THEN
2305	IF (vars_defined) THEN
2306	DO k=1, nlevSTD
2307	DO i=1, klon
2308	IF (O3daySTD(i,k).NE.missing_val) THEN
2309	zx_tmp_fi3d_STD(i,k) = O3daySTD(i,k) * 1.e+9
2310	ELSE
2311	zx_tmp_fi3d_STD(i,k) = missing_val
2312	ENDIF
2313	ENDDO
2314	ENDDO !k=1, nlevSTD
2315	ENDIF
2316	CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d_STD)
2317	ENDIF
2318	CALL histwrite_phy(o_uxv,uvSTD(:,:))
2319	CALL histwrite_phy(o_vxq,vqSTD(:,:))
2320	CALL histwrite_phy(o_vxT,vTSTD(:,:))
2321	CALL histwrite_phy(o_wxq,wqSTD(:,:))
2322	CALL histwrite_phy(o_vxphi,vphiSTD(:,:))
2323	CALL histwrite_phy(o_wxT,wTSTD(:,:))
2324	CALL histwrite_phy(o_uxu,u2STD(:,:))
2325	CALL histwrite_phy(o_vxv,v2STD(:,:))
2326	CALL histwrite_phy(o_TxT,T2STD(:,:))
2327	! ENDDO !nfiles
2328	ENDIF
2329	#endif
2330	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2331	IF (iflag_phytrac == 1 ) then
2332	IF (type_trac == 'lmdz' .OR. type_trac == 'coag') THEN
2333	DO iq=nqo+1, nqtot
2334	!--3D fields
2335	CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
2336	CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
2337	CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
2338	CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo))
2339	CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))
2340	CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))
2341	CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo))
2342	CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo))
2343	CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo))
2344	CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))
2345	CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))
2346	CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))
2347	CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))
2348	CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo))
2349	!--2D fields
2350	CALL histwrite_phy(o_dtr_dry(iq-nqo), flux_tr_dry(:,iq-nqo))
2351	zx_tmp_fi2d=0.
2352	IF (vars_defined) THEN
2353	DO k=1,klev
2354	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,iq-nqo)
2355	ENDDO
2356	ENDIF
2357	CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
2358	ENDDO !--iq
2359	ENDIF !--type_trac
2360	!
2361	IF (type_trac == 'co2i') THEN
2362	DO iq=nqo+1, nqtot
2363	!--3D fields
2364	CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
2365	CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
2366	CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
2367	CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo))
2368	!--2D fields
2369	!--CO2 burden
2370	zx_tmp_fi2d=0.
2371	IF (vars_defined) THEN
2372	DO k=1,klev
2373	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,iq-nqo)
2374	ENDDO
2375	ENDIF
2376	CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
2377	ENDDO !--iq
2378	!--CO2 net fluxes
2379	CALL histwrite_phy(o_flx_co2_land, fco2_land)
2380	CALL histwrite_phy(o_flx_co2_ocean, fco2_ocean)
2381	CALL histwrite_phy(o_flx_co2_ff, fco2_ff)
2382	CALL histwrite_phy(o_flx_co2_bb, fco2_bb)
2383	ENDIF !--type_trac co2i
2384
2385	ENDIF !(iflag_phytrac==1)
2386
2387	IF (.NOT.vars_defined) THEN
2388	!$OMP MASTER
2389	#ifndef CPP_IOIPSL_NO_OUTPUT
2390	DO iff=1,nfiles
2391	IF (clef_files(iff)) THEN
2392	CALL histend(nid_files(iff))
2393	ndex2d = 0
2394	ndex3d = 0
2395	ENDIF ! clef_files
2396	ENDDO ! iff
2397	#endif
2398	#ifdef CPP_XIOS
2399	!On finalise l'initialisation:
2400	CALL wxios_closedef()
2401	#endif
2402	!$OMP END MASTER
2403	!$OMP BARRIER
2404	vars_defined = .TRUE.
2405
2406	ENDIF !--.NOT.vars_defined
2407
2408	ENDDO
2409
2410	IF (vars_defined) THEN
2411	! On synchronise les fichiers pour IOIPSL
2412	#ifndef CPP_IOIPSL_NO_OUTPUT
2413	!$OMP MASTER
2414	DO iff=1,nfiles
2415	IF (ok_sync .AND. clef_files(iff)) THEN
2416	CALL histsync(nid_files(iff))
2417	ENDIF
2418	END DO
2419	!$OMP END MASTER
2420	#endif
2421	ENDIF
2422
2423	END SUBROUTINE phys_output_write
2424
2425	END MODULE phys_output_write_mod

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

Download in other formats:

Original Format