Context Navigation

phys_output_write_mod.F90 @ 2057

Last change on this file since 2057 was 2057, checked in by Ehouarn Millour, 10 years ago
Preparatory stuff to fix and improve the slab ocean model. FC
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
File size: 56.4 KB

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

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

Download in other formats:

Original Format