Context Navigation

phys_output_ctrlout_mod.F90 @ 1934

Last change on this file since 1934 was 1924, checked in by idelkadi, 11 years ago

Rajout de nouveaux diagnostiques pour la comparaison entre les fractions nuageuses sorties derictement dans LMDZ et celles sorties via le simulateur ISCCP. Routines concernees : physiq.F90, phys_output_ctrlout_mod.F90, phys_local_var_mod.F90 et phys_output_write_mod.F90

Rajout de la fraction d'ensoleillement par jour aux arguments de l'interface avec Cosp (entree pour le simulateur ISCCP). Routine concernée : physiq.F90

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.0 KB

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

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

Download in other formats:

Original Format