绘图行星边界层高度

Question

我想知道是否有人在fortran 90或python（或NCL）中绘制地图上的行星边界层高度（lat / lon中的2D）的简单程序。 我在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()

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