使用 ncks 的 4D netcdf 变量的 Hyperslab
Hyperslab of a 4D netcdf variable using ncks
我有一个很大的 netcdf 文件，我只需要其中的某些数据。因此，我想使用 ncks 创建此 netcdf 文件的细分。
netcdf文件如下：
Source:
F:\LECOB\Model091208_195356.nc
Format:
64bit
Global Attributes:
Model = 's26.bobshelf.20141113'
Title = 'S-NWM_BiP'
Dimensions:
ni_t = 682
nj_t = 712
nk_t = 29
time = 1 (UNLIMITED)
ni_w = 682
nj_w = 712
nk_w = 30
ni_u = 681
nj_u = 712
nk_u = 29
ni_v = 682
nj_v = 711
nk_v = 29
ni_f = 681
nj_f = 711
nk_f = 29
Variables:
time
Size: 1x1
Dimensions: time
Datatype: double
Attributes:
units = 'seconds from 2009-dec-08 17:00:14'
long_name = 'time'
standard_name = 'time'
time_origin = '2009-dec-08 17:00:14'
calendar = 'gregorian'
content = 'T'
axis = 'T'
associate = 'undefined'
ni_t
Size: 682x1
Dimensions: ni_t
Datatype: single
Attributes:
long_name = 'x_grid_index'
standard_name = 'x_grid_index'
content = 'X'
axis = 'X'
nj_t
Size: 712x1
Dimensions: nj_t
Datatype: single
Attributes:
long_name = 'y_grid_index'
standard_name = 'y_grid_index'
content = 'Y'
axis = 'Y'
nk_t
Size: 29x1
Dimensions: nk_t
Datatype: single
Attributes:
long_name = 'z_grid_index'
standard_name = 'z_grid_index'
content = 'Z'
axis = 'Z'
positive = 'up'
ni_w
Size: 682x1
Dimensions: ni_w
Datatype: single
Attributes:
long_name = 'x_grid_index_at_w_location'
standard_name = 'x_grid_index_at_w_location'
content = 'X'
axis = 'X'
nj_w
Size: 712x1
Dimensions: nj_w
Datatype: single
Attributes:
long_name = 'y_grid_index_at_w_location'
standard_name = 'y_grid_index_at_w_location'
content = 'Y'
axis = 'Y'
nk_w
Size: 30x1
Dimensions: nk_w
Datatype: single
Attributes:
long_name = 'z_grid_index_at_w_location'
standard_name = 'z_grid_index_at_w_location'
content = 'Z'
axis = 'Z'
positive = 'up'
ni_u
Size: 681x1
Dimensions: ni_u
Datatype: single
Attributes:
long_name = 'x_grid_index_at_u_location'
standard_name = 'x_grid_index_at_u_location'
content = 'X'
axis = 'X'
nj_u
Size: 712x1
Dimensions: nj_u
Datatype: single
Attributes:
long_name = 'y_grid_index_at_u_location'
standard_name = 'y_grid_index_at_u_location'
content = 'Y'
axis = 'Y'
nk_u
Size: 29x1
Dimensions: nk_u
Datatype: single
Attributes:
long_name = 'z_grid_index_at_u_location'
standard_name = 'z_grid_index_at_u_location'
content = 'Z'
axis = 'Z'
positive = 'up'
ni_v
Size: 682x1
Dimensions: ni_v
Datatype: single
Attributes:
long_name = 'x_grid_index_at_v_location'
standard_name = 'x_grid_index_at_v_location'
content = 'X'
axis = 'X'
nj_v
Size: 711x1
Dimensions: nj_v
Datatype: single
Attributes:
long_name = 'y_grid_index_at_v_location'
standard_name = 'y_grid_index_at_v_location'
content = 'Y'
axis = 'Y'
nk_v
Size: 29x1
Dimensions: nk_v
Datatype: single
Attributes:
long_name = 'z_grid_index_at_v_location'
standard_name = 'z_grid_index_at_v_location'
content = 'Z'
axis = 'Z'
positive = 'up'
ssh
Size: 682x712x1
Dimensions: ni_t,nj_t,time
Datatype: single
Attributes:
units = 'm'
long_name = 'sea surface height above geoid'
standard_name = 'sea_surface_height_above_geoid'
content = 'TYX'
associate = 'time latitude_t longitude_t'
coordinates = 'time latitude_t longitude_t'
_FillValue = -9999
CFL2D
Size: 682x712x1
Dimensions: ni_t,nj_t,time
Datatype: single
Attributes:
units = 's'
long_name = 'CFL2D'
standard_name = 'CFL2D'
content = 'TYX'
associate = 'time latitude_t longitude_t'
coordinates = 'time latitude_t longitude_t'
_FillValue = -9999
tem
Size: 682x712x29x1
Dimensions: ni_t,nj_t,nk_t,time
Datatype: single
Attributes:
units = 'degrees_Celsius'
long_name = 'sea_water_potential_temperature'
standard_name = 'sea_water_potential_temperature'
content = 'TZYX'
associate = 'time depth_t latitude_t longitude_t'
coordinates = 'time depth_t latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
sal
Size: 682x712x29x1
Dimensions: ni_t,nj_t,nk_t,time
Datatype: single
Attributes:
units = '1e-3'
long_name = 'sea water salinity'
standard_name = 'sea_water_salinity'
content = 'TZYX'
associate = 'time depth_t latitude_t longitude_t'
coordinates = 'time depth_t latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
vel_u
Size: 681x712x29x1
Dimensions: ni_u,nj_u,nk_u,time
Datatype: single
Attributes:
units = 'm s-1'
long_name = 'sea_water_x_velocity_at_u_location'
standard_name = 'sea_water_x_velocity_at_u_location'
content = 'TZYX'
associate = 'time depth_u latitude_u longitude_u'
coordinates = 'time depth_u latitude_u longitude_u'
_FillValue = -9999
positive = 'up'
vel_v
Size: 682x711x29x1
Dimensions: ni_v,nj_v,nk_v,time
Datatype: single
Attributes:
units = 'm s-1'
long_name = 'sea_water_y_velocity_at_v_location'
standard_name = 'sea_water_y_velocity_at_v_location'
content = 'TZYX'
associate = 'time depth_v latitude_v longitude_v'
coordinates = 'time depth_v latitude_v longitude_v'
_FillValue = -9999
positive = 'up'
kh
Size: 682x712x30x1
Dimensions: ni_w,nj_w,nk_w,time
Datatype: single
Attributes:
units = 'm2/s'
long_name = 'kh'
standard_name = 'kh'
content = 'TZYX'
associate = 'time depth_w latitude_t longitude_t'
coordinates = 'time depth_w latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
tken
Size: 682x712x30x1
Dimensions: ni_w,nj_w,nk_w,time
Datatype: single
Attributes:
units = '(m/s)2'
long_name = 'tken'
standard_name = 'tken'
content = 'TZYX'
associate = 'time depth_w latitude_t longitude_t'
coordinates = 'time depth_w latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
w
Size: 682x712x30x1
Dimensions: ni_w,nj_w,nk_w,time
Datatype: single
Attributes:
units = 'm s-1'
long_name = 'vertical_sea_water_velocity_at_w_location'
standard_name = 'vertical_sea_water_velocity_at_w_location'
content = 'TZYX'
associate = 'time depth_w latitude_t longitude_t'
coordinates = 'time depth_w latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
现在我对 4D 变量感兴趣 vel_u (ni_u,nj_u,nk_u,time)。
我想提取 ni_u 151 到 152，nj_u 234 到 235，nk_u 一切和 time 一切。
这个问题帮助了我 and so did the following link http://nco.sourceforge.net/nco.html#crd
使用这些，我在 linux 计算机上尝试了以下代码：
ncks -C -F -d vel_u,151,152,1 20091208_195356.nc test.nc
这给了我两个问题：
- 它复制了所有变量，而不仅仅是 vel_u，即使我按照问题
中的建议使用了 -C
- 我不知道如何指定只使用维度
nj_u 的 234 到 235
那么如何将 vel_u 变量的那些部分（ni_u 151 到 152，nj_u 234 到 235）放入我的 test.nc 文件中？
非常感谢任何答案！
我想我找到了答案：-d 是维度，-v 是变量，所以我的问题的答案如下：
ncks -C -F -d nj_u,234,235,1 -d ni_u,151,152,1 -v vel_u 20091208_195356.nc test.nc
-C 确保只复制变量 vel_u
-F 因为提示应该从 1 而不是 0 开始
-d 名称维度、最小值、最大值、步骤
-v 名称变量
input.nc
output.nc
这似乎正是我想要的。
我有一个很大的 netcdf 文件，我只需要其中的某些数据。因此，我想使用 ncks 创建此 netcdf 文件的细分。 netcdf文件如下：
Source:
F:\LECOB\Model091208_195356.nc
Format:
64bit
Global Attributes:
Model = 's26.bobshelf.20141113'
Title = 'S-NWM_BiP'
Dimensions:
ni_t = 682
nj_t = 712
nk_t = 29
time = 1 (UNLIMITED)
ni_w = 682
nj_w = 712
nk_w = 30
ni_u = 681
nj_u = 712
nk_u = 29
ni_v = 682
nj_v = 711
nk_v = 29
ni_f = 681
nj_f = 711
nk_f = 29
Variables:
time
Size: 1x1
Dimensions: time
Datatype: double
Attributes:
units = 'seconds from 2009-dec-08 17:00:14'
long_name = 'time'
standard_name = 'time'
time_origin = '2009-dec-08 17:00:14'
calendar = 'gregorian'
content = 'T'
axis = 'T'
associate = 'undefined'
ni_t
Size: 682x1
Dimensions: ni_t
Datatype: single
Attributes:
long_name = 'x_grid_index'
standard_name = 'x_grid_index'
content = 'X'
axis = 'X'
nj_t
Size: 712x1
Dimensions: nj_t
Datatype: single
Attributes:
long_name = 'y_grid_index'
standard_name = 'y_grid_index'
content = 'Y'
axis = 'Y'
nk_t
Size: 29x1
Dimensions: nk_t
Datatype: single
Attributes:
long_name = 'z_grid_index'
standard_name = 'z_grid_index'
content = 'Z'
axis = 'Z'
positive = 'up'
ni_w
Size: 682x1
Dimensions: ni_w
Datatype: single
Attributes:
long_name = 'x_grid_index_at_w_location'
standard_name = 'x_grid_index_at_w_location'
content = 'X'
axis = 'X'
nj_w
Size: 712x1
Dimensions: nj_w
Datatype: single
Attributes:
long_name = 'y_grid_index_at_w_location'
standard_name = 'y_grid_index_at_w_location'
content = 'Y'
axis = 'Y'
nk_w
Size: 30x1
Dimensions: nk_w
Datatype: single
Attributes:
long_name = 'z_grid_index_at_w_location'
standard_name = 'z_grid_index_at_w_location'
content = 'Z'
axis = 'Z'
positive = 'up'
ni_u
Size: 681x1
Dimensions: ni_u
Datatype: single
Attributes:
long_name = 'x_grid_index_at_u_location'
standard_name = 'x_grid_index_at_u_location'
content = 'X'
axis = 'X'
nj_u
Size: 712x1
Dimensions: nj_u
Datatype: single
Attributes:
long_name = 'y_grid_index_at_u_location'
standard_name = 'y_grid_index_at_u_location'
content = 'Y'
axis = 'Y'
nk_u
Size: 29x1
Dimensions: nk_u
Datatype: single
Attributes:
long_name = 'z_grid_index_at_u_location'
standard_name = 'z_grid_index_at_u_location'
content = 'Z'
axis = 'Z'
positive = 'up'
ni_v
Size: 682x1
Dimensions: ni_v
Datatype: single
Attributes:
long_name = 'x_grid_index_at_v_location'
standard_name = 'x_grid_index_at_v_location'
content = 'X'
axis = 'X'
nj_v
Size: 711x1
Dimensions: nj_v
Datatype: single
Attributes:
long_name = 'y_grid_index_at_v_location'
standard_name = 'y_grid_index_at_v_location'
content = 'Y'
axis = 'Y'
nk_v
Size: 29x1
Dimensions: nk_v
Datatype: single
Attributes:
long_name = 'z_grid_index_at_v_location'
standard_name = 'z_grid_index_at_v_location'
content = 'Z'
axis = 'Z'
positive = 'up'
ssh
Size: 682x712x1
Dimensions: ni_t,nj_t,time
Datatype: single
Attributes:
units = 'm'
long_name = 'sea surface height above geoid'
standard_name = 'sea_surface_height_above_geoid'
content = 'TYX'
associate = 'time latitude_t longitude_t'
coordinates = 'time latitude_t longitude_t'
_FillValue = -9999
CFL2D
Size: 682x712x1
Dimensions: ni_t,nj_t,time
Datatype: single
Attributes:
units = 's'
long_name = 'CFL2D'
standard_name = 'CFL2D'
content = 'TYX'
associate = 'time latitude_t longitude_t'
coordinates = 'time latitude_t longitude_t'
_FillValue = -9999
tem
Size: 682x712x29x1
Dimensions: ni_t,nj_t,nk_t,time
Datatype: single
Attributes:
units = 'degrees_Celsius'
long_name = 'sea_water_potential_temperature'
standard_name = 'sea_water_potential_temperature'
content = 'TZYX'
associate = 'time depth_t latitude_t longitude_t'
coordinates = 'time depth_t latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
sal
Size: 682x712x29x1
Dimensions: ni_t,nj_t,nk_t,time
Datatype: single
Attributes:
units = '1e-3'
long_name = 'sea water salinity'
standard_name = 'sea_water_salinity'
content = 'TZYX'
associate = 'time depth_t latitude_t longitude_t'
coordinates = 'time depth_t latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
vel_u
Size: 681x712x29x1
Dimensions: ni_u,nj_u,nk_u,time
Datatype: single
Attributes:
units = 'm s-1'
long_name = 'sea_water_x_velocity_at_u_location'
standard_name = 'sea_water_x_velocity_at_u_location'
content = 'TZYX'
associate = 'time depth_u latitude_u longitude_u'
coordinates = 'time depth_u latitude_u longitude_u'
_FillValue = -9999
positive = 'up'
vel_v
Size: 682x711x29x1
Dimensions: ni_v,nj_v,nk_v,time
Datatype: single
Attributes:
units = 'm s-1'
long_name = 'sea_water_y_velocity_at_v_location'
standard_name = 'sea_water_y_velocity_at_v_location'
content = 'TZYX'
associate = 'time depth_v latitude_v longitude_v'
coordinates = 'time depth_v latitude_v longitude_v'
_FillValue = -9999
positive = 'up'
kh
Size: 682x712x30x1
Dimensions: ni_w,nj_w,nk_w,time
Datatype: single
Attributes:
units = 'm2/s'
long_name = 'kh'
standard_name = 'kh'
content = 'TZYX'
associate = 'time depth_w latitude_t longitude_t'
coordinates = 'time depth_w latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
tken
Size: 682x712x30x1
Dimensions: ni_w,nj_w,nk_w,time
Datatype: single
Attributes:
units = '(m/s)2'
long_name = 'tken'
standard_name = 'tken'
content = 'TZYX'
associate = 'time depth_w latitude_t longitude_t'
coordinates = 'time depth_w latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
w
Size: 682x712x30x1
Dimensions: ni_w,nj_w,nk_w,time
Datatype: single
Attributes:
units = 'm s-1'
long_name = 'vertical_sea_water_velocity_at_w_location'
standard_name = 'vertical_sea_water_velocity_at_w_location'
content = 'TZYX'
associate = 'time depth_w latitude_t longitude_t'
coordinates = 'time depth_w latitude_t longitude_t'
_FillValue = -9999
positive = 'up'
现在我对 4D 变量感兴趣 vel_u (ni_u,nj_u,nk_u,time)。
我想提取 ni_u 151 到 152，nj_u 234 到 235，nk_u 一切和 time 一切。
这个问题帮助了我
ncks -C -F -d vel_u,151,152,1 20091208_195356.nc test.nc
这给了我两个问题：
- 它复制了所有变量，而不仅仅是 vel_u，即使我按照问题
- 我不知道如何指定只使用维度
nj_u 的 234 到 235
-C
那么如何将 vel_u 变量的那些部分（ni_u 151 到 152，nj_u 234 到 235）放入我的 test.nc 文件中？
非常感谢任何答案！
我想我找到了答案：-d 是维度，-v 是变量，所以我的问题的答案如下：
ncks -C -F -d nj_u,234,235,1 -d ni_u,151,152,1 -v vel_u 20091208_195356.nc test.nc
-C 确保只复制变量 vel_u
-F 因为提示应该从 1 而不是 0 开始
-d 名称维度、最小值、最大值、步骤
-v 名称变量
input.nc
output.nc
这似乎正是我想要的。