Context Navigation

phys_output_ctrlout_mod.F90 @ 2003

Last change on this file since 2003 was 2003, checked in by Laurent Fairhead, 10 years ago

Nouvelle version qui inclut les effets des aérosols et propose les mêmes diagnostics des effets
directs et indirects que l'ancienne version du rayonnement.
OB

New RRTM version that includes the effects of aerosols and outputs the same direct and indirect effects
diagnostics as the old version
OB

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: 92.5 KB

Line
1	MODULE phys_output_ctrlout_mod
2
3	USE phys_output_var_mod
4	USE indice_sol_mod
5	USE aero_mod, only : naero_tot,name_aero_tau
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	TYPE(ctrl_out), SAVE :: o_flat = ctrl_out((/ 5, 1, 10, 10, 5, 10, 11, 11, 11 /), &
44	'flat', 'Latent heat flux', 'W/m2', (/ ('', i=1, 9) /))
45	TYPE(ctrl_out), SAVE :: o_slp = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
46	'slp', 'Sea Level Pressure', 'Pa', (/ ('', i=1, 9) /))
47	TYPE(ctrl_out), SAVE :: o_tsol = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
48	'tsol', 'Surface Temperature', 'K', (/ ('', i=1, 9) /))
49	TYPE(ctrl_out), SAVE :: o_t2m = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
50	't2m', 'Temperature 2m', 'K', (/ ('', i=1, 9) /))
51	TYPE(ctrl_out), SAVE :: o_t2m_min = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
52	't2m_min', 'Temp 2m min', 'K', &
53	(/ "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)" /))
54	TYPE(ctrl_out), SAVE :: o_t2m_max = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
55	't2m_max', 'Temp 2m max', 'K', &
56	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
57	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /))
58
59	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_t2m_srf = (/ &
60	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /), &
61	't2m_ter', "Temp 2m "//clnsurf(1), "K", (/ ('', i=1, 9) /)), &
62	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /), &
63	't2m_lic', "Temp 2m "//clnsurf(2), "K", (/ ('', i=1, 9) /)), &
64	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /), &
65	't2m_oce', "Temp 2m "//clnsurf(3), "K", (/ ('', i=1, 9) /)), &
66	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /), &
67	't2m_sic', "Temp 2m "//clnsurf(4), "K", (/ ('', i=1, 9) /)) /)
68
69	TYPE(ctrl_out), SAVE :: o_wind10m = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
70	'wind10m', '10-m wind speed', 'm/s', (/ ('', i=1, 9) /))
71	TYPE(ctrl_out), SAVE :: o_wind10max = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
72	'wind10max', '10m wind speed max', 'm/s', &
73	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
74	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /))
75	TYPE(ctrl_out), SAVE :: o_sicf = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
76	'sicf', 'Sea-ice fraction', '-', (/ ('', i=1, 9) /))
77	TYPE(ctrl_out), SAVE :: o_q2m = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
78	'q2m', 'Specific humidity 2m', 'kg/kg', (/ ('', i=1, 9) /))
79	TYPE(ctrl_out), SAVE :: o_ustar = ctrl_out((/ 1, 1, 10, 5, 10, 10, 11, 11, 11 /), &
80	'ustar', 'Friction velocity', 'm/s', (/ ('', i=1, 9) /))
81	TYPE(ctrl_out), SAVE :: o_u10m = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
82	'u10m', 'Vent zonal 10m', 'm/s', (/ ('', i=1, 9) /))
83	TYPE(ctrl_out), SAVE :: o_v10m = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
84	'v10m', 'Vent meridien 10m', 'm/s', (/ ('', i=1, 9) /))
85	TYPE(ctrl_out), SAVE :: o_psol = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
86	'psol', 'Surface Pressure', 'Pa', (/ ('', i=1, 9) /))
87	TYPE(ctrl_out), SAVE :: o_qsurf = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
88	'qsurf', 'Surface Air humidity', 'kg/kg', (/ ('', i=1, 9) /))
89
90	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_ustar_srf = (/ &
91	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'ustar_ter', &
92	"Friction velocity "//clnsurf(1),"m/s", (/ ('', i=1, 9) /)), &
93	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'ustar_lic', &
94	"Friction velocity "//clnsurf(2),"m/s", (/ ('', i=1, 9) /)), &
95	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'ustar_oce', &
96	"Friction velocity "//clnsurf(3),"m/s", (/ ('', i=1, 9) /)), &
97	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'ustar_sic', &
98	"Friction velocity "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)) /)
99
100	TYPE(ctrl_out), SAVE, DIMENSION(5) :: o_wstar = (/ &
101	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'wstar_ter', &
102	"Friction velocity "//clnsurf(1),"m/s", (/ ('', i=1, 9) /)), &
103	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'wstar_lic', &
104	"Friction velocity "//clnsurf(2),"m/s", (/ ('', i=1, 9) /)), &
105	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'wstar_oce', &
106	"Friction velocity "//clnsurf(3),"m/s", (/ ('', i=1, 9) /)), &
107	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'wstar_sic', &
108	"Friction velocity "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)), &
109	ctrl_out((/ 5, 5, 10, 10, 10, 10, 11, 11, 11 /),'wstar', &
110	"w* convective velocity "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)) /)
111
112	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_u10m_srf = (/ &
113	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'u10m_ter', &
114	"Vent Zonal 10m "//clnsurf(1),"m/s", (/ ('', i=1, 9) /)), &
115	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'u10m_lic', &
116	"Vent Zonal 10m "//clnsurf(2),"m/s", (/ ('', i=1, 9) /)), &
117	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'u10m_oce', &
118	"Vent Zonal 10m "//clnsurf(3),"m/s", (/ ('', i=1, 9) /)), &
119	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'u10m_sic', &
120	"Vent Zonal 10m "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)) /)
121
122	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_v10m_srf = (/ &
123	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'v10m_ter', &
124	"Vent meredien 10m "//clnsurf(1),"m/s", (/ ('', i=1, 9) /)), &
125	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'v10m_lic', &
126	"Vent meredien 10m "//clnsurf(2),"m/s", (/ ('', i=1, 9) /)), &
127	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'v10m_oce', &
128	"Vent meredien 10m "//clnsurf(3),"m/s", (/ ('', i=1, 9) /)), &
129	ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11 /),'v10m_sic', &
130	"Vent meredien 10m "//clnsurf(4),"m/s", (/ ('', i=1, 9) /)) /)
131
132	TYPE(ctrl_out), SAVE :: o_qsol = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
133	'qsol', 'Soil watter content', 'mm', (/ ('', i=1, 9) /))
134	TYPE(ctrl_out), SAVE :: o_ndayrain = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
135	'ndayrain', 'Number of dayrain(liq+sol)', '-', &
136	(/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
137	TYPE(ctrl_out), SAVE :: o_precip = ctrl_out((/ 1, 1, 1, 10, 5, 10, 11, 11, 11 /), &
138	'precip', 'Precip Totale liq+sol', 'kg/(s*m2)', (/ ('', i=1, 9) /))
139	TYPE(ctrl_out), SAVE :: o_plul = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
140	'plul', 'Large-scale Precip.', 'kg/(s*m2)', (/ ('', i=1, 9) /))
141	TYPE(ctrl_out), SAVE :: o_pluc = ctrl_out((/ 1, 1, 1, 10, 5, 10, 11, 11, 11 /), &
142	'pluc', 'Convective Precip.', 'kg/(s*m2)', (/ ('', i=1, 9) /))
143	TYPE(ctrl_out), SAVE :: o_snow = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11 /), &
144	'snow', 'Snow fall', 'kg/(s*m2)', (/ ('', i=1, 9) /))
145	TYPE(ctrl_out), SAVE :: o_evap = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
146	'evap', 'Evaporat', 'kg/(s*m2)', (/ ('', i=1, 9) /))
147
148	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_evap_srf = (/ &
149	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evap_ter', &
150	"evaporation at surface "//clnsurf(1),"kg/(s*m2)", (/ ('', i=1, 9) /)), &
151	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evap_lic', &
152	"evaporation at surface "//clnsurf(2),"kg/(s*m2)", (/ ('', i=1, 9) /)), &
153	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evap_oce', &
154	"evaporation at surface "//clnsurf(3),"kg/(s*m2)", (/ ('', i=1, 9) /)), &
155	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evap_sic', &
156	"evaporation at surface "//clnsurf(4),"kg/(s*m2)", (/ ('', i=1, 9) /)) /)
157
158	TYPE(ctrl_out), SAVE :: o_msnow = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
159	'msnow', 'Surface snow amount', 'kg/m2', (/ ('', i=1, 9) /))
160	TYPE(ctrl_out), SAVE :: o_fsnow = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
161	'fsnow', 'Surface snow area fraction', '-', (/ ('', i=1, 9) /))
162	TYPE(ctrl_out), SAVE :: o_tops = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
163	'tops', 'Solar rad. at TOA', 'W/m2', (/ ('', i=1, 9) /))
164	TYPE(ctrl_out), SAVE :: o_tops0 = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
165	'tops0', 'CS Solar rad. at TOA', 'W/m2', (/ ('', i=1, 9) /))
166	TYPE(ctrl_out), SAVE :: o_topl = ctrl_out((/ 1, 1, 10, 5, 10, 10, 11, 11, 11 /), &
167	'topl', 'IR rad. at TOA', 'W/m2', (/ ('', i=1, 9) /))
168	TYPE(ctrl_out), SAVE :: o_topl0 = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
169	'topl0', 'IR rad. at TOA', 'W/m2', (/ ('', i=1, 9) /))
170	TYPE(ctrl_out), SAVE :: o_SWupTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
171	'SWupTOA', 'SWup at TOA', 'W/m2', (/ ('', i=1, 9) /))
172	TYPE(ctrl_out), SAVE :: o_SWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
173	'SWupTOAclr', 'SWup clear sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
174	TYPE(ctrl_out), SAVE :: o_SWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
175	'SWdnTOA', 'SWdn at TOA', 'W/m2', (/ ('', i=1, 9) /))
176	TYPE(ctrl_out), SAVE :: o_SWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
177	'SWdnTOAclr', 'SWdn clear sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
178	TYPE(ctrl_out), SAVE :: o_nettop = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /), &
179	'nettop', 'Net dn radiatif flux at TOA', 'W/m2', (/ ('', i=1, 9) /))
180	TYPE(ctrl_out), SAVE :: o_SWup200 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
181	'SWup200', 'SWup at 200mb', 'W/m2', (/ ('', i=1, 9) /))
182	TYPE(ctrl_out), SAVE :: o_SWup200clr = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
183	'SWup200clr', 'SWup clear sky at 200mb', 'W/m2', (/ ('', i=1, 9) /))
184	TYPE(ctrl_out), SAVE :: o_SWdn200 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
185	'SWdn200', 'SWdn at 200mb', 'W/m2', (/ ('', i=1, 9) /))
186	TYPE(ctrl_out), SAVE :: o_SWdn200clr = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
187	'SWdn200clr', 'SWdn clear sky at 200mb', 'W/m2', (/ ('', i=1, 9) /))
188
189	! arajouter
190	! type(ctrl_out),save :: o_LWupTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /),'LWupTOA', &
191	! (/ ('', i=1, 9) /))
192	! type(ctrl_out),save :: o_LWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /),'LWupTOAclr', &
193	! (/ ('', i=1, 9) /))
194	! type(ctrl_out),save :: o_LWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /),'LWdnTOA', &
195	! (/ ('', i=1, 9) /))
196	! type(ctrl_out),save :: o_LWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11 /),'LWdnTOAclr', &
197	! (/ ('', i=1, 9) /))
198	TYPE(ctrl_out), SAVE :: o_LWup200 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
199	'LWup200', 'LWup at 200mb', 'W/m2', (/ ('', i=1, 9) /))
200	TYPE(ctrl_out), SAVE :: o_LWup200clr = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
201	'LWup200clr', 'LWup clear sky at 200mb', 'W/m2', (/ ('', i=1, 9) /))
202	TYPE(ctrl_out), SAVE :: o_LWdn200 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
203	'LWdn200', 'LWdn at 200mb', 'W/m2', (/ ('', i=1, 9) /))
204	TYPE(ctrl_out), SAVE :: o_LWdn200clr = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
205	'LWdn200clr', 'LWdn clear sky at 200mb', 'W/m2', (/ ('', i=1, 9) /))
206	TYPE(ctrl_out), SAVE :: o_sols = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
207	'sols', 'Solar rad. at surf.', 'W/m2', (/ ('', i=1, 9) /))
208	TYPE(ctrl_out), SAVE :: o_sols0 = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
209	'sols0', 'Solar rad. at surf.', 'W/m2', (/ ('', i=1, 9) /))
210	TYPE(ctrl_out), SAVE :: o_soll = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
211	'soll', 'IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
212	TYPE(ctrl_out), SAVE :: o_soll0 = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
213	'soll0', 'IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
214	TYPE(ctrl_out), SAVE :: o_radsol = ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /), &
215	'radsol', 'Rayonnement au sol', 'W/m2', (/ ('', i=1, 9) /))
216	TYPE(ctrl_out), SAVE :: o_SWupSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
217	'SWupSFC', 'SWup at surface', 'W/m2', (/ ('', i=1, 9) /))
218	TYPE(ctrl_out), SAVE :: o_SWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
219	'SWupSFCclr', 'SWup clear sky at surface', 'W/m2', (/ ('', i=1, 9) /))
220	TYPE(ctrl_out), SAVE :: o_SWdnSFC = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11 /), &
221	'SWdnSFC', 'SWdn at surface', 'W/m2', (/ ('', i=1, 9) /))
222	TYPE(ctrl_out), SAVE :: o_SWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
223	'SWdnSFCclr', 'SWdn clear sky at surface', 'W/m2', (/ ('', i=1, 9) /))
224	TYPE(ctrl_out), SAVE :: o_LWupSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
225	'LWupSFC', 'Upwd. IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
226	TYPE(ctrl_out), SAVE :: o_LWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
227	'LWupSFCclr', 'CS Upwd. IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
228	TYPE(ctrl_out), SAVE :: o_LWdnSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
229	'LWdnSFC', 'Down. IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
230	TYPE(ctrl_out), SAVE :: o_LWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11 /), &
231	'LWdnSFCclr', 'Down. CS IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
232	TYPE(ctrl_out), SAVE :: o_bils = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
233	'bils', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
234	TYPE(ctrl_out), SAVE :: o_bils_tke = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
235	'bils_tke', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
236	TYPE(ctrl_out), SAVE :: o_bils_diss = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
237	'bils_diss', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
238	TYPE(ctrl_out), SAVE :: o_bils_ec = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
239	'bils_ec', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
240	TYPE(ctrl_out), SAVE :: o_bils_kinetic = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
241	'bils_kinetic', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
242	TYPE(ctrl_out), SAVE :: o_bils_enthalp = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
243	'bils_enthalp', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
244	TYPE(ctrl_out), SAVE :: o_bils_latent = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
245	'bils_latent', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
246	TYPE(ctrl_out), SAVE :: o_sens = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11 /), &
247	'sens', 'Sensible heat flux', 'W/m2', (/ ('', i=1, 9) /))
248	TYPE(ctrl_out), SAVE :: o_fder = ctrl_out((/ 1, 2, 10, 10, 10, 10, 11, 11, 11 /), &
249	'fder', 'Heat flux derivation', 'W/m2', (/ ('', i=1, 9) /))
250	TYPE(ctrl_out), SAVE :: o_ffonte = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
251	'ffonte', 'Thermal flux for snow melting', 'W/m2', (/ ('', i=1, 9) /))
252	TYPE(ctrl_out), SAVE :: o_fqcalving = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
253	'fqcalving', 'Ice Calving', 'kg/m2/s', (/ ('', i=1, 9) /))
254	TYPE(ctrl_out), SAVE :: o_fqfonte = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
255	'fqfonte', 'Land ice melt', 'kg/m2/s', (/ ('', i=1, 9) /))
256	TYPE(ctrl_out), SAVE :: o_taux = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
257	'taux', 'Zonal wind stress', 'Pa', (/ ('', i=1, 9) /))
258	TYPE(ctrl_out), SAVE :: o_tauy = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
259	'tauy', 'Meridional wind stress', 'Pa', (/ ('', i=1, 9) /))
260
261	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_taux_srf = (/ &
262	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'taux_ter', &
263	"Zonal wind stress"//clnsurf(1), "Pa", (/ ('', i=1, 9) /)), &
264	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'taux_lic', &
265	"Zonal wind stress"//clnsurf(2), "Pa", (/ ('', i=1, 9) /)), &
266	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'taux_oce', &
267	"Zonal wind stress"//clnsurf(3), "Pa", (/ ('', i=1, 9) /)), &
268	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'taux_sic', &
269	"Zonal wind stress"//clnsurf(4), "Pa", (/ ('', i=1, 9) /)) /)
270
271	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tauy_srf = (/ &
272	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tauy_ter', &
273	"Meridional wind stress "//clnsurf(1),"Pa", (/ ('', i=1, 9) /)), &
274	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tauy_lic', &
275	"Meridional wind stress "//clnsurf(2),"Pa", (/ ('', i=1, 9) /)), &
276	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tauy_oce', &
277	"Meridional wind stress "//clnsurf(3),"Pa", (/ ('', i=1, 9) /)), &
278	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tauy_sic', &
279	"Meridional wind stress "//clnsurf(4),"Pa", (/ ('', i=1, 9) /)) /)
280
281	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_pourc_srf = (/ &
282	ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /),'pourc_ter', &
283	"% "//clnsurf(1),"%", (/ ('', i=1, 9) /)), &
284	ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /),'pourc_lic', &
285	"% "//clnsurf(2),"%", (/ ('', i=1, 9) /)), &
286	ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /),'pourc_oce', &
287	"% "//clnsurf(3),"%", (/ ('', i=1, 9) /)), &
288	ctrl_out((/ 1, 7, 10, 10, 10, 10, 11, 11, 11 /),'pourc_sic', &
289	"% "//clnsurf(4),"%", (/ ('', i=1, 9) /)) /)
290
291	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_fract_srf = (/ &
292	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'fract_ter', &
293	"Fraction "//clnsurf(1),"1", (/ ('', i=1, 9) /)), &
294	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'fract_lic', &
295	"Fraction "//clnsurf(2),"1", (/ ('', i=1, 9) /)), &
296	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'fract_oce', &
297	"Fraction "//clnsurf(3),"1", (/ ('', i=1, 9) /)), &
298	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'fract_sic', &
299	"Fraction "//clnsurf(4),"1", (/ ('', i=1, 9) /)) /)
300
301	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tsol_srf = (/ &
302	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tsol_ter', &
303	"Temperature "//clnsurf(1),"K", (/ ('', i=1, 9) /)), &
304	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tsol_lic', &
305	"Temperature "//clnsurf(2),"K", (/ ('', i=1, 9) /)), &
306	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tsol_oce', &
307	"Temperature "//clnsurf(3),"K", (/ ('', i=1, 9) /)), &
308	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'tsol_sic', &
309	"Temperature "//clnsurf(4),"K", (/ ('', i=1, 9) /)) /)
310
311	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_evappot_srf = (/ &
312	ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11 /),'evappot_ter', &
313	"Temperature"//clnsurf(1),"K", (/ ('', i=1, 9) /)), &
314	ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /),'evappot_lic', &
315	"Temperature"//clnsurf(2),"K", (/ ('', i=1, 9) /)), &
316	ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /),'evappot_oce', &
317	"Temperature"//clnsurf(3),"K", (/ ('', i=1, 9) /)), &
318	ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /),'evappot_sic', &
319	"Temperature"//clnsurf(4),"K", (/ ('', i=1, 9) /)) /)
320
321	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_sens_srf = (/ &
322	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'sens_ter', &
323	"Sensible heat flux "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
324	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'sens_lic', &
325	"Sensible heat flux "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
326	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'sens_oce', &
327	"Sensible heat flux "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
328	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'sens_sic', &
329	"Sensible heat flux "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
330
331	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_lat_srf = (/ &
332	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'lat_ter', &
333	"Latent heat flux "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
334	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'lat_lic', &
335	"Latent heat flux "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
336	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'lat_oce', &
337	"Latent heat flux "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
338	ctrl_out((/ 1, 6, 10, 7, 10, 10, 11, 11, 11 /),'lat_sic', &
339	"Latent heat flux "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
340
341	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_flw_srf = (/ &
342	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'flw_ter', &
343	"LW "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
344	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'flw_lic', &
345	"LW "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
346	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'flw_oce', &
347	"LW "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
348	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'flw_sic', &
349	"LW "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
350
351	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_fsw_srf = (/ &
352	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'fsw_ter', &
353	"SW "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
354	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'fsw_lic', &
355	"SW "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
356	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'fsw_oce', &
357	"SW "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
358	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'fsw_sic', &
359	"SW "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
360
361	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_wbils_srf = (/ &
362	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbils_ter', &
363	"Bilan sol "//clnsurf(1),"W/m2", (/ ('', i=1, 9) /)), &
364	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbils_lic', &
365	"Bilan sol "//clnsurf(2),"W/m2", (/ ('', i=1, 9) /)), &
366	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbils_oce', &
367	"Bilan sol "//clnsurf(3),"W/m2", (/ ('', i=1, 9) /)), &
368	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbils_sic', &
369	"Bilan sol "//clnsurf(4),"W/m2", (/ ('', i=1, 9) /)) /)
370
371	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_wbilo_srf = (/ &
372	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbilo_ter', &
373	"Bilan eau "//clnsurf(1),"kg/(m2*s)", (/ ('', i=1, 9) /)), &
374	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbilo_lic', &
375	"Bilan eau "//clnsurf(2),"kg/(m2*s)", (/ ('', i=1, 9) /)), &
376	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbilo_oce', &
377	"Bilan eau "//clnsurf(3),"kg/(m2*s)", (/ ('', i=1, 9) /)), &
378	ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'wbilo_sic', &
379	"Bilan eau "//clnsurf(4),"kg/(m2*s)", (/ ('', i=1, 9) /)) /)
380
381	TYPE(ctrl_out), SAVE :: o_cdrm = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
382	'cdrm', 'Momentum drag coef.', '-', (/ ('', i=1, 9) /))
383	TYPE(ctrl_out), SAVE :: o_cdrh = ctrl_out((/ 1, 10, 10, 7, 10, 10, 11, 11, 11 /), &
384	'cdrh', 'Heat drag coef.', '-', (/ ('', i=1, 9) /))
385	TYPE(ctrl_out), SAVE :: o_cldl = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
386	'cldl', 'Low-level cloudiness', '-', (/ ('', i=1, 9) /))
387	TYPE(ctrl_out), SAVE :: o_cldm = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
388	'cldm', 'Mid-level cloudiness', '-', (/ ('', i=1, 9) /))
389	TYPE(ctrl_out), SAVE :: o_cldh = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
390	'cldh', 'High-level cloudiness', '-', (/ ('', i=1, 9) /))
391	TYPE(ctrl_out), SAVE :: o_cldt = ctrl_out((/ 1, 1, 2, 10, 5, 10, 11, 11, 11 /), &
392	'cldt', 'Total cloudiness', '-', (/ ('', i=1, 9) /))
393	TYPE(ctrl_out), SAVE :: o_JrNt = ctrl_out((/ 1, 1, 10, 7, 10, 10, 11, 11, 11 /), &
394	'JrNt', '1 if Day 0 if Night', '-', (/ ('', i=1, 9) /))
395	TYPE(ctrl_out), SAVE :: o_cldhjn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
396	'cldhjn', 'High-level cloudiness Day', '-', (/ ('', i=1, 9) /))
397	TYPE(ctrl_out), SAVE :: o_cldmjn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11/), &
398	'cldmjn', 'Mid-level cloudiness day', '-', (/ ('', i=1, 9) /))
399	TYPE(ctrl_out), SAVE :: o_cldljn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11/), &
400	'cldljn', 'Low-level cloudiness day', '-', (/ ('', i=1, 9) /))
401	TYPE(ctrl_out), SAVE :: o_cldtjn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11/), &
402	'cldtjn', 'Total cloudiness day', '-', (/ ('', i=1, 9) /))
403
404	TYPE(ctrl_out), SAVE :: o_cldq = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
405	'cldq', 'Cloud liquid water path', 'kg/m2', (/ ('', i=1, 9) /))
406	TYPE(ctrl_out), SAVE :: o_lwp = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
407	'lwp', 'Cloud water path', 'kg/m2', (/ ('', i=1, 9) /))
408	TYPE(ctrl_out), SAVE :: o_iwp = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
409	'iwp', 'Cloud ice water path', 'kg/m2', (/ ('', i=1, 9) /))
410	TYPE(ctrl_out), SAVE :: o_ue = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
411	'ue', 'Zonal energy transport', '-', (/ ('', i=1, 9) /))
412	TYPE(ctrl_out), SAVE :: o_ve = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
413	've', 'Merid energy transport', '-', (/ ('', i=1, 9) /))
414	TYPE(ctrl_out), SAVE :: o_uq = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
415	'uq', 'Zonal humidity transport', '-', (/ ('', i=1, 9) /))
416	TYPE(ctrl_out), SAVE :: o_vq = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
417	'vq', 'Merid humidity transport', '-', (/ ('', i=1, 9) /))
418	TYPE(ctrl_out), SAVE :: o_cape = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
419	'cape', 'Conv avlbl pot ener', 'J/kg', (/ ('', i=1, 9) /))
420	TYPE(ctrl_out), SAVE :: o_pbase = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
421	'pbase', 'Cld base pressure', 'Pa', (/ ('', i=1, 9) /))
422	TYPE(ctrl_out), SAVE :: o_ptop = ctrl_out((/ 1, 5, 10, 10, 10, 10, 11, 11, 11 /), &
423	'ptop', 'Cld top pressure', 'Pa', (/ ('', i=1, 9) /))
424	TYPE(ctrl_out), SAVE :: o_fbase = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
425	'fbase', 'Cld base mass flux', 'kg/m2/s', (/ ('', i=1, 9) /))
426	TYPE(ctrl_out), SAVE :: o_plcl = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
427	'plcl', 'Lifting Condensation Level', 'hPa', (/ ('', i=1, 9) /))
428	TYPE(ctrl_out), SAVE :: o_plfc = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
429	'plfc', 'Level of Free Convection', 'hPa', (/ ('', i=1, 9) /))
430	TYPE(ctrl_out), SAVE :: o_wbeff = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
431	'wbeff', 'Conv. updraft velocity at LFC (<100)', 'm/s', (/ ('', i=1, 9) /))
432	TYPE(ctrl_out), SAVE :: o_prw = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
433	'prw', 'Precipitable water', 'kg/m2', (/ ('', i=1, 9) /))
434	TYPE(ctrl_out), SAVE :: o_s_pblh = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
435	's_pblh', 'Boundary Layer Height', 'm', (/ ('', i=1, 9) /))
436	TYPE(ctrl_out), SAVE :: o_s_pblt = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
437	's_pblt', 't at Boundary Layer Height', 'K', (/ ('', i=1, 9) /))
438	TYPE(ctrl_out), SAVE :: o_s_lcl = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
439	's_lcl', 'Condensation level', 'm', (/ ('', i=1, 9) /))
440	TYPE(ctrl_out), SAVE :: o_s_therm = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
441	's_therm', 'Exces du thermique', 'K', (/ ('', i=1, 9) /))
442	!IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
443	! type(ctrl_out),save :: o_s_capCL = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_capCL', &
444	! (/ ('', i=1, 9) /))
445	! type(ctrl_out),save :: o_s_oliqCL = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_oliqCL', &
446	! (/ ('', i=1, 9) /))
447	! type(ctrl_out),save :: o_s_cteiCL = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_cteiCL', &
448	! (/ ('', i=1, 9) /))
449	! type(ctrl_out),save :: o_s_trmb1 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_trmb1', &
450	! (/ ('', i=1, 9) /))
451	! type(ctrl_out),save :: o_s_trmb2 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_trmb2', &
452	! (/ ('', i=1, 9) /))
453	! type(ctrl_out),save :: o_s_trmb3 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /),'s_trmb3', &
454	!(/ ('', i=1, 9) /))
455	TYPE(ctrl_out), SAVE :: o_slab_bils = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
456	'slab_bils_oce', 'Bilan au sol sur ocean slab', 'W/m2', (/ ('', i=1, 9) /))
457	TYPE(ctrl_out), SAVE :: o_ale_bl = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
458	'ale_bl', 'ALE BL', 'm2/s2', (/ ('', i=1, 9) /))
459	TYPE(ctrl_out), SAVE :: o_alp_bl = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
460	'alp_bl', 'ALP BL', 'm2/s2', (/ ('', i=1, 9) /))
461	TYPE(ctrl_out), SAVE :: o_ale_wk = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
462	'ale_wk', 'ALE WK', 'm2/s2', (/ ('', i=1, 9) /))
463	TYPE(ctrl_out), SAVE :: o_alp_wk = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
464	'alp_wk', 'ALP WK', 'm2/s2', (/ ('', i=1, 9) /))
465	TYPE(ctrl_out), SAVE :: o_ale = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
466	'ale', 'ALE', 'm2/s2', (/ ('', i=1, 9) /))
467	TYPE(ctrl_out), SAVE :: o_alp = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
468	'alp', 'ALP', 'W/m2', (/ ('', i=1, 9) /))
469	TYPE(ctrl_out), SAVE :: o_cin = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
470	'cin', 'Convective INhibition', 'm2/s2', (/ ('', i=1, 9) /))
471	TYPE(ctrl_out), SAVE :: o_wape = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
472	'wape', '', '', (/ ('', i=1, 9) /))
473
474	!!! nrlmd le 10/04/2012
475
476	!-------Spectre de thermiques de type 2 au LCL
477	TYPE(ctrl_out), SAVE :: o_n2 = ctrl_out((/ 1, 1, 1, 6, 10, 10, 11, 11, 11 /), &
478	'n2', 'Nombre de panaches de type 2', ' ', (/ ('', i=1, 9) /))
479	TYPE(ctrl_out), SAVE :: o_s2 = ctrl_out((/ 1, 1, 1, 6, 10, 10, 11, 11, 11 /), &
480	's2', 'Surface moyenne des panaches de type 2', 'm2', (/ ('', i=1, 9) /))
481
482	!-------Déclenchement stochastique
483	TYPE(ctrl_out), SAVE :: o_proba_notrig = ctrl_out((/ 1, 1, 1, 6, 10, 10, 11, 11, 11 /), &
484	'proba_notrig', &
485	'Probabilité de non-déclenchement', ' ', (/ ('', i=1, 9) /))
486	TYPE(ctrl_out), SAVE :: o_random_notrig = ctrl_out((/ 1, 1, 1, 6, 10, 10, 11, 11, 11 /), &
487	'random_notrig', &
488	'Tirage aléatoire de non-déclenchement', ' ', (/ ('', i=1, 9) /))
489	TYPE(ctrl_out), SAVE :: o_ale_bl_stat = ctrl_out((/ 1, 1, 1, 6, 10, 10, 11, 11, 11 /), &
490	'ale_bl_stat', &
491	'ALE_BL_STAT', 'm2/s2', (/ ('', i=1, 9) /))
492	TYPE(ctrl_out), SAVE :: o_ale_bl_trig = ctrl_out((/ 1, 1, 1, 6, 10, 10, 11, 11, 11 /), &
493	'ale_bl_trig', &
494	'ALE_BL_STAT + Condition P>Pseuil', 'm2/s2', (/ ('', i=1, 9) /))
495
496	!-------Fermeture statistique
497	TYPE(ctrl_out), SAVE :: o_alp_bl_det = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
498	'alp_bl_det', 'ALP_BL_DET', 'W/m2', (/ ('', i=1, 9) /))
499	TYPE(ctrl_out), SAVE :: o_alp_bl_fluct_m = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
500	'alp_bl_fluct_m', 'ALP_BL_FLUCT_M', 'W/m2', (/ ('', i=1, 9) /))
501	TYPE(ctrl_out), SAVE :: o_alp_bl_fluct_tke = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
502	'alp_bl_fluct_tke', 'ALP_BL_FLUCT_TKE', 'W/m2', (/ ('', i=1, 9) /))
503	TYPE(ctrl_out), SAVE :: o_alp_bl_conv = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
504	'alp_bl_conv', 'ALP_BL_CONV', 'W/m2', (/ ('', i=1, 9) /))
505	TYPE(ctrl_out), SAVE :: o_alp_bl_stat = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
506	'alp_bl_stat', 'ALP_BL_STAT', 'W/m2', (/ ('', i=1, 9) /))
507
508	!!! fin nrlmd le 10/04/2012
509
510	! Champs interpolles sur des niveaux de pression ??? a faire correctement
511
512	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_uSTDlevs = (/ &
513	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u850', "Zonal wind 850hPa", "m/s", &
514	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
515	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u700', "Zonal wind 700hPa", "m/s", &
516	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
517	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u500', "Zonal wind 500hPa", "m/s", &
518	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
519	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u200', "Zonal wind 200hPa", "m/s", &
520	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
521	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u100', "Zonal wind 100hPa", "m/s", &
522	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
523	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u50', "Zonal wind 50hPa", "m/s", &
524	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
525	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u10', "Zonal wind 10hPa", "m/s", &
526	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
527
528	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_vSTDlevs = (/ &
529	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v850', "Meridional wind 850hPa", "m/s", &
530	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
531	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v700', "Meridional wind 700hPa", "m/s", &
532	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
533	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v500', "Meridional wind 500hPa", "m/s", &
534	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
535	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v200', "Meridional wind 200hPa", "m/s", &
536	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
537	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v100', "Meridional wind 100hPa", "m/s", &
538	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
539	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v50', "Meridional wind 50hPa", "m/s", &
540	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
541	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v10', "Meridional wind 10hPa", "m/s", &
542	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
543
544	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_wSTDlevs = (/ &
545	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w850', "Vertical wind 1hPa", "Pa/s", &
546	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
547	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w700', "Vertical wind 700hPa", "Pa/s", &
548	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
549	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w500', "Vertical wind 500hPa", "Pa/s", &
550	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
551	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w200', "Vertical wind 200hPa", "Pa/s", &
552	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
553	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w100', "Vertical wind 100hPa", "Pa/s", &
554	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
555	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w50', "Vertical wind 50hPa", "Pa/s", &
556	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
557	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w10', "Vertical wind 10hPa", "Pa/s", &
558	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
559
560	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_tSTDlevs = (/ &
561	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t850', "Temperature 1hPa", "K", &
562	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
563	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t700', "Temperature 700hPa", "K", &
564	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
565	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t500', "Temperature 500hPa", "K", &
566	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
567	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t200', "Temperature 200hPa", "K", &
568	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
569	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t100', "Temperature 100hPa", "K", &
570	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
571	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t50', "Temperature 50hPa", "K", &
572	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
573	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t10', "Temperature 10hPa", "K", &
574	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
575
576	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_qSTDlevs = (/ &
577	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q850', "Specific humidity 1hPa", &
578	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
579	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q700', "Specific humidity 700hPa", &
580	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
581	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q500', "Specific humidity 500hPa", &
582	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
583	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q200', "Specific humidity 200hPa", &
584	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
585	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q100', "Specific humidity 100hPa", &
586	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
587	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q50', "Specific humidity 50hPa", &
588	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
589	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q10', "Specific humidity 10hPa", &
590	"kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
591
592	TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_zSTDlevs = (/ &
593	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z850', "Geopotential height 1hPa", &
594	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
595	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z700', "Geopotential height 700hPa", &
596	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
597	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z500', "Geopotential height 500hPa", &
598	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
599	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z200', "Geopotential height 200hPa", &
600	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
601	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z100', "Geopotential height 100hPa", &
602	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
603	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z50', "Geopotential height 50hPa", &
604	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
605	ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z10', "Geopotential height 10hPa", &
606	"m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
607
608	TYPE(ctrl_out), SAVE :: o_t_oce_sic = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
609	't_oce_sic', 'Temp mixte oce-sic', 'K', (/ ('', i=1, 9) /))
610	TYPE(ctrl_out), SAVE :: o_weakinv = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
611	'weakinv', 'Weak inversion', '-', (/ ('', i=1, 9) /))
612	TYPE(ctrl_out), SAVE :: o_dthmin = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
613	'dthmin', 'dTheta mini', 'K/m', (/ ('', i=1, 9) /))
614
615	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_u10_srf = (/ &
616	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'u10_ter', "", "", (/ ('', i=1, 9) /)), &
617	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'u10_lic', "", "", (/ ('', i=1, 9) /)), &
618	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'u10_oce', "", "", (/ ('', i=1, 9) /)), &
619	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'u10_sic', "", "", (/ ('', i=1, 9) /)) /)
620
621	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_v10_srf = (/ &
622	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'v10_ter', "", "", (/ ('', i=1, 9) /)), &
623	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'v10_lic', "", "", (/ ('', i=1, 9) /)), &
624	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'v10_oce', "", "", (/ ('', i=1, 9) /)), &
625	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'v10_sic', "", "", (/ ('', i=1, 9) /)) /)
626
627	TYPE(ctrl_out), SAVE :: o_cldtau = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
628	'cldtau', 'Cloud optical thickness', '1', (/ ('', i=1, 9) /))
629	TYPE(ctrl_out), SAVE :: o_cldemi = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
630	'cldemi', 'Cloud optical emissivity', '1', (/ ('', i=1, 9) /))
631	TYPE(ctrl_out), SAVE :: o_rh2m = ctrl_out((/ 5, 5, 10, 10, 10, 10, 11, 11, 11 /), &
632	'rh2m', 'Relative humidity at 2m', '%', (/ ('', i=1, 9) /))
633	TYPE(ctrl_out), SAVE :: o_rh2m_min = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
634	'rh2m_min', 'Min Relative humidity at 2m', '%', &
635	(/ '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)' /))
636	TYPE(ctrl_out), SAVE :: o_rh2m_max = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
637	'rh2m_max', 'Max Relative humidity at 2m', '%', &
638	(/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
639	't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
640	TYPE(ctrl_out), SAVE :: o_qsat2m = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
641	'qsat2m', 'Saturant humidity at 2m', '%', (/ ('', i=1, 9) /))
642	TYPE(ctrl_out), SAVE :: o_tpot = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
643	'tpot', 'Surface air potential temperature', 'K', (/ ('', i=1, 9) /))
644	TYPE(ctrl_out), SAVE :: o_tpote = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11 /), &
645	'tpote', &
646	'Surface air equivalent potential temperature', 'K', (/ ('', i=1, 9) /))
647	TYPE(ctrl_out), SAVE :: o_tke = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
648	'tke ', 'TKE', 'm2/s2', (/ ('', i=1, 9) /))
649	TYPE(ctrl_out), SAVE :: o_tke_max = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
650	'tke_max', 'TKE max', 'm2/s2', &
651	(/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
652	't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
653
654	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tke_srf = (/ &
655	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_ter', &
656	"Max Turb. Kinetic Energy "//clnsurf(1),"-", (/ ('', i=1, 9) /)), &
657	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_lic', &
658	"Max Turb. Kinetic Energy "//clnsurf(2),"-", (/ ('', i=1, 9) /)), &
659	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_oce', &
660	"Max Turb. Kinetic Energy "//clnsurf(3),"-", (/ ('', i=1, 9) /)), &
661	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_sic', &
662	"Max Turb. Kinetic Energy "//clnsurf(4),"-", (/ ('', i=1, 9) /)) /)
663
664	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tke_max_srf = (/ &
665	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_max_ter', &
666	"Max Turb. Kinetic Energy "//clnsurf(1),"-", &
667	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
668	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /)), &
669	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_max_lic', &
670	"Max Turb. Kinetic Energy "//clnsurf(2),"-", &
671	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
672	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /)), &
673	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_max_oce', &
674	"Max Turb. Kinetic Energy "//clnsurf(3),"-", &
675	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
676	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /)), &
677	ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'tke_max_sic', &
678	"Max Turb. Kinetic Energy "//clnsurf(4),"-", &
679	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
680	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /)) /)
681
682	TYPE(ctrl_out), SAVE :: o_kz = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
683	'kz', 'Kz melange', 'm2/s', (/ ('', i=1, 9) /))
684	TYPE(ctrl_out), SAVE :: o_kz_max = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
685	'kz_max', 'Kz melange max', 'm2/s', &
686	(/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
687	't_max(X)', "t_max(X)", "t_max(X)", "t_max(X)" /))
688	TYPE(ctrl_out), SAVE :: o_SWnetOR = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
689	'SWnetOR', 'Sfce net SW radiation OR', 'W/m2', (/ ('', i=1, 9) /))
690	TYPE(ctrl_out), SAVE :: o_SWdownOR = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
691	'SWdownOR', 'Sfce incident SW radiation OR', 'W/m2', (/ ('', i=1, 9) /))
692	TYPE(ctrl_out), SAVE :: o_LWdownOR = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
693	'LWdownOR', 'Sfce incident LW radiation OR', 'W/m2', (/ ('', i=1, 9) /))
694	TYPE(ctrl_out), SAVE :: o_snowl = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
695	'snowl', 'Solid Large-scale Precip.', 'kg/(m2*s)', (/ ('', i=1, 9) /))
696	TYPE(ctrl_out), SAVE :: o_cape_max = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
697	'cape_max', 'CAPE max.', 'J/kg', &
698	(/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
699	"t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /))
700	TYPE(ctrl_out), SAVE :: o_solldown = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11 /), &
701	'solldown', 'Down. IR rad. at surface', 'W/m2', (/ ('', i=1, 9) /))
702	TYPE(ctrl_out), SAVE :: o_dtsvdfo = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
703	'dtsvdfo', 'Boundary-layer dTs(o)', 'K/s', (/ ('', i=1, 9) /))
704	TYPE(ctrl_out), SAVE :: o_dtsvdft = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
705	'dtsvdft', 'Boundary-layer dTs(t)', 'K/s', (/ ('', i=1, 9) /))
706	TYPE(ctrl_out), SAVE :: o_dtsvdfg = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
707	'dtsvdfg', 'Boundary-layer dTs(g)', 'K/s', (/ ('', i=1, 9) /))
708	TYPE(ctrl_out), SAVE :: o_dtsvdfi = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
709	'dtsvdfi', 'Boundary-layer dTs(g)', 'K/s', (/ ('', i=1, 9) /))
710	TYPE(ctrl_out), SAVE :: o_rugs = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
711	'rugs', 'rugosity', '-', (/ ('', i=1, 9) /))
712	TYPE(ctrl_out), SAVE :: o_topswad = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
713	'topswad', 'ADE at TOA', 'W/m2', (/ ('', i=1, 9) /))
714	TYPE(ctrl_out), SAVE :: o_topswad0 = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
715	'topswad0', 'ADE clear-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
716	TYPE(ctrl_out), SAVE :: o_topswai = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
717	'topswai', 'AIE at TOA', 'W/m2', (/ ('', i=1, 9) /))
718	TYPE(ctrl_out), SAVE :: o_solswad = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
719	'solswad', 'ADE at SRF', 'W/m2', (/ ('', i=1, 9) /))
720	TYPE(ctrl_out), SAVE :: o_solswad0 = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
721	'solswad0', 'ADE clear-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
722	TYPE(ctrl_out), SAVE :: o_solswai = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
723	'solswai', 'AIE at SFR', 'W/m2', (/ ('', i=1, 9) /))
724
725	type(ctrl_out),save,dimension(naero_tot) :: o_tausumaero = &
726	(/ ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_ASBCM', &
727	"Aerosol Optical depth at 550 nm "//name_aero_tau(1),"1", (/ ('', i=1, 9) /)), &
728	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_ASPOMM', &
729	"Aerosol Optical depth at 550 nm "//name_aero_tau(2),"1", (/ ('', i=1, 9) /)), &
730	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_ASSO4M', &
731	"Aerosol Optical depth at 550 nm "//name_aero_tau(3),"1", (/ ('', i=1, 9) /)), &
732	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_CSSO4M', &
733	"Aerosol Optical depth at 550 nm "//name_aero_tau(4),"1", (/ ('', i=1, 9) /)), &
734	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_SSSSM', &
735	"Aerosol Optical depth at 550 nm "//name_aero_tau(5),"1", (/ ('', i=1, 9) /)), &
736	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_ASSSM', &
737	"Aerosol Optical depth at 550 nm "//name_aero_tau(6),"1", (/ ('', i=1, 9) /)), &
738	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_CSSSM', &
739	"Aerosol Optical depth at 550 nm "//name_aero_tau(7),"1", (/ ('', i=1, 9) /)), &
740	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_CIDUSTM', &
741	"Aerosol Optical depth at 550 nm "//name_aero_tau(8),"1", (/ ('', i=1, 9) /)), &
742	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_AIBCM', &
743	"Aerosol Optical depth at 550 nm "//name_aero_tau(9),"1", (/ ('', i=1, 9) /)), &
744	ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /),'OD550_AIPOMM', &
745	"Aerosol Optical depth at 550 nm "//name_aero_tau(10),"1", (/ ('', i=1, 9) /)),&
746	ctrl_out((/ 2, 2, 10, 10, 10, 10, 11, 11, 11 /),'OD550_STRAT', &
747	"Aerosol Optical depth at 550 nm "//name_aero_tau(11),"1", (/ ('', i=1, 9) /)) /)
748	!
749	TYPE(ctrl_out), SAVE :: o_od550aer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
750	'od550aer', 'Total aerosol optical depth at 550nm', '-', (/ ('', i=1, 9) /))
751	TYPE(ctrl_out), SAVE :: o_od865aer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
752	'od865aer', 'Total aerosol optical depth at 870nm', '-', (/ ('', i=1, 9) /))
753	TYPE(ctrl_out), SAVE :: o_absvisaer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
754	'absvisaer', 'Absorption aerosol visible optical depth', '-', (/ ('', i=1, 9) /))
755	TYPE(ctrl_out), SAVE :: o_od550lt1aer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
756	'od550lt1aer', 'Fine mode optical depth', '-', (/ ('', i=1, 9) /))
757	TYPE(ctrl_out), SAVE :: o_sconcso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
758	'sconcso4', 'Surface Concentration of Sulfate ', 'kg/m3', (/ ('', i=1, 9) /))
759	TYPE(ctrl_out), SAVE :: o_sconcoa = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
760	'sconcoa', 'Surface Concentration of Organic Aerosol ', 'kg/m3', (/ ('', i=1, 9) /))
761	TYPE(ctrl_out), SAVE :: o_sconcbc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
762	'sconcbc', 'Surface Concentration of Black Carbon ', 'kg/m3', (/ ('', i=1, 9) /))
763	TYPE(ctrl_out), SAVE :: o_sconcss = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
764	'sconcss', 'Surface Concentration of Sea Salt ', 'kg/m3', (/ ('', i=1, 9) /))
765	TYPE(ctrl_out), SAVE :: o_sconcdust = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
766	'sconcdust', 'Surface Concentration of Dust ', 'kg/m3', (/ ('', i=1, 9) /))
767	TYPE(ctrl_out), SAVE :: o_concso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
768	'concso4', 'Concentration of Sulfate ', 'kg/m3', (/ ('', i=1, 9) /))
769	TYPE(ctrl_out), SAVE :: o_concoa = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
770	'concoa', 'Concentration of Organic Aerosol ', 'kg/m3', (/ ('', i=1, 9) /))
771	TYPE(ctrl_out), SAVE :: o_concbc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
772	'concbc', 'Concentration of Black Carbon ', 'kg/m3', (/ ('', i=1, 9) /))
773	TYPE(ctrl_out), SAVE :: o_concss = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
774	'concss', 'Concentration of Sea Salt ', 'kg/m3', (/ ('', i=1, 9) /))
775	TYPE(ctrl_out), SAVE :: o_concdust = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
776	'concdust', 'Concentration of Dust ', 'kg/m3', (/ ('', i=1, 9) /))
777	TYPE(ctrl_out), SAVE :: o_loadso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
778	'loadso4', 'Column Load of Sulfate ', 'kg/m2', (/ ('', i=1, 9) /))
779	TYPE(ctrl_out), SAVE :: o_loadoa = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
780	'loadoa', 'Column Load of Organic Aerosol ', 'kg/m2', (/ ('', i=1, 9) /))
781	TYPE(ctrl_out), SAVE :: o_loadbc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
782	'loadbc', 'Column Load of Black Carbon ', 'kg/m2', (/ ('', i=1, 9) /))
783	TYPE(ctrl_out), SAVE :: o_loadss = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
784	'loadss', 'Column Load of Sea Salt ', 'kg/m2', (/ ('', i=1, 9) /))
785	TYPE(ctrl_out), SAVE :: o_loaddust = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
786	'loaddust', 'Column Load of Dust ', 'kg/m2', (/ ('', i=1, 9) /))
787	TYPE(ctrl_out), SAVE :: o_swtoaas_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
788	'swtoaas_nat', 'Natural aerosol radiative forcing all-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
789	TYPE(ctrl_out), SAVE :: o_swsrfas_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
790	'swsrfas_nat', 'Natural aerosol radiative forcing all-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
791	TYPE(ctrl_out), SAVE :: o_swtoacs_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
792	'swtoacs_nat', 'Natural aerosol radiative forcing clear-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
793	TYPE(ctrl_out), SAVE :: o_swsrfcs_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
794	'swsrfcs_nat', 'Natural aerosol radiative forcing clear-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
795	TYPE(ctrl_out), SAVE :: o_swtoaas_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
796	'swtoaas_ant', 'Anthropogenic aerosol radiative forcing all-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
797	TYPE(ctrl_out), SAVE :: o_swsrfas_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
798	'swsrfas_ant', 'Anthropogenic aerosol radiative forcing all-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
799	TYPE(ctrl_out), SAVE :: o_swtoacs_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
800	'swtoacs_ant', 'Anthropogenic aerosol radiative forcing clear-sky at TOA', 'W/m2', (/ ('', i=1, 9) /))
801	TYPE(ctrl_out), SAVE :: o_swsrfcs_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
802	'swsrfcs_ant', 'Anthropogenic aerosol radiative forcing clear-sky at SRF', 'W/m2', (/ ('', i=1, 9) /))
803	TYPE(ctrl_out), SAVE :: o_swtoacf_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
804	'swtoacf_nat', 'Natural aerosol impact on cloud radiative forcing at TOA', 'W/m2', (/ ('', i=1, 9) /))
805	TYPE(ctrl_out), SAVE :: o_swsrfcf_nat = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
806	'swsrfcf_nat', 'Natural aerosol impact on cloud radiative forcing at SRF', 'W/m2', (/ ('', i=1, 9) /))
807	TYPE(ctrl_out), SAVE :: o_swtoacf_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
808	'swtoacf_ant', 'Anthropogenic aerosol impact on cloud radiative forcing at TOA', 'W/m2', (/ ('', i=1, 9) /))
809	TYPE(ctrl_out), SAVE :: o_swsrfcf_ant = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
810	'swsrfcf_ant', 'Anthropogenic aerosol impact on cloud radiative forcing at SRF', 'W/m2', (/ ('', i=1, 9) /))
811	TYPE(ctrl_out), SAVE :: o_swtoacf_zero = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
812	'swtoacf_zero', 'Cloud radiative forcing (allsky-clearsky fluxes) at TOA', 'W/m2', (/ ('', i=1, 9) /))
813	TYPE(ctrl_out), SAVE :: o_swsrfcf_zero = ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11 /), &
814	'swsrfcf_zero', 'Cloud radiative forcing (allsky-clearsky fluxes) at SRF', 'W/m2', (/ ('', i=1, 9) /))
815	TYPE(ctrl_out), SAVE :: o_cldncl = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
816	'cldncl', 'CDNC at top of liquid water cloud', 'm-3', (/ ('', i=1, 9) /))
817	TYPE(ctrl_out), SAVE :: o_reffclwtop = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
818	'reffclwtop', 'Droplet effective radius at top of liquid water cloud', 'm', (/ ('', i=1, 9) /))
819	TYPE(ctrl_out), SAVE :: o_cldnvi = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
820	'cldnvi', 'Column Integrated Cloud Droplet Number', 'm-2', (/ ('', i=1, 9) /))
821	TYPE(ctrl_out), SAVE :: o_lcc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
822	'lcc', 'Cloud liquid fraction at top of cloud', '1', (/ ('', i=1, 9) /))
823
824
825	!!!!!!!!!!!!!!!!!!!!!! 3D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
826	TYPE(ctrl_out), SAVE :: o_ec550aer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
827	'ec550aer', 'Extinction at 550nm', 'm^-1', (/ ('', i=1, 9) /))
828	TYPE(ctrl_out), SAVE :: o_lwcon = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
829	'lwcon', 'Cloud liquid water content', 'kg/kg', (/ ('', i=1, 9) /))
830	TYPE(ctrl_out), SAVE :: o_iwcon = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
831	'iwcon', 'Cloud ice water content', 'kg/kg', (/ ('', i=1, 9) /))
832	TYPE(ctrl_out), SAVE :: o_temp = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11 /), &
833	'temp', 'Air temperature', 'K', (/ ('', i=1, 9) /))
834	TYPE(ctrl_out), SAVE :: o_theta = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11 /), &
835	'theta', 'Potential air temperature', 'K', (/ ('', i=1, 9) /))
836	TYPE(ctrl_out), SAVE :: o_ovap = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11 /), &
837	'ovap', 'Specific humidity', 'kg/kg', (/ ('', i=1, 9) /))
838	TYPE(ctrl_out), SAVE :: o_ovapinit = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
839	'ovapinit', 'Specific humidity (begin of timestep)', 'kg/kg', (/ ('', i=1, 9) /))
840	TYPE(ctrl_out), SAVE :: o_oliq = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11 /), &
841	'oliq', 'Condensed water', 'kg/kg', (/ ('', i=1, 9) /))
842	TYPE(ctrl_out), SAVE :: o_wvapp = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
843	'wvapp', '', '', (/ ('', i=1, 9) /))
844	TYPE(ctrl_out), SAVE :: o_geop = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
845	'geop', 'Geopotential height', 'm2/s2', (/ ('', i=1, 9) /))
846	TYPE(ctrl_out), SAVE :: o_vitu = ctrl_out((/ 2, 3, 4, 6, 10, 10, 11, 11, 11 /), &
847	'vitu', 'Zonal wind', 'm/s', (/ ('', i=1, 9) /))
848	TYPE(ctrl_out), SAVE :: o_vitv = ctrl_out((/ 2, 3, 4, 6, 10, 10, 11, 11, 11 /), &
849	'vitv', 'Meridional wind', 'm/s', (/ ('', i=1, 9) /))
850	TYPE(ctrl_out), SAVE :: o_vitw = ctrl_out((/ 2, 3, 10, 6, 10, 10, 11, 11, 11 /), &
851	'vitw', 'Vertical wind', 'Pa/s', (/ ('', i=1, 9) /))
852	TYPE(ctrl_out), SAVE :: o_pres = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
853	'pres', 'Air pressure', 'Pa', (/ ('', i=1, 9) /))
854	TYPE(ctrl_out), SAVE :: o_paprs = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
855	'paprs', 'Air pressure Inter-Couches', 'Pa', (/ ('', i=1, 9) /))
856	TYPE(ctrl_out), SAVE :: o_mass = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
857	'mass', 'Masse Couches', 'kg/m2', (/ ('', i=1, 9) /))
858	TYPE(ctrl_out), SAVE :: o_zfull = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
859	'zfull', 'Altitude of full pressure levels', 'm', (/ ('', i=1, 9) /))
860	TYPE(ctrl_out), SAVE :: o_zhalf = ctrl_out((/ 2, 3, 10, 10, 10, 10, 11, 11, 11 /), &
861	'zhalf', 'Altitude of half pressure levels', 'm', (/ ('', i=1, 9) /))
862	TYPE(ctrl_out), SAVE :: o_rneb = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
863	'rneb', 'Cloud fraction', '-', (/ ('', i=1, 9) /))
864	TYPE(ctrl_out), SAVE :: o_rnebjn = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11,11 /), &
865	'rnebjn', 'Cloud fraction in day', '-', (/ ('', i=1, 9) /))
866	TYPE(ctrl_out), SAVE :: o_rnebcon = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
867	'rnebcon', 'Convective Cloud Fraction', '-', (/ ('', i=1, 9) /))
868	TYPE(ctrl_out), SAVE :: o_rnebls = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
869	'rnebls', 'LS Cloud fraction', '-', (/ ('', i=1, 9) /))
870	TYPE(ctrl_out), SAVE :: o_rhum = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
871	'rhum', 'Relative humidity', '-', (/ ('', i=1, 9) /))
872	TYPE(ctrl_out), SAVE :: o_ozone = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
873	'ozone', 'Ozone mole fraction', '-', (/ ('', i=1, 9) /))
874	TYPE(ctrl_out), SAVE :: o_ozone_light = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
875	'ozone_daylight', 'Daylight ozone mole fraction', '-', (/ ('', i=1, 9) /))
876	TYPE(ctrl_out), SAVE :: o_upwd = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
877	'upwd', 'saturated updraft', 'kg/m2/s', (/ ('', i=1, 9) /))
878	TYPE(ctrl_out), SAVE :: o_dtphy = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
879	'dtphy', 'Physics dT', 'K/s', (/ ('', i=1, 9) /))
880	TYPE(ctrl_out), SAVE :: o_dqphy = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
881	'dqphy', 'Physics dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
882	TYPE(ctrl_out), SAVE :: o_pr_con_l = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
883	'pr_con_l', 'Convective precipitation lic', ' ', (/ ('', i=1, 9) /))
884	TYPE(ctrl_out), SAVE :: o_pr_con_i = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
885	'pr_con_i', 'Convective precipitation ice', ' ', (/ ('', i=1, 9) /))
886	TYPE(ctrl_out), SAVE :: o_pr_lsc_l = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
887	'pr_lsc_l', 'Large scale precipitation lic', ' ', (/ ('', i=1, 9) /))
888	TYPE(ctrl_out), SAVE :: o_pr_lsc_i = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11 /), &
889	'pr_lsc_i', 'Large scale precipitation ice', ' ', (/ ('', i=1, 9) /))
890	TYPE(ctrl_out), SAVE :: o_re = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
891	're', 'Cloud droplet effective radius', 'um', (/ ('', i=1, 9) /))
892	TYPE(ctrl_out), SAVE :: o_fl = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
893	'fl', 'Denominator of Cloud droplet effective radius', ' ', (/ ('', i=1, 9) /))
894	TYPE(ctrl_out), SAVE :: o_scdnc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
895	'scdnc', 'Cloud droplet number concentration', 'm-3', (/ ('', i=1, 9) /))
896	TYPE(ctrl_out), SAVE :: o_reffclws = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
897	'reffclws', 'Stratiform Cloud Droplet Effective Radius (aerosol diags.)', 'm', (/ ('', i=1, 9) /))
898	TYPE(ctrl_out), SAVE :: o_reffclwc = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
899	'reffclwc', 'Convective Cloud Droplet Effective Radius (aerosol diags.)', 'm', (/ ('', i=1, 9) /))
900	TYPE(ctrl_out), SAVE :: o_lcc3d = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
901	'lcc3d', 'Cloud liquid fraction', '1', (/ ('', i=1, 9) /))
902	TYPE(ctrl_out), SAVE :: o_lcc3dcon = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
903	'lcc3dcon', 'Convective cloud liquid fraction', '1', (/ ('', i=1, 9) /))
904	TYPE(ctrl_out), SAVE :: o_lcc3dstra = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11 /), &
905	'lcc3dstra', 'Stratiform cloud liquid fraction', '1', (/ ('', i=1, 9) /))
906	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
907
908	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_albe_srf = (/ &
909	ctrl_out((/ 3, 7, 10, 7, 10, 10, 11, 11, 11 /),'albe_ter', "Albedo VIS surf. "//clnsurf(1),"-", (/ ('', i=1, 9) /)), &
910	ctrl_out((/ 3, 7, 10, 7, 10, 10, 11, 11, 11 /),'albe_lic', "Albedo VIS surf. "//clnsurf(2),"-", (/ ('', i=1, 9) /)), &
911	ctrl_out((/ 3, 7, 10, 7, 10, 10, 11, 11, 11 /),'albe_oce', "Albedo VIS surf. "//clnsurf(3),"-", (/ ('', i=1, 9) /)), &
912	ctrl_out((/ 3, 7, 10, 7, 10, 10, 11, 11, 11 /),'albe_sic', "Albedo VIS surf. "//clnsurf(4),"-", (/ ('', i=1, 9) /)) /)
913
914	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_ages_srf = (/ &
915	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'ages_ter', "Snow age", "day", (/ ('', i=1, 9) /)), &
916	ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /),'ages_lic', "Snow age", "day", (/ ('', i=1, 9) /)), &
917	ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /),'ages_oce',"Snow age", "day", (/ ('', i=1, 9) /)), &
918	ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /),'ages_sic',"Snow age", "day", (/ ('', i=1, 9) /)) /)
919
920	TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_rugs_srf = (/ &
921	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'rugs_ter', "Surface roughness "//clnsurf(1),"m", (/ ('', i=1, 9) /)), &
922	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'rugs_lic', "Surface roughness "//clnsurf(2),"m", (/ ('', i=1, 9) /)), &
923	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'rugs_oce', "Surface roughness "//clnsurf(3),"m", (/ ('', i=1, 9) /)), &
924	ctrl_out((/ 3, 6, 10, 10, 10, 10, 11, 11, 11 /),'rugs_sic', "Surface roughness "//clnsurf(4),"m", (/ ('', i=1, 9) /)) /)
925
926	TYPE(ctrl_out), SAVE :: o_alb1 = ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /), &
927	'alb1', 'Surface VIS albedo', '-', (/ ('', i=1, 9) /))
928	TYPE(ctrl_out), SAVE :: o_alb2 = ctrl_out((/ 3, 10, 10, 10, 10, 10, 11, 11, 11 /), &
929	'alb2', 'Surface Near IR albedo', '-', (/ ('', i=1, 9) /))
930	TYPE(ctrl_out), SAVE :: o_clwcon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
931	'clwcon', 'Convective Cloud Liquid water content', 'kg/kg', (/ ('', i=1, 9) /))
932	TYPE(ctrl_out), SAVE :: o_Ma = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
933	'Ma', 'undilute adiab updraft', 'kg/m2/s', (/ ('', i=1, 9) /))
934	TYPE(ctrl_out), SAVE :: o_dnwd = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
935	'dnwd', 'saturated downdraft', 'kg/m2/s', (/ ('', i=1, 9) /))
936	TYPE(ctrl_out), SAVE :: o_dnwd0 = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
937	'dnwd0', 'unsat. downdraft', 'kg/m2/s', (/ ('', i=1, 9) /))
938	TYPE(ctrl_out), SAVE :: o_mc = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
939	'mc', 'Convective mass flux', 'kg/m2/s', (/ ('', i=1, 9) /))
940	TYPE(ctrl_out), SAVE :: o_ftime_con = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
941	'ftime_con', 'Fraction of time convection Occurs', ' ', &
942	(/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /))
943	TYPE(ctrl_out), SAVE :: o_dtdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
944	'dtdyn', 'Dynamics dT', 'K/s', (/ ('', i=1, 9) /))
945	TYPE(ctrl_out), SAVE :: o_dqdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
946	'dqdyn', 'Dynamics dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
947	TYPE(ctrl_out), SAVE :: o_dudyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
948	'dudyn', 'Dynamics dU', 'm/s2', (/ ('', i=1, 9) /))
949	TYPE(ctrl_out), SAVE :: o_dvdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
950	'dvdyn', 'Dynamics dV', 'm/s2', (/ ('', i=1, 9) /))
951	TYPE(ctrl_out), SAVE :: o_dtcon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
952	'dtcon', 'Convection dT', 'K/s', (/ ('', i=1, 9) /))
953	TYPE(ctrl_out), SAVE :: o_ducon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
954	'ducon', 'Convection du', 'm/s2', (/ ('', i=1, 9) /))
955	TYPE(ctrl_out), SAVE :: o_dvcon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
956	'dvcon', 'Convection dv', 'm/s2', (/ ('', i=1, 9) /))
957	TYPE(ctrl_out), SAVE :: o_dqcon = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
958	'dqcon', 'Convection dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
959	TYPE(ctrl_out), SAVE :: o_dtwak = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
960	'dtwak', 'Wake dT', 'K/s', (/ ('', i=1, 9) /))
961	TYPE(ctrl_out), SAVE :: o_dqwak = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
962	'dqwak', 'Wake dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
963	TYPE(ctrl_out), SAVE :: o_wake_h = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
964	'wake_h', 'wake_h', '-', (/ ('', i=1, 9) /))
965	TYPE(ctrl_out), SAVE :: o_wake_s = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
966	'wake_s', 'wake_s', '-', (/ ('', i=1, 9) /))
967	TYPE(ctrl_out), SAVE :: o_wake_deltat = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
968	'wake_deltat', 'wake_deltat', ' ', (/ ('', i=1, 9) /))
969	TYPE(ctrl_out), SAVE :: o_wake_deltaq = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
970	'wake_deltaq', 'wake_deltaq', ' ', (/ ('', i=1, 9) /))
971	TYPE(ctrl_out), SAVE :: o_wake_omg = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
972	'wake_omg', 'wake_omg', ' ', (/ ('', i=1, 9) /))
973	TYPE(ctrl_out), SAVE :: o_wdtrainA = ctrl_out((/ 4, 1, 10, 4, 1, 10, 11, 11, 110 /), &
974	'wdtrainA', 'precipitation from AA', '-', (/ ('', i=1, 9) /))
975	TYPE(ctrl_out), SAVE :: o_wdtrainM = ctrl_out((/ 4, 1, 10, 4, 1, 10, 11, 11, 110 /), &
976	'wdtrainM', 'precipitation from mixture', '-', (/ ('', i=1, 9) /))
977	TYPE(ctrl_out), SAVE :: o_Vprecip = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
978	'Vprecip', 'precipitation vertical profile', '-', (/ ('', i=1, 9) /))
979	TYPE(ctrl_out), SAVE :: o_ftd = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
980	'ftd', 'tend temp due aux descentes precip', '-', (/ ('', i=1, 9) /))
981	TYPE(ctrl_out), SAVE :: o_fqd = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11 /), &
982	'fqd', 'tend vap eau due aux descentes precip', '-', (/ ('', i=1, 9) /))
983	TYPE(ctrl_out), SAVE :: o_dtlsc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
984	'dtlsc', 'Condensation dT', 'K/s', (/ ('', i=1, 9) /))
985	TYPE(ctrl_out), SAVE :: o_dtlschr = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
986	'dtlschr', 'Large-scale condensational heating rate', 'K/s', (/ ('', i=1, 9) /))
987	TYPE(ctrl_out), SAVE :: o_dqlsc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
988	'dqlsc', 'Condensation dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
989	TYPE(ctrl_out), SAVE :: o_beta_prec = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
990	'beta_prec', 'LS Conversion rate to prec', '(kg/kg)/s', (/ ('', i=1, 9) /))
991	TYPE(ctrl_out), SAVE :: o_dtvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
992	'dtvdf', 'Boundary-layer dT', 'K/s', (/ ('', i=1, 9) /))
993	TYPE(ctrl_out), SAVE :: o_dtdis = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
994	'dtdis', 'TKE dissipation dT', 'K/s', (/ ('', i=1, 9) /))
995	TYPE(ctrl_out), SAVE :: o_dqvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
996	'dqvdf', 'Boundary-layer dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
997	TYPE(ctrl_out), SAVE :: o_dteva = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
998	'dteva', 'Reevaporation dT', 'K/s', (/ ('', i=1, 9) /))
999	TYPE(ctrl_out), SAVE :: o_dqeva = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1000	'dqeva', 'Reevaporation dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1001
1002	!!!!!!!!!!!!!!!! Specifique thermiques
1003	TYPE(ctrl_out), SAVE :: o_dqlscth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1004	'dqlscth', 'dQ therm.', '(kg/kg)/s', (/ ('', i=1, 9) /))
1005	TYPE(ctrl_out), SAVE :: o_dqlscst = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1006	'dqlscst', 'dQ strat.', '(kg/kg)/s', (/ ('', i=1, 9) /))
1007	TYPE(ctrl_out), SAVE :: o_dtlscth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1008	'dtlscth', 'dQ therm.', 'K/s', (/ ('', i=1, 9) /))
1009	TYPE(ctrl_out), SAVE :: o_dtlscst = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1010	'dtlscst', 'dQ strat.', 'K/s', (/ ('', i=1, 9) /))
1011	TYPE(ctrl_out), SAVE :: o_plulth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1012	'plulth', 'Rainfall therm.', 'K/s', (/ ('', i=1, 9) /))
1013	TYPE(ctrl_out), SAVE :: o_plulst = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1014	'plulst', 'Rainfall strat.', 'K/s', (/ ('', i=1, 9) /))
1015	TYPE(ctrl_out), SAVE :: o_lmaxth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1016	'lmaxth', "Upper level thermals", "", (/ ('', i=1, 9) /))
1017	TYPE(ctrl_out), SAVE :: o_ptconvth = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1018	'ptconvth', 'POINTS CONVECTIFS therm.', ' ', (/ ('', i=1, 9) /))
1019	!!!!!!!!!!!!!!!!!!!!!!!!
1020	TYPE(ctrl_out), SAVE :: o_ptconv = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1021	'ptconv', 'POINTS CONVECTIFS', ' ', (/ ('', i=1, 9) /))
1022	TYPE(ctrl_out), SAVE :: o_ratqs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1023	'ratqs', 'RATQS', ' ', (/ ('', i=1, 9) /))
1024	TYPE(ctrl_out), SAVE :: o_dtthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1025	'dtthe', 'Thermal dT', 'K/s', (/ ('', i=1, 9) /))
1026	TYPE(ctrl_out), SAVE :: o_duthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1027	'duthe', 'Thermal du', 'm/s2', (/ ('', i=1, 9) /))
1028	TYPE(ctrl_out), SAVE :: o_dvthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1029	'dvthe', 'Thermal dv', 'm/s2', (/ ('', i=1, 9) /))
1030	TYPE(ctrl_out), SAVE :: o_f_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1031	'f_th', 'Thermal plume mass flux', 'kg/(m2*s)', (/ ('', i=1, 9) /))
1032	TYPE(ctrl_out), SAVE :: o_e_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1033	'e_th', 'Thermal plume entrainment', 'K/s', (/ ('', i=1, 9) /))
1034	TYPE(ctrl_out), SAVE :: o_w_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1035	'w_th', 'Thermal plume vertical velocity', 'm/s', (/ ('', i=1, 9) /))
1036	TYPE(ctrl_out), SAVE :: o_lambda_th = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1037	'lambda_th', 'Thermal plume vertical velocity', 'm/s', (/ ('', i=1, 9) /))
1038	TYPE(ctrl_out), SAVE :: o_ftime_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1039	'ftime_th', 'Fraction of time Shallow convection occurs', ' ', (/ ('', i=1, 9) /))
1040	TYPE(ctrl_out), SAVE :: o_q_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1041	'q_th', 'Thermal plume total humidity', 'kg/kg', (/ ('', i=1, 9) /))
1042	TYPE(ctrl_out), SAVE :: o_a_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1043	'a_th', "Thermal plume fraction", "", (/ ('', i=1, 9) /))
1044	TYPE(ctrl_out), SAVE :: o_d_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1045	'd_th', 'Thermal plume detrainment', 'K/s', (/ ('', i=1, 9) /))
1046	TYPE(ctrl_out), SAVE :: o_f0_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1047	'f0_th', 'Thermal closure mass flux', 'K/s', (/ ('', i=1, 9) /))
1048	TYPE(ctrl_out), SAVE :: o_zmax_th = ctrl_out((/ 4, 4, 4, 5, 10, 10, 11, 11, 11 /), &
1049	'zmax_th', 'Thermal plume height', 'K/s', (/ ('', i=1, 9) /))
1050	TYPE(ctrl_out), SAVE :: o_dqthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1051	'dqthe', 'Thermal dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1052	TYPE(ctrl_out), SAVE :: o_dtajs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1053	'dtajs', 'Dry adjust. dT', 'K/s', (/ ('', i=1, 9) /))
1054	TYPE(ctrl_out), SAVE :: o_dqajs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1055	'dqajs', 'Dry adjust. dQ', '(kg/kg)/s', (/ ('', i=1, 9) /))
1056	TYPE(ctrl_out), SAVE :: o_dtswr = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1057	'dtswr', 'SW radiation dT', 'K/s', (/ ('', i=1, 9) /))
1058	TYPE(ctrl_out), SAVE :: o_dtsw0 = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1059	'dtsw0', 'CS SW radiation dT', 'K/s', (/ ('', i=1, 9) /))
1060	TYPE(ctrl_out), SAVE :: o_dtlwr = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1061	'dtlwr', 'LW radiation dT', 'K/s', (/ ('', i=1, 9) /))
1062	TYPE(ctrl_out), SAVE :: o_dtlw0 = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1063	'dtlw0', 'CS LW radiation dT', 'K/s', (/ ('', i=1, 9) /))
1064	TYPE(ctrl_out), SAVE :: o_dtec = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1065	'dtec', 'Cinetic dissip dT', 'K/s', (/ ('', i=1, 9) /))
1066	TYPE(ctrl_out), SAVE :: o_duvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1067	'duvdf', 'Boundary-layer dU', 'm/s2', (/ ('', i=1, 9) /))
1068	TYPE(ctrl_out), SAVE :: o_dvvdf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1069	'dvvdf', 'Boundary-layer dV', 'm/s2', (/ ('', i=1, 9) /))
1070	TYPE(ctrl_out), SAVE :: o_duoro = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1071	'duoro', 'Orography dU', 'm/s2', (/ ('', i=1, 9) /))
1072	TYPE(ctrl_out), SAVE :: o_dvoro = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1073	'dvoro', 'Orography dV', 'm/s2', (/ ('', i=1, 9) /))
1074	TYPE(ctrl_out), SAVE :: o_dulif = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1075	'dulif', 'Orography dU', 'm/s2', (/ ('', i=1, 9) /))
1076	TYPE(ctrl_out), SAVE :: o_dvlif = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1077	'dvlif', 'Orography dV', 'm/s2', (/ ('', i=1, 9) /))
1078	TYPE(ctrl_out), SAVE :: o_duhin = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1079	'duhin', 'Hines GWD dU', 'm/s2', (/ ('', i=1, 9) /))
1080	TYPE(ctrl_out), SAVE :: o_dvhin = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1081	'dvhin', 'Hines GWD dV', 'm/s2', (/ ('', i=1, 9) /))
1082	TYPE(ctrl_out), SAVE :: o_dtoro = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1083	'dtoro', 'Orography dT', 'K/s', (/ ('', i=1, 9) /))
1084	TYPE(ctrl_out), SAVE :: o_dtlif = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1085	'dtlif', 'Orography dT', 'K/s', (/ ('', i=1, 9) /))
1086	TYPE(ctrl_out), SAVE :: o_dthin = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1087	'dthin', 'Hines GWD dT', 'K/s', (/ ('', i=1, 9) /))
1088
1089	type(ctrl_out), save:: o_du_gwd_rando &
1090	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'du_gwd_rando', &
1091	"Random gravity waves dU/dt", "m/s2", (/ ('', i=1, 9) /))
1092	type(ctrl_out), save:: o_dv_gwd_rando &
1093	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'dv_gwd_rando', &
1094	"Random gravity waves dV/dt", "m/s2", (/ ('', i=1, 9) /))
1095	type(ctrl_out), save:: o_vstr_gwd_rando &
1096	= ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'vstr_gwd_rando', &
1097	"meridional wind stress random gravity waves", "Pa", (/ ('', i=1, 9) /))
1098
1099	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac(:)
1100	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac_cum(:)
1101	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_vdf(:)
1102	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_the(:)
1103	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_con(:)
1104	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_lessi_impa(:)
1105	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_lessi_nucl(:)
1106	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_insc(:)
1107	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_bcscav(:)
1108	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_evapls(:)
1109	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_ls(:)
1110	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_trsp(:)
1111	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_sscav(:)
1112	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_sat(:)
1113	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_uscav(:)
1114	TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_dry(:)
1115
1116	TYPE(ctrl_out), SAVE :: o_rsu = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1117	'rsu', 'SW upward radiation', 'W m-2', (/ ('', i=1, 9) /))
1118	TYPE(ctrl_out), SAVE :: o_rsd = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1119	'rsd', 'SW downward radiation', 'W m-2', (/ ('', i=1, 9) /))
1120	TYPE(ctrl_out), SAVE :: o_rlu = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1121	'rlu', 'LW upward radiation', 'W m-2', (/ ('', i=1, 9) /))
1122	TYPE(ctrl_out), SAVE :: o_rld = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1123	'rld', 'LW downward radiation', 'W m-2', (/ ('', i=1, 9) /))
1124	TYPE(ctrl_out), SAVE :: o_rsucs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1125	'rsucs', 'SW CS upward radiation', 'W m-2', (/ ('', i=1, 9) /))
1126	TYPE(ctrl_out), SAVE :: o_rsdcs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1127	'rsdcs', 'SW CS downward radiation', 'W m-2', (/ ('', i=1, 9) /))
1128	TYPE(ctrl_out), SAVE :: o_rlucs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1129	'rlucs', 'LW CS upward radiation', 'W m-2', (/ ('', i=1, 9) /))
1130	TYPE(ctrl_out), SAVE :: o_rldcs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1131	'rldcs', 'LW CS downward radiation', 'W m-2', (/ ('', i=1, 9) /))
1132	TYPE(ctrl_out), SAVE :: o_tnt = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1133	'tnt', 'Tendency of air temperature', 'K s-1', (/ ('', i=1, 9) /))
1134	TYPE(ctrl_out), SAVE :: o_tntc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1135	'tntc', 'Tendency of air temperature due to Moist Convection', 'K s-1', (/ ('', i=1, 9) /))
1136	TYPE(ctrl_out), SAVE :: o_tntr = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1137	'tntr', 'Air temperature tendency due to Radiative heating', 'K s-1', (/ ('', i=1, 9) /))
1138	TYPE(ctrl_out), SAVE :: o_tntscpbl = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1139	'tntscpbl', 'Air temperature tendency due to St cloud and precipitation and BL mixing', &
1140	'K s-1', (/ ('', i=1, 9) /))
1141	TYPE(ctrl_out), SAVE :: o_tnhus = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1142	'tnhus', 'Tendency of specific humidity', 's-1', (/ ('', i=1, 9) /))
1143	TYPE(ctrl_out), SAVE :: o_tnhusc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1144	'tnhusc', 'Tendency of specific humidity due to convection', 's-1', (/ ('', i=1, 9) /))
1145	TYPE(ctrl_out), SAVE :: o_tnhusscpbl = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1146	'tnhusscpbl', 'Tendency of Specific humidity due to ST cl, precip and BL mixing', 's-1', (/ ('', i=1, 9) /))
1147	TYPE(ctrl_out), SAVE :: o_evu = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1148	'evu', 'Eddy viscosity coefficient for Momentum Variables', 'm2 s-1', (/ ('', i=1, 9) /))
1149	TYPE(ctrl_out), SAVE :: o_h2o = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1150	'h2o', 'Mass Fraction of Water', '1', (/ ('', i=1, 9) /))
1151	TYPE(ctrl_out), SAVE :: o_mcd = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1152	'mcd', 'Downdraft COnvective Mass Flux', 'kg/(m2*s)', (/ ('', i=1, 9) /))
1153	TYPE(ctrl_out), SAVE :: o_dmc = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1154	'dmc', 'Deep COnvective Mass Flux', 'kg/(m2*s)', (/ ('', i=1, 9) /))
1155	TYPE(ctrl_out), SAVE :: o_ref_liq = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1156	'ref_liq', 'Effective radius of convective cloud liquid water particle', 'm', (/ ('', i=1, 9) /))
1157	TYPE(ctrl_out), SAVE :: o_ref_ice = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1158	'ref_ice', 'Effective radius of startiform cloud ice particle', 'm', (/ ('', i=1, 9) /))
1159	TYPE(ctrl_out), SAVE :: o_rsut4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1160	'rsut4co2', 'TOA Out SW in 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1161	TYPE(ctrl_out), SAVE :: o_rlut4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1162	'rlut4co2', 'TOA Out LW in 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1163	TYPE(ctrl_out), SAVE :: o_rsutcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1164	'rsutcs4co2', 'TOA Out CS SW in 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1165	TYPE(ctrl_out), SAVE :: o_rlutcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1166	'rlutcs4co2', 'TOA Out CS LW in 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1167	TYPE(ctrl_out), SAVE :: o_rsu4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1168	'rsu4co2', 'Upwelling SW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1169	TYPE(ctrl_out), SAVE :: o_rlu4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1170	'rlu4co2', 'Upwelling LW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1171	TYPE(ctrl_out), SAVE :: o_rsucs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1172	'rsucs4co2', 'Upwelling CS SW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1173	TYPE(ctrl_out), SAVE :: o_rlucs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1174	'rlucs4co2', 'Upwelling CS LW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1175	TYPE(ctrl_out), SAVE :: o_rsd4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1176	'rsd4co2', 'Downwelling SW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1177	TYPE(ctrl_out), SAVE :: o_rld4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1178	'rld4co2', 'Downwelling LW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1179	TYPE(ctrl_out), SAVE :: o_rsdcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1180	'rsdcs4co2', 'Downwelling CS SW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1181	TYPE(ctrl_out), SAVE :: o_rldcs4co2 = ctrl_out((/ 5, 10, 10, 10, 10, 10, 11, 11, 11 /), &
1182	'rldcs4co2', 'Downwelling CS LW 4xCO2 atmosphere', 'W/m2', (/ ('', i=1, 9) /))
1183	TYPE(ctrl_out), SAVE :: o_snowsrf = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1184	'snowsrf', 'Snow mass at surface', 'kg/m2', (/ ('', i=1, 9) /))
1185	TYPE(ctrl_out), SAVE :: o_qsnow = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1186	'qsnow', 'Water contained in snow', 'kg/m2', (/ ('', i=1, 9) /))
1187	TYPE(ctrl_out), SAVE :: o_snowhgt = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1188	'snowhgt', 'Snow height at surface', 'm', (/ ('', i=1, 9) /))
1189	TYPE(ctrl_out), SAVE :: o_toice = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1190	'to_ice', 'Snow passed to ice model', 'kg/m2', (/ ('', i=1, 9) /))
1191	TYPE(ctrl_out), SAVE :: o_sissnow = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1192	'sissnow', 'Snow in snow model', 'kg/m2', (/ ('', i=1, 9) /))
1193	TYPE(ctrl_out), SAVE :: o_runoff = ctrl_out((/ 1, 1, 10, 1, 10, 10, 11, 11, 11 /), &
1194	'runoff', 'Run-off rate land ice', 'kg/m2/s', (/ ('', i=1, 9) /))
1195	TYPE(ctrl_out), SAVE :: o_albslw3 = ctrl_out((/ 1, 1, 1, 1, 10, 10, 11, 11, 11 /), &
1196	'albslw3', 'Surface albedo LW3', '-', (/ ('', i=1, 9) /))
1197
1198	!!!!!!!!!!!!! Sorties niveaux standards de pression NMC
1199	TYPE(ctrl_out), SAVE :: o_tnondef = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1200	'tnondef', 'Undefined value of T', 'K', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1201	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1202	TYPE(ctrl_out), SAVE :: o_ta = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1203	'ta', 'Air temperature', 'K', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1204	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1205	TYPE(ctrl_out), SAVE :: o_zg = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1206	'zg', 'Geopotential height', 'm', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1207	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1208	TYPE(ctrl_out), SAVE :: o_hus = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1209	'hus', 'Specific humidity', '1', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1210	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1211	TYPE(ctrl_out), SAVE :: o_hur = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1212	'hur', 'Relative humidity', '%', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1213	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1214	TYPE(ctrl_out), SAVE :: o_ua = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1215	'ua', 'Eastward wind', 'm s-1', (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
1216	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1217	TYPE(ctrl_out), SAVE :: o_va = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1218	'va', 'Northward wind', 'm s-1', (/ ('', i=1, 9)/))
1219	TYPE(ctrl_out), SAVE :: o_wap = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1220	'wap', 'Lagrangian tendency of air pressure', 'Pa s-1', (/ "inst(X)", "inst(X)", "inst(X)", &
1221	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1222	TYPE(ctrl_out), SAVE :: o_psbg = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1223	'psbg', 'Pressure sfce below ground', '%', (/ "inst(X)", "inst(X)", "inst(X)", &
1224	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1225	TYPE(ctrl_out), SAVE :: o_tro3 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1226	'tro3', 'Ozone mole fraction', '1e-9', (/ "inst(X)", "inst(X)", "inst(X)", &
1227	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1228	TYPE(ctrl_out), SAVE :: o_tro3_daylight = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
1229	'tro3_daylight', 'Daylight ozone mole fraction', '1e-9', (/ "inst(X)", "inst(X)", "inst(X)", &
1230	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1231	TYPE(ctrl_out), SAVE :: o_uxv = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1232	'uv', 'uv', 'm2/s2', (/ "inst(X)", "inst(X)", "inst(X)", &
1233	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1234	TYPE(ctrl_out), SAVE :: o_vxq = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1235	'vq', 'vq', 'm/s * (kg/kg)', (/ "inst(X)", "inst(X)", "inst(X)", &
1236	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1237	TYPE(ctrl_out), SAVE :: o_vxT = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1238	'vT', 'vT', 'mK/s', (/ "inst(X)", "inst(X)", "inst(X)", &
1239	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1240	TYPE(ctrl_out), SAVE :: o_wxq = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1241	'wq', 'wq', '(Pa/s)*(kg/kg)', (/ "inst(X)", "inst(X)", "inst(X)", &
1242	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1243	TYPE(ctrl_out), SAVE :: o_vxphi = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1244	'vphi', 'vphi', 'm2/s', (/ "inst(X)", "inst(X)", "inst(X)", &
1245	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1246	TYPE(ctrl_out), SAVE :: o_wxT = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1247	'wT', 'wT', '"K*Pa/s', (/ "inst(X)", "inst(X)", "inst(X)", &
1248	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1249	TYPE(ctrl_out), SAVE :: o_uxu = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1250	'u2', 'u2', 'm2/s2', (/ "inst(X)", "inst(X)", "inst(X)", &
1251	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1252	TYPE(ctrl_out), SAVE :: o_vxv = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1253	'v2', 'v2', 'm2/s2', (/ "inst(X)", "inst(X)", "inst(X)", &
1254	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1255	TYPE(ctrl_out), SAVE :: o_TxT = ctrl_out((/ 11, 11, 11, 11, 11, 11, 6, 6, 6 /), &
1256	'T2', 'T2', 'K2', (/ "inst(X)", "inst(X)", "inst(X)", &
1257	"inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
1258
1259	END MODULE phys_output_ctrlout_mod

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

Download in other formats:

Original Format