Context Navigation

phys_output_mod.F90 @ 1628

Last change on this file since 1628 was 1628, checked in by Laurent Fairhead, 12 years ago

Modifications concerning ALE and ALP.

new default values for the prescribed boundary layer values alp_bl_prescr=0.1 and ale_bl_prescr=4.
iflag_coupl=0 -> ale/p_bl=ale/p_bl_prescr

phys_output_mod.F90 and phys_output_write.h

Slight changes in the control of ALE and ALP outputs
(in particular ale_bl is in the outputs even with iflag_coupl=0.)

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 106.4 KB

Line
1	! $Id: phys_output_mod.F90 1628 2012-05-16 09:46:34Z fairhead $
2	!
3	! Abderrahmane 12 2007
4	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
5	!!! Ecreture des Sorties du modele dans les fichiers Netcdf :
6	! histmth.nc : moyennes mensuelles
7	! histday.nc : moyennes journalieres
8	! histhf.nc : moyennes toutes les 3 heures
9	! histins.nc : valeurs instantanees
10	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
11
12	MODULE phys_output_mod
13
14	IMPLICIT NONE
15
16	private histdef2d, histdef3d, conf_physoutputs
17
18
19	integer, parameter :: nfiles = 6
20	logical, dimension(nfiles), save :: clef_files
21	logical, dimension(nfiles), save :: clef_stations
22	integer, dimension(nfiles), save :: lev_files
23	integer, dimension(nfiles), save :: nid_files
24	integer, dimension(nfiles), save :: nnid_files
25	!!$OMP THREADPRIVATE(clef_files, clef_stations, lev_files,nid_files,nnid_files)
26	integer, dimension(nfiles), private, save :: nnhorim
27
28	integer, dimension(nfiles), private, save :: nhorim, nvertm
29	integer, dimension(nfiles), private, save :: nvertap, nvertbp, nvertAlt
30	! integer, dimension(nfiles), private, save :: nvertp0
31	real, dimension(nfiles), private, save :: zoutm
32	real, private, save :: zdtime
33	CHARACTER(len=20), dimension(nfiles), private, save :: type_ecri
34	!$OMP THREADPRIVATE(nhorim, nvertm, zoutm,zdtime,type_ecri)
35
36	! integer, save :: nid_hf3d
37
38	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
39	!! Definition pour chaque variable du niveau d ecriture dans chaque fichier
40	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!/ histmth, histday, histhf, histins /),'!!!!!!!!!!!!
41	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
42
43	integer, private:: levmin(nfiles) = 1
44	integer, private:: levmax(nfiles)
45
46	TYPE ctrl_out
47	integer,dimension(6) :: flag
48	character(len=20) :: name
49	END TYPE ctrl_out
50
51	!!! Comosentes de la coordonnee sigma-hybride
52	!!! Ap et Bp
53	type(ctrl_out),save :: o_Ahyb = ctrl_out((/ 1, 1, 1, 1, 1, 1 /), 'Ap')
54	type(ctrl_out),save :: o_Bhyb = ctrl_out((/ 1, 1, 1, 1, 1, 1 /), 'Bp')
55	type(ctrl_out),save :: o_Alt = ctrl_out((/ 1, 1, 1, 1, 1, 1 /), 'Alt')
56
57	!!! 1D
58	type(ctrl_out),save :: o_phis = ctrl_out((/ 1, 1, 10, 5, 1, 1 /), 'phis')
59	type(ctrl_out),save :: o_aire = ctrl_out((/ 1, 1, 10, 10, 1, 1 /),'aire')
60	type(ctrl_out),save :: o_contfracATM = ctrl_out((/ 10, 1, 1, 10, 10, 10 /),'contfracATM')
61	type(ctrl_out),save :: o_contfracOR = ctrl_out((/ 10, 1, 1, 10, 10, 10 /),'contfracOR')
62	type(ctrl_out),save :: o_aireTER = ctrl_out((/ 10, 10, 1, 10, 10, 10 /),'aireTER')
63
64	!!! 2D
65	type(ctrl_out),save :: o_flat = ctrl_out((/ 5, 1, 10, 10, 5, 10 /),'flat')
66	type(ctrl_out),save :: o_slp = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'slp')
67	type(ctrl_out),save :: o_tsol = ctrl_out((/ 1, 1, 1, 5, 10, 10 /),'tsol')
68	type(ctrl_out),save :: o_t2m = ctrl_out((/ 1, 1, 1, 5, 10, 10 /),'t2m')
69	type(ctrl_out),save :: o_t2m_min = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'t2m_min')
70	type(ctrl_out),save :: o_t2m_max = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'t2m_max')
71	type(ctrl_out),save,dimension(4) :: o_t2m_srf = (/ ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'t2m_ter'), &
72	ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'t2m_lic'), &
73	ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'t2m_oce'), &
74	ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'t2m_sic') /)
75
76	type(ctrl_out),save :: o_wind10m = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'wind10m')
77	type(ctrl_out),save :: o_wind10max = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'wind10max')
78	type(ctrl_out),save :: o_sicf = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'sicf')
79	type(ctrl_out),save :: o_q2m = ctrl_out((/ 1, 1, 1, 5, 10, 10 /),'q2m')
80	type(ctrl_out),save :: o_u10m = ctrl_out((/ 1, 1, 1, 5, 10, 10 /),'u10m')
81	type(ctrl_out),save :: o_v10m = ctrl_out((/ 1, 1, 1, 5, 10, 10 /),'v10m')
82	type(ctrl_out),save :: o_psol = ctrl_out((/ 1, 1, 1, 5, 10, 10 /),'psol')
83	type(ctrl_out),save :: o_qsurf = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'qsurf')
84
85	type(ctrl_out),save,dimension(4) :: o_u10m_srf = (/ ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'u10m_ter'), &
86	ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'u10m_lic'), &
87	ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'u10m_oce'), &
88	ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'u10m_sic') /)
89
90	type(ctrl_out),save,dimension(4) :: o_v10m_srf = (/ ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'v10m_ter'), &
91	ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'v10m_lic'), &
92	ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'v10m_oce'), &
93	ctrl_out((/ 10, 6, 10, 10, 10, 10 /),'v10m_sic') /)
94
95	type(ctrl_out),save :: o_qsol = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'qsol')
96
97	type(ctrl_out),save :: o_ndayrain = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'ndayrain')
98	type(ctrl_out),save :: o_precip = ctrl_out((/ 1, 1, 1, 10, 5, 10 /),'precip')
99	type(ctrl_out),save :: o_plul = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'plul')
100
101	type(ctrl_out),save :: o_pluc = ctrl_out((/ 1, 1, 1, 10, 5, 10 /),'pluc')
102	type(ctrl_out),save :: o_snow = ctrl_out((/ 1, 1, 10, 10, 5, 10 /),'snow')
103	type(ctrl_out),save :: o_evap = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'evap')
104	type(ctrl_out),save,dimension(4) :: o_evap_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'evap_ter'), &
105	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'evap_lic'), &
106	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'evap_oce'), &
107	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'evap_sic') /)
108	type(ctrl_out),save :: o_msnow = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'msnow')
109	type(ctrl_out),save :: o_fsnow = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsnow')
110
111	type(ctrl_out),save :: o_tops = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'tops')
112	type(ctrl_out),save :: o_tops0 = ctrl_out((/ 1, 5, 10, 10, 10, 10 /),'tops0')
113	type(ctrl_out),save :: o_topl = ctrl_out((/ 1, 1, 10, 5, 10, 10 /),'topl')
114	type(ctrl_out),save :: o_topl0 = ctrl_out((/ 1, 5, 10, 10, 10, 10 /),'topl0')
115	type(ctrl_out),save :: o_SWupTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'SWupTOA')
116	type(ctrl_out),save :: o_SWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'SWupTOAclr')
117	type(ctrl_out),save :: o_SWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'SWdnTOA')
118	type(ctrl_out),save :: o_SWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'SWdnTOAclr')
119	type(ctrl_out),save :: o_nettop = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'nettop')
120
121	type(ctrl_out),save :: o_SWup200 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'SWup200')
122	type(ctrl_out),save :: o_SWup200clr = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'SWup200clr')
123	type(ctrl_out),save :: o_SWdn200 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'SWdn200')
124	type(ctrl_out),save :: o_SWdn200clr = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'SWdn200clr')
125
126	! arajouter
127	! type(ctrl_out),save :: o_LWupTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWupTOA')
128	! type(ctrl_out),save :: o_LWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWupTOAclr')
129	! type(ctrl_out),save :: o_LWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWdnTOA')
130	! type(ctrl_out),save :: o_LWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10 /),'LWdnTOAclr')
131
132	type(ctrl_out),save :: o_LWup200 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'LWup200')
133	type(ctrl_out),save :: o_LWup200clr = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'LWup200clr')
134	type(ctrl_out),save :: o_LWdn200 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'LWdn200')
135	type(ctrl_out),save :: o_LWdn200clr = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'LWdn200clr')
136	type(ctrl_out),save :: o_sols = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'sols')
137	type(ctrl_out),save :: o_sols0 = ctrl_out((/ 1, 5, 10, 10, 10, 10 /),'sols0')
138	type(ctrl_out),save :: o_soll = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'soll')
139	type(ctrl_out),save :: o_soll0 = ctrl_out((/ 1, 5, 10, 10, 10, 10 /),'soll0')
140	type(ctrl_out),save :: o_radsol = ctrl_out((/ 1, 7, 10, 10, 10, 10 /),'radsol')
141	type(ctrl_out),save :: o_SWupSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10 /),'SWupSFC')
142	type(ctrl_out),save :: o_SWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10 /),'SWupSFCclr')
143	type(ctrl_out),save :: o_SWdnSFC = ctrl_out((/ 1, 1, 10, 10, 5, 10 /),'SWdnSFC')
144	type(ctrl_out),save :: o_SWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10 /),'SWdnSFCclr')
145	type(ctrl_out),save :: o_LWupSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10 /),'LWupSFC')
146	type(ctrl_out),save :: o_LWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10 /),'LWupSFCclr')
147	type(ctrl_out),save :: o_LWdnSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10 /),'LWdnSFC')
148	type(ctrl_out),save :: o_LWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10 /),'LWdnSFCclr')
149	type(ctrl_out),save :: o_bils = ctrl_out((/ 1, 2, 10, 5, 10, 10 /),'bils')
150	type(ctrl_out),save :: o_sens = ctrl_out((/ 1, 1, 10, 10, 5, 10 /),'sens')
151	type(ctrl_out),save :: o_fder = ctrl_out((/ 1, 2, 10, 10, 10, 10 /),'fder')
152	type(ctrl_out),save :: o_ffonte = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'ffonte')
153	type(ctrl_out),save :: o_fqcalving = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fqcalving')
154	type(ctrl_out),save :: o_fqfonte = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fqfonte')
155
156	type(ctrl_out),save :: o_taux = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'taux')
157	type(ctrl_out),save :: o_tauy = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'tauy')
158	type(ctrl_out),save,dimension(4) :: o_taux_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'taux_ter'), &
159	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'taux_lic'), &
160	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'taux_oce'), &
161	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'taux_sic') /)
162
163	type(ctrl_out),save,dimension(4) :: o_tauy_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tauy_ter'), &
164	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tauy_lic'), &
165	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tauy_oce'), &
166	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tauy_sic') /)
167
168
169	type(ctrl_out),save,dimension(4) :: o_pourc_srf = (/ ctrl_out((/ 1, 7, 10, 10, 10, 10 /),'pourc_ter'), &
170	ctrl_out((/ 1, 7, 10, 10, 10, 10 /),'pourc_lic'), &
171	ctrl_out((/ 1, 7, 10, 10, 10, 10 /),'pourc_oce'), &
172	ctrl_out((/ 1, 7, 10, 10, 10, 10 /),'pourc_sic') /)
173
174	type(ctrl_out),save,dimension(4) :: o_fract_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'fract_ter'), &
175	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'fract_lic'), &
176	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'fract_oce'), &
177	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'fract_sic') /)
178
179	type(ctrl_out),save,dimension(4) :: o_tsol_srf = (/ ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tsol_ter'), &
180	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tsol_lic'), &
181	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tsol_oce'), &
182	ctrl_out((/ 1, 6, 10, 10, 10, 10 /),'tsol_sic') /)
183
184	type(ctrl_out),save,dimension(4) :: o_sens_srf = (/ ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'sens_ter'), &
185	ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'sens_lic'), &
186	ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'sens_oce'), &
187	ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'sens_sic') /)
188
189	type(ctrl_out),save,dimension(4) :: o_lat_srf = (/ ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'lat_ter'), &
190	ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'lat_lic'), &
191	ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'lat_oce'), &
192	ctrl_out((/ 1, 6, 10, 7, 10, 10 /),'lat_sic') /)
193
194	type(ctrl_out),save,dimension(4) :: o_flw_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'flw_ter'), &
195	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'flw_lic'), &
196	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'flw_oce'), &
197	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'flw_sic') /)
198
199	type(ctrl_out),save,dimension(4) :: o_fsw_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsw_ter'), &
200	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsw_lic'), &
201	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsw_oce'), &
202	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fsw_sic') /)
203
204	type(ctrl_out),save,dimension(4) :: o_wbils_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbils_ter'), &
205	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbils_lic'), &
206	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbils_oce'), &
207	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbils_sic') /)
208
209	type(ctrl_out),save,dimension(4) :: o_wbilo_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbilo_ter'), &
210	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbilo_lic'), &
211	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbilo_oce'), &
212	ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbilo_sic') /)
213
214
215	type(ctrl_out),save :: o_cdrm = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'cdrm')
216	type(ctrl_out),save :: o_cdrh = ctrl_out((/ 1, 10, 10, 7, 10, 10 /),'cdrh')
217	type(ctrl_out),save :: o_cldl = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'cldl')
218	type(ctrl_out),save :: o_cldm = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'cldm')
219	type(ctrl_out),save :: o_cldh = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'cldh')
220	type(ctrl_out),save :: o_cldt = ctrl_out((/ 1, 1, 2, 10, 5, 10 /),'cldt')
221	type(ctrl_out),save :: o_cldq = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'cldq')
222	type(ctrl_out),save :: o_lwp = ctrl_out((/ 1, 5, 10, 10, 10, 10 /),'lwp')
223	type(ctrl_out),save :: o_iwp = ctrl_out((/ 1, 5, 10, 10, 10, 10 /),'iwp')
224	type(ctrl_out),save :: o_ue = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'ue')
225	type(ctrl_out),save :: o_ve = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'ve')
226	type(ctrl_out),save :: o_uq = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'uq')
227	type(ctrl_out),save :: o_vq = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'vq')
228
229	type(ctrl_out),save :: o_cape = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'cape')
230	type(ctrl_out),save :: o_pbase = ctrl_out((/ 1, 5, 10, 10, 10, 10 /),'pbase')
231	type(ctrl_out),save :: o_ptop = ctrl_out((/ 1, 5, 10, 10, 10, 10 /),'ptop')
232	type(ctrl_out),save :: o_fbase = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'fbase')
233	type(ctrl_out),save :: o_plcl = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'plcl')
234	type(ctrl_out),save :: o_plfc = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'plfc')
235	type(ctrl_out),save :: o_wbeff = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'wbeff')
236	type(ctrl_out),save :: o_prw = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'prw')
237
238	type(ctrl_out),save :: o_s_pblh = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_pblh')
239	type(ctrl_out),save :: o_s_pblt = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_pblt')
240	type(ctrl_out),save :: o_s_lcl = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_lcl')
241	type(ctrl_out),save :: o_s_therm = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_therm')
242	!IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
243	! type(ctrl_out),save :: o_s_capCL = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_capCL')
244	! type(ctrl_out),save :: o_s_oliqCL = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_oliqCL')
245	! type(ctrl_out),save :: o_s_cteiCL = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_cteiCL')
246	! type(ctrl_out),save :: o_s_trmb1 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_trmb1')
247	! type(ctrl_out),save :: o_s_trmb2 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_trmb2')
248	! type(ctrl_out),save :: o_s_trmb3 = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'s_trmb3')
249
250	type(ctrl_out),save :: o_slab_bils = ctrl_out((/ 1, 1, 10, 10, 10, 10 /),'slab_bils_oce')
251
252	type(ctrl_out),save :: o_ale_bl = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'ale_bl')
253	type(ctrl_out),save :: o_alp_bl = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'alp_bl')
254	type(ctrl_out),save :: o_ale_wk = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'ale_wk')
255	type(ctrl_out),save :: o_alp_wk = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'alp_wk')
256
257	type(ctrl_out),save :: o_ale = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'ale')
258	type(ctrl_out),save :: o_alp = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'alp')
259	type(ctrl_out),save :: o_cin = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'cin')
260	type(ctrl_out),save :: o_wape = ctrl_out((/ 1, 1, 1, 10, 10, 10 /),'wape')
261
262
263	! Champs interpolles sur des niveaux de pression ??? a faire correctement
264
265	type(ctrl_out),save,dimension(7) :: o_uSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u850'), &
266	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u700'), &
267	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u500'), &
268	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u200'), &
269	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u100'), &
270	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u50'), &
271	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'u10') /)
272
273
274	type(ctrl_out),save,dimension(7) :: o_vSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v850'), &
275	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v700'), &
276	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v500'), &
277	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v200'), &
278	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v100'), &
279	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v50'), &
280	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'v10') /)
281
282	type(ctrl_out),save,dimension(7) :: o_wSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w850'), &
283	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w700'), &
284	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w500'), &
285	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w200'), &
286	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w100'), &
287	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w50'), &
288	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'w10') /)
289
290	type(ctrl_out),save,dimension(7) :: o_tSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t850'), &
291	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t700'), &
292	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t500'), &
293	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t200'), &
294	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t100'), &
295	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t50'), &
296	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'t10') /)
297
298	type(ctrl_out),save,dimension(7) :: o_qSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q850'), &
299	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q700'), &
300	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q500'), &
301	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q200'), &
302	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q100'), &
303	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q50'), &
304	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'q10') /)
305
306	type(ctrl_out),save,dimension(7) :: o_zSTDlevs = (/ ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z850'), &
307	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z700'), &
308	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z500'), &
309	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z200'), &
310	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z100'), &
311	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z50'), &
312	ctrl_out((/ 1, 7, 7, 10, 10, 10 /),'z10') /)
313
314
315	type(ctrl_out),save :: o_t_oce_sic = ctrl_out((/ 1, 10, 10, 10, 10, 10 /),'t_oce_sic')
316
317	type(ctrl_out),save :: o_weakinv = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'weakinv')
318	type(ctrl_out),save :: o_dthmin = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'dthmin')
319	type(ctrl_out),save,dimension(4) :: o_u10_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'u10_ter'), &
320	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'u10_lic'), &
321	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'u10_oce'), &
322	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'u10_sic') /)
323
324	type(ctrl_out),save,dimension(4) :: o_v10_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'v10_ter'), &
325	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'v10_lic'), &
326	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'v10_oce'), &
327	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'v10_sic') /)
328
329	type(ctrl_out),save :: o_cldtau = ctrl_out((/ 10, 5, 10, 10, 10, 10 /),'cldtau')
330	type(ctrl_out),save :: o_cldemi = ctrl_out((/ 10, 5, 10, 10, 10, 10 /),'cldemi')
331	type(ctrl_out),save :: o_rh2m = ctrl_out((/ 5, 5, 10, 10, 10, 10 /),'rh2m')
332	type(ctrl_out),save :: o_rh2m_min = ctrl_out((/ 10, 5, 10, 10, 10, 10 /),'rh2m_min')
333	type(ctrl_out),save :: o_rh2m_max = ctrl_out((/ 10, 5, 10, 10, 10, 10 /),'rh2m_max')
334	type(ctrl_out),save :: o_qsat2m = ctrl_out((/ 10, 5, 10, 10, 10, 10 /),'qsat2m')
335	type(ctrl_out),save :: o_tpot = ctrl_out((/ 10, 5, 10, 10, 10, 10 /),'tpot')
336	type(ctrl_out),save :: o_tpote = ctrl_out((/ 10, 5, 10, 10, 10, 10 /),'tpote')
337	type(ctrl_out),save :: o_tke = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'tke ')
338	type(ctrl_out),save :: o_tke_max = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'tke_max')
339
340	type(ctrl_out),save,dimension(4) :: o_tke_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_ter'), &
341	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_lic'), &
342	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_oce'), &
343	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_sic') /)
344
345	type(ctrl_out),save,dimension(4) :: o_tke_max_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_max_ter'), &
346	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_max_lic'), &
347	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_max_oce'), &
348	ctrl_out((/ 10, 4, 10, 10, 10, 10 /),'tke_max_sic') /)
349
350	type(ctrl_out),save :: o_kz = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'kz')
351	type(ctrl_out),save :: o_kz_max = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'kz_max')
352	type(ctrl_out),save :: o_SWnetOR = ctrl_out((/ 10, 10, 2, 10, 10, 10 /),'SWnetOR')
353	type(ctrl_out),save :: o_SWdownOR = ctrl_out((/ 10, 10, 2, 10, 10, 10 /),'SWdownOR')
354	type(ctrl_out),save :: o_LWdownOR = ctrl_out((/ 10, 10, 2, 10, 10, 10 /),'LWdownOR')
355
356	type(ctrl_out),save :: o_snowl = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'snowl')
357	type(ctrl_out),save :: o_cape_max = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'cape_max')
358	type(ctrl_out),save :: o_solldown = ctrl_out((/ 10, 1, 10, 10, 10, 10 /),'solldown')
359
360	type(ctrl_out),save :: o_dtsvdfo = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'dtsvdfo')
361	type(ctrl_out),save :: o_dtsvdft = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'dtsvdft')
362	type(ctrl_out),save :: o_dtsvdfg = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'dtsvdfg')
363	type(ctrl_out),save :: o_dtsvdfi = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'dtsvdfi')
364	type(ctrl_out),save :: o_rugs = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'rugs')
365
366	type(ctrl_out),save :: o_topswad = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'topswad')
367	type(ctrl_out),save :: o_topswai = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'topswai')
368	type(ctrl_out),save :: o_solswad = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'solswad')
369	type(ctrl_out),save :: o_solswai = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'solswai')
370
371	type(ctrl_out),save,dimension(10) :: o_tausumaero = (/ ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASBCM'), &
372	ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASPOMM'), &
373	ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASSO4M'), &
374	ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_CSSO4M'), &
375	ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_SSSSM'), &
376	ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_ASSSM'), &
377	ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_CSSSM'), &
378	ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_CIDUSTM'), &
379	ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_AIBCM'), &
380	ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'OD550_AIPOMM') /)
381
382	type(ctrl_out),save :: o_od550aer = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'od550aer')
383	type(ctrl_out),save :: o_od865aer = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'od865aer')
384	type(ctrl_out),save :: o_absvisaer = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'absvisaer')
385	type(ctrl_out),save :: o_od550lt1aer = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'od550lt1aer')
386
387	type(ctrl_out),save :: o_sconcso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'sconcso4')
388	type(ctrl_out),save :: o_sconcoa = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'sconcoa')
389	type(ctrl_out),save :: o_sconcbc = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'sconcbc')
390	type(ctrl_out),save :: o_sconcss = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'sconcss')
391	type(ctrl_out),save :: o_sconcdust = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'sconcdust')
392	type(ctrl_out),save :: o_concso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'concso4')
393	type(ctrl_out),save :: o_concoa = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'concoa')
394	type(ctrl_out),save :: o_concbc = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'concbc')
395	type(ctrl_out),save :: o_concss = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'concss')
396	type(ctrl_out),save :: o_concdust = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'concdust')
397	type(ctrl_out),save :: o_loadso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'loadso4')
398	type(ctrl_out),save :: o_loadoa = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'loadoa')
399	type(ctrl_out),save :: o_loadbc = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'loadbc')
400	type(ctrl_out),save :: o_loadss = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'loadss')
401	type(ctrl_out),save :: o_loaddust = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'loaddust')
402
403	type(ctrl_out),save :: o_swtoaas_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swtoaas_nat')
404	type(ctrl_out),save :: o_swsrfas_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swsrfas_nat')
405	type(ctrl_out),save :: o_swtoacs_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swtoacs_nat')
406	type(ctrl_out),save :: o_swsrfcs_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swsrfcs_nat')
407
408	type(ctrl_out),save :: o_swtoaas_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swtoaas_ant')
409	type(ctrl_out),save :: o_swsrfas_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swsrfas_ant')
410	type(ctrl_out),save :: o_swtoacs_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swtoacs_ant')
411	type(ctrl_out),save :: o_swsrfcs_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swsrfcs_ant')
412
413	type(ctrl_out),save :: o_swtoacf_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swtoacf_nat')
414	type(ctrl_out),save :: o_swsrfcf_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swsrfcf_nat')
415	type(ctrl_out),save :: o_swtoacf_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swtoacf_ant')
416	type(ctrl_out),save :: o_swsrfcf_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swsrfcf_ant')
417	type(ctrl_out),save :: o_swtoacf_zero = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swtoacf_zero')
418	type(ctrl_out),save :: o_swsrfcf_zero = ctrl_out((/ 4, 6, 10, 10, 10, 10 /),'swsrfcf_zero')
419
420	type(ctrl_out),save :: o_cldncl = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'cldncl')
421	type(ctrl_out),save :: o_reffclwtop = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'reffclwtop')
422	type(ctrl_out),save :: o_cldnvi = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'cldnvi')
423	type(ctrl_out),save :: o_lcc = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'lcc')
424
425
426	!!!!!!!!!!!!!!!!!!!!!! 3D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
427	type(ctrl_out),save :: o_ec550aer = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'ec550aer')
428	type(ctrl_out),save :: o_lwcon = ctrl_out((/ 2, 5, 10, 10, 10, 10 /),'lwcon')
429	type(ctrl_out),save :: o_iwcon = ctrl_out((/ 2, 5, 10, 10, 10, 10 /),'iwcon')
430	type(ctrl_out),save :: o_temp = ctrl_out((/ 2, 3, 4, 10, 10, 10 /),'temp')
431	type(ctrl_out),save :: o_theta = ctrl_out((/ 2, 3, 4, 10, 10, 10 /),'theta')
432	type(ctrl_out),save :: o_ovap = ctrl_out((/ 2, 3, 4, 10, 10, 10 /),'ovap')
433	type(ctrl_out),save :: o_ovapinit = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'ovapinit')
434	type(ctrl_out),save :: o_wvapp = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'wvapp')
435	type(ctrl_out),save :: o_geop = ctrl_out((/ 2, 3, 10, 10, 10, 10 /),'geop')
436	type(ctrl_out),save :: o_vitu = ctrl_out((/ 2, 3, 4, 6, 10, 10 /),'vitu')
437	type(ctrl_out),save :: o_vitv = ctrl_out((/ 2, 3, 4, 6, 10, 10 /),'vitv')
438	type(ctrl_out),save :: o_vitw = ctrl_out((/ 2, 3, 10, 6, 10, 10 /),'vitw')
439	type(ctrl_out),save :: o_pres = ctrl_out((/ 2, 3, 10, 10, 10, 10 /),'pres')
440	type(ctrl_out),save :: o_paprs = ctrl_out((/ 2, 3, 10, 10, 10, 10 /),'paprs')
441	type(ctrl_out),save :: o_mass = ctrl_out((/ 2, 3, 10, 10, 10, 10 /),'mass')
442	type(ctrl_out),save :: o_zfull = ctrl_out((/ 2, 3, 10, 10, 10, 10 /),'zfull')
443	type(ctrl_out),save :: o_zhalf = ctrl_out((/ 2, 3, 10, 10, 10, 10 /),'zhalf')
444	type(ctrl_out),save :: o_rneb = ctrl_out((/ 2, 5, 10, 10, 10, 10 /),'rneb')
445	type(ctrl_out),save :: o_rnebcon = ctrl_out((/ 2, 5, 10, 10, 10, 10 /),'rnebcon')
446	type(ctrl_out),save :: o_rhum = ctrl_out((/ 2, 5, 10, 10, 10, 10 /),'rhum')
447	type(ctrl_out),save :: o_ozone = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'ozone')
448	type(ctrl_out),save :: o_ozone_light = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'ozone_daylight')
449	type(ctrl_out),save :: o_upwd = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'upwd')
450	type(ctrl_out),save :: o_dtphy = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'dtphy')
451	type(ctrl_out),save :: o_dqphy = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'dqphy')
452	type(ctrl_out),save :: o_pr_con_l = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'pr_con_l')
453	type(ctrl_out),save :: o_pr_con_i = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'pr_con_i')
454	type(ctrl_out),save :: o_pr_lsc_l = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'pr_lsc_l')
455	type(ctrl_out),save :: o_pr_lsc_i = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'pr_lsc_i')
456	type(ctrl_out),save :: o_re = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'re')
457	type(ctrl_out),save :: o_fl = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'fl')
458	type(ctrl_out),save :: o_scdnc = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'scdnc')
459	type(ctrl_out),save :: o_reffclws = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'reffclws')
460	type(ctrl_out),save :: o_reffclwc = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'reffclwc')
461	type(ctrl_out),save :: o_lcc3d = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'lcc3d')
462	type(ctrl_out),save :: o_lcc3dcon = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'lcc3dcon')
463	type(ctrl_out),save :: o_lcc3dstra = ctrl_out((/ 2, 6, 10, 10, 10, 10 /),'lcc3dstra')
464	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
465
466	type(ctrl_out),save,dimension(4) :: o_albe_srf = (/ ctrl_out((/ 3, 7, 10, 7, 10, 10 /),'albe_ter'), &
467	ctrl_out((/ 3, 7, 10, 7, 10, 10 /),'albe_lic'), &
468	ctrl_out((/ 3, 7, 10, 7, 10, 10 /),'albe_oce'), &
469	ctrl_out((/ 3, 7, 10, 7, 10, 10 /),'albe_sic') /)
470
471	type(ctrl_out),save,dimension(4) :: o_ages_srf = (/ ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'ages_ter'), &
472	ctrl_out((/ 3, 10, 10, 10, 10, 10 /),'ages_lic'), &
473	ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'ages_oce'), &
474	ctrl_out((/ 3, 10, 10, 10, 10, 10 /),'ages_sic') /)
475
476	type(ctrl_out),save,dimension(4) :: o_rugs_srf = (/ ctrl_out((/ 3, 6, 10, 10, 10, 10 /),'rugs_ter'), &
477	ctrl_out((/ 3, 6, 10, 10, 10, 10 /),'rugs_lic'), &
478	ctrl_out((/ 3, 6, 10, 10, 10, 10 /),'rugs_oce'), &
479	ctrl_out((/ 3, 6, 10, 10, 10, 10 /),'rugs_sic') /)
480
481	type(ctrl_out),save :: o_alb1 = ctrl_out((/ 3, 10, 10, 10, 10, 10 /),'alb1')
482	type(ctrl_out),save :: o_alb2 = ctrl_out((/ 3, 10, 10, 10, 10, 10 /),'alb2')
483
484	type(ctrl_out),save :: o_clwcon = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'clwcon')
485	type(ctrl_out),save :: o_Ma = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'Ma')
486	type(ctrl_out),save :: o_dnwd = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dnwd')
487	type(ctrl_out),save :: o_dnwd0 = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dnwd0')
488	type(ctrl_out),save :: o_mc = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'mc')
489	type(ctrl_out),save :: o_ftime_con = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'ftime_con')
490	type(ctrl_out),save :: o_dtdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtdyn')
491	type(ctrl_out),save :: o_dqdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dqdyn')
492	type(ctrl_out),save :: o_dudyn = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dudyn') !AXC
493	type(ctrl_out),save :: o_dvdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dvdyn') !AXC
494	type(ctrl_out),save :: o_dtcon = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtcon')
495	type(ctrl_out),save :: o_ducon = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'ducon')
496	type(ctrl_out),save :: o_dqcon = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dqcon')
497	type(ctrl_out),save :: o_dtwak = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'dtwak')
498	type(ctrl_out),save :: o_dqwak = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'dqwak')
499	type(ctrl_out),save :: o_wake_h = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'wake_h')
500	type(ctrl_out),save :: o_wake_s = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'wake_s')
501	type(ctrl_out),save :: o_wake_deltat = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'wake_deltat')
502	type(ctrl_out),save :: o_wake_deltaq = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'wake_deltaq')
503	type(ctrl_out),save :: o_wake_omg = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'wake_omg')
504	type(ctrl_out),save :: o_Vprecip = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'Vprecip')
505	type(ctrl_out),save :: o_ftd = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'ftd')
506	type(ctrl_out),save :: o_fqd = ctrl_out((/ 4, 5, 10, 10, 10, 10 /),'fqd')
507	type(ctrl_out),save :: o_dtlsc = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtlsc')
508	type(ctrl_out),save :: o_dtlschr = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtlschr')
509	type(ctrl_out),save :: o_dqlsc = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dqlsc')
510	type(ctrl_out),save :: o_dtvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtvdf')
511	type(ctrl_out),save :: o_dqvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dqvdf')
512	type(ctrl_out),save :: o_dteva = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dteva')
513	type(ctrl_out),save :: o_dqeva = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dqeva')
514
515	!!!!!!!!!!!!!!!! Specifique thermiques
516	type(ctrl_out),save :: o_dqlscth = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'dqlscth')
517	type(ctrl_out),save :: o_dqlscst = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'dqlscst')
518	type(ctrl_out),save :: o_dtlscth = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'dtlscth')
519	type(ctrl_out),save :: o_dtlscst = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'dtlscst')
520	type(ctrl_out),save :: o_plulth = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'plulth')
521	type(ctrl_out),save :: o_plulst = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'plulst')
522	type(ctrl_out),save :: o_lmaxth = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'lmaxth')
523	type(ctrl_out),save :: o_ptconvth = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'ptconvth')
524	!!!!!!!!!!!!!!!!!!!!!!!!
525
526
527	type(ctrl_out),save :: o_ptconv = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'ptconv')
528	type(ctrl_out),save :: o_ratqs = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'ratqs')
529	type(ctrl_out),save :: o_dtthe = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtthe')
530	type(ctrl_out),save :: o_f_th = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'f_th')
531	type(ctrl_out),save :: o_e_th = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'e_th')
532	type(ctrl_out),save :: o_w_th = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'w_th')
533	type(ctrl_out),save :: o_lambda_th = ctrl_out((/ 10, 10, 10, 10, 10, 10 /),'lambda_th')
534	type(ctrl_out),save :: o_ftime_th = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'ftime_th')
535	type(ctrl_out),save :: o_q_th = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'q_th')
536	type(ctrl_out),save :: o_a_th = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'a_th')
537	type(ctrl_out),save :: o_d_th = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'d_th')
538	type(ctrl_out),save :: o_f0_th = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'f0_th')
539	type(ctrl_out),save :: o_zmax_th = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'zmax_th')
540	type(ctrl_out),save :: o_dqthe = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dqthe')
541	type(ctrl_out),save :: o_dtajs = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtajs')
542	type(ctrl_out),save :: o_dqajs = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dqajs')
543	type(ctrl_out),save :: o_dtswr = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtswr')
544	type(ctrl_out),save :: o_dtsw0 = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtsw0')
545	type(ctrl_out),save :: o_dtlwr = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtlwr')
546	type(ctrl_out),save :: o_dtlw0 = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtlw0')
547	type(ctrl_out),save :: o_dtec = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtec')
548	type(ctrl_out),save :: o_duvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'duvdf')
549	type(ctrl_out),save :: o_dvvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dvvdf')
550	type(ctrl_out),save :: o_duoro = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'duoro')
551	type(ctrl_out),save :: o_dvoro = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dvoro')
552	type(ctrl_out),save :: o_dulif = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dulif')
553	type(ctrl_out),save :: o_dvlif = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dvlif')
554	type(ctrl_out),save :: o_duhin = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'duhin')
555	type(ctrl_out),save :: o_dvhin = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dvhin')
556	type(ctrl_out),save :: o_dtoro = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtoro')
557	type(ctrl_out),save :: o_dtlif = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dtlif')
558	type(ctrl_out),save :: o_dthin = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dthin')
559
560	type(ctrl_out),save,allocatable :: o_trac(:)
561
562	type(ctrl_out),save :: o_rsu = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'rsu')
563	type(ctrl_out),save :: o_rsd = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'rsd')
564	type(ctrl_out),save :: o_rlu = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'rlu')
565	type(ctrl_out),save :: o_rld = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'rld')
566	type(ctrl_out),save :: o_rsucs = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'rsucs')
567	type(ctrl_out),save :: o_rsdcs = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'rsdcs')
568	type(ctrl_out),save :: o_rlucs = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'rlucs')
569	type(ctrl_out),save :: o_rldcs = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'rldcs')
570
571	type(ctrl_out),save :: o_tnt = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'tnt')
572	type(ctrl_out),save :: o_tntc = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'tntc')
573	type(ctrl_out),save :: o_tntr = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'tntr')
574	type(ctrl_out),save :: o_tntscpbl = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'tntscpbl')
575
576	type(ctrl_out),save :: o_tnhus = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'tnhus')
577	type(ctrl_out),save :: o_tnhusc = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'tnhusc')
578	type(ctrl_out),save :: o_tnhusscpbl = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'tnhusscpbl')
579
580	type(ctrl_out),save :: o_evu = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'evu')
581
582	type(ctrl_out),save :: o_h2o = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'h2o')
583
584	type(ctrl_out),save :: o_mcd = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'mcd')
585	type(ctrl_out),save :: o_dmc = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'dmc')
586	type(ctrl_out),save :: o_ref_liq = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'ref_liq')
587	type(ctrl_out),save :: o_ref_ice = ctrl_out((/ 4, 10, 10, 10, 10, 10 /),'ref_ice')
588
589	type(ctrl_out),save :: o_rsut4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rsut4co2')
590	type(ctrl_out),save :: o_rlut4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rlut4co2')
591	type(ctrl_out),save :: o_rsutcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rsutcs4co2')
592	type(ctrl_out),save :: o_rlutcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rlutcs4co2')
593
594	type(ctrl_out),save :: o_rsu4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rsu4co2')
595	type(ctrl_out),save :: o_rlu4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rlu4co2')
596	type(ctrl_out),save :: o_rsucs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rsucs4co2')
597	type(ctrl_out),save :: o_rlucs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rlucs4co2')
598	type(ctrl_out),save :: o_rsd4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rsd4co2')
599	type(ctrl_out),save :: o_rld4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rld4co2')
600	type(ctrl_out),save :: o_rsdcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rsdcs4co2')
601	type(ctrl_out),save :: o_rldcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10 /),'rldcs4co2')
602
603
604	CONTAINS
605
606	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
607	!!!!!!!!! Ouverture des fichier et definition des variable de sortie !!!!!!!!
608	!! histbeg, histvert et histdef
609	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
610
611	SUBROUTINE phys_output_open(rlon,rlat,pim,tabij,ipt,jpt,plon,plat, &
612	jjmp1,nlevSTD,clevSTD,nbteta, &
613	ctetaSTD, dtime, ok_veget, &
614	type_ocean, iflag_pbl,ok_mensuel,ok_journe, &
615	ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz, &
616	phys_out_filestations, &
617	new_aod, aerosol_couple)
618
619	USE iophy
620	USE dimphy
621	USE infotrac
622	USE ioipsl
623	USE mod_phys_lmdz_para
624	USE aero_mod, only : naero_spc,name_aero
625
626	IMPLICIT NONE
627	include "dimensions.h"
628	include "temps.h"
629	include "indicesol.h"
630	include "clesphys.h"
631	include "thermcell.h"
632	include "comvert.h"
633	include "iniprint.h"
634
635	real,dimension(klon),intent(in) :: rlon
636	real,dimension(klon),intent(in) :: rlat
637	integer, intent(in) :: pim
638	INTEGER, DIMENSION(pim) :: tabij
639	INTEGER,dimension(pim), intent(in) :: ipt, jpt
640	REAL,dimension(pim), intent(in) :: plat, plon
641	REAL,dimension(pim,2) :: plat_bounds, plon_bounds
642
643	integer :: jjmp1
644	integer :: nbteta, nlevSTD, radpas
645	logical :: ok_mensuel, ok_journe, ok_hf, ok_instan
646	logical :: ok_LES,ok_ade,ok_aie
647	logical :: new_aod, aerosol_couple
648	integer, intent(in):: read_climoz ! read ozone climatology
649	! Allowed values are 0, 1 and 2
650	! 0: do not read an ozone climatology
651	! 1: read a single ozone climatology that will be used day and night
652	! 2: read two ozone climatologies, the average day and night
653	! climatology and the daylight climatology
654
655	real :: dtime
656	integer :: idayref
657	real :: zjulian
658	real, dimension(klev) :: Ahyb, Bhyb, Alt
659	character(len=4), dimension(nlevSTD) :: clevSTD
660	integer :: nsrf, k, iq, iiq, iff, i, j, ilev
661	integer :: naero
662	logical :: ok_veget
663	integer :: iflag_pbl
664	CHARACTER(len=4) :: bb2
665	CHARACTER(len=2) :: bb3
666	character(len=6) :: type_ocean
667	CHARACTER(len=3) :: ctetaSTD(nbteta)
668	real, dimension(nfiles) :: ecrit_files
669	CHARACTER(len=20), dimension(nfiles) :: phys_out_filenames
670	INTEGER, dimension(iim*jjmp1) :: ndex2d
671	INTEGER, dimension(iimjjmp1klev) :: ndex3d
672	integer :: imin_ins, imax_ins
673	integer :: jmin_ins, jmax_ins
674	integer, dimension(nfiles) :: phys_out_levmin, phys_out_levmax
675	integer, dimension(nfiles) :: phys_out_filelevels
676	CHARACTER(len=20), dimension(nfiles) :: type_ecri_files, phys_out_filetypes
677	character(len=20), dimension(nfiles) :: chtimestep = (/ 'DefFreq', 'DefFreq','DefFreq', 'DefFreq', 'DefFreq', 'DefFreq' /)
678	logical, dimension(nfiles) :: phys_out_filekeys
679	logical, dimension(nfiles) :: phys_out_filestations
680
681	!!!!!!!!!! stockage dans une region limitee pour chaque fichier !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
682	! entre [phys_out_lonmin,phys_out_lonmax] et [phys_out_latmin,phys_out_latmax]
683
684	logical, dimension(nfiles), save :: phys_out_regfkey = (/ .false., .false., .false., .false., .false., .false. /)
685	real, dimension(nfiles), save :: phys_out_lonmin = (/ -180., -180., -180., -180., -180., -180. /)
686	real, dimension(nfiles), save :: phys_out_lonmax = (/ 180., 180., 180., 180., 180., 180. /)
687	real, dimension(nfiles), save :: phys_out_latmin = (/ -90., -90., -90., -90., -90., -90. /)
688	real, dimension(nfiles), save :: phys_out_latmax = (/ 90., 90., 90., 90., 90., 90. /)
689
690	write(lunout,*) 'Debut phys_output_mod.F90'
691	! Initialisations (Valeurs par defaut
692
693	if (.not. allocated(o_trac)) ALLOCATE(o_trac(nqtot))
694
695	levmax = (/ klev, klev, klev, klev, klev, klev /)
696
697	phys_out_filenames(1) = 'histmth'
698	phys_out_filenames(2) = 'histday'
699	phys_out_filenames(3) = 'histhf'
700	phys_out_filenames(4) = 'histins'
701	phys_out_filenames(5) = 'histLES'
702	phys_out_filenames(6) = 'histstn'
703
704	type_ecri(1) = 'ave(X)'
705	type_ecri(2) = 'ave(X)'
706	type_ecri(3) = 'ave(X)'
707	type_ecri(4) = 'inst(X)'
708	type_ecri(5) = 'ave(X)'
709	type_ecri(6) = 'inst(X)'
710
711	clef_files(1) = ok_mensuel
712	clef_files(2) = ok_journe
713	clef_files(3) = ok_hf
714	clef_files(4) = ok_instan
715	clef_files(5) = ok_LES
716	clef_files(6) = ok_instan
717
718	!sortir des fichiers "stations" si clef_stations(:)=.TRUE.
719	clef_stations(1) = .FALSE.
720	clef_stations(2) = .FALSE.
721	clef_stations(3) = .FALSE.
722	clef_stations(4) = .FALSE.
723	clef_stations(5) = .FALSE.
724	clef_stations(6) = .FALSE.
725
726	lev_files(1) = lev_histmth
727	lev_files(2) = lev_histday
728	lev_files(3) = lev_histhf
729	lev_files(4) = lev_histins
730	lev_files(5) = lev_histLES
731	lev_files(6) = lev_histins
732
733	ecrit_files(1) = ecrit_mth
734	ecrit_files(2) = ecrit_day
735	ecrit_files(3) = ecrit_hf
736	ecrit_files(4) = ecrit_ins
737	ecrit_files(5) = ecrit_LES
738	ecrit_files(6) = ecrit_ins
739
740	!! Lectures des parametres de sorties dans physiq.def
741
742	call getin('phys_out_regfkey',phys_out_regfkey)
743	call getin('phys_out_lonmin',phys_out_lonmin)
744	call getin('phys_out_lonmax',phys_out_lonmax)
745	call getin('phys_out_latmin',phys_out_latmin)
746	call getin('phys_out_latmax',phys_out_latmax)
747	phys_out_levmin(:)=levmin(:)
748	call getin('phys_out_levmin',levmin)
749	phys_out_levmax(:)=levmax(:)
750	call getin('phys_out_levmax',levmax)
751	call getin('phys_out_filenames',phys_out_filenames)
752	phys_out_filekeys(:)=clef_files(:)
753	call getin('phys_out_filekeys',clef_files)
754	phys_out_filestations(:)=clef_stations(:)
755	call getin('phys_out_filestations',clef_stations)
756	phys_out_filelevels(:)=lev_files(:)
757	call getin('phys_out_filelevels',lev_files)
758	call getin('phys_out_filetimesteps',chtimestep)
759	phys_out_filetypes(:)=type_ecri(:)
760	call getin('phys_out_filetypes',type_ecri)
761
762	type_ecri_files(:)=type_ecri(:)
763
764	write(lunout,*)'phys_out_lonmin=',phys_out_lonmin
765	write(lunout,*)'phys_out_lonmax=',phys_out_lonmax
766	write(lunout,*)'phys_out_latmin=',phys_out_latmin
767	write(lunout,*)'phys_out_latmax=',phys_out_latmax
768	write(lunout,*)'phys_out_filenames=',phys_out_filenames
769	write(lunout,*)'phys_out_filetypes=',type_ecri
770	write(lunout,*)'phys_out_filekeys=',clef_files
771	write(lunout,*)'phys_out_filestations=',clef_stations
772	write(lunout,*)'phys_out_filelevels=',lev_files
773
774	!!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
775	! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !!
776	! Appel des histbeg pour definir les variables (nom, moy ou inst, freq de sortie ..
777	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
778
779	zdtime = dtime ! Frequence ou l on moyenne
780
781	! Calcul des Ahyb, Bhyb et Alt
782	do k=1,klev
783	Ahyb(k)=(ap(k)+ap(k+1))/2.
784	Bhyb(k)=(bp(k)+bp(k+1))/2.
785	Alt(k)=log(preff/presnivs(k))*8.
786	enddo
787	! if(prt_level.ge.1) then
788	write(lunout,*)'Ap Hybrid = ',Ahyb(1:klev)
789	write(lunout,*)'Bp Hybrid = ',Bhyb(1:klev)
790	write(lunout,*)'Alt approx des couches pour une haut d echelle de 8km = ',Alt(1:klev)
791	! endif
792	DO iff=1,nfiles
793
794	IF (clef_files(iff)) THEN
795
796	if ( chtimestep(iff).eq.'DefFreq' ) then
797	! Par defaut ecrit_files = (ecrit_mensuel ecrit_jour ecrit_hf ...)*86400.
798	ecrit_files(iff)=ecrit_files(iff)*86400.
799	else
800	call convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff))
801	endif
802	write(lunout,*)'ecrit_files(',iff,')= ',ecrit_files(iff)
803
804	zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde
805
806	idayref = day_ref
807	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
808
809	!!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
810	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
811	if (phys_out_regfkey(iff)) then
812
813	imin_ins=1
814	imax_ins=iim
815	jmin_ins=1
816	jmax_ins=jjmp1
817
818	! correction abderr
819	do i=1,iim
820	write(lunout,*)'io_lon(i)=',io_lon(i)
821	if (io_lon(i).le.phys_out_lonmin(iff)) imin_ins=i
822	if (io_lon(i).le.phys_out_lonmax(iff)) imax_ins=i+1
823	enddo
824
825	do j=1,jjmp1
826	write(lunout,*)'io_lat(j)=',io_lat(j)
827	if (io_lat(j).ge.phys_out_latmin(iff)) jmax_ins=j+1
828	if (io_lat(j).ge.phys_out_latmax(iff)) jmin_ins=j
829	enddo
830
831	write(lunout,*)'On stoke le fichier histoire numero ',iff,' sur ', &
832	imin_ins,imax_ins,jmin_ins,jmax_ins
833	write(lunout,*)'longitudes : ', &
834	io_lon(imin_ins),io_lon(imax_ins), &
835	'latitudes : ', &
836	io_lat(jmax_ins),io_lat(jmin_ins)
837
838	CALL histbeg(phys_out_filenames(iff),iim,io_lon,jjmp1,io_lat, &
839	imin_ins,imax_ins-imin_ins+1, &
840	jmin_ins,jmax_ins-jmin_ins+1, &
841	itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
842	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
843	!IM fichiers stations
844	else if (clef_stations(iff)) THEN
845
846	write(lunout,*)'phys_output_mod phys_out_filenames=',phys_out_filenames(iff)
847
848	call histbeg_phy_all(rlon,rlat,pim,tabij,ipt,jpt,plon,plat,plon_bounds,plat_bounds, &
849	phys_out_filenames(iff), &
850	itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
851	else
852	CALL histbeg_phy(phys_out_filenames(iff),itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
853	endif
854
855	CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "Pa", &
856	levmax(iff) - levmin(iff) + 1, &
857	presnivs(levmin(iff):levmax(iff)), nvertm(iff),"down")
858
859	!!!!!!!!!!!!! Traitement des champs 3D pour histhf !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
860	!!!!!!!!!!!!!!! A Revoir plus tard !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
861	! IF (iff.eq.3.and.lev_files(iff).ge.4) THEN
862	! CALL histbeg_phy("histhf3d",itau_phy, &
863	! & zjulian, dtime, &
864	! & nhorim, nid_hf3d)
865
866	! CALL histvert(nid_hf3d, "presnivs", &
867	! & "Vertical levels", "mb", &
868	! & klev, presnivs/100., nvertm)
869	! ENDIF
870	!
871	!!!! Composentes de la coordonnee sigma-hybride
872	CALL histvert(nid_files(iff), "Ahyb","Ahyb comp of Hyb Cord ", "Pa", &
873	levmax(iff) - levmin(iff) + 1,Ahyb,nvertap(iff))
874
875	CALL histvert(nid_files(iff), "Bhyb","Bhyb comp of Hyb Cord", " ", &
876	levmax(iff) - levmin(iff) + 1,Bhyb,nvertbp(iff))
877
878	CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", &
879	levmax(iff) - levmin(iff) + 1,Alt,nvertAlt(iff))
880
881	! CALL histvert(nid_files(iff), "preff","Reference pressure", "Pa", &
882	! 1,preff,nvertp0(iff))
883	!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
884	IF (.NOT.clef_stations(iff)) THEN
885	!
886	!IM: there is no way to have one single value in a netcdf file
887	!
888	type_ecri(1) = 'once'
889	type_ecri(2) = 'once'
890	type_ecri(3) = 'once'
891	type_ecri(4) = 'once'
892	type_ecri(5) = 'once'
893	type_ecri(6) = 'once'
894	CALL histdef2d(iff,clef_stations(iff),o_aire%flag,o_aire%name,"Grid area", "-")
895	CALL histdef2d(iff,clef_stations(iff),o_contfracATM%flag,o_contfracATM%name,"% sfce ter+lic", "-")
896	ENDIF
897	type_ecri(:) = type_ecri_files(:)
898
899	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
900	CALL histdef2d(iff,clef_stations(iff),o_phis%flag,o_phis%name,"Surface geop.height", "m2/s2" )
901	CALL histdef2d(iff,clef_stations(iff),o_contfracOR%flag,o_contfracOR%name,"% sfce terre OR", "-" )
902	CALL histdef2d(iff,clef_stations(iff),o_aireTER%flag,o_aireTER%name,"Grid area CONT", "-" )
903	CALL histdef2d(iff,clef_stations(iff),o_flat%flag,o_flat%name, "Latent heat flux", "W/m2")
904	CALL histdef2d(iff,clef_stations(iff),o_slp%flag,o_slp%name, "Sea Level Pressure", "Pa" )
905	CALL histdef2d(iff,clef_stations(iff),o_tsol%flag,o_tsol%name, "Surface Temperature", "K")
906	CALL histdef2d(iff,clef_stations(iff),o_t2m%flag,o_t2m%name, "Temperature 2m", "K" )
907	IF (.NOT.clef_stations(iff)) THEN
908	!
909	!IM: there is no way to have one single value in a netcdf file
910	!
911	type_ecri(1) = 't_min(X)'
912	type_ecri(2) = 't_min(X)'
913	type_ecri(3) = 't_min(X)'
914	type_ecri(4) = 't_min(X)'
915	type_ecri(5) = 't_min(X)'
916	type_ecri(6) = 't_min(X)'
917	CALL histdef2d(iff,clef_stations(iff),o_t2m_min%flag,o_t2m_min%name, "Temp 2m min", "K" )
918	type_ecri(1) = 't_max(X)'
919	type_ecri(2) = 't_max(X)'
920	type_ecri(3) = 't_max(X)'
921	type_ecri(4) = 't_max(X)'
922	type_ecri(5) = 't_max(X)'
923	type_ecri(6) = 't_max(X)'
924	CALL histdef2d(iff,clef_stations(iff),o_t2m_max%flag,o_t2m_max%name, "Temp 2m max", "K" )
925	ENDIF
926	type_ecri(:) = type_ecri_files(:)
927	CALL histdef2d(iff,clef_stations(iff),o_wind10m%flag,o_wind10m%name, "10-m wind speed", "m/s")
928	CALL histdef2d(iff,clef_stations(iff),o_wind10max%flag,o_wind10max%name, "10m wind speed max", "m/s")
929	CALL histdef2d(iff,clef_stations(iff),o_sicf%flag,o_sicf%name, "Sea-ice fraction", "-" )
930	CALL histdef2d(iff,clef_stations(iff),o_q2m%flag,o_q2m%name, "Specific humidity 2m", "kg/kg")
931	CALL histdef2d(iff,clef_stations(iff),o_u10m%flag,o_u10m%name, "Vent zonal 10m", "m/s" )
932	CALL histdef2d(iff,clef_stations(iff),o_v10m%flag,o_v10m%name, "Vent meridien 10m", "m/s")
933	CALL histdef2d(iff,clef_stations(iff),o_psol%flag,o_psol%name, "Surface Pressure", "Pa" )
934	CALL histdef2d(iff,clef_stations(iff),o_qsurf%flag,o_qsurf%name, "Surface Air humidity", "kg/kg")
935
936	if (.not. ok_veget) then
937	CALL histdef2d(iff,clef_stations(iff),o_qsol%flag,o_qsol%name, "Soil watter content", "mm" )
938	endif
939
940	CALL histdef2d(iff,clef_stations(iff),o_ndayrain%flag,o_ndayrain%name, "Number of dayrain(liq+sol)", "-")
941	CALL histdef2d(iff,clef_stations(iff),o_precip%flag,o_precip%name, "Precip Totale liq+sol", "kg/(s*m2)" )
942	CALL histdef2d(iff,clef_stations(iff),o_plul%flag,o_plul%name, "Large-scale Precip.", "kg/(s*m2)")
943	CALL histdef2d(iff,clef_stations(iff),o_pluc%flag,o_pluc%name, "Convective Precip.", "kg/(s*m2)")
944	CALL histdef2d(iff,clef_stations(iff),o_snow%flag,o_snow%name, "Snow fall", "kg/(s*m2)" )
945	CALL histdef2d(iff,clef_stations(iff),o_msnow%flag,o_msnow%name, "Surface snow amount", "kg/m2" )
946	CALL histdef2d(iff,clef_stations(iff),o_fsnow%flag,o_fsnow%name, "Surface snow area fraction", "-" )
947	CALL histdef2d(iff,clef_stations(iff),o_evap%flag,o_evap%name, "Evaporat", "kg/(s*m2)" )
948	CALL histdef2d(iff,clef_stations(iff),o_tops%flag,o_tops%name, "Solar rad. at TOA", "W/m2")
949	CALL histdef2d(iff,clef_stations(iff),o_tops0%flag,o_tops0%name, "CS Solar rad. at TOA", "W/m2")
950	CALL histdef2d(iff,clef_stations(iff),o_topl%flag,o_topl%name, "IR rad. at TOA", "W/m2" )
951	CALL histdef2d(iff,clef_stations(iff),o_topl0%flag,o_topl0%name, "IR rad. at TOA", "W/m2")
952	CALL histdef2d(iff,clef_stations(iff),o_SWupTOA%flag,o_SWupTOA%name, "SWup at TOA", "W/m2")
953	CALL histdef2d(iff,clef_stations(iff),o_SWupTOAclr%flag,o_SWupTOAclr%name, "SWup clear sky at TOA", "W/m2")
954	CALL histdef2d(iff,clef_stations(iff),o_SWdnTOA%flag,o_SWdnTOA%name, "SWdn at TOA", "W/m2" )
955	CALL histdef2d(iff,clef_stations(iff),o_SWdnTOAclr%flag,o_SWdnTOAclr%name, "SWdn clear sky at TOA", "W/m2")
956	CALL histdef2d(iff,clef_stations(iff),o_nettop%flag,o_nettop%name, "Net dn radiatif flux at TOA", "W/m2")
957	CALL histdef2d(iff,clef_stations(iff),o_SWup200%flag,o_SWup200%name, "SWup at 200mb", "W/m2" )
958	CALL histdef2d(iff,clef_stations(iff),o_SWup200clr%flag,o_SWup200clr%name, "SWup clear sky at 200mb", "W/m2")
959	CALL histdef2d(iff,clef_stations(iff),o_SWdn200%flag,o_SWdn200%name, "SWdn at 200mb", "W/m2" )
960	CALL histdef2d(iff,clef_stations(iff),o_SWdn200clr%flag,o_SWdn200clr%name, "SWdn clear sky at 200mb", "W/m2")
961	CALL histdef2d(iff,clef_stations(iff),o_LWup200%flag,o_LWup200%name, "LWup at 200mb", "W/m2")
962	CALL histdef2d(iff,clef_stations(iff),o_LWup200clr%flag,o_LWup200clr%name, "LWup clear sky at 200mb", "W/m2")
963	CALL histdef2d(iff,clef_stations(iff),o_LWdn200%flag,o_LWdn200%name, "LWdn at 200mb", "W/m2")
964	CALL histdef2d(iff,clef_stations(iff),o_LWdn200clr%flag,o_LWdn200clr%name, "LWdn clear sky at 200mb", "W/m2")
965	CALL histdef2d(iff,clef_stations(iff),o_sols%flag,o_sols%name, "Solar rad. at surf.", "W/m2")
966	CALL histdef2d(iff,clef_stations(iff),o_sols0%flag,o_sols0%name, "Solar rad. at surf.", "W/m2")
967	CALL histdef2d(iff,clef_stations(iff),o_soll%flag,o_soll%name, "IR rad. at surface", "W/m2")
968	CALL histdef2d(iff,clef_stations(iff),o_radsol%flag,o_radsol%name, "Rayonnement au sol", "W/m2")
969	CALL histdef2d(iff,clef_stations(iff),o_soll0%flag,o_soll0%name, "IR rad. at surface", "W/m2")
970	CALL histdef2d(iff,clef_stations(iff),o_SWupSFC%flag,o_SWupSFC%name, "SWup at surface", "W/m2")
971	CALL histdef2d(iff,clef_stations(iff),o_SWupSFCclr%flag,o_SWupSFCclr%name, "SWup clear sky at surface", "W/m2")
972	CALL histdef2d(iff,clef_stations(iff),o_SWdnSFC%flag,o_SWdnSFC%name, "SWdn at surface", "W/m2")
973	CALL histdef2d(iff,clef_stations(iff),o_SWdnSFCclr%flag,o_SWdnSFCclr%name, "SWdn clear sky at surface", "W/m2")
974	CALL histdef2d(iff,clef_stations(iff),o_LWupSFC%flag,o_LWupSFC%name, "Upwd. IR rad. at surface", "W/m2")
975	CALL histdef2d(iff,clef_stations(iff),o_LWdnSFC%flag,o_LWdnSFC%name, "Down. IR rad. at surface", "W/m2")
976	CALL histdef2d(iff,clef_stations(iff),o_LWupSFCclr%flag,o_LWupSFCclr%name, "CS Upwd. IR rad. at surface", "W/m2")
977	CALL histdef2d(iff,clef_stations(iff),o_LWdnSFCclr%flag,o_LWdnSFCclr%name, "Down. CS IR rad. at surface", "W/m2")
978	CALL histdef2d(iff,clef_stations(iff),o_bils%flag,o_bils%name, "Surf. total heat flux", "W/m2")
979	CALL histdef2d(iff,clef_stations(iff),o_sens%flag,o_sens%name, "Sensible heat flux", "W/m2")
980	CALL histdef2d(iff,clef_stations(iff),o_fder%flag,o_fder%name, "Heat flux derivation", "W/m2")
981	CALL histdef2d(iff,clef_stations(iff),o_ffonte%flag,o_ffonte%name, "Thermal flux for snow melting", "W/m2")
982	CALL histdef2d(iff,clef_stations(iff),o_fqcalving%flag,o_fqcalving%name, "Ice Calving", "kg/m2/s")
983	CALL histdef2d(iff,clef_stations(iff),o_fqfonte%flag,o_fqfonte%name, "Land ice melt", "kg/m2/s")
984
985	CALL histdef2d(iff,clef_stations(iff),o_taux%flag,o_taux%name, "Zonal wind stress","Pa")
986	CALL histdef2d(iff,clef_stations(iff),o_tauy%flag,o_tauy%name, "Meridional wind stress","Pa")
987
988	DO nsrf = 1, nbsrf
989	CALL histdef2d(iff,clef_stations(iff),o_pourc_srf(nsrf)%flag,o_pourc_srf(nsrf)%name,"% "//clnsurf(nsrf),"%")
990	CALL histdef2d(iff,clef_stations(iff),o_fract_srf(nsrf)%flag,o_fract_srf(nsrf)%name,"Fraction "//clnsurf(nsrf),"1")
991	CALL histdef2d(iff,clef_stations(iff), &
992	o_taux_srf(nsrf)%flag,o_taux_srf(nsrf)%name,"Zonal wind stress"//clnsurf(nsrf),"Pa")
993	CALL histdef2d(iff,clef_stations(iff), &
994	o_tauy_srf(nsrf)%flag,o_tauy_srf(nsrf)%name,"Meridional wind stress "//clnsurf(nsrf),"Pa")
995	CALL histdef2d(iff,clef_stations(iff), &
996	o_tsol_srf(nsrf)%flag,o_tsol_srf(nsrf)%name,"Temperature "//clnsurf(nsrf),"K")
997	CALL histdef2d(iff,clef_stations(iff), &
998	o_u10m_srf(nsrf)%flag,o_u10m_srf(nsrf)%name,"Vent Zonal 10m "//clnsurf(nsrf),"m/s")
999	CALL histdef2d(iff,clef_stations(iff), &
1000	o_evap_srf(nsrf)%flag,o_evap_srf(nsrf)%name,"evaporation at surface "//clnsurf(nsrf),"kg/(s*m2)")
1001	CALL histdef2d(iff,clef_stations(iff), &
1002	o_v10m_srf(nsrf)%flag,o_v10m_srf(nsrf)%name,"Vent meredien 10m "//clnsurf(nsrf),"m/s")
1003	CALL histdef2d(iff,clef_stations(iff), &
1004	o_t2m_srf(nsrf)%flag,o_t2m_srf(nsrf)%name,"Temp 2m "//clnsurf(nsrf),"K")
1005	CALL histdef2d(iff,clef_stations(iff), &
1006	o_sens_srf(nsrf)%flag,o_sens_srf(nsrf)%name,"Sensible heat flux "//clnsurf(nsrf),"W/m2")
1007	CALL histdef2d(iff,clef_stations(iff), &
1008	o_lat_srf(nsrf)%flag,o_lat_srf(nsrf)%name,"Latent heat flux "//clnsurf(nsrf),"W/m2")
1009	CALL histdef2d(iff,clef_stations(iff), &
1010	o_flw_srf(nsrf)%flag,o_flw_srf(nsrf)%name,"LW "//clnsurf(nsrf),"W/m2")
1011	CALL histdef2d(iff,clef_stations(iff), &
1012	o_fsw_srf(nsrf)%flag,o_fsw_srf(nsrf)%name,"SW "//clnsurf(nsrf),"W/m2")
1013	CALL histdef2d(iff,clef_stations(iff), &
1014	o_wbils_srf(nsrf)%flag,o_wbils_srf(nsrf)%name,"Bilan sol "//clnsurf(nsrf),"W/m2" )
1015	CALL histdef2d(iff,clef_stations(iff), &
1016	o_wbilo_srf(nsrf)%flag,o_wbilo_srf(nsrf)%name,"Bilan eau "//clnsurf(nsrf),"kg/(m2*s)")
1017	if (iflag_pbl>1 .and. lev_files(iff).gt.10 ) then
1018	CALL histdef2d(iff,clef_stations(iff), &
1019	o_tke_srf(nsrf)%flag,o_tke_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-")
1020
1021	IF (.NOT.clef_stations(iff)) THEN
1022	!
1023	!IM: there is no way to have one single value in a netcdf file
1024	!
1025	type_ecri(1) = 't_max(X)'
1026	type_ecri(2) = 't_max(X)'
1027	type_ecri(3) = 't_max(X)'
1028	type_ecri(4) = 't_max(X)'
1029	type_ecri(5) = 't_max(X)'
1030	type_ecri(6) = 't_max(X)'
1031	CALL histdef2d(iff,clef_stations(iff), &
1032	o_tke_max_srf(nsrf)%flag,o_tke_max_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-")
1033	type_ecri(:) = type_ecri_files(:)
1034	ENDIF
1035
1036	endif
1037
1038	CALL histdef2d(iff,clef_stations(iff), &
1039	o_albe_srf(nsrf)%flag,o_albe_srf(nsrf)%name,"Albedo VIS surf. "//clnsurf(nsrf),"-")
1040	CALL histdef2d(iff,clef_stations(iff), &
1041	o_rugs_srf(nsrf)%flag,o_rugs_srf(nsrf)%name,"Surface roughness "//clnsurf(nsrf),"m")
1042	CALL histdef2d(iff,clef_stations(iff), &
1043	o_ages_srf(nsrf)%flag,o_ages_srf(nsrf)%name,"Snow age", "day")
1044	END DO
1045
1046	IF (new_aod .AND. (.NOT. aerosol_couple)) THEN
1047	IF (ok_ade.OR.ok_aie) THEN
1048
1049	CALL histdef2d(iff,clef_stations(iff), &
1050	o_od550aer%flag,o_od550aer%name, "Total aerosol optical depth at 550nm", "-")
1051	CALL histdef2d(iff,clef_stations(iff), &
1052	o_od865aer%flag,o_od865aer%name, "Total aerosol optical depth at 870nm", "-")
1053	CALL histdef2d(iff,clef_stations(iff), &
1054	o_absvisaer%flag,o_absvisaer%name, "Absorption aerosol visible optical depth", "-")
1055	CALL histdef2d(iff,clef_stations(iff), &
1056	o_od550lt1aer%flag,o_od550lt1aer%name, "Fine mode optical depth", "-")
1057
1058
1059	CALL histdef2d(iff,clef_stations(iff), &
1060	o_sconcso4%flag,o_sconcso4%name,"Surface Concentration of Sulfate ","kg/m3")
1061	CALL histdef2d(iff,clef_stations(iff), &
1062	o_sconcoa%flag,o_sconcoa%name,"Surface Concentration of Organic Aerosol ","kg/m3")
1063	CALL histdef2d(iff,clef_stations(iff), &
1064	o_sconcbc%flag,o_sconcbc%name,"Surface Concentration of Black Carbon ","kg/m3")
1065	CALL histdef2d(iff,clef_stations(iff), &
1066	o_sconcss%flag,o_sconcss%name,"Surface Concentration of Sea Salt ","kg/m3")
1067	CALL histdef2d(iff,clef_stations(iff), &
1068	o_sconcdust%flag,o_sconcdust%name,"Surface Concentration of Dust ","kg/m3")
1069	CALL histdef3d(iff,clef_stations(iff), &
1070	o_concso4%flag,o_concso4%name,"Concentration of Sulfate ","kg/m3")
1071	CALL histdef3d(iff,clef_stations(iff), &
1072	o_concoa%flag,o_concoa%name,"Concentration of Organic Aerosol ","kg/m3")
1073	CALL histdef3d(iff,clef_stations(iff), &
1074	o_concbc%flag,o_concbc%name,"Concentration of Black Carbon ","kg/m3")
1075	CALL histdef3d(iff,clef_stations(iff), &
1076	o_concss%flag,o_concss%name,"Concentration of Sea Salt ","kg/m3")
1077	CALL histdef3d(iff,clef_stations(iff), &
1078	o_concdust%flag,o_concdust%name,"Concentration of Dust ","kg/m3")
1079	CALL histdef2d(iff,clef_stations(iff), &
1080	o_loadso4%flag,o_loadso4%name,"Column Load of Sulfate ","kg/m2")
1081	CALL histdef2d(iff,clef_stations(iff), &
1082	o_loadoa%flag,o_loadoa%name,"Column Load of Organic Aerosol ","kg/m2")
1083	CALL histdef2d(iff,clef_stations(iff), &
1084	o_loadbc%flag,o_loadbc%name,"Column Load of Black Carbon ","kg/m2")
1085	CALL histdef2d(iff,clef_stations(iff), &
1086	o_loadss%flag,o_loadss%name,"Column Load of Sea Salt ","kg/m2")
1087	CALL histdef2d(iff,clef_stations(iff), &
1088	o_loaddust%flag,o_loaddust%name,"Column Load of Dust ","kg/m2")
1089
1090	DO naero = 1, naero_spc
1091	CALL histdef2d(iff,clef_stations(iff), &
1092	o_tausumaero(naero)%flag,o_tausumaero(naero)%name,"Aerosol Optical depth at 550 nm "//name_aero(naero),"1")
1093	END DO
1094	ENDIF
1095	ENDIF
1096
1097	IF (ok_ade) THEN
1098	CALL histdef2d(iff,clef_stations(iff), &
1099	o_topswad%flag,o_topswad%name, "ADE at TOA", "W/m2")
1100	CALL histdef2d(iff,clef_stations(iff), &
1101	o_solswad%flag,o_solswad%name, "ADE at SRF", "W/m2")
1102
1103	CALL histdef2d(iff,clef_stations(iff), &
1104	o_swtoaas_nat%flag,o_swtoaas_nat%name, "Natural aerosol radiative forcing all-sky at TOA", "W/m2")
1105	CALL histdef2d(iff,clef_stations(iff), &
1106	o_swsrfas_nat%flag,o_swsrfas_nat%name, "Natural aerosol radiative forcing all-sky at SRF", "W/m2")
1107	CALL histdef2d(iff,clef_stations(iff), &
1108	o_swtoacs_nat%flag,o_swtoacs_nat%name, "Natural aerosol radiative forcing clear-sky at TOA", "W/m2")
1109	CALL histdef2d(iff,clef_stations(iff), &
1110	o_swsrfcs_nat%flag,o_swsrfcs_nat%name, "Natural aerosol radiative forcing clear-sky at SRF", "W/m2")
1111
1112	CALL histdef2d(iff,clef_stations(iff), &
1113	o_swtoaas_ant%flag,o_swtoaas_ant%name, "Anthropogenic aerosol radiative forcing all-sky at TOA", "W/m2")
1114	CALL histdef2d(iff,clef_stations(iff), &
1115	o_swsrfas_ant%flag,o_swsrfas_ant%name, "Anthropogenic aerosol radiative forcing all-sky at SRF", "W/m2")
1116	CALL histdef2d(iff,clef_stations(iff), &
1117	o_swtoacs_ant%flag,o_swtoacs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at TOA", "W/m2")
1118	CALL histdef2d(iff,clef_stations(iff), &
1119	o_swsrfcs_ant%flag,o_swsrfcs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at SRF", "W/m2")
1120
1121	IF (.NOT. aerosol_couple) THEN
1122	CALL histdef2d(iff,clef_stations(iff), &
1123	o_swtoacf_nat%flag,o_swtoacf_nat%name, "Natural aerosol impact on cloud radiative forcing at TOA", "W/m2")
1124	CALL histdef2d(iff,clef_stations(iff), &
1125	o_swsrfcf_nat%flag,o_swsrfcf_nat%name, "Natural aerosol impact on cloud radiative forcing at SRF", "W/m2")
1126	CALL histdef2d(iff, clef_stations(iff), o_swtoacf_ant%flag, &
1127	o_swtoacf_ant%name, &
1128	"Anthropogenic aerosol impact on cloud radiative forcing at TOA", &
1129	"W/m2")
1130	CALL histdef2d(iff, clef_stations(iff), o_swsrfcf_ant%flag, &
1131	o_swsrfcf_ant%name, &
1132	"Anthropogenic aerosol impact on cloud radiative forcing at SRF", &
1133	"W/m2")
1134	CALL histdef2d(iff,clef_stations(iff), &
1135	o_swtoacf_zero%flag,o_swtoacf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at TOA", "W/m2")
1136	CALL histdef2d(iff,clef_stations(iff), &
1137	o_swsrfcf_zero%flag,o_swsrfcf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at SRF", "W/m2")
1138	ENDIF
1139	ENDIF
1140
1141	IF (ok_aie) THEN
1142	CALL histdef2d(iff,clef_stations(iff), &
1143	o_topswai%flag,o_topswai%name, "AIE at TOA", "W/m2")
1144	CALL histdef2d(iff,clef_stations(iff), &
1145	o_solswai%flag,o_solswai%name, "AIE at SFR", "W/m2")
1146	!Cloud droplet number concentration
1147	CALL histdef3d(iff,clef_stations(iff), &
1148	o_scdnc%flag,o_scdnc%name, "Cloud droplet number concentration","m-3")
1149	CALL histdef2d(iff,clef_stations(iff), &
1150	o_cldncl%flag,o_cldncl%name, "CDNC at top of liquid water cloud", "m-3")
1151	CALL histdef3d(iff,clef_stations(iff), &
1152	o_reffclws%flag,o_reffclws%name, "Stratiform Cloud Droplet Effective Radius (aerosol diags.)","m")
1153	CALL histdef3d(iff,clef_stations(iff), &
1154	o_reffclwc%flag,o_reffclwc%name, "Convective Cloud Droplet Effective Radius (aerosol diags.)","m")
1155	CALL histdef2d(iff,clef_stations(iff), &
1156	o_cldnvi%flag,o_cldnvi%name, "Column Integrated Cloud Droplet Number", "m-2")
1157	CALL histdef3d(iff,clef_stations(iff), &
1158	o_lcc3d%flag,o_lcc3d%name, "Cloud liquid fraction","1")
1159	CALL histdef3d(iff,clef_stations(iff), &
1160	o_lcc3dcon%flag,o_lcc3dcon%name, "Convective cloud liquid fraction","1")
1161	CALL histdef3d(iff,clef_stations(iff), &
1162	o_lcc3dstra%flag,o_lcc3dstra%name, "Stratiform cloud liquid fraction","1")
1163	CALL histdef2d(iff,clef_stations(iff), &
1164	o_lcc%flag,o_lcc%name, "Cloud liquid fraction at top of cloud","1")
1165	CALL histdef2d(iff,clef_stations(iff), &
1166	o_reffclwtop%flag,o_reffclwtop%name, "Droplet effective radius at top of liquid water cloud", "m")
1167	ENDIF
1168
1169
1170	CALL histdef2d(iff,clef_stations(iff), &
1171	o_alb1%flag,o_alb1%name, "Surface VIS albedo", "-")
1172	CALL histdef2d(iff,clef_stations(iff), &
1173	o_alb2%flag,o_alb2%name, "Surface Near IR albedo", "-")
1174	CALL histdef2d(iff,clef_stations(iff), &
1175	o_cdrm%flag,o_cdrm%name, "Momentum drag coef.", "-")
1176	CALL histdef2d(iff,clef_stations(iff), &
1177	o_cdrh%flag,o_cdrh%name, "Heat drag coef.", "-" )
1178	CALL histdef2d(iff,clef_stations(iff), &
1179	o_cldl%flag,o_cldl%name, "Low-level cloudiness", "-")
1180	CALL histdef2d(iff,clef_stations(iff), &
1181	o_cldm%flag,o_cldm%name, "Mid-level cloudiness", "-")
1182	CALL histdef2d(iff,clef_stations(iff), &
1183	o_cldh%flag,o_cldh%name, "High-level cloudiness", "-")
1184	CALL histdef2d(iff,clef_stations(iff), &
1185	o_cldt%flag,o_cldt%name, "Total cloudiness", "-")
1186	CALL histdef2d(iff,clef_stations(iff), &
1187	o_cldq%flag,o_cldq%name, "Cloud liquid water path", "kg/m2")
1188	CALL histdef2d(iff,clef_stations(iff), &
1189	o_lwp%flag,o_lwp%name, "Cloud water path", "kg/m2")
1190	CALL histdef2d(iff,clef_stations(iff), &
1191	o_iwp%flag,o_iwp%name, "Cloud ice water path", "kg/m2" )
1192	CALL histdef2d(iff,clef_stations(iff), &
1193	o_ue%flag,o_ue%name, "Zonal energy transport", "-")
1194	CALL histdef2d(iff,clef_stations(iff), &
1195	o_ve%flag,o_ve%name, "Merid energy transport", "-")
1196	CALL histdef2d(iff,clef_stations(iff), &
1197	o_uq%flag,o_uq%name, "Zonal humidity transport", "-")
1198	CALL histdef2d(iff,clef_stations(iff), &
1199	o_vq%flag,o_vq%name, "Merid humidity transport", "-")
1200
1201	IF(iflag_con.GE.3) THEN ! sb
1202	CALL histdef2d(iff,clef_stations(iff), &
1203	o_cape%flag,o_cape%name, "Conv avlbl pot ener", "J/kg")
1204	CALL histdef2d(iff,clef_stations(iff), &
1205	o_pbase%flag,o_pbase%name, "Cld base pressure", "Pa")
1206	CALL histdef2d(iff,clef_stations(iff), &
1207	o_ptop%flag,o_ptop%name, "Cld top pressure", "Pa")
1208	CALL histdef2d(iff,clef_stations(iff), &
1209	o_fbase%flag,o_fbase%name, "Cld base mass flux", "kg/m2/s")
1210	if (iflag_con /= 30) then
1211	CALL histdef2d(iff,clef_stations(iff), &
1212	o_plcl%flag,o_plcl%name, "Lifting Condensation Level", "hPa")
1213	CALL histdef2d(iff,clef_stations(iff), &
1214	o_plfc%flag,o_plfc%name, "Level of Free Convection", "hPa")
1215	CALL histdef2d(iff,clef_stations(iff), &
1216	o_wbeff%flag,o_wbeff%name, "Conv. updraft velocity at LFC (<100)", "m/s")
1217	end if
1218	CALL histdef2d(iff,clef_stations(iff), &
1219	o_prw%flag,o_prw%name, "Precipitable water", "kg/m2")
1220	IF (.NOT.clef_stations(iff)) THEN
1221	!
1222	!IM: there is no way to have one single value in a netcdf file
1223	!
1224	type_ecri(1) = 't_max(X)'
1225	type_ecri(2) = 't_max(X)'
1226	type_ecri(3) = 't_max(X)'
1227	type_ecri(4) = 't_max(X)'
1228	type_ecri(5) = 't_max(X)'
1229	type_ecri(6) = 't_max(X)'
1230	CALL histdef2d(iff,clef_stations(iff), &
1231	o_cape_max%flag,o_cape_max%name, "CAPE max.", "J/kg")
1232	ENDIF
1233	type_ecri(:) = type_ecri_files(:)
1234	CALL histdef3d(iff,clef_stations(iff), &
1235	o_upwd%flag,o_upwd%name, "saturated updraft", "kg/m2/s")
1236	CALL histdef3d(iff,clef_stations(iff), &
1237	o_Ma%flag,o_Ma%name, "undilute adiab updraft", "kg/m2/s")
1238	CALL histdef3d(iff,clef_stations(iff), &
1239	o_dnwd%flag,o_dnwd%name, "saturated downdraft", "kg/m2/s")
1240	CALL histdef3d(iff,clef_stations(iff), &
1241	o_dnwd0%flag,o_dnwd0%name, "unsat. downdraft", "kg/m2/s")
1242	CALL histdef3d(iff,clef_stations(iff), &
1243	o_mc%flag,o_mc%name, "Convective mass flux", "kg/m2/s")
1244	type_ecri(1) = 'inst(X)'
1245	type_ecri(2) = 'inst(X)'
1246	type_ecri(3) = 'inst(X)'
1247	type_ecri(4) = 'inst(X)'
1248	type_ecri(5) = 'inst(X)'
1249	type_ecri(6) = 'inst(X)'
1250	CALL histdef2d(iff,clef_stations(iff), &
1251	o_ftime_con%flag,o_ftime_con%name, "Fraction of time convection Occurs", " ")
1252	type_ecri(:) = type_ecri_files(:)
1253	ENDIF !iflag_con .GE. 3
1254
1255	CALL histdef2d(iff,clef_stations(iff), &
1256	o_s_pblh%flag,o_s_pblh%name, "Boundary Layer Height", "m")
1257	CALL histdef2d(iff,clef_stations(iff), &
1258	o_s_pblt%flag,o_s_pblt%name, "t at Boundary Layer Height", "K")
1259	CALL histdef2d(iff,clef_stations(iff), &
1260	o_s_lcl%flag,o_s_lcl%name, "Condensation level", "m")
1261	CALL histdef2d(iff,clef_stations(iff), &
1262	o_s_therm%flag,o_s_therm%name, "Exces du thermique", "K")
1263	!IM : Les champs suivants (s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
1264	!CALL histdef2d(iff,clef_stations(iff), &
1265	!o_s_capCL%flag,o_s_capCL%name, "Conv avlbl pot enerfor ABL", "J/m2" )
1266	!CALL histdef2d(iff,clef_stations(iff), &
1267	!o_s_oliqCL%flag,o_s_oliqCL%name, "Liq Water in BL", "kg/m2")
1268	!CALL histdef2d(iff,clef_stations(iff), &
1269	!o_s_cteiCL%flag,o_s_cteiCL%name, "Instability criteria(ABL)", "K")
1270	!CALL histdef2d(iff,clef_stations(iff), &
1271	!o_s_trmb1%flag,o_s_trmb1%name, "deep_cape(HBTM2)", "J/m2")
1272	!CALL histdef2d(iff,clef_stations(iff), &
1273	!o_s_trmb2%flag,o_s_trmb2%name, "inhibition (HBTM2)", "J/m2")
1274	!CALL histdef2d(iff,clef_stations(iff), &
1275	!o_s_trmb3%flag,o_s_trmb3%name, "Point Omega (HBTM2)", "m")
1276
1277	! Champs interpolles sur des niveaux de pression
1278
1279	type_ecri(1) = 'inst(X)'
1280	type_ecri(2) = 'inst(X)'
1281	type_ecri(3) = 'inst(X)'
1282	type_ecri(4) = 'inst(X)'
1283	type_ecri(5) = 'inst(X)'
1284	type_ecri(6) = 'inst(X)'
1285
1286	! Attention a reverifier
1287
1288	ilev=0
1289	DO k=1, nlevSTD
1290	bb2=clevSTD(k)
1291	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200" &
1292	.OR.bb2.EQ."100".OR.bb2.EQ."50".OR.bb2.EQ."10")THEN
1293	ilev=ilev+1
1294	! print*,'ilev k bb2 flag name ',ilev,k, bb2,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name
1295	CALL histdef2d(iff,clef_stations(iff), &
1296	o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name,"Zonal wind "//bb2//"hPa", "m/s")
1297	CALL histdef2d(iff,clef_stations(iff), &
1298	o_vSTDlevs(ilev)%flag,o_vSTDlevs(ilev)%name,"Meridional wind "//bb2//"hPa", "m/s")
1299	CALL histdef2d(iff,clef_stations(iff), &
1300	o_wSTDlevs(ilev)%flag,o_wSTDlevs(ilev)%name,"Vertical wind "//bb2//"hPa", "Pa/s")
1301	CALL histdef2d(iff,clef_stations(iff), &
1302	o_zSTDlevs(ilev)%flag,o_zSTDlevs(ilev)%name,"Geopotential height "//bb2//"hPa", "m")
1303	CALL histdef2d(iff,clef_stations(iff), &
1304	o_qSTDlevs(ilev)%flag,o_qSTDlevs(ilev)%name,"Specific humidity "//bb2//"hPa", "kg/kg" )
1305	CALL histdef2d(iff,clef_stations(iff), &
1306	o_tSTDlevs(ilev)%flag,o_tSTDlevs(ilev)%name,"Temperature "//bb2//"hPa", "K")
1307	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR."500".OR.bb2.EQ."200".OR.bb2.EQ."50".OR.bb2.EQ."10")
1308	ENDDO
1309	type_ecri(:) = type_ecri_files(:)
1310
1311	CALL histdef2d(iff,clef_stations(iff), &
1312	o_t_oce_sic%flag,o_t_oce_sic%name, "Temp mixte oce-sic", "K")
1313
1314	IF (type_ocean=='slab') &
1315	CALL histdef2d(iff,clef_stations(iff), &
1316	o_slab_bils%flag, o_slab_bils%name,"Bilan au sol sur ocean slab", "W/m2")
1317
1318	! Couplage conv-CL
1319	IF (iflag_con.GE.3) THEN
1320	CALL histdef2d(iff,clef_stations(iff), &
1321	o_ale_bl%flag,o_ale_bl%name, "ALE BL", "m2/s2")
1322	CALL histdef2d(iff,clef_stations(iff), &
1323	o_alp_bl%flag,o_alp_bl%name, "ALP BL", "m2/s2")
1324	ENDIF !(iflag_con.GE.3)
1325
1326	CALL histdef2d(iff,clef_stations(iff), &
1327	o_weakinv%flag,o_weakinv%name, "Weak inversion", "-")
1328	CALL histdef2d(iff,clef_stations(iff), &
1329	o_dthmin%flag,o_dthmin%name, "dTheta mini", "K/m")
1330	CALL histdef2d(iff,clef_stations(iff), &
1331	o_rh2m%flag,o_rh2m%name, "Relative humidity at 2m", "%" )
1332
1333	IF (.NOT.clef_stations(iff)) THEN
1334	!
1335	!IM: there is no way to have one single value in a netcdf file
1336	!
1337	type_ecri(1) = 't_min(X)'
1338	type_ecri(2) = 't_min(X)'
1339	type_ecri(3) = 't_min(X)'
1340	type_ecri(4) = 't_min(X)'
1341	type_ecri(5) = 't_min(X)'
1342	type_ecri(6) = 't_min(X)'
1343	CALL histdef2d(iff,clef_stations(iff),o_rh2m_min%flag,o_rh2m_min%name, "Min Relative humidity at 2m", "%" )
1344	type_ecri(1) = 't_max(X)'
1345	type_ecri(2) = 't_max(X)'
1346	type_ecri(3) = 't_max(X)'
1347	type_ecri(4) = 't_max(X)'
1348	type_ecri(5) = 't_max(X)'
1349	type_ecri(6) = 't_max(X)'
1350	CALL histdef2d(iff,clef_stations(iff),o_rh2m_max%flag,o_rh2m_max%name, "Max Relative humidity at 2m", "%" )
1351	ENDIF
1352
1353	type_ecri(:) = type_ecri_files(:)
1354	CALL histdef2d(iff,clef_stations(iff),o_qsat2m%flag,o_qsat2m%name, "Saturant humidity at 2m", "%")
1355	CALL histdef2d(iff,clef_stations(iff),o_tpot%flag,o_tpot%name, "Surface air potential temperature", "K")
1356	CALL histdef2d(iff,clef_stations(iff), &
1357	o_tpote%flag,o_tpote%name, "Surface air equivalent potential temperature", "K")
1358	CALL histdef2d(iff,clef_stations(iff),o_SWnetOR%flag,o_SWnetOR%name, "Sfce net SW radiation OR", "W/m2")
1359	CALL histdef2d(iff,clef_stations(iff),o_SWdownOR%flag,o_SWdownOR%name, "Sfce incident SW radiation OR", "W/m2")
1360	CALL histdef2d(iff,clef_stations(iff),o_LWdownOR%flag,o_LWdownOR%name, "Sfce incident LW radiation OR", "W/m2")
1361	CALL histdef2d(iff,clef_stations(iff),o_snowl%flag,o_snowl%name, "Solid Large-scale Precip.", "kg/(m2*s)")
1362
1363	CALL histdef2d(iff,clef_stations(iff),o_solldown%flag,o_solldown%name, "Down. IR rad. at surface", "W/m2")
1364	CALL histdef2d(iff,clef_stations(iff),o_dtsvdfo%flag,o_dtsvdfo%name, "Boundary-layer dTs(o)", "K/s")
1365	CALL histdef2d(iff,clef_stations(iff),o_dtsvdft%flag,o_dtsvdft%name, "Boundary-layer dTs(t)", "K/s")
1366	CALL histdef2d(iff,clef_stations(iff),o_dtsvdfg%flag,o_dtsvdfg%name, "Boundary-layer dTs(g)", "K/s")
1367	CALL histdef2d(iff,clef_stations(iff),o_dtsvdfi%flag,o_dtsvdfi%name, "Boundary-layer dTs(g)", "K/s")
1368	CALL histdef2d(iff,clef_stations(iff),o_rugs%flag,o_rugs%name, "rugosity", "-" )
1369
1370	! Champs 3D:
1371	CALL histdef3d(iff,clef_stations(iff),o_ec550aer%flag,o_ec550aer%name, "Extinction at 550nm", "m^-1")
1372	CALL histdef3d(iff,clef_stations(iff),o_lwcon%flag,o_lwcon%name, "Cloud liquid water content", "kg/kg")
1373	CALL histdef3d(iff,clef_stations(iff),o_iwcon%flag,o_iwcon%name, "Cloud ice water content", "kg/kg")
1374	CALL histdef3d(iff,clef_stations(iff),o_temp%flag,o_temp%name, "Air temperature", "K" )
1375	CALL histdef3d(iff,clef_stations(iff),o_theta%flag,o_theta%name, "Potential air temperature", "K" )
1376	CALL histdef3d(iff,clef_stations(iff),o_ovap%flag,o_ovap%name, "Specific humidity", "kg/kg" )
1377	CALL histdef3d(iff,clef_stations(iff), &
1378	o_ovapinit%flag,o_ovapinit%name, "Specific humidity (begin of timestep)", "kg/kg" )
1379	CALL histdef3d(iff,clef_stations(iff), &
1380	o_geop%flag,o_geop%name, "Geopotential height", "m2/s2")
1381	CALL histdef3d(iff,clef_stations(iff), &
1382	o_vitu%flag,o_vitu%name, "Zonal wind", "m/s" )
1383	CALL histdef3d(iff,clef_stations(iff), &
1384	o_vitv%flag,o_vitv%name, "Meridional wind", "m/s" )
1385	CALL histdef3d(iff,clef_stations(iff), &
1386	o_vitw%flag,o_vitw%name, "Vertical wind", "Pa/s" )
1387	CALL histdef3d(iff,clef_stations(iff), &
1388	o_pres%flag,o_pres%name, "Air pressure", "Pa" )
1389	CALL histdef3d(iff,clef_stations(iff), &
1390	o_paprs%flag,o_paprs%name, "Air pressure Inter-Couches", "Pa" )
1391	CALL histdef3d(iff,clef_stations(iff), &
1392	o_mass%flag,o_mass%name, "Masse Couches", "kg/m2" )
1393	CALL histdef3d(iff,clef_stations(iff), &
1394	o_zfull%flag,o_zfull%name, "Altitude of full pressure levels", "m" )
1395	CALL histdef3d(iff,clef_stations(iff), &
1396	o_zhalf%flag,o_zhalf%name, "Altitude of half pressure levels", "m" )
1397	CALL histdef3d(iff,clef_stations(iff), &
1398	o_rneb%flag,o_rneb%name, "Cloud fraction", "-")
1399	CALL histdef3d(iff,clef_stations(iff), &
1400	o_rnebcon%flag,o_rnebcon%name, "Convective Cloud Fraction", "-")
1401	CALL histdef3d(iff,clef_stations(iff), &
1402	o_rhum%flag,o_rhum%name, "Relative humidity", "-")
1403	CALL histdef3d(iff,clef_stations(iff), &
1404	o_ozone%flag,o_ozone%name, "Ozone mole fraction", "-")
1405	if (read_climoz == 2) &
1406	CALL histdef3d(iff,clef_stations(iff), &
1407	o_ozone_light%flag,o_ozone_light%name, &
1408	"Daylight ozone mole fraction", "-")
1409	CALL histdef3d(iff,clef_stations(iff), &
1410	o_dtphy%flag,o_dtphy%name, "Physics dT", "K/s")
1411	CALL histdef3d(iff,clef_stations(iff), &
1412	o_dqphy%flag,o_dqphy%name, "Physics dQ", "(kg/kg)/s")
1413	CALL histdef3d(iff,clef_stations(iff), &
1414	o_cldtau%flag,o_cldtau%name, "Cloud optical thickness", "1")
1415	CALL histdef3d(iff,clef_stations(iff), &
1416	o_cldemi%flag,o_cldemi%name, "Cloud optical emissivity", "1")
1417	!IM: bug ?? dimensionnement variables (klon,klev+1) pmflxr, pmflxs, prfl, psfl
1418	CALL histdef3d(iff,clef_stations(iff), &
1419	o_pr_con_l%flag,o_pr_con_l%name, "Convective precipitation lic", " ")
1420	CALL histdef3d(iff,clef_stations(iff), &
1421	o_pr_con_i%flag,o_pr_con_i%name, "Convective precipitation ice", " ")
1422	CALL histdef3d(iff,clef_stations(iff), &
1423	o_pr_lsc_l%flag,o_pr_lsc_l%name, "Large scale precipitation lic", " ")
1424	CALL histdef3d(iff,clef_stations(iff), &
1425	o_pr_lsc_i%flag,o_pr_lsc_i%name, "Large scale precipitation ice", " ")
1426	!Cloud droplet effective radius
1427	CALL histdef3d(iff,clef_stations(iff), &
1428	o_re%flag,o_re%name, "Cloud droplet effective radius","um")
1429	CALL histdef3d(iff,clef_stations(iff), &
1430	o_fl%flag,o_fl%name, "Denominator of Cloud droplet effective radius"," ")
1431	!FH Sorties pour la couche limite
1432	if (iflag_pbl>1) then
1433	CALL histdef3d(iff,clef_stations(iff), &
1434	o_tke%flag,o_tke%name, "TKE", "m2/s2")
1435	IF (.NOT.clef_stations(iff)) THEN
1436	!
1437	!IM: there is no way to have one single value in a netcdf file
1438	!
1439	type_ecri(1) = 't_max(X)'
1440	type_ecri(2) = 't_max(X)'
1441	type_ecri(3) = 't_max(X)'
1442	type_ecri(4) = 't_max(X)'
1443	type_ecri(5) = 't_max(X)'
1444	type_ecri(6) = 't_max(X)'
1445	CALL histdef3d(iff,clef_stations(iff), &
1446	o_tke_max%flag,o_tke_max%name, "TKE max", "m2/s2")
1447	ENDIF
1448	type_ecri(:) = type_ecri_files(:)
1449	endif
1450
1451	CALL histdef3d(iff,clef_stations(iff), &
1452	o_kz%flag,o_kz%name, "Kz melange", "m2/s")
1453	IF (.NOT.clef_stations(iff)) THEN
1454	!
1455	!IM: there is no way to have one single value in a netcdf file
1456	!
1457	type_ecri(1) = 't_max(X)'
1458	type_ecri(2) = 't_max(X)'
1459	type_ecri(3) = 't_max(X)'
1460	type_ecri(4) = 't_max(X)'
1461	type_ecri(5) = 't_max(X)'
1462	type_ecri(6) = 't_max(X)'
1463	CALL histdef3d(iff,clef_stations(iff), &
1464	o_kz_max%flag,o_kz_max%name, "Kz melange max", "m2/s" )
1465	ENDIF
1466	type_ecri(:) = type_ecri_files(:)
1467	CALL histdef3d(iff,clef_stations(iff), &
1468	o_clwcon%flag,o_clwcon%name, "Convective Cloud Liquid water content", "kg/kg")
1469	CALL histdef3d(iff,clef_stations(iff), &
1470	o_dtdyn%flag,o_dtdyn%name, "Dynamics dT", "K/s")
1471	CALL histdef3d(iff,clef_stations(iff), &
1472	o_dqdyn%flag,o_dqdyn%name, "Dynamics dQ", "(kg/kg)/s")
1473	CALL histdef3d(iff,clef_stations(iff), &
1474	o_dudyn%flag,o_dudyn%name, "Dynamics dU", "m/s2")
1475	CALL histdef3d(iff,clef_stations(iff), &
1476	o_dvdyn%flag,o_dvdyn%name, "Dynamics dV", "m/s2")
1477	CALL histdef3d(iff,clef_stations(iff), &
1478	o_dtcon%flag,o_dtcon%name, "Convection dT", "K/s")
1479	CALL histdef3d(iff,clef_stations(iff), &
1480	o_ducon%flag,o_ducon%name, "Convection du", "m/s2")
1481	CALL histdef3d(iff,clef_stations(iff), &
1482	o_dqcon%flag,o_dqcon%name, "Convection dQ", "(kg/kg)/s")
1483
1484	! Wakes
1485	IF(iflag_con.EQ.3) THEN
1486	IF (iflag_wake >= 1) THEN
1487	CALL histdef2d(iff,clef_stations(iff), &
1488	o_ale_wk%flag,o_ale_wk%name, "ALE WK", "m2/s2")
1489	CALL histdef2d(iff,clef_stations(iff), &
1490	o_alp_wk%flag,o_alp_wk%name, "ALP WK", "m2/s2")
1491	CALL histdef2d(iff,clef_stations(iff),o_wake_h%flag,o_wake_h%name, "wake_h", "-")
1492	CALL histdef2d(iff,clef_stations(iff),o_wake_s%flag,o_wake_s%name, "wake_s", "-")
1493	CALL histdef3d(iff,clef_stations(iff),o_dtwak%flag,o_dtwak%name, "Wake dT", "K/s")
1494	CALL histdef3d(iff,clef_stations(iff),o_dqwak%flag,o_dqwak%name, "Wake dQ", "(kg/kg)/s")
1495	CALL histdef3d(iff,clef_stations(iff),o_wake_deltat%flag,o_wake_deltat%name, "wake_deltat", " ")
1496	CALL histdef3d(iff,clef_stations(iff),o_wake_deltaq%flag,o_wake_deltaq%name, "wake_deltaq", " ")
1497	CALL histdef3d(iff,clef_stations(iff),o_wake_omg%flag,o_wake_omg%name, "wake_omg", " ")
1498	CALL histdef2d(iff,clef_stations(iff),o_wape%flag,o_WAPE%name, "WAPE", "m2/s2")
1499	ENDIF
1500	CALL histdef2d(iff,clef_stations(iff), &
1501	o_ale%flag,o_ale%name, "ALE", "m2/s2")
1502	CALL histdef2d(iff,clef_stations(iff), &
1503	o_alp%flag,o_alp%name, "ALP", "W/m2")
1504	CALL histdef2d(iff,clef_stations(iff),o_cin%flag,o_cin%name, "Convective INhibition", "m2/s2")
1505	CALL histdef3d(iff,clef_stations(iff),o_Vprecip%flag,o_Vprecip%name, "precipitation vertical profile", "-")
1506	CALL histdef3d(iff,clef_stations(iff),o_ftd%flag,o_ftd%name, "tend temp due aux descentes precip", "-")
1507	CALL histdef3d(iff,clef_stations(iff),o_fqd%flag,o_fqd%name,"tend vap eau due aux descentes precip", "-")
1508	ENDIF !(iflag_con.EQ.3)
1509
1510	CALL histdef3d(iff,clef_stations(iff),o_dtlsc%flag,o_dtlsc%name, "Condensation dT", "K/s")
1511	CALL histdef3d(iff,clef_stations(iff),o_dtlschr%flag,o_dtlschr%name,"Large-scale condensational heating rate","K/s")
1512	CALL histdef3d(iff,clef_stations(iff),o_dqlsc%flag,o_dqlsc%name, "Condensation dQ", "(kg/kg)/s")
1513	CALL histdef3d(iff,clef_stations(iff),o_dtvdf%flag,o_dtvdf%name, "Boundary-layer dT", "K/s")
1514	CALL histdef3d(iff,clef_stations(iff),o_dqvdf%flag,o_dqvdf%name, "Boundary-layer dQ", "(kg/kg)/s")
1515	CALL histdef3d(iff,clef_stations(iff),o_dteva%flag,o_dteva%name, "Reevaporation dT", "K/s")
1516	CALL histdef3d(iff,clef_stations(iff),o_dqeva%flag,o_dqeva%name, "Reevaporation dQ", "(kg/kg)/s")
1517	CALL histdef3d(iff,clef_stations(iff),o_ptconv%flag,o_ptconv%name, "POINTS CONVECTIFS", " ")
1518	CALL histdef3d(iff,clef_stations(iff),o_ratqs%flag,o_ratqs%name, "RATQS", " ")
1519	CALL histdef3d(iff,clef_stations(iff),o_dtthe%flag,o_dtthe%name, "Thermal dT", "K/s")
1520
1521	if(iflag_thermals.gt.1) THEN
1522	CALL histdef3d(iff,clef_stations(iff),o_dqlscth%flag,o_dqlscth%name, "dQ therm.", "(kg/kg)/s")
1523	CALL histdef3d(iff,clef_stations(iff),o_dqlscst%flag,o_dqlscst%name, "dQ strat.", "(kg/kg)/s")
1524	CALL histdef3d(iff,clef_stations(iff),o_dtlscth%flag,o_dtlscth%name, "dQ therm.", "K/s")
1525	CALL histdef3d(iff,clef_stations(iff),o_dtlscst%flag,o_dtlscst%name, "dQ strat.", "K/s")
1526	CALL histdef2d(iff,clef_stations(iff),o_plulth%flag,o_plulth%name, "Rainfall therm.", "K/s")
1527	CALL histdef2d(iff,clef_stations(iff),o_plulst%flag,o_plulst%name, "Rainfall strat.", "K/s")
1528	CALL histdef2d(iff,clef_stations(iff),o_lmaxth%flag,o_lmaxth%name, "Upper level thermals", "")
1529	CALL histdef3d(iff,clef_stations(iff),o_ptconvth%flag,o_ptconvth%name, "POINTS CONVECTIFS therm.", " ")
1530	CALL histdef3d(iff,clef_stations(iff),o_f_th%flag,o_f_th%name, "Thermal plume mass flux", "kg/(m2*s)")
1531	CALL histdef3d(iff,clef_stations(iff),o_e_th%flag,o_e_th%name,"Thermal plume entrainment","K/s")
1532	CALL histdef3d(iff,clef_stations(iff),o_w_th%flag,o_w_th%name,"Thermal plume vertical velocity","m/s")
1533	CALL histdef3d(iff,clef_stations(iff), &
1534	o_lambda_th%flag,o_lambda_th%name,"Thermal plume vertical velocity","m/s")
1535	CALL histdef2d(iff,clef_stations(iff), &
1536	o_ftime_th%flag,o_ftime_th%name,"Fraction of time Shallow convection occurs"," ")
1537	CALL histdef3d(iff,clef_stations(iff), &
1538	o_q_th%flag,o_q_th%name, "Thermal plume total humidity", "kg/kg")
1539	CALL histdef3d(iff,clef_stations(iff), &
1540	o_a_th%flag,o_a_th%name, "Thermal plume fraction", "")
1541	CALL histdef3d(iff,clef_stations(iff), &
1542	o_d_th%flag,o_d_th%name, "Thermal plume detrainment", "K/s")
1543
1544	CALL histdef2d(iff,clef_stations(iff), &
1545	o_f0_th%flag,o_f0_th%name, "Thermal closure mass flux", "K/s")
1546	CALL histdef2d(iff,clef_stations(iff), &
1547	o_zmax_th%flag,o_zmax_th%name, "Thermal plume height", "K/s")
1548	CALL histdef3d(iff,clef_stations(iff), &
1549	o_dqthe%flag,o_dqthe%name, "Thermal dQ", "(kg/kg)/s")
1550	endif !iflag_thermals.gt.1
1551	CALL histdef3d(iff,clef_stations(iff), &
1552	o_dtajs%flag,o_dtajs%name, "Dry adjust. dT", "K/s")
1553	CALL histdef3d(iff,clef_stations(iff), &
1554	o_dqajs%flag,o_dqajs%name, "Dry adjust. dQ", "(kg/kg)/s")
1555	CALL histdef3d(iff,clef_stations(iff), &
1556	o_dtswr%flag,o_dtswr%name, "SW radiation dT", "K/s")
1557	CALL histdef3d(iff,clef_stations(iff), &
1558	o_dtsw0%flag,o_dtsw0%name, "CS SW radiation dT", "K/s")
1559	CALL histdef3d(iff,clef_stations(iff), &
1560	o_dtlwr%flag,o_dtlwr%name, "LW radiation dT", "K/s")
1561	CALL histdef3d(iff,clef_stations(iff), &
1562	o_dtlw0%flag,o_dtlw0%name, "CS LW radiation dT", "K/s")
1563	CALL histdef3d(iff,clef_stations(iff), &
1564	o_dtec%flag,o_dtec%name, "Cinetic dissip dT", "K/s")
1565	CALL histdef3d(iff,clef_stations(iff), &
1566	o_duvdf%flag,o_duvdf%name, "Boundary-layer dU", "m/s2")
1567	CALL histdef3d(iff,clef_stations(iff), &
1568	o_dvvdf%flag,o_dvvdf%name, "Boundary-layer dV", "m/s2")
1569
1570	IF (ok_orodr) THEN
1571	CALL histdef3d(iff,clef_stations(iff), &
1572	o_duoro%flag,o_duoro%name, "Orography dU", "m/s2")
1573	CALL histdef3d(iff,clef_stations(iff), &
1574	o_dvoro%flag,o_dvoro%name, "Orography dV", "m/s2")
1575	CALL histdef3d(iff,clef_stations(iff), &
1576	o_dtoro%flag,o_dtoro%name, "Orography dT", "K/s")
1577	ENDIF
1578
1579	IF (ok_orolf) THEN
1580	CALL histdef3d(iff,clef_stations(iff), &
1581	o_dulif%flag,o_dulif%name, "Orography dU", "m/s2")
1582	CALL histdef3d(iff,clef_stations(iff), &
1583	o_dvlif%flag,o_dvlif%name, "Orography dV", "m/s2")
1584	CALL histdef3d(iff,clef_stations(iff), &
1585	o_dtlif%flag,o_dtlif%name, "Orography dT", "K/s")
1586	ENDIF
1587
1588	IF (ok_hines) then
1589	CALL histdef3d(iff,clef_stations(iff), &
1590	o_duhin%flag,o_duhin%name, "Hines GWD dU", "m/s2")
1591	CALL histdef3d(iff,clef_stations(iff), &
1592	o_dvhin%flag,o_dvhin%name, "Hines GWD dV", "m/s2")
1593
1594	CALL histdef3d(iff,clef_stations(iff), &
1595	o_dthin%flag,o_dthin%name, "Hines GWD dT", "K/s")
1596	ENDIF
1597
1598	CALL histdef3d(iff,clef_stations(iff), &
1599	o_rsu%flag,o_rsu%name, "SW upward radiation", "W m-2")
1600	CALL histdef3d(iff,clef_stations(iff), &
1601	o_rsd%flag,o_rsd%name, "SW downward radiation", "W m-2")
1602	CALL histdef3d(iff,clef_stations(iff), &
1603	o_rlu%flag,o_rlu%name, "LW upward radiation", "W m-2")
1604	CALL histdef3d(iff,clef_stations(iff), &
1605	o_rld%flag,o_rld%name, "LW downward radiation", "W m-2")
1606
1607	CALL histdef3d(iff,clef_stations(iff), &
1608	o_rsucs%flag,o_rsucs%name, "SW CS upward radiation", "W m-2")
1609	CALL histdef3d(iff,clef_stations(iff), &
1610	o_rsdcs%flag,o_rsdcs%name, "SW CS downward radiation", "W m-2")
1611	CALL histdef3d(iff,clef_stations(iff), &
1612	o_rlucs%flag,o_rlucs%name, "LW CS upward radiation", "W m-2")
1613	CALL histdef3d(iff,clef_stations(iff), &
1614	o_rldcs%flag,o_rldcs%name, "LW CS downward radiation", "W m-2")
1615
1616	CALL histdef3d(iff,clef_stations(iff), &
1617	o_tnt%flag,o_tnt%name, "Tendency of air temperature", "K s-1")
1618
1619	CALL histdef3d(iff,clef_stations(iff), &
1620	o_tntc%flag,o_tntc%name, "Tendency of air temperature due to Moist Convection", &
1621	"K s-1")
1622
1623	CALL histdef3d(iff,clef_stations(iff), &
1624	o_tntr%flag,o_tntr%name, "Air temperature tendency due to Radiative heating", &
1625	"K s-1")
1626
1627	CALL histdef3d(iff,clef_stations(iff), &
1628	o_tntscpbl%flag,o_tntscpbl%name, "Air temperature tendency due to St cloud and precipitation and BL mixing", &
1629	"K s-1")
1630
1631	CALL histdef3d(iff,clef_stations(iff), &
1632	o_tnhus%flag,o_tnhus%name, "Tendency of specific humidity", "s-1")
1633
1634	CALL histdef3d(iff,clef_stations(iff), &
1635	o_tnhusc%flag,o_tnhusc%name, "Tendency of specific humidity due to convection", "s-1")
1636
1637	CALL histdef3d(iff,clef_stations(iff), &
1638	o_tnhusscpbl%flag,o_tnhusscpbl%name, "Tendency of Specific humidity due to ST cl, precip and BL mixing", &
1639	"s-1")
1640
1641	CALL histdef3d(iff,clef_stations(iff), &
1642	o_evu%flag,o_evu%name, "Eddy viscosity coefficient for Momentum Variables", "m2 s-1")
1643
1644	CALL histdef3d(iff,clef_stations(iff), &
1645	o_h2o%flag,o_h2o%name, "Mass Fraction of Water", "1")
1646
1647	CALL histdef3d(iff,clef_stations(iff), &
1648	o_mcd%flag,o_mcd%name, "Downdraft COnvective Mass Flux", "kg/(m2*s)")
1649
1650	CALL histdef3d(iff,clef_stations(iff), &
1651	o_dmc%flag,o_dmc%name, "Deep COnvective Mass Flux", "kg/(m2*s)")
1652
1653	CALL histdef3d(iff,clef_stations(iff), &
1654	o_ref_liq%flag,o_ref_liq%name, "Effective radius of convective cloud liquid water particle", "m")
1655
1656	CALL histdef3d(iff,clef_stations(iff), &
1657	o_ref_ice%flag,o_ref_ice%name, "Effective radius of startiform cloud ice particle", "m")
1658
1659	if (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR. &
1660	RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR. &
1661	RCFC12_per.NE.RCFC12_act) THEN
1662
1663	CALL histdef2d(iff,clef_stations(iff),o_rsut4co2%flag,o_rsut4co2%name, &
1664	"TOA Out SW in 4xCO2 atmosphere", "W/m2")
1665	CALL histdef2d(iff,clef_stations(iff),o_rlut4co2%flag,o_rlut4co2%name, &
1666	"TOA Out LW in 4xCO2 atmosphere", "W/m2")
1667	CALL histdef2d(iff,clef_stations(iff),o_rsutcs4co2%flag,o_rsutcs4co2%name, &
1668	"TOA Out CS SW in 4xCO2 atmosphere", "W/m2")
1669	CALL histdef2d(iff,clef_stations(iff),o_rlutcs4co2%flag,o_rlutcs4co2%name, &
1670	"TOA Out CS LW in 4xCO2 atmosphere", "W/m2")
1671
1672	CALL histdef3d(iff,clef_stations(iff),o_rsu4co2%flag,o_rsu4co2%name, &
1673	"Upwelling SW 4xCO2 atmosphere", "W/m2")
1674	CALL histdef3d(iff,clef_stations(iff),o_rlu4co2%flag,o_rlu4co2%name, &
1675	"Upwelling LW 4xCO2 atmosphere", "W/m2")
1676	CALL histdef3d(iff,clef_stations(iff),o_rsucs4co2%flag,o_rsucs4co2%name, &
1677	"Upwelling CS SW 4xCO2 atmosphere", "W/m2")
1678	CALL histdef3d(iff,clef_stations(iff),o_rlucs4co2%flag,o_rlucs4co2%name, &
1679	"Upwelling CS LW 4xCO2 atmosphere", "W/m2")
1680
1681	CALL histdef3d(iff,clef_stations(iff),o_rsd4co2%flag,o_rsd4co2%name, &
1682	"Downwelling SW 4xCO2 atmosphere", "W/m2")
1683	CALL histdef3d(iff,clef_stations(iff),o_rld4co2%flag,o_rld4co2%name, &
1684	"Downwelling LW 4xCO2 atmosphere", "W/m2")
1685	CALL histdef3d(iff,clef_stations(iff),o_rsdcs4co2%flag,o_rsdcs4co2%name, &
1686	"Downwelling CS SW 4xCO2 atmosphere", "W/m2")
1687	CALL histdef3d(iff,clef_stations(iff),o_rldcs4co2%flag,o_rldcs4co2%name, &
1688	"Downwelling CS LW 4xCO2 atmosphere", "W/m2")
1689
1690	endif
1691
1692
1693	IF (nqtot>=3) THEN
1694	DO iq=3,nqtot
1695	iiq=niadv(iq)
1696	o_trac(iq-2) = ctrl_out((/ 4, 5, 1, 1, 1, 10 /),tname(iiq))
1697	CALL histdef3d (iff,clef_stations(iff), &
1698	o_trac(iq-2)%flag,o_trac(iq-2)%name,'Tracer '//ttext(iiq), "-" )
1699	ENDDO
1700	ENDIF
1701
1702	CALL histend(nid_files(iff))
1703
1704	ndex2d = 0
1705	ndex3d = 0
1706
1707	ENDIF ! clef_files
1708
1709	ENDDO ! iff
1710	write(lunout,*)'Fin phys_output_mod.F90'
1711	end subroutine phys_output_open
1712
1713	SUBROUTINE histdef2d (iff,lpoint,flag_var,nomvar,titrevar,unitvar)
1714
1715	use ioipsl
1716	USE dimphy
1717	USE mod_phys_lmdz_para
1718	USE iophy
1719
1720	IMPLICIT NONE
1721
1722	include "dimensions.h"
1723	include "temps.h"
1724	include "indicesol.h"
1725	include "clesphys.h"
1726
1727	integer :: iff
1728	logical :: lpoint
1729	integer, dimension(nfiles) :: flag_var
1730	character(len=20) :: nomvar
1731	character(len=*) :: titrevar
1732	character(len=*) :: unitvar
1733
1734	real zstophym
1735
1736	if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then
1737	zstophym=zoutm(iff)
1738	else
1739	zstophym=zdtime
1740	endif
1741
1742	! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def
1743	call conf_physoutputs(nomvar,flag_var)
1744
1745	if(.NOT.lpoint) THEN
1746	if ( flag_var(iff)<=lev_files(iff) ) then
1747	call histdef (nid_files(iff),nomvar,titrevar,unitvar, &
1748	iim,jj_nb,nhorim(iff), 1,1,1, -99, 32, &
1749	type_ecri(iff), zstophym,zoutm(iff))
1750	endif
1751	else
1752	if ( flag_var(iff)<=lev_files(iff) ) then
1753	call histdef (nid_files(iff),nomvar,titrevar,unitvar, &
1754	npstn,1,nhorim(iff), 1,1,1, -99, 32, &
1755	type_ecri(iff), zstophym,zoutm(iff))
1756	endif
1757	endif
1758	end subroutine histdef2d
1759
1760	SUBROUTINE histdef3d (iff,lpoint,flag_var,nomvar,titrevar,unitvar)
1761
1762	use ioipsl
1763	USE dimphy
1764	USE mod_phys_lmdz_para
1765	USE iophy
1766
1767	IMPLICIT NONE
1768
1769	include "dimensions.h"
1770	include "temps.h"
1771	include "indicesol.h"
1772	include "clesphys.h"
1773
1774	integer :: iff
1775	logical :: lpoint
1776	integer, dimension(nfiles) :: flag_var
1777	character(len=20) :: nomvar
1778	character(len=*) :: titrevar
1779	character(len=*) :: unitvar
1780
1781	real zstophym
1782
1783	! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def
1784	call conf_physoutputs(nomvar,flag_var)
1785
1786	if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then
1787	zstophym=zoutm(iff)
1788	else
1789	zstophym=zdtime
1790	endif
1791
1792	if(.NOT.lpoint) THEN
1793	if ( flag_var(iff)<=lev_files(iff) ) then
1794	call histdef (nid_files(iff), nomvar, titrevar, unitvar, &
1795	iim, jj_nb, nhorim(iff), klev, levmin(iff), &
1796	levmax(iff)-levmin(iff)+1, nvertm(iff), 32, type_ecri(iff), &
1797	zstophym, zoutm(iff))
1798	endif
1799	else
1800	if ( flag_var(iff)<=lev_files(iff) ) then
1801	call histdef (nid_files(iff), nomvar, titrevar, unitvar, &
1802	npstn,1,nhorim(iff), klev, levmin(iff), &
1803	levmax(iff)-levmin(iff)+1, nvertm(iff), 32, &
1804	type_ecri(iff), zstophym,zoutm(iff))
1805	endif
1806	endif
1807	end subroutine histdef3d
1808
1809	SUBROUTINE conf_physoutputs(nam_var,flag_var)
1810	!!! Lecture des noms et niveau de sortie des variables dans output.def
1811	! en utilisant les routines getin de IOIPSL
1812	use ioipsl
1813
1814	IMPLICIT NONE
1815
1816	include 'iniprint.h'
1817
1818	character(len=20) :: nam_var
1819	integer, dimension(nfiles) :: flag_var
1820
1821	IF(prt_level>10) WRITE(lunout,*)'Avant getin: nam_var flag_var ',nam_var,flag_var(:)
1822	call getin('flag_'//nam_var,flag_var)
1823	call getin('name_'//nam_var,nam_var)
1824	IF(prt_level>10) WRITE(lunout,*)'Apres getin: nam_var flag_var ',nam_var,flag_var(:)
1825
1826	END SUBROUTINE conf_physoutputs
1827
1828	SUBROUTINE convers_timesteps(str,dtime,timestep)
1829
1830	use ioipsl
1831	USE phys_cal_mod
1832
1833	IMPLICIT NONE
1834
1835	character(len=20) :: str
1836	character(len=10) :: type
1837	integer :: ipos,il
1838	real :: ttt,xxx,timestep,dayseconde,dtime
1839	parameter (dayseconde=86400.)
1840	include "temps.h"
1841	include "comconst.h"
1842	include "iniprint.h"
1843
1844	ipos=scan(str,'0123456789.',.true.)
1845	!
1846	il=len_trim(str)
1847	write(lunout,*)ipos,il
1848	read(str(1:ipos),*) ttt
1849	write(lunout,*)ttt
1850	type=str(ipos+1:il)
1851
1852
1853	if ( il == ipos ) then
1854	type='day'
1855	endif
1856
1857	if ( type == 'day'.or.type == 'days'.or.type == 'jours'.or.type == 'jour' ) timestep = ttt * dayseconde
1858	if ( type == 'mounths'.or.type == 'mth'.or.type == 'mois' ) then
1859	write(lunout,*)'annee_ref,day_ref mon_len',annee_ref,day_ref,mth_len
1860	timestep = ttt * dayseconde * mth_len
1861	endif
1862	if ( type == 'hours'.or.type == 'hr'.or.type == 'heurs') timestep = ttt * dayseconde / 24.
1863	if ( type == 'mn'.or.type == 'minutes' ) timestep = ttt * 60.
1864	if ( type == 's'.or.type == 'sec'.or.type == 'secondes' ) timestep = ttt
1865	if ( type == 'TS' ) timestep = ttt * dtime
1866
1867	write(lunout,*)'type = ',type
1868	write(lunout,*)'nb j/h/m = ',ttt
1869	write(lunout,*)'timestep(s)=',timestep
1870
1871	END SUBROUTINE convers_timesteps
1872
1873	END MODULE phys_output_mod
1874

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

Download in other formats:

Original Format