from pyiem.util import ncopen
from pyiem import era5land
import numpy as np
from tqdm import tqdm

dom = era5land.DOMAINS[""]

progress = tqdm(list(range(1950,2025)))
for year in progress:
    progress.set_description(f"{year}")
    with ncopen(f"{year}_era5land_hourly.nc", "a") as nc:
        nc.variables['lat'][:] = dom["YAXIS"]
        nc.variables['lon'][:] = dom["XAXIS"]
        nc.createDimension('nv', 2)
        nc.variables['lat'].bounds = "lat_bnds"
        nc.variables['lon'].bounds = "lon_bnds"
        lat_bnds = nc.createVariable("lat_bnds", float, ("lat", "nv"))
        lat_bnds[:, 0] = dom["YAXIS"] - 0.05
        lat_bnds[:, 1] = dom["YAXIS"] + 0.05
        lon_bnds = nc.createVariable("lon_bnds", float, ("lon", "nv"))
        lon_bnds[:, 0] = dom["XAXIS"] - 0.05
        lon_bnds[:, 1] = dom["XAXIS"] + 0.05