#[derive(Debug, Clone)]
pub struct Grib1Parameter {
pub number: u8,
pub abbreviation: &'static str,
pub name: &'static str,
pub units: &'static str,
}
pub fn get_parameter(center_id: u8, parameter: u8) -> Option<Grib1Parameter> {
match center_id {
98 => get_ecmwf_parameter(parameter), 7 => get_ncep_parameter(parameter), _ => get_wmo_standard_parameter(parameter), }
}
fn get_ecmwf_parameter(parameter: u8) -> Option<Grib1Parameter> {
let param = match parameter {
1 => Grib1Parameter {
number: 1,
abbreviation: "sp",
name: "Surface pressure",
units: "Pa",
},
2 => Grib1Parameter {
number: 2,
abbreviation: "prmsl",
name: "Pressure reduced to MSL",
units: "Pa",
},
11 => Grib1Parameter {
number: 11,
abbreviation: "t",
name: "Temperature",
units: "K",
},
29 => Grib1Parameter {
number: 29,
abbreviation: "sst",
name: "Sea surface temperature",
units: "K",
},
31 => Grib1Parameter {
number: 31,
abbreviation: "u",
name: "U component of wind",
units: "m s-1",
},
32 => Grib1Parameter {
number: 32,
abbreviation: "v",
name: "V component of wind",
units: "m s-1",
},
33 => Grib1Parameter {
number: 33,
abbreviation: "u",
name: "U velocity",
units: "m s-1",
},
34 => Grib1Parameter {
number: 34,
abbreviation: "v",
name: "V velocity",
units: "m s-1",
},
39 => Grib1Parameter {
number: 39,
abbreviation: "w",
name: "Vertical velocity",
units: "Pa s-1",
},
51 => Grib1Parameter {
number: 51,
abbreviation: "q",
name: "Specific humidity",
units: "kg kg-1",
},
52 => Grib1Parameter {
number: 52,
abbreviation: "r",
name: "Relative humidity",
units: "%",
},
53 => Grib1Parameter {
number: 53,
abbreviation: "q",
name: "Humidity mixing ratio",
units: "kg kg-1",
},
54 => Grib1Parameter {
number: 54,
abbreviation: "pwat",
name: "Precipitable water",
units: "kg m-2",
},
59 => Grib1Parameter {
number: 59,
abbreviation: "prate",
name: "Precipitation rate",
units: "kg m-2 s-1",
},
61 => Grib1Parameter {
number: 61,
abbreviation: "tp",
name: "Total precipitation",
units: "m",
},
71 => Grib1Parameter {
number: 71,
abbreviation: "tcc",
name: "Total cloud cover",
units: "%",
},
129 => Grib1Parameter {
number: 129,
abbreviation: "z",
name: "Geopotential",
units: "m2 s-2",
},
130 => Grib1Parameter {
number: 130,
abbreviation: "t",
name: "Temperature",
units: "K",
},
131 => Grib1Parameter {
number: 131,
abbreviation: "u",
name: "U component of wind",
units: "m s-1",
},
132 => Grib1Parameter {
number: 132,
abbreviation: "v",
name: "V component of wind",
units: "m s-1",
},
133 => Grib1Parameter {
number: 133,
abbreviation: "q",
name: "Specific humidity",
units: "kg kg-1",
},
134 => Grib1Parameter {
number: 134,
abbreviation: "sp",
name: "Surface pressure",
units: "Pa",
},
135 => Grib1Parameter {
number: 135,
abbreviation: "w",
name: "Vertical velocity",
units: "Pa s-1",
},
141 => Grib1Parameter {
number: 141,
abbreviation: "sd",
name: "Snow depth",
units: "m of water equivalent",
},
151 => Grib1Parameter {
number: 151,
abbreviation: "prmsl",
name: "Pressure reduced to MSL",
units: "Pa",
},
157 => Grib1Parameter {
number: 157,
abbreviation: "r",
name: "Relative humidity",
units: "%",
},
165 => Grib1Parameter {
number: 165,
abbreviation: "u10",
name: "10 metre U wind component",
units: "m s-1",
},
166 => Grib1Parameter {
number: 166,
abbreviation: "v10",
name: "10 metre V wind component",
units: "m s-1",
},
167 => Grib1Parameter {
number: 167,
abbreviation: "t2m",
name: "2 metre temperature",
units: "K",
},
168 => Grib1Parameter {
number: 168,
abbreviation: "d2m",
name: "2 metre dewpoint temperature",
units: "K",
},
_ => return None,
};
Some(param)
}
fn get_ncep_parameter(parameter: u8) -> Option<Grib1Parameter> {
get_wmo_standard_parameter(parameter)
}
fn get_wmo_standard_parameter(parameter: u8) -> Option<Grib1Parameter> {
let param = match parameter {
1 => Grib1Parameter {
number: 1,
abbreviation: "pres",
name: "Pressure",
units: "Pa",
},
2 => Grib1Parameter {
number: 2,
abbreviation: "prmsl",
name: "Pressure reduced to MSL",
units: "Pa",
},
7 => Grib1Parameter {
number: 7,
abbreviation: "gh",
name: "Geopotential height",
units: "gpm",
},
11 => Grib1Parameter {
number: 11,
abbreviation: "t",
name: "Temperature",
units: "K",
},
33 => Grib1Parameter {
number: 33,
abbreviation: "u",
name: "U-component of wind",
units: "m s-1",
},
34 => Grib1Parameter {
number: 34,
abbreviation: "v",
name: "V-component of wind",
units: "m s-1",
},
39 => Grib1Parameter {
number: 39,
abbreviation: "w",
name: "Vertical velocity",
units: "Pa s-1",
},
51 => Grib1Parameter {
number: 51,
abbreviation: "q",
name: "Specific humidity",
units: "kg kg-1",
},
52 => Grib1Parameter {
number: 52,
abbreviation: "r",
name: "Relative humidity",
units: "%",
},
61 => Grib1Parameter {
number: 61,
abbreviation: "tp",
name: "Total precipitation",
units: "kg m-2",
},
_ => return None,
};
Some(param)
}
pub fn get_level_type_info(level_type: u8) -> (&'static str, &'static str) {
match level_type {
1 => ("surface", ""),
2 => ("cloud_base", ""),
3 => ("cloud_top", ""),
4 => ("isotherm_zero", "m"),
100 => ("isobaric", "hPa"),
102 => ("mean_sea_level", ""),
103 => ("fixed_height", "m"),
105 => ("fixed_height_above_ground", "m"),
106 => ("sigma", "sigma"),
107 => ("sigma_height", "sigma"),
108 => ("sigma_pressure", "sigma"),
109 => ("hybrid", "hybrid"),
111 => ("depth_below_surface", "m"),
112 => ("layer_between_depths", "m"),
113 => ("isentropic", "K"),
114 => ("layer_between_isentropic", "K"),
200 => ("entire_atmosphere", ""),
201 => ("entire_ocean", ""),
_ => ("unknown", ""),
}
}
#[cfg(test)]
mod tests {
use super::*;
#[test]
fn test_ecmwf_parameters() {
let param = get_parameter(98, 11).unwrap();
assert_eq!(param.abbreviation, "t");
assert_eq!(param.name, "Temperature");
let param = get_parameter(98, 131).unwrap();
assert_eq!(param.abbreviation, "u");
}
#[test]
fn test_level_types() {
let (name, units) = get_level_type_info(100);
assert_eq!(name, "isobaric");
assert_eq!(units, "hPa");
let (name, _) = get_level_type_info(1);
assert_eq!(name, "surface");
}
}