python - 惑星境界層の高さをプロット

Question

Fortran 90またはpython（またはNCL）で惑星境界層の高さを地図（緯度/経度の2D）にプロットする簡単なプログラムがあるかどうか知りたいです。F.90 のプログラムを使用していますが、機能しないため、2 番目のプログラムと比較したいと思います。ありがとうございました

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

ファイル名、変数名などのデータに一致するように、この例を編集する必要がある場合があります。