Parameter list and GRIB definitions

HARMONIE system output

The HARMONIE system writes its primary output, in FA format, to the upper air history files ICMSHHARM+llll and the SURFEX history files ICMSHHARM+llll.sfx, where HARM is the four-character experiment identifier set in the configuration file config_exp.h, and llll is normally the current timestep in hours. The files are designed to be complete snapshots of respective model state described by the system for a particular time point. In addition more model output including post-processing/diagnostic fields can be written out during the forecast model integration, such as those model diagnostics or pressure level diagnostics, also in FA format, as PFHARMDOMAIN+llll. The FA files can be considered to be internal format files. All of them can be converted to GRIB files during the run for external usage. The name convention is as follows:

Forecast upper air history files: ICMSHHARM+llll -> fcYYYYMMDDHH+lll_grib (GRIB1) or fcYYYYMMDDHH+lll_grib2 (GRIB2)
Forecast Surfex history files: ICMSHHARM+llll.sfx -> fcYYYYMMDDHH+lll_grib_sfx (GRIB1 only)
Forecast Surfex selected output: ICMSSELE+llll.sfx -> fcYYYYMMDDHH+lll_grib_sfxs (GRIB1 only)
Postprocess files: PFHARMDOMAIN+llll.hfp -> fcYYYYMMDDHH+lllgrib_fp (GRIB1) or fcYYYYMMDDHH+lllgrib2_fp (GRIB2)
Analysis upper air history files: ICMSHANAL+0000 -> anYYYYMMDDHH+000grib (GRIB1) or anYYYYMMDDHH+000grib2 (GRIB2) (1)
Analysis SURFEX history files: ICMSHANAL+0000.sfx -> sa2019041600+000grib_sfx (only GRIB1 for the time being)

GRIB1 table 2 version in HARMONIE

To avoid conflicts with archived HIRLAM data HARMONIE uses version 253 of table 2. The table is based on the standard WMO version 3 of table 2 and postion 000-127 is kept the same as in the WMO. Note that accumulated and instantaneous versions of the same parameter differ only by the time range indicator. It is thus not sufficient to specify parameter, type and level when you refer to an accumulated parameter, but the time range indicator has to be included as well.

The translation of SURFEX files to GRIB1 is still incomplete and contains several WMO violations. This is not changed in the current release but will revised later. However, the upper air history file also includes the most common surface parameters and should be sufficient for most users.

The current table 2 version 253 definition files for gribapi can be found in `util/glgrib_api/definitions/`. These local definition files assume centre=233 (Dublin) and should be copied to your own GRIB-API installation. You are strongly recommended to set your own code for generating centre fore operational usage of the data.

GRIB2 in HARMONIE

The possibility to convert to GRIB2 has been introduced in release-43h2. So far the conversion is restricted to atmospheric history and fullpos files only. To get the output in GRIB2 set ARCHIVE_FORMAT=GRIB2 in ecf/config_exp.h. Please notice that if ARCHIVE_FORMAT=GRIB2 is selected, SURFEX files will be converted to GRIB1 anyway (for the time being). To convert from GRIB1 with GRIB2 using grib_filter we have to tell EcCodes how to translate the parameters. This is done by using the internal HARMONIE tables and setting

export ECCODES_DEFINITION_PATH=$SOME_PATH_TO_GL/gl/definitions:$SOME_PATH_TO_ECCODES/share/eccodes/definitions

Note that there are a few parameters that are not translated to GRIB2 to and those has to be excluded explicitly.

List of parameters

header abbreviations in the tables:

abbr.	description	see table
lvT	levelType	level types
iOP	indicatorOfParameter	indicator of parameter
d	discipline
pC	parameterCategory
pN	parameterNumber
lev	level
sT	stepType	time range indicator

3D model state variables on model levels (1-NLEV), levelType=hybrid

FA name	shortName	NCname	lvT	iOP	pC	pN	lev	sT	units	description
SNNNHUMI.SPECIFI	q	hus	hy	51	1	0	1	ins	kg/kg	Specific humidity
SNNNLIQUID_WATER	cwat_cond	clw	hy	76	1	83	1	ins	kg/kg	Specific cloud liquid water content
SNNNSOLID_WATER	ciwc_cond	cli	hy	58	1	84	1	ins	kg/kg	Specific cloud ice water content
SNNNSNOW	snow_cond	#	hy	184	1	86	1	ins	kg/kg	Specific snow water content
SNNNRAIN	rain_cond	#	hy	181	1	85	1	ins	kg/kg	Specific rain water content
SNNNGRAUPEL	grpl_cond	#	hy	201	1	32	1	ins	kg/kg	Specific graupel
SNNNTKE	tke	tke	hy	200	19	11	1	ins	J/kg	Turbulent Kinetic Energy
SNNNCLOUD_FRACTI	tcc	clt	hy	71	6	192	1	ins	0-1	Total cloud cover
SNNNPRESS.DEPART	pdep	#	hy	212	3	8	1	ins	Pa	Pressure departure
SNNNTEMPERATURE	t	ta	hy	11	0	0	1	ins	K	Temperature
SNNNVERTIC.DIVER	vdiv	#	hy	213	2	192	1	ins	s-1	Vertical Divergence
SNNNWIND.U.PHYS	u	ua	hy	33	2	2	1	ins	m/s	u-component of wind
SNNNWIND.V.PHYS	v	va	hy	34	2	3	1	ins	m/s	v-component of wind

2D Surface, prognostic/diagnostic near-surface and soil variables, levelType=heightAboveGround

