lib/numru/gphys/gphys_grads_io.rb in gphys-1.1.1 vs lib/numru/gphys/gphys_grads_io.rb in gphys-1.2.2

- old
+ new

@@ -78,20 +78,21 @@ file = GrADS_Gridded.open(file) elsif ! file.is_a?(GrADS_Gridded) raise ArgumentError, "1st arg must be a GrADS_Gridded or a file name" end - grvar = file.var(varname) + grvar = GrADSVar.new(file,varname,true,false) data = VArrayGrADS.new(grvar) -# axposnames = [ "lon", "lat", "lev", "time" ] - axposnames = [ "x", "y", "z", "t" ] - rank = 4 + shape_current = grvar.shape_current + dimnames = grvar.dim_names + rank = dimnames.length bare_index = [ false ] * rank # will be true if coord var is not found axes = Array.new for i in 0...rank - axpos = VArrayGrADS.new( file.var(axposnames[i]) )[0...data.shape_current[i]] + grvar = GrADSVar.new(file,dimnames[i],false,true) + axpos = VArrayGrADS.new(grvar)[0...shape_current[i]] cell_center = true cell = false axis = Axis.new(cell,bare_index[i])