netcdf using xarray

xr_nc.py

"""
Created on Wed May 19 17:55:14 2021

@author: mathew

"""
import xarray as xr
import numpy as np


# making a random spatial dataset
# 12 time-step and 30x30 grid points


vals = np.random.rand(12,30,30)
time = np.arange(12)
lats = np.linspace(10,20,30)
lon = np.linspace(70,80,30)

# creating xarray dataset

dset = xr.Dataset({
    "values" : (["time","lon","lat"],vals),},
    coords = {"time":(["time",],time),
              "lon":(["lon",],lon),
              "lat":(["lat",],lats)})


dset.to_netcdf("file.nc")