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