FA name	shortName	NCname	lvT	iOP	pC	pN	lev	sT	units	description
SURFPRESSION	pres	ps	hag	1	3	0	0	ins	Pa	Surface pressure
SURFTEMPERATURE	t	ts_rad	hag	11	0	0	0	ins	K	Surface temperature
CLSTEMPERATURE	t	tas	hag	11	0	0	2	ins	K	Temperature at 2m
CLSMAXI.TEMPERAT	tmax	tasmax	hag	15	0	0	2	max	K	Maximum temperature (`FREQ_RESET_TEMP`)
CLSMINI.TEMPERAT	tmin	tasmin	hag	16	0	0	2	min	K	Minimum temperature (`FREQ_RESET_TEMP`)
CLSVENT.ZONAL	u	uas	hag	33	2	2	10	ins	m/s	u-component of wind at 10m, relative to model coordinates
CLSVENT.MERIDIEN	v	vas	hag	34	2	3	10	ins	m/s	v-component of wind at 10m, relative to model coordinates
CLSHUMI.SPECIFIQ	q	huss	hag	51	1	0	2	ins	kg/kg	Specific humidity at 2m
CLSHUMI.RELATIVE	r	hurs	hag	52	1	192	2	ins	0-1	Relative humidity at 2m
SURFRESERV.NEIGE	sdwe	snw	hag	65	1	60	0	ins	kg/m2	Snow depth water equivalent
CLPMHAUT.MOD.XFU	mld	zmla	hag	67	19	3	0	ins	m	Height (in meters) of the PBL out of the model
SURFNEBUL.TOTALE	tcc	clt_ins	hag	71	6	192	0	ins	0-1	Total cloud cover
SURFNEBUL.CONVEC	ccc	clc_ins	hag	72	6	193	0	ins	0-1	Convective cloud cover
SURFNEBUL.BASSE	lcc	cll_ins	hag	73	6	194	0	ins	0-1	Low cloud cover
SURFNEBUL.MOYENN	mcc	clm_ins	hag	74	6	195	0	ins	0-1	Medium cloud cover
SURFNEBUL.HAUTE	hcc	clh_ins	hag	75	6	196	0	ins	0-1	High cloud cover
SURFRAYT.SOLAIRE	swavr	#	hag	116	4	7	0	ins	W/m2	Instantaneous surface solar radiation (SW down global) Parameter identifier was 116, again is???
SURFRAYT.TERREST	lwavr	#	hag	115	5	4	0	ins	W/m2	Instantaneous longwave radiation flux
SURFCAPE.MOD.XFU	cape	cape	hag	160	7	6	0	ins	J/kg	Model output CAPE (not calculated by AROME physics)
SURFDIAGHAIL	xhail	#	hag	161	1	203	0	ins	0-1	AROME hail diagnostic, `LXXDIAGH = .TRUE.`
CLSU.RAF.MOD.XFU	ugst	ugs	hag	162	2	23	10	max	m/s	U-momentum of gusts from the model. `LXXGST = .TRUE.` in `NAMXFU`. gives gust between current and previous output time step (`FREQ_RESET_GUST`)
CLSV.RAF.MOD.XFU	vgst	vgs	hag	163	2	24	10	max	m/s	V-momentum of gusts from the model. `LXXGST = .TRUE.` in `NAMXFU`. gives gust between current and previous output time step (`FREQ_RESET_GUST`)
SURFINSPLUIE	rain	#	hag	181	1	65	0	ins	kg/m2	Instantaneous rain
SURFINSNEIGE	snow	#	hag	184	1	53	0	ins	kg/m2	Instantaneous snow
SURFINSGRAUPEL	grpl	#	hag	201	1	75	0	ins	kg/m2	Instantaneous graupel
CLSMINI.HUMI.REL	rmn2m	#	hag	241	1	1	2	min	0-1	Minimum relative moisture at 2m over 3h
CLSMAXI.HUMI.REL	rmx2m	#	hag	242	1	1	2	max	0-1	Maximum relative moisture at 2m over 3h
CLSRAFALES.POS	fg	wsgsmax	hag	228	2	22	10	max	m/s	Gust wind speed

2D Surface, accumulated near-surface and soil variables

Note that all these are coded with stepType=accum

FA name	shortName	NCname	lvT	iOP	d	pC	pN	lev	sT	units	description
S065RAYT SOL CL	cssw	#	hy	130	0	4	11	65	acc	J/m2	SW net clear sky rad
S065RAYT THER CL	cslw	#	hy	131	0	5	6	65	acc	J/m2	LW net clear sky rad
SURFACCGRAUPEL	grpl	prgrpl	hag	201	0	1	75	0	acc	kg/m2	Accumulated graupel
SURFACCNEIGE	snow	prsn	hag	184	0	1	53	0	acc	kg/m2	Accumulated snowfall
SURFACCPLUIE	rain	prrain	hag	181	0	1	65	0	acc	kg/m2	Accumulated rain
SURFDIR NORM IRR	dne	rids	hag	140	3	6	3	0	acc	J/m2	Direct normal exposure
SURFFLU.CHA.SENS	shf	hfss	hag	122	0	0	11	0	acc	J/m2	Sensible heat flux
SURFFLU.LAT.MEVA	lhe	hfls_eva	hag	132	0	1	193	0	acc	J/m2	Latent heat flux through evaporation
SURFFLU.LAT.MSUB	lhsub	hfls_sbl	hag	244	0	1	202	0	acc	J/kg	Latent Heat Sublimation
SURFFLU.MEVAP.EA	wevap	evspsbl	hag	245	0	1	6	0	acc	kg/m2	Water evaporation
SURFFLU.MSUBL.NE	snsub	sbl_snow	hag	246	0	1	62	0	acc	kg/m2	Snow sublimation
SURFFLU.RAY.SOLA	nswrs	rsns	hag	111	0	4	9	0	acc	J/m2	Net shortwave radiation flux (surface)
SURFFLU.RAY.THER	nlwrs	rlns	hag	112	0	5	5	0	acc	J/m2	Net longwave radiation flux (surface)
SURFRAYT DIR SUR	swavr	rsdsdir	hag	116	0	4	7	0	acc	J/m2	Shortwave radiation flux
SURFRAYT SOLA DE	grad	rsds	hag	117	0	4	3	0	acc	J/m2	Global radiation flux
SURFRAYT THER DE	lwavr	rlds	hag	115	0	5	4	0	acc	J/m2	Longwave radiation flux
SURFTENS.TURB.ME	vflx	tauv	hag	125	0	2	199	0	acc	N/m2	Momentum flux, v-component
SURFTENS.TURB.ZO	uflx	tauu	hag	124	0	2	198	0	acc	N/m2	Momentum flux, u-component

2D TOA, diagnostic and accumulated variables, levelType=nominalTop

