关于CDO守恒插值

cdo无法处理坐标系（指griddes中xval和yval）中含有missing value的格点。NCL同样不行（见这个ppt）。这种情况往往出现在曲线网格curvilinear中。

对于曲线网格的守恒插值，往往需要提供xbound和ybound。对于一些网格，特别是卫星数据，出现cdo remapcon (Abort): Source grid cell corner coordinates missing! 的原因往往出在这里。虽然可以计算得到（例如见How to calculate the corner points of a grid if centre points are know），但也可能会出现意料不到的问题。

好在cdo本身提供了投影网格（即gridtype = projection），为此往往需要以下变量（更具体的参见cdo user guide 1.5.2.4节）

gridtype  = projection
gridsize  = 
datatype  = float
xsize     = 
ysize     = 
xname     = x
xunits    = "meters"
yname     = y
yunits    = "meters"
xfirst    = 
xinc      = 
yfirst    = 
yinc      = 

其中，最后四个常数非常重要，需要针对不同的资料寻找对应的格点值。

对于写好的网格描述文件，可以使用setgrid将描述文件赋给netCDF，之后就可以使用remapcon插值了。

除此以外，cdo支持PROJ库语法，写在proj_params的内容会被直接传递到PROJ库。对于以下五种投影，也可以选择不显式地写出proj_params：

• rotated_latitude_longitude
• lambert_conformal_conic
• lambert_azimuthal_equal_area
• sinusoidal
• polar_stereographic