GuoKeCCC的个人博客分享 http://blog.sciencenet.cn/u/GuoKeCCC

博文

CDO将曲线网格(curve linear)重采样为经纬网格(lonlat)

已有 5250 次阅读 2022-4-29 16:52 |个人分类:CDO|系统分类:科研笔记

一些气候模式的输出数据的格式并非规则的经纬网格,而是网格大小不一的曲线网格。比如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



https://m.sciencenet.cn/blog-3486942-1336245.html

上一篇:CDO (Climate Data Operators)安装遇到的问题及解决
下一篇:CDO每隔n段时间(小时,天,月份等)统计变量(求和,平均,最大最小等)

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-27 07:22

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部