Python 绘制行星边界层高度图

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

---

谢谢

这里是一个带有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()

---

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

您有什么类型的数据文件？有无数种方法可以做到这一点。你目前的做法是什么？是的，对不起，我没有提供足够的信息。我在netcdf中有一个文件，我使用理查森号码…你试了什么？您是否需要计算算法或可视化库？我希望这些算法能够从netcdf文件中绘制PBL高度图（内部有空气温度、表面压力、a和b混合坐标…）。我不知道是否有办法下载这个页面上的文件？在这种情况下，我很难在不展示的情况下解释我所做的尝试，所以这不是一个提出此类问题的好地方。PBL高度有多种定义，你应该查阅动力气象学教科书和相关论文。我的文件中没有pblh。。。我得计算一下！我正在寻找算法。。。