我正在尝试为靠近挪威 sauda 的 lwe_percipitation_rate 创建一个小子集，我目前正在努力使用 OPeNDAP 中的坐标系。
latitude_bounds = [57, 61]
longitude_bounds = [3.5,7]
lats = d.variables['lat'][:]
lons = d.variables['lon'][:]
#alt1
# latitude lower and upper index
latli = np.argmin( np.abs( lats - latitude_bounds[0] ) )
latui = np.argmin( np.abs( lats - latitude_bounds[1] ) )
# longitude lower and upper index
lonli = np.argmin( np.abs( lons - longitude_bounds[0] ) )
lonui = np.argmin( np.abs( lons - longitude_bounds[1] ) )
lwe_subset = d.variables['lwe_precipitation_rate'][ : , latli:latui , lonli:lonui ]
#alt2
lat_inds = np.where((lats > latitude_bounds[0]) & (lats < latitude_bounds[1]))
lon_inds = np.where((lons > longitude_bounds[0]) & (lons < longitude_bounds[1]))
lwe_subset = d.variables['lwe_precipitation_rate'][:,lat_inds,lon_inds]
这是目前我正在使用的两个替代方案，这是我感兴趣的区域： 红色方块标记了感兴趣的区域，靠近指定的坐标
