Context Navigation

phys_output_mod.F90 @ 1707

Last change on this file since 1707 was 1707, checked in by Laurent Fairhead, 11 years ago

Version testing basée sur la r1706

Testing release based on r1706

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

File size: 112.1 KB

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

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

Download in other formats:

Original Format