FA name	shortName	NCname	lvT	iOP	d	pC	pN	sT	units	description
SOMMFLU.RAY.SOLA	nswrt	rsnt	nt	113	0	4	9	acc	J/m2	Net shortwave radiation flux(atmosph.top)
SOMMFLU.RAY.THER	nlwrt	rlnt	nt	114	0	5	5	acc	J/m2	Net longwave radiation flux(atmosph.top)
SOMMRAYT.SOLAIRE	nswrt	#	nt	113	0	4	9	ins	W/m2	Net shortwave radiation flux(atmosph.top)
SOMMRAYT.TERREST	nlwrt	#	nt	114	0	5	5	ins	W/m2	Net longwave radiation flux(atmosph.top)
SOMMRAYT SOL CL	cssw	rsntcs	nt	130	0	4	11	acc	J/m2	TOA Net shortwave clear sky radiation(atmosph.top)
SOMMRAYT THER CL	cslw	rlntcs	nt	131	0	5	6	acc	J/m2	TOA Net longwave clear sky radiation(atmosph.top)
TOPRAYT DIR SOM	swavr	rsdt	nt	116	0	4	7	acc	J/m2	TOA Accumulated SW down radiation Parameter identifier was 117
SOMMTBOZCLEAR	btozcs	#	nt	170	-1	-1	-1	-	K	Brightness temperature OZ clear
SOMMTBOZCLOUD	btozcl	#	nt	171	-1	-1	-1	-	K	Brightness temperature OZ cloud
SOMMTBIRCLEAR	btircs	#	nt	172	-1	-1	-1	-	K	Brightness temperature IR clear
SOMMTBIRCLOUD	btircl	#	nt	173	-1	-1	-1	-	K	Brightness temperature IR cloud
SOMMTBWVCLEAR	btwvcs	#	nt	174	-1	-1	-1	-	K	Brightness temperature WV clear
SOMMTBWVCLOUD	btwvcl	#	nt	175	-1	-1	-1	-	K	Brightness temperature WV cloud

2D Surface, Postprocessed variables (fullpos)

FA name	shortName	NCname	lvT	iOP	pC	pN	sT	units	description
SURFCAPE.POS.F00	cape	cape	hag	160	7	6	ins	J/kg	Convective available potential energy (CAPE)
SURFCIEN.POS.F00	cin	cin	hag	165	7	7	ins	J/kg	Convective inhibition (CIN)
SURFLIFTCONDLEV	lcl	#	ac	167	3	6	ins	m	Lifting condensation level (LCL)
SURFFREECONVLEV	lfc	#	lfc	168	3	6	ins	m	Level of free convection (LFC)
SURFEQUILIBRLEV	lnb	#	lnb	169	3	6	ins	m	Level of neutral buoyancy (LNB)

2D Surface, constant near-surface and soil variables

FA name	shortName	NCname	lvT	iOP	d	pC	pN	lev	sT	units	description
SPECSURFGEOPOTEN	z	phis_s	hag	6	0	3	4	0	ins	m2/s2	Geopotential relative to mean sea level. "... contains a GRID POINT orography which is the interpolation of the departure orography"
SURFIND.TERREMER	lsm	lsm	hag	81	2	0	0	0	ins	0-1	Land-sea mask
SURFAEROS.SEA	aers	#	hag	251	0	13	192	0	ins	kg/kg	Surface aerosol sea (Marine aerosols, locally defined GRIB)
SURFAEROS.LAND	aerl	#	hag	252	0	13	193	0	ins	kg/kg	Surface aerosol land (Continental aerosols, locally defined GRIB)
SURFAEROS.SOOT	aerc	#	hag	253	0	13	194	0	ins	kg/kg	Surface carbon aerosol (Carbone aerosols, locally defined GRIB)
SURFAEROS.DESERT	aerd	#	hag	254	0	13	195	0	ins	kg/kg	Surface aerosol desert (Desert aerosols, locally defined GRIB)
SURFAEROS.VOLCAN	#	#	hag	197	-1	-1	-1	-1			Surface aerosol volcan (Stratospheric ash, to be locally defined GRIB)
SURFAEROS.SULFAT	#	#	hag	198	-1	-1	-1	-1			Surface aerosol sulfate (Stratospheric sulfate, to be locally defined GRIB)
SURFA.OF.OZONE	ao	#	hag	248	0	14	192	0	ins	kg/kg	A Ozone, First ozone profile (A), locally defined GRIB
SURFB.OF.OZONE	bo	#	hag	249	0	14	193	0	ins	kg/kg	B Ozone, Second ozone profile (B), locally defined GRIB
SURFC.OF.OZONE	co	#	hag	250	0	14	194	0	ins	kg/kg	C Ozone, Third ozone profile (C), locally defined GRIB
PROFTEMPERATURE	slt	#	dbl	85	2	3	18	0	ins	K	Soil Temperature
PROFRESERV.EAU	sm	#	dbl	86	2	3	20	0	ins	kg/m2	Deep Soil Wetness
PROFPROP.RMAX.EA	swv	#	dbl	238	2	3	25	0	ins	kg/m2	Climate relaxed deep soil wetness
PROFRESERV.GLACE	wsoice	#	dbl	193	2	3	22	0	ins	kg/m2	Deep soil ice

2D variables on special surfaces

FA name	shortName	NCname	lvT	iOP	pC	pN	lev	sT	units	description
KT273ISOT_ALTIT	h	#	isot	8	3	6	27315	ins	m	Altitude of 0-degree isotherm
KT263ISOT_ALTIT	h	#	isot	8	3	6	26315	ins	m	Altitude of -10-degree isotherm
SURFISOTPW0.MALT	h	#	isot0wb	8	3	6	0	ins	m	Altitude of iso-tprimw=0
SURFTOT.WAT.VAPO	wvint	prw	ea	54	1	64	0	ins	kg/m2	Total column integral water vapour
WFPOWERINS	wfpower_ins	wfpower_ins	ea	211	2	39	0	ins	MW	Wind power production, instantaneous (`LWINDFARM=.TRUE.` in `NAMPHY`)
WFPOWERACC	wfpower_acc	wfpower_acc	ea	211	2	39	0	acc	MJ	Wind power production, accumulated (`LWINDFARM=.TRUE.` in `NAMPHY`)

Postprocessed variables on different surface types

Through the postprocessing sofware fullpos HARMONIE offers a number of variables postprocessed on different surface types. For the current choice of variables, surfaces and levels please see scr/Select_postp.pl.

State variables and diagnostics on pressure levels, leveltype=isobaricInhPa

