Context Navigation

phys_output_mod.F90 @ 925

Last change on this file since 925 was 922, checked in by Laurent Fairhead, 17 years ago
Pour cvs sorry LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 77.9 KB

Line
1	!
2	! $Header$
3	!
4	! Abderrahmane 12 2007
5	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
6	!!! Ecreture des Sorties du modele dans les fichiers Netcdf :
7	! histmth.nc : moyennes mensuelles
8	! histday.nc : moyennes journalieres
9	! histhf.nc : moyennes toutes les 3 heures
10	! histins.nc : valeurs instantanees
11	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
12
13	MODULE phys_output_mod
14
15	IMPLICIT NONE
16
17	integer, parameter :: nfiles = 4
18	logical, dimension(nfiles), save :: clef_files
19	integer, dimension(nfiles), save :: lev_files
20	integer, dimension(nfiles), save :: nid_files
21	! integer, save :: nid_hf3d
22
23	!!!!!!!!!! stockage dans une region limitee pour chaque fichier !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
24	! entre [lonmin_reg,lonmax_reg] et [latmin_reg,latmax_reg]
25
26	logical, dimension(nfiles), save :: ok_reglim = (/ .false., .true., .true., .false. /)
27	real, dimension(nfiles), save :: lonmin_reg = (/ 0., 0., 0., 0. /)
28	real, dimension(nfiles), save :: lonmax_reg = (/ 90., 90., 90., 90. /)
29	real, dimension(nfiles), save :: latmin_reg = (/ 0., 0., 0., 0. /)
30	real, dimension(nfiles), save :: latmax_reg = (/ 90., 90., 90., 90. /)
31
32	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
33	!! Definition pour chaque variable du niveau d ecriture dans chaque fichier
34	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!/ histmth, histday, histhf, histins /)!!!!!!!!!!!!
35	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
36
37	!!! 1D
38	integer, dimension(nfiles) , save :: flag_phis = (/ 1, 1, 10, 1 /)
39	integer, dimension(nfiles) , save :: flag_aire = (/ 1, 1, 10, 1 /)
40	integer, dimension(nfiles) , save :: flag_contfracATM = (/ 10, 1, 1, 10 /)
41	integer, dimension(nfiles) , save :: flag_contfracOR = (/ 10, 1, 1, 10 /)
42	integer, dimension(nfiles) , save :: flag_aireTER = (/ 10, 10, 1, 10 /)
43
44	!!! 2D
45	integer, dimension(nfiles) , save :: flag_flat = (/ 10, 1, 10, 10 /)
46	integer, dimension(nfiles) , save :: flag_slp = (/ 1, 1, 1, 10 /)
47	integer, dimension(nfiles) , save :: flag_tsol = (/ 1, 1, 1, 1 /)
48	integer, dimension(nfiles) , save :: flag_t2m = (/ 1, 1, 1, 1 /)
49	integer, dimension(nfiles) , save :: flag_t2m_min = (/ 1, 1, 10, 10 /)
50	integer, dimension(nfiles) , save :: flag_t2m_max = (/ 1, 1, 10, 10 /)
51	integer, dimension(nfiles) , save :: flag_t2m_sol = (/ 10, 4, 10, 10 /)
52	integer, dimension(nfiles) , save :: flag_wind10m = (/ 1, 1, 1, 10 /)
53	integer, dimension(nfiles) , save :: flag_wind10max = (/ 10, 1, 10, 10 /)
54	integer, dimension(nfiles) , save :: flag_sicf = (/ 1, 1, 10, 10 /)
55	integer, dimension(nfiles) , save :: flag_q2m = (/ 1, 1, 1, 1 /)
56	integer, dimension(nfiles) , save :: flag_u10m = (/ 1, 1, 1, 1 /)
57	integer, dimension(nfiles) , save :: flag_v10m = (/ 1, 1, 1, 1 /)
58	integer, dimension(nfiles) , save :: flag_psol = (/ 1, 1, 1, 1 /)
59	integer, dimension(nfiles) , save :: flag_qsurf = (/ 1, 10, 10, 10 /)
60
61	integer, dimension(nfiles) , save :: flag_u10m_sol = (/ 10, 4, 10, 10 /)
62	integer, dimension(nfiles) , save :: flag_v10m_sol = (/ 10, 4, 10, 10 /)
63
64	integer, dimension(nfiles) , save :: flag_qsol = (/ 1, 10, 10, 1 /)
65
66	integer, dimension(nfiles),save :: flag_ndayrain = (/ 1, 10, 10, 10 /)
67	integer, dimension(nfiles),save :: flag_precip(nfiles)= (/ 1, 1, 1, 1 /)
68	integer, dimension(nfiles), save :: flag_plul(nfiles) = (/ 1, 1, 10, 1 /)
69
70	integer, dimension(nfiles) , save :: flag_pluc = (/ 1, 1, 10, 1 /)
71	integer, dimension(nfiles) , save :: flag_snow = (/ 1, 1, 10, 1 /)
72	integer, dimension(nfiles) , save :: flag_evap = (/ 1, 1, 10, 1 /)
73	integer, dimension(nfiles) , save :: flag_tops = (/ 1, 1, 10, 10 /)
74	integer, dimension(nfiles) , save :: flag_tops0 = (/ 1, 5, 10, 10 /)
75	integer, dimension(nfiles) , save :: flag_topl = (/ 1, 1, 10, 1 /)
76	integer, dimension(nfiles) , save :: flag_topl0 = (/ 1, 5, 10, 10 /)
77	integer, dimension(nfiles) , save :: flag_SWupTOA = (/ 1, 4, 10, 10 /)
78	integer, dimension(nfiles) , save :: flag_SWupTOAclr = (/ 1, 4, 10, 10 /)
79	integer, dimension(nfiles) , save :: flag_SWdnTOA = (/ 1, 4, 10, 10 /)
80	integer, dimension(nfiles) , save :: flag_SWdnTOAclr = (/ 1, 4, 10, 10 /)
81	integer, dimension(nfiles) , save :: flag_SWup200 = (/ 1, 10, 10, 10 /)
82	integer, dimension(nfiles) , save :: flag_SWup200clr = (/ 10, 1, 10, 10 /)
83	integer, dimension(nfiles) , save :: flag_SWdn200 = (/ 1, 10, 10, 10 /)
84	integer, dimension(nfiles) , save :: flag_SWdn200clr = (/ 10, 1, 10, 10 /)
85
86	! arajouter
87	! integer, dimension(nfiles) , save :: flag_LWupTOA = (/ 1, 4, 10, 10 /)
88	! integer, dimension(nfiles) , save :: flag_LWupTOAclr = (/ 1, 4, 10, 10 /)
89	! integer, dimension(nfiles) , save :: flag_LWdnTOA = (/ 1, 4, 10, 10 /)
90	! integer, dimension(nfiles) , save :: flag_LWdnTOAclr = (/ 1, 4, 10, 10 /)
91
92	integer, dimension(nfiles) , save :: flag_LWup200 = (/ 1, 10, 10, 10 /)
93	integer, dimension(nfiles) , save :: flag_LWup200clr = (/ 1, 10, 10, 10 /)
94	integer, dimension(nfiles) , save :: flag_LWdn200 = (/ 1, 10, 10, 10 /)
95	integer, dimension(nfiles) , save :: flag_LWdn200clr = (/ 1, 10, 10, 10 /)
96	integer, dimension(nfiles) , save :: flag_sols = (/ 1, 1, 10, 1 /)
97	integer, dimension(nfiles) , save :: flag_sols0 = (/ 1, 5, 10, 10 /)
98	integer, dimension(nfiles) , save :: flag_soll = (/ 1, 1, 10, 1 /)
99	integer, dimension(nfiles) , save :: flag_soll0 = (/ 1, 5, 10, 10 /)
100	integer, dimension(nfiles) , save :: flag_radsol = (/ 1, 1, 10, 10 /)
101	integer, dimension(nfiles) , save :: flag_SWupSFC = (/ 1, 4, 10, 10 /)
102	integer, dimension(nfiles) , save :: flag_SWupSFCclr = (/ 1, 4, 10, 10 /)
103	integer, dimension(nfiles) , save :: flag_SWdnSFC = (/ 1, 1, 10, 10 /)
104	integer, dimension(nfiles) , save :: flag_SWdnSFCclr = (/ 1, 4, 10, 10 /)
105	integer, dimension(nfiles) , save :: flag_LWupSFC = (/ 1, 4, 10, 10 /)
106	integer, dimension(nfiles) , save :: flag_LWupSFCclr = (/ 1, 4, 10, 10 /)
107	integer, dimension(nfiles) , save :: flag_LWdnSFC = (/ 1, 4, 10, 10 /)
108	integer, dimension(nfiles) , save :: flag_LWdnSFCclr = (/ 1, 4, 10, 10 /)
109	integer, dimension(nfiles) , save :: flag_bils = (/ 1, 2, 10, 1 /)
110	integer, dimension(nfiles) , save :: flag_sens = (/ 1, 1, 10, 1 /)
111	integer, dimension(nfiles) , save :: flag_fder = (/ 1, 2, 10, 1 /)
112	integer, dimension(nfiles) , save :: flag_ffonte = (/ 1, 10, 10, 10 /)
113	integer, dimension(nfiles) , save :: flag_fqcalving = (/ 1, 10, 10, 10 /)
114	integer, dimension(nfiles) , save :: flag_fqfonte = (/ 1, 10, 10, 10 /)
115
116	integer, dimension(nfiles) , save :: flag_taux_sol = (/ 1, 4, 10, 1 /)
117	integer, dimension(nfiles) , save :: flag_tauy_sol = (/ 1, 4, 10, 1 /)
118
119	integer, dimension(nfiles) , save :: flag_pourc_sol = (/ 1, 4, 10, 1 /)
120	integer, dimension(nfiles) , save :: flag_fract_sol = (/ 1, 4, 10, 1 /)
121	integer, dimension(nfiles) , save :: flag_tsol_sol = (/ 1, 4, 10, 1 /)
122	integer, dimension(nfiles) , save :: flag_sens_sol = (/ 1, 4, 10, 1 /)
123	integer, dimension(nfiles) , save :: flag_lat_sol = (/ 1, 4, 10, 1 /)
124	integer, dimension(nfiles) , save :: flag_flw_sol = (/ 1, 10, 10, 10 /)
125	integer, dimension(nfiles) , save :: flag_fsw_sol = (/ 1, 10, 10, 10 /)
126	integer, dimension(nfiles) , save :: flag_wbils_sol = (/ 1, 10, 10, 10 /)
127	integer, dimension(nfiles) , save :: flag_wbilo_sol = (/ 1, 10, 10, 10 /)
128
129	integer, dimension(nfiles) , save :: flag_cdrm = (/ 1, 10, 10, 1 /)
130	integer, dimension(nfiles) , save :: flag_cdrh = (/ 1, 10, 10, 1 /)
131	integer, dimension(nfiles) , save :: flag_cldl = (/ 1, 1, 10, 10 /)
132	integer, dimension(nfiles) , save :: flag_cldm = (/ 1, 1, 10, 10 /)
133	integer, dimension(nfiles) , save :: flag_cldh = (/ 1, 1, 10, 10 /)
134	integer, dimension(nfiles) , save :: flag_cldt = (/ 1, 1, 2, 10 /)
135	integer, dimension(nfiles) , save :: flag_cldq = (/ 1, 1, 10, 10 /)
136	integer, dimension(nfiles) , save :: flag_lwp = (/ 1, 5, 10, 10 /)
137	integer, dimension(nfiles) , save :: flag_iwp = (/ 1, 5, 10, 10 /)
138	integer, dimension(nfiles) , save :: flag_ue = (/ 1, 10, 10, 10 /)
139	integer, dimension(nfiles) , save :: flag_ve = (/ 1, 10, 10, 10 /)
140	integer, dimension(nfiles) , save :: flag_uq = (/ 1, 10, 10, 10 /)
141	integer, dimension(nfiles) , save :: flag_vq = (/ 1, 10, 10, 10 /)
142
143	integer, dimension(nfiles) , save :: flag_cape = (/ 1, 10, 10, 10 /)
144	integer, dimension(nfiles) , save :: flag_pbase = (/ 1, 10, 10, 10 /)
145	integer, dimension(nfiles) , save :: flag_ptop = (/ 1, 4, 10, 10 /)
146	integer, dimension(nfiles) , save :: flag_fbase = (/ 1, 10, 10, 10 /)
147	integer, dimension(nfiles) , save :: flag_prw = (/ 1, 1, 10, 10 /)
148
149	integer, dimension(nfiles) , save :: flag_s_pblh = (/ 1, 10, 10, 1 /)
150	integer, dimension(nfiles) , save :: flag_s_pblt = (/ 1, 10, 10, 1 /)
151	integer, dimension(nfiles) , save :: flag_s_lclt = (/ 1, 10, 10, 1 /)
152	integer, dimension(nfiles) , save :: flag_s_capCL = (/ 1, 10, 10, 1 /)
153	integer, dimension(nfiles) , save :: flag_s_oliqCLL = (/ 1, 10, 10, 1 /)
154	integer, dimension(nfiles) , save :: flag_s_cteiCL = (/ 1, 10, 10, 1 /)
155	integer, dimension(nfiles) , save :: flag_s_therm = (/ 1, 10, 10, 1 /)
156	integer, dimension(nfiles) , save :: flag_s_trmb1 = (/ 1, 10, 10, 1 /)
157	integer, dimension(nfiles) , save :: flag_s_trmb2 = (/ 1, 10, 10, 1 /)
158	integer, dimension(nfiles) , save :: flag_s_trmb3 = (/ 1, 10, 10, 1 /)
159
160	integer, dimension(nfiles) , save :: flag_slab_bils = (/ 1, 1, 10, 10 /)
161
162	! Champs interpolles sur des niveaux de pression ??? a faire correctement
163	! if=1 on ecrit u v w phi sur 850 700 500 200 au niv 1
164	! if=2 on ecrit w et ph 500 seulement au niv 1
165	! et u v sur 850 700 500 200
166	! if=3 on ecrit ph a 500 seulement au niv 1
167	! on ecrit u v t q a 850 700 500 200 au niv 3
168	! on ecrit ph a 500 au niv 3
169
170	integer, dimension(nfiles) , save :: flag_ulevsSTD = (/ 1, 1, 3, 10 /)
171	integer, dimension(nfiles) , save :: flag_vlevsSTD = (/ 1, 1, 3, 10 /)
172	integer, dimension(nfiles) , save :: flag_wlevsSTD = (/ 1, 1, 10, 10 /)
173	integer, dimension(nfiles) , save :: flag_tlevsSTD = (/ 10, 10, 3, 10 /)
174	integer, dimension(nfiles) , save :: flag_qlevsSTD = (/ 10, 10, 3, 10 /)
175	integer, dimension(nfiles) , save :: flag_philevsSTD = (/ 1, 1, 1, 10 /)
176
177	integer, dimension(nfiles) , save :: flag_fluxo = (/ 1, 1, 10, 10 /)
178	integer, dimension(nfiles) , save :: flag_fluxg = (/ 1, 1, 10, 10 /)
179	integer, dimension(nfiles) , save :: flag_t_oce_sic = (/ 1, 10, 10, 10 /)
180
181	integer, dimension(nfiles) , save :: flag_lmt_bils = (/ 1, 1, 10, 10 /)
182	integer, dimension(nfiles) , save :: flag_tslab = (/ 1, 1, 10, 10 /)
183	integer, dimension(nfiles) , save :: flag_seaice = (/ 1, 1, 10, 10 /)
184	integer, dimension(nfiles) , save :: flag_siceh = (/ 1, 1, 10, 10 /)
185
186	integer, dimension(nfiles) , save :: flag_weakinv = (/ 10, 1, 10, 10 /)
187	integer, dimension(nfiles) , save :: flag_dthmin = (/ 10, 1, 10, 10 /)
188	integer, dimension(nfiles) , save :: flag_u10_sol = (/ 10, 4, 10, 10 /)
189	integer, dimension(nfiles) , save :: flag_v10_sol = (/ 10, 4, 10, 10 /)
190	integer, dimension(nfiles) , save :: flag_meantaucld = (/ 10, 5, 10, 10 /)
191	integer, dimension(nfiles) , save :: flag_cldtau = (/ 10, 5, 10, 10 /)
192	integer, dimension(nfiles) , save :: flag_cldemi = (/ 10, 5, 10, 10 /)
193	integer, dimension(nfiles) , save :: flag_rh2m = (/ 10, 5, 10, 10 /)
194	integer, dimension(nfiles) , save :: flag_qsat2m = (/ 10, 5, 10, 10 /)
195	integer, dimension(nfiles) , save :: flag_tpot = (/ 10, 5, 10, 10 /)
196	integer, dimension(nfiles) , save :: flag_tpote = (/ 10, 5, 10, 10 /)
197	integer, dimension(nfiles) , save :: flag_tke = (/ 4, 10, 10, 10 /)
198	integer, dimension(nfiles) , save :: flag_tke_max = (/ 4, 10, 10, 10 /)
199	integer, dimension(nfiles) , save :: flag_tke_sol = (/ 10, 4, 10, 10 /)
200	integer, dimension(nfiles) , save :: flag_tke_max_sol = (/ 10, 4, 10, 10 /)
201	integer, dimension(nfiles) , save :: flag_kz = (/ 4, 10, 10, 10 /)
202	integer, dimension(nfiles) , save :: flag_kz_max = (/ 4, 10, 10, 10 /)
203	integer, dimension(nfiles) , save :: flag_SWnetOR = (/ 10, 10, 2, 10 /)
204	integer, dimension(nfiles) , save :: flag_SWdownOR = (/ 10, 10, 2, 10 /)
205	integer, dimension(nfiles) , save :: flag_LWdownOR = (/ 10, 10, 2, 10 /)
206
207	integer, dimension(nfiles) , save :: flag_snowl = (/ 10, 1, 10, 10 /)
208	integer, dimension(nfiles) , save :: flag_cape_max = (/ 10, 1, 10, 10 /)
209	integer, dimension(nfiles) , save :: flag_solldown = (/ 10, 1, 10, 1 /)
210
211	integer, dimension(nfiles) , save :: flag_dtsvdfo = (/ 10, 10, 10, 1 /)
212	integer, dimension(nfiles) , save :: flag_dtsvdft = (/ 10, 10, 10, 1 /)
213	integer, dimension(nfiles) , save :: flag_dtsvdfg = (/ 10, 10, 10, 1 /)
214	integer, dimension(nfiles) , save :: flag_dtsvdfi = (/ 10, 10, 10, 1 /)
215	integer, dimension(nfiles) , save :: flag_rugs = (/ 10, 10, 10, 1 /)
216
217	!!!!!!!!!!!!!!!!!!!!!! 3D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
218	integer, dimension(nfiles) , save :: flag_lwcon = (/ 2, 5, 10, 10 /)
219	integer, dimension(nfiles) , save :: flag_iwcon = (/ 2, 5, 10, 10 /)
220	integer, dimension(nfiles) , save :: flag_temp = (/ 2, 3, 4, 1 /)
221	integer, dimension(nfiles) , save :: flag_ovap = (/ 2, 3, 4, 10 /)
222	integer, dimension(nfiles) , save :: flag_wvapp = (/ 2, 10, 10, 10 /)
223	integer, dimension(nfiles) , save :: flag_geop = (/ 2, 3, 10, 1 /)
224	integer, dimension(nfiles) , save :: flag_vitu = (/ 2, 3, 4, 1 /)
225	integer, dimension(nfiles) , save :: flag_vitv = (/ 2, 3, 4, 1 /)
226	integer, dimension(nfiles) , save :: flag_vitw = (/ 2, 3, 10, 10 /)
227	integer, dimension(nfiles) , save :: flag_pres = (/ 2, 3, 10, 1 /)
228	integer, dimension(nfiles) , save :: flag_rneb = (/ 2, 5, 10, 10 /)
229	integer, dimension(nfiles) , save :: flag_rnebcon = (/ 2, 5, 10, 10 /)
230	integer, dimension(nfiles) , save :: flag_rhum = (/ 2, 10, 10, 10 /)
231	integer, dimension(nfiles) , save :: flag_ozone = (/ 2, 10, 10, 10 /)
232	integer, dimension(nfiles) , save :: flag_upwd = (/ 2, 10, 10, 10 /)
233	integer, dimension(nfiles) , save :: flag_dtphy = (/ 2, 10, 10, 10 /)
234	integer, dimension(nfiles) , save :: flag_dqphy = (/ 2, 10, 10, 10 /)
235	! arevoir probleme noms (geop, ovap a ?????)
236
237	integer, dimension(nfiles) , save :: flag_albe_sol = (/ 3, 4, 10, 1 /)
238	integer, dimension(nfiles) , save :: flag_ages_sol = (/ 3, 10, 10, 10 /)
239	integer, dimension(nfiles) , save :: flag_rugs_sol = (/ 3, 4, 10, 1 /)
240
241	integer, dimension(nfiles) , save :: flag_albs = (/ 3, 10, 10, 1 /)
242	integer, dimension(nfiles) , save :: flag_albslw = (/ 3, 10, 10, 1 /)
243
244	integer, dimension(nfiles) , save :: flag_clwcon = (/ 4, 10, 10, 10 /)
245	integer, dimension(nfiles) , save :: flag_Ma = (/ 4, 10, 10, 10 /)
246	integer, dimension(nfiles) , save :: flag_dnwd = (/ 4, 10, 10, 10 /)
247	integer, dimension(nfiles) , save :: flag_dnwd0 = (/ 4, 10, 10, 10 /)
248	integer, dimension(nfiles) , save :: flag_dtdyn = (/ 4, 10, 10, 10 /)
249	integer, dimension(nfiles) , save :: flag_dqdyn = (/ 4, 10, 10, 10 /)
250	integer, dimension(nfiles) , save :: flag_dtcon = (/ 4, 5, 10, 10 /)
251	integer, dimension(nfiles) , save :: flag_ducon = (/ 4, 10, 10, 10 /)
252	integer, dimension(nfiles) , save :: flag_dqcon = (/ 4, 5, 10, 10 /)
253	integer, dimension(nfiles) , save :: flag_dtlsc = (/ 4, 10, 10, 10 /)
254	integer, dimension(nfiles) , save :: flag_dtlschr = (/ 4, 10, 10, 10 /)
255	integer, dimension(nfiles) , save :: flag_dqlsc = (/ 4, 10, 10, 10 /)
256	integer, dimension(nfiles) , save :: flag_dtvdf = (/ 4, 10, 10, 1 /)
257	integer, dimension(nfiles) , save :: flag_dqvdf = (/ 4, 10, 10, 1 /)
258	integer, dimension(nfiles) , save :: flag_dteva = (/ 4, 10, 10, 10 /)
259	integer, dimension(nfiles) , save :: flag_dqeva = (/ 4, 10, 10, 10 /)
260	integer, dimension(nfiles) , save :: flag_ptconv = (/ 4, 10, 10, 10 /)
261	integer, dimension(nfiles) , save :: flag_ratqs = (/ 4, 10, 10, 10 /)
262	integer, dimension(nfiles) , save :: flag_dtthe = (/ 4, 10, 10, 10 /)
263	integer, dimension(nfiles) , save :: flag_dqthe = (/ 4, 10, 10, 10 /)
264	integer, dimension(nfiles) , save :: flag_dtajs = (/ 4, 10, 10, 10 /)
265	integer, dimension(nfiles) , save :: flag_dqajs = (/ 4, 10, 10, 10 /)
266	integer, dimension(nfiles) , save :: flag_dtswr = (/ 4, 10, 10, 10 /)
267	integer, dimension(nfiles) , save :: flag_dtsw0 = (/ 4, 10, 10, 10 /)
268	integer, dimension(nfiles) , save :: flag_dtlwr = (/ 4, 10, 10, 10 /)
269	integer, dimension(nfiles) , save :: flag_dtlw0 = (/ 4, 10, 10, 10 /)
270	integer, dimension(nfiles) , save :: flag_dtec = (/ 4, 10, 10, 10 /)
271	integer, dimension(nfiles) , save :: flag_duvdf = (/ 4, 10, 10, 10 /)
272	integer, dimension(nfiles) , save :: flag_dvvdf = (/ 4, 10, 10, 10 /)
273	integer, dimension(nfiles) , save :: flag_duoro = (/ 4, 10, 10, 10 /)
274	integer, dimension(nfiles) , save :: flag_dvoro = (/ 4, 10, 10, 10 /)
275	integer, dimension(nfiles) , save :: flag_dulif = (/ 4, 10, 10, 10 /)
276	integer, dimension(nfiles) , save :: flag_dvlif = (/ 4, 10, 10, 10 /)
277	integer, dimension(nfiles) , save :: flag_trac = (/ 4, 10, 10, 10 /)
278
279	CONTAINS
280
281	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
282	!!!!!!!!! Ouverture des fichier et definition des variable de sortie !!!!!!!!
283	!! histbeg, histvert et histdef
284	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
285
286	SUBROUTINE phys_output_open ( jj_nb, jjmp1, nqmax, &
287	& nlevSTD, clevSTD, nbteta, ctetaSTD, &
288	& dtime, presnivs, ok_veget, ocean, iflag_pbl, &
289	& ok_mensuel,ok_journe,ok_hf,ok_instan, &
290	& rlon, rlat, nid_fil)
291
292
293
294	USE iophy
295	USE dimphy
296	USE ioipsl
297
298	IMPLICIT NONE
299	include "dimensions.h"
300	include "temps.h"
301	include "indicesol.h"
302	include "advtrac.h"
303	include "clesphys.h"
304
305	integer :: jj_nb, jjmp1, nqmax
306	integer :: nbteta, nlevSTD, radpas
307	logical :: ok_mensuel, ok_journe, ok_hf, ok_instan
308	real :: dtime, zstophy, zout
309	integer :: nhori, nvert
310	integer :: idayref
311	real :: zjulian
312	real, dimension(klev) :: presnivs
313	character(len=4), dimension(nlevSTD) :: clevSTD
314	integer :: nsrf, k, iq, iiq, iff, i, j
315	logical :: ok_veget
316	integer :: iflag_pbl
317	CHARACTER(len=3) :: bb2
318	CHARACTER(len=2) :: bb3
319	character(len=6) :: ocean
320	CHARACTER(len=3) :: ctetaSTD(nbteta)
321	integer, dimension(nfiles) :: nid_fil
322	real, dimension(nfiles) :: ecrit_files
323	CHARACTER(len=20), dimension(nfiles) :: name_files
324	CHARACTER(len=20), dimension(nfiles) :: type_ecri
325	INTEGER, dimension(iim*jjmp1) :: ndex2d
326	INTEGER, dimension(iimjjmp1klev) :: ndex3d
327	REAL, dimension(klon) :: rlon, rlat
328	integer :: imin_ins, imax_ins
329	integer :: jmin_ins, jmax_ins
330	REAL, dimension(iim,jjmp1) :: zx_lon, zx_lat
331
332	name_files(1) = 'histmth'
333	name_files(2) = 'histday'
334	name_files(3) = 'histhf'
335	name_files(4) = 'histins'
336
337	type_ecri(1) = 'ave(X)'
338	type_ecri(2) = 'ave(X)'
339	type_ecri(3) = 'ave(X)'
340	type_ecri(4) = 'inst(X)'
341
342	clef_files(1) = ok_mensuel
343	clef_files(2) = ok_journe
344	clef_files(3) = ok_hf
345	clef_files(4) = ok_instan
346
347	lev_files(1) = lev_histmth
348	lev_files(2) = lev_histday
349	lev_files(3) = lev_histhf
350	lev_files(4) = 1
351
352	ecrit_files(1) = ecrit_mth
353	ecrit_files(2) = ecrit_day
354	ecrit_files(3) = ecrit_hf
355	ecrit_files(4) = ecrit_ins
356
357	!!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
358	! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !!
359	! Appel des histbeg pour definir les variables (nom, moy ou inst, freq de sortie ..
360	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
361
362	DO iff=1,nfiles
363
364	IF (clef_files(iff)) THEN
365
366	zstophy = dtime ! Frequence ou l on moyenne
367	zout = ecrit_files(iff) ! Frequence ou l on ecrit
368
369	idayref = day_ref
370
371	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
372
373	! print*,'Dans hist_def, name_files(iff) ',name_files(iff)
374
375	!!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
376	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
377	if (ok_reglim(iff)) then
378
379	! a reecrire plus simplement en appelant io_lat et io_lon de iophy
380
381	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
382	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
383
384	DO i = 1, iim
385	zx_lon(i,1) = rlon(i+1)
386	zx_lon(i,jjmp1) = rlon(i+1)
387	ENDDO
388
389
390	imin_ins=1
391	imax_ins=iim
392	jmin_ins=1
393	jmax_ins=jjmp1
394
395	do i=1,iim-1
396	if (zx_lon(i,1).lt.lonmin_reg(iff)) imin_ins=i
397	if (zx_lon(i,1).le.lonmax_reg(iff)) imax_ins=i+1
398	enddo
399
400	do j=1,jjmp1-1
401	if (zx_lat(1,j).ge.latmin_reg(iff)) jmax_ins=j
402	if (zx_lat(1,j).gt.latmax_reg(iff)) jmin_ins=j
403	enddo
404
405
406	print*,'On stoke le fichier hist sur, ', &
407	& imin_ins,imax_ins,jmin_ins,jmax_ins
408	print*,'On stoke le fichier instantanne sur, ', &
409	& zx_lon(imin_ins,1),zx_lon(imax_ins,1), &
410	& zx_lat(1,jmin_ins),zx_lat(1,jmax_ins)
411
412
413	CALL histbeg(name_files(iff),iim,zx_lon(:,1), &
414	& jjmp1,zx_lat(1,:), &
415	& imin_ins,imax_ins-imin_ins+1, &
416	& jmin_ins,jmax_ins-jmin_ins+1, &
417	& itau_phy, zjulian, dtime, &
418	& nhori, nid_fil(iff))
419
420	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
421	else
422
423	CALL histbeg_phy(name_files(iff),itau_phy, &
424	& zjulian, dtime, &
425	& nhori, nid_fil(iff))
426	endif
427
428	nid_files(iff)=nid_fil(iff)
429	! print*,'Dans hist_def, nid_files(iff) ',nid_files(iff)
430
431	CALL histvert(nid_files(iff),"presnivs", &
432	& "Vertical levels","mb", &
433	& klev, presnivs/100., nvert)
434
435	!!!!!!!!!!!!! Traitement des champs 3D pour histhf !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
436	!!!!!!!!!!!!!!! A Revoir plus tard !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
437	! IF (iff.eq.3.and.lev_files(iff).ge.4) THEN
438	! CALL histbeg_phy("histhf3d",itau_phy, &
439	! & zjulian, dtime, &
440	! & nhori, nid_hf3d)
441
442	! CALL histvert(nid_hf3d,"presnivs", &
443	! & "Vertical levels","mb", &
444	! & klev, presnivs/100., nvert)
445	! ENDIF
446
447
448	!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
449	IF (flag_phis(iff)<=lev_files(iff)) THEN
450	CALL histdef(nid_files(iff),"phis", &
451	& "Surface geop. height", "m2/s2", &
452	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
453	& "once", zstophy,zout)
454	ENDIF
455
456	IF (flag_aire(iff)<=lev_files(iff)) THEN
457	CALL histdef(nid_files(iff), "aire", "Grid area", "-", &
458	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
459	& "once", zstophy,zout)
460	ENDIF
461
462
463	IF (flag_contfracATM(iff)<=lev_files(iff)) THEN
464	CALL histdef(nid_files(iff), &
465	& "contfracATM","% sfce ter+lic ", &
466	& "-", iim,jj_nb,nhori, 1,1,1, -99, 32, &
467	& "once", zstophy,zout)
468	ENDIF
469
470	IF (flag_contfracOR(iff)<=lev_files(iff)) THEN
471	CALL histdef(nid_files(iff), &
472	& "contfracOR","% sfce terre OR", &
473	& "-", iim,jj_nb,nhori, 1,1,1, -99, 32, &
474	& "once", zstophy,zout)
475	ENDIF
476
477	IF (flag_aireTER(iff)<=lev_files(iff)) THEN
478	CALL histdef(nid_files(iff), &
479	& "aireTER","Grid area CONT","-", &
480	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
481	& "once", zstophy,zout)
482	ENDIF
483
484	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
485
486	IF (flag_flat(iff)<=lev_files(iff)) THEN
487	CALL histdef(nid_files(iff), "flat", "Latent heat flux", &
488	& "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, &
489	& type_ecri(iff), zstophy,zout)
490	ENDIF
491
492	IF (flag_slp(iff)<=lev_files(iff)) THEN
493	CALL histdef(nid_files(iff), "slp", &
494	& "Sea Level Pressure", "Pa", &
495	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
496	& type_ecri(iff), zstophy,zout)
497	ENDIF
498
499	IF (flag_tsol(iff)<=lev_files(iff)) THEN
500	CALL histdef(nid_files(iff), "tsol", &
501	& "Surface Temperature", "K", &
502	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
503	& type_ecri(iff), zstophy,zout)
504	ENDIF
505
506	IF (flag_t2m(iff)<=lev_files(iff)) THEN
507	CALL histdef(nid_files(iff), "t2m", &
508	& "Temperature 2m", "K", &
509	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
510	& type_ecri(iff), zstophy,zout)
511	ENDIF
512
513	IF (flag_t2m_min(iff)<=lev_files(iff)) THEN
514	CALL histdef(nid_files(iff), "t2m_min", &
515	& "Temp 2m min", "K", &
516	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
517	& "t_min(X)", zstophy,zout)
518	ENDIF
519
520	IF (flag_t2m_max(iff)<=lev_files(iff)) THEN
521	CALL histdef(nid_files(iff), "t2m_max", &
522	& "Temp 2m max", "K", &
523	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
524	& "t_max(X)", zstophy,zout)
525	ENDIF
526
527	IF (flag_wind10m(iff)<=lev_files(iff)) THEN
528	CALL histdef(nid_files(iff), "wind10m", &
529	& "10-m wind speed","m/s", &
530	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
531	& type_ecri(iff), zstophy,zout)
532	ENDIF
533
534	IF (flag_wind10max(iff)<=lev_files(iff)) THEN
535	CALL histdef(nid_files(iff), "wind10max", &
536	& "10-m wind speed max.", "m/s", &
537	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
538	& "t_max(X)", zstophy,zout)
539	ENDIF
540
541	IF (flag_sicf(iff)<=lev_files(iff)) THEN
542	CALL histdef(nid_files(iff), "sicf", &
543	& "Sea-ice fraction", "-", &
544	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
545	& type_ecri(iff), zstophy,zout)
546	ENDIF
547
548	IF (flag_q2m(iff)<=lev_files(iff)) THEN
549	CALL histdef(nid_files(iff), "q2m", &
550	& "Specific humidity 2m", &
551	& "kg/kg", iim,jj_nb,nhori, &
552	& 1,1,1, -99, 32, &
553	& type_ecri(iff), zstophy,zout)
554	ENDIF
555
556	IF (flag_u10m(iff)<=lev_files(iff)) THEN
557	CALL histdef(nid_files(iff), "u10m", &
558	& "Vent zonal 10m", "m/s", &
559	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
560	& type_ecri(iff), zstophy,zout)
561	ENDIF
562
563	IF (flag_v10m(iff)<=lev_files(iff)) THEN
564	CALL histdef(nid_files(iff), "v10m", &
565	& "Vent meridien 10m", &
566	& "m/s", iim,jj_nb,nhori, &
567	& 1,1,1, -99, 32, &
568	& type_ecri(iff), zstophy,zout)
569	ENDIF
570
571	IF (flag_psol(iff)<=lev_files(iff)) THEN
572	CALL histdef(nid_files(iff),"psol", &
573	& "Surface Pressure","Pa", &
574	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
575	& type_ecri(iff), zstophy,zout)
576	ENDIF
577
578	IF (flag_qsurf(iff)<=lev_files(iff)) THEN
579	CALL histdef(nid_files(iff),"qsurf", &
580	& "Surface Air humidity", &
581	& "kg/kg",iim,jj_nb,nhori, 1,1,1, -99, 32, &
582	& type_ecri(iff), zstophy,zout)
583	ENDIF
584
585	if (.not. ok_veget) then
586	IF (flag_qsol(iff)<=lev_files(iff)) THEN
587	CALL histdef(nid_files(iff), "qsol", &
588	& "Soil watter content", &
589	& "mm", iim,jj_nb,nhori, &
590	& 1,1,1, -99, 32, &
591	& type_ecri(iff), zstophy,zout)
592	ENDIF
593	endif
594
595	IF (flag_ndayrain(iff)<=lev_files(iff)) THEN
596	CALL histdef(nid_files(iff), "ndayrain", &
597	& "Number of day with rain (liq+sol)", "-", &
598	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
599	& "inst(X)", zout,zout)
600	ENDIF
601
602	IF (flag_precip(iff)<=lev_files(iff)) THEN
603	CALL histdef(nid_files(iff),"precip", &
604	& "Precip Totale liq+sol","kg/(s*m2)", &
605	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
606	& type_ecri(iff), zstophy,zout)
607	ENDIF
608
609	IF (flag_plul(iff)<=lev_files(iff)) THEN
610	CALL histdef(nid_files(iff), "plul", &
611	& "Large-scale Precip.","kg/(s*m2)", &
612	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
613	& type_ecri(iff), zstophy,zout)
614	ENDIF
615
616	IF (flag_pluc(iff)<=lev_files(iff)) THEN
617	CALL histdef(nid_files(iff), "pluc", &
618	& "Convective Precip.", "kg/(s*m2)", &
619	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
620	& type_ecri(iff), zstophy,zout)
621	ENDIF
622
623	IF (flag_snow(iff)<=lev_files(iff)) THEN
624	CALL histdef(nid_files(iff),"snow", &
625	& "Snow fall","kg/(s*m2)", &
626	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
627	& type_ecri(iff), zstophy,zout)
628	ENDIF
629
630	IF (flag_evap(iff)<=lev_files(iff)) THEN
631	CALL histdef(nid_files(iff),"evap", &
632	& "Evaporat", "kg/(s*m2)", &
633	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
634	& type_ecri(iff), zstophy,zout)
635	ENDIF
636
637	IF (flag_tops(iff)<=lev_files(iff)) THEN
638	CALL histdef(nid_files(iff),"tops", &
639	& "Solar rad. at TOA","W/m2", &
640	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
641	& type_ecri(iff), zstophy,zout)
642	ENDIF
643
644	IF (flag_tops0(iff)<=lev_files(iff)) THEN
645	CALL histdef(nid_files(iff),"tops0", &
646	& "CS Solar rad. at TOA", &
647	& "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, &
648	& type_ecri(iff), zstophy,zout)
649	ENDIF
650
651	IF (flag_topl(iff)<=lev_files(iff)) THEN
652	CALL histdef(nid_files(iff), "topl", &
653	& "IR rad. at TOA", &
654	& "W/m2",iim,jj_nb,nhori, 1,1,1, -99, 32, &
655	& type_ecri(iff), zstophy,zout)
656	ENDIF
657
658	IF (flag_topl0(iff)<=lev_files(iff)) THEN
659	CALL histdef(nid_files(iff),"topl0", &
660	& "IR rad. at TOA","W/m2", &
661	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
662	& type_ecri(iff), zstophy,zout)
663	ENDIF
664
665	IF (flag_SWupTOA(iff)<=lev_files(iff)) THEN
666	CALL histdef(nid_files(iff),"SWupTOA", &
667	& "SWup at TOA","W/m2", &
668	& iim,jj_nb,nhori, 1,1,1,-99, &
669	& 32, type_ecri(iff), zstophy,zout)
670	ENDIF
671
672	IF (flag_SWupTOAclr(iff)<=lev_files(iff)) THEN
673	CALL histdef(nid_files(iff), "SWupTOAclr", &
674	& "SWup clear sky at TOA","W/m2", &
675	& iim,jj_nb,nhori, 1,1,1,-99, &
676	& 32, type_ecri(iff), zstophy,zout)
677	ENDIF
678
679	IF (flag_SWdnTOA(iff)<=lev_files(iff)) THEN
680	CALL histdef(nid_files(iff),"SWdnTOA", &
681	& "SWdn at TOA","W/m2", &
682	& iim,jj_nb,nhori, 1,1,1,-99, &
683	& 32, type_ecri(iff), zstophy,zout)
684	ENDIF
685
686	IF (flag_SWdnTOAclr(iff)<=lev_files(iff)) THEN
687	CALL histdef(nid_files(iff), "SWdnTOAclr", &
688	& "SWdn clear sky at TOA","W/m2", &
689	& iim,jj_nb,nhori, 1,1,1,-99, &
690	& 32, type_ecri(iff), zstophy,zout)
691	ENDIF
692
693	IF (flag_SWup200(iff)<=lev_files(iff)) THEN
694	CALL histdef(nid_files(iff),"SWup200", &
695	& "SWup at 200mb","W/m2", &
696	& iim,jj_nb,nhori, 1,1,1,-99, &
697	& 32, type_ecri(iff), zstophy,zout)
698	ENDIF
699
700	IF (flag_SWup200clr(iff)<=lev_files(iff)) THEN
701	CALL histdef(nid_files(iff), "SWup200clr", &
702	& "SWup clear sky at 200mb","W/m2", &
703	& iim,jj_nb,nhori, 1,1,1,-99, &
704	& 32, type_ecri(iff), zstophy,zout)
705	ENDIF
706
707	IF (flag_SWdn200(iff)<=lev_files(iff)) THEN
708	CALL histdef(nid_files(iff),"SWdn200", &
709	& "SWdn at 200mb","W/m2", &
710	& iim,jj_nb,nhori, 1,1,1,-99, &
711	& 32, type_ecri(iff), zstophy,zout)
712	ENDIF
713
714	IF (flag_SWdn200clr(iff)<=lev_files(iff)) THEN
715	CALL histdef(nid_files(iff), "SWdn200clr", &
716	& "SWdn clear sky at 200mb","W/m2", &
717	& iim,jj_nb,nhori, 1,1,1,-99, &
718	& 32, type_ecri(iff), zstophy,zout)
719	ENDIF
720
721	IF (flag_LWup200(iff)<=lev_files(iff)) THEN
722	CALL histdef(nid_files(iff),"LWup200", &
723	& "LWup at 200mb","W/m2", &
724	& iim,jj_nb,nhori, 1,1,1,-99, &
725	& 32, type_ecri(iff), zstophy,zout)
726	ENDIF
727
728	IF (flag_LWup200clr(iff)<=lev_files(iff)) THEN
729	CALL histdef(nid_files(iff), "LWup200clr", &
730	& "LWup clear sky at 200mb","W/m2", &
731	& iim,jj_nb,nhori, 1,1,1,-99, &
732	& 32, type_ecri(iff), zstophy,zout)
733	ENDIF
734
735	IF (flag_LWdn200(iff)<=lev_files(iff)) THEN
736	CALL histdef(nid_files(iff),"LWdn200", &
737	& "LWdn at 200mb","W/m2", &
738	& iim,jj_nb,nhori, 1,1,1,-99, &
739	& 32, type_ecri(iff), zstophy,zout)
740	ENDIF
741
742	IF (flag_LWdn200clr(iff)<=lev_files(iff)) THEN
743	CALL histdef(nid_files(iff), "LWdn200clr", &
744	& "LWdn clear sky at 200mb","W/m2", &
745	& iim,jj_nb,nhori, 1,1,1,-99, &
746	& 32, type_ecri(iff), zstophy,zout)
747	ENDIF
748
749	IF (flag_sols(iff)<=lev_files(iff)) THEN
750	CALL histdef(nid_files(iff), "sols", &
751	& "Solar rad. at surf.", &
752	& "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, &
753	& type_ecri(iff), zstophy,zout)
754	ENDIF
755
756	IF (flag_sols0(iff)<=lev_files(iff)) THEN
757	CALL histdef(nid_files(iff), "sols0", &
758	& "Solar rad. at surf.", &
759	& "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, &
760	& type_ecri(iff), zstophy,zout)
761	ENDIF
762
763	IF (flag_soll(iff)<=lev_files(iff)) THEN
764	CALL histdef(nid_files(iff),"soll", &
765	& "IR rad. at surface","W/m2", &
766	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
767	& type_ecri(iff), zstophy,zout)
768	ENDIF
769
770	IF (flag_radsol(iff)<=lev_files(iff)) THEN
771	CALL histdef(nid_files(iff),"radsol", &
772	& "Rayonnement au sol", &
773	& "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, &
774	& type_ecri(iff), zstophy,zout)
775	ENDIF
776
777	IF (flag_soll0(iff)<=lev_files(iff)) THEN
778	CALL histdef(nid_files(iff), "soll0", &
779	& "IR rad. at surface", &
780	& "W/m2",iim,jj_nb,nhori, 1,1,1, -99, 32, &
781	& type_ecri(iff), zstophy,zout)
782	ENDIF
783
784	IF (flag_SWupSFC(iff)<=lev_files(iff)) THEN
785	CALL histdef(nid_files(iff), "SWupSFC", &
786	& "SWup at surface", &
787	& "W/m2",iim,jj_nb,nhori, 1,1,1,-99, &
788	& 32, type_ecri(iff), zstophy,zout)
789	ENDIF
790
791	IF (flag_SWupSFCclr(iff)<=lev_files(iff)) THEN
792	CALL histdef(nid_files(iff), "SWupSFCclr", &
793	& "SWup clear sky at surface","W/m2", &
794	& iim,jj_nb,nhori, 1,1,1,-99, &
795	& 32, type_ecri(iff), zstophy,zout)
796	ENDIF
797
798	IF (flag_SWdnSFC(iff)<=lev_files(iff)) THEN
799	CALL histdef(nid_files(iff),"SWdnSFC", &
800	& "SWdn at surface","W/m2", &
801	& iim,jj_nb,nhori, 1,1,1,-99, &
802	& 32, type_ecri(iff), zstophy,zout)
803	ENDIF
804
805	IF (flag_SWdnSFCclr(iff)<=lev_files(iff)) THEN
806	CALL histdef(nid_files(iff), "SWdnSFCclr", &
807	& "SWdn clear sky at surface","W/m2", &
808	& iim,jj_nb,nhori, 1,1,1,-99, &
809	& 32, type_ecri(iff), zstophy,zout)
810	ENDIF
811
812	IF (flag_LWupSFC(iff)<=lev_files(iff)) THEN
813	CALL histdef(nid_files(iff),"LWupSFC", &
814	& "Upwd. IR rad. at surface", "W/m2", &
815	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
816	& type_ecri(iff), zstophy,zout)
817	ENDIF
818
819	IF (flag_LWdnSFC(iff)<=lev_files(iff)) THEN
820	CALL histdef(nid_files(iff),"LWdnSFC", &
821	& "Down. IR rad. at surface", &
822	& "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, &
823	& type_ecri(iff), zstophy,zout)
824	ENDIF
825
826	IF (flag_LWupSFCclr(iff)<=lev_files(iff)) THEN
827	CALL histdef(nid_files(iff),"LWupSFCclr", &
828	& "CS Upwd. IR rad. at surface", &
829	& "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, &
830	& type_ecri(iff), zstophy,zout)
831	ENDIF
832
833	IF (flag_LWdnSFCclr(iff)<=lev_files(iff)) THEN
834	CALL histdef(nid_files(iff),"LWdnSFCclr", &
835	& "Down. CS IR rad. at surface", &
836	& "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, &
837	& type_ecri(iff), zstophy,zout)
838	ENDIF
839
840	IF (flag_bils(iff)<=lev_files(iff)) THEN
841	CALL histdef(nid_files(iff), &
842	& "bils","Surf. total heat flux","W/m2", &
843	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
844	& type_ecri(iff), zstophy,zout)
845	ENDIF
846
847	IF (flag_sens(iff)<=lev_files(iff)) THEN
848	CALL histdef(nid_files(iff),"sens", &
849	& "Sensible heat flux","W/m2", &
850	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
851	& type_ecri(iff), zstophy,zout)
852	ENDIF
853
854	IF (flag_fder(iff)<=lev_files(iff)) THEN
855	CALL histdef(nid_files(iff),"fder", &
856	& "Heat flux derivation","W/m2", &
857	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
858	& type_ecri(iff), zstophy,zout)
859	ENDIF
860
861	IF (flag_ffonte(iff)<=lev_files(iff)) THEN
862	CALL histdef(nid_files(iff),"ffonte", &
863	& "Thermal flux for snow melting", &
864	& "W/m2",iim,jj_nb,nhori, 1,1,1, -99, 32, &
865	& type_ecri(iff), zstophy,zout)
866	ENDIF
867
868	IF (flag_fqcalving(iff)<=lev_files(iff)) THEN
869	CALL histdef(nid_files(iff), "fqcalving", &
870	& "Ice Calving", &
871	& "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32, &
872	& type_ecri(iff), zstophy,zout)
873	ENDIF
874
875	IF (flag_fqfonte(iff)<=lev_files(iff)) THEN
876	CALL histdef(nid_files(iff), "fqfonte", &
877	& "Land ice melt", &
878	& "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32, &
879	& type_ecri(iff), zstophy,zout)
880	ENDIF
881
882	DO nsrf = 1, nbsrf
883	IF (flag_pourc_sol(iff)<=lev_files(iff)) THEN
884	call histdef(nid_files(iff), "pourc_"//clnsurf(nsrf), &
885	& "% "//clnsurf(nsrf), "%", &
886	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
887	& type_ecri(iff), zstophy,zout)
888	ENDIF
889
890	IF (flag_fract_sol(iff)<=lev_files(iff)) THEN
891	call histdef(nid_files(iff), "fract_"//clnsurf(nsrf), &
892	& "Fraction "//clnsurf(nsrf), "1", &
893	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
894	& type_ecri(iff), zstophy,zout)
895	ENDIF
896
897	IF (flag_taux_sol(iff)<=lev_files(iff)) THEN
898	call histdef(nid_files(iff), "taux_"//clnsurf(nsrf), &
899	& "Zonal wind stress"//clnsurf(nsrf), "Pa", &
900	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
901	& type_ecri(iff), zstophy,zout)
902	ENDIF
903
904	IF (flag_tauy_sol(iff)<=lev_files(iff)) THEN
905	call histdef(nid_files(iff), "tauy_"//clnsurf(nsrf), &
906	& "Meridional wind stress "//clnsurf(nsrf), "Pa", &
907	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
908	& type_ecri(iff), zstophy,zout)
909	ENDIF
910
911	IF (flag_tsol_sol(iff)<=lev_files(iff)) THEN
912	call histdef(nid_files(iff), "tsol_"//clnsurf(nsrf), &
913	& "Temperature "//clnsurf(nsrf), "K", &
914	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
915	& type_ecri(iff), zstophy,zout)
916	ENDIF
917
918	IF (flag_u10m_sol(iff)<=lev_files(iff)) THEN
919	call histdef(nid_files(iff), "u10m_"//clnsurf(nsrf), &
920	& "Vent Zonal 10m "//clnsurf(nsrf), "m/s", &
921	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
922	& type_ecri(iff), zstophy,zout)
923	ENDIF
924
925	IF (flag_v10m_sol(iff)<=lev_files(iff)) THEN
926	call histdef(nid_files(iff), "v10m_"//clnsurf(nsrf), &
927	& "Vent meredien 10m "//clnsurf(nsrf), "m/s", &
928	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
929	& type_ecri(iff), zstophy,zout)
930	ENDIF
931
932	IF (flag_t2m_sol(iff)<=lev_files(iff)) THEN
933	call histdef(nid_files(iff), "t2m_"//clnsurf(nsrf), &
934	& "Temp 2m "//clnsurf(nsrf), "K", &
935	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
936	& type_ecri(iff), zstophy,zout)
937	ENDIF
938
939	IF (flag_sens_sol(iff)<=lev_files(iff)) THEN
940	call histdef(nid_files(iff), "sens_"//clnsurf(nsrf), &
941	& "Sensible heat flux "//clnsurf(nsrf), "W/m2", &
942	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
943	& type_ecri(iff), zstophy,zout)
944	ENDIF
945
946	IF (flag_lat_sol(iff)<=lev_files(iff)) THEN
947	call histdef(nid_files(iff), "lat_"//clnsurf(nsrf), &
948	& "Latent heat flux "//clnsurf(nsrf), "W/m2", &
949	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
950	& type_ecri(iff), zstophy,zout)
951	ENDIF
952
953	IF (flag_flw_sol(iff)<=lev_files(iff)) THEN
954	call histdef(nid_files(iff), "flw_"//clnsurf(nsrf), &
955	& "LW "//clnsurf(nsrf), "W/m2", &
956	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
957	& type_ecri(iff), zstophy,zout)
958	ENDIF
959
960	IF (flag_fsw_sol(iff)<=lev_files(iff)) THEN
961	call histdef(nid_files(iff), "fsw_"//clnsurf(nsrf), &
962	& "SW "//clnsurf(nsrf), "W/m2", &
963	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
964	& type_ecri(iff), zstophy,zout)
965	ENDIF
966
967	IF (flag_wbils_sol(iff)<=lev_files(iff)) THEN
968	call histdef(nid_files(iff), "wbils_"//clnsurf(nsrf), &
969	& "Bilan sol "//clnsurf(nsrf), "W/m2", &
970	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
971	& type_ecri(iff), zstophy,zout)
972	ENDIF
973
974	IF (flag_wbilo_sol(iff)<=lev_files(iff)) THEN
975	call histdef(nid_files(iff), "wbilo_"//clnsurf(nsrf), &
976	& "Bilan eau "//clnsurf(nsrf), "kg/(m2*s)", &
977	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
978	& type_ecri(iff), zstophy,zout)
979	ENDIF
980
981	if (iflag_pbl>1 .and. lev_files(iff).gt.10 ) then
982	IF (flag_tke_sol(iff)<=lev_files(iff)) THEN
983	call histdef(nid_files(iff), "tke_"//clnsurf(nsrf), &
984	& "Max Turb. Kinetic Energy "//clnsurf(nsrf), "-", &
985	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
986	& type_ecri(iff), zstophy,zout)
987	ENDIF
988
989	IF (flag_tke_max_sol(iff)<=lev_files(iff)) THEN
990	call histdef(nid_files(iff), "tke_max_"//clnsurf(nsrf), &
991	& "Max Turb. Kinetic Energy "//clnsurf(nsrf), "-", &
992	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
993	& "t_max(X)", zstophy,zout)
994	ENDIF
995	endif
996
997	END DO
998
999	IF (flag_cdrm(iff)<=lev_files(iff)) THEN
1000	CALL histdef(nid_files(iff), "cdrm", &
1001	& "Momentum drag coef.", "-", &
1002	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1003	& type_ecri(iff), zstophy,zout)
1004	ENDIF
1005
1006	IF (flag_cdrh(iff)<=lev_files(iff)) THEN
1007	CALL histdef(nid_files(iff), "cdrh", &
1008	& "Heat drag coef.", "-", &
1009	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1010	& type_ecri(iff), zstophy,zout)
1011	ENDIF
1012
1013	IF (flag_cldl(iff)<=lev_files(iff)) THEN
1014	CALL histdef(nid_files(iff), "cldl", &
1015	& "Low-level cloudiness", "-", &
1016	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1017	& type_ecri(iff), zstophy,zout)
1018	ENDIF
1019
1020	IF (flag_cldm(iff)<=lev_files(iff)) THEN
1021	CALL histdef(nid_files(iff), "cldm", &
1022	& "Mid-level cloudiness", "-", &
1023	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1024	& type_ecri(iff), zstophy,zout)
1025	ENDIF
1026
1027	IF (flag_cldh(iff)<=lev_files(iff)) THEN
1028	CALL histdef(nid_files(iff), "cldh", &
1029	& "High-level cloudiness", "-", &
1030	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1031	& type_ecri(iff), zstophy,zout)
1032	ENDIF
1033
1034	IF (flag_cldt(iff)<=lev_files(iff)) THEN
1035	CALL histdef(nid_files(iff), "cldt", &
1036	& "Total cloudiness", "%", &
1037	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1038	& type_ecri(iff), zstophy,zout)
1039	ENDIF
1040
1041	IF (flag_cldq(iff)<=lev_files(iff)) THEN
1042	CALL histdef(nid_files(iff),"cldq", &
1043	& "Cloud liquid water path","kg/m2", &
1044	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1045	& type_ecri(iff), zstophy,zout)
1046	ENDIF
1047
1048	IF (flag_lwp(iff)<=lev_files(iff)) THEN
1049	CALL histdef(nid_files(iff),"lwp", &
1050	& "Cloud water path","kg/m2", &
1051	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1052	& type_ecri(iff), zstophy,zout)
1053	ENDIF
1054
1055	IF (flag_iwp(iff)<=lev_files(iff)) THEN
1056	CALL histdef(nid_files(iff),"iwp", &
1057	& "Cloud ice water path","kg/m2", &
1058	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1059	& type_ecri(iff), zstophy,zout)
1060	ENDIF
1061
1062	IF (flag_ue(iff)<=lev_files(iff)) THEN
1063	CALL histdef(nid_files(iff), "ue", &
1064	& "Zonal energy transport", "-", &
1065	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1066	& type_ecri(iff), zstophy,zout)
1067	ENDIF
1068
1069	IF (flag_ve(iff)<=lev_files(iff)) THEN
1070	CALL histdef(nid_files(iff), "ve", &
1071	& "Merid energy transport", "-", &
1072	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1073	& type_ecri(iff), zstophy,zout)
1074	ENDIF
1075
1076	IF (flag_uq(iff)<=lev_files(iff)) THEN
1077	CALL histdef(nid_files(iff), "uq", &
1078	& "Zonal humidity transport", "-", &
1079	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1080	& type_ecri(iff), zstophy,zout)
1081	ENDIF
1082
1083	IF (flag_vq(iff)<=lev_files(iff)) THEN
1084	CALL histdef(nid_files(iff), "vq", &
1085	& "Merid humidity transport", "-", &
1086	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1087	& type_ecri(iff), zstophy,zout)
1088	ENDIF
1089
1090	IF(iflag_con.GE.3) THEN ! sb
1091	IF (flag_cape(iff)<=lev_files(iff)) THEN
1092	CALL histdef(nid_files(iff), "cape", &
1093	& "Conv avlbl pot ener", "J/kg", &
1094	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1095	& type_ecri(iff), zstophy,zout)
1096	ENDIF
1097
1098	IF (flag_pbase(iff)<=lev_files(iff)) THEN
1099	CALL histdef(nid_files(iff), "pbase", &
1100	& "Cld base pressure", "mb", &
1101	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1102	& type_ecri(iff), zstophy,zout)
1103	ENDIF
1104
1105	IF (flag_ptop(iff)<=lev_files(iff)) THEN
1106	CALL histdef(nid_files(iff), "ptop", &
1107	& "Cld top pressure", "mb", &
1108	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1109	& type_ecri(iff), zstophy,zout)
1110	ENDIF
1111
1112	IF (flag_fbase(iff)<=lev_files(iff)) THEN
1113	CALL histdef(nid_files(iff), "fbase", &
1114	& "Cld base mass flux", "kg/m2/s", &
1115	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1116	& type_ecri(iff), zstophy,zout)
1117	ENDIF
1118
1119	IF (flag_prw(iff)<=lev_files(iff)) THEN
1120	CALL histdef(nid_files(iff), "prw", &
1121	& "Precipitable water", "kg/m2", &
1122	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1123	& type_ecri(iff), zstophy,zout)
1124	ENDIF
1125
1126	ENDIF !iflag_con .GE. 3
1127
1128	IF (flag_s_pblh(iff)<=lev_files(iff)) THEN
1129	CALL histdef(nid_files(iff), "s_pblh", &
1130	& "Boundary Layer Height", "m", &
1131	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1132	& type_ecri(iff), zstophy,zout)
1133	ENDIF
1134
1135	IF (flag_s_pblt(iff)<=lev_files(iff)) THEN
1136	CALL histdef(nid_files(iff), "s_pblt", &
1137	& "t at Boundary Layer Height","K", &
1138	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1139	& type_ecri(iff), zstophy,zout)
1140	ENDIF
1141
1142	IF (flag_s_lclt(iff)<=lev_files(iff)) THEN
1143	CALL histdef(nid_files(iff), "s_lcl", &
1144	& "Condensation level", "m", &
1145	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1146	& type_ecri(iff), zstophy,zout)
1147	ENDIF
1148
1149	IF (flag_s_capCL(iff)<=lev_files(iff)) THEN
1150	CALL histdef(nid_files(iff), "s_capCL", &
1151	& "Conv avlbl pot enerfor ABL", "J/m2", &
1152	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1153	& type_ecri(iff), zstophy,zout)
1154	ENDIF
1155
1156	IF (flag_s_oliqCLL(iff)<=lev_files(iff)) THEN
1157	CALL histdef(nid_files(iff), "s_oliqCL", &
1158	& "Liq Water in BL", "kg/m2", &
1159	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1160	& type_ecri(iff), zstophy,zout)
1161	ENDIF
1162
1163	IF (flag_s_cteiCL(iff)<=lev_files(iff)) THEN
1164	CALL histdef(nid_files(iff), "s_cteiCL", &
1165	& "Instability criteria(ABL)", "K", &
1166	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1167	& type_ecri(iff), zstophy,zout)
1168	ENDIF
1169
1170	IF (flag_s_therm(iff)<=lev_files(iff)) THEN
1171	CALL histdef(nid_files(iff), "s_therm", &
1172	& "Exces du thermique", "K", &
1173	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1174	& type_ecri(iff), zstophy,zout)
1175	ENDIF
1176
1177	IF (flag_s_trmb1(iff)<=lev_files(iff)) THEN
1178	CALL histdef(nid_files(iff), "s_trmb1", &
1179	& "deep_cape(HBTM2)", "J/m2", &
1180	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1181	& type_ecri(iff), zstophy,zout)
1182	ENDIF
1183
1184	IF (flag_s_trmb2(iff)<=lev_files(iff)) THEN
1185	CALL histdef(nid_files(iff), "s_trmb2", &
1186	& "inhibition (HBTM2)", "J/m2", &
1187	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1188	& type_ecri(iff), zstophy,zout)
1189	ENDIF
1190
1191	IF (flag_s_trmb3(iff)<=lev_files(iff)) THEN
1192	CALL histdef(nid_files(iff), "s_trmb3", &
1193	& "Point Omega (HBTM2)", "m", &
1194	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1195	& type_ecri(iff), zstophy,zout)
1196	ENDIF
1197
1198	! Champs interpolles sur des niveaux de pression
1199	! iif=1 on ecrit u v w phi sur 850 700 500 200 au niv 1
1200	! iif=2 on ecrit w et ph 500 seulement au niv 1
1201	! et u v sur 850 700 500 200
1202	! iif=3 on ecrit ph a 500 seulement au niv 1
1203	! on ecrit u v t q a 850 700 500 200 au niv 3
1204
1205	DO k=1, nlevSTD
1206	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
1207	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
1208	IF(bb2.EQ."850".OR.bb2.EQ."700".OR. &
1209	& bb2.EQ."500".OR.bb2.EQ."200") THEN
1210
1211	IF (flag_ulevsSTD(iff)<=lev_files(iff)) THEN
1212	CALL histdef(nid_files(iff), "u"//bb2, &
1213	& "Zonal wind "//bb2//"mb","m/s", &
1214	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1215	& "inst(X)", zout,zout)
1216	ENDIF
1217
1218	IF (flag_vlevsSTD(iff)<=lev_files(iff)) THEN
1219	CALL histdef(nid_files(iff), "v"//bb2, &
1220	& "Meridional wind "//bb2//"mb","m/s", &
1221	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1222	& "inst(X)", zout,zout)
1223	ENDIF
1224
1225	IF (flag_wlevsSTD(iff)<=lev_files(iff)) THEN
1226	CALL histdef(nid_files(iff), "w"//bb2, &
1227	& "Vertical wind "//bb2//"mb","m/s", &
1228	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1229	& "inst(X)", zout,zout)
1230	ENDIF
1231
1232	IF (flag_philevsSTD(iff)<=lev_files(iff)) THEN
1233	CALL histdef(nid_files(iff), "phi"//bb2, &
1234	& "Geopotential "//bb2//"mb","m", &
1235	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1236	& "inst(X)", zout,zout)
1237	ENDIF
1238
1239	IF (flag_qlevsSTD(iff)<=lev_files(iff)) THEN
1240	CALL histdef(nid_files(iff), "q"//bb2, &
1241	& "Specific humidity "//bb2//"mb","kg/kg", &
1242	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1243	& "inst(X)", zout,zout)
1244	ENDIF
1245
1246	IF (flag_tlevsSTD(iff)<=lev_files(iff)) THEN
1247	CALL histdef(nid_files(iff), "t"//bb2, &
1248	& "Temperature "//bb2//"mb","K", &
1249	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1250	& "inst(X)", zout,zout)
1251	ENDIF
1252
1253	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
1254	ENDDO
1255
1256	!IM diagnostiques flux ocean-atm ou ocean-glace de mer
1257	!IM pour utilisation dans un modele de "slab" ocean
1258
1259	IF (flag_fluxo(iff)<=lev_files(iff)) THEN
1260	CALL histdef(nid_files(iff), "fluxo", &
1261	& "Flux turbulents ocean-atmosphere", "W/m2", &
1262	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1263	& type_ecri(iff), zstophy,zout)
1264	ENDIF
1265
1266	IF (flag_fluxg(iff)<=lev_files(iff)) THEN
1267	CALL histdef(nid_files(iff), "fluxg", &
1268	& "Flux turbulents ocean-glace de mer","W/m2", &
1269	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1270	& type_ecri(iff), zstophy,zout)
1271	ENDIF
1272
1273	IF (flag_t_oce_sic(iff)<=lev_files(iff)) THEN
1274	CALL histdef(nid_files(iff), "t_oce_sic", &
1275	& "Temp. mixte oce-sic","K", &
1276	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1277	& type_ecri(iff), zstophy,zout)
1278	ENDIF
1279
1280	IF (OCEAN.EQ.'force ') THEN
1281	IF (flag_lmt_bils(iff)<=lev_files(iff)) THEN
1282	CALL histdef(nid_files(iff), "lmt_bils", &
1283	& "Bilan au sol atmosphere forcee", "W/m2", &
1284	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1285	& type_ecri(iff), zstophy,zout)
1286	ENDIF
1287
1288	ELSE IF (OCEAN.EQ.'slab ') THEN
1289
1290	IF (flag_slab_bils(iff)<=lev_files(iff)) THEN
1291	CALL histdef(nid_files(iff), "slab_bils", &
1292	& "Bilan au sol Slab", "W/m2", &
1293	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1294	& type_ecri(iff), zstophy,zout)
1295	ENDIF
1296
1297	IF (flag_tslab(iff)<=lev_files(iff)) THEN
1298	CALL histdef(nid_files(iff), "tslab", "Slab SST ", "K", &
1299	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1300	& type_ecri(iff), zstophy,zout)
1301	ENDIF
1302
1303	IF (flag_seaice(iff)<=lev_files(iff)) THEN
1304	CALL histdef(nid_files(iff),"seaice","Slab seaice","kg/m2", &
1305	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1306	& type_ecri(iff), zstophy,zout)
1307	ENDIF
1308
1309	IF (flag_siceh(iff)<=lev_files(iff)) THEN
1310	CALL histdef(nid_files(iff),"siceh", &
1311	& "Slab seaice height","m", &
1312	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1313	& type_ecri(iff), zstophy,zout)
1314	ENDIF
1315	ENDIF
1316
1317	IF (flag_weakinv(iff)<=lev_files(iff)) THEN
1318	CALL histdef(nid_files(iff), "weakinv", &
1319	& "Weak inversion", "-", &
1320	& iim,jjmp1,nhori, 1,1,1, -99, 32, &
1321	& type_ecri(iff), zstophy,zout)
1322	ENDIF
1323
1324	IF (flag_dthmin(iff)<=lev_files(iff)) THEN
1325	CALL histdef(nid_files(iff), "dthmin", &
1326	& "dTheta mini", "K/m", &
1327	& iim,jjmp1,nhori, 1,1,1, -99, 32, &
1328	& type_ecri(iff), zstophy,zout)
1329	ENDIF
1330
1331	IF (flag_meantaucld(iff)<=lev_files(iff)) THEN
1332	CALL histdef(nid_files(iff),"meantaucld", &
1333	& "ISCCP mean cloud optical thickness","1", &
1334	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1335	& type_ecri(iff), zstophy,zout)
1336	ENDIF
1337
1338	IF (flag_cldtau(iff)<=lev_files(iff)) THEN
1339	CALL histdef(nid_files(iff),"cldtau", &
1340	& "Cloud optical thickness","1", &
1341	& iim,jj_nb,nhori, klev,1,klev, nvert, 32, &
1342	& type_ecri(iff), zstophy,zout)
1343	ENDIF
1344
1345	IF (flag_cldemi(iff)<=lev_files(iff)) THEN
1346	CALL histdef(nid_files(iff),"cldemi", &
1347	& "Cloud optical emissivity","1", &
1348	& iim,jj_nb,nhori, klev,1,klev, nvert, 32, &
1349	& type_ecri(iff), zstophy,zout)
1350	ENDIF
1351
1352	IF (flag_rh2m(iff)<=lev_files(iff)) THEN
1353	CALL histdef(nid_files(iff), "rh2m", &
1354	& "Relative humidity at 2m", "%", &
1355	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1356	& type_ecri(iff), zstophy,zout)
1357	ENDIF
1358
1359	IF (flag_qsat2m(iff)<=lev_files(iff)) THEN
1360	CALL histdef(nid_files(iff), "qsat2m", &
1361	& "Saturant humidity at 2m", "%", &
1362	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1363	& type_ecri(iff), zstophy,zout)
1364	ENDIF
1365
1366	IF (flag_tpot(iff)<=lev_files(iff)) THEN
1367	CALL histdef(nid_files(iff), "tpot", &
1368	& "Surface air potential temperature", "K", &
1369	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1370	& type_ecri(iff), zstophy,zout)
1371	ENDIF
1372
1373	IF (flag_tpote(iff)<=lev_files(iff)) THEN
1374	CALL histdef(nid_files(iff), "tpote", &
1375	& "Surface air equivalent potential temperature", "K", &
1376	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1377	& type_ecri(iff), zstophy,zout)
1378	ENDIF
1379
1380	IF (flag_SWnetOR(iff)<=lev_files(iff)) THEN
1381	CALL histdef(nid_files(iff), "SWnetOR", &
1382	& "Sfce net SW radiation OR", "W/m2", &
1383	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1384	& type_ecri(iff), zstophy,zout)
1385	ENDIF
1386
1387	IF (flag_SWdownOR(iff)<=lev_files(iff)) THEN
1388	CALL histdef(nid_files(iff), "SWdownOR", &
1389	& "Sfce incident SW radiation OR", "W/m2", &
1390	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1391	& type_ecri(iff), zstophy,zout)
1392	ENDIF
1393
1394	IF (flag_LWdownOR(iff)<=lev_files(iff)) THEN
1395	CALL histdef(nid_files(iff), "LWdownOR", &
1396	& "Sfce incident LW radiation OR", "W/m2", &
1397	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1398	& type_ecri(iff), zstophy,zout)
1399	ENDIF
1400
1401	IF (flag_snowl(iff)<=lev_files(iff)) THEN
1402	CALL histdef(nid_files(iff), "snowl", &
1403	& "Solid Large-scale Precip.","kg/(m2*s)", &
1404	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1405	& type_ecri(iff), zstophy,zout)
1406	ENDIF
1407
1408	IF (flag_cape_max(iff)<=lev_files(iff)) THEN
1409	CALL histdef(nid_files(iff), "cape_max", "CAPE max.", &
1410	& "J/kg", &
1411	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1412	& 't_max(X)', zstophy,zout)
1413	ENDIF
1414
1415	IF (flag_solldown(iff)<=lev_files(iff)) THEN
1416	CALL histdef(nid_files(iff), "solldown", &
1417	& "Down. IR rad. at surface", &
1418	& "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, &
1419	& type_ecri(iff), zstophy,zout)
1420	ENDIF
1421
1422	IF (flag_dtsvdfo(iff)<=lev_files(iff)) THEN
1423	CALL histdef(nid_files(iff), "dtsvdfo", &
1424	& "Boundary-layer dTs(o)", &
1425	& "K/s", iim,jj_nb,nhori, 1,1,1, -99, 32, &
1426	& type_ecri(iff), zout,zout)
1427	ENDIF
1428
1429	IF (flag_dtsvdft(iff)<=lev_files(iff)) THEN
1430	CALL histdef(nid_files(iff), "dtsvdft", &
1431	& "Boundary-layer dTs(t)", "K/s", &
1432	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1433	& type_ecri(iff), zout,zout)
1434	ENDIF
1435
1436	IF (flag_dtsvdfg(iff)<=lev_files(iff)) THEN
1437	CALL histdef(nid_files(iff), "dtsvdfg", &
1438	& "Boundary-layer dTs(g)", "K/s", &
1439	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1440	& type_ecri(iff), zout,zout)
1441	ENDIF
1442
1443	IF (flag_dtsvdfi(iff)<=lev_files(iff)) THEN
1444	CALL histdef(nid_files(iff), "dtsvdfi", &
1445	& "Boundary-layer dTs(g)", "K/s", &
1446	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1447	& type_ecri(iff), zout,zout)
1448	ENDIF
1449
1450	IF (flag_rugs(iff)<=lev_files(iff)) THEN
1451	CALL histdef(nid_files(iff), "rugs", "rugosity", "-", &
1452	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1453	& "inst(X)", zout,zout)
1454	ENDIF
1455
1456	! Champs 3D:
1457
1458	IF (flag_lwcon(iff)<=lev_files(iff)) THEN
1459	CALL histdef(nid_files(iff),"lwcon", &
1460	& "Cloud liquid water content","kg/kg", &
1461	& iim,jj_nb,nhori, klev,1,klev, nvert, 32, &
1462	& type_ecri(iff), zstophy,zout)
1463	ENDIF
1464
1465	IF (flag_iwcon(iff)<=lev_files(iff)) THEN
1466	CALL histdef(nid_files(iff),"iwcon", &
1467	& "Cloud ice water content","kg/kg", &
1468	& iim,jj_nb,nhori, klev,1,klev, nvert, 32, &
1469	& type_ecri(iff), zstophy,zout)
1470	ENDIF
1471
1472	IF (flag_temp(iff)<=lev_files(iff)) THEN
1473	CALL histdef(nid_files(iff),"temp","Air temperature","K", &
1474	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1475	& type_ecri(iff), zstophy,zout)
1476	ENDIF
1477
1478	IF (flag_ovap(iff)<=lev_files(iff)) THEN
1479	CALL histdef(nid_files(iff), "ovap", &
1480	& "Specific humidity","kg/kg", &
1481	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1482	& type_ecri(iff), zstophy,zout)
1483	ENDIF
1484
1485
1486	IF (flag_geop(iff)<=lev_files(iff)) THEN
1487	CALL histdef(nid_files(iff), "geop", &
1488	& "Geopotential height", "m2/s2", &
1489	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1490	& type_ecri(iff), zstophy,zout)
1491	ENDIF
1492
1493	IF (flag_vitu(iff)<=lev_files(iff)) THEN
1494	CALL histdef(nid_files(iff), "vitu", "Zonal wind", "m/s", &
1495	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1496	& type_ecri(iff), zstophy,zout)
1497	ENDIF
1498
1499	IF (flag_vitv(iff)<=lev_files(iff)) THEN
1500	CALL histdef(nid_files(iff),"vitv","Meridional wind","m/s", &
1501	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1502	& type_ecri(iff), zstophy,zout)
1503	ENDIF
1504
1505	IF (flag_vitw(iff)<=lev_files(iff)) THEN
1506	CALL histdef(nid_files(iff),"vitw","Vertical wind","m/s", &
1507	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1508	& type_ecri(iff), zstophy,zout)
1509	ENDIF
1510
1511	IF (flag_pres(iff)<=lev_files(iff)) THEN
1512	CALL histdef(nid_files(iff), "pres", "Air pressure", "Pa", &
1513	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1514	& type_ecri(iff), zstophy,zout)
1515	ENDIF
1516
1517	IF (flag_rneb(iff)<=lev_files(iff)) THEN
1518	CALL histdef(nid_files(iff), "rneb", "Cloud fraction", "-", &
1519	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1520	& type_ecri(iff), zstophy,zout)
1521	ENDIF
1522
1523	IF (flag_rnebcon(iff)<=lev_files(iff)) THEN
1524	CALL histdef(nid_files(iff), "rnebcon", &
1525	& "Convective Cloud Fraction", "-", &
1526	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1527	& type_ecri(iff), zstophy,zout)
1528	ENDIF
1529
1530	IF (flag_rhum(iff)<=lev_files(iff)) THEN
1531	CALL histdef(nid_files(iff), "rhum", "Relative humidity", "-", &
1532	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1533	& type_ecri(iff), zstophy,zout)
1534	ENDIF
1535
1536	IF (flag_ozone(iff)<=lev_files(iff)) THEN
1537	CALL histdef(nid_files(iff), "ozone", &
1538	& "Ozone concentration", "ppmv", &
1539	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1540	& type_ecri(iff), zstophy,zout)
1541	ENDIF
1542
1543	IF (flag_upwd(iff)<=lev_files(iff)) THEN
1544	CALL histdef(nid_files(iff), "upwd", &
1545	& "saturated updraft", "kg/m2/s", &
1546	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1547	& type_ecri(iff), zstophy,zout)
1548	ENDIF
1549
1550	IF (flag_dtphy(iff)<=lev_files(iff)) THEN
1551	CALL histdef(nid_files(iff), "dtphy", "Physics dT", "K/s", &
1552	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1553	& type_ecri(iff), zstophy,zout)
1554	ENDIF
1555
1556	IF (flag_dqphy(iff)<=lev_files(iff)) THEN
1557	CALL histdef(nid_files(iff), "dqphy", &
1558	& "Physics dQ", "(kg/kg)/s", &
1559	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1560	& type_ecri(iff), zstophy,zout)
1561	ENDIF
1562
1563	DO nsrf=1, nbsrf
1564	IF (flag_albe_sol(iff)<=lev_files(iff)) THEN
1565	call histdef(nid_files(iff), "albe_"//clnsurf(nsrf), &
1566	& "Albedo surf. "//clnsurf(nsrf), "-", &
1567	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1568	& type_ecri(iff), zstophy,zout)
1569	ENDIF
1570
1571	IF (flag_rugs_sol(iff)<=lev_files(iff)) THEN
1572	call histdef(nid_files(iff), "rugs_"//clnsurf(nsrf), &
1573	& "Latent heat flux "//clnsurf(nsrf), "W/m2", &
1574	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1575	& type_ecri(iff), zstophy,zout)
1576	ENDIF
1577
1578	IF (flag_ages_sol(iff)<=lev_files(iff)) THEN
1579	CALL histdef(nid_files(iff), &
1580	& "ages_"//clnsurf(nsrf), "Snow age","day", &
1581	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1582	& type_ecri(iff), zstophy,zout)
1583	ENDIF
1584
1585	ENDDO !nsrf=1, nbsrf
1586
1587	IF (flag_albs(iff)<=lev_files(iff)) THEN
1588	CALL histdef(nid_files(iff),"albs","Surface albedo","-", &
1589	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1590	& type_ecri(iff), zstophy,zout)
1591	ENDIF
1592
1593	IF (flag_albslw(iff)<=lev_files(iff)) THEN
1594	CALL histdef(nid_files(iff),"albslw", &
1595	& "Surface albedo LW", "-", &
1596	& iim,jj_nb,nhori, 1,1,1, -99, 32, &
1597	& type_ecri(iff), zstophy,zout)
1598	ENDIF
1599
1600	!FH Sorties pour la couche limite
1601	if (iflag_pbl>1) then
1602	IF (flag_tke(iff)<=lev_files(iff)) THEN
1603	CALL histdef(nid_files(iff), "tke","TKE","m2/s2", &
1604	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1605	& type_ecri(iff), zstophy,zout)
1606	ENDIF
1607
1608	IF (flag_tke_max(iff)<=lev_files(iff)) THEN
1609	CALL histdef(nid_files(iff), "tke_max","TKE max","m2/s2", &
1610	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1611	& "t_max(X)", zstophy,zout)
1612	ENDIF
1613	endif
1614
1615	IF (flag_kz(iff)<=lev_files(iff)) THEN
1616	CALL histdef(nid_files(iff), "kz","Kz melange","m2/s", &
1617	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1618	& type_ecri(iff), zstophy,zout)
1619	ENDIF
1620
1621	IF (flag_kz_max(iff)<=lev_files(iff)) THEN
1622	CALL histdef(nid_files(iff), "kz_max", &
1623	& "Kz melange max","m2/s", &
1624	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1625	& "t_max(X)", zstophy,zout)
1626	ENDIF
1627
1628	IF (flag_clwcon(iff)<=lev_files(iff)) THEN
1629	CALL histdef(nid_files(iff), "clwcon", &
1630	& "Convective Cloud Liquid water content", &
1631	& "kg/kg", &
1632	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1633	& type_ecri(iff), zstophy,zout)
1634	ENDIF
1635
1636	IF (flag_Ma(iff)<=lev_files(iff)) THEN
1637	CALL histdef(nid_files(iff),"Ma", &
1638	& "undilute adiab updraft","kg/m2/s", &
1639	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1640	& type_ecri(iff), zstophy,zout)
1641	ENDIF
1642
1643	IF (flag_dnwd(iff)<=lev_files(iff)) THEN
1644	CALL histdef(nid_files(iff), "dnwd", &
1645	& "saturated downdraft","kg/m2/s", &
1646	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1647	& type_ecri(iff), zstophy,zout)
1648	ENDIF
1649
1650	IF (flag_dnwd0(iff)<=lev_files(iff)) THEN
1651	CALL histdef(nid_files(iff), "dnwd0", &
1652	& "unsat. downdraft", "kg/m2/s", &
1653	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1654	& type_ecri(iff), zstophy,zout)
1655	ENDIF
1656
1657	IF (flag_dtdyn(iff)<=lev_files(iff)) THEN
1658	CALL histdef(nid_files(iff),"dtdyn","Dynamics dT","K/s", &
1659	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1660	& type_ecri(iff), zstophy,zout)
1661	ENDIF
1662
1663	IF (flag_dqdyn(iff)<=lev_files(iff)) THEN
1664	CALL histdef(nid_files(iff), "dqdyn", &
1665	& "Dynamics dQ", "(kg/kg)/s", &
1666	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1667	& type_ecri(iff), zstophy,zout)
1668	ENDIF
1669
1670	IF (flag_dtcon(iff)<=lev_files(iff)) THEN
1671	CALL histdef(nid_files(iff),"dtcon","Convection dT","K/s", &
1672	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1673	& type_ecri(iff), zstophy,zout)
1674	ENDIF
1675
1676	IF (flag_ducon(iff)<=lev_files(iff)) THEN
1677	CALL histdef(nid_files(iff),"ducon","Convection du","m/s2", &
1678	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1679	& type_ecri(iff), zstophy,zout)
1680	ENDIF
1681
1682	IF (flag_dqcon(iff)<=lev_files(iff)) THEN
1683	CALL histdef(nid_files(iff), "dqcon", &
1684	& "Convection dQ", "(kg/kg)/s", &
1685	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1686	& type_ecri(iff), zstophy,zout)
1687	ENDIF
1688
1689	IF (flag_dtlsc(iff)<=lev_files(iff)) THEN
1690	CALL histdef(nid_files(iff), "dtlsc", &
1691	& "Condensation dT", "K/s", &
1692	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1693	& type_ecri(iff), zstophy,zout)
1694	ENDIF
1695
1696	IF (flag_dtlschr(iff)<=lev_files(iff)) THEN
1697	CALL histdef(nid_files(iff), "dtlschr", &
1698	& "Large-scale condensational heating rate", &
1699	& "K/s",iim,jj_nb, nhori, klev,1,klev,nvert, &
1700	& 32,type_ecri(iff), zstophy,zout)
1701	ENDIF
1702
1703	IF (flag_dqlsc(iff)<=lev_files(iff)) THEN
1704	CALL histdef(nid_files(iff), "dqlsc", &
1705	& "Condensation dQ", "(kg/kg)/s", &
1706	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1707	& type_ecri(iff), zstophy,zout)
1708	ENDIF
1709
1710	IF (flag_dtvdf(iff)<=lev_files(iff)) THEN
1711	CALL histdef(nid_files(iff), "dtvdf", &
1712	& "Boundary-layer dT", "K/s", &
1713	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1714	& type_ecri(iff), zstophy,zout)
1715	ENDIF
1716
1717	IF (flag_dqvdf(iff)<=lev_files(iff)) THEN
1718	CALL histdef(nid_files(iff), "dqvdf", &
1719	& "Boundary-layer dQ","(kg/kg)/s", &
1720	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1721	& type_ecri(iff), zstophy,zout)
1722	ENDIF
1723
1724	IF (flag_dteva(iff)<=lev_files(iff)) THEN
1725	CALL histdef(nid_files(iff), "dteva", &
1726	& "Reevaporation dT", "K/s", &
1727	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1728	& type_ecri(iff), zstophy,zout)
1729	ENDIF
1730
1731	IF (flag_dqeva(iff)<=lev_files(iff)) THEN
1732	CALL histdef(nid_files(iff), "dqeva", &
1733	& "Reevaporation dQ","(kg/kg)/s", &
1734	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1735	& type_ecri(iff), zstophy,zout)
1736	ENDIF
1737
1738	IF (flag_ptconv(iff)<=lev_files(iff)) THEN
1739	CALL histdef(nid_files(iff), "ptconv", &
1740	& "POINTS CONVECTIFS"," ", &
1741	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1742	& type_ecri(iff), zstophy,zout)
1743	ENDIF
1744
1745	IF (flag_ratqs(iff)<=lev_files(iff)) THEN
1746	CALL histdef(nid_files(iff), "ratqs", "RATQS"," ", &
1747	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1748	& type_ecri(iff), zstophy,zout)
1749	ENDIF
1750
1751	IF (flag_dtthe(iff)<=lev_files(iff)) THEN
1752	CALL histdef(nid_files(iff), "dtthe", &
1753	& "Dry adjust. dT", "K/s", &
1754	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1755	& type_ecri(iff), zstophy,zout)
1756	ENDIF
1757
1758	IF (flag_dqthe(iff)<=lev_files(iff)) THEN
1759	CALL histdef(nid_files(iff),"dqthe", &
1760	& "Dry adjust. dQ","(kg/kg)/s", &
1761	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1762	& type_ecri(iff), zstophy,zout)
1763	ENDIF
1764
1765	IF (flag_dtajs(iff)<=lev_files(iff)) THEN
1766	CALL histdef(nid_files(iff), "dtajs", &
1767	& "Dry adjust. dT", "K/s", &
1768	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1769	& type_ecri(iff), zstophy,zout)
1770	ENDIF
1771
1772	IF (flag_dqajs(iff)<=lev_files(iff)) THEN
1773	CALL histdef(nid_files(iff),"dqajs", &
1774	& "Dry adjust. dQ","(kg/kg)/s", &
1775	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1776	& type_ecri(iff), zstophy,zout)
1777	ENDIF
1778
1779	IF (flag_dtswr(iff)<=lev_files(iff)) THEN
1780	CALL histdef(nid_files(iff), "dtswr", &
1781	& "SW radiation dT", "K/s", &
1782	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1783	& type_ecri(iff), zstophy,zout)
1784	ENDIF
1785
1786	IF (flag_dtsw0(iff)<=lev_files(iff)) THEN
1787	CALL histdef(nid_files(iff), "dtsw0", &
1788	& "CS SW radiation dT", "K/s", &
1789	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1790	& type_ecri(iff), zstophy,zout)
1791	ENDIF
1792
1793	IF (flag_dtlwr(iff)<=lev_files(iff)) THEN
1794	CALL histdef(nid_files(iff), "dtlwr", &
1795	& "LW radiation dT", "K/s", &
1796	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1797	& type_ecri(iff), zstophy,zout)
1798	ENDIF
1799
1800	IF (flag_dtlw0(iff)<=lev_files(iff)) THEN
1801	CALL histdef(nid_files(iff),"dtlw0", &
1802	& "CS LW radiation dT","K/s", &
1803	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1804	& type_ecri(iff), zstophy,zout)
1805	ENDIF
1806
1807	IF (flag_dtec(iff)<=lev_files(iff)) THEN
1808	CALL histdef(nid_files(iff), "dtec", &
1809	& "Cinetic dissip dT", "K/s", &
1810	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1811	& type_ecri(iff), zstophy,zout)
1812	ENDIF
1813
1814	IF (flag_duvdf(iff)<=lev_files(iff)) THEN
1815	CALL histdef(nid_files(iff), "duvdf", &
1816	& "Boundary-layer dU", "m/s2", &
1817	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1818	& type_ecri(iff), zstophy,zout)
1819	ENDIF
1820
1821	IF (flag_dvvdf(iff)<=lev_files(iff)) THEN
1822	CALL histdef(nid_files(iff), "dvvdf", &
1823	& "Boundary-layer dV", "m/s2", &
1824	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1825	& type_ecri(iff), zstophy,zout)
1826	ENDIF
1827
1828	IF (ok_orodr) THEN
1829	IF (flag_duoro(iff)<=lev_files(iff)) THEN
1830	CALL histdef(nid_files(iff), "duoro", &
1831	& "Orography dU", "m/s2", &
1832	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1833	& type_ecri(iff), zstophy,zout)
1834	ENDIF
1835
1836	IF (flag_dvoro(iff)<=lev_files(iff)) THEN
1837	CALL histdef(nid_files(iff), "dvoro", &
1838	& "Orography dV", "m/s2", &
1839	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1840	& type_ecri(iff), zstophy,zout)
1841	ENDIF
1842
1843	ENDIF
1844
1845	IF (ok_orolf) THEN
1846	IF (flag_dulif(iff)<=lev_files(iff)) THEN
1847	CALL histdef(nid_files(iff), "dulif", &
1848	& "Orography dU", "m/s2", &
1849	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1850	& type_ecri(iff), zstophy,zout)
1851	ENDIF
1852
1853	IF (flag_dvlif(iff)<=lev_files(iff)) THEN
1854	CALL histdef(nid_files(iff), "dvlif", &
1855	& "Orography dV", "m/s2", &
1856	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1857	& type_ecri(iff), zstophy,zout)
1858	ENDIF
1859	ENDIF
1860
1861	IF (flag_trac(iff)<=lev_files(iff)) THEN
1862	if (nqmax>=3) THEN
1863	DO iq=3,nqmax
1864	iiq=niadv(iq)
1865	CALL histdef(nid_files(iff),tnom(iq),ttext(iiq), "-", &
1866	& iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
1867	& type_ecri(iff), zstophy,zout)
1868	ENDDO
1869	endif
1870	ENDIF
1871
1872	CALL histend(nid_files(iff))
1873
1874	ndex2d = 0
1875	ndex3d = 0
1876
1877	ENDIF ! clef_files
1878
1879	ENDDO !
1880	end subroutine phys_output_open
1881
1882	END MODULE phys_output_mod
1883

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

Download in other formats:

Original Format