Context Navigation

phys_output_ctrlout_mod.F90 @ 2668

Last change on this file since 2668 was 2665, checked in by dcugnet, 8 years ago
A (re)startphy.nc file (standard name: "startphy0.nc") can be read by ce0l to get land mask, so mask can be defined (in decreasing priority order) from: 1) "o2a.nc file" if this file is found 2) "startphy0.nc" if this file is found 3) "Relief.nc" otherwise Sub-cell scales parameters for orographic gravity waves can be read from file "oro_params.nc" if the configuration key "read_orop" is TRUE. The effect is to bypass the "grid_noro" routine in ce0l, so that any pre-defined mask (from o2a.nc or startphy0.nc) is then overwritten. The gcm stops if the "limit.nc" records number differs from the current year number of days. A warning is issued in case the gcm calendar does not match the time axis attribute "calendar" (if available) from the "limit.nc" file. This attribute is now added to the "limit.nc" time axis. Few simplifications in grid_noro Few parameters changes in acama_gwd and flott_gwd. Variable d_u can be saved in the outputs.
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory
File size: 117.1 KB

Line
1	MODULE phys_output_ctrlout_mod
2
3	USE phys_output_var_mod
4	USE indice_sol_mod
5	USE aero_mod
6
7
8
9	IMPLICIT NONE
10	INTEGER, PRIVATE :: i
11
12	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
13	!! Definition pour chaque variable du niveau d ecriture dans chaque fichier,
14	!! de son nom, de sa description, de son unité et du type d'écriture.
15	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!/ histmth, histday, histhf, histins /),'!!!!!!!!!!!!
16	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
17	! CHARACTER(len=20), dimension(nfiles) :: TEF = type_ecri_files
18
19	!!! Comosantes de la coordonnee sigma-hybride
20	!!! Ap et Bp
21	TYPE(ctrl_out), SAVE :: o_Ahyb = ctrl_out((/ 1, 1, 1, 1, 1, 1, 11, 11, 11 /), &
22	'Ap', '', '', (/ ('', i=1, 9) /))
23	TYPE(ctrl_out), SAVE :: o_Bhyb = ctrl_out((/ 1, 1, 1, 1, 1, 1, 11, 11, 11 /), &
24	'Bp', '', '', (/ ('', i=1, 9) /))
25	TYPE(ctrl_out), SAVE :: o_Alt = ctrl_out((/ 1, 1, 1, 1, 1, 1, 11, 11, 11 /), &
26	'Alt', '', '', (/ ('', i=1, 9) /))
27
28	!!! 1D
29	TYPE(ctrl_out), SAVE :: o_phis = ctrl_out((/ 1, 1, 10, 5, 1, 1, 11, 11, 11 /), &
30	'phis', 'Surface geop.height', 'm2/s2', (/ ('', i=1, 9) /))
31	TYPE(ctrl_out), SAVE :: o_aire = ctrl_out((/ 1, 1, 10, 10, 1, 1, 11, 11, 11 /), &
32	'aire', 'Grid area', '-', (/ 'once', 'once', 'once', 'once', 'once', 'once', &
33	'once', 'once', 'once' /))
34	TYPE(ctrl_out), SAVE :: o_contfracATM = ctrl_out((/ 10, 1, 1, 10, 10, 10, 11, 11, 11 /), &
35	'contfracATM', '% sfce ter+lic', '-', &
36	(/ 'once', 'once', 'once', 'once', 'once', 'once', 'once', 'once', 'once' /))
37	TYPE(ctrl_out), SAVE :: o_contfracOR = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
38	'contfracOR', '% sfce terre OR', '-', (/ ('', i=1, 9) /))
39	TYPE(ctrl_out), SAVE :: o_aireTER = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
40	'aireTER', 'Grid area CONT', '-', (/ ('', i=1, 9) /))
41
42	!!! 2D
43
44	! Marine
45
46	TYPE(ctrl_out), SAVE :: o_alt_tropo = ctrl_out((/1,1,1,1,1,10,10,10,10/),&
47	'alt_tropo','Tropopause pressure','hPa',&
48	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
49	"inst(X)", "inst(X)","inst(X)" /))
50
51	TYPE(ctrl_out), SAVE :: o_map_prop_hc = ctrl_out((/1,1,1,1,1,10,10,10,10/),&
52	'map_prop_hc','Proportion of high clouds',' ',&
53	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
54	"inst(X)", "inst(X)","inst(X)" /))
55
56	TYPE(ctrl_out), SAVE :: o_map_prop_hist = &
57	ctrl_out((/1,1,1,1,1,1,10,10,10/),&
58	'map_prop_hist','Proportion of high ice semi-transp clouds',' ',&
59	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
60	"inst(X)", "inst(X)","inst(X)" /))
61
62	TYPE(ctrl_out), SAVE :: o_map_emis_hc = &
63	ctrl_out((/1,1,1,1,1,1,10,10,10/),&
64	'map_emis_hc','Emissivity of high clouds',' ',&
65	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
66	"inst(X)", "inst(X)","inst(X)" /))
67
68	TYPE(ctrl_out), SAVE :: o_map_iwp_hc = &
69	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
70	'map_iwp_hc','Ice water path of high clouds','g/m2',&
71	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
72	"inst(X)", "inst(X)","inst(X)" /))
73
74	TYPE(ctrl_out), SAVE :: o_map_deltaz_hc = &
75	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
76	'map_deltaz_hc','geom thickness of high clouds','m',&
77	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
78	"inst(X)", "inst(X)","inst(X)" /))
79
80	TYPE(ctrl_out), SAVE :: o_map_pcld_hc = &
81	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
82	'map_pcld_hc','cloud pressure of high clouds','hPa',&
83	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
84	"inst(X)", "inst(X)","inst(X)" /))
85
86	TYPE(ctrl_out), SAVE :: o_map_tcld_hc = &
87	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
88	'map_tcld_hc','cloud temperature of high clouds','K',&
89	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
90	"inst(X)", "inst(X)","inst(X)" /))
91
92
93	TYPE(ctrl_out), SAVE :: o_map_emis_hist = &
94	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
95	'map_emis_hist','Emissivity of high ice st clouds',' ',&
96	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
97	"inst(X)", "inst(X)","inst(X)" /))
98
99	TYPE(ctrl_out), SAVE :: o_map_iwp_hist = &
100	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
101	'map_iwp_hist','Ice water path of high ice st clouds','g/m2',&
102	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
103	"inst(X)", "inst(X)","inst(X)" /))
104
105	TYPE(ctrl_out), SAVE :: o_map_deltaz_hist = &
106	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
107	'map_deltaz_hist','geom thickness of high ice st clouds','m',&
108	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
109	"inst(X)", "inst(X)","inst(X)" /))
110
111	TYPE(ctrl_out), SAVE :: o_map_rad_hist = &
112	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
113	'map_rad_hist','ice crystals radius in high ice st clouds','µm',&
114	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
115	"inst(X)", "inst(X)","inst(X)" /))
116
117
118	TYPE(ctrl_out), SAVE :: o_map_emis_Cb = &
119	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
120	'map_emis_Cb','Emissivity of high Cb clouds',' ',&
121	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
122	"inst(X)", "inst(X)","inst(X)" /))
123
124	TYPE(ctrl_out), SAVE :: o_map_pcld_Cb = &
125	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
126	'map_pcld_Cb','cloud pressure of high Cb clouds','hPa',&
127	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
128	"inst(X)", "inst(X)","inst(X)" /))
129
130	TYPE(ctrl_out), SAVE :: o_map_tcld_Cb = &
131	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
132	'map_tcld_Cb','cloud temperature of high Cb clouds','K',&
133	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
134	"inst(X)", "inst(X)","inst(X)" /))
135
136
137	TYPE(ctrl_out), SAVE :: o_map_emis_Anv = &
138	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
139	'map_emis_Anv','Emissivity of high Anv clouds',' ',&
140	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
141	"inst(X)", "inst(X)","inst(X)" /))
142
143	TYPE(ctrl_out), SAVE :: o_map_pcld_Anv = &
144	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
145	'map_pcld_Anv','cloud pressure of high Anv clouds','hPa',&
146	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
147	"inst(X)", "inst(X)","inst(X)" /))
148
149	TYPE(ctrl_out), SAVE :: o_map_tcld_Anv = &
150	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
151	'map_tcld_Anv','cloud temperature of high Anv clouds','K',&
152	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
153	"inst(X)", "inst(X)","inst(X)" /))
154
155	TYPE(ctrl_out), SAVE :: o_map_emis_ThCi = &
156	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
157	'map_emis_ThCi','Emissivity of high ThCi clouds',' ',&
158	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
159	"inst(X)", "inst(X)","inst(X)" /))
160
161	TYPE(ctrl_out), SAVE :: o_map_pcld_ThCi = &
162	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
163	'map_pcld_ThCi','cloud pressure of high ThCi clouds','hPa',&
164	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
165	"inst(X)", "inst(X)","inst(X)" /))
166
167	TYPE(ctrl_out), SAVE :: o_map_tcld_ThCi = &
168	ctrl_out((/10,10,1,10,10,10,10,10,10/),&
169	'map_tcld_ThCi','cloud temperature of high ThCi clouds','K',&
170	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
171	"inst(X)", "inst(X)","inst(X)" /))
172
173	TYPE(ctrl_out), SAVE :: o_map_ntot = &
174	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
175	'map_ntot','total AIRS cloud fraction',' ',&
176	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
177	"inst(X)", "inst(X)","inst(X)" /))
178
179	TYPE(ctrl_out), SAVE :: o_map_hc = &
180	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
181	'map_hc','high clouds AIRS cloud fraction',' ',&
182	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
183	"inst(X)", "inst(X)","inst(X)" /))
184
185	TYPE(ctrl_out), SAVE :: o_map_hist = &
186	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
187	'map_hist','high clouds ice st AIRS cloud fraction',' ',&
188	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
189	"inst(X)", "inst(X)","inst(X)" /))
190
191	TYPE(ctrl_out), SAVE :: o_map_Cb = &
192	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
193	'map_Cb','high clouds Cb AIRS cloud fraction',' ',&
194	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
195	"inst(X)", "inst(X)","inst(X)" /))
196
197	TYPE(ctrl_out), SAVE :: o_map_ThCi = &
198	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
199	'map_ThCi','high clouds ThCi AIRS cloud fraction',' ',&
200	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
201	"inst(X)", "inst(X)","inst(X)" /))
202
203	TYPE(ctrl_out), SAVE :: o_map_Anv = &
204	ctrl_out((/1,1,1,1,1,10,10,10,10/),&
205	'map_Anv','high clouds Anv AIRS cloud fraction',' ',&
206	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
207	"inst(X)", "inst(X)","inst(X)" /))
208
209
210	! Fin Marine
211
212	TYPE(ctrl_out), SAVE :: o_flat = ctrl_out((/ 5, 1, 10, 10, 5, 10, 11, 11, 11 /), &
213	'flat', 'Latent heat flux', 'W/m2', (/ ('', i=1, 9) /))
214	TYPE(ctrl_out), SAVE :: o_ptstar = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
215	'ptstar', 'Air Surface Temperature', 'K', (/ ('', i=1, 9) /))
216	TYPE(ctrl_out), SAVE :: o_pt0 = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
217	'pt0', 'Standard Air Surface Temperature', 'K', (/ ('', i=1, 9) /))
218	TYPE(ctrl_out), SAVE :: o_slp = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
219	'slp', 'Sea Level Pressure', 'Pa', (/ ('', i=1, 9) /))
220	TYPE(ctrl_out), SAVE :: o_tsol = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
221	'tsol', 'Surface Temperature', 'K', (/ ('', i=1, 9) /))
222	TYPE(ctrl_out), SAVE :: o_t2m = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
223	't2m', 'Temperature 2m', 'K', (/ ('', i=1, 9) /))
224	TYPE(ctrl_out), SAVE :: o_t2m_min = ctrl_out((/ 20, 1, 10, 10, 10, 10, 11, 11, 11 /), &
225	't2m_min', 'Temp 2m min', 'K', &
226	(/ "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)" /))
227	TYPE(ctrl_out), SAVE :: o_t2m_max = ctrl_out((/ 20, 1, 10, 10, 10, 10, 11, 11, 11 /), &
228	't2m_max', 'Temp 2m max', 'K', &
229	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
230	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /))
231
232	TYPE(ctrl_out), SAVE :: o_t2m_min_mon = ctrl_out((/ 1, 20, 20, 20, 20, 20, 20, 20, 20 /), &
233	't2m_min_mon', 'Monthly average min 2m temperature', 'K', &
234	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
235	TYPE(ctrl_out), SAVE :: o_t2m_max_mon = ctrl_out((/ 1, 20, 20, 20, 20, 20, 20, 20, 20 /), &
236	't2m_max_mon', 'Monthly average max 2m temperature', 'K', &
237	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
238	"inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
239
240	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_t2m_srf = (/ &
241	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /), &
242	't2m_ter', "Temp 2m "//clnsurf(1), "K", (/ ('', i=1, 9) /)), &
243	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /), &
244	't2m_lic', "Temp 2m "//clnsurf(2), "K", (/ ('', i=1, 9) /)), &
245	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /), &
246	't2m_oce', "Temp 2m "//clnsurf(3), "K", (/ ('', i=1, 9) /)), &
247	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /), &
248	't2m_sic', "Temp 2m "//clnsurf(4), "K", (/ ('', i=1, 9) /)) /)
249
250	TYPE(ctrl_out), SAVE :: o_gusts = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
251	'gusts', 'surface gustiness', 'm2/s2', (/ ('', i=1, 9) /))
252
253	TYPE(ctrl_out), SAVE :: o_wind10m = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
254	'wind10m', '10-m wind speed', 'm/s', (/ ('', i=1, 9) /))
255	TYPE(ctrl_out), SAVE :: o_wind10max = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
256	'wind10max', '10m wind speed max', 'm/s', &
257	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
258	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /))
259	TYPE(ctrl_out), SAVE :: o_sicf = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
260	'sicf', 'Sea-ice fraction', '-', (/ ('', i=1, 9) /))
261	TYPE(ctrl_out), SAVE :: o_q2m = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
262	'q2m', 'Specific humidity 2m', 'kg/kg', (/ ('', i=1, 9) /))
263	TYPE(ctrl_out), SAVE :: o_ustar = ctrl_out((/ 1, 1, 10, 5, 10, 10, 11, 11, 11 /), &
264	'ustar', 'Friction velocity', 'm/s', (/ ('', i=1, 9) /))
265	TYPE(ctrl_out), SAVE :: o_u10m = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
266	'u10m', 'Vent zonal 10m', 'm/s', (/ ('', i=1, 9) /))
267	TYPE(ctrl_out), SAVE :: o_v10m = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
268	'v10m', 'Vent meridien 10m', 'm/s', (/ ('', i=1, 9) /))
269	TYPE(ctrl_out), SAVE :: o_psol = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
270	'psol', 'Surface Pressure', 'Pa', (/ ('', i=1, 9) /))
271	TYPE(ctrl_out), SAVE :: o_qsurf = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
272	'qsurf', 'Surface Air humidity', 'kg/kg', (/ ('', i=1, 9) /))
273
274	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_ustar_srf = (/ &
275	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'ustar_ter', &
276	"Friction velocity "//clnsurf(1),"m/s", (/ ('', i=1, 9) /)), &
277	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'ustar_lic', &
278	"Friction velocity "//clnsurf(2),"m/s", (/ ('', i=1, 9) /)), &
279	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'ustar_oce', &
280	"Friction velocity "//clnsurf(3),"m/s", (/ ('', i=1, 9) /)), &
281	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'ustar_sic', &
282	"Friction velocity "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)) /)
283
284	TYPE(ctrl_out), SAVE, DIMENSION(5) :: o_wstar = (/ &
285	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'wstar_ter', &
286	"Friction velocity "//clnsurf(1),"m/s", (/ ('', i=1, 9) /)), &
287	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'wstar_lic', &
288	"Friction velocity "//clnsurf(2),"m/s", (/ ('', i=1, 9) /)), &
289	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'wstar_oce', &
290	"Friction velocity "//clnsurf(3),"m/s", (/ ('', i=1, 9) /)), &
291	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'wstar_sic', &
292	"Friction velocity "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)), &
293	ctrl_out((/ 5, 5, 10, 10, 10, 10, 11, 11, 11 /),'wstar', &
294	"w* convective velocity "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)) /)
295
296	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_u10m_srf = (/ &
297	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'u10m_ter', &
298	"Vent Zonal 10m "//clnsurf(1),"m/s", (/ ('', i=1, 9) /)), &
299	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'u10m_lic', &
300	"Vent Zonal 10m "//clnsurf(2),"m/s", (/ ('', i=1, 9) /)), &
301	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'u10m_oce', &
302	"Vent Zonal 10m "//clnsurf(3),"m/s", (/ ('', i=1, 9) /)), &
303	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'u10m_sic', &
304	"Vent Zonal 10m "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)) /)
305
306	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_v10m_srf = (/ &
307	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'v10m_ter', &
308	"Vent meredien 10m "//clnsurf(1),"m/s", (/ ('', i=1, 9) /)), &
309	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'v10m_lic', &
310	"Vent meredien 10m "//clnsurf(2),"m/s", (/ ('', i=1, 9) /)), &
311	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'v10m_oce', &
312	"Vent meredien 10m "//clnsurf(3),"m/s", (/ ('', i=1, 9) /)), &
313	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'v10m_sic', &
314	"Vent meredien 10m "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)) /)
315
316	TYPE(ctrl_out), SAVE :: o_qsol = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
317	'qsol', 'Soil watter content', 'mm', (/ ('', i=1, 9) /))
318	TYPE(ctrl_out), SAVE :: o_ndayrain = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
319	'ndayrain', 'Number of dayrain(liq+sol)', '-', &
320	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
321	TYPE(ctrl_out), SAVE :: o_precip = ctrl_out((/ 1, 1, 1, 10, 5, 10, 11, 11, 11 /), &
322	'precip', 'Precip Totale liq+sol', 'kg/(s*m2)', (/ ('', i=1, 9) /))
323	TYPE(ctrl_out), SAVE :: o_plul = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
324	'plul', 'Large-scale Precip.', 'kg/(s*m2)', (/ ('', i=1, 9) /))
325	TYPE(ctrl_out), SAVE :: o_plun = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
326	'plun', 'Numerical Precip.', 'kg/(s*m2)', (/ ('', i=1, 9) /))
327	TYPE(ctrl_out), SAVE :: o_pluc = ctrl_out((/ 1, 1, 1, 10, 5, 10, 11, 11, 11 /), &
328	'pluc', 'Convective Precip.', 'kg/(s*m2)', (/ ('', i=1, 9) /))
329	TYPE(ctrl_out), SAVE :: o_snow = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11 /), &
330	'snow', 'Snow fall', 'kg/(s*m2)', (/ ('', i=1, 9) /))
331	TYPE(ctrl_out), SAVE :: o_evap = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
332	'evap', 'Evaporat', 'kg/(s*m2)', (/ ('', i=1, 9) /))
333
334	TYPE(ctrl_out), SAVE :: o_sens_prec_liq_oce = ctrl_out((/ 5, 5, 10, 10, 5, 10, 11, 11, 11 /), &
335	'sens_rain_oce', 'Sensible heat flux of liquid prec. over ocean', 'W/m2', (/ ('', i=1, 9) /))
336	TYPE(ctrl_out), SAVE :: o_sens_prec_liq_sic = ctrl_out((/ 5, 5, 10, 10, 5, 10, 11, 11, 11 /), &
337	'sens_rain_sic', 'Sensible heat flux of liquid prec. over seaice', 'W/m2', (/ ('', i=1, 9) /))
338	TYPE(ctrl_out), SAVE :: o_sens_prec_sol_oce = ctrl_out((/ 5, 5, 10, 10, 5, 10, 11, 11, 11 /), &
339	'sens_snow_oce', 'Sensible heat flux of solid prec. over ocean', 'W/m2', (/ ('', i=1, 9) /))
340	TYPE(ctrl_out), SAVE :: o_sens_prec_sol_sic = ctrl_out((/ 5, 5, 10, 10, 5, 10, 11, 11, 11 /), &
341	'sens_snow_sic', 'Sensible heat flux of solid prec. over seaice', 'W/m2', (/ ('', i=1, 9) /))
342	TYPE(ctrl_out), SAVE :: o_lat_prec_liq_oce = ctrl_out((/ 5, 5, 10, 10, 5, 10, 11, 11, 11 /), &
343	'lat_rain_oce', 'Latent heat flux of liquid prec. over ocean', 'W/m2', (/ ('', i=1, 9) /))
344	TYPE(ctrl_out), SAVE :: o_lat_prec_liq_sic = ctrl_out((/ 5, 5, 10, 10, 5, 10, 11, 11, 11 /), &
345	'lat_rain_sic', 'Latent heat flux of liquid prec. over seaice', 'W/m2', (/ ('', i=1, 9) /))
346	TYPE(ctrl_out), SAVE :: o_lat_prec_sol_oce = ctrl_out((/ 5, 5, 10, 10, 5, 10, 11, 11, 11 /), &
347	'lat_snow_oce', 'Latent heat flux of solid prec. over ocean', 'W/m2', (/ ('', i=1, 9) /))
348	TYPE(ctrl_out), SAVE :: o_lat_prec_sol_sic = ctrl_out((/ 5, 5, 10, 10, 5, 10, 11, 11, 11 /), &
349	'lat_snow_sic', 'Latent heat flux of solid prec. over seaice', 'W/m2', (/ ('', i=1, 9) /))
350
351
352	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_evap_srf = (/ &
353	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evap_ter', &
354	"evaporation at surface "//clnsurf(1),"kg/(s*m2)", (/ ('', i=1, 9) /)), &
355	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evap_lic', &
356	"evaporation at surface "//clnsurf(2),"kg/(s*m2)", (/ ('', i=1, 9) /)), &
357	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evap_oce', &
358	"evaporation at surface "//clnsurf(3),"kg/(s*m2)", (/ ('', i=1, 9) /)), &
359	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evap_sic', &
360	"evaporation at surface "//clnsurf(4),"kg/(s*m2)", (/ ('', i=1, 9) /)) /)
361
362	TYPE(ctrl_out), SAVE :: o_msnow = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
363	'msnow', 'Surface snow amount', 'kg/m2', (/ ('', i=1, 9) /))
364	TYPE(ctrl_out), SAVE :: o_fsnow = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
365	'fsnow', 'Surface snow area fraction', '-', (/ ('', i=1, 9) /))
366	TYPE(ctrl_out), SAVE :: o_tops = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
367	'tops', 'Solar rad. at TOA', 'W/m2', (/ ('', i=1, 9) /))
368	TYPE(ctrl_out), SAVE :: o_tops0 = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
369	'tops0', 'CS Solar rad. at TOA', 'W/m2', (/ ('', i=1, 9) /))
370	TYPE(ctrl_out), SAVE :: o_topl = ctrl_out((/ 1, 1, 10, 5, 10, 10, 11, 11, 11 /), &
371	'topl', 'IR rad. at TOA', 'W/m2', (/ ('', i=1, 9) /))
372	TYPE(ctrl_out), SAVE :: o_topl0 = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
373	'topl0', 'IR rad. at TOA', 'W/m2', (/ ('', i=1, 9) /))
374	TYPE(ctrl_out), SAVE :: o_SWupTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
375	'SWupTOA', 'SWup at TOA', 'W/m2', (/ ('', i=1, 9) /))
376	TYPE(ctrl_out), SAVE :: o_SWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
377	'SWupTOAclr', 'SWup clear sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
378	TYPE(ctrl_out), SAVE :: o_SWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
379	'SWdnTOA', 'SWdn at TOA', 'W/m2', (/ ('', i=1, 9) /))
380	TYPE(ctrl_out), SAVE :: o_SWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
381	'SWdnTOAclr', 'SWdn clear sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
382	TYPE(ctrl_out), SAVE :: o_nettop = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
383	'nettop', 'Net dn radiatif flux at TOA', 'W/m2', (/ ('', i=1, 9) /))
384	TYPE(ctrl_out), SAVE :: o_SWup200 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
385	'SWup200', 'SWup at 200mb', 'W/m2', (/ ('', i=1, 9) /))
386	TYPE(ctrl_out), SAVE :: o_SWup200clr = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
387	'SWup200clr', 'SWup clear sky at 200mb', 'W/m2', (/ ('', i=1, 9) /))
388	TYPE(ctrl_out), SAVE :: o_SWdn200 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
389	'SWdn200', 'SWdn at 200mb', 'W/m2', (/ ('', i=1, 9) /))
390	TYPE(ctrl_out), SAVE :: o_SWdn200clr = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
391	'SWdn200clr', 'SWdn clear sky at 200mb', 'W/m2', (/ ('', i=1, 9) /))
392
393	! arajouter
394	! type(ctrl_out),save :: o_LWupTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /),'LWupTOA', &
395	! (/ ('', i=1, 9) /))
396	! type(ctrl_out),save :: o_LWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /),'LWupTOAclr', &
397	! (/ ('', i=1, 9) /))
398	! type(ctrl_out),save :: o_LWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /),'LWdnTOA', &
399	! (/ ('', i=1, 9) /))
400	! type(ctrl_out),save :: o_LWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /),'LWdnTOAclr', &
401	! (/ ('', i=1, 9) /))
402	TYPE(ctrl_out), SAVE :: o_LWup200 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
403	'LWup200', 'LWup at 200mb', 'W/m2', (/ ('', i=1, 9) /))
404	TYPE(ctrl_out), SAVE :: o_LWup200clr = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
405	'LWup200clr', 'LWup clear sky at 200mb', 'W/m2', (/ ('', i=1, 9) /))
406	TYPE(ctrl_out), SAVE :: o_LWdn200 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
407	'LWdn200', 'LWdn at 200mb', 'W/m2', (/ ('', i=1, 9) /))
408	TYPE(ctrl_out), SAVE :: o_LWdn200clr = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
409	'LWdn200clr', 'LWdn clear sky at 200mb', 'W/m2', (/ ('', i=1, 9) /))
410	TYPE(ctrl_out), SAVE :: o_sols = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
411	'sols', 'Solar rad. at surf.', 'W/m2', (/ ('', i=1, 9) /))
412	TYPE(ctrl_out), SAVE :: o_sols0 = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
413	'sols0', 'Solar rad. at surf.', 'W/m2', (/ ('', i=1, 9) /))
414	TYPE(ctrl_out), SAVE :: o_soll = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
415	'soll', 'IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
416	TYPE(ctrl_out), SAVE :: o_soll0 = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
417	'soll0', 'IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
418	TYPE(ctrl_out), SAVE :: o_radsol = ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /), &
419	'radsol', 'Rayonnement au sol', 'W/m2', (/ ('', i=1, 9) /))
420	TYPE(ctrl_out), SAVE :: o_SWupSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
421	'SWupSFC', 'SWup at surface', 'W/m2', (/ ('', i=1, 9) /))
422	TYPE(ctrl_out), SAVE :: o_SWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
423	'SWupSFCclr', 'SWup clear sky at surface', 'W/m2', (/ ('', i=1, 9) /))
424	TYPE(ctrl_out), SAVE :: o_SWdnSFC = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11 /), &
425	'SWdnSFC', 'SWdn at surface', 'W/m2', (/ ('', i=1, 9) /))
426	TYPE(ctrl_out), SAVE :: o_SWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
427	'SWdnSFCclr', 'SWdn clear sky at surface', 'W/m2', (/ ('', i=1, 9) /))
428	TYPE(ctrl_out), SAVE :: o_LWupSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
429	'LWupSFC', 'Upwd. IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
430	TYPE(ctrl_out), SAVE :: o_LWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
431	'LWupSFCclr', 'CS Upwd. IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
432	TYPE(ctrl_out), SAVE :: o_LWdnSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
433	'LWdnSFC', 'Down. IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
434	TYPE(ctrl_out), SAVE :: o_LWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
435	'LWdnSFCclr', 'Down. CS IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
436	TYPE(ctrl_out), SAVE :: o_bils = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
437	'bils', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
438	TYPE(ctrl_out), SAVE :: o_bils_tke = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
439	'bils_tke', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
440	TYPE(ctrl_out), SAVE :: o_bils_diss = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
441	'bils_diss', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
442	TYPE(ctrl_out), SAVE :: o_bils_ec = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
443	'bils_ec', 'Surf. total heat flux correction', 'W/m2', (/ ('', i=1, 9) /))
444	TYPE(ctrl_out), SAVE :: o_bils_ech = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
445	'bils_ech', 'Surf. total heat flux correction', 'W/m2', (/ ('', i=1, 9) /))
446	TYPE(ctrl_out), SAVE :: o_bils_kinetic = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
447	'bils_kinetic', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
448	TYPE(ctrl_out), SAVE :: o_bils_enthalp = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
449	'bils_enthalp', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
450	TYPE(ctrl_out), SAVE :: o_bils_latent = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
451	'bils_latent', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
452	TYPE(ctrl_out), SAVE :: o_sens = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11 /), &
453	'sens', 'Sensible heat flux', 'W/m2', (/ ('', i=1, 9) /))
454	TYPE(ctrl_out), SAVE :: o_fder = ctrl_out((/ 1, 2, 10, 10, 10, 10, 11, 11, 11 /), &
455	'fder', 'Heat flux derivation', 'W/m2', (/ ('', i=1, 9) /))
456	TYPE(ctrl_out), SAVE :: o_ffonte = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
457	'ffonte', 'Thermal flux for snow melting', 'W/m2', (/ ('', i=1, 9) /))
458	TYPE(ctrl_out), SAVE :: o_fqcalving = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
459	'fqcalving', 'Ice Calving', 'kg/m2/s', (/ ('', i=1, 9) /))
460	TYPE(ctrl_out), SAVE :: o_fqfonte = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
461	'fqfonte', 'Land ice melt', 'kg/m2/s', (/ ('', i=1, 9) /))
462	TYPE(ctrl_out), SAVE :: o_runofflic = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
463	'runofflic', 'Land ice melt to ocean', 'kg/m2/s', (/ ('', i=1, 9) /))
464	TYPE(ctrl_out), SAVE :: o_taux = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
465	'taux', 'Zonal wind stress', 'Pa', (/ ('', i=1, 9) /))
466	TYPE(ctrl_out), SAVE :: o_tauy = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
467	'tauy', 'Meridional wind stress', 'Pa', (/ ('', i=1, 9) /))
468
469	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_taux_srf = (/ &
470	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'taux_ter', &
471	"Zonal wind stress"//clnsurf(1), "Pa", (/ ('', i=1, 9) /)), &
472	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'taux_lic', &
473	"Zonal wind stress"//clnsurf(2), "Pa", (/ ('', i=1, 9) /)), &
474	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'taux_oce', &
475	"Zonal wind stress"//clnsurf(3), "Pa", (/ ('', i=1, 9) /)), &
476	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'taux_sic', &
477	"Zonal wind stress"//clnsurf(4), "Pa", (/ ('', i=1, 9) /)) /)
478
479	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tauy_srf = (/ &
480	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tauy_ter', &
481	"Meridional wind stress "//clnsurf(1),"Pa", (/ ('', i=1, 9) /)), &
482	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tauy_lic', &
483	"Meridional wind stress "//clnsurf(2),"Pa", (/ ('', i=1, 9) /)), &
484	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tauy_oce', &
485	"Meridional wind stress "//clnsurf(3),"Pa", (/ ('', i=1, 9) /)), &
486	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tauy_sic', &
487	"Meridional wind stress "//clnsurf(4),"Pa", (/ ('', i=1, 9) /)) /)
488
489	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_pourc_srf = (/ &
490	ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /),'pourc_ter', &
491	"% "//clnsurf(1),"%", (/ ('', i=1, 9) /)), &
492	ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /),'pourc_lic', &
493	"% "//clnsurf(2),"%", (/ ('', i=1, 9) /)), &
494	ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /),'pourc_oce', &
495	"% "//clnsurf(3),"%", (/ ('', i=1, 9) /)), &
496	ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /),'pourc_sic', &
497	"% "//clnsurf(4),"%", (/ ('', i=1, 9) /)) /)
498
499	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_fract_srf = (/ &
500	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'fract_ter', &
501	"Fraction "//clnsurf(1),"1", (/ ('', i=1, 9) /)), &
502	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'fract_lic', &
503	"Fraction "//clnsurf(2),"1", (/ ('', i=1, 9) /)), &
504	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'fract_oce', &
505	"Fraction "//clnsurf(3),"1", (/ ('', i=1, 9) /)), &
506	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'fract_sic', &
507	"Fraction "//clnsurf(4),"1", (/ ('', i=1, 9) /)) /)
508
509	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tsol_srf = (/ &
510	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tsol_ter', &
511	"Temperature "//clnsurf(1),"K", (/ ('', i=1, 9) /)), &
512	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tsol_lic', &
513	"Temperature "//clnsurf(2),"K", (/ ('', i=1, 9) /)), &
514	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tsol_oce', &
515	"Temperature "//clnsurf(3),"K", (/ ('', i=1, 9) /)), &
516	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tsol_sic', &
517	"Temperature "//clnsurf(4),"K", (/ ('', i=1, 9) /)) /)
518
519	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_evappot_srf = (/ &
520	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evappot_ter', &
521	"Temperature"//clnsurf(1),"K", (/ ('', i=1, 9) /)), &
522	ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /),'evappot_lic', &
523	"Temperature"//clnsurf(2),"K", (/ ('', i=1, 9) /)), &
524	ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /),'evappot_oce', &
525	"Temperature"//clnsurf(3),"K", (/ ('', i=1, 9) /)), &
526	ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /),'evappot_sic', &
527	"Temperature"//clnsurf(4),"K", (/ ('', i=1, 9) /)) /)
528
529	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_sens_srf = (/ &
530	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'sens_ter', &
531	"Sensible heat flux "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
532	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'sens_lic', &
533	"Sensible heat flux "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
534	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'sens_oce', &
535	"Sensible heat flux "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
536	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'sens_sic', &
537	"Sensible heat flux "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
538
539	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_lat_srf = (/ &
540	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'lat_ter', &
541	"Latent heat flux "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
542	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'lat_lic', &
543	"Latent heat flux "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
544	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'lat_oce', &
545	"Latent heat flux "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
546	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'lat_sic', &
547	"Latent heat flux "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
548
549	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_flw_srf = (/ &
550	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'flw_ter', &
551	"LW "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
552	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'flw_lic', &
553	"LW "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
554	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'flw_oce', &
555	"LW "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
556	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'flw_sic', &
557	"LW "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
558
559	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_fsw_srf = (/ &
560	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'fsw_ter', &
561	"SW "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
562	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'fsw_lic', &
563	"SW "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
564	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'fsw_oce', &
565	"SW "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
566	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'fsw_sic', &
567	"SW "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
568
569	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_wbils_srf = (/ &
570	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbils_ter', &
571	"Bilan sol "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
572	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbils_lic', &
573	"Bilan sol "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
574	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbils_oce', &
575	"Bilan sol "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
576	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbils_sic', &
577	"Bilan sol "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
578
579	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_wbilo_srf = (/ &
580	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbilo_ter', &
581	"Bilan eau "//clnsurf(1),"kg/(m2*s)", (/ ('', i=1, 9) /)), &
582	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbilo_lic', &
583	"Bilan eau "//clnsurf(2),"kg/(m2*s)", (/ ('', i=1, 9) /)), &
584	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbilo_oce', &
585	"Bilan eau "//clnsurf(3),"kg/(m2*s)", (/ ('', i=1, 9) /)), &
586	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbilo_sic', &
587	"Bilan eau "//clnsurf(4),"kg/(m2*s)", (/ ('', i=1, 9) /)) /)
588
589	TYPE(ctrl_out), SAVE :: o_cdrm = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
590	'cdrm', 'Momentum drag coef.', '-', (/ ('', i=1, 9) /))
591	TYPE(ctrl_out), SAVE :: o_cdrh = ctrl_out((/ 1, 10, 10, 7, 10, 10, 11, 11, 11 /), &
592	'cdrh', 'Heat drag coef.', '-', (/ ('', i=1, 9) /))
593	TYPE(ctrl_out), SAVE :: o_cldl = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
594	'cldl', 'Low-level cloudiness', '-', (/ ('', i=1, 9) /))
595	TYPE(ctrl_out), SAVE :: o_cldm = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
596	'cldm', 'Mid-level cloudiness', '-', (/ ('', i=1, 9) /))
597	TYPE(ctrl_out), SAVE :: o_cldh = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
598	'cldh', 'High-level cloudiness', '-', (/ ('', i=1, 9) /))
599	TYPE(ctrl_out), SAVE :: o_cldt = ctrl_out((/ 1, 1, 2, 10, 5, 10, 11, 11, 11 /), &
600	'cldt', 'Total cloudiness', '-', (/ ('', i=1, 9) /))
601	TYPE(ctrl_out), SAVE :: o_JrNt = ctrl_out((/ 1, 1, 10, 7, 10, 10, 11, 11, 11 /), &
602	'JrNt', '1 if Day 0 if Night', '-', (/ ('', i=1, 9) /))
603	TYPE(ctrl_out), SAVE :: o_cldhjn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
604	'cldhjn', 'High-level cloudiness Day', '-', (/ ('', i=1, 9) /))
605	TYPE(ctrl_out), SAVE :: o_cldmjn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11/), &
606	'cldmjn', 'Mid-level cloudiness day', '-', (/ ('', i=1, 9) /))
607	TYPE(ctrl_out), SAVE :: o_cldljn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11/), &
608	'cldljn', 'Low-level cloudiness day', '-', (/ ('', i=1, 9) /))
609	TYPE(ctrl_out), SAVE :: o_cldtjn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11/), &
610	'cldtjn', 'Total cloudiness day', '-', (/ ('', i=1, 9) /))
611
612	TYPE(ctrl_out), SAVE :: o_cldq = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
613	'cldq', 'Cloud liquid water path', 'kg/m2', (/ ('', i=1, 9) /))
614	TYPE(ctrl_out), SAVE :: o_lwp = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
615	'lwp', 'Cloud water path', 'kg/m2', (/ ('', i=1, 9) /))
616	TYPE(ctrl_out), SAVE :: o_iwp = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
617	'iwp', 'Cloud ice water path', 'kg/m2', (/ ('', i=1, 9) /))
618	TYPE(ctrl_out), SAVE :: o_ue = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
619	'ue', 'Zonal energy transport', '-', (/ ('', i=1, 9) /))
620	TYPE(ctrl_out), SAVE :: o_ve = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
621	've', 'Merid energy transport', '-', (/ ('', i=1, 9) /))
622	TYPE(ctrl_out), SAVE :: o_uq = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
623	'uq', 'Zonal humidity transport', '-', (/ ('', i=1, 9) /))
624	TYPE(ctrl_out), SAVE :: o_vq = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
625	'vq', 'Merid humidity transport', '-', (/ ('', i=1, 9) /))
626	TYPE(ctrl_out), SAVE :: o_cape = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
627	'cape', 'Conv avlbl pot ener', 'J/kg', (/ ('', i=1, 9) /))
628	TYPE(ctrl_out), SAVE :: o_pbase = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
629	'pbase', 'Cld base pressure', 'Pa', (/ ('', i=1, 9) /))
630	TYPE(ctrl_out), SAVE :: o_ptop = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
631	'ptop', 'Cld top pressure', 'Pa', (/ ('', i=1, 9) /))
632	TYPE(ctrl_out), SAVE :: o_fbase = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
633	'fbase', 'Cld base mass flux', 'kg/m2/s', (/ ('', i=1, 9) /))
634	TYPE(ctrl_out), SAVE :: o_plcl = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
635	'plcl', 'Lifting Condensation Level', 'hPa', (/ ('', i=1, 9) /))
636	TYPE(ctrl_out), SAVE :: o_plfc = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
637	'plfc', 'Level of Free Convection', 'hPa', (/ ('', i=1, 9) /))
638	TYPE(ctrl_out), SAVE :: o_wbeff = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
639	'wbeff', 'Conv. updraft velocity at LFC (<100)', 'm/s', (/ ('', i=1, 9) /))
640	TYPE(ctrl_out), SAVE :: o_prw = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
641	'prw', 'Precipitable water', 'kg/m2', (/ ('', i=1, 9) /))
642	TYPE(ctrl_out), SAVE :: o_prlw = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
643	'prlw', 'Precipitable liquid water', 'kg/m2', (/ ('', i=1, 9) /))
644	TYPE(ctrl_out), SAVE :: o_prsw = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
645	'prsw', 'Precipitable solid water', 'kg/m2', (/ ('', i=1, 9) /))
646	TYPE(ctrl_out), SAVE :: o_s_pblh = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
647	's_pblh', 'Boundary Layer Height', 'm', (/ ('', i=1, 9) /))
648	TYPE(ctrl_out), SAVE :: o_s_pblt = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
649	's_pblt', 't at Boundary Layer Height', 'K', (/ ('', i=1, 9) /))
650	TYPE(ctrl_out), SAVE :: o_s_lcl = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
651	's_lcl', 'Condensation level', 'm', (/ ('', i=1, 9) /))
652	TYPE(ctrl_out), SAVE :: o_s_therm = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
653	's_therm', 'Exces du thermique', 'K', (/ ('', i=1, 9) /))
654	!IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
655	! type(ctrl_out),save :: o_s_capCL = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_capCL', &
656	! (/ ('', i=1, 9) /))
657	! type(ctrl_out),save :: o_s_oliqCL = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_oliqCL', &
658	! (/ ('', i=1, 9) /))
659	! type(ctrl_out),save :: o_s_cteiCL = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_cteiCL', &
660	! (/ ('', i=1, 9) /))
661	! type(ctrl_out),save :: o_s_trmb1 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_trmb1', &
662	! (/ ('', i=1, 9) /))
663	! type(ctrl_out),save :: o_s_trmb2 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_trmb2', &
664	! (/ ('', i=1, 9) /))
665	! type(ctrl_out),save :: o_s_trmb3 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_trmb3', &
666	!(/ ('', i=1, 9) /))
667	TYPE(ctrl_out), SAVE :: o_slab_bils = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
668	'slab_bils', 'flux atmos - slab ponderes foce', 'W/m2', (/ ('', i=1, 9) /))
669	TYPE(ctrl_out), SAVE :: o_slab_bilg = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
670	'slab_bilg', 'flux glace - slab ponderes fsic', 'W/m2', (/ ('', i=1, 9) /))
671	TYPE(ctrl_out), SAVE :: o_slab_qflux = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
672	'slab_qflux', 'Correction flux slab', 'W/m2', (/ ('', i=1, 9) /))
673	TYPE(ctrl_out), SAVE :: o_tslab = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
674	'tslab', 'Temperature ocean slab', 'K', (/ ('', i=1, 9) /))
675	TYPE(ctrl_out), SAVE :: o_slab_tice = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
676	'slab_tice', 'Temperature banquise slab', 'K', (/ ('', i=1, 9) /))
677	TYPE(ctrl_out), SAVE :: o_slab_sic = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
678	'seaice', 'Epaisseur banquise slab', 'kg/m2', (/ ('', i=1, 9) /))
679	TYPE(ctrl_out), SAVE :: o_slab_hdiff = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
680	'slab_hdiff', 'Horizontal diffusion', 'W/m2', (/ ('', i=1, 9) /))
681	TYPE(ctrl_out), SAVE :: o_slab_ekman = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
682	'slab_ekman', 'Ekman heat transport', 'W/m2', (/ ('', i=1, 9) /))
683	TYPE(ctrl_out), SAVE :: o_ale_bl = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
684	'ale_bl', 'ALE BL', 'm2/s2', (/ ('', i=1, 9) /))
685	TYPE(ctrl_out), SAVE :: o_alp_bl = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
686	'alp_bl', 'ALP BL', 'W/m2', (/ ('', i=1, 9) /))
687	TYPE(ctrl_out), SAVE :: o_ale_wk = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
688	'ale_wk', 'ALE WK', 'm2/s2', (/ ('', i=1, 9) /))
689	TYPE(ctrl_out), SAVE :: o_alp_wk = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
690	'alp_wk', 'ALP WK', 'W/m2', (/ ('', i=1, 9) /))
691	!!!
692	!nrlmd+jyg<
693	type(ctrl_out),save :: o_dtvdf_x = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
694	'dtvdf_x', ' dtvdf off_wake','K/s', (/ ('', i=1, 9) /))
695	type(ctrl_out),save :: o_dtvdf_w = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
696	'dtvdf_w', ' dtvdf within_wake','K/s', (/ ('', i=1, 9) /))
697	type(ctrl_out),save :: o_dqvdf_x = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
698	'dqvdf_x', ' dqvdf off_wake','kg/kg/s', (/ ('', i=1, 9) /))
699	type(ctrl_out),save :: o_dqvdf_w = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
700	'dqvdf_w', ' dqvdf within_wake','kg/kg/s', (/ ('', i=1, 9) /))
701	!!
702	type(ctrl_out),save :: o_sens_x = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
703	'sens_x', 'sens off_wake', 'W/m2', (/ ('', i=1, 9) /))
704	type(ctrl_out),save :: o_sens_w = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
705	'sens_w', 'sens within_wake', 'W/m2', (/ ('', i=1, 9) /))
706	type(ctrl_out),save :: o_flat_x = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
707	'flat_x', 'flat off_wake', 'W/m2', (/ ('', i=1, 9) /))
708	type(ctrl_out),save :: o_flat_w = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
709	'flat_w', 'flat within_wake', 'W/m2', (/ ('', i=1, 9) /))
710	!!
711	type(ctrl_out),save :: o_delta_tsurf = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
712	'delta_tsurf', 'Temperature difference (w-x)', 'K', (/ ('', i=1, 9) /))
713	type(ctrl_out),save :: o_cdragh_x = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
714	'cdragh_x', 'cdragh off-wake', '', (/ ('', i=1, 9) /))
715	type(ctrl_out),save :: o_cdragh_w = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
716	'cdragh_w', 'cdragh within-wake', '', (/ ('', i=1, 9) /))
717	type(ctrl_out),save :: o_cdragm_x = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
718	'cdragm_x', 'cdragm off-wake', '', (/ ('', i=1, 9) /))
719	type(ctrl_out),save :: o_cdragm_w = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
720	'cdragm_w', 'cdrgam within-wake', '', (/ ('', i=1, 9) /))
721	type(ctrl_out),save :: o_kh = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
722	'kh', 'Kh', 'kg/s/m2', (/ ('', i=1, 9) /))
723	type(ctrl_out),save :: o_kh_x = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
724	'kh_x', 'Kh off-wake', 'kg/s/m2', (/ ('', i=1, 9) /))
725	type(ctrl_out),save :: o_kh_w = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
726	'kh_w', 'Kh within-wake', 'kg/s/m2', (/ ('', i=1, 9) /))
727	!>nrlmd+jyg
728	!!!
729	TYPE(ctrl_out), SAVE :: o_ale = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
730	'ale', 'ALE', 'm2/s2', (/ ('', i=1, 9) /))
731	TYPE(ctrl_out), SAVE :: o_alp = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
732	'alp', 'ALP', 'W/m2', (/ ('', i=1, 9) /))
733	TYPE(ctrl_out), SAVE :: o_cin = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
734	'cin', 'Convective INhibition', 'm2/s2', (/ ('', i=1, 9) /))
735	TYPE(ctrl_out), SAVE :: o_wape = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
736	'wape', '', 'm2/s2', (/ ('', i=1, 9) /))
737
738	!!! nrlmd le 10/04/2012
739
740	!-------Spectre de thermiques de type 2 au LCL
741	TYPE(ctrl_out), SAVE :: o_n2 = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11 /), &
742	'n2', 'Nombre de panaches de type 2', ' ', (/ ('', i=1, 9) /))
743	TYPE(ctrl_out), SAVE :: o_s2 = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11 /), &
744	's2', 'Surface moyenne des panaches de type 2', 'm2', (/ ('', i=1, 9) /))
745
746	!-------Déclenchement stochastique
747	TYPE(ctrl_out), SAVE :: o_proba_notrig = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11 /), &
748	'proba_notrig', &
749	'Probabilite de non-declenchement', ' ', (/ ('', i=1, 9) /))
750	TYPE(ctrl_out), SAVE :: o_random_notrig = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11 /), &
751	'random_notrig', &
752	'Tirage aleatoire de non-declenchement', ' ', (/ ('', i=1, 9) /))
753	TYPE(ctrl_out), SAVE :: o_ale_bl_stat = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11 /), &
754	'ale_bl_stat', &
755	'ALE_BL_STAT', 'm2/s2', (/ ('', i=1, 9) /))
756	TYPE(ctrl_out), SAVE :: o_ale_bl_trig = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11 /), &
757	'ale_bl_trig', &
758	'ALE_BL_STAT + Condition S>Sthreshold', 'm2/s2', (/ ('', i=1, 9) /))
759
760	!-------Fermeture statistique
761	TYPE(ctrl_out), SAVE :: o_alp_bl_det = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
762	'alp_bl_det', 'ALP_BL_DET', 'W/m2', (/ ('', i=1, 9) /))
763	TYPE(ctrl_out), SAVE :: o_alp_bl_fluct_m = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
764	'alp_bl_fluct_m', 'ALP_BL_FLUCT_M', 'W/m2', (/ ('', i=1, 9) /))
765	TYPE(ctrl_out), SAVE :: o_alp_bl_fluct_tke = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
766	'alp_bl_fluct_tke', 'ALP_BL_FLUCT_TKE', 'W/m2', (/ ('', i=1, 9) /))
767	TYPE(ctrl_out), SAVE :: o_alp_bl_conv = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
768	'alp_bl_conv', 'ALP_BL_CONV', 'W/m2', (/ ('', i=1, 9) /))
769	TYPE(ctrl_out), SAVE :: o_alp_bl_stat = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
770	'alp_bl_stat', 'ALP_BL_STAT', 'W/m2', (/ ('', i=1, 9) /))
771
772	!!! fin nrlmd le 10/04/2012
773
774	! Champs interpolles sur des niveaux de pression ??? a faire correctement
775
776	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_uSTDlevs = (/ &
777	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u850', "Zonal wind 850hPa", "m/s", &
778	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
779	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u700', "Zonal wind 700hPa", "m/s", &
780	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
781	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u500', "Zonal wind 500hPa", "m/s", &
782	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
783	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u200', "Zonal wind 200hPa", "m/s", &
784	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
785	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u100', "Zonal wind 100hPa", "m/s", &
786	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
787	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u50', "Zonal wind 50hPa", "m/s", &
788	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
789	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u10', "Zonal wind 10hPa", "m/s", &
790	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
791
792	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_vSTDlevs = (/ &
793	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v850', "Meridional wind 850hPa", "m/s", &
794	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
795	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v700', "Meridional wind 700hPa", "m/s", &
796	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
797	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v500', "Meridional wind 500hPa", "m/s", &
798	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
799	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v200', "Meridional wind 200hPa", "m/s", &
800	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
801	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v100', "Meridional wind 100hPa", "m/s", &
802	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
803	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v50', "Meridional wind 50hPa", "m/s", &
804	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
805	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v10', "Meridional wind 10hPa", "m/s", &
806	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
807
808	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_wSTDlevs = (/ &
809	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w850', "Vertical wind 850hPa", "Pa/s", &
810	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
811	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w700', "Vertical wind 700hPa", "Pa/s", &
812	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
813	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w500', "Vertical wind 500hPa", "Pa/s", &
814	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
815	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w200', "Vertical wind 200hPa", "Pa/s", &
816	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
817	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w100', "Vertical wind 100hPa", "Pa/s", &
818	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
819	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w50', "Vertical wind 50hPa", "Pa/s", &
820	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
821	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w10', "Vertical wind 10hPa", "Pa/s", &
822	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
823
824	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_tSTDlevs = (/ &
825	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t850', "Temperature 850hPa", "K", &
826	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
827	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t700', "Temperature 700hPa", "K", &
828	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
829	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t500', "Temperature 500hPa", "K", &
830	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
831	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t200', "Temperature 200hPa", "K", &
832	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
833	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t100', "Temperature 100hPa", "K", &
834	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
835	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t50', "Temperature 50hPa", "K", &
836	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
837	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t10', "Temperature 10hPa", "K", &
838	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
839
840	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_qSTDlevs = (/ &
841	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q850', "Specific humidity 850hPa", &
842	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
843	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q700', "Specific humidity 700hPa", &
844	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
845	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q500', "Specific humidity 500hPa", &
846	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
847	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q200', "Specific humidity 200hPa", &
848	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
849	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q100', "Specific humidity 100hPa", &
850	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
851	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q50', "Specific humidity 50hPa", &
852	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
853	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q10', "Specific humidity 10hPa", &
854	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
855
856	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_zSTDlevs = (/ &
857	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z850', "Geopotential height 850hPa", &
858	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
859	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z700', "Geopotential height 700hPa", &
860	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
861	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z500', "Geopotential height 500hPa", &
862	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
863	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z200', "Geopotential height 200hPa", &
864	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
865	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z100', "Geopotential height 100hPa", &
866	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
867	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z50', "Geopotential height 50hPa", &
868	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
869	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z10', "Geopotential height 10hPa", &
870	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
871
872	TYPE(ctrl_out), SAVE :: o_t_oce_sic = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
873	't_oce_sic', 'Temp mixte oce-sic', 'K', (/ ('', i=1, 9) /))
874	TYPE(ctrl_out), SAVE :: o_weakinv = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
875	'weakinv', 'Weak inversion', '-', (/ ('', i=1, 9) /))
876	TYPE(ctrl_out), SAVE :: o_dthmin = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
877	'dthmin', 'dTheta mini', 'K/m', (/ ('', i=1, 9) /))
878
879	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_u10_srf = (/ &
880	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'u10_ter', "", "", (/ ('', i=1, 9) /)), &
881	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'u10_lic', "", "", (/ ('', i=1, 9) /)), &
882	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'u10_oce', "", "", (/ ('', i=1, 9) /)), &
883	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'u10_sic', "", "", (/ ('', i=1, 9) /)) /)
884
885	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_v10_srf = (/ &
886	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'v10_ter', "", "", (/ ('', i=1, 9) /)), &
887	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'v10_lic', "", "", (/ ('', i=1, 9) /)), &
888	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'v10_oce', "", "", (/ ('', i=1, 9) /)), &
889	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'v10_sic', "", "", (/ ('', i=1, 9) /)) /)
890
891	TYPE(ctrl_out), SAVE :: o_cldtau = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
892	'cldtau', 'Cloud optical thickness', '1', (/ ('', i=1, 9) /))
893	TYPE(ctrl_out), SAVE :: o_cldemi = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
894	'cldemi', 'Cloud optical emissivity', '1', (/ ('', i=1, 9) /))
895	TYPE(ctrl_out), SAVE :: o_rh2m = ctrl_out((/ 5, 5, 10, 10, 10, 10, 11, 11, 11 /), &
896	'rh2m', 'Relative humidity at 2m', '%', (/ ('', i=1, 9) /))
897	TYPE(ctrl_out), SAVE :: o_rh2m_min = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
898	'rh2m_min', 'Min Relative humidity at 2m', '%', &
899	(/ 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)' /))
900	TYPE(ctrl_out), SAVE :: o_rh2m_max = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
901	'rh2m_max', 'Max Relative humidity at 2m', '%', &
902	(/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
903	't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
904	TYPE(ctrl_out), SAVE :: o_qsat2m = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
905	'qsat2m', 'Saturant humidity at 2m', '%', (/ ('', i=1, 9) /))
906	TYPE(ctrl_out), SAVE :: o_tpot = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
907	'tpot', 'Surface air potential temperature', 'K', (/ ('', i=1, 9) /))
908	TYPE(ctrl_out), SAVE :: o_tpote = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
909	'tpote', &
910	'Surface air equivalent potential temperature', 'K', (/ ('', i=1, 9) /))
911	TYPE(ctrl_out), SAVE :: o_tke = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
912	'tke ', 'TKE', 'm2/s2', (/ ('', i=1, 9) /))
913	TYPE(ctrl_out), SAVE :: o_tke_max = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
914	'tke_max', 'TKE max', 'm2/s2', &
915	(/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
916	't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
917
918	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tke_srf = (/ &
919	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_ter', &
920	"Max Turb. Kinetic Energy "//clnsurf(1),"m2/s2", (/ ('', i=1, 9) /)), &
921	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_lic', &
922	"Max Turb. Kinetic Energy "//clnsurf(2),"m2/s2", (/ ('', i=1, 9) /)), &
923	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_oce', &
924	"Max Turb. Kinetic Energy "//clnsurf(3),"m2/s2", (/ ('', i=1, 9) /)), &
925	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_sic', &
926	"Max Turb. Kinetic Energy "//clnsurf(4),"m2/s2", (/ ('', i=1, 9) /)) /)
927
928	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_l_mixmin = (/ &
929	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'l_mixmin_ter', &
930	"PBL mixing length "//clnsurf(1),"m", (/ ('', i=1, 9) /)), &
931	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'l_mixmin_lic', &
932	"PBL mixing length "//clnsurf(2),"m", (/ ('', i=1, 9) /)), &
933	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'l_mixmin_oce', &
934	"PBL mixing length "//clnsurf(3),"m", (/ ('', i=1, 9) /)), &
935	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'l_mixmin_sic', &
936	"PBL mixing length "//clnsurf(4),"m", (/ ('', i=1, 9) /)) /)
937
938	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_l_mix = (/ &
939	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'l_mix_ter', &
940	"min PBL mixing length "//clnsurf(1),"m", (/ ('', i=1, 9) /)), &
941	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'l_mix_lic', &
942	"min PBL mixing length "//clnsurf(2),"m", (/ ('', i=1, 9) /)), &
943	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'l_mix_oce', &
944	"min PBL mixing length "//clnsurf(3),"m", (/ ('', i=1, 9) /)), &
945	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'l_mix_sic', &
946	"min PBL mixing length "//clnsurf(4),"m", (/ ('', i=1, 9) /)) /)
947
948	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tke_max_srf = (/ &
949	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_max_ter', &
950	"Max Turb. Kinetic Energy "//clnsurf(1),"-", &
951	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
952	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /)), &
953	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_max_lic', &
954	"Max Turb. Kinetic Energy "//clnsurf(2),"-", &
955	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
956	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /)), &
957	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_max_oce', &
958	"Max Turb. Kinetic Energy "//clnsurf(3),"-", &
959	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
960	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /)), &
961	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_max_sic', &
962	"Max Turb. Kinetic Energy "//clnsurf(4),"-", &
963	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
964	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /)) /)
965
966	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_dltpbltke_srf = (/ &
967	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'dltpbltke_ter', &
968	"TKE difference (w - x) "//clnsurf(1),"-", (/ ('', i=1, 9) /)), &
969	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'dltpbltke_lic', &
970	"TKE difference (w - x) "//clnsurf(2),"-", (/ ('', i=1, 9) /)), &
971	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'dltpbltke_oce', &
972	"TKE difference (w - x) "//clnsurf(3),"-", (/ ('', i=1, 9) /)), &
973	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'dltpbltke_sic', &
974	"TKE difference (w - x) "//clnsurf(4),"-", (/ ('', i=1, 9) /)) /)
975
976	TYPE(ctrl_out), SAVE :: o_kz = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
977	'kz', 'Kz melange', 'm2/s', (/ ('', i=1, 9) /))
978	TYPE(ctrl_out), SAVE :: o_kz_max = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
979	'kz_max', 'Kz melange max', 'm2/s', &
980	(/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
981	't_max(X)', "t_max(X)", "t_max(X)", "t_max(X)" /))
982	TYPE(ctrl_out), SAVE :: o_SWnetOR = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
983	'SWnetOR', 'Sfce net SW radiation OR', 'W/m2', (/ ('', i=1, 9) /))
984	TYPE(ctrl_out), SAVE :: o_SWdownOR = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
985	'SWdownOR', 'Sfce incident SW radiation OR', 'W/m2', (/ ('', i=1, 9) /))
986	TYPE(ctrl_out), SAVE :: o_LWdownOR = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
987	'LWdownOR', 'Sfce incident LW radiation OR', 'W/m2', (/ ('', i=1, 9) /))
988	TYPE(ctrl_out), SAVE :: o_snowl = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
989	'snowl', 'Solid Large-scale Precip.', 'kg/(m2*s)', (/ ('', i=1, 9) /))
990	TYPE(ctrl_out), SAVE :: o_cape_max = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
991	'cape_max', 'CAPE max.', 'J/kg', &
992	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
993	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /))
994	TYPE(ctrl_out), SAVE :: o_solldown = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
995	'solldown', 'Down. IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
996	TYPE(ctrl_out), SAVE :: o_dtsvdfo = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
997	'dtsvdfo', 'Boundary-layer dTs(o)', 'K/s', (/ ('', i=1, 9) /))
998	TYPE(ctrl_out), SAVE :: o_dtsvdft = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
999	'dtsvdft', 'Boundary-layer dTs(t)', 'K/s', (/ ('', i=1, 9) /))
1000	TYPE(ctrl_out), SAVE :: o_dtsvdfg = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1001	'dtsvdfg', 'Boundary-layer dTs(g)', 'K/s', (/ ('', i=1, 9) /))
1002	TYPE(ctrl_out), SAVE :: o_dtsvdfi = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1003	'dtsvdfi', 'Boundary-layer dTs(g)', 'K/s', (/ ('', i=1, 9) /))
1004	TYPE(ctrl_out), SAVE :: o_z0m = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1005	'z0m', 'roughness length, momentum', '-', (/ ('', i=1, 9) /))
1006	TYPE(ctrl_out), SAVE :: o_z0h = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1007	'z0h', 'roughness length, enthalpy', '-', (/ ('', i=1, 9) /))
1008	TYPE(ctrl_out), SAVE :: o_topswad = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1009	'topswad', 'ADE at TOA', 'W/m2', (/ ('', i=1, 9) /))
1010	TYPE(ctrl_out), SAVE :: o_topswad0 = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1011	'topswad0', 'ADE clear-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
1012	TYPE(ctrl_out), SAVE :: o_topswai = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1013	'topswai', 'AIE at TOA', 'W/m2', (/ ('', i=1, 9) /))
1014	TYPE(ctrl_out), SAVE :: o_solswad = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1015	'solswad', 'ADE at SRF', 'W/m2', (/ ('', i=1, 9) /))
1016	TYPE(ctrl_out), SAVE :: o_solswad0 = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1017	'solswad0', 'ADE clear-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
1018	TYPE(ctrl_out), SAVE :: o_solswai = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1019	'solswai', 'AIE at SFR', 'W/m2', (/ ('', i=1, 9) /))
1020	TYPE(ctrl_out), SAVE :: o_toplwad = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1021	'toplwad', 'LW-ADE at TOA', 'W/m2', (/ ('', i=1, 9) /))
1022	TYPE(ctrl_out), SAVE :: o_toplwad0 = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1023	'toplwad0', 'LW-ADE clear-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
1024	TYPE(ctrl_out), SAVE :: o_toplwai = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1025	'toplwai', 'LW-AIE at TOA', 'W/m2', (/ ('', i=1, 9) /))
1026	TYPE(ctrl_out), SAVE :: o_sollwad = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1027	'sollwad', 'LW-ADE at SRF', 'W/m2', (/ ('', i=1, 9) /))
1028	TYPE(ctrl_out), SAVE :: o_sollwad0 = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1029	'sollwad0', 'LW-ADE clear-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
1030	TYPE(ctrl_out), SAVE :: o_sollwai = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1031	'sollwai', 'LW-AIE at SFR', 'W/m2', (/ ('', i=1, 9) /))
1032
1033
1034	TYPE(ctrl_out),SAVE,DIMENSION(naero_tot) :: o_tausumaero = &
1035	(/ ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(1), &
1036	"Aerosol Optical depth at 550 nm "//name_aero_tau(1),"1", (/ ('', i=1, 9) /)), &
1037	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(2), &
1038	"Aerosol Optical depth at 550 nm "//name_aero_tau(2),"2", (/ ('', i=1, 9) /)), &
1039	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(3), &
1040	"Aerosol Optical depth at 550 nm "//name_aero_tau(3),"3", (/ ('', i=1, 9) /)), &
1041	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(4), &
1042	"Aerosol Optical depth at 550 nm "//name_aero_tau(4),"4", (/ ('', i=1, 9) /)), &
1043	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(5), &
1044	"Aerosol Optical depth at 550 nm "//name_aero_tau(5),"5", (/ ('', i=1, 9) /)), &
1045	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(6), &
1046	"Aerosol Optical depth at 550 nm "//name_aero_tau(6),"6", (/ ('', i=1, 9) /)), &
1047	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(7), &
1048	"Aerosol Optical depth at 550 nm "//name_aero_tau(7),"7", (/ ('', i=1, 9) /)), &
1049	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(8), &
1050	"Aerosol Optical depth at 550 nm "//name_aero_tau(8),"8", (/ ('', i=1, 9) /)), &
1051	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(9), &
1052	"Aerosol Optical depth at 550 nm "//name_aero_tau(9),"9", (/ ('', i=1, 9) /)), &
1053	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(10), &
1054	"Aerosol Optical depth at 550 nm "//name_aero_tau(10),"10", (/ ('', i=1, 9) /)), &
1055	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(11), &
1056	"Aerosol Optical depth at 550 nm "//name_aero_tau(11),"11", (/ ('', i=1, 9) /)), &
1057	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(12), &
1058	"Aerosol Optical depth at 550 nm "//name_aero_tau(12),"12", (/ ('', i=1, 9) /)), &
1059	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(13), &
1060	"Aerosol Optical depth at 550 nm "//name_aero_tau(13),"13", (/ ('', i=1, 9) /)), &
1061	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_'//name_aero_tau(14), &
1062	"Aerosol Optical depth at 550 nm "//name_aero_tau(14),"14", (/ ('', i=1, 9) /)) /)
1063	!
1064	TYPE(ctrl_out), SAVE :: o_tausumaero_lw = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1065	'OD_10um_STRAT', 'Stratospheric Aerosol Optical depth at 10 um ', '1', (/ ('', i=1, 9) /))
1066	!
1067	TYPE(ctrl_out), SAVE :: o_od550aer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1068	'od550aer', 'Total aerosol optical depth at 550nm', '-', (/ ('', i=1, 9) /))
1069	TYPE(ctrl_out), SAVE :: o_od865aer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1070	'od865aer', 'Total aerosol optical depth at 870nm', '-', (/ ('', i=1, 9) /))
1071	TYPE(ctrl_out), SAVE :: o_absvisaer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1072	'absvisaer', 'Absorption aerosol visible optical depth', '-', (/ ('', i=1, 9) /))
1073	TYPE(ctrl_out), SAVE :: o_od550lt1aer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1074	'od550lt1aer', 'Fine mode optical depth', '-', (/ ('', i=1, 9) /))
1075	TYPE(ctrl_out), SAVE :: o_sconcso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1076	'sconcso4', 'Surface Concentration of Sulfate ', 'kg/m3', (/ ('', i=1, 9) /))
1077	TYPE(ctrl_out), SAVE :: o_sconcno3 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1078	'sconcno3', 'Surface Concentration of Nitrate ', 'kg/m3', (/ ('', i=1, 9) /))
1079	TYPE(ctrl_out), SAVE :: o_sconcoa = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1080	'sconcoa', 'Surface Concentration of Organic Aerosol ', 'kg/m3', (/ ('', i=1, 9) /))
1081	TYPE(ctrl_out), SAVE :: o_sconcbc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1082	'sconcbc', 'Surface Concentration of Black Carbon ', 'kg/m3', (/ ('', i=1, 9) /))
1083	TYPE(ctrl_out), SAVE :: o_sconcss = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1084	'sconcss', 'Surface Concentration of Sea Salt ', 'kg/m3', (/ ('', i=1, 9) /))
1085	TYPE(ctrl_out), SAVE :: o_sconcdust = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1086	'sconcdust', 'Surface Concentration of Dust ', 'kg/m3', (/ ('', i=1, 9) /))
1087	TYPE(ctrl_out), SAVE :: o_concso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1088	'concso4', 'Concentration of Sulfate ', 'kg/m3', (/ ('', i=1, 9) /))
1089	TYPE(ctrl_out), SAVE :: o_concno3 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1090	'concno3', 'Concentration of Nitrate ', 'kg/m3', (/ ('', i=1, 9) /))
1091	TYPE(ctrl_out), SAVE :: o_concoa = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1092	'concoa', 'Concentration of Organic Aerosol ', 'kg/m3', (/ ('', i=1, 9) /))
1093	TYPE(ctrl_out), SAVE :: o_concbc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1094	'concbc', 'Concentration of Black Carbon ', 'kg/m3', (/ ('', i=1, 9) /))
1095	TYPE(ctrl_out), SAVE :: o_concss = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1096	'concss', 'Concentration of Sea Salt ', 'kg/m3', (/ ('', i=1, 9) /))
1097	TYPE(ctrl_out), SAVE :: o_concdust = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1098	'concdust', 'Concentration of Dust ', 'kg/m3', (/ ('', i=1, 9) /))
1099	TYPE(ctrl_out), SAVE :: o_loadso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1100	'loadso4', 'Column Load of Sulfate ', 'kg/m2', (/ ('', i=1, 9) /))
1101	TYPE(ctrl_out), SAVE :: o_loadoa = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1102	'loadoa', 'Column Load of Organic Aerosol ', 'kg/m2', (/ ('', i=1, 9) /))
1103	TYPE(ctrl_out), SAVE :: o_loadbc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1104	'loadbc', 'Column Load of Black Carbon ', 'kg/m2', (/ ('', i=1, 9) /))
1105	TYPE(ctrl_out), SAVE :: o_loadss = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1106	'loadss', 'Column Load of Sea Salt ', 'kg/m2', (/ ('', i=1, 9) /))
1107	TYPE(ctrl_out), SAVE :: o_loaddust = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1108	'loaddust', 'Column Load of Dust ', 'kg/m2', (/ ('', i=1, 9) /))
1109	TYPE(ctrl_out), SAVE :: o_swtoaas_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1110	'swtoaas_nat', 'Natural aerosol radiative forcing all-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
1111	TYPE(ctrl_out), SAVE :: o_swsrfas_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1112	'swsrfas_nat', 'Natural aerosol radiative forcing all-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
1113	TYPE(ctrl_out), SAVE :: o_swtoacs_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1114	'swtoacs_nat', 'Natural aerosol radiative forcing clear-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
1115	TYPE(ctrl_out), SAVE :: o_swsrfcs_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1116	'swsrfcs_nat', 'Natural aerosol radiative forcing clear-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
1117	TYPE(ctrl_out), SAVE :: o_swtoaas_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1118	'swtoaas_ant', 'Anthropogenic aerosol radiative forcing all-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
1119	TYPE(ctrl_out), SAVE :: o_swsrfas_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1120	'swsrfas_ant', 'Anthropogenic aerosol radiative forcing all-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
1121	TYPE(ctrl_out), SAVE :: o_swtoacs_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1122	'swtoacs_ant', 'Anthropogenic aerosol radiative forcing clear-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
1123	TYPE(ctrl_out), SAVE :: o_swsrfcs_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1124	'swsrfcs_ant', 'Anthropogenic aerosol radiative forcing clear-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
1125	TYPE(ctrl_out), SAVE :: o_swtoacf_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1126	'swtoacf_nat', 'Natural aerosol impact on cloud radiative forcing at TOA', 'W/m2', (/ ('', i=1, 9) /))
1127	TYPE(ctrl_out), SAVE :: o_swsrfcf_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1128	'swsrfcf_nat', 'Natural aerosol impact on cloud radiative forcing at SRF', 'W/m2', (/ ('', i=1, 9) /))
1129	TYPE(ctrl_out), SAVE :: o_swtoacf_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1130	'swtoacf_ant', 'Anthropogenic aerosol impact on cloud radiative forcing at TOA', 'W/m2', (/ ('', i=1, 9) /))
1131	TYPE(ctrl_out), SAVE :: o_swsrfcf_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1132	'swsrfcf_ant', 'Anthropogenic aerosol impact on cloud radiative forcing at SRF', 'W/m2', (/ ('', i=1, 9) /))
1133	TYPE(ctrl_out), SAVE :: o_swtoacf_zero = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1134	'swtoacf_zero', 'Cloud radiative forcing (allsky-clearsky fluxes) at TOA', 'W/m2', (/ ('', i=1, 9) /))
1135	TYPE(ctrl_out), SAVE :: o_swsrfcf_zero = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1136	'swsrfcf_zero', 'Cloud radiative forcing (allsky-clearsky fluxes) at SRF', 'W/m2', (/ ('', i=1, 9) /))
1137	TYPE(ctrl_out), SAVE :: o_cldncl = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1138	'cldncl', 'CDNC at top of liquid water cloud', 'm-3', (/ ('', i=1, 9) /))
1139	TYPE(ctrl_out), SAVE :: o_reffclwtop = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1140	'reffclwtop', 'Droplet effective radius at top of liquid water cloud', 'm', (/ ('', i=1, 9) /))
1141	TYPE(ctrl_out), SAVE :: o_cldnvi = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1142	'cldnvi', 'Column Integrated Cloud Droplet Number', 'm-2', (/ ('', i=1, 9) /))
1143	TYPE(ctrl_out), SAVE :: o_lcc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1144	'lcc', 'Cloud liquid fraction at top of cloud', '1', (/ ('', i=1, 9) /))
1145
1146
1147	!!!!!!!!!!!!!!!!!!!!!! 3D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1148	TYPE(ctrl_out), SAVE :: o_ec550aer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1149	'ec550aer', 'Extinction at 550nm', 'm^-1', (/ ('', i=1, 9) /))
1150	TYPE(ctrl_out), SAVE :: o_lwcon = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1151	'lwcon', 'Cloud liquid water content', 'kg/kg', (/ ('', i=1, 9) /))
1152	TYPE(ctrl_out), SAVE :: o_iwcon = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1153	'iwcon', 'Cloud ice water content', 'kg/kg', (/ ('', i=1, 9) /))
1154	TYPE(ctrl_out), SAVE :: o_temp = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11 /), &
1155	'temp', 'Air temperature', 'K', (/ ('', i=1, 9) /))
1156	TYPE(ctrl_out), SAVE :: o_theta = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11 /), &
1157	'theta', 'Potential air temperature', 'K', (/ ('', i=1, 9) /))
1158	TYPE(ctrl_out), SAVE :: o_ovap = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11 /), &
1159	'ovap', 'Specific humidity', 'kg/kg', (/ ('', i=1, 9) /))
1160	TYPE(ctrl_out), SAVE :: o_ovapinit = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1161	'ovapinit', 'Specific humidity (begin of timestep)', 'kg/kg', (/ ('', i=1, 9) /))
1162	TYPE(ctrl_out), SAVE :: o_oliq = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11 /), &
1163	'oliq', 'Liquid water', 'kg/kg', (/ ('', i=1, 9) /))
1164	TYPE(ctrl_out), SAVE :: o_ocond = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11 /), &
1165	'ocond', 'Condensed water', 'kg/kg', (/ ('', i=1, 9) /))
1166	TYPE(ctrl_out), SAVE :: o_wvapp = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1167	'wvapp', '', '', (/ ('', i=1, 9) /))
1168	TYPE(ctrl_out), SAVE :: o_geop = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
1169	'geop', 'Geopotential height', 'm2/s2', (/ ('', i=1, 9) /))
1170	TYPE(ctrl_out), SAVE :: o_vitu = ctrl_out((/ 2, 3, 4, 6, 10, 10, 11, 11, 11 /), &
1171	'vitu', 'Zonal wind', 'm/s', (/ ('', i=1, 9) /))
1172	TYPE(ctrl_out), SAVE :: o_vitv = ctrl_out((/ 2, 3, 4, 6, 10, 10, 11, 11, 11 /), &
1173	'vitv', 'Meridional wind', 'm/s', (/ ('', i=1, 9) /))
1174	TYPE(ctrl_out), SAVE :: o_vitw = ctrl_out((/ 2, 3, 10, 6, 10, 10, 11, 11, 11 /), &
1175	'vitw', 'Vertical wind', 'Pa/s', (/ ('', i=1, 9) /))
1176	TYPE(ctrl_out), SAVE :: o_pres = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
1177	'pres', 'Air pressure', 'Pa', (/ ('', i=1, 9) /))
1178	TYPE(ctrl_out), SAVE :: o_paprs = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
1179	'paprs', 'Air pressure Inter-Couches', 'Pa', (/ ('', i=1, 9) /))
1180	TYPE(ctrl_out), SAVE :: o_mass = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
1181	'mass', 'Masse Couches', 'kg/m2', (/ ('', i=1, 9) /))
1182	TYPE(ctrl_out), SAVE :: o_zfull = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
1183	'zfull', 'Altitude of full pressure levels', 'm', (/ ('', i=1, 9) /))
1184	TYPE(ctrl_out), SAVE :: o_zhalf = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
1185	'zhalf', 'Altitude of half pressure levels', 'm', (/ ('', i=1, 9) /))
1186	TYPE(ctrl_out), SAVE :: o_rneb = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1187	'rneb', 'Cloud fraction', '-', (/ ('', i=1, 9) /))
1188	TYPE(ctrl_out), SAVE :: o_rnebjn = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11,11 /), &
1189	'rnebjn', 'Cloud fraction in day', '-', (/ ('', i=1, 9) /))
1190	TYPE(ctrl_out), SAVE :: o_rnebcon = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1191	'rnebcon', 'Convective Cloud Fraction', '-', (/ ('', i=1, 9) /))
1192	TYPE(ctrl_out), SAVE :: o_rnebls = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1193	'rnebls', 'LS Cloud fraction', '-', (/ ('', i=1, 9) /))
1194	TYPE(ctrl_out), SAVE :: o_rhum = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1195	'rhum', 'Relative humidity', '-', (/ ('', i=1, 9) /))
1196	TYPE(ctrl_out), SAVE :: o_ozone = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1197	'ozone', 'Ozone mole fraction', '-', (/ ('', i=1, 9) /))
1198	TYPE(ctrl_out), SAVE :: o_ozone_light = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1199	'ozone_daylight', 'Daylight ozone mole fraction', '-', (/ ('', i=1, 9) /))
1200	TYPE(ctrl_out), SAVE :: o_upwd = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1201	'upwd', 'saturated updraft', 'kg/m2/s', (/ ('', i=1, 9) /))
1202	TYPE(ctrl_out), SAVE :: o_epmax_diag = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1203	'epmax', 'epmax en fn cape', 'su', (/ ('', i=1, 9) /))
1204	TYPE(ctrl_out), SAVE :: o_ep = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1205	'ep', 'ep', 'su', (/ ('', i=1, 9) /))
1206	TYPE(ctrl_out), SAVE :: o_duphy = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1207	'duphy', 'Physics du', 'm/s2', (/ ('', i=1, 9) /))
1208	TYPE(ctrl_out), SAVE :: o_dtphy = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1209	'dtphy', 'Physics dT', 'K/s', (/ ('', i=1, 9) /))
1210	TYPE(ctrl_out), SAVE :: o_dqphy = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1211	'dqphy', 'Physics dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1212	TYPE(ctrl_out), SAVE :: o_dqphy2d = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1213	'dqphy2d', 'Physics dQ', '(kg/m2)/s', (/ ('', i=1, 9) /))
1214	TYPE(ctrl_out), SAVE :: o_dqlphy = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1215	'dqlphy', 'Physics dQL', '(kg/kg)/s', (/ ('', i=1, 9) /))
1216	TYPE(ctrl_out), SAVE :: o_dqlphy2d = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1217	'dqlphy2d', 'Physics dQL', '(kg/m2)/s', (/ ('', i=1, 9) /))
1218	TYPE(ctrl_out), SAVE :: o_dqsphy = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1219	'dqsphy', 'Physics dQS', '(kg/kg)/s', (/ ('', i=1, 9) /))
1220	TYPE(ctrl_out), SAVE :: o_dqsphy2d = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1221	'dqsphy2d', 'Physics dQS', '(kg/m2)/s', (/ ('', i=1, 9) /))
1222	TYPE(ctrl_out), SAVE :: o_pr_con_l = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1223	'pr_con_l', 'Convective precipitation lic', ' ', (/ ('', i=1, 9) /))
1224	TYPE(ctrl_out), SAVE :: o_pr_con_i = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1225	'pr_con_i', 'Convective precipitation ice', ' ', (/ ('', i=1, 9) /))
1226	TYPE(ctrl_out), SAVE :: o_pr_lsc_l = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1227	'pr_lsc_l', 'Large scale precipitation lic', ' ', (/ ('', i=1, 9) /))
1228	TYPE(ctrl_out), SAVE :: o_pr_lsc_i = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1229	'pr_lsc_i', 'Large scale precipitation ice', ' ', (/ ('', i=1, 9) /))
1230	TYPE(ctrl_out), SAVE :: o_re = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1231	're', 'Cloud droplet effective radius', 'um', (/ ('', i=1, 9) /))
1232	TYPE(ctrl_out), SAVE :: o_fl = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1233	'fl', 'Denominator of Cloud droplet effective radius', ' ', (/ ('', i=1, 9) /))
1234	TYPE(ctrl_out), SAVE :: o_scdnc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1235	'scdnc', 'Cloud droplet number concentration', 'm-3', (/ ('', i=1, 9) /))
1236	TYPE(ctrl_out), SAVE :: o_reffclws = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1237	'reffclws', 'Stratiform Cloud Droplet Effective Radius (aerosol diags.)', 'm', (/ ('', i=1, 9) /))
1238	TYPE(ctrl_out), SAVE :: o_reffclwc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1239	'reffclwc', 'Convective Cloud Droplet Effective Radius (aerosol diags.)', 'm', (/ ('', i=1, 9) /))
1240	TYPE(ctrl_out), SAVE :: o_lcc3d = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1241	'lcc3d', 'Cloud liquid fraction', '1', (/ ('', i=1, 9) /))
1242	TYPE(ctrl_out), SAVE :: o_lcc3dcon = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1243	'lcc3dcon', 'Convective cloud liquid fraction', '1', (/ ('', i=1, 9) /))
1244	TYPE(ctrl_out), SAVE :: o_lcc3dstra = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1245	'lcc3dstra', 'Stratiform cloud liquid fraction', '1', (/ ('', i=1, 9) /))
1246	TYPE(ctrl_out), SAVE :: o_stratomask = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
1247	'stratomask', 'Stratospheric fraction', '1', (/ ('', i=1, 9) /))
1248	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1249
1250	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_albe_srf = (/ &
1251	ctrl_out((/ 3, 7, 10, 7, 10, 10, 11, 11, 11 /),'albe_ter', "Albedo VIS surf. "//clnsurf(1),"-", (/ ('', i=1, 9) /)), &
1252	ctrl_out((/ 3, 7, 10, 7, 10, 10, 11, 11, 11 /),'albe_lic', "Albedo VIS surf. "//clnsurf(2),"-", (/ ('', i=1, 9) /)), &
1253	ctrl_out((/ 3, 7, 10, 7, 10, 10, 11, 11, 11 /),'albe_oce', "Albedo VIS surf. "//clnsurf(3),"-", (/ ('', i=1, 9) /)), &
1254	ctrl_out((/ 3, 7, 10, 7, 10, 10, 11, 11, 11 /),'albe_sic', "Albedo VIS surf. "//clnsurf(4),"-", (/ ('', i=1, 9) /)) /)
1255
1256	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_ages_srf = (/ &
1257	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'ages_ter', "Snow age", "day", (/ ('', i=1, 9) /)), &
1258	ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /),'ages_lic', "Snow age", "day", (/ ('', i=1, 9) /)), &
1259	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'ages_oce',"Snow age", "day", (/ ('', i=1, 9) /)), &
1260	ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /),'ages_sic',"Snow age", "day", (/ ('', i=1, 9) /)) /)
1261
1262	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_snow_srf = (/ &
1263	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'snow_ter', "Snow", "kg/m2", (/ ('', i=1, 9) /)), &
1264	ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /),'snow_lic', "Snow", "kg/m2", (/ ('', i=1, 9) /)), &
1265	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'snow_oce',"Snow", "kg/m2", (/ ('', i=1, 9) /)), &
1266	ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /),'snow_sic',"Snow", "kg/m2", (/ ('', i=1, 9) /)) /)
1267
1268	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_z0m_srf = (/ &
1269	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'z0m_ter', "Surface roughness "//clnsurf(1),"m", (/ ('', i=1, 9) /)), &
1270	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'z0m_lic', "Surface roughness "//clnsurf(2),"m", (/ ('', i=1, 9) /)), &
1271	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'z0m_oce', "Surface roughness "//clnsurf(3),"m", (/ ('', i=1, 9) /)), &
1272	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'z0m_sic', "Surface roughness "//clnsurf(4),"m", (/ ('', i=1, 9) /)) /)
1273
1274	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_z0h_srf = (/ &
1275	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'z0h_ter', "Surface roughness "//clnsurf(1),"m", (/ ('', i=1, 9) /)), &
1276	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'z0h_lic', "Surface roughness "//clnsurf(2),"m", (/ ('', i=1, 9) /)), &
1277	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'z0h_oce', "Surface roughness "//clnsurf(3),"m", (/ ('', i=1, 9) /)), &
1278	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'z0h_sic', "Surface roughness "//clnsurf(4),"m", (/ ('', i=1, 9) /)) /)
1279
1280	TYPE(ctrl_out), SAVE :: o_alb1 = ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1281	'alb1', 'Surface VIS albedo', '-', (/ ('', i=1, 9) /))
1282	TYPE(ctrl_out), SAVE :: o_alb2 = ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1283	'alb2', 'Surface Near IR albedo', '-', (/ ('', i=1, 9) /))
1284	TYPE(ctrl_out), SAVE :: o_clwcon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1285	'clwcon', 'Convective Cloud Liquid water content', 'kg/kg', (/ ('', i=1, 9) /))
1286	TYPE(ctrl_out), SAVE :: o_Ma = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1287	'Ma', 'undilute adiab updraft', 'kg/m2/s', (/ ('', i=1, 9) /))
1288	TYPE(ctrl_out), SAVE :: o_dnwd = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1289	'dnwd', 'saturated downdraft', 'kg/m2/s', (/ ('', i=1, 9) /))
1290	TYPE(ctrl_out), SAVE :: o_dnwd0 = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1291	'dnwd0', 'unsat. downdraft', 'kg/m2/s', (/ ('', i=1, 9) /))
1292	TYPE(ctrl_out), SAVE :: o_mc = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1293	'mc', 'Convective mass flux', 'kg/m2/s', (/ ('', i=1, 9) /))
1294	TYPE(ctrl_out), SAVE :: o_ftime_con = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1295	'ftime_con', 'Fraction of time convection Occurs', ' ', &
1296	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /))
1297	TYPE(ctrl_out), SAVE :: o_dtdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1298	'dtdyn', 'Dynamics dT', 'K/s', (/ ('', i=1, 9) /))
1299	TYPE(ctrl_out), SAVE :: o_dqdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1300	'dqdyn', 'Dynamics dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1301	TYPE(ctrl_out), SAVE :: o_dqdyn2d = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1302	'dqdyn2d', 'Dynamics dQ', '(kg/m2)/s', (/ ('', i=1, 9) /))
1303	TYPE(ctrl_out), SAVE :: o_dqldyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1304	'dqldyn', 'Dynamics dQL', '(kg/kg)/s', (/ ('', i=1, 9) /))
1305	TYPE(ctrl_out), SAVE :: o_dqldyn2d = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1306	'dqldyn2d', 'Dynamics dQL', '(kg/m2)/s', (/ ('', i=1, 9) /))
1307	TYPE(ctrl_out), SAVE :: o_dqsdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1308	'dqsdyn', 'Dynamics dQS', '(kg/kg)/s', (/ ('', i=1, 9) /))
1309	TYPE(ctrl_out), SAVE :: o_dqsdyn2d = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1310	'dqsdyn2d', 'Dynamics dQS', '(kg/m2)/s', (/ ('', i=1, 9) /))
1311	TYPE(ctrl_out), SAVE :: o_dudyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1312	'dudyn', 'Dynamics dU', 'm/s2', (/ ('', i=1, 9) /))
1313	TYPE(ctrl_out), SAVE :: o_dvdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1314	'dvdyn', 'Dynamics dV', 'm/s2', (/ ('', i=1, 9) /))
1315	TYPE(ctrl_out), SAVE :: o_dtcon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1316	'dtcon', 'Convection dT', 'K/s', (/ ('', i=1, 9) /))
1317	TYPE(ctrl_out), SAVE :: o_ducon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1318	'ducon', 'Convection du', 'm/s2', (/ ('', i=1, 9) /))
1319	TYPE(ctrl_out), SAVE :: o_dvcon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1320	'dvcon', 'Convection dv', 'm/s2', (/ ('', i=1, 9) /))
1321	TYPE(ctrl_out), SAVE :: o_dqcon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1322	'dqcon', 'Convection dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1323	TYPE(ctrl_out), SAVE :: o_dqcon2d = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1324	'dqcon2d', 'Convection dQ', '(kg/m2)/s', (/ ('', i=1, 9) /))
1325	TYPE(ctrl_out), SAVE :: o_dtwak = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1326	'dtwak', 'Wake dT', 'K/s', (/ ('', i=1, 9) /))
1327	TYPE(ctrl_out), SAVE :: o_dqwak = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1328	'dqwak', 'Wake dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1329	TYPE(ctrl_out), SAVE :: o_dqwak2d = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1330	'dqwak2d', 'Wake dQ', '(kg/m2)/s', (/ ('', i=1, 9) /))
1331	TYPE(ctrl_out), SAVE :: o_wake_h = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1332	'wake_h', 'wake_h', '-', (/ ('', i=1, 9) /))
1333	TYPE(ctrl_out), SAVE :: o_wake_s = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1334	'wake_s', 'wake_s', '-', (/ ('', i=1, 9) /))
1335	TYPE(ctrl_out), SAVE :: o_wake_deltat = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1336	'wake_deltat', 'wake_deltat', ' ', (/ ('', i=1, 9) /))
1337	TYPE(ctrl_out), SAVE :: o_wake_deltaq = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1338	'wake_deltaq', 'wake_deltaq', ' ', (/ ('', i=1, 9) /))
1339	TYPE(ctrl_out), SAVE :: o_wake_omg = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1340	'wake_omg', 'wake_omg', ' ', (/ ('', i=1, 9) /))
1341	TYPE(ctrl_out), SAVE :: o_wdtrainA = ctrl_out((/ 4, 5, 10, 4, 10, 10, 11, 11, 110 /), &
1342	'wdtrainA', 'precipitation from AA', '-', (/ ('', i=1, 9) /))
1343	TYPE(ctrl_out), SAVE :: o_wdtrainM = ctrl_out((/ 4, 5, 10, 4, 10, 10, 11, 11, 110 /), &
1344	'wdtrainM', 'precipitation from mixture', '-', (/ ('', i=1, 9) /))
1345	TYPE(ctrl_out), SAVE :: o_Vprecip = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1346	'Vprecip', 'precipitation vertical profile', '-', (/ ('', i=1, 9) /))
1347	TYPE(ctrl_out), SAVE :: o_ftd = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1348	'ftd', 'tend temp due aux descentes precip', '-', (/ ('', i=1, 9) /))
1349	TYPE(ctrl_out), SAVE :: o_fqd = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
1350	'fqd', 'tend vap eau due aux descentes precip', '-', (/ ('', i=1, 9) /))
1351	TYPE(ctrl_out), SAVE :: o_dtlsc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1352	'dtlsc', 'Condensation dT', 'K/s', (/ ('', i=1, 9) /))
1353	TYPE(ctrl_out), SAVE :: o_dtlschr = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1354	'dtlschr', 'Large-scale condensational heating rate', 'K/s', (/ ('', i=1, 9) /))
1355	TYPE(ctrl_out), SAVE :: o_dqlsc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1356	'dqlsc', 'Condensation dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1357	TYPE(ctrl_out), SAVE :: o_dqlsc2d = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1358	'dqlsc2d', 'Condensation dQ', '(kg/m2)/s', (/ ('', i=1, 9) /))
1359	TYPE(ctrl_out), SAVE :: o_beta_prec = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1360	'beta_prec', 'LS Conversion rate to prec', '(kg/kg)/s', (/ ('', i=1, 9) /))
1361	TYPE(ctrl_out), SAVE :: o_dtvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1362	'dtvdf', 'Boundary-layer dT', 'K/s', (/ ('', i=1, 9) /))
1363	TYPE(ctrl_out), SAVE :: o_dtdis = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1364	'dtdis', 'TKE dissipation dT', 'K/s', (/ ('', i=1, 9) /))
1365	TYPE(ctrl_out), SAVE :: o_dqvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1366	'dqvdf', 'Boundary-layer dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1367	TYPE(ctrl_out), SAVE :: o_dqvdf2d = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1368	'dqvdf2d', 'Boundary-layer dQ', '(kg/m2)/s', (/ ('', i=1, 9) /))
1369	TYPE(ctrl_out), SAVE :: o_dteva = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1370	'dteva', 'Reevaporation dT', 'K/s', (/ ('', i=1, 9) /))
1371	TYPE(ctrl_out), SAVE :: o_dqeva = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1372	'dqeva', 'Reevaporation dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1373	TYPE(ctrl_out), SAVE :: o_dqeva2d = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1374	'dqeva2d', 'Reevaporation dQ', '(kg/m2)/s', (/ ('', i=1, 9) /))
1375
1376	!!!!!!!!!!!!!!!! Specifique thermiques
1377	TYPE(ctrl_out), SAVE :: o_dqlscth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1378	'dqlscth', 'dQ therm.', '(kg/kg)/s', (/ ('', i=1, 9) /))
1379	TYPE(ctrl_out), SAVE :: o_dqlscth2d = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1380	'dqlscth2d', 'dQ therm.', '(kg/m2)/s', (/ ('', i=1, 9) /))
1381	TYPE(ctrl_out), SAVE :: o_dqlscst = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1382	'dqlscst', 'dQ strat.', '(kg/kg)/s', (/ ('', i=1, 9) /))
1383	TYPE(ctrl_out), SAVE :: o_dqlscst2d = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1384	'dqlscst2d', 'dQ strat.', '(kg/m2)/s', (/ ('', i=1, 9) /))
1385	TYPE(ctrl_out), SAVE :: o_dtlscth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1386	'dtlscth', 'dQ therm.', 'K/s', (/ ('', i=1, 9) /))
1387	TYPE(ctrl_out), SAVE :: o_dtlscst = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1388	'dtlscst', 'dQ strat.', 'K/s', (/ ('', i=1, 9) /))
1389	TYPE(ctrl_out), SAVE :: o_plulth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1390	'plulth', 'Rainfall therm.', 'K/s', (/ ('', i=1, 9) /))
1391	TYPE(ctrl_out), SAVE :: o_plulst = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1392	'plulst', 'Rainfall strat.', 'K/s', (/ ('', i=1, 9) /))
1393	TYPE(ctrl_out), SAVE :: o_lmaxth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1394	'lmaxth', "Upper level thermals", "", (/ ('', i=1, 9) /))
1395	TYPE(ctrl_out), SAVE :: o_ptconvth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1396	'ptconvth', 'POINTS CONVECTIFS therm.', ' ', (/ ('', i=1, 9) /))
1397	!!!!!!!!!!!!!!!!!!!!!!!!
1398	TYPE(ctrl_out), SAVE :: o_ptconv = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1399	'ptconv', 'POINTS CONVECTIFS', ' ', (/ ('', i=1, 9) /))
1400	TYPE(ctrl_out), SAVE :: o_ratqs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1401	'ratqs', 'RATQS', ' ', (/ ('', i=1, 9) /))
1402	TYPE(ctrl_out), SAVE :: o_dtthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1403	'dtthe', 'Thermal dT', 'K/s', (/ ('', i=1, 9) /))
1404	TYPE(ctrl_out), SAVE :: o_duthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1405	'duthe', 'Thermal du', 'm/s2', (/ ('', i=1, 9) /))
1406	TYPE(ctrl_out), SAVE :: o_dvthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1407	'dvthe', 'Thermal dv', 'm/s2', (/ ('', i=1, 9) /))
1408	TYPE(ctrl_out), SAVE :: o_f_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1409	'f_th', 'Thermal plume mass flux', 'kg/(m2*s)', (/ ('', i=1, 9) /))
1410	TYPE(ctrl_out), SAVE :: o_e_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1411	'e_th', 'Thermal plume entrainment', 'K/s', (/ ('', i=1, 9) /))
1412	TYPE(ctrl_out), SAVE :: o_w_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1413	'w_th', 'Thermal plume vertical velocity', 'm/s', (/ ('', i=1, 9) /))
1414	TYPE(ctrl_out), SAVE :: o_lambda_th = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1415	'lambda_th', 'Thermal plume vertical velocity', 'm/s', (/ ('', i=1, 9) /))
1416	TYPE(ctrl_out), SAVE :: o_ftime_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1417	'ftime_th', 'Fraction of time Shallow convection occurs', ' ', (/ ('', i=1, 9) /))
1418	TYPE(ctrl_out), SAVE :: o_q_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1419	'q_th', 'Thermal plume total humidity', 'kg/kg', (/ ('', i=1, 9) /))
1420	TYPE(ctrl_out), SAVE :: o_a_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1421	'a_th', "Thermal plume fraction", "", (/ ('', i=1, 9) /))
1422	TYPE(ctrl_out), SAVE :: o_d_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1423	'd_th', 'Thermal plume detrainment', 'K/s', (/ ('', i=1, 9) /))
1424	TYPE(ctrl_out), SAVE :: o_f0_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1425	'f0_th', 'Thermal closure mass flux', 'K/s', (/ ('', i=1, 9) /))
1426	TYPE(ctrl_out), SAVE :: o_zmax_th = ctrl_out((/ 4, 4, 4, 5, 10, 10, 11, 11, 11 /), &
1427	'zmax_th', 'Thermal plume height', 'K/s', (/ ('', i=1, 9) /))
1428	TYPE(ctrl_out), SAVE :: o_dqthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1429	'dqthe', 'Thermal dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1430	TYPE(ctrl_out), SAVE :: o_dqthe2d = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1431	'dqthe2d', 'Thermal dQ', '(kg/m2)/s', (/ ('', i=1, 9) /))
1432	TYPE(ctrl_out), SAVE :: o_dtajs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1433	'dtajs', 'Dry adjust. dT', 'K/s', (/ ('', i=1, 9) /))
1434	TYPE(ctrl_out), SAVE :: o_dqajs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1435	'dqajs', 'Dry adjust. dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1436	TYPE(ctrl_out), SAVE :: o_dqajs2d = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1437	'dqajs2d', 'Dry adjust. dQ', '(kg/m2)/s', (/ ('', i=1, 9) /))
1438	TYPE(ctrl_out), SAVE :: o_dtswr = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1439	'dtswr', 'SW radiation dT', 'K/s', (/ ('', i=1, 9) /))
1440	TYPE(ctrl_out), SAVE :: o_dtsw0 = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1441	'dtsw0', 'CS SW radiation dT', 'K/s', (/ ('', i=1, 9) /))
1442	TYPE(ctrl_out), SAVE :: o_dtlwr = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1443	'dtlwr', 'LW radiation dT', 'K/s', (/ ('', i=1, 9) /))
1444	TYPE(ctrl_out), SAVE :: o_dtlw0 = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1445	'dtlw0', 'CS LW radiation dT', 'K/s', (/ ('', i=1, 9) /))
1446	TYPE(ctrl_out), SAVE :: o_dtec = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1447	'dtec', 'Cinetic dissip dT', 'K/s', (/ ('', i=1, 9) /))
1448	TYPE(ctrl_out), SAVE :: o_duvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1449	'duvdf', 'Boundary-layer dU', 'm/s2', (/ ('', i=1, 9) /))
1450	TYPE(ctrl_out), SAVE :: o_dvvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1451	'dvvdf', 'Boundary-layer dV', 'm/s2', (/ ('', i=1, 9) /))
1452	TYPE(ctrl_out), SAVE :: o_duoro = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1453	'duoro', 'Orography dU', 'm/s2', (/ ('', i=1, 9) /))
1454	TYPE(ctrl_out), SAVE :: o_dvoro = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1455	'dvoro', 'Orography dV', 'm/s2', (/ ('', i=1, 9) /))
1456	TYPE(ctrl_out), SAVE :: o_dulif = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1457	'dulif', 'Orography dU', 'm/s2', (/ ('', i=1, 9) /))
1458	TYPE(ctrl_out), SAVE :: o_dvlif = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1459	'dvlif', 'Orography dV', 'm/s2', (/ ('', i=1, 9) /))
1460	TYPE(ctrl_out), SAVE :: o_du_gwd_hines = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1461	'du_gwd_hines', 'Hines GWD dU', 'm/s2', (/ ('', i=1, 9) /))
1462	TYPE(ctrl_out), SAVE :: o_dv_gwd_hines = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1463	'dv_gwd_hines', 'Hines GWD dV', 'm/s2', (/ ('', i=1, 9) /))
1464	TYPE(ctrl_out), SAVE :: o_du_gwd_front = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1465	'du_gwd_front', 'Fronts GWD dU', 'm/s2', (/ ('', i=1, 9) /))
1466	TYPE(ctrl_out), SAVE :: o_dv_gwd_front = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1467	'dv_gwd_front', 'Fronts GWD dV', 'm/s2', (/ ('', i=1, 9) /))
1468	TYPE(ctrl_out), SAVE :: o_east_gwstress = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1469	'east_gwstress', 'Eastward GW Stress', 'Pa', (/ ('', i=1, 9) /))
1470	TYPE(ctrl_out), SAVE :: o_west_gwstress = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1471	'west_gwstress', 'Westward GW Stress', 'Pa', (/ ('', i=1, 9) /))
1472	TYPE(ctrl_out), SAVE :: o_dtoro = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1473	'dtoro', 'Orography dT', 'K/s', (/ ('', i=1, 9) /))
1474	TYPE(ctrl_out), SAVE :: o_dtlif = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1475	'dtlif', 'Orography dT', 'K/s', (/ ('', i=1, 9) /))
1476	TYPE(ctrl_out), SAVE :: o_dthin = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1477	'dthin', 'Hines GWD dT', 'K/s', (/ ('', i=1, 9) /))
1478	TYPE(ctrl_out), SAVE :: o_dqch4 = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1479	'dqch4', 'H2O due to CH4 oxidation & photolysis', '(kg/kg)/s', (/ ('', i=1, 9) /))
1480
1481	type(ctrl_out), save:: o_du_gwd_rando &
1482	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'du_gwd_rando', &
1483	"Random gravity waves dU/dt", "m/s2", (/ ('', i=1, 9) /))
1484	type(ctrl_out), save:: o_dv_gwd_rando &
1485	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'dv_gwd_rando', &
1486	"Random gravity waves dV/dt", "m/s2", (/ ('', i=1, 9) /))
1487	type(ctrl_out), save:: o_ustr_gwd_hines &
1488	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'ustr_gwd_hines', &
1489	"zonal wind stress Hines gravity waves", "Pa", (/ ('', i=1, 9) /))
1490	type(ctrl_out), save:: o_vstr_gwd_hines &
1491	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'vstr_gwd_hines', &
1492	"meridional wind stress Hines gravity waves", "Pa", (/ ('', i=1, 9) /))
1493	type(ctrl_out), save:: o_ustr_gwd_front &
1494	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'ustr_gwd_front', &
1495	"zonal wind stress fronts gravity waves", "Pa", (/ ('', i=1, 9) /))
1496	type(ctrl_out), save:: o_vstr_gwd_front &
1497	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'vstr_gwd_front', &
1498	"meridional wind stress fronts gravity waves", "Pa", (/ ('', i=1, 9) /))
1499	type(ctrl_out), save:: o_ustr_gwd_rando &
1500	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'ustr_gwd_rando', &
1501	"zonal wind stress random gravity waves", "Pa", (/ ('', i=1, 9) /))
1502	type(ctrl_out), save:: o_vstr_gwd_rando &
1503	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'vstr_gwd_rando', &
1504	"meridional wind stress random gravity waves", "Pa", (/ ('', i=1, 9) /))
1505
1506	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac(:)
1507	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac_cum(:)
1508	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_vdf(:)
1509	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_the(:)
1510	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_con(:)
1511	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_lessi_impa(:)
1512	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_lessi_nucl(:)
1513	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_insc(:)
1514	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_bcscav(:)
1515	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_evapls(:)
1516	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_ls(:)
1517	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_trsp(:)
1518	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_sscav(:)
1519	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_sat(:)
1520	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_uscav(:)
1521	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_dry(:)
1522
1523	TYPE(ctrl_out), SAVE :: o_rsu = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1524	'rsu', 'SW upward radiation', 'W m-2', (/ ('', i=1, 9) /))
1525	TYPE(ctrl_out), SAVE :: o_rsd = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1526	'rsd', 'SW downward radiation', 'W m-2', (/ ('', i=1, 9) /))
1527	TYPE(ctrl_out), SAVE :: o_rlu = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1528	'rlu', 'LW upward radiation', 'W m-2', (/ ('', i=1, 9) /))
1529	TYPE(ctrl_out), SAVE :: o_rld = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1530	'rld', 'LW downward radiation', 'W m-2', (/ ('', i=1, 9) /))
1531	TYPE(ctrl_out), SAVE :: o_rsucs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1532	'rsucs', 'SW CS upward radiation', 'W m-2', (/ ('', i=1, 9) /))
1533	TYPE(ctrl_out), SAVE :: o_rsdcs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1534	'rsdcs', 'SW CS downward radiation', 'W m-2', (/ ('', i=1, 9) /))
1535	TYPE(ctrl_out), SAVE :: o_rlucs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1536	'rlucs', 'LW CS upward radiation', 'W m-2', (/ ('', i=1, 9) /))
1537	TYPE(ctrl_out), SAVE :: o_rldcs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1538	'rldcs', 'LW CS downward radiation', 'W m-2', (/ ('', i=1, 9) /))
1539	TYPE(ctrl_out), SAVE :: o_tnt = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1540	'tnt', 'Tendency of air temperature', 'K s-1', (/ ('', i=1, 9) /))
1541	TYPE(ctrl_out), SAVE :: o_tntc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1542	'tntc', 'Tendency of air temperature due to Moist Convection', 'K s-1', (/ ('', i=1, 9) /))
1543	TYPE(ctrl_out), SAVE :: o_tntr = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1544	'tntr', 'Air temperature tendency due to Radiative heating', 'K s-1', (/ ('', i=1, 9) /))
1545	TYPE(ctrl_out), SAVE :: o_tntscpbl = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1546	'tntscpbl', 'Air temperature tendency due to St cloud and precipitation and BL mixing', &
1547	'K s-1', (/ ('', i=1, 9) /))
1548	TYPE(ctrl_out), SAVE :: o_tnhus = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1549	'tnhus', 'Tendency of specific humidity', 's-1', (/ ('', i=1, 9) /))
1550	TYPE(ctrl_out), SAVE :: o_tnhusc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1551	'tnhusc', 'Tendency of specific humidity due to convection', 's-1', (/ ('', i=1, 9) /))
1552	TYPE(ctrl_out), SAVE :: o_tnhusscpbl = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1553	'tnhusscpbl', 'Tendency of Specific humidity due to ST cl, precip and BL mixing', 's-1', (/ ('', i=1, 9) /))
1554	TYPE(ctrl_out), SAVE :: o_evu = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1555	'evu', 'Eddy viscosity coefficient for Momentum Variables', 'm2 s-1', (/ ('', i=1, 9) /))
1556	TYPE(ctrl_out), SAVE :: o_h2o = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1557	'h2o', 'Mass Fraction of Water', '1', (/ ('', i=1, 9) /))
1558	TYPE(ctrl_out), SAVE :: o_mcd = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1559	'mcd', 'Downdraft COnvective Mass Flux', 'kg/(m2*s)', (/ ('', i=1, 9) /))
1560	TYPE(ctrl_out), SAVE :: o_dmc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1561	'dmc', 'Deep COnvective Mass Flux', 'kg/(m2*s)', (/ ('', i=1, 9) /))
1562	TYPE(ctrl_out), SAVE :: o_ref_liq = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1563	'ref_liq', 'Effective radius of convective cloud liquid water particle', 'm', (/ ('', i=1, 9) /))
1564	TYPE(ctrl_out), SAVE :: o_ref_ice = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1565	'ref_ice', 'Effective radius of startiform cloud ice particle', 'm', (/ ('', i=1, 9) /))
1566	TYPE(ctrl_out), SAVE :: o_rsut4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1567	'rsut4co2', 'TOA Out SW in 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1568	TYPE(ctrl_out), SAVE :: o_rlut4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1569	'rlut4co2', 'TOA Out LW in 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1570	TYPE(ctrl_out), SAVE :: o_rsutcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1571	'rsutcs4co2', 'TOA Out CS SW in 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1572	TYPE(ctrl_out), SAVE :: o_rlutcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1573	'rlutcs4co2', 'TOA Out CS LW in 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1574	TYPE(ctrl_out), SAVE :: o_rsu4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1575	'rsu4co2', 'Upwelling SW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1576	TYPE(ctrl_out), SAVE :: o_rlu4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1577	'rlu4co2', 'Upwelling LW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1578	TYPE(ctrl_out), SAVE :: o_rsucs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1579	'rsucs4co2', 'Upwelling CS SW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1580	TYPE(ctrl_out), SAVE :: o_rlucs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1581	'rlucs4co2', 'Upwelling CS LW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1582	TYPE(ctrl_out), SAVE :: o_rsd4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1583	'rsd4co2', 'Downwelling SW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1584	TYPE(ctrl_out), SAVE :: o_rld4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1585	'rld4co2', 'Downwelling LW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1586	TYPE(ctrl_out), SAVE :: o_rsdcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1587	'rsdcs4co2', 'Downwelling CS SW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1588	TYPE(ctrl_out), SAVE :: o_rldcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1589	'rldcs4co2', 'Downwelling CS LW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1590	TYPE(ctrl_out), SAVE :: o_snowsrf = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1591	'snowsrf', 'Snow mass at surface', 'kg/m2', (/ ('', i=1, 9) /))
1592	TYPE(ctrl_out), SAVE :: o_qsnow = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1593	'qsnow', 'Water contained in snow', 'kg/m2', (/ ('', i=1, 9) /))
1594	TYPE(ctrl_out), SAVE :: o_snowhgt = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1595	'snowhgt', 'Snow height at surface', 'm', (/ ('', i=1, 9) /))
1596	TYPE(ctrl_out), SAVE :: o_toice = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1597	'to_ice', 'Snow passed to ice model', 'kg/m2', (/ ('', i=1, 9) /))
1598	TYPE(ctrl_out), SAVE :: o_sissnow = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1599	'sissnow', 'Snow in snow model', 'kg/m2', (/ ('', i=1, 9) /))
1600	TYPE(ctrl_out), SAVE :: o_runoff = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1601	'runoff', 'Run-off rate land ice', 'kg/m2/s', (/ ('', i=1, 9) /))
1602	TYPE(ctrl_out), SAVE :: o_albslw3 = ctrl_out((/ 1, 1, 1, 1, 10, 10, 11, 11, 11 /), &
1603	'albslw3', 'Surface albedo LW3', '-', (/ ('', i=1, 9) /))
1604
1605	!!!!!!!!!!!!! Sorties niveaux standards de pression NMC
1606	TYPE(ctrl_out), SAVE :: o_tnondef = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1607	'tnondef', 'Undefined value of T', 'K', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1608	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1609	TYPE(ctrl_out), SAVE :: o_ta = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1610	'ta', 'Air temperature', 'K', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1611	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1612	TYPE(ctrl_out), SAVE :: o_zg = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1613	'zg', 'Geopotential height', 'm', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1614	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1615	TYPE(ctrl_out), SAVE :: o_hus = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1616	'hus', 'Specific humidity', '1', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1617	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1618	TYPE(ctrl_out), SAVE :: o_hur = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1619	'hur', 'Relative humidity', '%', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1620	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1621	TYPE(ctrl_out), SAVE :: o_ua = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1622	'ua', 'Eastward wind', 'm s-1', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1623	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1624	TYPE(ctrl_out), SAVE :: o_va = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1625	'va', 'Northward wind', 'm s-1', (/ ('', i=1, 9)/))
1626	TYPE(ctrl_out), SAVE :: o_wap = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1627	'wap', 'Lagrangian tendency of air pressure', 'Pa s-1', (/ "inst(X)", "inst(X)", "inst(X)", &
1628	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1629	TYPE(ctrl_out), SAVE :: o_psbg = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1630	'psbg', 'Pressure sfce below ground', '%', (/ "inst(X)", "inst(X)", "inst(X)", &
1631	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1632	TYPE(ctrl_out), SAVE :: o_tro3 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1633	'tro3', 'Ozone mole fraction', '1e-9', (/ "inst(X)", "inst(X)", "inst(X)", &
1634	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1635	TYPE(ctrl_out), SAVE :: o_tro3_daylight = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1636	'tro3_daylight', 'Daylight ozone mole fraction', '1e-9', (/ "inst(X)", "inst(X)", "inst(X)", &
1637	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1638	TYPE(ctrl_out), SAVE :: o_uxv = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1639	'uv', 'uv', 'm2/s2', (/ "inst(X)", "inst(X)", "inst(X)", &
1640	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1641	TYPE(ctrl_out), SAVE :: o_vxq = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1642	'vxq', 'vxq', 'm/s * (kg/kg)', (/ "inst(X)", "inst(X)", "inst(X)", &
1643	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1644	TYPE(ctrl_out), SAVE :: o_vxT = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1645	'vT', 'vT', 'mK/s', (/ "inst(X)", "inst(X)", "inst(X)", &
1646	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1647	TYPE(ctrl_out), SAVE :: o_wxq = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1648	'wq', 'wq', '(Pa/s)*(kg/kg)', (/ "inst(X)", "inst(X)", "inst(X)", &
1649	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1650	TYPE(ctrl_out), SAVE :: o_vxphi = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1651	'vphi', 'vphi', 'm2/s', (/ "inst(X)", "inst(X)", "inst(X)", &
1652	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1653	TYPE(ctrl_out), SAVE :: o_wxT = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1654	'wT', 'wT', '"K*Pa/s', (/ "inst(X)", "inst(X)", "inst(X)", &
1655	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1656	TYPE(ctrl_out), SAVE :: o_uxu = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1657	'u2', 'u2', 'm2/s2', (/ "inst(X)", "inst(X)", "inst(X)", &
1658	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1659	TYPE(ctrl_out), SAVE :: o_vxv = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1660	'v2', 'v2', 'm2/s2', (/ "inst(X)", "inst(X)", "inst(X)", &
1661	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1662	TYPE(ctrl_out), SAVE :: o_TxT = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1663	'T2', 'T2', 'K2', (/ "inst(X)", "inst(X)", "inst(X)", &
1664	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1665
1666	#ifdef CPP_Dust
1667	#include "Dust/spla_output_dat.h"
1668	#endif
1669
1670	END MODULE phys_output_ctrlout_mod

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

Download in other formats:

Original Format