JiaweiZhuang / xESMF

Universal Regridder for Geospatial Data
http://xesmf.readthedocs.io/
MIT License
269 stars 49 forks source link

Kernel Died Error during Weight Creation with Finer Resolutions #92

Open AhmadAmmar opened 4 years ago

AhmadAmmar commented 4 years ago

I am able to run the following code successfully if I specify a grid resolution size of 0.5, 0.25 and even 0.1. However, I get a kernel died error if I specify finer resolutions like 0.099, 0.066 or 0.033. What am I doing wrong?

ds = xr.open_dataset('/home/splash/data/metoc_output_ammar/metoc_output/ncep/grtofs/rtofs_glo_2ds_3hrly_uv_20200416.nc').isel(time=0)
dr = ds['uo']

ds['lat'].values, ds['lon'].values

ds_out = xr.Dataset({'lat': (['lat'], np.arange(-90, 90, 0.066)),
                     'lon': (['lon'], np.arange(-180, 180, 0.066)),
                    }
                   )

regridder = xe.Regridder(ds, ds_out, 'bilinear', reuse_weights=True)
dr_out = regridder(dr)
dr_out.to_netcdf('rtofs_uo_0_066.nc')