Metview
import metview as mv
IO
import zipfile
import cdsapi
c = cdsapi.Client()
c.retrieve(
"cems-fire-historical",
{
"format": "zip",
"product_type": "reanalysis",
"variable": [
"danger_risk",
],
"version": "3.1",
"dataset": "Consolidated dataset",
"year": "2020",
"month": "05",
"day": "15",
},
"fire.zip",
)
with zipfile.ZipFile("fire.zip", "r") as zip_ref:
zip_ref.extractall(".")
fire = mv.read("ECMWF_FWI_DANGER_RISK_20200515_1200_hr_v3.1_con.nc")
Write grib file
fs.write("file.grib")
Read grib file
fieldset = mv.read("file.grib")
Fieldset (grib object)
fieldset.describe()
fieldset.ls()
ds = fieldset.to_dataset()
fieldset = mv.dataset_to_fieldset(ds)
t2m = fs.select(shortName="2t")
Regrid
t2m_0x25 = mv.read(data=t2m, grid=[0.25, 0.25])
t2m_N400 = mv.read(data=t2m_0x25, grid="N400")
fs = mv.read("T2D01260600012606001")
fs_N640 = mv.read(data=fs, grid="N640")
fs_N640.write("T2D01260600012606001_N640")