import xarray as xr import numpy as np # Set xarray DataArray metadata for the NetCDF file ds_dims = ['lat', 'lon'] ds_coords = {'lat': np.arange(-90, 90+1, 1), 'lon': np.arange(0, 359+1, 1)} # Create a DataArray and write it to a NetCDF file (assuming data has a shape of (num_lats x num_lons)) ds = xr.DataArray(data, coords=ds_coords, dims=ds_dims, attrs={'long_name': '2m Temperature'}, name='tmean').to_dataset() ds.to_netcdf('out.nc', encoding={var: dict(zlib=True, complevel=9) for var in ds.data_vars})