如何使用 NCO 提取单个变量。关联变量问题答案

【问题标题】：How to extract a single variable using NCO. Associated variable problem如何使用 NCO 提取单个变量。关联变量问题
【发布时间】：2021-02-04 10:05:28
【问题描述】：

我正在尝试使用 NCO 从 NetCDF 文件（超过 thredds）中提取单个变量 kd_490。

我的代码如下：

ncks -v kd_490 -d lat,40.0,70.0 -d lon,-20.0,15.0 https://rsg.pml.ac.uk/thredds/dodsC/cci/v4.2-release/geographic/daily/kd/1998/ESACCI-OC-L3S-K_490-MERGED-1D_DAILY_4km_GEO_PML_KD490_Lee-19980102-fv4.2.nc out.nc

但是，与 kd_490 一起，它还提取 kd_490_bias 和 kd_490_rmsd。我知道 ncks 提取“关联”变量。见here。但是，我不清楚为什么 NCO 将这些标识为“关联”变量。

我想不出只选择变量 kd_490 的方法。在代码中添加“-C”会导致网格错误。有人知道怎么做吗？

注意：CDO 可以解决这个问题。但是，CDO 在代码中提取空间子集的效率较低，因此 NCO 更合适。

【问题讨论】：

标签： netcdf nco

【解决方案1】：

由于kd_490 中的ancillary_variables 属性，NCO 将kd_490_bias 和kd_490_rmsd 识别为关联变量：

   float kd_490(time,lat,lon) ;
      kd_490:_FillValue = 9.96921e+36f ;
      kd_490:long_name = "Downwelling attenuation coefficient at 490nm, derived using Lee 2005 equation and bbw from Zhang 2009 (following the SeaDAS Kd_lee algorithm)" ;
      kd_490:units = "m-1" ;
      kd_490:ancillary_variables = "kd_490_rmsd kd_490_bias" ;
      kd_490:grid_mapping = "crs" ;
      kd_490:parameter_vocab_uri = "http://vocab.ndg.nerc.ac.uk/term/P071/19/CFSN0064" ;
      kd_490:standard_name = "volume_attenuation_coefficient_of_downwelling_radiative_flux_in_sea_water" ;
      kd_490:units_nonstandard = "m^-1" ;
      kd_490:_ChunkSizes = 1, 270, 270 ;

如记录的here。要提取 kd_490 而不使用辅助变量，但使用网格变量（这是其他“关联”变量），这对我有用：

ncks -O -C -v kd_490,crs,lat,lon -d lat,40.0,70.0 -d lon,-20.0,15.0 https://rsg.pml.ac.uk/thredds/dodsC/cci/v4.2-release/geographic/daily/kd/1998/ESACCI-OC-L3S-K_490-MERGED-1D_DAILY_4km_GEO_PML_KD490_Lee-19980102-fv4.2.nc ~/foo2.nc

【讨论】：