gphys 1.5.6 → 1.5.7

data/lib/numru/gphys/grads_gridded.rb

@@ -340,6 +340,7 @@ module NumRu
       "cray_32bit_ieee"=>nil,
       "365_day_calendar"=>nil,
       "pdef"=>nil,
+      "edef"=>nil,
     }
 
     case(@mode)
@@ -380,7 +381,7 @@ module NumRu
   end
 
   def ndims
-    4
+    @options["edef"] ? 5 : 4
   end
 
   def nvars
@@ -480,7 +481,7 @@ module NumRu
     end
     @title = "<no title>" if @title==""
     @undef = -999.0 if @undef==nil
-    return <<EOS
+    ctl = <<EOS
 DSET    #{if @dset[0]=="/" then @dset else "^"+@dset end}
 TITLE   #{@title}
 UNDEF   #{@undef}
@@ -489,15 +490,23 @@ XDEF    #{@dimensions[0][:len]} #{@dimensions[0][:flag]} #{@dimensions[0][:spec]
 YDEF    #{@dimensions[1][:len]} #{@dimensions[1][:flag]} #{@dimensions[1][:spec]}
 ZDEF    #{@dimensions[2][:len]} #{@dimensions[2][:flag]} #{@dimensions[2][:spec]}
 TDEF    #{@dimensions[3][:len]} #{@dimensions[3][:flag]} #{@dimensions[3][:spec]}
+EOS
+    if ( @options["edef"] )
+    ctl += <<EOS
+EDEF    #{@dimensions[4][:len]} #{@dimensions[4][:flag]} #{@dimensions[4][:spec]}
+EOS
+    end
+    ctl += <<EOS
 VARS    #{nvars}
 #{@variables.collect{|i| i[:name]+" "+i[:nlev].to_s+" "+i[:option].to_s+
        " "+i[:description]}.join("\n")}
 ENDVARS
 EOS
+    return ctl
   end
 
   def inspect
-    return <<EOS
+    str = <<EOS
 #{self.class}
 file: #{@ctlfile.path}
 DSET #{@dset}
@@ -506,11 +515,19 @@ XDIM #{@dimensions[0].inspect}
 YDIM #{@dimensions[1].inspect}
 ZDIM #{@dimensions[2].inspect}
 TDIM #{@dimensions[3].inspect}
+EOS
+    if ( @options["edef"] )
+      str += <<EOS
+EDEF #{@dimensions[4].inspect}
+EOS
+    end
+    str += <<EOS
 VARS
 #{@variables.collect{|i| "  "+i[:name]+"\t"+i[:nlev].to_s+"\t"+i[:option].to_s+
        "\t"+i[:description]}.join("\n")}
 ENDVARS
 EOS
+    return str
   end
 
   def put(ary)
@@ -530,7 +547,7 @@ EOS
 
   end
 
-  def get(name, z, t, lonlat=nil)
+  def get(name, z, t, lonlat=nil, ens=nil)
 
 # t:    [0,1,2,..] : record number
 # time: [0,1/24,2/24,...] days since 00:00Z01jan2000 : Numeric with units
@@ -593,6 +610,9 @@ EOS
         init_min = target_date.min
         @dset_r.gsub!( /%n2/, sprintf("%02d",init_min) )
       end
+      if ( ens && @dset_r =~ /%e(\d+)/ )
+        @dset_r.gsub!( /%e(\d+)/, sprintf("%0#{$1}d",ens) )
+      end
 
       # find t in the first record in the target file (="init_t")
       # and determine the record number in the target file
@@ -609,10 +629,10 @@ EOS
 
       t_in_target_file = t - init_t
 
-      start_byte = start_byte(name, z, t_in_target_file)
+      start_byte = start_byte(name, z, t_in_target_file, 0)
       @datafile = File.open(@dset_r,"rb")
     else
-      start_byte = start_byte(name, z, t)
+      start_byte = start_byte(name, z, t, ens)
       @datafile = File.open(@dset,"rb")
     end
 
@@ -777,9 +797,8 @@ EOS
         else
           raise "Invalid line: "+line
         end
-      when /^\s*[XYZT]DEF/i
-#        /^\s*([XYZT])DEF\s+(\d+)\s+(\S+)\s+(.*)$/i =~ line
-        /^\s*([XYZT])DEF\s+(\d+)\s+(\S+)+(.*)$/i =~ line
+      when /^\s*[XYZTE]DEF/i
+        /^\s*([XYZTE])DEF\s+(\d+)\s+(\S+)+(.*)$/i =~ line
         if ( (len=$2) && (flag=$3) && (spec=$4))
           dim = {:len=>len.to_i, :flag=>flag, :spec=>spec}
           case $1
