Context Navigation

phys_output_write_mod.F90 @ 2371

Last change on this file since 2371 was 2366, checked in by oboucher, 10 years ago
Second batch of changes for diagnosing SW radiative diagnostics at the physics timestep resolution radsol is now computed in physiq.F90
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: 65.6 KB

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

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

Download in other formats:

Original Format