require "numru/gphys"
include NumRu

ncfn = "../test_a3_L080km/lee-wave_test_a3_VelX_rank??????.nc"
vname = "VelX"
outfn = "data_tmp/a3_L080km/lee-wave_test_a3_VelX.nc"

#ncfn = "../test_a3_L080km/lee-wave_test_a3_VelZ_rank??????.nc"
#vname = "VelZ"
#outfn = "data_tmp/a3_L080km/lee-wave_test_a3_VelZ.nc"

#vname = "PTemp"
#ncfn = "../test_a3_L080km/lee-wave_test_a3_PTemp_rank??????.nc"
#outfn = "data_tmp/a3_L080km/lee-wave_test_a3_PTemp.nc"

vname = "VelX"
sym   = "a3_L080km"
ncfn  = "../test_"+sym+"/lee-wave_test_a3_"+vname+"_rank??????.nc"
outfn = "data_tmp/"+sym+"/lee-wave_test_a3_"+vname+".nc"



url = ncfn + "@" + vname
gphys = GPhys::IO.open_gturl( url )

#< iterate over the last dimension >
outfile = NetCDF.create( outfn )
GPhys::NetCDF_IO.each_along_dims_write(gphys, outfile, -1) do |sub|
  [sub]
end
outfile.close