require "numru/gphys"
include NumRu


gphys = GPhys::NetCDF_IO.open('intr_SPRGC_glevel5_alpha05PI.nc', 'h')
outfile = NetCDF.create('tmp.nc')

t_axis = gphys.axis('t')
t_a = VArray.new(t_axis.pos.val, {"long_name"=>"time", "units"=>"seconds since 2010-11-9 00:00:00 +0:00"}, "t")
t = Axis.new.set_pos(t_a)

level_a = VArray.new(NArray.sfloat(1).indgen(1,1), {"long_name"=>"level", "units"=>"m"}, "Z")
level = Axis.new.set_pos(level_a)

shape = gphys.data.val.shape
#data = NArray.float(shape[0],shape[1],shape[2], 1)
#p data
#data[true,true,true,0] = gphys.data.val

data = gphys.data.copy
data.reshape!(shape[0],shape[1], 1, shape[2])
p data.val[true,true,0,2]

outgphys = GPhys.new(Grid.new(gphys.axis('lon'), gphys.axis('lat'), level, t), data)

GPhys::NetCDF_IO.write(outfile, outgphys)
outfile.close