Context Navigation

phys_output_mod.F90 @ 1259

Last change on this file since 1259 was 1258, checked in by idelkadi, 15 years ago
-Rajout des champs des precipitations verticales lsc (liq et ice) et conv (liq et ice) -Appel dans physiq.F de l'interface avec le simulateur COSP sous la cles CPP_COSP avec une frequence=freq_COSP(s) -Declaration de freq_COSP dans clesphys.h et lecture dans physiq.def via conf_phys.F90 -Sorie dans newmicro.F des rayons effectifs des particules liq et ice des nuages necessaires pour cosp
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 70.4 KB

Line
1	!
2	! $Id: phys_output_mod.F90 1258 2009-11-09 15:24:36Z idelkadi $
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	private histdef2d, histdef3d, conf_physoutputs
18
19
20	integer, parameter :: nfiles = 5
21	logical, dimension(nfiles), save :: clef_files
22	integer, dimension(nfiles), save :: lev_files
23	integer, dimension(nfiles), save :: nid_files
24	!!$OMP THREADPRIVATE(clef_files, lev_files,nid_files)
25
26	integer, dimension(nfiles), private, save :: nhorim, nvertm
27	integer, dimension(nfiles), private, save :: nvertap, nvertbp, nvertAlt
28	! integer, dimension(nfiles), private, save :: nvertp0
29	real, dimension(nfiles), private, save :: zoutm
30	real, private, save :: zdtime
31	CHARACTER(len=20), dimension(nfiles), private, save :: type_ecri
32	!$OMP THREADPRIVATE(nhorim, nvertm, zoutm,zdtime,type_ecri)
33
34	! integer, save :: nid_hf3d
35
36	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
37	!! Definition pour chaque variable du niveau d ecriture dans chaque fichier
38	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!/ histmth, histday, histhf, histins /),'!!!!!!!!!!!!
39	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
40
41	integer, private:: levmin(nfiles) = 1
42	integer, private:: levmax(nfiles)
43
44	TYPE ctrl_out
45	integer,dimension(5) :: flag
46	character(len=20) :: name
47	END TYPE ctrl_out
48
49	!!! Comosentes de la coordonnee sigma-hybride
50	!!! Ap et Bp
51	type(ctrl_out),save :: o_Ahyb = ctrl_out((/ 1, 1, 1, 1, 1 /), 'Ap')
52	type(ctrl_out),save :: o_Bhyb = ctrl_out((/ 1, 1, 1, 1, 1 /), 'Bp')
53	type(ctrl_out),save :: o_Alt = ctrl_out((/ 1, 1, 1, 1, 1 /), 'Alt')
54
55	!!! 1D
56	type(ctrl_out),save :: o_phis = ctrl_out((/ 1, 1, 10, 1, 1 /), 'phis')
57	type(ctrl_out),save :: o_aire = ctrl_out((/ 1, 1, 10, 1, 1 /),'aire')
58	type(ctrl_out),save :: o_contfracATM = ctrl_out((/ 10, 1, 1, 10, 10 /),'contfracATM')
59	type(ctrl_out),save :: o_contfracOR = ctrl_out((/ 10, 1, 1, 10, 10 /),'contfracOR')
60	type(ctrl_out),save :: o_aireTER = ctrl_out((/ 10, 10, 1, 10, 10 /),'aireTER')
61
62	!!! 2D
63	type(ctrl_out),save :: o_flat = ctrl_out((/ 10, 1, 10, 10, 1 /),'flat')
64	type(ctrl_out),save :: o_slp = ctrl_out((/ 1, 1, 1, 10, 1 /),'slp')
65	type(ctrl_out),save :: o_tsol = ctrl_out((/ 1, 1, 1, 1, 1 /),'tsol')
66	type(ctrl_out),save :: o_t2m = ctrl_out((/ 1, 1, 1, 1, 1 /),'t2m')
67	type(ctrl_out),save :: o_t2m_min = ctrl_out((/ 1, 1, 10, 10, 10 /),'t2m_min')
68	type(ctrl_out),save :: o_t2m_max = ctrl_out((/ 1, 1, 10, 10, 10 /),'t2m_max')
69	type(ctrl_out),save,dimension(4) :: o_t2m_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_ter'), &
70	ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_lic'), &
71	ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_oce'), &
72	ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_sic') /)
73
74	type(ctrl_out),save :: o_wind10m = ctrl_out((/ 1, 1, 1, 10, 10 /),'wind10m')
75	type(ctrl_out),save :: o_wind10max = ctrl_out((/ 10, 1, 10, 10, 10 /),'wind10max')
76	type(ctrl_out),save :: o_sicf = ctrl_out((/ 1, 1, 10, 10, 10 /),'sicf')
77	type(ctrl_out),save :: o_q2m = ctrl_out((/ 1, 1, 1, 1, 1 /),'q2m')
78	type(ctrl_out),save :: o_u10m = ctrl_out((/ 1, 1, 1, 1, 1 /),'u10m')
79	type(ctrl_out),save :: o_v10m = ctrl_out((/ 1, 1, 1, 1, 1 /),'v10m')
80	type(ctrl_out),save :: o_psol = ctrl_out((/ 1, 1, 1, 1, 1 /),'psol')
81	type(ctrl_out),save :: o_qsurf = ctrl_out((/ 1, 10, 10, 10, 10 /),'qsurf')
82
83	type(ctrl_out),save,dimension(4) :: o_u10m_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_ter'), &
84	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_lic'), &
85	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_oce'), &
86	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_sic') /)
87
88	type(ctrl_out),save,dimension(4) :: o_v10m_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_ter'), &
89	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_lic'), &
90	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_oce'), &
91	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_sic') /)
92
93	type(ctrl_out),save :: o_qsol = ctrl_out((/ 1, 10, 10, 1, 1 /),'qsol')
94
95	type(ctrl_out),save :: o_ndayrain = ctrl_out((/ 1, 10, 10, 10, 10 /),'ndayrain')
96	type(ctrl_out),save :: o_precip = ctrl_out((/ 1, 1, 1, 1, 1 /),'precip')
97	type(ctrl_out),save :: o_plul = ctrl_out((/ 1, 1, 1, 1, 10 /),'plul')
98
99	type(ctrl_out),save :: o_pluc = ctrl_out((/ 1, 1, 1, 1, 10 /),'pluc')
100	type(ctrl_out),save :: o_snow = ctrl_out((/ 1, 1, 10, 1, 10 /),'snow')
101	type(ctrl_out),save :: o_evap = ctrl_out((/ 1, 1, 10, 1, 10 /),'evap')
102	type(ctrl_out),save :: o_tops = ctrl_out((/ 1, 1, 10, 10, 10 /),'tops')
103	type(ctrl_out),save :: o_tops0 = ctrl_out((/ 1, 5, 10, 10, 10 /),'tops0')
104	type(ctrl_out),save :: o_topl = ctrl_out((/ 1, 1, 10, 1, 10 /),'topl')
105	type(ctrl_out),save :: o_topl0 = ctrl_out((/ 1, 5, 10, 10, 10 /),'topl0')
106	type(ctrl_out),save :: o_SWupTOA = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWupTOA')
107	type(ctrl_out),save :: o_SWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWupTOAclr')
108	type(ctrl_out),save :: o_SWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWdnTOA')
109	type(ctrl_out),save :: o_SWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWdnTOAclr')
110	type(ctrl_out),save :: o_SWup200 = ctrl_out((/ 1, 10, 10, 10, 10 /),'SWup200')
111	type(ctrl_out),save :: o_SWup200clr = ctrl_out((/ 10, 1, 10, 10, 10 /),'SWup200clr')
112	type(ctrl_out),save :: o_SWdn200 = ctrl_out((/ 1, 10, 10, 10, 10 /),'SWdn200')
113	type(ctrl_out),save :: o_SWdn200clr = ctrl_out((/ 10, 1, 10, 10, 10 /),'SWdn200clr')
114
115	! arajouter
116	! type(ctrl_out),save :: o_LWupTOA = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWupTOA')
117	! type(ctrl_out),save :: o_LWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWupTOAclr')
118	! type(ctrl_out),save :: o_LWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWdnTOA')
119	! type(ctrl_out),save :: o_LWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWdnTOAclr')
120
121	type(ctrl_out),save :: o_LWup200 = ctrl_out((/ 1, 10, 10, 10, 10 /),'LWup200')
122	type(ctrl_out),save :: o_LWup200clr = ctrl_out((/ 1, 10, 10, 10, 10 /),'LWup200clr')
123	type(ctrl_out),save :: o_LWdn200 = ctrl_out((/ 1, 10, 10, 10, 10 /),'LWdn200')
124	type(ctrl_out),save :: o_LWdn200clr = ctrl_out((/ 1, 10, 10, 10, 10 /),'LWdn200clr')
125	type(ctrl_out),save :: o_sols = ctrl_out((/ 1, 1, 10, 1, 10 /),'sols')
126	type(ctrl_out),save :: o_sols0 = ctrl_out((/ 1, 5, 10, 10, 10 /),'sols0')
127	type(ctrl_out),save :: o_soll = ctrl_out((/ 1, 1, 10, 1, 10 /),'soll')
128	type(ctrl_out),save :: o_soll0 = ctrl_out((/ 1, 5, 10, 10, 10 /),'soll0')
129	type(ctrl_out),save :: o_radsol = ctrl_out((/ 1, 1, 10, 10, 10 /),'radsol')
130	type(ctrl_out),save :: o_SWupSFC = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWupSFC')
131	type(ctrl_out),save :: o_SWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWupSFCclr')
132	type(ctrl_out),save :: o_SWdnSFC = ctrl_out((/ 1, 1, 10, 10, 10 /),'SWdnSFC')
133	type(ctrl_out),save :: o_SWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWdnSFCclr')
134	type(ctrl_out),save :: o_LWupSFC = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWupSFC')
135	type(ctrl_out),save :: o_LWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWupSFCclr')
136	type(ctrl_out),save :: o_LWdnSFC = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWdnSFC')
137	type(ctrl_out),save :: o_LWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWdnSFCclr')
138	type(ctrl_out),save :: o_bils = ctrl_out((/ 1, 2, 10, 1, 10 /),'bils')
139	type(ctrl_out),save :: o_sens = ctrl_out((/ 1, 1, 10, 1, 1 /),'sens')
140	type(ctrl_out),save :: o_fder = ctrl_out((/ 1, 2, 10, 1, 10 /),'fder')
141	type(ctrl_out),save :: o_ffonte = ctrl_out((/ 1, 10, 10, 10, 10 /),'ffonte')
142	type(ctrl_out),save :: o_fqcalving = ctrl_out((/ 1, 10, 10, 10, 10 /),'fqcalving')
143	type(ctrl_out),save :: o_fqfonte = ctrl_out((/ 1, 10, 10, 10, 10 /),'fqfonte')
144
145	type(ctrl_out),save,dimension(4) :: o_taux_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'taux_ter'), &
146	ctrl_out((/ 1, 4, 10, 1, 10 /),'taux_lic'), &
147	ctrl_out((/ 1, 4, 10, 1, 10 /),'taux_oce'), &
148	ctrl_out((/ 1, 4, 10, 1, 10 /),'taux_sic') /)
149
150	type(ctrl_out),save,dimension(4) :: o_tauy_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'tauy_ter'), &
151	ctrl_out((/ 1, 4, 10, 1, 10 /),'tauy_lic'), &
152	ctrl_out((/ 1, 4, 10, 1, 10 /),'tauy_oce'), &
153	ctrl_out((/ 1, 4, 10, 1, 10 /),'tauy_sic') /)
154
155
156	type(ctrl_out),save,dimension(4) :: o_pourc_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'pourc_ter'), &
157	ctrl_out((/ 1, 4, 10, 1, 10 /),'pourc_lic'), &
158	ctrl_out((/ 1, 4, 10, 1, 10 /),'pourc_oce'), &
159	ctrl_out((/ 1, 4, 10, 1, 10 /),'pourc_sic') /)
160
161	type(ctrl_out),save,dimension(4) :: o_fract_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'fract_ter'), &
162	ctrl_out((/ 1, 4, 10, 1, 10 /),'fract_lic'), &
163	ctrl_out((/ 1, 4, 10, 1, 10 /),'fract_oce'), &
164	ctrl_out((/ 1, 4, 10, 1, 10 /),'fract_sic') /)
165
166	type(ctrl_out),save,dimension(4) :: o_tsol_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'tsol_ter'), &
167	ctrl_out((/ 1, 4, 10, 1, 10 /),'tsol_lic'), &
168	ctrl_out((/ 1, 4, 10, 1, 10 /),'tsol_oce'), &
169	ctrl_out((/ 1, 4, 10, 1, 10 /),'tsol_sic') /)
170
171	type(ctrl_out),save,dimension(4) :: o_sens_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_ter'), &
172	ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_lic'), &
173	ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_oce'), &
174	ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_sic') /)
175
176	type(ctrl_out),save,dimension(4) :: o_lat_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_ter'), &
177	ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_lic'), &
178	ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_oce'), &
179	ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_sic') /)
180
181	type(ctrl_out),save,dimension(4) :: o_flw_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_ter'), &
182	ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_lic'), &
183	ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_oce'), &
184	ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_sic') /)
185
186	type(ctrl_out),save,dimension(4) :: o_fsw_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'fsw_ter'), &
187	ctrl_out((/ 1, 10, 10, 10, 10 /),'fsw_lic'), &
188	ctrl_out((/ 1, 10, 10, 10, 10 /),'fsw_oce'), &
189	ctrl_out((/ 1, 10, 10, 10, 10 /),'fsw_sic') /)
190
191	type(ctrl_out),save,dimension(4) :: o_wbils_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'wbils_ter'), &
192	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbils_lic'), &
193	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbils_oce'), &
194	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbils_sic') /)
195
196	type(ctrl_out),save,dimension(4) :: o_wbilo_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'wbilo_ter'), &
197	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbilo_lic'), &
198	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbilo_oce'), &
199	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbilo_sic') /)
200
201
202	type(ctrl_out),save :: o_cdrm = ctrl_out((/ 1, 10, 10, 1, 10 /),'cdrm')
203	type(ctrl_out),save :: o_cdrh = ctrl_out((/ 1, 10, 10, 1, 10 /),'cdrh')
204	type(ctrl_out),save :: o_cldl = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldl')
205	type(ctrl_out),save :: o_cldm = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldm')
206	type(ctrl_out),save :: o_cldh = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldh')
207	type(ctrl_out),save :: o_cldt = ctrl_out((/ 1, 1, 2, 10, 10 /),'cldt')
208	type(ctrl_out),save :: o_cldq = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldq')
209	type(ctrl_out),save :: o_lwp = ctrl_out((/ 1, 5, 10, 10, 10 /),'lwp')
210	type(ctrl_out),save :: o_iwp = ctrl_out((/ 1, 5, 10, 10, 10 /),'iwp')
211	type(ctrl_out),save :: o_ue = ctrl_out((/ 1, 10, 10, 10, 10 /),'ue')
212	type(ctrl_out),save :: o_ve = ctrl_out((/ 1, 10, 10, 10, 10 /),'ve')
213	type(ctrl_out),save :: o_uq = ctrl_out((/ 1, 10, 10, 10, 10 /),'uq')
214	type(ctrl_out),save :: o_vq = ctrl_out((/ 1, 10, 10, 10, 10 /),'vq')
215
216	type(ctrl_out),save :: o_cape = ctrl_out((/ 1, 10, 10, 10, 10 /),'cape')
217	type(ctrl_out),save :: o_pbase = ctrl_out((/ 1, 10, 10, 10, 10 /),'pbase')
218	type(ctrl_out),save :: o_ptop = ctrl_out((/ 1, 4, 10, 10, 10 /),'ptop')
219	type(ctrl_out),save :: o_fbase = ctrl_out((/ 1, 10, 10, 10, 10 /),'fbase')
220	type(ctrl_out),save :: o_prw = ctrl_out((/ 1, 1, 10, 10, 10 /),'prw')
221
222	type(ctrl_out),save :: o_s_pblh = ctrl_out((/ 1, 10, 10, 1, 1 /),'s_pblh')
223	type(ctrl_out),save :: o_s_pblt = ctrl_out((/ 1, 10, 10, 1, 1 /),'s_pblt')
224	type(ctrl_out),save :: o_s_lcl = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_lcl')
225	type(ctrl_out),save :: o_s_capCL = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_capCL')
226	type(ctrl_out),save :: o_s_oliqCL = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_oliqCL')
227	type(ctrl_out),save :: o_s_cteiCL = ctrl_out((/ 1, 10, 10, 1, 1 /),'s_cteiCL')
228	type(ctrl_out),save :: o_s_therm = ctrl_out((/ 1, 10, 10, 1, 1 /),'s_therm')
229	type(ctrl_out),save :: o_s_trmb1 = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_trmb1')
230	type(ctrl_out),save :: o_s_trmb2 = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_trmb2')
231	type(ctrl_out),save :: o_s_trmb3 = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_trmb3')
232
233	type(ctrl_out),save :: o_slab_bils = ctrl_out((/ 1, 1, 10, 10, 10 /),'slab_bils_oce')
234
235	type(ctrl_out),save :: o_ale_bl = ctrl_out((/ 1, 1, 1, 1, 10 /),'ale_bl')
236	type(ctrl_out),save :: o_alp_bl = ctrl_out((/ 1, 1, 1, 1, 10 /),'alp_bl')
237	type(ctrl_out),save :: o_ale_wk = ctrl_out((/ 1, 1, 1, 1, 10 /),'ale_wk')
238	type(ctrl_out),save :: o_alp_wk = ctrl_out((/ 1, 1, 1, 1, 10 /),'alp_wk')
239
240	type(ctrl_out),save :: o_ale = ctrl_out((/ 1, 1, 1, 1, 10 /),'ale')
241	type(ctrl_out),save :: o_alp = ctrl_out((/ 1, 1, 1, 1, 10 /),'alp')
242	type(ctrl_out),save :: o_cin = ctrl_out((/ 1, 1, 1, 1, 10 /),'cin')
243	type(ctrl_out),save :: o_wape = ctrl_out((/ 1, 1, 1, 1, 10 /),'wape')
244
245
246	! Champs interpolles sur des niveaux de pression ??? a faire correctement
247
248	type(ctrl_out),save,dimension(6) :: o_uSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'u850'), &
249	ctrl_out((/ 1, 1, 3, 10, 10 /),'u700'), &
250	ctrl_out((/ 1, 1, 3, 10, 10 /),'u500'), &
251	ctrl_out((/ 1, 1, 3, 10, 10 /),'u200'), &
252	ctrl_out((/ 1, 1, 3, 10, 10 /),'u50'), &
253	ctrl_out((/ 1, 1, 3, 10, 10 /),'u10') /)
254
255
256	type(ctrl_out),save,dimension(6) :: o_vSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'v850'), &
257	ctrl_out((/ 1, 1, 3, 10, 10 /),'v700'), &
258	ctrl_out((/ 1, 1, 3, 10, 10 /),'v500'), &
259	ctrl_out((/ 1, 1, 3, 10, 10 /),'v200'), &
260	ctrl_out((/ 1, 1, 3, 10, 10 /),'v50'), &
261	ctrl_out((/ 1, 1, 3, 10, 10 /),'v10') /)
262
263	type(ctrl_out),save,dimension(6) :: o_wSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'w850'), &
264	ctrl_out((/ 1, 1, 3, 10, 10 /),'w700'), &
265	ctrl_out((/ 1, 1, 3, 10, 10 /),'w500'), &
266	ctrl_out((/ 1, 1, 3, 10, 10 /),'w200'), &
267	ctrl_out((/ 1, 1, 3, 10, 10 /),'w50'), &
268	ctrl_out((/ 1, 1, 3, 10, 10 /),'w10') /)
269
270	type(ctrl_out),save,dimension(6) :: o_tSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'t850'), &
271	ctrl_out((/ 1, 1, 3, 10, 10 /),'t700'), &
272	ctrl_out((/ 1, 1, 3, 10, 10 /),'t500'), &
273	ctrl_out((/ 1, 1, 3, 10, 10 /),'t200'), &
274	ctrl_out((/ 1, 1, 3, 10, 10 /),'t50'), &
275	ctrl_out((/ 1, 1, 3, 10, 10 /),'t10') /)
276
277	type(ctrl_out),save,dimension(6) :: o_qSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'q850'), &
278	ctrl_out((/ 1, 1, 3, 10, 10 /),'q700'), &
279	ctrl_out((/ 1, 1, 3, 10, 10 /),'q500'), &
280	ctrl_out((/ 1, 1, 3, 10, 10 /),'q200'), &
281	ctrl_out((/ 1, 1, 3, 10, 10 /),'q50'), &
282	ctrl_out((/ 1, 1, 3, 10, 10 /),'q10') /)
283
284	type(ctrl_out),save,dimension(6) :: o_phiSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'phi850'), &
285	ctrl_out((/ 1, 1, 3, 10, 10 /),'phi700'), &
286	ctrl_out((/ 1, 1, 3, 10, 10 /),'phi500'), &
287	ctrl_out((/ 1, 1, 3, 10, 10 /),'phi200'), &
288	ctrl_out((/ 1, 1, 3, 10, 10 /),'phi50'), &
289	ctrl_out((/ 1, 1, 3, 10, 10 /),'phi10') /)
290
291
292	type(ctrl_out),save :: o_t_oce_sic = ctrl_out((/ 1, 10, 10, 10, 10 /),'t_oce_sic')
293
294	type(ctrl_out),save :: o_weakinv = ctrl_out((/ 10, 1, 10, 10, 10 /),'weakinv')
295	type(ctrl_out),save :: o_dthmin = ctrl_out((/ 10, 1, 10, 10, 10 /),'dthmin')
296	type(ctrl_out),save,dimension(4) :: o_u10_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'u10_ter'), &
297	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10_lic'), &
298	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10_oce'), &
299	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10_sic') /)
300
301	type(ctrl_out),save,dimension(4) :: o_v10_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'v10_ter'), &
302	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10_lic'), &
303	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10_oce'), &
304	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10_sic') /)
305
306	type(ctrl_out),save :: o_cldtau = ctrl_out((/ 10, 5, 10, 10, 10 /),'cldtau')
307	type(ctrl_out),save :: o_cldemi = ctrl_out((/ 10, 5, 10, 10, 10 /),'cldemi')
308	type(ctrl_out),save :: o_rh2m = ctrl_out((/ 10, 5, 10, 10, 10 /),'rh2m')
309	type(ctrl_out),save :: o_qsat2m = ctrl_out((/ 10, 5, 10, 10, 10 /),'qsat2m')
310	type(ctrl_out),save :: o_tpot = ctrl_out((/ 10, 5, 10, 10, 10 /),'tpot')
311	type(ctrl_out),save :: o_tpote = ctrl_out((/ 10, 5, 10, 10, 10 /),'tpote')
312	type(ctrl_out),save :: o_tke = ctrl_out((/ 4, 10, 10, 10, 10 /),'tke ')
313	type(ctrl_out),save :: o_tke_max = ctrl_out((/ 4, 10, 10, 10, 10 /),'tke_max')
314
315	type(ctrl_out),save,dimension(4) :: o_tke_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_ter'), &
316	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_lic'), &
317	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_oce'), &
318	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_sic') /)
319
320	type(ctrl_out),save,dimension(4) :: o_tke_max_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_max_ter'), &
321	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_max_lic'), &
322	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_max_oce'), &
323	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_max_sic') /)
324
325	type(ctrl_out),save :: o_kz = ctrl_out((/ 4, 10, 10, 10, 10 /),'kz')
326	type(ctrl_out),save :: o_kz_max = ctrl_out((/ 4, 10, 10, 10, 10 /),'kz_max')
327	type(ctrl_out),save :: o_SWnetOR = ctrl_out((/ 10, 10, 2, 10, 10 /),'SWnetOR')
328	type(ctrl_out),save :: o_SWdownOR = ctrl_out((/ 10, 10, 2, 10, 10 /),'SWdownOR')
329	type(ctrl_out),save :: o_LWdownOR = ctrl_out((/ 10, 10, 2, 10, 10 /),'LWdownOR')
330
331	type(ctrl_out),save :: o_snowl = ctrl_out((/ 10, 1, 10, 10, 10 /),'snowl')
332	type(ctrl_out),save :: o_cape_max = ctrl_out((/ 10, 1, 10, 10, 10 /),'cape_max')
333	type(ctrl_out),save :: o_solldown = ctrl_out((/ 10, 1, 10, 1, 10 /),'solldown')
334
335	type(ctrl_out),save :: o_dtsvdfo = ctrl_out((/ 10, 10, 10, 1, 10 /),'dtsvdfo')
336	type(ctrl_out),save :: o_dtsvdft = ctrl_out((/ 10, 10, 10, 1, 10 /),'dtsvdft')
337	type(ctrl_out),save :: o_dtsvdfg = ctrl_out((/ 10, 10, 10, 1, 10 /),'dtsvdfg')
338	type(ctrl_out),save :: o_dtsvdfi = ctrl_out((/ 10, 10, 10, 1, 10 /),'dtsvdfi')
339	type(ctrl_out),save :: o_rugs = ctrl_out((/ 10, 10, 10, 1, 1 /),'rugs')
340
341	type(ctrl_out),save :: o_topswad = ctrl_out((/ 4, 10, 10, 10, 10 /),'topswad')
342	type(ctrl_out),save :: o_topswai = ctrl_out((/ 4, 10, 10, 10, 10 /),'topswai')
343	type(ctrl_out),save :: o_solswad = ctrl_out((/ 4, 10, 10, 10, 10 /),'solswad')
344	type(ctrl_out),save :: o_solswai = ctrl_out((/ 4, 10, 10, 10, 10 /),'solswai')
345
346	type(ctrl_out),save,dimension(10) :: o_tausumaero = (/ ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_ASBCM'), &
347	ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_ASPOMM'), &
348	ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_ASSO4M'), &
349	ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_CSSO4M'), &
350	ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_SSSSM'), &
351	ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_ASSSM'), &
352	ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_CSSSM'), &
353	ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_CIDUSTM'), &
354	ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_AIBCM'), &
355	ctrl_out((/ 4, 10, 10, 10, 10 /),'OD550_AIPOMM') /)
356
357	type(ctrl_out),save :: o_swtoaas_nat = ctrl_out((/ 4, 10, 10, 10, 10 /),'swtoaas_nat')
358	type(ctrl_out),save :: o_swsrfas_nat = ctrl_out((/ 4, 10, 10, 10, 10 /),'swsrfas_nat')
359	type(ctrl_out),save :: o_swtoacs_nat = ctrl_out((/ 4, 10, 10, 10, 10 /),'swtoacs_nat')
360	type(ctrl_out),save :: o_swsrfcs_nat = ctrl_out((/ 4, 10, 10, 10, 10 /),'swsrfcs_nat')
361
362	type(ctrl_out),save :: o_swtoaas_ant = ctrl_out((/ 4, 10, 10, 10, 10 /),'swtoaas_ant')
363	type(ctrl_out),save :: o_swsrfas_ant = ctrl_out((/ 4, 10, 10, 10, 10 /),'swsrfas_ant')
364	type(ctrl_out),save :: o_swtoacs_ant = ctrl_out((/ 4, 10, 10, 10, 10 /),'swtoacs_ant')
365	type(ctrl_out),save :: o_swsrfcs_ant = ctrl_out((/ 4, 10, 10, 10, 10 /),'swsrfcs_ant')
366
367	type(ctrl_out),save :: o_swtoacf_nat = ctrl_out((/ 4, 10, 10, 10, 10 /),'swtoacf_nat')
368	type(ctrl_out),save :: o_swsrfcf_nat = ctrl_out((/ 4, 10, 10, 10, 10 /),'swsrfcf_nat')
369	type(ctrl_out),save :: o_swtoacf_ant = ctrl_out((/ 4, 10, 10, 10, 10 /),'swtoacf_ant')
370	type(ctrl_out),save :: o_swsrfcf_ant = ctrl_out((/ 4, 10, 10, 10, 10 /),'swsrfcf_ant')
371	type(ctrl_out),save :: o_swtoacf_zero = ctrl_out((/ 4, 10, 10, 10, 10 /),'swtoacf_zero')
372	type(ctrl_out),save :: o_swsrfcf_zero = ctrl_out((/ 4, 10, 10, 10, 10 /),'swsrfcf_zero')
373
374
375	!!!!!!!!!!!!!!!!!!!!!! 3D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
376	type(ctrl_out),save :: o_lwcon = ctrl_out((/ 2, 5, 10, 10, 1 /),'lwcon')
377	type(ctrl_out),save :: o_iwcon = ctrl_out((/ 2, 5, 10, 10, 10 /),'iwcon')
378	type(ctrl_out),save :: o_temp = ctrl_out((/ 2, 3, 4, 1, 1 /),'temp')
379	type(ctrl_out),save :: o_theta = ctrl_out((/ 2, 3, 4, 1, 1 /),'theta')
380	type(ctrl_out),save :: o_ovap = ctrl_out((/ 2, 3, 4, 1, 1 /),'ovap')
381	type(ctrl_out),save :: o_ovapinit = ctrl_out((/ 2, 3, 4, 1, 1 /),'ovapinit')
382	type(ctrl_out),save :: o_wvapp = ctrl_out((/ 2, 10, 10, 10, 10 /),'wvapp')
383	type(ctrl_out),save :: o_geop = ctrl_out((/ 2, 3, 10, 1, 1 /),'geop')
384	type(ctrl_out),save :: o_vitu = ctrl_out((/ 2, 3, 4, 1, 1 /),'vitu')
385	type(ctrl_out),save :: o_vitv = ctrl_out((/ 2, 3, 4, 1, 1 /),'vitv')
386	type(ctrl_out),save :: o_vitw = ctrl_out((/ 2, 3, 10, 10, 1 /),'vitw')
387	type(ctrl_out),save :: o_pres = ctrl_out((/ 2, 3, 10, 1, 1 /),'pres')
388	type(ctrl_out),save :: o_rneb = ctrl_out((/ 2, 5, 10, 10, 1 /),'rneb')
389	type(ctrl_out),save :: o_rnebcon = ctrl_out((/ 2, 5, 10, 10, 1 /),'rnebcon')
390	type(ctrl_out),save :: o_rhum = ctrl_out((/ 2, 10, 10, 10, 10 /),'rhum')
391	type(ctrl_out),save :: o_ozone = ctrl_out((/ 2, 10, 10, 10, 10 /),'ozone')
392	type(ctrl_out),save :: o_upwd = ctrl_out((/ 2, 10, 10, 10, 10 /),'upwd')
393	type(ctrl_out),save :: o_dtphy = ctrl_out((/ 2, 10, 10, 10, 1 /),'dtphy')
394	type(ctrl_out),save :: o_dqphy = ctrl_out((/ 2, 10, 10, 10, 1 /),'dqphy')
395	type(ctrl_out),save :: o_pr_con_l = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_con_l')
396	type(ctrl_out),save :: o_pr_con_i = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_con_i')
397	type(ctrl_out),save :: o_pr_lsc_l = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_lsc_l')
398	type(ctrl_out),save :: o_pr_lsc_i = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_lsc_i')
399	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
400
401	type(ctrl_out),save,dimension(4) :: o_albe_srf = (/ ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_ter'), &
402	ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_lic'), &
403	ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_oce'), &
404	ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_sic') /)
405
406	type(ctrl_out),save,dimension(4) :: o_ages_srf = (/ ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_ter'), &
407	ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_lic'), &
408	ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_oce'), &
409	ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_sic') /)
410
411	type(ctrl_out),save,dimension(4) :: o_rugs_srf = (/ ctrl_out((/ 3, 4, 10, 1, 10 /),'rugs_ter'), &
412	ctrl_out((/ 3, 4, 10, 1, 10 /),'rugs_lic'), &
413	ctrl_out((/ 3, 4, 10, 1, 10 /),'rugs_oce'), &
414	ctrl_out((/ 3, 4, 10, 1, 10 /),'rugs_sic') /)
415
416	type(ctrl_out),save :: o_albs = ctrl_out((/ 3, 10, 10, 1, 10 /),'albs')
417	type(ctrl_out),save :: o_albslw = ctrl_out((/ 3, 10, 10, 1, 10 /),'albslw')
418
419	type(ctrl_out),save :: o_clwcon = ctrl_out((/ 4, 10, 10, 10, 10 /),'clwcon')
420	type(ctrl_out),save :: o_Ma = ctrl_out((/ 4, 10, 10, 10, 10 /),'Ma')
421	type(ctrl_out),save :: o_dnwd = ctrl_out((/ 4, 10, 10, 10, 10 /),'dnwd')
422	type(ctrl_out),save :: o_dnwd0 = ctrl_out((/ 4, 10, 10, 10, 10 /),'dnwd0')
423	type(ctrl_out),save :: o_dtdyn = ctrl_out((/ 4, 10, 10, 10, 1 /),'dtdyn')
424	type(ctrl_out),save :: o_dqdyn = ctrl_out((/ 4, 10, 10, 10, 1 /),'dqdyn')
425	type(ctrl_out),save :: o_dudyn = ctrl_out((/ 4, 10, 10, 10, 1 /),'dudyn') !AXC
426	type(ctrl_out),save :: o_dvdyn = ctrl_out((/ 4, 10, 10, 10, 1 /),'dvdyn') !AXC
427	type(ctrl_out),save :: o_dtcon = ctrl_out((/ 4, 5, 10, 10, 10 /),'dtcon')
428	type(ctrl_out),save :: o_ducon = ctrl_out((/ 4, 10, 10, 10, 10 /),'ducon')
429	type(ctrl_out),save :: o_dqcon = ctrl_out((/ 4, 5, 10, 10, 10 /),'dqcon')
430	type(ctrl_out),save :: o_dtwak = ctrl_out((/ 4, 5, 10, 10, 10 /),'dtwak')
431	type(ctrl_out),save :: o_dqwak = ctrl_out((/ 4, 5, 10, 10, 10 /),'dqwak')
432	type(ctrl_out),save :: o_wake_h = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_h')
433	type(ctrl_out),save :: o_wake_s = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_s')
434	type(ctrl_out),save :: o_wake_deltat = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_deltat')
435	type(ctrl_out),save :: o_wake_deltaq = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_deltaq')
436	type(ctrl_out),save :: o_wake_omg = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_omg')
437	type(ctrl_out),save :: o_Vprecip = ctrl_out((/ 10, 10, 10, 10, 10 /),'Vprecip')
438	type(ctrl_out),save :: o_ftd = ctrl_out((/ 4, 5, 10, 10, 10 /),'ftd')
439	type(ctrl_out),save :: o_fqd = ctrl_out((/ 4, 5, 10, 10, 10 /),'fqd')
440	type(ctrl_out),save :: o_dtlsc = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtlsc')
441	type(ctrl_out),save :: o_dtlschr = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtlschr')
442	type(ctrl_out),save :: o_dqlsc = ctrl_out((/ 4, 10, 10, 10, 10 /),'dqlsc')
443	type(ctrl_out),save :: o_dtvdf = ctrl_out((/ 4, 10, 10, 1, 10 /),'dtvdf')
444	type(ctrl_out),save :: o_dqvdf = ctrl_out((/ 4, 10, 10, 1, 10 /),'dqvdf')
445	type(ctrl_out),save :: o_dteva = ctrl_out((/ 4, 10, 10, 10, 10 /),'dteva')
446	type(ctrl_out),save :: o_dqeva = ctrl_out((/ 4, 10, 10, 10, 10 /),'dqeva')
447	type(ctrl_out),save :: o_ptconv = ctrl_out((/ 4, 10, 10, 10, 10 /),'ptconv')
448	type(ctrl_out),save :: o_ratqs = ctrl_out((/ 4, 10, 10, 10, 10 /),'ratqs')
449	type(ctrl_out),save :: o_dtthe = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtthe')
450	type(ctrl_out),save :: o_f_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'f_th')
451	type(ctrl_out),save :: o_e_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'e_th')
452	type(ctrl_out),save :: o_w_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'w_th')
453	type(ctrl_out),save :: o_lambda_th = ctrl_out((/ 10, 10, 10, 10, 10 /),'lambda_th')
454	type(ctrl_out),save :: o_q_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'q_th')
455	type(ctrl_out),save :: o_a_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'a_th')
456	type(ctrl_out),save :: o_d_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'d_th')
457	type(ctrl_out),save :: o_f0_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'f0_th')
458	type(ctrl_out),save :: o_zmax_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'zmax_th')
459	type(ctrl_out),save :: o_dqthe = ctrl_out((/ 4, 10, 10, 10, 1 /),'dqthe')
460	type(ctrl_out),save :: o_dtajs = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtajs')
461	type(ctrl_out),save :: o_dqajs = ctrl_out((/ 4, 10, 10, 10, 10 /),'dqajs')
462	type(ctrl_out),save :: o_dtswr = ctrl_out((/ 4, 10, 10, 10, 1 /),'dtswr')
463	type(ctrl_out),save :: o_dtsw0 = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtsw0')
464	type(ctrl_out),save :: o_dtlwr = ctrl_out((/ 4, 10, 10, 10, 1 /),'dtlwr')
465	type(ctrl_out),save :: o_dtlw0 = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtlw0')
466	type(ctrl_out),save :: o_dtec = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtec')
467	type(ctrl_out),save :: o_duvdf = ctrl_out((/ 4, 10, 10, 10, 10 /),'duvdf')
468	type(ctrl_out),save :: o_dvvdf = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvvdf')
469	type(ctrl_out),save :: o_duoro = ctrl_out((/ 4, 10, 10, 10, 10 /),'duoro')
470	type(ctrl_out),save :: o_dvoro = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvoro')
471	type(ctrl_out),save :: o_dulif = ctrl_out((/ 4, 10, 10, 10, 10 /),'dulif')
472	type(ctrl_out),save :: o_dvlif = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvlif')
473
474	! Attention a refaire correctement
475	type(ctrl_out),save,dimension(2) :: o_trac = (/ ctrl_out((/ 4, 10, 10, 10, 10 /),'trac01'), &
476	ctrl_out((/ 4, 10, 10, 10, 10 /),'trac02') /)
477	CONTAINS
478
479	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
480	!!!!!!!!! Ouverture des fichier et definition des variable de sortie !!!!!!!!
481	!! histbeg, histvert et histdef
482	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
483
484	SUBROUTINE phys_output_open(jjmp1,nlevSTD,clevSTD,nbteta, &
485	ctetaSTD,dtime, ok_veget, &
486	type_ocean, iflag_pbl,ok_mensuel,ok_journe, &
487	ok_hf,ok_instan,ok_LES,ok_ade,ok_aie)
488
489	USE iophy
490	USE dimphy
491	USE infotrac
492	USE ioipsl
493	USE mod_phys_lmdz_para
494	USE aero_mod, only : naero_spc,name_aero
495
496	IMPLICIT NONE
497	include "dimensions.h"
498	include "temps.h"
499	include "indicesol.h"
500	include "clesphys.h"
501	include "thermcell.h"
502	include "comvert.h"
503
504	integer :: jjmp1
505	integer :: nbteta, nlevSTD, radpas
506	logical :: ok_mensuel, ok_journe, ok_hf, ok_instan
507	logical :: ok_LES,ok_ade,ok_aie
508	real :: dtime
509	integer :: idayref
510	real :: zjulian
511	real, dimension(klev) :: Ahyb, Bhyb, Alt
512	character(len=4), dimension(nlevSTD) :: clevSTD
513	integer :: nsrf, k, iq, iiq, iff, i, j, ilev
514	integer :: naero
515	logical :: ok_veget
516	integer :: iflag_pbl
517	CHARACTER(len=4) :: bb2
518	CHARACTER(len=2) :: bb3
519	character(len=6) :: type_ocean
520	CHARACTER(len=3) :: ctetaSTD(nbteta)
521	real, dimension(nfiles) :: ecrit_files
522	CHARACTER(len=20), dimension(nfiles) :: phys_out_filenames
523	INTEGER, dimension(iim*jjmp1) :: ndex2d
524	INTEGER, dimension(iimjjmp1klev) :: ndex3d
525	integer :: imin_ins, imax_ins
526	integer :: jmin_ins, jmax_ins
527	integer, dimension(nfiles) :: phys_out_levmin, phys_out_levmax
528	integer, dimension(nfiles) :: phys_out_filelevels
529	CHARACTER(len=20), dimension(nfiles) :: type_ecri_files, phys_out_filetypes
530	character(len=20), dimension(nfiles) :: chtimestep = (/ 'DefFreq', 'DefFreq','DefFreq', 'DefFreq', 'DefFreq' /)
531	logical, dimension(nfiles) :: phys_out_filekeys
532
533	!!!!!!!!!! stockage dans une region limitee pour chaque fichier !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
534	! entre [phys_out_lonmin,phys_out_lonmax] et [phys_out_latmin,phys_out_latmax]
535
536	logical, dimension(nfiles), save :: phys_out_regfkey = (/ .false., .false., .false., .false., .false. /)
537	real, dimension(nfiles), save :: phys_out_lonmin = (/ -180., -180., -180., -180., -180. /)
538	real, dimension(nfiles), save :: phys_out_lonmax = (/ 180., 180., 180., 180., 180. /)
539	real, dimension(nfiles), save :: phys_out_latmin = (/ -90., -90., -90., -90., -90. /)
540	real, dimension(nfiles), save :: phys_out_latmax = (/ 90., 90., 90., 90., 90. /)
541
542
543
544	!
545	print*,'Debut phys_output_mod.F90'
546	! Initialisations (Valeurs par defaut
547	levmax = (/ klev, klev, klev, klev, klev /)
548
549	phys_out_filenames(1) = 'histmth'
550	phys_out_filenames(2) = 'histday'
551	phys_out_filenames(3) = 'histhf'
552	phys_out_filenames(4) = 'histins'
553	phys_out_filenames(5) = 'histLES'
554
555	type_ecri(1) = 'ave(X)'
556	type_ecri(2) = 'ave(X)'
557	type_ecri(3) = 'ave(X)'
558	type_ecri(4) = 'inst(X)'
559	type_ecri(5) = 'inst(X)'
560
561	clef_files(1) = ok_mensuel
562	clef_files(2) = ok_journe
563	clef_files(3) = ok_hf
564	clef_files(4) = ok_instan
565	clef_files(5) = ok_LES
566
567	lev_files(1) = lev_histmth
568	lev_files(2) = lev_histday
569	lev_files(3) = lev_histhf
570	lev_files(4) = lev_histins
571	lev_files(5) = lev_histLES
572
573
574	ecrit_files(1) = ecrit_mth
575	ecrit_files(2) = ecrit_day
576	ecrit_files(3) = ecrit_hf
577	ecrit_files(4) = ecrit_ins
578	ecrit_files(5) = ecrit_LES
579
580	!! Lectures des parametres de sorties dans physiq.def
581
582	call getin('phys_out_regfkey',phys_out_regfkey)
583	call getin('phys_out_lonmin',phys_out_lonmin)
584	call getin('phys_out_lonmax',phys_out_lonmax)
585	call getin('phys_out_latmin',phys_out_latmin)
586	call getin('phys_out_latmax',phys_out_latmax)
587	phys_out_levmin(:)=levmin(:)
588	call getin('phys_out_levmin',levmin)
589	phys_out_levmax(:)=levmax(:)
590	call getin('phys_out_levmax',levmax)
591	call getin('phys_out_filenames',phys_out_filenames)
592	phys_out_filekeys(:)=clef_files(:)
593	call getin('phys_out_filekeys',clef_files)
594	phys_out_filelevels(:)=lev_files(:)
595	call getin('phys_out_filelevels',lev_files)
596	call getin('phys_out_filetimesteps',chtimestep)
597	phys_out_filetypes(:)=type_ecri(:)
598	call getin('phys_out_filetypes',type_ecri)
599
600	type_ecri_files(:)=type_ecri(:)
601
602	print*,'phys_out_lonmin=',phys_out_lonmin
603	print*,'phys_out_lonmax=',phys_out_lonmax
604	print*,'phys_out_latmin=',phys_out_latmin
605	print*,'phys_out_latmax=',phys_out_latmax
606	print*,'phys_out_filenames=',phys_out_filenames
607	print*,'phys_out_filetypes=',type_ecri
608	print*,'phys_out_filekeys=',clef_files
609	print*,'phys_out_filelevels=',lev_files
610
611	!!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
612	! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !!
613	! Appel des histbeg pour definir les variables (nom, moy ou inst, freq de sortie ..
614	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
615
616	zdtime = dtime ! Frequence ou l on moyenne
617
618	! Calcul des Ahyb, Bhyb et Alt
619	do k=1,klev
620	Ahyb(k)=(ap(k)+ap(k+1))/2.
621	Bhyb(k)=(bp(k)+bp(k+1))/2.
622	Alt(k)=log(preff/presnivs(k))*8.
623	enddo
624	! if(prt_level.ge.1) then
625	print*,'Ap Hybrid = ',Ahyb(1:klev)
626	print*,'Bp Hybrid = ',Bhyb(1:klev)
627	print*,'Alt approx des couches pour une haut d echelle de 8km = ',Alt(1:klev)
628	! endif
629	DO iff=1,nfiles
630
631	IF (clef_files(iff)) THEN
632
633	if ( chtimestep(iff).eq.'DefFreq' ) then
634	! Par defaut ecrit_files = (ecrit_mensuel ecrit_jour ecrit_hf ...)*86400.
635	ecrit_files(iff)=ecrit_files(iff)*86400.
636	else
637	call convers_timesteps(chtimestep(iff),ecrit_files(iff))
638	endif
639	print*,'ecrit_files(',iff,')= ',ecrit_files(iff)
640
641	zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde
642
643	idayref = day_ref
644	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
645
646	!!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
647	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
648	if (phys_out_regfkey(iff)) then
649
650	imin_ins=1
651	imax_ins=iim
652	jmin_ins=1
653	jmax_ins=jjmp1
654
655	! correction abderr
656	do i=1,iim
657	print*,'io_lon(i)=',io_lon(i)
658	if (io_lon(i).le.phys_out_lonmin(iff)) imin_ins=i
659	if (io_lon(i).le.phys_out_lonmax(iff)) imax_ins=i+1
660	enddo
661
662	do j=1,jjmp1
663	print*,'io_lat(j)=',io_lat(j)
664	if (io_lat(j).ge.phys_out_latmin(iff)) jmax_ins=j+1
665	if (io_lat(j).ge.phys_out_latmax(iff)) jmin_ins=j
666	enddo
667
668	print*,'On stoke le fichier histoire numero ',iff,' sur ', &
669	imin_ins,imax_ins,jmin_ins,jmax_ins
670	print*,'longitudes : ', &
671	io_lon(imin_ins),io_lon(imax_ins), &
672	'latitudes : ', &
673	io_lat(jmax_ins),io_lat(jmin_ins)
674
675	CALL histbeg(phys_out_filenames(iff),iim,io_lon,jjmp1,io_lat, &
676	imin_ins,imax_ins-imin_ins+1, &
677	jmin_ins,jmax_ins-jmin_ins+1, &
678	itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
679	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
680	else
681	CALL histbeg_phy(phys_out_filenames(iff),itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
682	endif
683
684	CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "mb", &
685	levmax(iff) - levmin(iff) + 1, &
686	presnivs(levmin(iff):levmax(iff))/100., nvertm(iff))
687
688	!!!!!!!!!!!!! Traitement des champs 3D pour histhf !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
689	!!!!!!!!!!!!!!! A Revoir plus tard !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
690	! IF (iff.eq.3.and.lev_files(iff).ge.4) THEN
691	! CALL histbeg_phy("histhf3d",itau_phy, &
692	! & zjulian, dtime, &
693	! & nhorim, nid_hf3d)
694
695	! CALL histvert(nid_hf3d, "presnivs", &
696	! & "Vertical levels", "mb", &
697	! & klev, presnivs/100., nvertm)
698	! ENDIF
699	!
700	!!!! Composentes de la coordonnee sigma-hybride
701	CALL histvert(nid_files(iff), "Ahyb","Ahyb comp of Hyb Cord ", "Pa", &
702	levmax(iff) - levmin(iff) + 1,Ahyb,nvertap(iff))
703
704	CALL histvert(nid_files(iff), "Bhyb","Bhyb comp of Hyb Cord", " ", &
705	levmax(iff) - levmin(iff) + 1,Bhyb,nvertbp(iff))
706
707	CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", &
708	levmax(iff) - levmin(iff) + 1,Alt,nvertAlt(iff))
709
710	! CALL histvert(nid_files(iff), "preff","Reference pressure", "Pa", &
711	! 1,preff,nvertp0(iff))
712	!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
713	CALL histdef2d(iff,o_phis%flag,o_phis%name,"Surface geop.height", "m2/s2")
714	type_ecri(1) = 'once'
715	type_ecri(2) = 'once'
716	type_ecri(3) = 'once'
717	type_ecri(4) = 'once'
718	type_ecri(5) = 'once'
719	CALL histdef2d(iff,o_aire%flag,o_aire%name,"Grid area", "-")
720	CALL histdef2d(iff,o_contfracATM%flag,o_contfracATM%name,"% sfce ter+lic", "-")
721	type_ecri(:) = type_ecri_files(:)
722
723	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
724	CALL histdef2d(iff,o_contfracOR%flag,o_contfracOR%name,"% sfce terre OR", "-" )
725	CALL histdef2d(iff,o_aireTER%flag,o_aireTER%name,"Grid area CONT", "-" )
726	CALL histdef2d(iff,o_flat%flag,o_flat%name, "Latent heat flux", "W/m2")
727	CALL histdef2d(iff,o_slp%flag,o_slp%name, "Sea Level Pressure", "Pa" )
728	CALL histdef2d(iff,o_tsol%flag,o_tsol%name, "Surface Temperature", "K")
729	CALL histdef2d(iff,o_t2m%flag,o_t2m%name, "Temperature 2m", "K" )
730	type_ecri(1) = 't_min(X)'
731	type_ecri(2) = 't_min(X)'
732	type_ecri(3) = 't_min(X)'
733	type_ecri(4) = 't_min(X)'
734	type_ecri(5) = 't_min(X)'
735	CALL histdef2d(iff,o_t2m_min%flag,o_t2m_min%name, "Temp 2m min", "K" )
736	type_ecri(1) = 't_max(X)'
737	type_ecri(2) = 't_max(X)'
738	type_ecri(3) = 't_max(X)'
739	type_ecri(4) = 't_max(X)'
740	type_ecri(5) = 't_max(X)'
741	CALL histdef2d(iff,o_t2m_max%flag,o_t2m_max%name, "Temp 2m max", "K" )
742	type_ecri(:) = type_ecri_files(:)
743	CALL histdef2d(iff,o_wind10m%flag,o_wind10m%name, "10-m wind speed", "m/s")
744	CALL histdef2d(iff,o_wind10max%flag,o_wind10max%name, "10m wind speed max", "m/s")
745	CALL histdef2d(iff,o_sicf%flag,o_sicf%name, "Sea-ice fraction", "-" )
746	CALL histdef2d(iff,o_q2m%flag,o_q2m%name, "Specific humidity 2m", "kg/kg")
747	CALL histdef2d(iff,o_u10m%flag,o_u10m%name, "Vent zonal 10m", "m/s" )
748	CALL histdef2d(iff,o_v10m%flag,o_v10m%name, "Vent meridien 10m", "m/s")
749	CALL histdef2d(iff,o_psol%flag,o_psol%name, "Surface Pressure", "Pa" )
750	CALL histdef2d(iff,o_qsurf%flag,o_qsurf%name, "Surface Air humidity", "kg/kg")
751
752	if (.not. ok_veget) then
753	CALL histdef2d(iff,o_qsol%flag,o_qsol%name, "Soil watter content", "mm" )
754	endif
755
756	CALL histdef2d(iff,o_ndayrain%flag,o_ndayrain%name, "Number of dayrain(liq+sol)", "-")
757	CALL histdef2d(iff,o_precip%flag,o_precip%name, "Precip Totale liq+sol", "kg/(s*m2)" )
758	CALL histdef2d(iff,o_plul%flag,o_plul%name, "Large-scale Precip.", "kg/(s*m2)")
759	CALL histdef2d(iff,o_pluc%flag,o_pluc%name, "Convective Precip.", "kg/(s*m2)")
760	CALL histdef2d(iff,o_snow%flag,o_snow%name, "Snow fall", "kg/(s*m2)" )
761	CALL histdef2d(iff,o_evap%flag,o_evap%name, "Evaporat", "kg/(s*m2)" )
762	CALL histdef2d(iff,o_tops%flag,o_tops%name, "Solar rad. at TOA", "W/m2")
763	CALL histdef2d(iff,o_tops0%flag,o_tops0%name, "CS Solar rad. at TOA", "W/m2")
764	CALL histdef2d(iff,o_topl%flag,o_topl%name, "IR rad. at TOA", "W/m2" )
765	CALL histdef2d(iff,o_topl0%flag,o_topl0%name, "IR rad. at TOA", "W/m2")
766	CALL histdef2d(iff,o_SWupTOA%flag,o_SWupTOA%name, "SWup at TOA", "W/m2")
767	CALL histdef2d(iff,o_SWupTOAclr%flag,o_SWupTOAclr%name, "SWup clear sky at TOA", "W/m2")
768	CALL histdef2d(iff,o_SWdnTOA%flag,o_SWdnTOA%name, "SWdn at TOA", "W/m2" )
769	CALL histdef2d(iff,o_SWdnTOAclr%flag,o_SWdnTOAclr%name, "SWdn clear sky at TOA", "W/m2")
770	CALL histdef2d(iff,o_SWup200%flag,o_SWup200%name, "SWup at 200mb", "W/m2" )
771	CALL histdef2d(iff,o_SWup200clr%flag,o_SWup200clr%name, "SWup clear sky at 200mb", "W/m2")
772	CALL histdef2d(iff,o_SWdn200%flag,o_SWdn200%name, "SWdn at 200mb", "W/m2" )
773	CALL histdef2d(iff,o_SWdn200clr%flag,o_SWdn200clr%name, "SWdn clear sky at 200mb", "W/m2")
774	CALL histdef2d(iff,o_LWup200%flag,o_LWup200%name, "LWup at 200mb", "W/m2")
775	CALL histdef2d(iff,o_LWup200clr%flag,o_LWup200clr%name, "LWup clear sky at 200mb", "W/m2")
776	CALL histdef2d(iff,o_LWdn200%flag,o_LWdn200%name, "LWdn at 200mb", "W/m2")
777	CALL histdef2d(iff,o_LWdn200clr%flag,o_LWdn200clr%name, "LWdn clear sky at 200mb", "W/m2")
778	CALL histdef2d(iff,o_sols%flag,o_sols%name, "Solar rad. at surf.", "W/m2")
779	CALL histdef2d(iff,o_sols0%flag,o_sols0%name, "Solar rad. at surf.", "W/m2")
780	CALL histdef2d(iff,o_soll%flag,o_soll%name, "IR rad. at surface", "W/m2")
781	CALL histdef2d(iff,o_radsol%flag,o_radsol%name, "Rayonnement au sol", "W/m2")
782	CALL histdef2d(iff,o_soll0%flag,o_soll0%name, "IR rad. at surface", "W/m2")
783	CALL histdef2d(iff,o_SWupSFC%flag,o_SWupSFC%name, "SWup at surface", "W/m2")
784	CALL histdef2d(iff,o_SWupSFCclr%flag,o_SWupSFCclr%name, "SWup clear sky at surface", "W/m2")
785	CALL histdef2d(iff,o_SWdnSFC%flag,o_SWdnSFC%name, "SWdn at surface", "W/m2")
786	CALL histdef2d(iff,o_SWdnSFCclr%flag,o_SWdnSFCclr%name, "SWdn clear sky at surface", "W/m2")
787	CALL histdef2d(iff,o_LWupSFC%flag,o_LWupSFC%name, "Upwd. IR rad. at surface", "W/m2")
788	CALL histdef2d(iff,o_LWdnSFC%flag,o_LWdnSFC%name, "Down. IR rad. at surface", "W/m2")
789	CALL histdef2d(iff,o_LWupSFCclr%flag,o_LWupSFCclr%name, "CS Upwd. IR rad. at surface", "W/m2")
790	CALL histdef2d(iff,o_LWdnSFCclr%flag,o_LWdnSFCclr%name, "Down. CS IR rad. at surface", "W/m2")
791	CALL histdef2d(iff,o_bils%flag,o_bils%name, "Surf. total heat flux", "W/m2")
792	CALL histdef2d(iff,o_sens%flag,o_sens%name, "Sensible heat flux", "W/m2")
793	CALL histdef2d(iff,o_fder%flag,o_fder%name, "Heat flux derivation", "W/m2")
794	CALL histdef2d(iff,o_ffonte%flag,o_ffonte%name, "Thermal flux for snow melting", "W/m2")
795	CALL histdef2d(iff,o_fqcalving%flag,o_fqcalving%name, "Ice Calving", "kg/m2/s")
796	CALL histdef2d(iff,o_fqfonte%flag,o_fqfonte%name, "Land ice melt", "kg/m2/s")
797
798	DO nsrf = 1, nbsrf
799	CALL histdef2d(iff,o_pourc_srf(nsrf)%flag,o_pourc_srf(nsrf)%name,"% "//clnsurf(nsrf),"%")
800	CALL histdef2d(iff,o_fract_srf(nsrf)%flag,o_fract_srf(nsrf)%name,"Fraction "//clnsurf(nsrf),"1")
801	CALL histdef2d(iff,o_taux_srf(nsrf)%flag,o_taux_srf(nsrf)%name,"Zonal wind stress"//clnsurf(nsrf),"Pa")
802	CALL histdef2d(iff,o_tauy_srf(nsrf)%flag,o_tauy_srf(nsrf)%name,"Meridional wind stress "//clnsurf(nsrf),"Pa")
803	CALL histdef2d(iff,o_tsol_srf(nsrf)%flag,o_tsol_srf(nsrf)%name,"Temperature "//clnsurf(nsrf),"K")
804	CALL histdef2d(iff,o_u10m_srf(nsrf)%flag,o_u10m_srf(nsrf)%name,"Vent Zonal 10m "//clnsurf(nsrf),"m/s")
805	CALL histdef2d(iff,o_v10m_srf(nsrf)%flag,o_v10m_srf(nsrf)%name,"Vent meredien 10m "//clnsurf(nsrf),"m/s")
806	CALL histdef2d(iff,o_t2m_srf(nsrf)%flag,o_t2m_srf(nsrf)%name,"Temp 2m "//clnsurf(nsrf),"K")
807	CALL histdef2d(iff,o_sens_srf(nsrf)%flag,o_sens_srf(nsrf)%name,"Sensible heat flux "//clnsurf(nsrf),"W/m2")
808	CALL histdef2d(iff,o_lat_srf(nsrf)%flag,o_lat_srf(nsrf)%name,"Latent heat flux "//clnsurf(nsrf),"W/m2")
809	CALL histdef2d(iff,o_flw_srf(nsrf)%flag,o_flw_srf(nsrf)%name,"LW "//clnsurf(nsrf),"W/m2")
810	CALL histdef2d(iff,o_fsw_srf(nsrf)%flag,o_fsw_srf(nsrf)%name,"SW "//clnsurf(nsrf),"W/m2")
811	CALL histdef2d(iff,o_wbils_srf(nsrf)%flag,o_wbils_srf(nsrf)%name,"Bilan sol "//clnsurf(nsrf),"W/m2" )
812	CALL histdef2d(iff,o_wbilo_srf(nsrf)%flag,o_wbilo_srf(nsrf)%name,"Bilan eau "//clnsurf(nsrf),"kg/(m2*s)")
813	if (iflag_pbl>1 .and. lev_files(iff).gt.10 ) then
814	CALL histdef2d(iff,o_tke_srf(nsrf)%flag,o_tke_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-")
815	type_ecri(1) = 't_max(X)'
816	type_ecri(2) = 't_max(X)'
817	type_ecri(3) = 't_max(X)'
818	type_ecri(4) = 't_max(X)'
819	type_ecri(5) = 't_max(X)'
820	CALL histdef2d(iff,o_tke_max_srf(nsrf)%flag,o_tke_max_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-")
821	type_ecri(:) = type_ecri_files(:)
822	endif
823	CALL histdef2d(iff,o_albe_srf(nsrf)%flag,o_albe_srf(nsrf)%name,"Albedo surf. "//clnsurf(nsrf),"-")
824	CALL histdef2d(iff,o_rugs_srf(nsrf)%flag,o_rugs_srf(nsrf)%name,"Latent heat flux "//clnsurf(nsrf),"W/m2")
825	CALL histdef2d(iff,o_ages_srf(nsrf)%flag,o_ages_srf(nsrf)%name,"Snow age", "day")
826	END DO
827
828	DO naero = 1, naero_spc
829	CALL histdef2d(iff,o_tausumaero(naero)%flag,o_tausumaero(naero)%name,"Aerosol Optical depth at 550 nm "//name_aero(naero),"1")
830	END DO
831
832	IF (ok_ade) THEN
833	CALL histdef2d(iff,o_topswad%flag,o_topswad%name, "ADE at TOA", "W/m2")
834	CALL histdef2d(iff,o_solswad%flag,o_solswad%name, "ADE at SRF", "W/m2")
835
836	CALL histdef2d(iff,o_swtoaas_nat%flag,o_swtoaas_nat%name, "Natural aerosol radiative forcing all-sky at TOA", "W/m2")
837	CALL histdef2d(iff,o_swsrfas_nat%flag,o_swsrfas_nat%name, "Natural aerosol radiative forcing all-sky at SRF", "W/m2")
838	CALL histdef2d(iff,o_swtoacs_nat%flag,o_swtoacs_nat%name, "Natural aerosol radiative forcing clear-sky at TOA", "W/m2")
839	CALL histdef2d(iff,o_swsrfcs_nat%flag,o_swsrfcs_nat%name, "Natural aerosol radiative forcing clear-sky at SRF", "W/m2")
840
841	CALL histdef2d(iff,o_swtoaas_ant%flag,o_swtoaas_ant%name, "Anthropogenic aerosol radiative forcing all-sky at TOA", "W/m2")
842	CALL histdef2d(iff,o_swsrfas_ant%flag,o_swsrfas_ant%name, "Anthropogenic aerosol radiative forcing all-sky at SRF", "W/m2")
843	CALL histdef2d(iff,o_swtoacs_ant%flag,o_swtoacs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at TOA", "W/m2")
844	CALL histdef2d(iff,o_swsrfcs_ant%flag,o_swsrfcs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at SRF", "W/m2")
845
846	CALL histdef2d(iff,o_swtoacf_nat%flag,o_swtoacf_nat%name, "Natural aerosol impact on cloud radiative forcing at TOA", "W/m2")
847	CALL histdef2d(iff,o_swsrfcf_nat%flag,o_swsrfcf_nat%name, "Natural aerosol impact on cloud radiative forcing at SRF", "W/m2")
848	CALL histdef2d(iff,o_swtoacf_ant%flag,o_swtoacf_ant%name, "Anthropogenic aerosol impact on cloud radiative forcing at TOA", "W/m2")
849	CALL histdef2d(iff,o_swsrfcf_ant%flag,o_swsrfcf_ant%name, "Anthropogenic aerosol impact on cloud radiative forcing at SRF", "W/m2")
850	CALL histdef2d(iff,o_swtoacf_zero%flag,o_swtoacf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at TOA", "W/m2")
851	CALL histdef2d(iff,o_swsrfcf_zero%flag,o_swsrfcf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at SRF", "W/m2")
852
853	ENDIF
854
855	IF (ok_aie) THEN
856	CALL histdef2d(iff,o_topswai%flag,o_topswai%name, "AIE at TOA", "W/m2")
857	CALL histdef2d(iff,o_solswai%flag,o_solswai%name, "AIE at SFR", "W/m2")
858	ENDIF
859
860
861	CALL histdef2d(iff,o_albs%flag,o_albs%name, "Surface albedo", "-")
862	CALL histdef2d(iff,o_albslw%flag,o_albslw%name, "Surface albedo LW", "-")
863	CALL histdef2d(iff,o_cdrm%flag,o_cdrm%name, "Momentum drag coef.", "-")
864	CALL histdef2d(iff,o_cdrh%flag,o_cdrh%name, "Heat drag coef.", "-" )
865	CALL histdef2d(iff,o_cldl%flag,o_cldl%name, "Low-level cloudiness", "-")
866	CALL histdef2d(iff,o_cldm%flag,o_cldm%name, "Mid-level cloudiness", "-")
867	CALL histdef2d(iff,o_cldh%flag,o_cldh%name, "High-level cloudiness", "-")
868	CALL histdef2d(iff,o_cldt%flag,o_cldt%name, "Total cloudiness", "%")
869	CALL histdef2d(iff,o_cldq%flag,o_cldq%name, "Cloud liquid water path", "kg/m2")
870	CALL histdef2d(iff,o_lwp%flag,o_lwp%name, "Cloud water path", "kg/m2")
871	CALL histdef2d(iff,o_iwp%flag,o_iwp%name, "Cloud ice water path", "kg/m2" )
872	CALL histdef2d(iff,o_ue%flag,o_ue%name, "Zonal energy transport", "-")
873	CALL histdef2d(iff,o_ve%flag,o_ve%name, "Merid energy transport", "-")
874	CALL histdef2d(iff,o_uq%flag,o_uq%name, "Zonal humidity transport", "-")
875	CALL histdef2d(iff,o_vq%flag,o_vq%name, "Merid humidity transport", "-")
876
877	IF(iflag_con.GE.3) THEN ! sb
878	CALL histdef2d(iff,o_cape%flag,o_cape%name, "Conv avlbl pot ener", "J/kg")
879	CALL histdef2d(iff,o_pbase%flag,o_pbase%name, "Cld base pressure", "mb")
880	CALL histdef2d(iff,o_ptop%flag,o_ptop%name, "Cld top pressure", "mb")
881	CALL histdef2d(iff,o_fbase%flag,o_fbase%name, "Cld base mass flux", "kg/m2/s")
882	CALL histdef2d(iff,o_prw%flag,o_prw%name, "Precipitable water", "kg/m2")
883	type_ecri(1) = 't_max(X)'
884	type_ecri(2) = 't_max(X)'
885	type_ecri(3) = 't_max(X)'
886	type_ecri(4) = 't_max(X)'
887	type_ecri(5) = 't_max(X)'
888	CALL histdef2d(iff,o_cape_max%flag,o_cape_max%name, "CAPE max.", "J/kg")
889	type_ecri(:) = type_ecri_files(:)
890	CALL histdef3d(iff,o_upwd%flag,o_upwd%name, "saturated updraft", "kg/m2/s")
891	CALL histdef3d(iff,o_Ma%flag,o_Ma%name, "undilute adiab updraft", "kg/m2/s")
892	CALL histdef3d(iff,o_dnwd%flag,o_dnwd%name, "saturated downdraft", "kg/m2/s")
893	CALL histdef3d(iff,o_dnwd0%flag,o_dnwd0%name, "unsat. downdraft", "kg/m2/s")
894	ENDIF !iflag_con .GE. 3
895
896	CALL histdef2d(iff,o_s_pblh%flag,o_s_pblh%name, "Boundary Layer Height", "m")
897	CALL histdef2d(iff,o_s_pblt%flag,o_s_pblt%name, "t at Boundary Layer Height", "K")
898	CALL histdef2d(iff,o_s_lcl%flag,o_s_lcl%name, "Condensation level", "m")
899	CALL histdef2d(iff,o_s_capCL%flag,o_s_capCL%name, "Conv avlbl pot enerfor ABL", "J/m2" )
900	CALL histdef2d(iff,o_s_oliqCL%flag,o_s_oliqCL%name, "Liq Water in BL", "kg/m2")
901	CALL histdef2d(iff,o_s_cteiCL%flag,o_s_cteiCL%name, "Instability criteria(ABL)", "K")
902	CALL histdef2d(iff,o_s_therm%flag,o_s_therm%name, "Exces du thermique", "K")
903	CALL histdef2d(iff,o_s_trmb1%flag,o_s_trmb1%name, "deep_cape(HBTM2)", "J/m2")
904	CALL histdef2d(iff,o_s_trmb2%flag,o_s_trmb2%name, "inhibition (HBTM2)", "J/m2")
905	CALL histdef2d(iff,o_s_trmb3%flag,o_s_trmb3%name, "Point Omega (HBTM2)", "m")
906
907	! Champs interpolles sur des niveaux de pression
908
909	type_ecri(1) = 'inst(X)'
910	type_ecri(2) = 'inst(X)'
911	type_ecri(3) = 'inst(X)'
912	type_ecri(4) = 'inst(X)'
913	type_ecri(5) = 'inst(X)'
914
915	! Attention a reverifier
916
917	ilev=0
918	DO k=1, nlevSTD
919	! IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
920	bb2=clevSTD(k)
921	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200".OR.bb2.EQ."50".OR.bb2.EQ."10")THEN
922	ilev=ilev+1
923	print*,'ilev k bb2 flag name ',ilev,k, bb2,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name
924	CALL histdef2d(iff,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name,"Zonal wind "//bb2//"mb", "m/s")
925	CALL histdef2d(iff,o_vSTDlevs(ilev)%flag,o_vSTDlevs(ilev)%name,"Meridional wind "//bb2//"mb", "m/s")
926	CALL histdef2d(iff,o_wSTDlevs(ilev)%flag,o_wSTDlevs(ilev)%name,"Vertical wind "//bb2//"mb", "Pa/s")
927	CALL histdef2d(iff,o_phiSTDlevs(ilev)%flag,o_phiSTDlevs(ilev)%name,"Geopotential "//bb2//"mb", "m")
928	CALL histdef2d(iff,o_qSTDlevs(ilev)%flag,o_qSTDlevs(ilev)%name,"Specific humidity "//bb2//"mb", "kg/kg" )
929	CALL histdef2d(iff,o_tSTDlevs(ilev)%flag,o_tSTDlevs(ilev)%name,"Temperature "//bb2//"mb", "K")
930	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR."500".OR.bb2.EQ."200".OR.bb2.EQ."50".OR.bb2.EQ."10")
931	ENDDO
932	type_ecri(:) = type_ecri_files(:)
933
934	CALL histdef2d(iff,o_t_oce_sic%flag,o_t_oce_sic%name, "Temp mixte oce-sic", "K")
935
936	IF (type_ocean=='slab') &
937	CALL histdef2d(iff,o_slab_bils%flag, o_slab_bils%name,"Bilan au sol sur ocean slab", "W/m2")
938
939	! Couplage conv-CL
940	IF (iflag_con.GE.3) THEN
941	IF (iflag_coupl.EQ.1) THEN
942	CALL histdef2d(iff,o_ale_bl%flag,o_ale_bl%name, "ALE BL", "m2/s2")
943	CALL histdef2d(iff,o_alp_bl%flag,o_alp_bl%name, "ALP BL", "m2/s2")
944	ENDIF
945	ENDIF !(iflag_con.GE.3)
946
947	CALL histdef2d(iff,o_weakinv%flag,o_weakinv%name, "Weak inversion", "-")
948	CALL histdef2d(iff,o_dthmin%flag,o_dthmin%name, "dTheta mini", "K/m")
949	CALL histdef2d(iff,o_rh2m%flag,o_rh2m%name, "Relative humidity at 2m", "%" )
950	CALL histdef2d(iff,o_qsat2m%flag,o_qsat2m%name, "Saturant humidity at 2m", "%")
951	CALL histdef2d(iff,o_tpot%flag,o_tpot%name, "Surface air potential temperature", "K")
952	CALL histdef2d(iff,o_tpote%flag,o_tpote%name, "Surface air equivalent potential temperature", "K")
953	CALL histdef2d(iff,o_SWnetOR%flag,o_SWnetOR%name, "Sfce net SW radiation OR", "W/m2")
954	CALL histdef2d(iff,o_SWdownOR%flag,o_SWdownOR%name, "Sfce incident SW radiation OR", "W/m2")
955	CALL histdef2d(iff,o_LWdownOR%flag,o_LWdownOR%name, "Sfce incident LW radiation OR", "W/m2")
956	CALL histdef2d(iff,o_snowl%flag,o_snowl%name, "Solid Large-scale Precip.", "kg/(m2*s)")
957
958	CALL histdef2d(iff,o_solldown%flag,o_solldown%name, "Down. IR rad. at surface", "W/m2")
959	CALL histdef2d(iff,o_dtsvdfo%flag,o_dtsvdfo%name, "Boundary-layer dTs(o)", "K/s")
960	CALL histdef2d(iff,o_dtsvdft%flag,o_dtsvdft%name, "Boundary-layer dTs(t)", "K/s")
961	CALL histdef2d(iff,o_dtsvdfg%flag,o_dtsvdfg%name, "Boundary-layer dTs(g)", "K/s")
962	CALL histdef2d(iff,o_dtsvdfi%flag,o_dtsvdfi%name, "Boundary-layer dTs(g)", "K/s")
963	CALL histdef2d(iff,o_rugs%flag,o_rugs%name, "rugosity", "-" )
964
965	! Champs 3D:
966	CALL histdef3d(iff,o_lwcon%flag,o_lwcon%name, "Cloud liquid water content", "kg/kg")
967	CALL histdef3d(iff,o_iwcon%flag,o_iwcon%name, "Cloud ice water content", "kg/kg")
968	CALL histdef3d(iff,o_temp%flag,o_temp%name, "Air temperature", "K" )
969	CALL histdef3d(iff,o_theta%flag,o_theta%name, "Potential air temperature", "K" )
970	CALL histdef3d(iff,o_ovap%flag,o_ovap%name, "Specific humidity + dqphy", "kg/kg" )
971	CALL histdef3d(iff,o_ovapinit%flag,o_ovapinit%name, "Specific humidity", "kg/kg" )
972	CALL histdef3d(iff,o_geop%flag,o_geop%name, "Geopotential height", "m2/s2")
973	CALL histdef3d(iff,o_vitu%flag,o_vitu%name, "Zonal wind", "m/s" )
974	CALL histdef3d(iff,o_vitv%flag,o_vitv%name, "Meridional wind", "m/s" )
975	CALL histdef3d(iff,o_vitw%flag,o_vitw%name, "Vertical wind", "Pa/s" )
976	CALL histdef3d(iff,o_pres%flag,o_pres%name, "Air pressure", "Pa" )
977	CALL histdef3d(iff,o_rneb%flag,o_rneb%name, "Cloud fraction", "-")
978	CALL histdef3d(iff,o_rnebcon%flag,o_rnebcon%name, "Convective Cloud Fraction", "-")
979	CALL histdef3d(iff,o_rhum%flag,o_rhum%name, "Relative humidity", "-")
980	CALL histdef3d(iff,o_ozone%flag,o_ozone%name, "Ozone mole fraction", "-")
981	CALL histdef3d(iff,o_dtphy%flag,o_dtphy%name, "Physics dT", "K/s")
982	CALL histdef3d(iff,o_dqphy%flag,o_dqphy%name, "Physics dQ", "(kg/kg)/s")
983	CALL histdef3d(iff,o_cldtau%flag,o_cldtau%name, "Cloud optical thickness", "1")
984	CALL histdef3d(iff,o_cldemi%flag,o_cldemi%name, "Cloud optical emissivity", "1")
985	!IM: bug ?? dimensionnement variables (klon,klev+1) pmflxr, pmflxs, prfl, psfl
986	CALL histdef3d(iff,o_pr_con_l%flag,o_pr_con_l%name, "Convective precipitation lic", " ")
987	CALL histdef3d(iff,o_pr_con_i%flag,o_pr_con_i%name, "Convective precipitation ice", " ")
988	CALL histdef3d(iff,o_pr_lsc_l%flag,o_pr_lsc_l%name, "Large scale precipitation lic", " ")
989	CALL histdef3d(iff,o_pr_lsc_i%flag,o_pr_lsc_i%name, "Large scale precipitation ice", " ")
990	!FH Sorties pour la couche limite
991	if (iflag_pbl>1) then
992	CALL histdef3d(iff,o_tke%flag,o_tke%name, "TKE", "m2/s2")
993	type_ecri(1) = 't_max(X)'
994	type_ecri(2) = 't_max(X)'
995	type_ecri(3) = 't_max(X)'
996	type_ecri(4) = 't_max(X)'
997	type_ecri(5) = 't_max(X)'
998	CALL histdef3d(iff,o_tke_max%flag,o_tke_max%name, "TKE max", "m2/s2")
999	type_ecri(:) = type_ecri_files(:)
1000	endif
1001
1002	CALL histdef3d(iff,o_kz%flag,o_kz%name, "Kz melange", "m2/s")
1003	type_ecri(1) = 't_max(X)'
1004	type_ecri(2) = 't_max(X)'
1005	type_ecri(3) = 't_max(X)'
1006	type_ecri(4) = 't_max(X)'
1007	type_ecri(5) = 't_max(X)'
1008	CALL histdef3d(iff,o_kz_max%flag,o_kz_max%name, "Kz melange max", "m2/s" )
1009	type_ecri(:) = type_ecri_files(:)
1010	CALL histdef3d(iff,o_clwcon%flag,o_clwcon%name, "Convective Cloud Liquid water content", "kg/kg")
1011	CALL histdef3d(iff,o_dtdyn%flag,o_dtdyn%name, "Dynamics dT", "K/s")
1012	CALL histdef3d(iff,o_dqdyn%flag,o_dqdyn%name, "Dynamics dQ", "(kg/kg)/s")
1013	CALL histdef3d(iff,o_dudyn%flag,o_dudyn%name, "Dynamics dU", "m/s2")
1014	CALL histdef3d(iff,o_dvdyn%flag,o_dvdyn%name, "Dynamics dV", "m/s2")
1015	CALL histdef3d(iff,o_dtcon%flag,o_dtcon%name, "Convection dT", "K/s")
1016	CALL histdef3d(iff,o_ducon%flag,o_ducon%name, "Convection du", "m/s2")
1017	CALL histdef3d(iff,o_dqcon%flag,o_dqcon%name, "Convection dQ", "(kg/kg)/s")
1018
1019	! Wakes
1020	IF(iflag_con.EQ.3) THEN
1021	IF (iflag_wake == 1) THEN
1022	CALL histdef2d(iff,o_ale_wk%flag,o_ale_wk%name, "ALE WK", "m2/s2")
1023	CALL histdef2d(iff,o_alp_wk%flag,o_alp_wk%name, "ALP WK", "m2/s2")
1024	CALL histdef2d(iff,o_ale%flag,o_ale%name, "ALE", "m2/s2")
1025	CALL histdef2d(iff,o_alp%flag,o_alp%name, "ALP", "W/m2")
1026	CALL histdef2d(iff,o_cin%flag,o_cin%name, "Convective INhibition", "m2/s2")
1027	CALL histdef2d(iff,o_wape%flag,o_WAPE%name, "WAPE", "m2/s2")
1028	CALL histdef2d(iff,o_wake_h%flag,o_wake_h%name, "wake_h", "-")
1029	CALL histdef2d(iff,o_wake_s%flag,o_wake_s%name, "wake_s", "-")
1030	CALL histdef3d(iff,o_dtwak%flag,o_dtwak%name, "Wake dT", "K/s")
1031	CALL histdef3d(iff,o_dqwak%flag,o_dqwak%name, "Wake dQ", "(kg/kg)/s")
1032	CALL histdef3d(iff,o_wake_deltat%flag,o_wake_deltat%name, "wake_deltat", " ")
1033	CALL histdef3d(iff,o_wake_deltaq%flag,o_wake_deltaq%name, "wake_deltaq", " ")
1034	CALL histdef3d(iff,o_wake_omg%flag,o_wake_omg%name, "wake_omg", " ")
1035	ENDIF
1036	CALL histdef3d(iff,o_Vprecip%flag,o_Vprecip%name, "precipitation vertical profile", "-")
1037	CALL histdef3d(iff,o_ftd%flag,o_ftd%name, "tend temp due aux descentes precip", "-")
1038	CALL histdef3d(iff,o_fqd%flag,o_fqd%name,"tend vap eau due aux descentes precip", "-")
1039	ENDIF !(iflag_con.EQ.3)
1040
1041	CALL histdef3d(iff,o_dtlsc%flag,o_dtlsc%name, "Condensation dT", "K/s")
1042	CALL histdef3d(iff,o_dtlschr%flag,o_dtlschr%name,"Large-scale condensational heating rate","K/s")
1043	CALL histdef3d(iff,o_dqlsc%flag,o_dqlsc%name, "Condensation dQ", "(kg/kg)/s")
1044	CALL histdef3d(iff,o_dtvdf%flag,o_dtvdf%name, "Boundary-layer dT", "K/s")
1045	CALL histdef3d(iff,o_dqvdf%flag,o_dqvdf%name, "Boundary-layer dQ", "(kg/kg)/s")
1046	CALL histdef3d(iff,o_dteva%flag,o_dteva%name, "Reevaporation dT", "K/s")
1047	CALL histdef3d(iff,o_dqeva%flag,o_dqeva%name, "Reevaporation dQ", "(kg/kg)/s")
1048	CALL histdef3d(iff,o_ptconv%flag,o_ptconv%name, "POINTS CONVECTIFS", " ")
1049	CALL histdef3d(iff,o_ratqs%flag,o_ratqs%name, "RATQS", " ")
1050	CALL histdef3d(iff,o_dtthe%flag,o_dtthe%name, "Dry adjust. dT", "K/s")
1051
1052	if(iflag_thermals.gt.1) THEN
1053	CALL histdef3d(iff,o_f_th%flag,o_f_th%name, "Thermal plume mass flux", "K/s")
1054	CALL histdef3d(iff,o_e_th%flag,o_e_th%name,"Thermal plume entrainment","K/s")
1055	CALL histdef3d(iff,o_w_th%flag,o_w_th%name,"Thermal plume vertical velocity","m/s")
1056	CALL histdef3d(iff,o_lambda_th%flag,o_lambda_th%name,"Thermal plume vertical velocity","m/s")
1057	CALL histdef3d(iff,o_q_th%flag,o_q_th%name, "Thermal plume total humidity", "kg/kg")
1058	CALL histdef3d(iff,o_a_th%flag,o_a_th%name, "Thermal plume fraction", "")
1059	CALL histdef3d(iff,o_d_th%flag,o_d_th%name, "Thermal plume detrainment", "K/s")
1060	endif !iflag_thermals.gt.1
1061	CALL histdef2d(iff,o_f0_th%flag,o_f0_th%name, "Thermal closure mass flux", "K/s")
1062	CALL histdef2d(iff,o_zmax_th%flag,o_zmax_th%name, "Thermal plume height", "K/s")
1063	CALL histdef3d(iff,o_dqthe%flag,o_dqthe%name, "Dry adjust. dQ", "(kg/kg)/s")
1064	CALL histdef3d(iff,o_dtajs%flag,o_dtajs%name, "Dry adjust. dT", "K/s")
1065	CALL histdef3d(iff,o_dqajs%flag,o_dqajs%name, "Dry adjust. dQ", "(kg/kg)/s")
1066	CALL histdef3d(iff,o_dtswr%flag,o_dtswr%name, "SW radiation dT", "K/s")
1067	CALL histdef3d(iff,o_dtsw0%flag,o_dtsw0%name, "CS SW radiation dT", "K/s")
1068	CALL histdef3d(iff,o_dtlwr%flag,o_dtlwr%name, "LW radiation dT", "K/s")
1069	CALL histdef3d(iff,o_dtlw0%flag,o_dtlw0%name, "CS LW radiation dT", "K/s")
1070	CALL histdef3d(iff,o_dtec%flag,o_dtec%name, "Cinetic dissip dT", "K/s")
1071	CALL histdef3d(iff,o_duvdf%flag,o_duvdf%name, "Boundary-layer dU", "m/s2")
1072	CALL histdef3d(iff,o_dvvdf%flag,o_dvvdf%name, "Boundary-layer dV", "m/s2")
1073
1074	IF (ok_orodr) THEN
1075	CALL histdef3d(iff,o_duoro%flag,o_duoro%name, "Orography dU", "m/s2")
1076	CALL histdef3d(iff,o_dvoro%flag,o_dvoro%name, "Orography dV", "m/s2")
1077	ENDIF
1078
1079	IF (ok_orolf) THEN
1080	CALL histdef3d(iff,o_dulif%flag,o_dulif%name, "Orography dU", "m/s2")
1081	CALL histdef3d(iff,o_dvlif%flag,o_dvlif%name, "Orography dV", "m/s2")
1082	ENDIF
1083
1084	if (nqtot>=3) THEN
1085	!Attention DO iq=3,nqtot
1086	DO iq=3,4
1087	iiq=niadv(iq)
1088	! CALL histdef3d (iff, o_trac%flag,'o_'//tnom(iq)%name,ttext(iiq), "-" )
1089	CALL histdef3d (iff, o_trac(iq-2)%flag,o_trac(iq-2)%name,ttext(iiq), "-" )
1090	ENDDO
1091	endif
1092
1093	CALL histend(nid_files(iff))
1094
1095	ndex2d = 0
1096	ndex3d = 0
1097
1098	ENDIF ! clef_files
1099
1100	ENDDO ! iff
1101	print*,'Fin phys_output_mod.F90'
1102	end subroutine phys_output_open
1103
1104	SUBROUTINE histdef2d (iff,flag_var,nomvar,titrevar,unitvar)
1105
1106	use ioipsl
1107	USE dimphy
1108	USE mod_phys_lmdz_para
1109
1110	IMPLICIT NONE
1111
1112	include "dimensions.h"
1113	include "temps.h"
1114	include "indicesol.h"
1115	include "clesphys.h"
1116
1117	integer :: iff
1118	integer, dimension(nfiles) :: flag_var
1119	character(len=20) :: nomvar
1120	character(len=*) :: titrevar
1121	character(len=*) :: unitvar
1122
1123	real zstophym
1124
1125	if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then
1126	zstophym=zoutm(iff)
1127	else
1128	zstophym=zdtime
1129	endif
1130
1131	! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def
1132	call conf_physoutputs(nomvar,flag_var)
1133
1134	if ( flag_var(iff)<=lev_files(iff) ) then
1135	call histdef (nid_files(iff),nomvar,titrevar,unitvar, &
1136	iim,jj_nb,nhorim(iff), 1,1,1, -99, 32, &
1137	type_ecri(iff), zstophym,zoutm(iff))
1138	endif
1139	end subroutine histdef2d
1140
1141	SUBROUTINE histdef3d (iff,flag_var,nomvar,titrevar,unitvar)
1142
1143	use ioipsl
1144	USE dimphy
1145	USE mod_phys_lmdz_para
1146
1147	IMPLICIT NONE
1148
1149	include "dimensions.h"
1150	include "temps.h"
1151	include "indicesol.h"
1152	include "clesphys.h"
1153
1154	integer :: iff
1155	integer, dimension(nfiles) :: flag_var
1156	character(len=20) :: nomvar
1157	character(len=*) :: titrevar
1158	character(len=*) :: unitvar
1159
1160	real zstophym
1161
1162	! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def
1163	call conf_physoutputs(nomvar,flag_var)
1164
1165	if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then
1166	zstophym=zoutm(iff)
1167	else
1168	zstophym=zdtime
1169	endif
1170
1171	if ( flag_var(iff)<=lev_files(iff) ) then
1172	call histdef (nid_files(iff), nomvar, titrevar, unitvar, &
1173	iim, jj_nb, nhorim(iff), klev, levmin(iff), &
1174	levmax(iff)-levmin(iff)+1, nvertm(iff), 32, type_ecri(iff), &
1175	zstophym, zoutm(iff))
1176	endif
1177	end subroutine histdef3d
1178
1179	SUBROUTINE conf_physoutputs(nam_var,flag_var)
1180	!!! Lecture des noms et niveau de sortie des variables dans output.def
1181	! en utilisant les routines getin de IOIPSL
1182	use ioipsl
1183
1184	IMPLICIT NONE
1185
1186	include 'iniprint.h'
1187
1188	character(len=20) :: nam_var
1189	integer, dimension(nfiles) :: flag_var
1190
1191	IF(prt_level>10) WRITE(lunout,*)'Avant getin: nam_var flag_var ',nam_var,flag_var(:)
1192	call getin('flag_'//nam_var,flag_var)
1193	call getin('name_'//nam_var,nam_var)
1194	IF(prt_level>10) WRITE(lunout,*)'Apres getin: nam_var flag_var ',nam_var,flag_var(:)
1195
1196	END SUBROUTINE conf_physoutputs
1197
1198	SUBROUTINE convers_timesteps(str,timestep)
1199
1200	use ioipsl
1201
1202	IMPLICIT NONE
1203
1204	character(len=20) :: str
1205	character(len=10) :: type
1206	integer :: ipos,il
1207	real :: ttt,xxx,timestep,dayseconde
1208	parameter (dayseconde=86400.)
1209	include "temps.h"
1210	include "comconst.h"
1211
1212	ipos=scan(str,'0123456789.',.true.)
1213	!
1214	il=len_trim(str)
1215	print*,ipos,il
1216	read(str(1:ipos),*) ttt
1217	print*,ttt
1218	type=str(ipos+1:il)
1219
1220
1221	if ( il == ipos ) then
1222	type='day'
1223	endif
1224
1225	if ( type == 'day'.or.type == 'days'.or.type == 'jours'.or.type == 'jour' ) timestep = ttt * dayseconde
1226	if ( type == 'mounths'.or.type == 'mth'.or.type == 'mois' ) then
1227	print*,'annee_ref,day_ref mon_len',annee_ref,day_ref,ioget_mon_len(annee_ref,day_ref)
1228	timestep = ttt * dayseconde * ioget_mon_len(annee_ref,day_ref)
1229	endif
1230	if ( type == 'hours'.or.type == 'hr'.or.type == 'heurs') timestep = ttt * dayseconde / 24.
1231	if ( type == 'mn'.or.type == 'minutes' ) timestep = ttt * 60.
1232	if ( type == 's'.or.type == 'sec'.or.type == 'secondes' ) timestep = ttt
1233	if ( type == 'TS' ) timestep = dtphys
1234
1235	print*,'type = ',type
1236	print*,'nb j/h/m = ',ttt
1237	print*,'timestep(s)=',timestep
1238
1239	END SUBROUTINE convers_timesteps
1240
1241	END MODULE phys_output_mod
1242

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

Download in other formats:

Original Format