FA name	shortName	NCname	lvT	iOP	pC	pN	lev	sT	units	description
PNNNNNWIND.U.PHY	u	ua	pl	33	2	2	NNNNN	ins	m/s	u-component of wind
PNNNNNWIND.V.PHY	v	va	pl	34	2	3	NNNNN	ins	m/s	v-component of wind
PNNNNNTEMPERATUR	t	ta	pl	11	0	0	NNNNN	ins	K	Temperature
PNNNNNHUMI.SPECI	q	hus	pl	51	1	0	NNNNN	ins	kg/kg	Specific humidity
PNNNNNLIQUID_WAT	cwat_cond	clw	pl	76	1	83	NNNNN	ins	kg/kg	Specific cloud liquid water content
PNNNNNSOLID_WATE	ciwc_cond	cli	pl	58	1	84	NNNNN	ins	kg/kg	Specific cloud ice water content
PNNNNNCLOUD_FRAC	tcc	#	pl	71	6	192	NNNNN	ins	0-1	Total cloud cover
PNNNNNSNOW	snow_cond	#	pl	184	1	86	NNNNN	ins	kg/kg	Specific snow water content
PNNNNNRAIN	rain_cond	#	pl	181	1	85	NNNNN	ins	kg/kg	Specific rain water content
PNNNNNGRAUPEL	grpl_cond	#	pl	201	1	32	NNNNN	ins	kg/kg	Specific graupel
PNNNNNGEOPOTENTI	z	phi	pl	6	3	4	NNNNN	ins	m2/s2	Geopotential
PNNNNNHUMI_RELAT	r	hur	pl	52	1	192	NNNNN	ins	0-1	Relative humidity
PNNNNNTHETA_PRIM	papt	thetaE	pl	14	0	3	NNNNN	ins	K	Pseudo-adiabatic potential temperature
PNNNNNTHETA_VIRT	vptmp	#	pl	176	0	15	NNNNN	ins	K	Virtual potential temperature
PNNNNNVERT.VELOC	w	wa	pl	40	2	9	NNNNN	ins	m/s	Geometrical vertical velocity
PNNNNNPOT_VORTIC	pv	pv	pl	4	2	14	NNNNN	ins	K m2/kg/s	Potential vorticity
PNNNNNABS_VORTIC	absv	#	pl	41	2	10	NNNNN	ins	s-1	Absolute vorticity
PNNNNNDIVERGENCE	d	#	pl	44	2	13	NNNNN	ins	s-1	Relative divergence

NNNNN is in Pascals.
From FullPos documentation: "Warning: fields on pressure levels bigger or equal to 1000 hPa are written out with truncated names; for example, temperature at 1000 hPa is P00000TEMPERATURE while P00500TEMPERATURE could be as well the temperature at 5 hPa or the temperature at 1005 hPa!"

State variables and diagnostics on height levels, levelType=heightAboveGround

FA name	shortName	NCname	lvT	iOP	pC	pN	lev	sT	units	description
HNNNNNWIND.U.PHY	u	ua	hag	33	2	2	NNNNN	ins	m/s	u-component of wind
HNNNNNWIND.V.PHY	v	va	hag	34	2	3	NNNNN	ins	m/s	v-component of wind
HNNNNNTEMPERATUR	t	ta	hag	11	0	0	NNNNN	ins	K	Temperature
HNNNNNLIQUID_WAT	cwat_cond	clw	hag	76	1	83	NNNNN	ins	kg/kg	Specific cloud liquid water content
HNNNNNSOLID_WATE	ciwc_cond	cli	hag	58	1	84	NNNNN	ins	kg/kg	Specific cloud ice water content
HNNNNNCLOUD_FRAC	tcc	clt	hag	71	6	192	NNNNN	ins	0-1	Total cloud cover
HNNNNNSNOW	snow_cond	#	hag	184	1	86	NNNNN	ins	kg/kg	Specific snow water content
HNNNNNRAIN	rain_cond	#	hag	181	1	85	NNNNN	ins	kg/kg	Specific rain water content
HNNNNNGRAUPEL	grpl_cond	#	hag	201	1	32	NNNNN	ins	kg/kg	Specific graupel
HNNNNNHUMI_RELAT	r	hur	hag	52	1	192	NNNNN	ins	0-1	Relative humidity
HNNNNNPRESSURE	pres	p	hag	1	3	0	NNNNN	ins	Pa	Pressure

NNNNN is in meters.

State variables and diagnostics on PV levels, GRIB1 level type 117, levelType=potentialVorticity

FA name	shortName	NCname	lvT	iOP	pC	pN	lev	sT	units	description
VNNNGEOPOTENTIEL	z	#	pv	6	3	4	NNN	ins	m2/s2	Geopotential
VNNNTEMPERATURE	t	#	pv	11	0	0	NNN	ins	K	Temperature
VNNNPRESSURE	pres	#	pv	1	3	0	NNN	ins	Pa	Pressure
VNNNHUMI_RELATIV	r	#	pv	52	1	192	NNN	ins	0-1	Relative humidity
VNNNHUMI.SPECIFI	q	#	pv	51	1	0	NNN	ins	kg/kg	Specific humidity
VNNNWIND.U.PHYS	u	#	pv	33	2	2	NNN	ins	m/s	u-component of wind
VNNNWIND.V.PHYS	v	#	pv	34	2	3	NNN	ins	m/s	v-component of wind
VNNNVITESSE_VERT	omega	#	pv	39	2	8	NNN	ins	Pa/s	Pressure vertical velocity (`DYNAMICS=h`)
VNNNVERT.VELOCIT	w	#	pv	40	2	9	NNN	ins	m/s	Geometrical vertical velocity (`DYNAMICS=nh`)
VNNNTEMPE_POTENT	pt	#	pv	13	0	2	NNN	ins	K	Potential temperature
VNNNABS_VORTICIT	absv	#	pv	41	2	10	NNN	ins	s-1	Absolute vorticity
VNNNDIVERGENCE	d	#	pv	44	2	13	NNN	ins	s-1	Relative divergence
VNNNTHETAPRIMW	papt	#	pv	14	0	3	NNN	ins	K	Pseudo-adiabatic potential temperature

"pv" stream is not available by default
NNN is in deci-PVU (1PVU = 1x10-6 K m2 kg-1 s-1) in FA files. PV levels must be in SI units in namelists
GRIB1 levels are in milli-PVU. Currently gl does not convert devi-PVU (FA) to milli-PVU (GRIB1)

State variables and diagnostics on Theta levels, GRIB1 level type 113, levelType=theta