@@ -802,8 +821,13 @@ EOS
             dim[:name] = 't'
             dim[:description] = 'time'
             idim=3
+          when /E/i
+            dim[:name] = 'e'
+            dim[:description] = 'ensemble'
+            idim=4
+            @options["edef"] = true
           end
-          if (idim!=3)
+          if (idim<3)
             if (dim[:flag] =~ /LINEAR/i)
               begin
                 dim[:start],dim[:increment] = spec.split.collect!{|i| i.to_f}
@@ -833,7 +857,7 @@ EOS
             else
               raise "invalid or not-yet-supported dimension flag: "+dim[:flag]
             end
-          else
+          elsif idim==3
             # idim = 3 --- time
             if (dim[:flag] =~ /LINEAR/i)
               start,increment= spec.split
@@ -845,6 +869,8 @@ EOS
             else
               raise "invalid dimension flag(only LINEAR is available for time)"
             end
+          else # idim==4 ensemble
+            dim[:levels] = NArray.sint(dim[:len]).indgen
           end
           @dimensions[idim]=dim
         else
@@ -944,6 +970,7 @@ EOS
         sy = sy*PI/180
         sy2 = sy2*PI/180
         xi = NArray.sfloat(imax).indgen(1) - bi
+        bj = jmax - bj if @options['yrev']
         yj = NArray.sfloat(jmax).indgen(1) - bj
         x = xi*dx
         #y = -yj*dy 
@@ -980,7 +1007,6 @@ EOS
         lon = NArray.sfloat(isize.to_i).indgen(1)
         lat = NArray.sfloat(jsize.to_i).indgen(1)
       end
-      p lon, lat
 
       dim = {:len=>isize.to_i, :flag=>"levels"}
       dim[:name] = 'x'
@@ -1012,7 +1038,7 @@ EOS
 
   end
 
-  def start_byte(name, level, time)
+  def start_byte(name, level, time, ens=0)
     # offset to read an xy section of the variable with NAME
     # at LEVEL(counted from 0) at TIME(conted from 0)
     if (map = @map[name] )
@@ -1023,6 +1049,12 @@ EOS
         raise "Time #{time} is not in the data period"
       end
       iblock = map[:start]+level*map[:zstep]+time*map[:tstep]
+      if @options["edef"]
+        if (ens<0 || ens>=@dimensions[4][:len])
+          raise "Ens #{ens} is not in the data"
+        end
+        iblock += ens * map[:tstep] * @dimensions[3][:len]
+      end
       str_byte = map[:offset] + @dimensions[0][:len] * @dimensions[1][:len] * iblock
       if( @options["sequential"] )
         str_byte += iblock*2 + 4
@@ -1243,6 +1275,7 @@ EOS
           @shape = [@dimensions[0][:len],@dimensions[1][:len],
                     @dimensions[2][:len],@dimensions[3][:len]]
           @shape[2] = var[:nlev].to_i if( var[:nlev] )
+          @shape[4] = @dimensions[4][:len] if @rank==5
 
           @attr[:missing_value] = NArray.sfloat(1).fill!(@ctl.undef)
           found = true
@@ -1414,8 +1447,8 @@ EOS
         end
       end
     else
-      h_sta = [0, 0, 0, 0]
-      h_end = [-1, -1, -1, -1]
+      h_sta = [0, 0, 0, 0, 0]
+      h_end = [-1, -1, -1, -1, -1]
       if hash == nil
 #      elsif hash.key?("index")==true
       else
@@ -1439,12 +1472,27 @@ EOS
       str_z = idx_z[h_sta[2]]; end_z = idx_z[h_end[2]]
       str_t = idx_t[h_sta[3]]; end_t = idx_t[h_end[3]]
 
