# 基礎的なpy3gradsの例
# Basic sample of py3grads.
from py3grads import Grads
ga = Grads(verbose=False)
indir='/data/REANALYSES/JRA55/org_Monthly/anl_p125/'
ctl_fn='anl_p125_hgt.monthly.ctl'
ctl_fn_fl=indir + ctl_fn
ga('open ' + ctl_fn_fl)
ga('set lat 20 30') # 緯度経度の設定
ga('set lon 180 200')
hgt=ga.exp('hgtprs')
lons=ga.exp('lon')
lats=ga.exp('lat')
from matplotlib import pyplot as plt
plt.contourf(lons,lats,hgt)
plt.show()