跳到主要内容

前言

以前用 grads 和 Fortran 时候,经常会用到 grd 二进制文件,并拿 ctl 文件画图

可以用 Python 的 xgrads 包来实现这些操作

下载 xgrads

pip install xgrads

官方网站:https://xgrads.readthedocs.io/en/latest/

注意事项

在下载 xgrads 之后,就可以使用了,使用之前,请确保你的 ctl 文件和 grd 文件正确。

示例

  1. 打开多个空间场相同,时间不同的 ctl 文件时,可以批量读取

    类似于 xarray 的 open_mfdataset

    from xgrads import open_mfdataset
    ds=open_mfdataset('/../*.ctl')
    ds
  2. 转 nc 格式

    from xgrads import open_mfdataset
    ds=open_dataset('/../xxxx.ctl')
    ds.attrs['pdef']='None' # 原来的 ds 的 attr 是没有 pdef 这个属性的,所以如果报错,可以加上这个,大致报错的原因都是因为属性缺少有些东西,缺啥补啥
    ds.to_netcdf('/../xx.nc')

更多的使用方法,具体参考官网