Context Navigation

phys_output_write_mod.F90 @ 3550

Last change on this file since 3550 was 3550, checked in by oboucher, 5 years ago
reordering loops for tracer output
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: 97.1 KB

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

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

Download in other formats:

Original Format