FA name	shortName	NCname	lvT	iOP	pC	pN	lev	sT	units	description
TNNNGEOPOTENTIEL	z	#	th	6	3	4	NNN	ins	m2/s2	Geopotential
TNNNTEMPERATURE	t	#	th	11	0	0	NNN	ins	K	Temperature
TNNNPRESSURE	pres	#	th	1	3	0	NNN	ins	Pa	Pressure
TNNNHUMI_RELATIV	r	#	th	52	1	192	NNN	ins	0-1	Relative humidity
TNNNHUMI.SPECIFI	q	#	th	51	1	0	NNN	ins	kg/kg	Specific humidity
TNNNWIND.U.PHYS	u	#	th	33	2	2	NNN	ins	m/s	u-component of wind
TNNNWIND.V.PHYS	v	#	th	34	2	3	NNN	ins	m/s	v-component of wind
TNNNVITESSE_VERT	omega	#	th	39	2	8	NNN	ins	Pa/s	Pressure vertical velocity (`DYNAMICS=h`)
TNNNVERT.VELOCIT	w	#	th	40	2	9	NNN	ins	m/s	Geometrical vertical velocity (`DYNAMICS=nh`)
TNNNABS_VORTICIT	absv	#	th	41	2	10	NNN	ins	s-1	Absolute vorticity
TNNNPOT_VORTICIT	pv	#	th	4	2	14	NNN	ins	K m2/kg/s	Potential vorticity
TNNNDIVERGENCE	d	#	th	44	2	13	NNN	ins	s-1	Relative divergence

"th" stream is not available by default
NNN is in Kelvin.

FA fields without any default GRIB1 translation

Some very special fields are left without any default translation. Please see in the gl documentation on how to add you own translation.

FA name	Unit	Comment
CUF1PRESSURE		Coupling error field.
THETAPWP_FLUX	K m-4 s-1	Instantaneous thetaprimwprim surface flux
CLPMOCON.MOD.XFU	kg kg-1 s-1	MOCON model output
ATMONEBUL.TOTALE		Accumulated Total cloud cover.
ATMONEBUL.CONVEC		Accumulated Convective cloud cover.
ATMONEBUL.BASSE		Accumulated Low cloud cover.
ATMONEBUL.MOYENN		Accumulated Medium cloud cover.
ATMONEBUL.HAUTE		Accumulated High cloud cover.
SURFCFU.Q.TURBUL		Accumulated contribution of Turbulence to Q.
SURFCFU.CT.TURBUL		Accumulated contribution of Turbulence to CpT
SUNSHI. DURATION		Sunshine duration.
SURFFL.U TURBUL		Contribution of Turbulence to U.
SURFFL.V TURBUL		Contribution of Turbulence to V.
SURFFL.Q TURBUL		Contribution of Turbulence to Q.
SURFFL.CT TURBUL		Contribution of Turbulence to CpT
SNNNSRC		Second order flux.

Variables postprocessed by gl

The following fields are can be generated by gl from a history file and are thus not necessarily available as FA fields in Harmonie's FA output. When calculating these post-processed fields, make sure the required fields to derive them are in the input files! For details, check util/gl/grb/postprocess.f90 & the routines called therein.

Single level fields

FA name	shortName	NCname	lvT	iOP	d	pC	pN	lev	sT	units	description
MSLPRESSURE	pres	psl	has	1	0	3	0	0	ins	Pa	MSLP. gl calculates MSLP independent of AROME/!FullPos
#	td	td	hag	17	0	0	6	2	ins	K	Dew point temperature
#	vis	#	hag	20	0	19	0	0	ins	m	Visibility
#	wdir	#	ttt	31	0	2	0	lll	ins	Deg. true	Wind direction. gl calculates based on u[33,ttt,lll] and v[34,ttt,lll] wind components
#	ws	#	ttt	32	0	2	1	lll	ins	m/s	Wind speed. gl calculates based on u[33,ttt,lll] and v[34,ttt,lll] wind components
TOT.WATER.PRECIP	tp	pr	hag	61	0	1	8	0	acc	kg/m2	Total precipitation, gl calculates TP![61,105,0]=rain![181,105,0]+snow![184,105,0]+graupel![201,105,0]+hail![204,105,0]
TOT.SOLID.PRECIP	tpsolid	prsolid	hag	185	0	1	200	0	acc	kg/m2	Total solid precipitation, gl calculates ![185,105,0]=snow![184,105,0]+graupel![201,105,0]+hail![204,105,0]
#	mld	zmla	hag	67	0	19	3	0	ins	m	Mixed layer depth/boundary layer height
#	tcc	#	hag	71	0	6	192	2	ins	0-1	Fog, cloud fraction of lowest model level
#	icei	#	hag	135				0	ins	-	Icing index
#	atmiceg	#	hy?	?	0	1	205		ins	m/s	Icing index, Atmospheric ice growth rate
#	icei2	#	hag/?	134	0	1	194	0	ins	-	Icing index version 2
#	psct	#	hag/ct?	136	0	4	0	0	ins	K	Pseudo satellite image, cloud top temperature (infrared)
#	pstb	#	hag	137	0	4	198	0	ins	K	Pseudo satellite image, water vapour brightness temperature
#	pstbc	#	hag	138	0	4	199	0	ins	K	Pseudo satellite image, water vapour br. temp. + correction for clouds
#	pscw	#	hag	139	0	4	200	0	ins	-	Pseudo satellite image, cloud water reflectivity (visible)
#	prtp	#	hag	144	0	1	19	0	ins	code	Precipitation type, 0:drizzle, 1:rain, 2:sleet, 3:snow, 4:freezing drizzle, 5:freezing rain, 6:graupel, 7:hail
#	fg	#	ttt	228	0	2	22	lll	max	m/s	Gust wind speed, calculated from `ugst` & `vgst` on corresponding level & levelType
#	hti	#	hag	148	0	17	193	0	ins	-	Helicopter Triggered lightning Index
#	transmit	#	hag	149	0	6	199	0	ins	-	Transmittance
#	cat	#	hag	145	0	19	22	0	ins	-\|%	CAT (clear air turbulence) index
#	bvf	#	hag	159	0	19	202	0	ins	s-1	Brunt Vaisala frequency

Integrated quantities

