使用GDAL python读取WRF NetCDF文件中的1-D变量

Question

我的问题很简单 例如，使用wrfout文件“out.nc” 该文件包含Geo2D，Geo3D和1D变量。

在Python 2.7中使用 GDAL 包，我可以像这样轻松地提取Geo2D变量：

## T2 is 2-d variable means temperature 2 m above the ground
temp = gdal.Open('NETCDF:"'+"out.nc"+'":T2')          

但是，当我想使用此代码提取1d数组时，它失败了。

## Time is 1-d array represent the timeseries throught the simulation period
time = gdal.Open('NETCDF:"'+"out.nc"+'":Time')       

什么都没发生！希望有人提供一些建议，轻松阅读任意维度的WRF输出变量！

Answer 1

您也可以在scipy.io中使用NetCDF阅读器：

import scipy.io.netcdf as nc

# Open a netcdf file object and assign the data values to a variable
time = nc.netcdf_file('out.nc', 'r').variables['Time'][:]

它的好处是scipy是一个非常受欢迎且广泛安装的软件包，而在某些方面类似于打开文件。