绘制行星边界层高度

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()

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