Source code for pycif.plugins.domains.lmdz_ico.read_domain

import numpy as np
import xarray as xr




def read_grid(self, **kwargs):
    """Reads a grid from NetCDF files"""
    self.unstructured_domain = True

    with xr.open_dataset(self.grid_path) as ds:
        self.nlat = 1
        self.nlon = ds.sizes["cell"]

        lon_vertices = ds[ds.lon.attrs["bounds"]].values
        lat_vertices = ds[ds.lat.attrs["bounds"]].values

        self.zlat = ds.lat.values[np.newaxis, :]
        self.zlon = ds.lon.values[np.newaxis, :]
        self.zlatc = lat_vertices.T
        self.zlonc = lon_vertices.T

        if "area" in ds:
            self.areas = ds.area.values

    with xr.open_dataset(self.vcoord_path) as ds:
        self.nlev = ds.sizes["lev"] - 1
        self.sigma_a = ds.ap.values
        self.sigma_b = ds.bp.values
        self.pressure_unit = ds.ap.attrs["units"]