PyGrADSは、GrADSの変数をPythonメモリに簡単にインポートできます.

from pylab import *
from grads import GrADS

# Start GrADS and open the data file
# ----------------------------------
ga = GrADS(Bin='grads',Echo=False,Port=True,Window=False)
ga.open('../data/model.ctl')

# XY slices
# ---------
ts1 = ga.exp('ts')   #ts1.shape=(46,73)
ts2 = ga.expr('ts')  #ts2.shape=(46,73)
print("XY Skin temperature: ")
print(ts1.data-ts2.data)

# XYT slices
# ----------
ga('set t 2 3')
ts1 = ga.exp('ts') #ts1.shape=(2,46,73)
ts2 = ga.expr('ts')#ts2.shape=(2,46,73)
print("XYT Skin temperature: ")
print(ts1.data-ts2.data)

# XYZ slices
# ----------
ga('set t 3')
ga('set z 5 6')
xx1 = ga.exp('ua')#xx1.shape=(2, 46, 73)
xx2 = ga.expr('ua')#xx2.shape=(2, 46, 73)
print("XYZ Zonal Wind")
print(xx1.data-xx2.data)

# XZT slices
# ----------
ga('set t 2 4')
ga('set z 1 5')
ga('set y 20')
xx2 = ga.expr('ua')#xx2.shape=(3,5,73)
print("XZT Zonal Wind")
print(xx2)