Context Navigation

phys_output_write_mod.F90 @ 4190

Last change on this file since 4190 was 4020, checked in by lguez, 3 years ago

Send 3 more fields to the ocean

Send 3 more fields to the ocean to compute CO2 flux at
ocean-atmosphere interface. The three fields are dter and dser, which
already existed, and a newly created field: dt_ds. So dter and dser
have to become state variables. The variable dt_ds of module
phys_state_var_mod is only allocated and defined if
activate_ocean_skin == 2 and type_ocean == "couple".

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: 104.7 KB

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

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

Download in other formats:

Original Format