读取nc文件绘制固定时间和经纬度的位势高度图代码

以下是一份Python代码，可以读取nc文件中的数据，绘制指定时间和经纬度的位势高度图：

import netCDF4 as nc
import matplotlib.pyplot as plt
import numpy as np

# 读取nc文件
file = nc.Dataset('filename.nc', 'r')
lat = file.variables['lat'][:]
lon = file.variables['lon'][:]
time = file.variables['time'][:]
z = file.variables['z'][:]

# 根据需要选择时间和经纬度
t = 0 # 时间索引
lat_index = 10 # 纬度索引
lon_index = 20 # 经度索引

# 获取数据
Z = z[t, :, :]
X, Y = np.meshgrid(lon, lat)

# 绘制位势高度图
fig = plt.figure(figsize=(8, 6))
ax = fig.add_subplot(111, projection='3d')
ax.plot_surface(X, Y, Z, cmap='jet')
ax.set_xlabel('Longitude')
ax.set_ylabel('Latitude')
ax.set_zlabel('Geopotential Height')
plt.show()

在这份代码中，我们首先使用NetCDF4库读取nc文件中的数据，然后根据需要选择一个时间和一个经纬度的位置。再根据选择的位置获取相应的数据，最后使用Matplotlib库绘制位势高度图。需要注意的是，数据的维度顺序可能会因文件而异，需要根据具体文件的情况对代码进行修改。