||
一些气候模式的输出数据的格式并非规则的经纬网格,而是网格大小不一的曲线网格。比如CESM模式的海洋模块使用的网格类型(grid type)为displaced-pole grid,意为把北极点移到了格陵兰岛。此网格的代号gx1v6,水平分辨率1.125° × 0.27~0.54°。
现在我们需要把数据重采样为0.25° × 0.25°的规则网格数据,使用CDO的步骤如下:
【1】生成权重文件
cdo genbil,r1440x720 Hs_mon_185001-201412.nc bilwgt_025deg.nc &
genbil表示使用bilinear双线性插值,Hs_mon_185001-201412.nc和bilwgt_025deg.nc分别是输入文件和生成的权重文件。
注意:执行这一步时,可能会报错:CDO genbil unsupported generic coordinates。这是说输入文件中的变量没有格网信息或者CDO不支持当前格网。这时需要把输入文件里保存坐标信息的变量(lon和lat)赋给要插值的变量(Hs)。命令如下:ncatted -a coordinates,Hs,c,c,"lon lat" Hs_mon_185001-201412.nc。ncatted需要安装NCO(NetCDF Operators)
才能执行。
另一种方法是使用CDO中的setgrid命令,具体可参考:https://stackoverflow.com/questions/67986451/converting-generic-grid-type-to-a-lonlat-gridtype-using-cdo
【2】使用权重文件,重采样为规则网格数据
cdo remap,r1440x720,bilwgt_025deg.nc Hs_mon_185001-201412.nc Hs_mon_185001-201412.025deg.nc
注意:这里使用remap命令,而非remapbil,因为权重文件bilwgt_025deg.nc中已经包含了插值方法。
参考文献:
https://www.ncl.ucar.edu/Support/talk_archives/2011/2634.html
https://code.mpimet.mpg.de/boards/1/topics/55
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-4-27 07:22
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社