basic sample

# 基礎的な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()