FA name	shortName	NCname	lvT	iOP	pC	pN	lev	sT	units	description
TOT.NEB.ICE	ciwc_vint	clivi	ea	58	1	70	0	ins	kg/m2	Vertical integral of cloud ice
TOT.NEB.WAT	cwat_vint	clqvi	ea	76	1	69	0	ins	kg/m2	Vertical integral of cloud liquid water
#	msca	#	ea	133	6	197	0	ins	0-1	Mask of significant cloud amount
#	cape	#	hag	160	7	6	0	ins	J/kg	Convective Available Potential Energy, comes in two flavours, `cape_version=1\|2`, where the second is compatible with the ECMWF version
#	cin	#	hag	165	7	7	0	ins	J/kg	Convective inhibition, , comes in two flavours, `cape_version=1\|2`, where the second is compatible with the ECMWF version
#	rain_vint	clrvi	ea	181	1	45	0	ins	kg/m2	Vertical integral of rain
#	snow_vint	clsvi	ea	184	1	46	0	ins	kg/m2	Vertical integral of snow
#	grpl_vint	clgvi	ea	201	1	74	0	ins	kg/m2	Vertical integral of garupel
#	cb	#	ea	186	6	11	0	ins	m	Cloud base
#	ct	#	ea	187	6	12	0	ins	m	Cloud top
#	cb38	#	hag?	189	6	198	3	ins	m	Cloud base >3/8
#	lgt	#	ea	209	17	192	0	ins	flash/h	Lightning intensity
#	lmxws	#	ea/mw?	142	3	6	0	ins	m	Level of max wind speed
#	maxucol	#	ea	164	2	2	0	ins	m/s	Max u-component in column
#	maxvcol	#	ea	177	2	3	0	ins	m/s	Max v-component in column
#	lmxice	#	ea	143	1	199	0	ins	m	Icing index, Level of max icing
#	mxicegr	#	ea	141	1	204	0	ins	1	Icing index, Max icing growth index (0-4)
#	blice	#	ea	146	19	6	0	ins	m	Icing index, Bottom level of icing
#	tlice	#	ea	147	19	5	0	ins	m	Icing index, Top level of icing
#	cat_maxlev	#	ea	150	19	201	0	ins	m	Level of max CAT index
#	cat_max	#	ea	151	19	197	0	ins	-	Max CAT index
#	cat_b	#	ea	152	19	198	0	ins	m	Bottom level of CAT
#	cat_t	#	ea	153	19	199	0	ins	m	Top level of CAT

GRIB encoding information

Time units, WMO code table 4

The following time units are used to encode GRIB edition 1 data

Code	Unit
0	Minute
1	Hour
13	15 minutes
14	30 minutes

Time range indicator, WMO code TABLE 5

Code	abbr	Definition
0	ins	Forecast product valid for reference time + P1 (P1 > 0), or Uninitialized analysis product for reference time (P1 = 0)
2	min/max	Product with a valid time ranging between reference time + P1 and reference time + P2. Used for min/max values
3	avg	Average (reference time + P1 to reference time + P2)
4	acc	Accumulation (reference time + P1 to reference time + P2) product considered valid at reference time + P2

Note that fields available as both instanteous and accumulated values like e.g. rain has the same parameter values and can only be distinguished by the time range indicator.

Level types, WMO Code table 3

level type	name	abbr	WMO/HIRLAM type definition	Units	notes
001	surface	sfc	Ground or water surface		WMO
002	cloudBase	cb	Cloud base level		WMO
003	cloudTop	ct	Level of cloud tops		WMO
004	isothermZero	isot0	Level of 0°C isotherm		WMO
005	adiabaticCondensation	ac	Level of Adiabatic Condensation Lifted from the Surface		WMO
006	maxWind	mw	Maximum wind level		WMO
007	tropopause	tp	Tropopause		WMO
008	nominalTop	nt	Top-of-atmosphere		WMO
020	isothermal	isot	Isothermal level	Temperature in 1/100 K	WMO
100	isobaricInhPa	pl	Isobaric level	hPa	WMO
102	meanSea	ms	At mean sea level
103	heightAboveSea	has	Specified altitude above mean sea level	Altitude in m	WMO
105	heightAboveGround	hag	Specified height above ground	Altitude in m	WMO
107	sigma		Sigma level	Sigma value in 1/10000	WMO
109	hybrid	hy	Hybrid level		WMO
112	depthBelowLandLayer	dbl
113	theta	th	Isentropic (theta) level	Potential temperature in K	WMO
117	potentialVorticity	pv	Potential vorticity surface	10-9 K m2 kg-1 s-1	WMO
192	isothermZeroWetBulb	isot0wb
200	entireAtmosphere	ea	Entire atmosphere (considered as a single layer)		WMO, vertically integrated
	levelFreeConvection	lfc			as heightAboveGround in GRIB1
	levelNeutralBuoyancy	lnb			as heightAboveGround in GRIB1

Harmonie GRIB1 code table 2 version 253 - Indicator of parameter

Below the indicator of parameter code table for the Harmonie model. It is based on the WMO code table 2 version 3 with local parameters added. Parameter indicators 128-254 are reserved for originating center use. Parameter indicators 000-127 should not be altered. In HARMONIE, radiation fluxes are assumed positive downwards (against the recommendation by WMO).

