Regrid using xesmf latlon
Fri May 27 2022 11:35:53 GMT+0000 (Coordinated Universal Time)
Saved by
@diptish
#python
import xarray as xr
import xesmf as xe
import numpy as np
def regrid(dr_in):
"""Convenience function for one-time regridding"""
ds_out = xr.Dataset(
{
"lat": (["lat"], np.arange(-90, 91, 1.0)),
"lon": (["lon"], np.arange(-180, 182, 2.0)),
}
)
#ds_out
regridder = xe.Regridder(dr_in, ds_out, 'bilinear', periodic=True)
dr_out = regridder(dr_in)
return dr_out
ERAin=xr.open_dataset('ERA5_monthly_SST.nc')
ERAOut=regrid(ERAin)
ERAOut
ERAOut.to_netcdf(path='ERA5_SST.nc',mode='w',format='NETCDF4')
content_copyCOPY
Comments