Context Navigation

phys_output_mod.F90 @ 1724

Last change on this file since 1724 was 1724, checked in by musat, 11 years ago

Ajout coefficient beta (=Evap/Evap? potentielle) fonction de la latitude
et de la saison (voir these AbdulKhadre? Traore).

Frederic Hourdin
IM

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

File size: 112.6 KB

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

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

Download in other formats:

Original Format