python - 绘制行星边界层高

Question

我想知道是否有人有一个简单的程序可以在 fortran 90 或 python（或 NCL）的地图（纬度/经度中的二维）上绘制行星边界层高度。我在 F.90 中使用了一个程序，但它不起作用，所以我想与第二个程序进行比较。谢谢

Answer 1

这是一个带有 matplotlib 和 netCDF4 模块的 Python 示例：

import matplotlib.pyplot as plt
from netCDF4 import Dataset

nc = Dataset('mydatafile.nc','r')
lon = nc.variables['lon'][:] 
lat = nc.variables['lat'][:] 
pblh = nc.variables['pblh'][:]
nc.close() 

plt.contourf(lon,lat,pblh)
plt.colorbar()
plt.savefig('pblh.png')
ptl.clf()

您可能需要编辑此示例以匹配您的数据，例如文件名、变量名等。