-      na = NArrayMiss.new(vartype,end_x-str_x+1,end_y-str_y+1,end_z-str_z+1,end_t-str_t+1)
-      for t in str_t..end_t
-        for z in str_z..end_z
-          na_xy = @ctl.get(@varname,z,t,[str_x,end_x,str_y,end_y])
-          mask = ( na_xy.ne @ctl.undef )
-          na[true,true,z-str_z,t-str_t] = NArrayMiss.to_nam(na_xy,mask)
+      if @rank == 5 # has ensemble dimension
+        idx_e = NArray.int(@shape[4]).indgen!
+        str_e = idx_e[h_sta[4]]; end_e = idx_e[h_end[4]]
+        na = NArrayMiss.new(vartype,end_x-str_x+1,end_y-str_y+1,end_z-str_z+1,end_t-str_t+1,end_e-str_e+1)
+        for e in str_e..end_e
+          for t in str_t..end_t
+            for z in str_z..end_z
+              na_xy = @ctl.get(@varname,z,t,[str_x,end_x,str_y,end_y],e)
+              mask = ( na_xy.ne @ctl.undef )
+              na[true,true,z-str_z,t-str_t,e-str_e] = NArrayMiss.to_nam(na_xy,mask)
+            end
+          end
+        end
+      else
+        na = NArrayMiss.new(vartype,end_x-str_x+1,end_y-str_y+1,end_z-str_z+1,end_t-str_t+1)
+        for t in str_t..end_t
+          for z in str_z..end_z
+            na_xy = @ctl.get(@varname,z,t,[str_x,end_x,str_y,end_y])
+            mask = ( na_xy.ne @ctl.undef )
+            na[true,true,z-str_z,t-str_t] = NArrayMiss.to_nam(na_xy,mask)
+          end
         end
       end
     end
@@ -1479,7 +1527,7 @@ EOS
      stride = Array.new
      set_stride = false
      a.each{|i|
-       if(i.is_a?(Fixnum))
+       if(i.is_a?(Integer))
          first.push(i)
          last.push(i)
          stride.push(1)
@@ -1511,7 +1559,7 @@ EOS
            if n==0 then
              k = at
              if at > 0
-               a[0..at-1].each{|x| if x.is_a?(Fixnum) then k -= 1 end}
+               a[0..at-1].each{|x| if x.is_a?(Integer) then k -= 1 end}
              end
              shape_tmp = na_tmp.shape
              shape_tmp[k] = i.length
@@ -1523,7 +1571,7 @@ EOS
          end
          return na
        else
-         raise TypeError, "argument must be Fixnum, Range, Hash, TrueClass, Array, or NArray"
+         raise TypeError, "argument must be Integer, Range, Hash, TrueClass, Array, or NArray"
        end
      }
 
@@ -1534,7 +1582,7 @@ EOS
      end
      shape = na.shape
      (a.length-1).downto(0){ |i|
-         shape.delete_at(i) if a[i].is_a?(Fixnum)
+         shape.delete_at(i) if a[i].is_a?(Integer)
       }
       na.reshape!( *shape )
      na
@@ -1548,7 +1596,7 @@ EOS
      stride = Array.new
      set_stride = false
      a.each{|i|
-       if(i.is_a?(Fixnum))
+       if(i.is_a?(Integer))
          first.push(i)
          last.push(i)
          stride.push(1)
@@ -1575,13 +1623,13 @@ EOS
          at = a.index(i)
          i = NArray.to_na(i) if i.is_a?(Array)
          val = NArray.to_na(val) if val.is_a?(Array)
-         rank_of_subset = a.dup.delete_if{|v| v.is_a?(Fixnum)}.length
+         rank_of_subset = a.dup.delete_if{|v| v.is_a?(Integer)}.length
          if val.rank != rank_of_subset
            raise "rank of the rhs (#{val.rank}) is not equal to the rank "+
                  "of the subset specified by #{a.inspect} (#{rank_of_subset})"
          end
          k = at
-         a[0..at-1].each{|x| if x.is_a?(Fixnum) then k -= 1 end}
+         a[0..at-1].each{|x| if x.is_a?(Integer) then k -= 1 end}
          if i.length != val.shape[k]
            raise "length of the #{k+1}-th dim of rhs is incorrect "+
                  "(#{i.length} for #{val.shape[k]})"
@@ -1598,7 +1646,7 @@ EOS
          end
          return self
        else
-         raise TypeError, "argument must be Fixnum, Range, Hash, TrueClass, Array, or NArray"
+         raise TypeError, "argument must be Integer, Range, Hash, TrueClass, Array, or NArray"
        end
      }
 

data/lib/numru/gphys/grib.rb

@@ -1229,7 +1229,7 @@ class GPhys
           nlon,nlat = shape
           (nlon*nlat)==nlen+nmiss || raise("length is not collect")
           if index
-            index = index.collect{|el| Fixnum===el ? el..el : el}
+            index = index.collect{|el| Integer===el ? el..el : el}
             il = index[0]
             i0 = il.first
             i1 = il.last
