Context Navigation

phys_output_write_mod.F90 @ 5227

Last change on this file since 5227 was 3483, checked in by acozic, 6 years ago
toplwad diagnostic works with INCA aeNP aerosols as well [ merge with rev 3381 of trunk ]
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: 91.5 KB

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

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

Download in other formats:

Original Format