Context Navigation

phys_output_write_mod.F90 @ 2580

Last change on this file since 2580 was 2580, checked in by idelkadi, 8 years ago
Implementation du simulateur AIRS: Le but du simulateur est de permettre la comparaison de proprietes macro-et microphysiques des nuages hauts de LMDZ avec celles restituees par les observations du satellite AIRS (Atmospheric IR Sounder). La methode est decrite dans Hendricks et al. Meteorol. Z., 2010
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: 73.8 KB

Line
1	!
2	! $Id: phys_output_write_mod.F90 2580 2016-06-29 12:23:19Z idelkadi $
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
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 (los 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, ivap, iliq, isol, new_aod, ok_sync, &
20	ptconv, read_climoz, clevSTD, ptconvth, &
21	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, nslay
27	USE mod_phys_lmdz_para, ONLY: is_north_pole_phy,is_south_pole_phy
28	USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
29	USE time_phylmdz_mod, only: day_step_phy, start_time, itau_phy
30	USE phys_output_ctrlout_mod, only: o_phis, o_aire, is_ter, is_lic, is_oce, &
31	is_ave, is_sic, o_contfracATM, o_contfracOR, &
32	o_aireTER, o_flat, o_slp, o_ptstar, o_pt0, o_tsol, &
33	o_t2m, o_t2m_min, o_t2m_max, &
34	o_t2m_min_mon, o_t2m_max_mon, &
35	o_q2m, o_ustar, o_u10m, o_v10m, &
36	o_wind10m, o_wind10max, o_gusts, o_sicf, &
37	o_psol, o_mass, o_qsurf, o_qsol, &
38	o_precip, o_ndayrain, o_plul, o_pluc, o_plun, &
39	o_snow, o_msnow, o_fsnow, o_evap, o_ep,o_epmax_diag, & ! epmax_cape
40	o_tops, o_tops0, o_topl, o_topl0, &
41	o_SWupTOA, o_SWupTOAclr, o_SWdnTOA, &
42	o_SWdnTOAclr, o_nettop, o_SWup200, &
43	o_SWup200clr, o_SWdn200, o_SWdn200clr, &
44	o_LWup200, o_LWup200clr, o_LWdn200, &
45	o_LWdn200clr, o_sols, o_sols0, &
46	o_soll, o_radsol, o_soll0, o_SWupSFC, &
47	o_SWupSFCclr, o_SWdnSFC, o_SWdnSFCclr, &
48	o_LWupSFC, o_LWdnSFC, o_LWupSFCclr, &
49	o_LWdnSFCclr, o_bils, o_bils_diss, &
50	o_bils_ec,o_bils_ech, o_bils_tke, o_bils_kinetic, &
51	o_bils_latent, o_bils_enthalp, o_sens, &
52	o_fder, o_ffonte, o_fqcalving, o_fqfonte, o_runofflic, &
53	o_taux, o_tauy, o_snowsrf, o_qsnow, &
54	o_snowhgt, o_toice, o_sissnow, o_runoff, &
55	o_albslw3, o_pourc_srf, o_fract_srf, &
56	o_taux_srf, o_tauy_srf, o_tsol_srf, &
57	o_evappot_srf, o_ustar_srf, o_u10m_srf, &
58	o_v10m_srf, o_t2m_srf, o_evap_srf, &
59	o_sens_srf, o_lat_srf, o_flw_srf, &
60	o_fsw_srf, o_wbils_srf, o_wbilo_srf, &
61	o_tke_srf, o_tke_max_srf,o_dltpbltke_srf, o_wstar, &
62	o_l_mixmin,o_l_mix, &
63	o_cdrm, o_cdrh, o_cldl, o_cldm, o_cldh, &
64	o_cldt, o_JrNt, o_cldljn, o_cldmjn, &
65	o_cldhjn, o_cldtjn, o_cldq, o_lwp, o_iwp, &
66	o_ue, o_ve, o_uq, o_vq, o_cape, o_pbase, &
67	o_ptop, o_fbase, o_plcl, o_plfc, &
68	o_wbeff, o_cape_max, o_upwd, o_ep,o_epmax_diag, o_Ma, &
69	o_dnwd, o_dnwd0, o_ftime_con, o_mc, &
70	o_prw, o_prlw, o_prsw, o_s_pblh, o_s_pblt, o_s_lcl, &
71	o_s_therm, o_uSTDlevs, o_vSTDlevs, &
72	o_wSTDlevs, o_zSTDlevs, o_qSTDlevs, &
73	o_tSTDlevs, epsfra, o_t_oce_sic, &
74	o_ale_bl, o_alp_bl, o_ale_wk, o_alp_wk, &
75	o_dtvdf_x , o_dtvdf_w , o_dqvdf_x , o_dqvdf_w , &
76	o_sens_x , o_sens_w , o_flat_x , o_flat_w , &
77	o_delta_tsurf, &
78	o_cdragh_x , o_cdragh_w , o_cdragm_x , o_cdragm_w , &
79	o_kh , o_kh_x , o_kh_w , &
80	o_ale, o_alp, o_cin, o_WAPE, o_wake_h, &
81	o_wake_s, o_wake_deltat, o_wake_deltaq, &
82	o_wake_omg, o_dtwak, o_dqwak, o_dqwak2d, o_Vprecip, &
83	o_ftd, o_fqd, o_wdtrainA, o_wdtrainM, &
84	o_n2, o_s2, o_proba_notrig, &
85	o_random_notrig, o_ale_bl_stat, &
86	o_ale_bl_trig, o_alp_bl_det, &
87	o_alp_bl_fluct_m, o_alp_bl_fluct_tke, &
88	o_alp_bl_conv, o_alp_bl_stat, &
89	o_slab_qflux, o_tslab, o_slab_bils, &
90	o_slab_bilg, o_slab_sic, o_slab_tice, &
91	o_weakinv, o_dthmin, o_cldtau, &
92	o_cldemi, o_pr_con_l, o_pr_con_i, &
93	o_pr_lsc_l, o_pr_lsc_i, o_re, o_fl, &
94	o_rh2m, o_rh2m_min, o_rh2m_max, &
95	o_qsat2m, o_tpot, o_tpote, o_SWnetOR, &
96	o_SWdownOR, o_LWdownOR, o_snowl, &
97	o_solldown, o_dtsvdfo, o_dtsvdft, &
98	o_dtsvdfg, o_dtsvdfi, o_z0m, o_z0h, o_od550aer, &
99	o_od865aer, o_absvisaer, o_od550lt1aer, &
100	o_sconcso4, o_sconcno3, o_sconcoa, o_sconcbc, &
101	o_sconcss, o_sconcdust, o_concso4, o_concno3, &
102	o_concoa, o_concbc, o_concss, o_concdust, &
103	o_loadso4, o_loadoa, o_loadbc, o_loadss, &
104	o_loaddust, o_tausumaero, o_tausumaero_lw, &
105	o_topswad, o_topswad0, o_solswad, o_solswad0, &
106	o_toplwad, o_toplwad0, o_sollwad, o_sollwad0, &
107	o_swtoaas_nat, o_swsrfas_nat, &
108	o_swtoacs_nat, o_swtoaas_ant, &
109	o_swsrfas_ant, o_swtoacs_ant, &
110	o_swsrfcs_ant, o_swtoacf_nat, &
111	o_swsrfcf_nat, o_swtoacf_ant, &
112	o_swsrfcs_nat, o_swsrfcf_ant, &
113	o_swtoacf_zero, o_swsrfcf_zero, &
114	o_topswai, o_solswai, o_scdnc, &
115	o_cldncl, o_reffclws, o_reffclwc, o_stratomask,&
116	o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
117	o_lcc3dstra, o_reffclwtop, o_ec550aer, &
118	o_lwcon, o_iwcon, o_temp, o_theta, &
119	o_ovapinit, o_ovap, o_oliq, o_ocond, o_geop, &
120	o_vitu, o_vitv, o_vitw, o_pres, o_paprs, &
121	o_zfull, o_zhalf, o_rneb, o_rnebjn, o_rnebcon, &
122	o_rnebls, o_rhum, o_ozone, o_ozone_light, &
123	o_dtphy, o_dqphy, o_dqphy2d, o_dqlphy, o_dqlphy2d, &
124	o_dqsphy, o_dqsphy2d, o_albe_srf, o_z0m_srf, o_z0h_srf, &
125	o_ages_srf, o_snow_srf, o_alb1, o_alb2, o_tke, &
126	o_tke_max, o_kz, o_kz_max, o_clwcon, &
127	o_dtdyn, o_dqdyn, o_dqdyn2d, o_dqldyn, o_dqldyn2d, &
128	o_dqsdyn, o_dqsdyn2d, o_dudyn, o_dvdyn, &
129	o_dtcon, o_tntc, o_ducon, o_dvcon, &
130	o_dqcon, o_dqcon2d, o_tnhusc, o_tnhusc, o_dtlsc, &
131	o_dtlschr, o_dqlsc, o_dqlsc2d, o_beta_prec, &
132	o_dtlscth, o_dtlscst, o_dqlscth, o_dqlscth2d, &
133	o_dqlscst, o_dqlscst2d, o_plulth, o_plulst, &
134	o_ptconvth, o_lmaxth, o_dtvdf, &
135	o_dtdis, o_dqvdf, o_dqvdf2d, o_dteva, o_dqeva, o_dqeva2d, &
136	o_ptconv, o_ratqs, o_dtthe, &
137	o_duthe, o_dvthe, o_ftime_th, &
138	o_f_th, o_e_th, o_w_th, o_q_th, &
139	o_a_th, o_d_th, o_f0_th, o_zmax_th, &
140	o_dqthe, o_dqthe2d, o_dtajs, o_dqajs, o_dqajs2d, o_dtswr, &
141	o_dtsw0, o_dtlwr, o_dtlw0, o_dtec, &
142	o_duvdf, o_dvvdf, o_duoro, o_dvoro, &
143	o_dtoro, o_dulif, o_dvlif, o_dtlif, &
144	o_du_gwd_hines, o_dv_gwd_hines, o_dthin, o_dqch4, o_rsu, &
145	o_du_gwd_front, o_dv_gwd_front, &
146	o_east_gwstress, o_west_gwstress, &
147	o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, &
148	o_rlucs, o_rldcs, o_tnt, o_tntr, &
149	o_tntscpbl, o_tnhus, o_tnhusscpbl, &
150	o_evu, o_h2o, o_mcd, o_dmc, o_ref_liq, &
151	o_ref_ice, o_rsut4co2, o_rlut4co2, &
152	o_rsutcs4co2, o_rlutcs4co2, o_rsu4co2, &
153	o_rlu4co2, o_rsucs4co2, o_rlucs4co2, &
154	o_rsd4co2, o_rld4co2, o_rsdcs4co2, &
155	o_rldcs4co2, o_tnondef, o_ta, o_zg, &
156	o_hus, o_hur, o_ua, o_va, o_wap, &
157	o_psbg, o_tro3, o_tro3_daylight, &
158	o_uxv, o_vxq, o_vxT, o_wxq, o_vxphi, &
159	o_wxT, o_uxu, o_vxv, o_TxT, o_trac, &
160	o_dtr_vdf, o_dtr_the, o_dtr_con, &
161	o_dtr_lessi_impa, o_dtr_lessi_nucl, &
162	o_dtr_insc, o_dtr_bcscav, o_dtr_evapls, &
163	o_dtr_ls, o_dtr_trsp, o_dtr_sscav, &
164	o_dtr_sat, o_dtr_uscav, o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
165	o_ustr_gwd_hines,o_vstr_gwd_hines,o_ustr_gwd_rando,o_vstr_gwd_rando, &
166	o_ustr_gwd_front,o_vstr_gwd_front, &
167	o_sens_prec_liq_oce, o_sens_prec_liq_sic, &
168	o_sens_prec_sol_oce, o_sens_prec_sol_sic, &
169	o_lat_prec_liq_oce, o_lat_prec_liq_sic, &
170	o_lat_prec_sol_oce, o_lat_prec_sol_sic, &
171	! Marine
172	o_map_prop_hc, o_map_prop_hist, o_map_emis_hc, o_map_iwp_hc, &
173	o_map_deltaz_hc, o_map_pcld_hc, o_map_tcld_hc, &
174	o_map_emis_hist, o_map_iwp_hist, o_map_deltaz_hist, &
175	o_map_rad_hist, &
176	o_map_emis_Cb, o_map_pcld_Cb, o_map_tcld_Cb, &
177	o_map_emis_ThCi, o_map_pcld_ThCi, o_map_tcld_ThCi, &
178	o_map_emis_Anv, o_map_pcld_Anv, o_map_tcld_Anv, &
179	o_map_ntot, o_map_hc,o_map_hist,o_map_Cb,o_map_ThCi,o_map_Anv, &
180	o_alt_tropo
181
182
183	USE phys_state_var_mod, only: pctsrf, paire_ter, rain_fall, snow_fall, &
184	qsol, z0m, z0h, fevap, agesno, &
185	nday_rain, rain_con, snow_con, &
186	topsw, toplw, toplw0, swup, swdn, &
187	topsw0, swup0, swdn0, SWup200, SWup200clr, &
188	SWdn200, SWdn200clr, LWup200, LWup200clr, &
189	LWdn200, LWdn200clr, solsw, solsw0, sollw, &
190	radsol, swradcorr, sollw0, sollwdown, sollw, gustiness, &
191	sollwdownclr, lwdn0, ftsol, ustar, u10m, &
192	v10m, pbl_tke, wake_delta_pbl_TKE, &
193	delta_tsurf, &
194	wstar, cape, ema_pcb, ema_pct, &
195	ema_cbmf, Ma, fm_therm, ale_bl, alp_bl, ale, &
196	alp, cin, wake_pe, wake_s, wake_deltat, &
197	wake_deltaq, ftd, fqd, ale_bl_trig, albsol1, &
198	rnebcon, wo, falb1, albsol2, coefh, clwcon0, &
199	ratqs, entr_therm, zqasc, detr_therm, f0, &
200	lwup, lwdn, lwup0, coefm, &
201	swupp, lwupp, swup0p, lwup0p, swdnp, lwdnp, &
202	swdn0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
203	vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &
204	vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &
205	T2sumSTD, nlevSTD, du_gwd_rando, du_gwd_front, &
206	ulevSTD, vlevSTD, wlevSTD, philevSTD, qlevSTD, tlevSTD, &
207	rhlevSTD, O3STD, O3daySTD, uvSTD, vqSTD, vTSTD, wqSTD, &
208	vphiSTD, wTSTD, u2STD, v2STD, T2STD, missing_val_nf90
209
210	USE phys_local_var_mod, only: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, &
211	t2m_min_mon, t2m_max_mon, evap, &
212	l_mixmin,l_mix, &
213	zu10m, zv10m, zq2m, zustar, zxqsurf, &
214	rain_lsc, rain_num, snow_lsc, bils, sens, fder, &
215	zxffonte, zxfqcalving, zxfqfonte, zxrunofflic, fluxu, &
216	fluxv, zxsnow, qsnow, snowhgt, to_ice, &
217	sissnow, runoff, albsol3_lic, evap_pot, &
218	t2m, fluxt, fluxlat, fsollw, fsolsw, &
219	wfbils, wfbilo, cdragm, cdragh, cldl, cldm, &
220	cldh, cldt, JrNt, cldljn, cldmjn, cldhjn, &
221	cldtjn, cldq, flwp, fiwp, ue, ve, uq, vq, &
222	plcl, plfc, wbeff, upwd, dnwd, dnwd0, prw, prlw, prsw, &
223	s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
224	vwriteSTD, wwriteSTD, phiwriteSTD, qwriteSTD, &
225	twriteSTD, ale_wake, alp_wake, &
226	dtvdf_x ,dtvdf_w ,dqvdf_x ,dqvdf_w , &
227	sens_x ,sens_w ,zxfluxlat_x,zxfluxlat_w, &
228	cdragh_x ,cdragh_w ,cdragm_x ,cdragm_w , &
229	kh ,kh_x ,kh_w , &
230	wake_h, &
231	wake_omg, d_t_wake, d_q_wake, Vprecip, &
232	wdtrainA, wdtrainM, n2, s2, proba_notrig, &
233	random_notrig, ale_bl_stat, &
234	alp_bl_det, alp_bl_fluct_m, alp_bl_conv, &
235	alp_bl_stat, alp_bl_fluct_tke, slab_wfbils, &
236	weak_inversion, dthmin, cldtau, cldemi, &
237	pmflxr, pmflxs, prfl, psfl, re, fl, rh2m, &
238	qsat2m, tpote, tpot, d_ts, od550aer, &
239	od865aer, absvisaer, od550lt1aer, sconcso4, sconcno3, &
240	sconcoa, sconcbc, sconcss, sconcdust, concso4, concno3, &
241	concoa, concbc, concss, concdust, loadso4, &
242	loadoa, loadbc, loadss, loaddust, tausum_aero, &
243	topswad_aero, topswad0_aero, solswad_aero, &
244	solswad0_aero, topsw_aero, solsw_aero, &
245	topsw0_aero, solsw0_aero, topswcf_aero, &
246	solswcf_aero, topswai_aero, solswai_aero, &
247	toplwad_aero, toplwad0_aero, sollwad_aero, &
248	sollwad0_aero, toplwai_aero, sollwai_aero, &
249	scdnc, cldncl, reffclws, reffclwc, cldnvi, stratomask,&
250	lcc, lcc3d, lcc3dcon, lcc3dstra, reffclwtop, &
251	ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
252	ql_seri, qs_seri, tr_seri, &
253	zphi, u_seri, v_seri, omega, cldfra, &
254	rneb, rnebjn, zx_rh, d_t_dyn, &
255	d_q_dyn, d_ql_dyn, d_qs_dyn, &
256	d_q_dyn2d, d_ql_dyn2d, d_qs_dyn2d, &
257	d_u_dyn, d_v_dyn, d_t_con, d_t_ajsb, d_t_ajs, &
258	d_u_ajs, d_v_ajs, &
259	d_u_con, d_v_con, d_q_con, d_q_ajs, d_t_lsc, &
260	d_t_lwr,d_t_lw0,d_t_swr,d_t_sw0, &
261	d_t_eva, d_q_lsc, beta_prec, d_t_lscth, &
262	d_t_lscst, d_q_lscth, d_q_lscst, plul_th, &
263	plul_st, d_t_vdf, d_t_diss, d_q_vdf, d_q_eva, &
264	zw2, fraca, zmax_th, d_q_ajsb, d_t_ec, d_u_vdf, &
265	d_v_vdf, d_u_oro, d_v_oro, d_t_oro, d_u_lif, &
266	d_v_lif, d_t_lif, du_gwd_hines, dv_gwd_hines, d_t_hin, &
267	dv_gwd_rando, dv_gwd_front, &
268	east_gwstress, west_gwstress, &
269	d_q_ch4, pmfd, pmfu, ref_liq, ref_ice, rhwriteSTD, &
270	ep, epmax_diag ! epmax_cape
271
272	USE phys_output_var_mod, only: vars_defined, snow_o, zfra_o, bils_diss, &
273	bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
274	itau_con, nfiles, clef_files, nid_files, &
275	zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
276	zustr_gwd_front, zvstr_gwd_front, &
277	sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
278	! Marine
279	map_prop_hc, map_prop_hist, &
280	map_emis_hc,map_iwp_hc,map_deltaz_hc,&
281	map_pcld_hc,map_tcld_hc,&
282	map_emis_hist,map_iwp_hist,map_deltaz_hist,&
283	map_rad_hist,&
284	map_ntot,map_hc,map_hist,&
285	map_Cb,map_ThCi,map_Anv,&
286	map_emis_Cb,map_pcld_Cb,map_tcld_Cb,&
287	map_emis_ThCi,map_pcld_ThCi,map_tcld_ThCi,&
288	map_emis_Anv,map_pcld_Anv,map_tcld_Anv, &
289	alt_tropo
290
291
292
293	USE ocean_slab_mod, only: tslab, slab_bils, slab_bilg, tice, seaice
294	USE pbl_surface_mod, only: snow
295	USE indice_sol_mod, only: nbsrf
296	USE infotrac_phy, only: nqtot, nqo, type_trac
297	USE geometry_mod, only: cell_area
298	USE surface_data, only: type_ocean, version_ocean, ok_veget, ok_snow
299	! USE aero_mod, only: naero_spc
300	USE aero_mod, only: naero_tot, id_STRAT_phy
301	USE ioipsl, only: histend, histsync
302	USE iophy, only: set_itau_iophy, histwrite_phy
303	USE netcdf, only: nf90_fill_real
304	USE print_control_mod, ONLY: prt_level,lunout
305
306
307	#ifdef CPP_XIOS
308	! ug Pour les sorties XIOS
309	USE xios, ONLY: xios_update_calendar
310	USE wxios, only: wxios_closedef, missing_val
311	#endif
312	USE phys_cal_mod, only : mth_len
313
314
315	IMPLICIT NONE
316
317
318	INCLUDE "clesphys.h"
319	INCLUDE "thermcell.h"
320	INCLUDE "compbl.h"
321	INCLUDE "YOMCST.h"
322
323	! Input
324	INTEGER :: itap, ivap, iliq, isol, read_climoz
325	INTEGER, DIMENSION(klon) :: lmax_th
326	LOGICAL :: aerosol_couple, ok_sync
327	LOGICAL :: ok_ade, ok_aie, new_aod
328	LOGICAL, DIMENSION(klon, klev) :: ptconv, ptconvth
329	REAL :: pdtphys
330	CHARACTER (LEN=4), DIMENSION(nlevSTD) :: clevSTD
331	REAL, DIMENSION(klon,nlevSTD) :: zx_tmp_fi3d_STD
332	REAL, DIMENSION(klon) :: pphis
333	REAL, DIMENSION(klon, klev) :: pplay, d_t
334	REAL, DIMENSION(klon, klev+1) :: paprs
335	REAL, DIMENSION(klon,klev,nqtot) :: qx, d_qx
336	REAL, DIMENSION(klon, klev) :: zmasse
337	INTEGER :: flag_aerosol_strat
338	INTEGER :: flag_aerosol
339	LOGICAL :: ok_cdnc
340	REAL, DIMENSION(3) :: freq_moyNMC
341
342	! Local
343	INTEGER :: itau_w
344	INTEGER :: i, iinit, iinitend=1, iff, iq, nsrf, k, ll, naero
345	REAL, DIMENSION (klon) :: zx_tmp_fi2d
346	REAL, DIMENSION (klon,klev) :: zx_tmp_fi3d, zpt_conv
347	REAL, DIMENSION (klon,klev+1) :: zx_tmp_fi3d1
348	CHARACTER (LEN=4) :: bb2
349	INTEGER, DIMENSION(nbp_lon*nbp_lat) :: ndex2d
350	INTEGER, DIMENSION(nbp_lonnbp_latklev) :: ndex3d
351	REAL, PARAMETER :: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
352	! REAL, PARAMETER :: missing_val=nf90_fill_real
353	#ifndef CPP_XIOS
354	REAL :: missing_val
355	#endif
356	REAL, PARAMETER :: un_jour=86400.
357
358	! On calcul le nouveau tau:
359	itau_w = itau_phy + itap
360	! On le donne à iophy pour que les histwrite y aient accès:
361	CALL set_itau_iophy(itau_w)
362
363	IF(.NOT.vars_defined) THEN
364	iinitend = 2
365	ELSE
366	iinitend = 1
367	ENDIF
368
369	! ug la boucle qui suit ne sert qu'une fois, pour l'initialisation, sinon il n'y a toujours qu'un seul passage:
370	DO iinit=1, iinitend
371	#ifdef CPP_XIOS
372	!$OMP MASTER
373	IF (vars_defined) THEN
374	if (prt_level >= 10) then
375	write(lunout,*)"phys_output_write: call xios_update_calendar, itau_w=",itau_w
376	endif
377	! CALL xios_update_calendar(itau_w)
378	CALL xios_update_calendar(itap)
379	END IF
380	!$OMP END MASTER
381	!$OMP BARRIER
382	#endif
383	! On procède à l'écriture ou à la définition des nombreuses variables:
384	!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
385	CALL histwrite_phy(o_phis, pphis)
386
387	zx_tmp_fi2d = cell_area
388	if (is_north_pole_phy) then
389	zx_tmp_fi2d(1) = cell_area(1)/nbp_lon
390	endif
391	if (is_south_pole_phy) then
392	zx_tmp_fi2d(klon) = cell_area(klon)/nbp_lon
393	endif
394	CALL histwrite_phy(o_aire, zx_tmp_fi2d)
395
396	IF (vars_defined) THEN
397	DO i=1, klon
398	zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
399	ENDDO
400	ENDIF
401
402	CALL histwrite_phy(o_contfracATM, zx_tmp_fi2d)
403	CALL histwrite_phy(o_contfracOR, pctsrf(:,is_ter))
404	CALL histwrite_phy(o_aireTER, paire_ter)
405	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
406	! Simulateur AIRS
407	IF (ok_airs) then
408	CALL histwrite_phy(o_alt_tropo,alt_tropo)
409
410	CALL histwrite_phy(o_map_prop_hc,map_prop_hc)
411	CALL histwrite_phy(o_map_prop_hist,map_prop_hist)
412
413	CALL histwrite_phy(o_map_emis_hc,map_emis_hc)
414	CALL histwrite_phy(o_map_iwp_hc,map_iwp_hc)
415	CALL histwrite_phy(o_map_deltaz_hc,map_deltaz_hc)
416	CALL histwrite_phy(o_map_pcld_hc,map_pcld_hc)
417	CALL histwrite_phy(o_map_tcld_hc,map_tcld_hc)
418
419	CALL histwrite_phy(o_map_emis_hist,map_emis_hist)
420	CALL histwrite_phy(o_map_iwp_hist,map_iwp_hist)
421	CALL histwrite_phy(o_map_deltaz_hist,map_deltaz_hist)
422
423	CALL histwrite_phy(o_map_ntot,map_ntot)
424	CALL histwrite_phy(o_map_hc,map_hc)
425	CALL histwrite_phy(o_map_hist,map_hist)
426
427	CALL histwrite_phy(o_map_Cb,map_Cb)
428	CALL histwrite_phy(o_map_ThCi,map_ThCi)
429	CALL histwrite_phy(o_map_Anv,map_Anv)
430
431	CALL histwrite_phy(o_map_emis_Cb,map_emis_Cb)
432	CALL histwrite_phy(o_map_pcld_Cb,map_pcld_Cb)
433	CALL histwrite_phy(o_map_tcld_Cb,map_tcld_Cb)
434
435	CALL histwrite_phy(o_map_emis_ThCi,map_emis_ThCi)
436	CALL histwrite_phy(o_map_pcld_ThCi,map_pcld_ThCi)
437	CALL histwrite_phy(o_map_tcld_ThCi,map_tcld_ThCi)
438
439	CALL histwrite_phy(o_map_emis_Anv,map_emis_Anv)
440	CALL histwrite_phy(o_map_pcld_Anv,map_pcld_Anv)
441	CALL histwrite_phy(o_map_tcld_Anv,map_tcld_Anv)
442	ENDIF
443
444	CALL histwrite_phy(o_flat, zxfluxlat)
445	CALL histwrite_phy(o_ptstar, ptstar)
446	CALL histwrite_phy(o_pt0, pt0)
447	CALL histwrite_phy(o_slp, slp)
448	CALL histwrite_phy(o_tsol, zxtsol)
449	CALL histwrite_phy(o_t2m, zt2m)
450	CALL histwrite_phy(o_t2m_min, zt2m)
451	CALL histwrite_phy(o_t2m_max, zt2m)
452	CALL histwrite_phy(o_t2m_max_mon, t2m_max_mon)
453	CALL histwrite_phy(o_t2m_min_mon, t2m_min_mon)
454
455	IF (vars_defined) THEN
456	DO i=1, klon
457	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
458	ENDDO
459	ENDIF
460	CALL histwrite_phy(o_wind10m, zx_tmp_fi2d)
461
462	IF (vars_defined) THEN
463	DO i=1, klon
464	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
465	ENDDO
466	ENDIF
467	CALL histwrite_phy(o_wind10max, zx_tmp_fi2d)
468
469	CALL histwrite_phy(o_gusts, gustiness)
470
471	IF (vars_defined) THEN
472	DO i = 1, klon
473	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
474	ENDDO
475	ENDIF
476	CALL histwrite_phy(o_sicf, zx_tmp_fi2d)
477	CALL histwrite_phy(o_q2m, zq2m)
478	CALL histwrite_phy(o_ustar, zustar)
479	CALL histwrite_phy(o_u10m, zu10m)
480	CALL histwrite_phy(o_v10m, zv10m)
481
482	IF (vars_defined) THEN
483	DO i = 1, klon
484	zx_tmp_fi2d(i) = paprs(i,1)
485	ENDDO
486	ENDIF
487	CALL histwrite_phy(o_psol, zx_tmp_fi2d)
488	CALL histwrite_phy(o_mass, zmasse)
489	CALL histwrite_phy(o_qsurf, zxqsurf)
490
491	IF (.NOT. ok_veget) THEN
492	CALL histwrite_phy(o_qsol, qsol)
493	ENDIF
494
495	IF (vars_defined) THEN
496	DO i = 1, klon
497	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
498	ENDDO
499	ENDIF
500
501	CALL histwrite_phy(o_precip, zx_tmp_fi2d)
502	CALL histwrite_phy(o_ndayrain, nday_rain)
503
504	! epmax_cape:
505	! CALL histwrite_phy(o_epmax_diag, epmax_diag)
506	CALL histwrite_phy(o_ep, ep)
507
508	IF (vars_defined) THEN
509	DO i = 1, klon
510	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
511	ENDDO
512	ENDIF
513	CALL histwrite_phy(o_plul, zx_tmp_fi2d)
514	CALL histwrite_phy(o_plun, rain_num)
515
516	IF (vars_defined) THEN
517	DO i = 1, klon
518	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
519	ENDDO
520	ENDIF
521	CALL histwrite_phy(o_pluc, zx_tmp_fi2d)
522	CALL histwrite_phy(o_snow, snow_fall)
523	CALL histwrite_phy(o_msnow, zxsnow)
524	CALL histwrite_phy(o_fsnow, zfra_o)
525	CALL histwrite_phy(o_evap, evap)
526	CALL histwrite_phy(o_tops, topsw*swradcorr)
527	CALL histwrite_phy(o_tops0, topsw0*swradcorr)
528	CALL histwrite_phy(o_topl, toplw)
529	CALL histwrite_phy(o_topl0, toplw0)
530
531	IF (vars_defined) THEN
532	zx_tmp_fi2d(:) = swup(:,klevp1)*swradcorr(:)
533	ENDIF
534	CALL histwrite_phy(o_SWupTOA, zx_tmp_fi2d)
535
536	IF (vars_defined) THEN
537	zx_tmp_fi2d(:) = swup0(:,klevp1)*swradcorr(:)
538	ENDIF
539	CALL histwrite_phy(o_SWupTOAclr, zx_tmp_fi2d)
540
541	IF (vars_defined) THEN
542	zx_tmp_fi2d(:) = swdn(:,klevp1)*swradcorr(:)
543	ENDIF
544	CALL histwrite_phy(o_SWdnTOA, zx_tmp_fi2d)
545
546	IF (vars_defined) THEN
547	zx_tmp_fi2d(:) = swdn0(:,klevp1)*swradcorr(:)
548	ENDIF
549	CALL histwrite_phy(o_SWdnTOAclr, zx_tmp_fi2d)
550
551	IF (vars_defined) THEN
552	zx_tmp_fi2d(:) = topsw(:)*swradcorr(:)-toplw(:)
553	ENDIF
554	CALL histwrite_phy(o_nettop, zx_tmp_fi2d)
555	CALL histwrite_phy(o_SWup200, SWup200*swradcorr)
556	CALL histwrite_phy(o_SWup200clr, SWup200clr*swradcorr)
557	CALL histwrite_phy(o_SWdn200, SWdn200*swradcorr)
558	CALL histwrite_phy(o_SWdn200clr, SWdn200clr*swradcorr)
559	CALL histwrite_phy(o_LWup200, LWup200)
560	CALL histwrite_phy(o_LWup200clr, LWup200clr)
561	CALL histwrite_phy(o_LWdn200, LWdn200)
562	CALL histwrite_phy(o_LWdn200clr, LWdn200clr)
563	CALL histwrite_phy(o_sols, solsw*swradcorr)
564	CALL histwrite_phy(o_sols0, solsw0*swradcorr)
565	CALL histwrite_phy(o_soll, sollw)
566	CALL histwrite_phy(o_soll0, sollw0)
567	CALL histwrite_phy(o_radsol, radsol)
568
569	IF (vars_defined) THEN
570	zx_tmp_fi2d(:) = swup(:,1)*swradcorr(:)
571	ENDIF
572	CALL histwrite_phy(o_SWupSFC, zx_tmp_fi2d)
573
574	IF (vars_defined) THEN
575	zx_tmp_fi2d(:) = swup0(:,1)*swradcorr(:)
576	ENDIF
577	CALL histwrite_phy(o_SWupSFCclr, zx_tmp_fi2d)
578
579	IF (vars_defined) THEN
580	zx_tmp_fi2d(:) = swdn(:,1)*swradcorr(:)
581	ENDIF
582	CALL histwrite_phy(o_SWdnSFC, zx_tmp_fi2d)
583
584	IF (vars_defined) THEN
585	zx_tmp_fi2d(:) = swdn0(:,1)*swradcorr(:)
586	ENDIF
587	CALL histwrite_phy(o_SWdnSFCclr, zx_tmp_fi2d)
588
589	IF (vars_defined) THEN
590	zx_tmp_fi2d(:)=sollwdown(:)-sollw(:)
591	ENDIF
592	CALL histwrite_phy(o_LWupSFC, zx_tmp_fi2d)
593	CALL histwrite_phy(o_LWdnSFC, sollwdown)
594
595	IF (vars_defined) THEN
596	sollwdownclr(1:klon) = -1.*lwdn0(1:klon,1)
597	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
598	ENDIF
599	CALL histwrite_phy(o_LWupSFCclr, zx_tmp_fi2d)
600	CALL histwrite_phy(o_LWdnSFCclr, sollwdownclr)
601	CALL histwrite_phy(o_bils, bils)
602	CALL histwrite_phy(o_bils_diss, bils_diss)
603	CALL histwrite_phy(o_bils_ec, bils_ec)
604	IF (iflag_ener_conserv>=1) THEN
605	CALL histwrite_phy(o_bils_ech, bils_ech)
606	ENDIF
607	CALL histwrite_phy(o_bils_tke, bils_tke)
608	CALL histwrite_phy(o_bils_kinetic, bils_kinetic)
609	CALL histwrite_phy(o_bils_latent, bils_latent)
610	CALL histwrite_phy(o_bils_enthalp, bils_enthalp)
611
612	IF (vars_defined) THEN
613	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
614	ENDIF
615	CALL histwrite_phy(o_sens, zx_tmp_fi2d)
616	CALL histwrite_phy(o_fder, fder)
617	CALL histwrite_phy(o_ffonte, zxffonte)
618	CALL histwrite_phy(o_fqcalving, zxfqcalving)
619	CALL histwrite_phy(o_fqfonte, zxfqfonte)
620	CALL histwrite_phy(o_runofflic, zxrunofflic)
621	IF (vars_defined) THEN
622	zx_tmp_fi2d=0.
623	DO nsrf=1,nbsrf
624	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf)
625	ENDDO
626	ENDIF
627	CALL histwrite_phy(o_taux, zx_tmp_fi2d)
628
629	IF (vars_defined) THEN
630	zx_tmp_fi2d=0.
631	DO nsrf=1,nbsrf
632	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf)
633	ENDDO
634	ENDIF
635	CALL histwrite_phy(o_tauy, zx_tmp_fi2d)
636
637	IF (ok_snow) THEN
638	CALL histwrite_phy(o_snowsrf, snow_o)
639	CALL histwrite_phy(o_qsnow, qsnow)
640	CALL histwrite_phy(o_snowhgt,snowhgt)
641	CALL histwrite_phy(o_toice,to_ice)
642	CALL histwrite_phy(o_sissnow,sissnow)
643	CALL histwrite_phy(o_runoff,runoff)
644	CALL histwrite_phy(o_albslw3,albsol3_lic)
645	ENDIF
646
647	DO nsrf = 1, nbsrf
648	IF (vars_defined) zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
649	CALL histwrite_phy(o_pourc_srf(nsrf), zx_tmp_fi2d)
650	IF (vars_defined) zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
651	CALL histwrite_phy(o_fract_srf(nsrf), zx_tmp_fi2d)
652	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
653	CALL histwrite_phy(o_taux_srf(nsrf), zx_tmp_fi2d)
654	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
655	CALL histwrite_phy(o_tauy_srf(nsrf), zx_tmp_fi2d)
656	IF (vars_defined) zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
657	CALL histwrite_phy(o_tsol_srf(nsrf), zx_tmp_fi2d)
658	IF (vars_defined) zx_tmp_fi2d(1 : klon) = evap_pot( 1 : klon, nsrf)
659	CALL histwrite_phy(o_evappot_srf(nsrf), zx_tmp_fi2d)
660	IF (vars_defined) zx_tmp_fi2d(1 : klon) = ustar(1 : klon, nsrf)
661	CALL histwrite_phy(o_ustar_srf(nsrf), zx_tmp_fi2d)
662	IF (vars_defined) zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf)
663	CALL histwrite_phy(o_u10m_srf(nsrf), zx_tmp_fi2d)
664	IF (vars_defined) zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf)
665	CALL histwrite_phy(o_v10m_srf(nsrf), zx_tmp_fi2d)
666	IF (vars_defined) zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf)
667	CALL histwrite_phy(o_t2m_srf(nsrf), zx_tmp_fi2d)
668	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf)
669	CALL histwrite_phy(o_evap_srf(nsrf), zx_tmp_fi2d)
670	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
671	CALL histwrite_phy(o_sens_srf(nsrf), zx_tmp_fi2d)
672	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
673	CALL histwrite_phy(o_lat_srf(nsrf), zx_tmp_fi2d)
674	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
675	CALL histwrite_phy(o_flw_srf(nsrf), zx_tmp_fi2d)
676	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
677	CALL histwrite_phy(o_fsw_srf(nsrf), zx_tmp_fi2d)
678	IF (vars_defined) zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
679	CALL histwrite_phy(o_wbils_srf(nsrf), zx_tmp_fi2d)
680	IF (vars_defined) zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
681	CALL histwrite_phy(o_wbilo_srf(nsrf), zx_tmp_fi2d)
682
683	IF (iflag_pbl > 1) THEN
684	CALL histwrite_phy(o_tke_srf(nsrf), pbl_tke(:,1:klev,nsrf))
685	CALL histwrite_phy(o_l_mix(nsrf), l_mix(:,1:klev,nsrf))
686	CALL histwrite_phy(o_l_mixmin(nsrf), l_mixmin(:,1:klev,nsrf))
687	CALL histwrite_phy(o_tke_max_srf(nsrf), pbl_tke(:,1:klev,nsrf))
688	ENDIF
689	!jyg<
690	IF (iflag_pbl > 1 .AND. iflag_wake>=1 .AND. iflag_pbl_split >=1) THEN
691	CALL histwrite_phy(o_dltpbltke_srf(nsrf), wake_delta_pbl_TKE(:,1:klev,nsrf))
692	ENDIF
693	!>jyg
694
695	ENDDO
696
697	IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_liq_o(1 : klon, 1)
698	CALL histwrite_phy(o_sens_prec_liq_oce, zx_tmp_fi2d)
699	IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_liq_o(1 : klon, 2)
700	CALL histwrite_phy(o_sens_prec_liq_sic, zx_tmp_fi2d)
701	IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_sol_o(1 : klon, 1)
702	CALL histwrite_phy(o_sens_prec_sol_oce, zx_tmp_fi2d)
703	IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_sol_o(1 : klon, 2)
704	CALL histwrite_phy(o_sens_prec_sol_sic, zx_tmp_fi2d)
705
706	IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_liq_o(1 : klon, 1)
707	CALL histwrite_phy(o_lat_prec_liq_oce, zx_tmp_fi2d)
708	IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_liq_o(1 : klon, 2)
709	CALL histwrite_phy(o_lat_prec_liq_sic, zx_tmp_fi2d)
710	IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_sol_o(1 : klon, 1)
711	CALL histwrite_phy(o_lat_prec_sol_oce, zx_tmp_fi2d)
712	IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_sol_o(1 : klon, 2)
713	CALL histwrite_phy(o_lat_prec_sol_sic, zx_tmp_fi2d)
714
715	DO nsrf=1,nbsrf+1
716	CALL histwrite_phy(o_wstar(nsrf), wstar(1 : klon, nsrf))
717	ENDDO
718
719	CALL histwrite_phy(o_cdrm, cdragm)
720	CALL histwrite_phy(o_cdrh, cdragh)
721	CALL histwrite_phy(o_cldl, cldl)
722	CALL histwrite_phy(o_cldm, cldm)
723	CALL histwrite_phy(o_cldh, cldh)
724	CALL histwrite_phy(o_cldt, cldt)
725	CALL histwrite_phy(o_JrNt, JrNt)
726	CALL histwrite_phy(o_cldljn, cldl*JrNt)
727	CALL histwrite_phy(o_cldmjn, cldm*JrNt)
728	CALL histwrite_phy(o_cldhjn, cldh*JrNt)
729	CALL histwrite_phy(o_cldtjn, cldt*JrNt)
730	CALL histwrite_phy(o_cldq, cldq)
731	IF (vars_defined) zx_tmp_fi2d(1:klon) = flwp(1:klon)
732	CALL histwrite_phy(o_lwp, zx_tmp_fi2d)
733	IF (vars_defined) zx_tmp_fi2d(1:klon) = fiwp(1:klon)
734	CALL histwrite_phy(o_iwp, zx_tmp_fi2d)
735	CALL histwrite_phy(o_ue, ue)
736	CALL histwrite_phy(o_ve, ve)
737	CALL histwrite_phy(o_uq, uq)
738	CALL histwrite_phy(o_vq, vq)
739	IF(iflag_con.GE.3) THEN ! sb
740	CALL histwrite_phy(o_cape, cape)
741	CALL histwrite_phy(o_pbase, ema_pcb)
742	CALL histwrite_phy(o_ptop, ema_pct)
743	CALL histwrite_phy(o_fbase, ema_cbmf)
744	if (iflag_con /= 30) then
745	CALL histwrite_phy(o_plcl, plcl)
746	CALL histwrite_phy(o_plfc, plfc)
747	CALL histwrite_phy(o_wbeff, wbeff)
748	end if
749
750	CALL histwrite_phy(o_cape_max, cape)
751
752	CALL histwrite_phy(o_upwd, upwd)
753	CALL histwrite_phy(o_Ma, Ma)
754	CALL histwrite_phy(o_dnwd, dnwd)
755	CALL histwrite_phy(o_dnwd0, dnwd0)
756	IF (vars_defined) zx_tmp_fi2d=float(itau_con)/float(itap)
757	CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
758	IF (vars_defined) THEN
759	IF(iflag_thermals>=1)THEN
760	zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm(:,1:klev)
761	ELSE
762	zx_tmp_fi3d=dnwd+dnwd0+upwd
763	ENDIF
764	ENDIF
765	CALL histwrite_phy(o_mc, zx_tmp_fi3d)
766	ENDIF !iflag_con .GE. 3
767	CALL histwrite_phy(o_prw, prw)
768	CALL histwrite_phy(o_prlw, prlw)
769	CALL histwrite_phy(o_prsw, prsw)
770	CALL histwrite_phy(o_s_pblh, s_pblh)
771	CALL histwrite_phy(o_s_pblt, s_pblt)
772	CALL histwrite_phy(o_s_lcl, s_lcl)
773	CALL histwrite_phy(o_s_therm, s_therm)
774	!IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
775	! IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN
776	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
777	! $o_s_capCL%name,itau_w,s_capCL)
778	! ENDIF
779	! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
780	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
781	! $o_s_oliqCL%name,itau_w,s_oliqCL)
782	! ENDIF
783	! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
784	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
785	! $o_s_cteiCL%name,itau_w,s_cteiCL)
786	! ENDIF
787	! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
788	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
789	! $o_s_trmb1%name,itau_w,s_trmb1)
790	! ENDIF
791	! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
792	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
793	! $o_s_trmb2%name,itau_w,s_trmb2)
794	! ENDIF
795	! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
796	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
797	! $o_s_trmb3%name,itau_w,s_trmb3)
798	! ENDIF
799
800	#ifdef CPP_IOIPSL
801	#ifndef CPP_XIOS
802	IF (.NOT.ok_all_xml) THEN
803	! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX:
804	! Champs interpolles sur des niveaux de pression
805	missing_val=missing_val_nf90
806	DO iff=1, nfiles
807	ll=0
808	DO k=1, nlevSTD
809	bb2=clevSTD(k)
810	IF(bb2.EQ."850".OR.bb2.EQ."700".OR. &
811	bb2.EQ."500".OR.bb2.EQ."200".OR. &
812	bb2.EQ."100".OR. &
813	bb2.EQ."50".OR.bb2.EQ."10") THEN
814
815	! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
816	ll=ll+1
817	CALL histwrite_phy(o_uSTDlevs(ll),uwriteSTD(:,k,iff), iff)
818	CALL histwrite_phy(o_vSTDlevs(ll),vwriteSTD(:,k,iff), iff)
819	CALL histwrite_phy(o_wSTDlevs(ll),wwriteSTD(:,k,iff), iff)
820	CALL histwrite_phy(o_zSTDlevs(ll),phiwriteSTD(:,k,iff), iff)
821	CALL histwrite_phy(o_qSTDlevs(ll),qwriteSTD(:,k,iff), iff)
822	CALL histwrite_phy(o_tSTDlevs(ll),twriteSTD(:,k,iff), iff)
823
824	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
825	ENDDO
826	ENDDO
827	ENDIF
828	#endif
829	#endif
830	#ifdef CPP_XIOS
831	IF(ok_all_xml) THEN
832	!XIOS CALL xios_get_field_attr("u850",default_value=missing_val)
833	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
834	ll=0
835	DO k=1, nlevSTD
836	bb2=clevSTD(k)
837	IF(bb2.EQ."850".OR.bb2.EQ."700".OR. &
838	bb2.EQ."500".OR.bb2.EQ."200".OR. &
839	bb2.EQ."100".OR. &
840	bb2.EQ."50".OR.bb2.EQ."10") THEN
841	ll=ll+1
842	CALL histwrite_phy(o_uSTDlevs(ll),ulevSTD(:,k))
843	CALL histwrite_phy(o_vSTDlevs(ll),vlevSTD(:,k))
844	CALL histwrite_phy(o_wSTDlevs(ll),wlevSTD(:,k))
845	CALL histwrite_phy(o_zSTDlevs(ll),philevSTD(:,k))
846	CALL histwrite_phy(o_qSTDlevs(ll),qlevSTD(:,k))
847	CALL histwrite_phy(o_tSTDlevs(ll),tlevSTD(:,k))
848	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
849	ENDDO
850	ENDIF
851	#endif
852	IF (vars_defined) THEN
853	DO i=1, klon
854	IF (pctsrf(i,is_oce).GT.epsfra.OR. &
855	pctsrf(i,is_sic).GT.epsfra) THEN
856	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ &
857	ftsol(i, is_sic) * pctsrf(i,is_sic))/ &
858	(pctsrf(i,is_oce)+pctsrf(i,is_sic))
859	ELSE
860	zx_tmp_fi2d(i) = 273.15
861	ENDIF
862	ENDDO
863	ENDIF
864	CALL histwrite_phy(o_t_oce_sic, zx_tmp_fi2d)
865
866	! Couplage convection-couche limite
867	IF (iflag_con.GE.3) THEN
868	IF (iflag_coupl>=1) THEN
869	CALL histwrite_phy(o_ale_bl, ale_bl)
870	CALL histwrite_phy(o_alp_bl, alp_bl)
871	ENDIF !iflag_coupl>=1
872	ENDIF !(iflag_con.GE.3)
873	! Wakes
874	IF (iflag_con.EQ.3) THEN
875	IF (iflag_wake>=1) THEN
876	CALL histwrite_phy(o_ale_wk, ale_wake)
877	CALL histwrite_phy(o_alp_wk, alp_wake)
878	IF (iflag_pbl_split>=1) THEN
879	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_x(1:klon,1:klev)/pdtphys
880	CALL histwrite_phy(o_dtvdf_x ,zx_tmp_fi3d)
881	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_w(1:klon,1:klev)/pdtphys
882	CALL histwrite_phy(o_dtvdf_w ,zx_tmp_fi3d)
883	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_x(1:klon,1:klev)/pdtphys
884	CALL histwrite_phy(o_dqvdf_x ,zx_tmp_fi3d)
885	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_w(1:klon,1:klev)/pdtphys
886	CALL histwrite_phy(o_dqvdf_w ,zx_tmp_fi3d)
887	CALL histwrite_phy(o_sens_x ,sens_x )
888	CALL histwrite_phy(o_sens_w ,sens_w )
889	CALL histwrite_phy(o_flat_x ,zxfluxlat_x)
890	CALL histwrite_phy(o_flat_w ,zxfluxlat_w)
891	CALL histwrite_phy(o_delta_tsurf,delta_tsurf)
892	CALL histwrite_phy(o_cdragh_x ,cdragh_x )
893	CALL histwrite_phy(o_cdragh_w ,cdragh_w )
894	CALL histwrite_phy(o_cdragm_x ,cdragm_x )
895	CALL histwrite_phy(o_cdragm_w ,cdragm_w )
896	CALL histwrite_phy(o_kh ,kh )
897	CALL histwrite_phy(o_kh_x ,kh_x )
898	CALL histwrite_phy(o_kh_w ,kh_w )
899	ENDIF ! (iflag_pbl_split>=1)
900	CALL histwrite_phy(o_ale, ale)
901	CALL histwrite_phy(o_alp, alp)
902	CALL histwrite_phy(o_cin, cin)
903	CALL histwrite_phy(o_WAPE, wake_pe)
904	CALL histwrite_phy(o_wake_h, wake_h)
905	CALL histwrite_phy(o_wake_s, wake_s)
906	CALL histwrite_phy(o_wake_deltat, wake_deltat)
907	CALL histwrite_phy(o_wake_deltaq, wake_deltaq)
908	CALL histwrite_phy(o_wake_omg, wake_omg)
909	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev) &
910	/pdtphys
911	CALL histwrite_phy(o_dtwak, zx_tmp_fi3d)
912	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
913	CALL histwrite_phy(o_dqwak, zx_tmp_fi3d)
914	CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
915	CALL histwrite_phy(o_dqwak2d, zx_tmp_fi2d)
916	ENDIF ! iflag_wake>=1
917	CALL histwrite_phy(o_ftd, ftd)
918	CALL histwrite_phy(o_fqd, fqd)
919	ENDIF !(iflag_con.EQ.3)
920	IF (iflag_con.EQ.3.OR.iflag_con.EQ.30) THEN
921	! sortie RomP convection descente insaturee iflag_con=30
922	! etendue a iflag_con=3 (jyg)
923	CALL histwrite_phy(o_Vprecip, Vprecip)
924	CALL histwrite_phy(o_wdtrainA, wdtrainA)
925	CALL histwrite_phy(o_wdtrainM, wdtrainM)
926	ENDIF !(iflag_con.EQ.3.or.iflag_con.EQ.30)
927	!!! nrlmd le 10/04/2012
928	IF (iflag_trig_bl>=1) THEN
929	CALL histwrite_phy(o_n2, n2)
930	CALL histwrite_phy(o_s2, s2)
931	CALL histwrite_phy(o_proba_notrig, proba_notrig)
932	CALL histwrite_phy(o_random_notrig, random_notrig)
933	CALL histwrite_phy(o_ale_bl_stat, ale_bl_stat)
934	CALL histwrite_phy(o_ale_bl_trig, ale_bl_trig)
935	ENDIF !(iflag_trig_bl>=1)
936	IF (iflag_clos_bl>=1) THEN
937	CALL histwrite_phy(o_alp_bl_det, alp_bl_det)
938	CALL histwrite_phy(o_alp_bl_fluct_m, alp_bl_fluct_m)
939	CALL histwrite_phy(o_alp_bl_fluct_tke, &
940	alp_bl_fluct_tke)
941	CALL histwrite_phy(o_alp_bl_conv, alp_bl_conv)
942	CALL histwrite_phy(o_alp_bl_stat, alp_bl_stat)
943	ENDIF !(iflag_clos_bl>=1)
944	!!! fin nrlmd le 10/04/2012
945	! Output of slab ocean variables
946	IF (type_ocean=='slab ') THEN
947	CALL histwrite_phy(o_slab_qflux, slab_wfbils)
948	CALL histwrite_phy(o_slab_bils, slab_bils)
949	IF (nslay.EQ.1) THEN
950	zx_tmp_fi2d(:)=tslab(:,1)
951	CALL histwrite_phy(o_tslab, zx_tmp_fi2d)
952	ELSE
953	CALL histwrite_phy(o_tslab, tslab)
954	END IF
955	IF (version_ocean=='sicINT') THEN
956	CALL histwrite_phy(o_slab_bilg, slab_bilg)
957	CALL histwrite_phy(o_slab_tice, tice)
958	CALL histwrite_phy(o_slab_sic, seaice)
959	END IF
960	ENDIF !type_ocean == force/slab
961	CALL histwrite_phy(o_weakinv, weak_inversion)
962	CALL histwrite_phy(o_dthmin, dthmin)
963	CALL histwrite_phy(o_cldtau, cldtau)
964	CALL histwrite_phy(o_cldemi, cldemi)
965	CALL histwrite_phy(o_pr_con_l, pmflxr(:,1:klev))
966	CALL histwrite_phy(o_pr_con_i, pmflxs(:,1:klev))
967	CALL histwrite_phy(o_pr_lsc_l, prfl(:,1:klev))
968	CALL histwrite_phy(o_pr_lsc_i, psfl(:,1:klev))
969	CALL histwrite_phy(o_re, re)
970	CALL histwrite_phy(o_fl, fl)
971	IF (vars_defined) THEN
972	DO i=1, klon
973	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
974	ENDDO
975	ENDIF
976	CALL histwrite_phy(o_rh2m, zx_tmp_fi2d)
977
978	IF (vars_defined) THEN
979	DO i=1, klon
980	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
981	ENDDO
982	ENDIF
983	CALL histwrite_phy(o_rh2m_min, zx_tmp_fi2d)
984
985	IF (vars_defined) THEN
986	DO i=1, klon
987	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
988	ENDDO
989	ENDIF
990	CALL histwrite_phy(o_rh2m_max, zx_tmp_fi2d)
991
992	CALL histwrite_phy(o_qsat2m, qsat2m)
993	CALL histwrite_phy(o_tpot, tpot)
994	CALL histwrite_phy(o_tpote, tpote)
995	IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
996	CALL histwrite_phy(o_SWnetOR, zx_tmp_fi2d)
997	IF (vars_defined) zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon))
998	CALL histwrite_phy(o_SWdownOR, zx_tmp_fi2d)
999	CALL histwrite_phy(o_LWdownOR, sollwdown)
1000	CALL histwrite_phy(o_snowl, snow_lsc)
1001	CALL histwrite_phy(o_solldown, sollwdown)
1002	CALL histwrite_phy(o_dtsvdfo, d_ts(:,is_oce))
1003	CALL histwrite_phy(o_dtsvdft, d_ts(:,is_ter))
1004	CALL histwrite_phy(o_dtsvdfg, d_ts(:,is_lic))
1005	CALL histwrite_phy(o_dtsvdfi, d_ts(:,is_sic))
1006	CALL histwrite_phy(o_z0m, z0m(:,nbsrf+1))
1007	CALL histwrite_phy(o_z0h, z0h(:,nbsrf+1))
1008	! OD550 per species
1009	!--OLIVIER
1010	!This is warranted by treating INCA aerosols as offline aerosols
1011	! IF (new_aod .and. (.not. aerosol_couple)) THEN
1012	IF (new_aod) THEN
1013	IF (flag_aerosol.GT.0) THEN
1014	CALL histwrite_phy(o_od550aer, od550aer)
1015	CALL histwrite_phy(o_od865aer, od865aer)
1016	CALL histwrite_phy(o_absvisaer, absvisaer)
1017	CALL histwrite_phy(o_od550lt1aer, od550lt1aer)
1018	CALL histwrite_phy(o_sconcso4, sconcso4)
1019	CALL histwrite_phy(o_sconcno3, sconcno3)
1020	CALL histwrite_phy(o_sconcoa, sconcoa)
1021	CALL histwrite_phy(o_sconcbc, sconcbc)
1022	CALL histwrite_phy(o_sconcss, sconcss)
1023	CALL histwrite_phy(o_sconcdust, sconcdust)
1024	CALL histwrite_phy(o_concso4, concso4)
1025	CALL histwrite_phy(o_concno3, concno3)
1026	CALL histwrite_phy(o_concoa, concoa)
1027	CALL histwrite_phy(o_concbc, concbc)
1028	CALL histwrite_phy(o_concss, concss)
1029	CALL histwrite_phy(o_concdust, concdust)
1030	CALL histwrite_phy(o_loadso4, loadso4)
1031	CALL histwrite_phy(o_loadoa, loadoa)
1032	CALL histwrite_phy(o_loadbc, loadbc)
1033	CALL histwrite_phy(o_loadss, loadss)
1034	CALL histwrite_phy(o_loaddust, loaddust)
1035	!--STRAT AER
1036	ENDIF
1037	IF (flag_aerosol.GT.0.OR.flag_aerosol_strat.GT.0) THEN
1038	! DO naero = 1, naero_spc
1039	!--correction mini bug OB
1040	DO naero = 1, naero_tot
1041	CALL histwrite_phy(o_tausumaero(naero), &
1042	tausum_aero(:,2,naero) )
1043	END DO
1044	ENDIF
1045	IF (flag_aerosol_strat.GT.0) THEN
1046	CALL histwrite_phy(o_tausumaero_lw, &
1047	tausum_aero(:,6,id_STRAT_phy) )
1048	ENDIF
1049	ENDIF
1050	IF (ok_ade) THEN
1051	CALL histwrite_phy(o_topswad, topswad_aero*swradcorr)
1052	CALL histwrite_phy(o_topswad0, topswad0_aero*swradcorr)
1053	CALL histwrite_phy(o_solswad, solswad_aero*swradcorr)
1054	CALL histwrite_phy(o_solswad0, solswad0_aero*swradcorr)
1055	CALL histwrite_phy(o_toplwad, toplwad_aero)
1056	CALL histwrite_phy(o_toplwad0, toplwad0_aero)
1057	CALL histwrite_phy(o_sollwad, sollwad_aero)
1058	CALL histwrite_phy(o_sollwad0, sollwad0_aero)
1059	!====MS forcing diagnostics
1060	if (new_aod) then
1061	zx_tmp_fi2d(:)=topsw_aero(:,1)*swradcorr(:)
1062	CALL histwrite_phy(o_swtoaas_nat,zx_tmp_fi2d)
1063	zx_tmp_fi2d(:)=solsw_aero(:,1)*swradcorr(:)
1064	CALL histwrite_phy(o_swsrfas_nat,zx_tmp_fi2d)
1065	zx_tmp_fi2d(:)=topsw0_aero(:,1)*swradcorr(:)
1066	CALL histwrite_phy(o_swtoacs_nat,zx_tmp_fi2d)
1067	zx_tmp_fi2d(:)=solsw0_aero(:,1)*swradcorr(:)
1068	CALL histwrite_phy(o_swsrfcs_nat,zx_tmp_fi2d)
1069	!ant
1070	zx_tmp_fi2d(:)=topsw_aero(:,2)*swradcorr(:)
1071	CALL histwrite_phy(o_swtoaas_ant,zx_tmp_fi2d)
1072	zx_tmp_fi2d(:)=solsw_aero(:,2)*swradcorr(:)
1073	CALL histwrite_phy(o_swsrfas_ant,zx_tmp_fi2d)
1074	zx_tmp_fi2d(:)=topsw0_aero(:,2)*swradcorr(:)
1075	CALL histwrite_phy(o_swtoacs_ant,zx_tmp_fi2d)
1076	zx_tmp_fi2d(:)=solsw0_aero(:,2)*swradcorr(:)
1077	CALL histwrite_phy(o_swsrfcs_ant,zx_tmp_fi2d)
1078	!cf
1079	if (.not. aerosol_couple) then
1080	zx_tmp_fi2d(:)=topswcf_aero(:,1)*swradcorr(:)
1081	CALL histwrite_phy(o_swtoacf_nat,zx_tmp_fi2d)
1082	zx_tmp_fi2d(:)=solswcf_aero(:,1)*swradcorr(:)
1083	CALL histwrite_phy(o_swsrfcf_nat,zx_tmp_fi2d)
1084	zx_tmp_fi2d(:)=topswcf_aero(:,2)*swradcorr(:)
1085	CALL histwrite_phy(o_swtoacf_ant,zx_tmp_fi2d)
1086	zx_tmp_fi2d(:)=solswcf_aero(:,2)*swradcorr(:)
1087	CALL histwrite_phy(o_swsrfcf_ant,zx_tmp_fi2d)
1088	zx_tmp_fi2d(:)=topswcf_aero(:,3)*swradcorr(:)
1089	CALL histwrite_phy(o_swtoacf_zero,zx_tmp_fi2d)
1090	zx_tmp_fi2d(:)=solswcf_aero(:,3)*swradcorr(:)
1091	CALL histwrite_phy(o_swsrfcf_zero,zx_tmp_fi2d)
1092	endif
1093	endif ! new_aod
1094	!====MS forcing diagnostics
1095	ENDIF
1096	IF (ok_aie) THEN
1097	CALL histwrite_phy(o_topswai, topswai_aero*swradcorr)
1098	CALL histwrite_phy(o_solswai, solswai_aero*swradcorr)
1099	ENDIF
1100	IF (flag_aerosol.GT.0.AND.ok_cdnc) THEN
1101	CALL histwrite_phy(o_scdnc, scdnc)
1102	CALL histwrite_phy(o_cldncl, cldncl)
1103	CALL histwrite_phy(o_reffclws, reffclws)
1104	CALL histwrite_phy(o_reffclwc, reffclwc)
1105	CALL histwrite_phy(o_cldnvi, cldnvi)
1106	CALL histwrite_phy(o_lcc, lcc)
1107	CALL histwrite_phy(o_lcc3d, lcc3d)
1108	CALL histwrite_phy(o_lcc3dcon, lcc3dcon)
1109	CALL histwrite_phy(o_lcc3dstra, lcc3dstra)
1110	CALL histwrite_phy(o_reffclwtop, reffclwtop)
1111	ENDIF
1112	! Champs 3D:
1113	IF (ok_ade .OR. ok_aie) then
1114	CALL histwrite_phy(o_ec550aer, ec550aer)
1115	ENDIF
1116	CALL histwrite_phy(o_lwcon, flwc)
1117	CALL histwrite_phy(o_iwcon, fiwc)
1118	CALL histwrite_phy(o_temp, t_seri)
1119	CALL histwrite_phy(o_theta, theta)
1120	CALL histwrite_phy(o_ovapinit, qx(:,:,ivap))
1121	CALL histwrite_phy(o_ovap, q_seri)
1122	CALL histwrite_phy(o_oliq, ql_seri)
1123	CALL histwrite_phy(o_ocond, ql_seri+qs_seri)
1124	CALL histwrite_phy(o_geop, zphi)
1125	CALL histwrite_phy(o_vitu, u_seri)
1126	CALL histwrite_phy(o_vitv, v_seri)
1127	CALL histwrite_phy(o_vitw, omega)
1128	CALL histwrite_phy(o_pres, pplay)
1129	CALL histwrite_phy(o_paprs, paprs(:,1:klev))
1130	CALL histwrite_phy(o_zfull,zphi/RG)
1131
1132	IF (flag_aerosol_strat.EQ.2) THEN
1133	CALL histwrite_phy(o_stratomask, stratomask)
1134	ENDIF
1135
1136	IF (vars_defined) THEN
1137	zx_tmp_fi3d(:,1)= pphis(:)/RG
1138	DO k = 2, klev
1139	DO i = 1, klon
1140	zx_tmp_fi3d(i,k) = zphi(i,k)/RG + &
1141	(zphi(i,k)-zphi(i,k-1))/RG * &
1142	(paprs(i,k)-pplay(i,k))/(pplay(i,k)-pplay(i,k-1))
1143	ENDDO
1144	ENDDO
1145	ENDIF
1146	CALL histwrite_phy(o_zhalf, zx_tmp_fi3d)
1147	CALL histwrite_phy(o_rneb, cldfra)
1148	CALL histwrite_phy(o_rnebcon, rnebcon)
1149	CALL histwrite_phy(o_rnebls, rneb)
1150	IF (vars_defined) THEN
1151	DO k=1, klev
1152	DO i=1, klon
1153	zx_tmp_fi3d(i,k)=cldfra(i,k)*JrNt(i)
1154	ENDDO
1155	ENDDO
1156	ENDIF
1157	CALL histwrite_phy(o_rnebjn, zx_tmp_fi3d)
1158	CALL histwrite_phy(o_rhum, zx_rh)
1159	CALL histwrite_phy(o_ozone, &
1160	wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
1161
1162	IF (read_climoz == 2) THEN
1163	CALL histwrite_phy(o_ozone_light, &
1164	wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
1165	ENDIF
1166
1167	CALL histwrite_phy(o_dtphy, d_t)
1168
1169	CALL histwrite_phy(o_dqphy, d_qx(:,:,ivap))
1170	CALL water_int(klon,klev,d_qx(:,:,ivap),zmasse,zx_tmp_fi2d)
1171	CALL histwrite_phy(o_dqphy2d, zx_tmp_fi2d)
1172
1173	CALL histwrite_phy(o_dqlphy, d_qx(:,:,iliq))
1174	CALL water_int(klon,klev,d_qx(:,:,iliq),zmasse,zx_tmp_fi2d)
1175	CALL histwrite_phy(o_dqlphy2d, zx_tmp_fi2d)
1176
1177	IF (nqo.EQ.3) THEN
1178	CALL histwrite_phy(o_dqsphy, d_qx(:,:,isol))
1179	CALL water_int(klon,klev,d_qx(:,:,isol),zmasse,zx_tmp_fi2d)
1180	CALL histwrite_phy(o_dqsphy2d, zx_tmp_fi2d)
1181	ELSE
1182	zx_tmp_fi3d=0.0
1183	CALL histwrite_phy(o_dqsphy, zx_tmp_fi3d)
1184	zx_tmp_fi2d=0.0
1185	CALL histwrite_phy(o_dqsphy2d, zx_tmp_fi2d)
1186	ENDIF
1187
1188	DO nsrf=1, nbsrf
1189	IF (vars_defined) zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
1190	CALL histwrite_phy(o_albe_srf(nsrf), zx_tmp_fi2d)
1191	IF (vars_defined) zx_tmp_fi2d(1 : klon) = z0m( 1 : klon, nsrf)
1192	CALL histwrite_phy(o_z0m_srf(nsrf), zx_tmp_fi2d)
1193	IF (vars_defined) zx_tmp_fi2d(1 : klon) = z0h( 1 : klon, nsrf)
1194	CALL histwrite_phy(o_z0h_srf(nsrf), zx_tmp_fi2d)
1195	IF (vars_defined) zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1196	CALL histwrite_phy(o_ages_srf(nsrf), zx_tmp_fi2d)
1197	IF (vars_defined) zx_tmp_fi2d(1 : klon) = snow( 1 : klon, nsrf)
1198	CALL histwrite_phy(o_snow_srf(nsrf), zx_tmp_fi2d)
1199	ENDDO !nsrf=1, nbsrf
1200	CALL histwrite_phy(o_alb1, albsol1)
1201	CALL histwrite_phy(o_alb2, albsol2)
1202	!FH Sorties pour la couche limite
1203	if (iflag_pbl>1) then
1204	zx_tmp_fi3d=0.
1205	IF (vars_defined) THEN
1206	do nsrf=1,nbsrf
1207	do k=1,klev
1208	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) &
1209	+pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
1210	enddo
1211	enddo
1212	ENDIF
1213	CALL histwrite_phy(o_tke, zx_tmp_fi3d)
1214
1215	CALL histwrite_phy(o_tke_max, zx_tmp_fi3d)
1216	ENDIF
1217
1218	CALL histwrite_phy(o_kz, coefh(:,:,is_ave))
1219
1220	CALL histwrite_phy(o_kz_max, coefh(:,:,is_ave))
1221
1222	CALL histwrite_phy(o_clwcon, clwcon0)
1223	CALL histwrite_phy(o_dtdyn, d_t_dyn)
1224
1225	CALL histwrite_phy(o_dqdyn, d_q_dyn)
1226
1227	CALL histwrite_phy(o_dqdyn2d,d_q_dyn2d)
1228
1229	CALL histwrite_phy(o_dqldyn, d_ql_dyn)
1230
1231	CALL histwrite_phy(o_dqldyn2d, d_ql_dyn2d)
1232
1233	CALL histwrite_phy(o_dqsdyn, d_qs_dyn)
1234
1235	CALL histwrite_phy(o_dqsdyn2d, d_qs_dyn2d)
1236
1237	CALL histwrite_phy(o_dudyn, d_u_dyn)
1238	CALL histwrite_phy(o_dvdyn, d_v_dyn)
1239
1240	IF (vars_defined) THEN
1241	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
1242	ENDIF
1243	CALL histwrite_phy(o_dtcon, zx_tmp_fi3d)
1244	if(iflag_thermals.eq.0)then
1245	IF (vars_defined) THEN
1246	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
1247	d_t_ajsb(1:klon,1:klev)/pdtphys
1248	ENDIF
1249	CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
1250	else if(iflag_thermals.ge.1.and.iflag_wake.EQ.1)then
1251	IF (vars_defined) THEN
1252	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
1253	d_t_ajs(1:klon,1:klev)/pdtphys + &
1254	d_t_wake(1:klon,1:klev)/pdtphys
1255	ENDIF
1256	CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
1257	endif
1258	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1259	CALL histwrite_phy(o_ducon, zx_tmp_fi3d)
1260	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
1261	CALL histwrite_phy(o_dvcon, zx_tmp_fi3d)
1262	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1263	CALL histwrite_phy(o_dqcon, zx_tmp_fi3d)
1264	CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1265	CALL histwrite_phy(o_dqcon2d, zx_tmp_fi2d)
1266
1267	IF(iflag_thermals.EQ.0) THEN
1268	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1269	CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d)
1270	ELSE IF(iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
1271	IF (vars_defined) THEN
1272	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys + &
1273	d_q_ajs(1:klon,1:klev)/pdtphys + &
1274	d_q_wake(1:klon,1:klev)/pdtphys
1275	ENDIF
1276	CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d)
1277	ENDIF
1278
1279	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
1280	CALL histwrite_phy(o_dtlsc, zx_tmp_fi3d)
1281	IF (vars_defined) zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+ &
1282	d_t_eva(1:klon,1:klev))/pdtphys
1283	CALL histwrite_phy(o_dtlschr, zx_tmp_fi3d)
1284	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1285	CALL histwrite_phy(o_dqlsc, zx_tmp_fi3d)
1286	CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1287	CALL histwrite_phy(o_dqlsc2d, zx_tmp_fi2d)
1288	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=beta_prec(1:klon,1:klev)
1289	CALL histwrite_phy(o_beta_prec, zx_tmp_fi3d)
1290	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1291	! Sorties specifiques a la separation thermiques/non thermiques
1292	if (iflag_thermals>=1) then
1293	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys
1294	CALL histwrite_phy(o_dtlscth, zx_tmp_fi3d)
1295	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys
1296	CALL histwrite_phy(o_dtlscst, zx_tmp_fi3d)
1297	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys
1298	CALL histwrite_phy(o_dqlscth, zx_tmp_fi3d)
1299	CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1300	CALL histwrite_phy(o_dqlscth2d, zx_tmp_fi2d)
1301	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys
1302	CALL histwrite_phy(o_dqlscst, zx_tmp_fi3d)
1303	CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1304	CALL histwrite_phy(o_dqlscst2d, zx_tmp_fi2d)
1305	CALL histwrite_phy(o_plulth, plul_th)
1306	CALL histwrite_phy(o_plulst, plul_st)
1307	IF (vars_defined) THEN
1308	do k=1,klev
1309	do i=1,klon
1310	if (ptconvth(i,k)) then
1311	zx_tmp_fi3d(i,k)=1.
1312	else
1313	zx_tmp_fi3d(i,k)=0.
1314	endif
1315	enddo
1316	enddo
1317	ENDIF
1318	CALL histwrite_phy(o_ptconvth, zx_tmp_fi3d)
1319	IF (vars_defined) THEN
1320	do i=1,klon
1321	zx_tmp_fi2d(1:klon)=lmax_th(:)
1322	enddo
1323	ENDIF
1324	CALL histwrite_phy(o_lmaxth, zx_tmp_fi2d)
1325	endif ! iflag_thermals>=1
1326	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1327	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1328	CALL histwrite_phy(o_dtvdf, zx_tmp_fi3d)
1329	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_diss(1:klon,1:klev)/pdtphys
1330	CALL histwrite_phy(o_dtdis, zx_tmp_fi3d)
1331	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1332	CALL histwrite_phy(o_dqvdf, zx_tmp_fi3d)
1333	CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1334	CALL histwrite_phy(o_dqvdf2d, zx_tmp_fi2d)
1335	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1336	CALL histwrite_phy(o_dteva, zx_tmp_fi3d)
1337	IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1338	CALL histwrite_phy(o_dqeva, zx_tmp_fi3d)
1339	CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1340	CALL histwrite_phy(o_dqeva2d, zx_tmp_fi2d)
1341	zpt_conv = 0.
1342	WHERE (ptconv) zpt_conv = 1.
1343	CALL histwrite_phy(o_ptconv, zpt_conv)
1344	CALL histwrite_phy(o_ratqs, ratqs)
1345	IF (vars_defined) THEN
1346	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys - &
1347	d_t_ajsb(1:klon,1:klev)/pdtphys
1348	ENDIF
1349	CALL histwrite_phy(o_dtthe, zx_tmp_fi3d)
1350	IF (vars_defined) THEN
1351	zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys
1352	ENDIF
1353	CALL histwrite_phy(o_duthe, zx_tmp_fi3d)
1354	IF (vars_defined) THEN
1355	zx_tmp_fi3d(1:klon,1:klev)=d_v_ajs(1:klon,1:klev)/pdtphys
1356	ENDIF
1357	CALL histwrite_phy(o_dvthe, zx_tmp_fi3d)
1358
1359	IF (iflag_thermals>=1) THEN
1360	! Pour l instant 0 a y reflichir pour les thermiques
1361	zx_tmp_fi2d=0.
1362	CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
1363	CALL histwrite_phy(o_f_th, fm_therm)
1364	CALL histwrite_phy(o_e_th, entr_therm)
1365	CALL histwrite_phy(o_w_th, zw2)
1366	CALL histwrite_phy(o_q_th, zqasc)
1367	CALL histwrite_phy(o_a_th, fraca)
1368	CALL histwrite_phy(o_d_th, detr_therm)
1369	CALL histwrite_phy(o_f0_th, f0)
1370	CALL histwrite_phy(o_zmax_th, zmax_th)
1371	IF (vars_defined) THEN
1372	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys - &
1373	d_q_ajsb(1:klon,1:klev)/pdtphys
1374	ENDIF
1375	CALL histwrite_phy(o_dqthe, zx_tmp_fi3d)
1376	CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1377	CALL histwrite_phy(o_dqthe2d, zx_tmp_fi2d)
1378	ENDIF !iflag_thermals
1379	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
1380	CALL histwrite_phy(o_dtajs, zx_tmp_fi3d)
1381	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
1382	CALL histwrite_phy(o_dqajs, zx_tmp_fi3d)
1383	CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1384	CALL histwrite_phy(o_dqajs2d, zx_tmp_fi2d)
1385	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys
1386	CALL histwrite_phy(o_dtswr, zx_tmp_fi3d)
1387	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_sw0(1:klon,1:klev)/pdtphys
1388	CALL histwrite_phy(o_dtsw0, zx_tmp_fi3d)
1389	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lwr(1:klon,1:klev)/pdtphys
1390	CALL histwrite_phy(o_dtlwr, zx_tmp_fi3d)
1391	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lw0(1:klon,1:klev)/pdtphys
1392	CALL histwrite_phy(o_dtlw0, zx_tmp_fi3d)
1393	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)/pdtphys
1394	CALL histwrite_phy(o_dtec, zx_tmp_fi3d)
1395	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1396	CALL histwrite_phy(o_duvdf, zx_tmp_fi3d)
1397	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1398	CALL histwrite_phy(o_dvvdf, zx_tmp_fi3d)
1399	IF (ok_orodr) THEN
1400	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1401	CALL histwrite_phy(o_duoro, zx_tmp_fi3d)
1402	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1403	CALL histwrite_phy(o_dvoro, zx_tmp_fi3d)
1404	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys
1405	CALL histwrite_phy(o_dtoro, zx_tmp_fi3d)
1406	ENDIF
1407	IF (ok_orolf) THEN
1408	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
1409	CALL histwrite_phy(o_dulif, zx_tmp_fi3d)
1410
1411	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
1412	CALL histwrite_phy(o_dvlif, zx_tmp_fi3d)
1413
1414	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys
1415	CALL histwrite_phy(o_dtlif, zx_tmp_fi3d)
1416	ENDIF
1417
1418	IF (ok_hines) THEN
1419	CALL histwrite_phy(o_du_gwd_hines, du_gwd_hines/pdtphys)
1420	CALL histwrite_phy(o_dv_gwd_hines, dv_gwd_hines/pdtphys)
1421	IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys
1422	CALL histwrite_phy(o_dthin, zx_tmp_fi3d)
1423	CALL histwrite_phy(o_ustr_gwd_hines, zustr_gwd_hines)
1424	CALL histwrite_phy(o_vstr_gwd_hines, zvstr_gwd_hines)
1425	end IF
1426
1427	if (.not. ok_hines .and. ok_gwd_rando) then
1428	CALL histwrite_phy(o_du_gwd_front, du_gwd_front / pdtphys)
1429	CALL histwrite_phy(o_dv_gwd_front, dv_gwd_front / pdtphys)
1430	CALL histwrite_phy(o_ustr_gwd_front, zustr_gwd_front)
1431	CALL histwrite_phy(o_vstr_gwd_front, zvstr_gwd_front)
1432	ENDIF
1433
1434	IF (ok_gwd_rando) then
1435	CALL histwrite_phy(o_du_gwd_rando, du_gwd_rando / pdtphys)
1436	CALL histwrite_phy(o_dv_gwd_rando, dv_gwd_rando / pdtphys)
1437	CALL histwrite_phy(o_ustr_gwd_rando, zustr_gwd_rando)
1438	CALL histwrite_phy(o_vstr_gwd_rando, zvstr_gwd_rando)
1439	CALL histwrite_phy(o_east_gwstress, east_gwstress )
1440	CALL histwrite_phy(o_west_gwstress, west_gwstress )
1441	end IF
1442
1443	IF (ok_qch4) then
1444	CALL histwrite_phy(o_dqch4, d_q_ch4 / pdtphys)
1445	ENDIF
1446
1447	DO k=1, klevp1
1448	zx_tmp_fi3d1(:,k)=swup(:,k)*swradcorr(:)
1449	ENDDO
1450	CALL histwrite_phy(o_rsu, zx_tmp_fi3d1)
1451	DO k=1, klevp1
1452	zx_tmp_fi3d1(:,k)=swdn(:,k)*swradcorr(:)
1453	ENDDO
1454	CALL histwrite_phy(o_rsd, zx_tmp_fi3d1)
1455	DO k=1, klevp1
1456	zx_tmp_fi3d1(:,k)=swup0(:,k)*swradcorr(:)
1457	ENDDO
1458	CALL histwrite_phy(o_rsucs, zx_tmp_fi3d1)
1459	DO k=1, klevp1
1460	zx_tmp_fi3d1(:,k)=swdn0(:,k)*swradcorr(:)
1461	ENDDO
1462	CALL histwrite_phy(o_rsdcs, zx_tmp_fi3d1)
1463
1464	CALL histwrite_phy(o_rlu, lwup)
1465	CALL histwrite_phy(o_rld, lwdn)
1466	CALL histwrite_phy(o_rlucs, lwup0)
1467	CALL histwrite_phy(o_rldcs, lwdn0)
1468
1469	IF(vars_defined) THEN
1470	zx_tmp_fi3d(1:klon,1:klev)=d_t(1:klon,1:klev)+ &
1471	d_t_dyn(1:klon,1:klev)
1472	ENDIF
1473	CALL histwrite_phy(o_tnt, zx_tmp_fi3d)
1474
1475	IF(vars_defined) THEN
1476	zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys + &
1477	d_t_lwr(1:klon,1:klev)/pdtphys
1478	ENDIF
1479	CALL histwrite_phy(o_tntr, zx_tmp_fi3d)
1480	IF(vars_defined) THEN
1481	zx_tmp_fi3d(1:klon,1:klev)= (d_t_lsc(1:klon,1:klev)+ &
1482	d_t_eva(1:klon,1:klev)+ &
1483	d_t_vdf(1:klon,1:klev))/pdtphys
1484	ENDIF
1485	CALL histwrite_phy(o_tntscpbl, zx_tmp_fi3d)
1486	IF(vars_defined) THEN
1487	zx_tmp_fi3d(1:klon,1:klev)=d_qx(1:klon,1:klev,ivap)+ &
1488	d_q_dyn(1:klon,1:klev)
1489	ENDIF
1490	CALL histwrite_phy(o_tnhus, zx_tmp_fi3d)
1491	IF(vars_defined) THEN
1492	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys+ &
1493	d_q_eva(1:klon,1:klev)/pdtphys
1494	ENDIF
1495	CALL histwrite_phy(o_tnhusscpbl, zx_tmp_fi3d)
1496	CALL histwrite_phy(o_evu, coefm(:,:,is_ave))
1497	IF(vars_defined) THEN
1498	zx_tmp_fi3d(1:klon,1:klev)=q_seri(1:klon,1:klev)+ &
1499	ql_seri(1:klon,1:klev)
1500	ENDIF
1501	CALL histwrite_phy(o_h2o, zx_tmp_fi3d)
1502	if (iflag_con >= 3) then
1503	IF(vars_defined) THEN
1504	zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+ &
1505	dnwd0(1:klon,1:klev))
1506	ENDIF
1507	CALL histwrite_phy(o_mcd, zx_tmp_fi3d)
1508	IF(vars_defined) THEN
1509	zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) + &
1510	dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev)
1511	ENDIF
1512	CALL histwrite_phy(o_dmc, zx_tmp_fi3d)
1513	else if (iflag_con == 2) then
1514	CALL histwrite_phy(o_mcd, pmfd)
1515	CALL histwrite_phy(o_dmc, pmfu + pmfd)
1516	end if
1517	CALL histwrite_phy(o_ref_liq, ref_liq)
1518	CALL histwrite_phy(o_ref_ice, ref_ice)
1519	if (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR. &
1520	RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR. &
1521	RCFC12_per.NE.RCFC12_act) THEN
1522	IF(vars_defined) zx_tmp_fi2d(:) = swupp(:,klevp1)*swradcorr(:)
1523	CALL histwrite_phy(o_rsut4co2, zx_tmp_fi2d)
1524	IF(vars_defined) zx_tmp_fi2d(:) = lwupp(:,klevp1)
1525	CALL histwrite_phy(o_rlut4co2, zx_tmp_fi2d)
1526	IF(vars_defined) zx_tmp_fi2d(:) = swup0p(:,klevp1)*swradcorr(:)
1527	CALL histwrite_phy(o_rsutcs4co2, zx_tmp_fi2d)
1528	IF(vars_defined) zx_tmp_fi2d(:) = lwup0p(:,klevp1)
1529	CALL histwrite_phy(o_rlutcs4co2, zx_tmp_fi2d)
1530	DO k=1, klevp1
1531	zx_tmp_fi3d1(:,k)=swupp(:,k)*swradcorr(:)
1532	ENDDO
1533	CALL histwrite_phy(o_rsu4co2, zx_tmp_fi3d1)
1534	DO k=1, klevp1
1535	zx_tmp_fi3d1(:,k)=swup0p(:,k)*swradcorr(:)
1536	ENDDO
1537	CALL histwrite_phy(o_rsucs4co2, zx_tmp_fi3d1)
1538	DO k=1, klevp1
1539	zx_tmp_fi3d1(:,k)=swdnp(:,k)*swradcorr(:)
1540	ENDDO
1541	CALL histwrite_phy(o_rsd4co2, zx_tmp_fi3d1)
1542	DO k=1, klevp1
1543	zx_tmp_fi3d1(:,k)=swdn0p(:,k)*swradcorr(:)
1544	ENDDO
1545	CALL histwrite_phy(o_rsdcs4co2, zx_tmp_fi3d1)
1546	CALL histwrite_phy(o_rlu4co2, lwupp)
1547	CALL histwrite_phy(o_rlucs4co2, lwup0p)
1548	CALL histwrite_phy(o_rld4co2, lwdnp)
1549	CALL histwrite_phy(o_rldcs4co2, lwdn0p)
1550	ENDIF
1551	!!!!!!!!!!!! Sorties niveaux de pression NMC !!!!!!!!!!!!!!!!!!!!
1552	#ifdef CPP_IOIPSL
1553	#ifndef CPP_XIOS
1554	IF (.NOT.ok_all_xml) THEN
1555	! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX:
1556	! Champs interpolles sur des niveaux de pression
1557	missing_val=missing_val_nf90
1558	DO iff=7, nfiles
1559
1560	CALL histwrite_phy(o_tnondef,tnondef(:,:,iff-6),iff)
1561	CALL histwrite_phy(o_ta,twriteSTD(:,:,iff-6),iff)
1562	CALL histwrite_phy(o_zg,phiwriteSTD(:,:,iff-6),iff)
1563	CALL histwrite_phy(o_hus,qwriteSTD(:,:,iff-6),iff)
1564	CALL histwrite_phy(o_hur,rhwriteSTD(:,:,iff-6),iff)
1565	CALL histwrite_phy(o_ua,uwriteSTD(:,:,iff-6),iff)
1566	CALL histwrite_phy(o_va,vwriteSTD(:,:,iff-6),iff)
1567	CALL histwrite_phy(o_wap,wwriteSTD(:,:,iff-6),iff)
1568	IF(vars_defined) THEN
1569	DO k=1, nlevSTD
1570	DO i=1, klon
1571	IF(tnondef(i,k,iff-6).NE.missing_val) THEN
1572	IF(freq_outNMC(iff-6).LT.0) THEN
1573	freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6)
1574	ELSE
1575	freq_moyNMC(iff-6)=freq_outNMC(iff-6)/freq_calNMC(iff-6)
1576	ENDIF
1577	zx_tmp_fi3d_STD(i,k) = (100.*tnondef(i,k,iff-6))/freq_moyNMC(iff-6)
1578	ELSE
1579	zx_tmp_fi3d_STD(i,k) = missing_val
1580	ENDIF
1581	ENDDO
1582	ENDDO
1583	ENDIF
1584	CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD,iff)
1585	IF(vars_defined) THEN
1586	DO k=1, nlevSTD
1587	DO i=1, klon
1588	IF(O3sumSTD(i,k,iff-6).NE.missing_val) THEN
1589	zx_tmp_fi3d_STD(i,k) = O3sumSTD(i,k,iff-6) * 1.e+9
1590	ELSE
1591	zx_tmp_fi3d_STD(i,k) = missing_val
1592	ENDIF
1593	ENDDO
1594	ENDDO !k=1, nlevSTD
1595	ENDIF
1596	CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD,iff)
1597	if (read_climoz == 2) THEN
1598	IF(vars_defined) THEN
1599	DO k=1, nlevSTD
1600	DO i=1, klon
1601	IF(O3daysumSTD(i,k,iff-6).NE.missing_val) THEN
1602	zx_tmp_fi3d_STD(i,k) = O3daysumSTD(i,k,iff-6) * 1.e+9
1603	ELSE
1604	zx_tmp_fi3d_STD(i,k) = missing_val
1605	ENDIF
1606	ENDDO
1607	ENDDO !k=1, nlevSTD
1608	ENDIF
1609	CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d_STD,iff)
1610	endif
1611	CALL histwrite_phy(o_uxv,uvsumSTD(:,:,iff-6),iff)
1612	CALL histwrite_phy(o_vxq,vqsumSTD(:,:,iff-6),iff)
1613	CALL histwrite_phy(o_vxT,vTsumSTD(:,:,iff-6),iff)
1614	CALL histwrite_phy(o_wxq,wqsumSTD(:,:,iff-6),iff)
1615	CALL histwrite_phy(o_vxphi,vphisumSTD(:,:,iff-6),iff)
1616	CALL histwrite_phy(o_wxT,wTsumSTD(:,:,iff-6),iff)
1617	CALL histwrite_phy(o_uxu,u2sumSTD(:,:,iff-6),iff)
1618	CALL histwrite_phy(o_vxv,v2sumSTD(:,:,iff-6),iff)
1619	CALL histwrite_phy(o_TxT,T2sumSTD(:,:,iff-6),iff)
1620	ENDDO !nfiles
1621	ENDIF
1622	#endif
1623	#endif
1624	#ifdef CPP_XIOS
1625	IF(ok_all_xml) THEN
1626	! DO iff=7, nfiles
1627
1628	! CALL histwrite_phy(o_tnondef,tnondef(:,:,3))
1629	CALL histwrite_phy(o_ta,tlevSTD(:,:))
1630	CALL histwrite_phy(o_zg,philevSTD(:,:))
1631	CALL histwrite_phy(o_hus,qlevSTD(:,:))
1632	CALL histwrite_phy(o_hur,rhlevSTD(:,:))
1633	CALL histwrite_phy(o_ua,ulevSTD(:,:))
1634	CALL histwrite_phy(o_va,vlevSTD(:,:))
1635	CALL histwrite_phy(o_wap,wlevSTD(:,:))
1636	! IF(vars_defined) THEN
1637	! DO k=1, nlevSTD
1638	! DO i=1, klon
1639	! IF(tnondef(i,k,3).NE.missing_val) THEN
1640	! IF(freq_outNMC(iff-6).LT.0) THEN
1641	! freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6)
1642	! ELSE
1643	! freq_moyNMC(iff-6)=freq_outNMC(iff-6)/freq_calNMC(iff-6)
1644	! ENDIF
1645	! zx_tmp_fi3d_STD(i,k) = (100.*tnondef(i,k,3))/freq_moyNMC(iff-6)
1646	! ELSE
1647	! zx_tmp_fi3d_STD(i,k) = missing_val
1648	! ENDIF
1649	! ENDDO
1650	! ENDDO
1651	! ENDIF
1652	! CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD)
1653	IF(vars_defined) THEN
1654	DO k=1, nlevSTD
1655	DO i=1, klon
1656	IF(O3STD(i,k).NE.missing_val) THEN
1657	zx_tmp_fi3d_STD(i,k) = O3STD(i,k) * 1.e+9
1658	ELSE
1659	zx_tmp_fi3d_STD(i,k) = missing_val
1660	ENDIF
1661	ENDDO
1662	ENDDO !k=1, nlevSTD
1663	ENDIF
1664	CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD)
1665	if (read_climoz == 2) THEN
1666	IF(vars_defined) THEN
1667	DO k=1, nlevSTD
1668	DO i=1, klon
1669	IF(O3daySTD(i,k).NE.missing_val) THEN
1670	zx_tmp_fi3d_STD(i,k) = O3daySTD(i,k) * 1.e+9
1671	ELSE
1672	zx_tmp_fi3d_STD(i,k) = missing_val
1673	ENDIF
1674	ENDDO
1675	ENDDO !k=1, nlevSTD
1676	ENDIF
1677	CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d_STD)
1678	endif
1679	CALL histwrite_phy(o_uxv,uvSTD(:,:))
1680	CALL histwrite_phy(o_vxq,vqSTD(:,:))
1681	CALL histwrite_phy(o_vxT,vTSTD(:,:))
1682	CALL histwrite_phy(o_wxq,wqSTD(:,:))
1683	CALL histwrite_phy(o_vxphi,vphiSTD(:,:))
1684	CALL histwrite_phy(o_wxT,wTSTD(:,:))
1685	CALL histwrite_phy(o_uxu,u2STD(:,:))
1686	CALL histwrite_phy(o_vxv,v2STD(:,:))
1687	CALL histwrite_phy(o_TxT,T2STD(:,:))
1688	! ENDDO !nfiles
1689	ENDIF
1690	#endif
1691	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1692	IF (nqtot.GE.nqo+1) THEN
1693	DO iq=nqo+1,nqtot
1694	IF (type_trac == 'lmdz' .OR. type_trac == 'repr') THEN
1695
1696	!jyg<
1697	!! CALL histwrite_phy(o_trac(iq-nqo), qx(:,:,iq))
1698	CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
1699	!>jyg
1700	CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
1701	CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
1702	CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo))
1703	CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))
1704	CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))
1705	CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo))
1706	CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo))
1707	CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo))
1708	CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))
1709	CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))
1710	CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))
1711	CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))
1712	CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo))
1713	zx_tmp_fi2d=0.
1714	IF(vars_defined) THEN
1715	DO k=1,klev
1716	!jyg<
1717	!! zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*qx(:,k,iq)
1718	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,iq-nqo)
1719	!>jyg
1720	ENDDO
1721	ENDIF
1722	CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
1723	endif
1724	ENDDO
1725	ENDIF
1726
1727	IF(.NOT.vars_defined) THEN
1728	!$OMP MASTER
1729	#ifndef CPP_IOIPSL_NO_OUTPUT
1730	DO iff=1,nfiles
1731	IF (clef_files(iff)) THEN
1732	CALL histend(nid_files(iff))
1733	ndex2d = 0
1734	ndex3d = 0
1735
1736	ENDIF ! clef_files
1737	ENDDO ! iff
1738	#endif
1739	#ifdef CPP_XIOS
1740	!On finalise l'initialisation:
1741	CALL wxios_closedef()
1742	#endif
1743
1744	!$OMP END MASTER
1745	!$OMP BARRIER
1746	vars_defined = .TRUE.
1747
1748	END IF
1749
1750	END DO
1751
1752	IF(vars_defined) THEN
1753	! On synchronise les fichiers pour IOIPSL
1754	#ifndef CPP_IOIPSL_NO_OUTPUT
1755	!$OMP MASTER
1756	DO iff=1,nfiles
1757	IF (ok_sync .AND. clef_files(iff)) THEN
1758	CALL histsync(nid_files(iff))
1759	ENDIF
1760	END DO
1761	!$OMP END MASTER
1762	#endif
1763	ENDIF
1764
1765	END SUBROUTINE phys_output_write
1766
1767	END MODULE phys_output_write_mod

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

Download in other formats:

Original Format