Par	Description	SI Units
000	Reserved	n/a
001	Pressure	Pa
002	Pressure reduced to MSL	Pa
003	Pressure tendency	Pa s-1
004	Potential vorticity	K m2 kg-1 s-1
005	ICAO Standard Atmosphere reference height	m
006	Geopotential	m2 s-2
007	Geopotential height	gpm
008	Geometrical height	m
009	Standard deviation of height	m
010	Total ozone	Dobson
011	Temperature	K
012	Virtual temperature	K
013	Potential temperature	K
014	Pseudo-adiabatic potential temperature	K
015	Maximum temperature	K
016	Minimum temperature	K
017	Dew-point temperature	K
018	Dew-point depression (or deficit)	K
019	Lapse rate	K m-1
020	Visibility	m
021	Radar spectra (1)	-
022	Radar spectra (2)	-
023	Radar spectra (3)	-
024	Parcel lifted index (to 500 hPa)	K
025	Temperature anomaly	K
026	Pressure anomaly	Pa
027	Geopotential height anomaly	gpm
028	Wave spectra (1)	-
029	Wave spectra (2)	-
030	Wave spectra (3)	-
031	Wind direction	Degree true
032	Wind speed	m s-1
033	u-component of wind	m s-1
034	v-component of wind	m s-1
035	Stream function	m2 s-1
036	Velocity potential	m2 s-1
037	Montgomery stream function	m2 s-1
038	Sigma coordinate vertical velocity	s-1
039	Vertical velocity	Pa s-1
040	Vertical velocity	m s-1
041	Absolute vorticity	s-1
042	Absolute divergence	s-1
043	Relative vorticity	s-1
044	Relative divergence	s-1
045	Vertical u-component shear	s-1
046	Vertical v-component shear	s-1
047	Direction of current	Degree true
048	Speed of current	m s-1
049	u-component of current	m s-1
050	v-component of current	m s-1
051	Specific humidity	kg kg-1
052	Relative humidity	%
053	Humidity mixing ratio	kg kg-1
054	Precipitable water	kg m-2
055	Vapor pressure	Pa
056	Saturation deficit	Pa
057	Evaporation	kg m-2
058	Cloud ice	kg m-2
059	Precipitation rate	kg m-2 s-1
060	Thunderstorm probability	%
061	Total precipitation	kg m-2
062	Large scale precipitation	kg m-2
063	Convective precipitation	kg m-2
064	Snowfall rate water equivalent	kg m-2 s-1
065	Water equivalent of accumulated snow depth	kg m-2
066	Snow depth	m
067	Mixed layer depth	m
068	Transient thermocline depth	m
069	Main thermocline depth	m
070	Main thermocline anomaly	m
071	Total cloud cover	%
072	Convective cloud cover	%
073	Low cloud cover	%
074	Medium cloud cover	%
075	High cloud cover	%
076	Cloud water	kg m-2
077	Best lifted index (to 500 hPa)	K
078	Convective snow	kg m-2
079	Large scale snow	kg m-2
080	Water temperature	K
081	Land cover (1 = land, 0 = sea)	Proportion
082	Deviation of sea level from mean	m
083	Surface roughness	m
084	Albedo	%
085	Soil temperature	K
086	Soil moisture content	kg m-2
087	Vegetation	%
088	Salinity	kg kg-1
089	Density	kg m-3
090	Water run-off	kg m-2
091	Ice cover (1 = ice, 0 = no ice)	Proportion
092	Ice thickness	m
093	Direction of ice drift	Degree true
094	Speed of ice drift	m s-1
095	u-component of ice drift	m s-1
096	v-component of ice drift	m s-1
097	Ice growth rate	m s-1
098	Ice divergence	s-1
099	Snow melt	kg m-2
100	Significant height of combined wind waves and swell	m
101	Direction of wind waves	Degree true
102	Significant height of wind waves	m
103	Mean period of wind waves	s
104	Direction of swell waves	Degree true
105	Significant height of swell waves	m
106	Mean period of swell waves	s
107	Primary wave direction	Degree true
108	Primary wave mean period	s
109	Secondary wave direction	Degree true
110	Secondary wave mean period	s
111	Net short-wave radiation flux (surface)	W m-2
112	Net long-wave radiation flux (surface)	W m-2
113	Net short-wave radiation flux (top of atmosphere)	W m-2
114	Net long-wave radiation flux (top of atmosphere)	W m-2
115	Long-wave radiation flux	W m-2
116	Short-wave radiation flux	W m-2
117	Global radiation flux	W m-2
118	Brightness temperature	K
119	Radiance (with respect to wave number)	W m-1 sr-1
120	Radiance (with respect to wave length)	W m-3 sr-1
121	Latent heat flux	W m-2
122	Sensible heat flux	W m-2
123	Boundary layer dissipation	W m-2
124	Momentum flux, u-component	N m-2
125	Momentum flux, v-component	N m-2
126	Wind mixing energy	J
127	Image data	-
128	Analysed RMS of PHI (CANARI)	m2 s-2
129	Forecasted RMS of PHI (CANARI)	m2 s-2
130	SW net clear sky rad	W m-2
131	LW net clear sky rad	W m-2
132	Latent heat flux through evaporation	W m-2
133	Mask of significant cloud amount	0-1
134	Icing index version 2	-
135	Icing index	Code table
136	Pseudo satellite image, cloud top temperature (infrared)	K
137	Pseudo satellite image, water vapour brightness temperature	K
138	Pseudo satellite image, water vapour br. temp. + correction for clouds	K
139	Pseudo satellite image, cloud water reflectivity (visible)	?
140	Direct normal irradiance	J m-2
141	Max icing growth index	-
142	Level of max wind speed	m
143	Level of max icing	m
144	Precipition Type	Code table
145	CAT index	- / %
146	Bottom level of icing	m
147	Top level of icing	m
148	Helicopter Triggered ligthning Index	-
149	Transmittance	-
150	Level of max CAT index	m
151	Max CAT index	-
152	Bottom level of CAT	m
153	Top level of CAT	m
154	Max Wind speed	m s-1
155	Available	#
156	Available	#
157	Available	#
158	Surface downward moon radiation	W m-2
159	ABrunt Vaisala frequency	s-1
160	CAPE	J kg-1
161	AROME hail diagnostic	%
162	U-momentum of gusts out of the model	m s-1
163	V-momentum of gusts out of the model	m s-1
164	Max u-component in column	m s-1
165	Convective inhibition (CIN)	J kg-1
166	MOCON out of the model	kg/kg s-1
167	Lifting condensation level (LCL)	m
168	Level of free convection (LFC)	m
169	Level of neutral boyancy (LNB)	m
170	Brightness temperature OZ clear	K
171	Brightness temperature OZ cloud	K
172	Brightness temperature IR clear	K
173	Brightness temperature IR cloud	K
174	Brightness temperature WV clear	K
175	Brightness temperature WV cloud	K
176	Virtual potential temperature	K
177	Max v-component in column	m s-1
178	Available	#
179	Available	#
180	Available	#
181	Rain	kg m-2
182	Stratiform Rain	kg m-2
183	Convective Rain	kg m-2
184	Snow	kg m-2
185	Total solid precipitation	kg m-2
186	Cloud base	m
187	Cloud top	m
188	Fraction of urban land	Proportion
189	Cloud base >3/8	m
190	Snow Albedo	Proportion
191	Snow density	kg/m3
192	Water on canopy	kg/m2
193	Soil ice	kg/m2
194	Available	#
195	Gravity wave stress U-comp	N/m2
196	Gravity wave stress V-comp	N/m2
197	Available	#
198	Available	#
199	Vegetation type	-
200	TKE	m2 s-2
201	Graupel	kg m-2
202	Stratiform Graupel	kg m-2
203	Convective Graupel	kg m-2
204	Hail	kg m-2
205	Stratiform Hail	kg m-2
206	Convective Hail	kg m-2
207	Available	#
208	Available	#
209	Lightning	flash h-1
210	Simulated reflectivity	dBz
211	Wind power production	MW or MJ
212	Pressure departure	Pa
213	Vertical divergence	s-1
214	UD_OMEGA	ms-1?
215	DD_OMEGA	ms-1?
216	UDMESHFRAC	-
217	DDMESHFRAC	-
218	PSHICONVCL	-
219	Surface albedo for non snow covered areas	Proportion
220	Standard deviation of orography * g	m2 s-2
221	Anisotropy coeff of topography	-
222	Direction of main axis of topography	rad
223	Roughness length of bare surface * g	m2 s-2
224	Roughness length for vegetation * g	m2 s-2
225	Fraction of clay within soil	Proportion
226	Fraction of sand within soil	Proportion
227	Maximum proportion of vegetation	Proportion
228	Gust wind speed	m s-1
229	Albedo of bare ground	Proportion
230	Albedo of vegetation	Proportion
231	Stomatal minimum resistance	s/m
232	Leaf area index	m2/m2
233	Thetaprimwprim surface flux	Km/s
234	Dominant vegetation index	-
235	Surface emissivity	-
236	Maximum soil depth	m
237	Soil depth	m
238	Soil wetness	kg/m2
239	Thermal roughness length * g	m2 s-2
240	Resistance to evapotransiration	s/m
241	Minimum relative moisture at 2 meters	%
242	Maximum relative moisture at 2 meters	%
243	Duration of total precipitations	s
244	Latent Heat Sublimation	W/m2
245	Water evaporation	kg/m2
246	Snow sublimation	kg/m2
247	Snow history	???
248	A OZONE	kg kg-1
249	B OZONE	kg kg-1
250	C OZONE	kg kg-1
251	Surface aerosol sea	kg kg-1
252	Surface aerosol land	kg kg-1
253	Surface aerosol soot	kg kg-1
254	Surface aerosol desert	kg kg-1
255	Missing value	n/a