@@ -1980,7 +1980,7 @@ even when the time is unique (i.e., length of the time coordinate is 1).
                 (0..dlen-1).each{|i| imap[n-2][i] = i}
               end
               next
-            when Fixnum
+            when Integer
               sha[n] = 1
               imap[n-2][ind] = 0
               dimdel.push(n)

data/lib/numru/gphys/grib_params.rb

@@ -706,8 +706,8 @@ NAMES_UNITS = {
 #150	COVMZ	m2/s2		Covariance between meridional and zonal components of the wind. Defined as [uv]-[u][v], where "[]" indicates the mean over the indicated time span.		
 #151	COVTZ1	K*m/s		Covariance between temperature and zonal component of the wind. Defined as [uT]-[u][T], where "[]" indicates the mean over the indicated time span.		
 #152	COVTM	K*m/s		Covariance between temperature and meridional component of the wind. Defined as [vT]-[v][T], where "[]" indicates the mean over the indicated time span.		
-  PARAM_CLWMR => ["Cloud water", "CLWMR", "Kg/kg", nil],
-  PARAM_O3MR => ["Ozone mixing ratio", "O3MR", "Kg/kg", nil],
+  PARAM_CLWMR => ["Cloud water", "CLWMR", "kg/kg", nil],
+  PARAM_O3MR => ["Ozone mixing ratio", "O3MR", "kg/kg", nil],
 #155	GFLUX	W/m2		Ground Heat Flux		
   PARAM_CIN => ["Convective inhibition", "CIN", "J/kg", nil],
   PARAM_CAPE => ["Convective Available Potential Energy", "CAPE", "J/kg", nil],		
@@ -941,7 +941,7 @@ PARAM_CWORK => ["Cloud work function", "CWORK", "J/kg", nil],
 PARAM_FGLU => ["U sfc momentum flux by long gravity wave", "FGLU", "N/m^2", nil],
 PARAM_FGLV => ["V sfc momentum flux by long gravity wave", "FGLV", "N/m^2", nil],
 PARAM_ADVUA => ["Advective zonal acceleration", "ADVUA", "m/s/day", nil],
-PARAM_VWV => ["Column total of meridional water vapor flux", "VWV", "Kg/m/s", nil],
+PARAM_VWV => ["Column total of meridional water vapor flux", "VWV", "kg/m/s", nil],
 PARAM_FGSV => ["V sfc momentum flux by short gravity wave", "FGSV", "N/m^2", nil],
 PARAM_GFLX => ["Ground heat flux", "GFLX", "W/m^2", nil],
 PARAM_UWV => ["Column total of zonal water vapor flux", "UWV", "kg/m/s", nil],

data/lib/numru/gphys/interpolate.rb

@@ -210,6 +210,85 @@ module NumRu
       interpolate(*coords)
     end
 
+    # Regrid for equally spaced grid (or index based regrid)
+    # 
+    # ARGUMENTS
+    # * axes [Array of Axis (or 1D VArray)] : if 1D VArray, an Axis is created
+    # * poses [Array of multi-D VArray (or GPhys)] : specifies where to
+    #   regrid. lengths of poses and axes must agree. The shape of each of pos
+    #   must be equal to the shape suggested by axes (i.e., corrected lengths
+    #   of axes). poses must have non-overwrapping names.
+    #   Each pos (an element of poses) must have a name equal to
+    #   a main coordinate varible (or a coordinate variable plus "_index").
+    #   E.g., if self has axes named ["x", "y", "t"], the name of a pos
+    #   must be one of "x", "y", "t", "x_index", "y_index", "t_index".
+    #   No two pos must share the coordinate name (w or wo _index).
+    def regrid2( axes, poses )
+      if axes.length != poses.length
+        raise(ArgumentError, "lengths of axes and poses do not agree")
+      end
+      axes = axes.map{|x| x.is_a?(VArray) ? Axis.new.set_pos(x) : x }
+      newshape = axes.map{|x| x.length}
+      dims = Array.new
+      cyclics = Array.new
+      axnms = axnames
+      eps = 1e-4
+      ndims = poses.length
+      scidxs = poses.map do |pos|
+        if pos.shape != newshape
+          raise(ArgumentError, "lengths of axes and poses do not agree")
+        end
+        nm = pos.name
+        index_type = false  # initialization
+        d = axnms.index(nm)
+        if !d
+          d = axnms.index(nm.sub(/_index$/,''))
+          raise("pos #{nm} does not match any of the axis names") if !d
+          index_type = true
+        end
+        dims.push(d)
+        cyclics.push( axis(d).cyclic_extendible? )
+        v = pos.val
+        if v.is_a?(NArrayMiss)
+          if v.get_mask.count_false == 0
+            v = v.to_na
+          else
+            raise ArgumentError, "pos cannot have data missing"
+          end
+        end
+        if index_type
+          v
+        else
+          cval = coord(d).val
+          c0 = cval[0]
+          dc = (cval[-1] - cval[0])/(shape[d]-1.0)
+          if (cval[1]-cval[0] - dc).abs > eps*dc
+            raise("coordinate #{nm} is not equally spaced." +
+                  " Use index based spacification if possible.")
+          end
+          (v - c0)/dc
+        end
+      end
+      na = val
+      if na.is_a?(NArrayMiss)
+        mask = na.get_mask
+        na = na.to_na
+      else
+        mask = NArray.byte( *na.shape )
+        mask = mask.fill!(1)
+      end
+      
+      rval, rmask = c_regrid2_w_idx(na, mask, dims, scidxs, cyclics)
+
+      newdata = VArray.new( NArrayMiss.to_nam(rval, rmask), @data, name )
+      newaxes = (0...rank).map{|d| axis(d)}
+      (0...dims.length).each do |d|
+        newaxes[ dims[d] ] = axes[d]
+      end
+      newgrid = Grid.new( *newaxes )
+      GPhys.new( newgrid, newdata )
+    end
+    
     # Reverse the main data (i.e., the dependent variable) and one of the
     # coordinates (an independent variable) through interpolation.
     # 
@@ -424,6 +503,50 @@ module NumRu
       ret
     end
 
+    # interpolate data missing alon dim's dimention
+    #
+    # ARGUMENTS:
+    # * dim [Integer or String] specifies the dimension
+    #
+    # * force_cyclic [nil, true, or Numeric] : if nil (default),
+    #   cyclic interpolation is made only if the cyclic_extendible?
+    #   applied to the dimension returns true (that is, the axis
+    #   is cyclic and only one grid-point extension is needed).
+    #   By setting true, cyclic extension is made as long as modulo
+    #   is defined, irrespective of the remaining distances. By
+    #   setting a Numeric, cyclic extension is forced by using its
+    #   value as the modulo.
+    def interpolate_missing(dim, force_cyclic: nil)
+      dim = dim_index(dim)
+      ax = axis(dim)
+      modulo = nil
+      case force_cyclic
+      when nil
+        cyclic = ax.cyclic_extendible?
+        modulo = ax.modulo if cyclic
+      when true
+        modulo = ax.modulo
+        if modulo
+          cyclic = true
+          modulo = modulo[0]
+        else
+          cyclic = false
+        end
+      when Numeric
+        cyclic = true
+        modulo = force_cyclic
+      else
+        raise(ArgumentError, "force_cyclic must be nil, true, or Numeric")
+      end
+      v = self.val
+      na = v.to_na
+      mask = v.get_mask
+      x = ax.pos.val
+      ve, mask = c_interpo_missing(na,mask,dim,x,cyclic,modulo)
+      dat = VArray.new( NArrayMiss.to_nam(ve,mask), self.data, self.name )
+      GPhys.new(grid, dat)
+    end
+    
     private
 
     def _interpo_find_position(crdmap)
@@ -987,7 +1110,36 @@ if $0 == __FILE__
   GGraph.tone gpcut,true,"int"=>0.15,"min"=>-1.05,"max"=>1.05
   GGraph.color_bar
   
-
+  GGraph::tone gdm,true,"color_bar"=>true, "title"=>"orig with missing"
+  p gdm.coord(1).att_names, gdm.coord(1).val.to_a
+  gde = gdm.interpolate_missing(1)
+  GGraph::tone gde,true,"color_bar"=>true, "title"=>"missing interpolated (1, non-cyclic)"
+  gde = gdm.interpolate_missing(1, force_cyclic: 2*Math::PI)
+  GGraph::tone gde,true,"color_bar"=>true, "title"=>"missing interpolated (1, cyclic)"
+  gde = gdm.interpolate_missing(0, force_cyclic: 2*Math::PI)
+  GGraph::tone gde,true,"color_bar"=>true, "title"=>"missing interpolated (0, cyclic)"
+
+  nth = 40
+  nr = 40
+  theta = NArray.float(nth).indgen! * (Math::PI/2/(nth-1))
+  #p "@@@@",gde.coord(0).val, gde.coord(1).val
+  r = NArray.float(nr).indgen! * (8.0/nr)
+  thax = Axis.new.set_pos( VArray.new(theta,
+                 {"long_name"=>"azimuth","units"=>"radian"}, "theta") )
+  rax = Axis.new.set_pos( VArray.new(r,
+                 {"long_name"=>"r","units"=>"m"}, "r") )
+  x = r.newdim(-1) * NMath.cos( theta.newdim(0) ) + gde.coord(0).val[0]
+  y = r.newdim(-1) * NMath.sin( theta.newdim(0) ) + gde.coord(1).val[0]
+  vx = VArray.new(x, {"long_name"=>"x","units"=>"m"}, "x")
+  vy = VArray.new(y, {"long_name"=>"y","units"=>"m"}, "y")
+
+  GGraph::tone gdd[false,1],true,"color_bar"=>true
+  gdr = gdd.regrid2( [rax, thax], [vx, vy] )
+  GGraph::tone gdr[false,1],true,"color_bar"=>true, "keep"=>true, #"tonc"=>true,
+               "title"=>"polar projected"
+  #GGraph::contour gdr,false
+  p "%%%", gdd, gdr
+  
   #< finish >
   DCL.grcls
 end

data/lib/numru/gphys/varraycomposite.rb

@@ -149,7 +149,9 @@ module NumRu
       # p @bound_idx, @varrays
     end
 
-    undef initialize_mapping
+    [:initialize_mapping, :mapping, :varray, :ary, :reshape].each do |method|
+      undef_method(method) if method_defined?(method)
+    end
 
     def set_att(name,val)
       @varrays.each{|va|
@@ -158,8 +160,6 @@ module NumRu
       self
     end
 
-    undef :mapping, :varray, :ary
-
     def inspect
       "<#{self.class.to_s} shape=#{shape.inspect} #_of_tiles=#{@bound_idx.collect{|b| b.length-1}.inspect}>"
       # "bounds=#{@bound_idx.collect{|b| b[1..-2]}.inspect}"
@@ -203,6 +203,9 @@ module NumRu
 	    raise TypeError, "Unexpected type #{v.class}"
 	  end
 	end
+        if val.is_a?(NArray) && v.is_a?(NArrayMiss)
+          val = NArrayMiss.to_nam(val)
+        end
 	val[*vidx] = v
       }
       val
@@ -286,7 +289,7 @@ module NumRu
     #def reshape!( *shape )
     #  raise "cannot reshape a #{class}. Use copy first to make it a VArray with NArray"
     #end
-    undef reshape!
+    #undef reshape!
 
     ## Not needed, so far:
     #def convention

data/lib/numru/gphys/version.rb

@@ -1,6 +1,6 @@
 module NumRu
   class GPhys
     # Add subsubminor version while under development; remove it and advance version to release
-    VERSION = "1.5.6"
+    VERSION = "1.5.7"
   end
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: gphys
 version: !ruby/object:Gem::Version
-  version: 1.5.6
+  version: 1.5.7
 platform: ruby
 authors:
 - Takeshi Horinouchi
@@ -12,7 +12,7 @@ authors:
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2019-03-19 00:00:00.000000000 Z
+date: 2025-10-13 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: narray
@@ -562,6 +562,7 @@ files:
 - lib/numru/ganalysis/fitting.rb
 - lib/numru/ganalysis/histogram.rb
 - lib/numru/ganalysis/log_p.rb
+- lib/numru/ganalysis/lomb_scargle.rb
 - lib/numru/ganalysis/met.rb
 - lib/numru/ganalysis/met_z.rb
 - lib/numru/ganalysis/narray_ext.rb
@@ -630,6 +631,7 @@ files:
 - testdata/UV.jan.nc
 - testdata/assoc_crds.nc
 - testdata/cira86.dat
+- testdata/pres.jan.nc
 homepage: http://www.gfd-dennou.org/arch/ruby/products/gphys/
 licenses:
 - BSD-2-Clause
@@ -651,8 +653,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubyforge_project: 
-rubygems_version: 2.5.2.1
+rubygems_version: 3.2.5
 signing_key: 
 specification_version: 4
 summary: A multi-purpose class to handle Gridded Physical quantities
@@ -671,3 +672,4 @@ test_files:
 - testdata/UV.jan.nc
 - testdata/assoc_crds.nc
 - testdata/cira86.dat
+- testdata/pres.jan.nc