我正在尝试使用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更合适。
发布于 2020-10-22 02:54:37
由于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.nchttps://stackoverflow.com/questions/64464455
复制相似问题