SURFEX output Harmonie GRIB1 code table 2 version 001

Levels are used in the conversion of SURFEX output to GRIB to indicate tile/patch/type/level:

level	description
300	Extra yet unknown SURFEX variables
301	Fraction of each vegetation types on PATCH 1
302	Fraction of each vegetation types on PATCH 2
303	Fraction of each vegetation types cy43 (ECOCLIMAP-SG)
600	Physiography fields?
720	Sea ice
730	Sea ice (TICE_LL)
755	Precip
760	Sea
770	in addition to FLake (or instead of it)
780	Flake
790	Patch (*_P fields)
800	ISBA
810	Gridpoint average
820	Surface boundary multi layer fields
830	ISBA - patch 1 (X001*, open land)
840	ISBA - patch 2 (X002*, forest)
950	Town energy balance model (TEB)

A small selection of fields available in the SURFEX output files is shown below.

FA name	shortName	NCname	lvT	iOP	lev	sT	units	description
FRAC_SEA	#	sftof	hag	32	300	ins	0-1	Fraction of sea
FRAC_WATER	#	sftlaf	hag	33	300	ins	0-1	Fraction of water
FRAC_NATURE	#	sftnf	hag	34	300	ins	0-1	Fraction of nature
FRAC_TOWN	#	sfturf	hag	35	300	ins	0-1	Fraction of town
COVER001	#	lsm		1	0	ins		LAND SEA MASK
COVER002-COVER243	#	#		002-243	0	ins		ECOCLIMAP I cover types
COVER255	#	#		255	0	ins		ECOCLIMAP I MY_COVER type
COVER301-COVER573	#	#		001-254 & 001-019	0	ins		ECOCLIMAP II cover types
ZS	#	orog	hag	8	0	ins	m	Oro hgt.
SST	#	tos	ms	11	0	ins	K	SST
SIC	#	siconca	ms	91	0	ins	0-1	SIC
T2M_SEA	#	tas_sea	hag	11	760	ins	K	T2m sea
Q2M_SEA	#	huss_sea	hag	51	760	ins	kg kg-1	Q2m sea
MER10M_SEA	#	vas_sea	hag	34	760	ins	m s-1	V10m sea
ZON10M_SEA	#	uas_sea	hag	33	760	ins	m s-1	U10m sea
T2M_WAT	#	tas_water	hag	11	772	ins	K	T2m water
Q2M_WAT	#	huss_water	hag	51	770	ins	kg kg-1	Q2m water
MER10M_WAT	#	vas_water	hag	34	770	ins	m s-1	V10m water
ZON10M_WAT	#	uas_water	hag	33	770	ins	m s-1	U10m water
DSNTISBA	#	snd	hag	66	0	ins	m	Snow depth
WSNTISBA	#	snw	hag	13	0	ins	kg m-2	Total snow reservoir
T2M_ISBA	#	tas_nature	hag	11	802	ins	K	T2m isba
Q2M_ISBA	#	huss_nature	hag	51	802	ins	kg kg-1	Q2m isba
X001T2M_P	#	tas	hag	11	832	ins	K	T2m of patch 1
X002T2M_P	#	tas	hag	11	842	ins	K	T2m of patch 2
T2M_TEB	#	tas_town	hag	11	950	ins	K	T2m town
T2MMAX_TEB	#	tasmax_town	hag	15	950	max	K	Max Temp for town
T2MMIN_TEB	#	tasmin_town	hag	16	950	min	K	Min Temp for town
TG`L`	#	tg_L`LL`	hag	11	800+	ins	K	Temperature of soil layer `L`(isba)
WG`L`	#	wsa_L`LL`	hag	86	800+	ins	m3 m-3	Liquid volumetric water content of soil layer `L`
WGI`L`	#	isa_L`LL`	hag	193	800+	ins	m3 m-3	Frozen volumetric water content of soil layer `L`
WR	#	wr	hag	12	800	ins	kg m-2	Liquid water retained by foliage (isba)
DG`L`	#	dsoil_L`LL`	hag	23	300	ins	m	Soil depth of soil layer `L`