Regridding introduces missing values at wraparound longitude

JiaweiZhuang / xESMF

Universal Regridder for Geospatial Data

http://xesmf.readthedocs.io/

MIT License

269 stars 49 forks source link

Closed spencerahill closed 6 years ago

spencerahill commented 6 years ago

Entirely possible I'm making a silly mistake, but I can't figure it out so far at least. Thanks in advance for help.

JiaweiZhuang commented 6 years ago

Oops, xESMF has some issue with periodic boundary conditions right now. Sorry for the trouble.

There are two quick ways to fix that:

Use conservative regridding. It uses cell corner information and won't miss any cells.
Manually extend one more grid point in the longitude direction. For example, let lon = [-180, -175..., 175, 180] and use the same data at -180 and 180

spencerahill commented 6 years ago

Thanks. I will try one or both of these and let you know how it goes.

JiaweiZhuang commented 6 years ago

Periodic boundary option is now included in v0.1.1. Just set regridder = xe.Regridder(..., periodic=True)

The option is used in this demo.