1

我想将纬度范围添加到现有的 netCDF 文件中。见http://www.cgd.ucar.edu/cms/eaton/netcdf/CF-20010629.htm#bnds边界的定义。

有没有可以自动添加边界的命令？我对使用 nco 或 netCDF4 或 cdms 工具包的任何解决方案持开放态度。但是，我更希望它是一个后处理解决方案。

4

1 回答 1

3

将边界添加到现有坐标，如下所示：

ncap2 -O -s 'crd@bounds="crd_bnds";defdim("bnds",2);crd_bnds[$crd,$bnds]=0.0;*crd_dff=0.5*(crd(1)-crd(0));crd_bnds(:,0)=crd-crd_dff;crd_bnds(:,1)=crd+crd_dff;' in.nc out.nc

这产生

zender@aerosol:~$ ncks --cdl -v crd ~/foo2.nc
netcdf foo2 {
  dimensions:
    bnds = 2 ;
    crd = 10 ;

  variables:
    int crd(crd) ;
      crd:bounds = "crd_bnds" ;

    double crd_bnds(crd,bnds) ;

  data:
    crd = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 ;

    crd_bnds = -0.5, 0.5, 0.5, 1.5, 1.5, 2.5, 2.5, 3.5, 3.5, 4.5, 4.5, 5.5, 5.5, 6.5, 6.5, 7.5, 7.5, 8.5, 8.5, 9.5 ;

} // group /

这里坐标被命名为“crd”，边界被命名为“crd_bnds”。当坐标改变一个固定值（如上）时，边界值的公式可以很简单，或者更多地涉及到可变间距（读者练习）。我们将在不久的将来添加一个 ncap2 函数来自动执行此操作。

cz

于 2016-04-17T18:04:13.317 回答