Context Navigation

phys_output_write_mod.F90 @ 3779

Last change on this file since 3779 was 3779, checked in by asima, 4 years ago
Code cleaning : cld*jn are removed as variables, because they only serve as output names.
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.0 KB

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

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

Download in other formats:

Original Format