require "numru/gphys"
include NumRu

gp = GPhys::IO.open "air.2012-01.nc","air"
gp = gp.cut("level"=>925..20, "lat"=>40)[true,{0..-1=>2},true]
     #^サブセット指定 (この時点ではまだ読み込まれない)

print "gp :\n"
p gp
print "pressure levels:\n"
p gp.coord("level").val
print "\ngp.mean(0) :\n"
p gp.mean(0)
print "\ngp.mean('lon','time') :\n"
p gp.mean('lon','time')
print "\ngp.mean :\n"
mn = gp.mean
p mn.class, mn
print "\nmean,max,min,stddev values:\n"
p gp.mean('lon','time').val, gp.max('lon','time').val, 
  gp.min('lon','time').val, gp.stddev('lon','time').val