plot行星邊界層高度

Question

我想知道是否有人有一個簡單的程序在fortran 90或python（或NCL）的地圖上繪製行星邊界層高度（2D經緯度）。 我正在使用F.90中的程序，但它不工作，所以我想與第二個程序進行比較。 謝謝

Answer 1

下面是一個Python的例子與matplotlib和netCDF4模塊：

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

您可能需要修改這個例子來匹配數據，例如文件名，變量名等