rb-grib 0.1.4 → 0.2.0

data/lib/numru/grib/setenv.rb

@@ -0,0 +1,12 @@
+def_path = File.join(File.expand_path(File.dirname(__FILE__)),"grib","definitions")
+if File.exist?(def_path)
+  env_name = 'GRIB_DEFINITION_PATH'
+  if gdp = ENV[env_name]
+    ENV[env_name] = "#{gdp}:#{def_path}"
+  else
+    path = `grib_info -d`
+    if File.exist?(path)
+      ENV[env_name] = "#{path}:#{def_path}"
+    end
+  end
+end

data/lib/numru/grib/version.rb

@@ -1,5 +1,5 @@
 module NumRu
   class Grib
-    VERSION = "0.1.4"
+    VERSION = "0.2.0"
   end
 end

data/rb-grib.gemspec

@@ -21,5 +21,6 @@ Gem::Specification.new do |s|
 
   # specify any dependencies here; for example:
   s.add_development_dependency "rspec"
-  # s.add_runtime_dependency "rest-client"
+  s.add_runtime_dependency "narray"
+  s.add_runtime_dependency "narray_miss"
 end

data/spec/grib_read_spec.rb

@@ -1,83 +1,104 @@
 require File.expand_path(File.join('.','spec_helper'), File.dirname(__FILE__))
 
 dir = File.expand_path(File.join('..','data'), File.dirname(__FILE__))
-fname_reg_gaus_surf_g1 = File.join(dir, 'regular_gaussian_surface.grib1')
-fname_reg_gaus_surf_g2 = File.join(dir, 'regular_gaussian_surface.grib2')
-fname_reg_gaus_ml_g1 = File.join(dir, 'regular_gaussian_model_level.grib1')
-fname_reg_gaus_ml_g2 = File.join(dir, 'regular_gaussian_model_level.grib2')
-fname_reg_gaus_pl_g1 = File.join(dir, 'regular_gaussian_pressure_level.grib1')
-fname_reg_gaus_pl_g2 = File.join(dir, 'regular_gaussian_pressure_level.grib2')
-fname_reg_gaus_pl_const_g1 = File.join(dir, 'regular_gaussian_pressure_level_constant.grib1')
-fname_reg_gaus_pl_const_g2 = File.join(dir, 'regular_gaussian_pressure_level_constant.grib2')
-fname_reg_ll_surf_g1 = File.join(dir, 'regular_latlon_surface.grib1')
-fname_reg_ll_surf_g2 = File.join(dir, 'regular_latlon_surface.grib2')
-fname_reg_ll_surf_const_g1 = File.join(dir, 'regular_latlon_surface_constant.grib1')
-fname_reg_ll_surf_const_g2 = File.join(dir, 'regular_latlon_surface_constant.grib2')
+GribData = {
+  File.join(dir, 'regular_gaussian_surface.grib1') =>
+  {"2t"=>{:dims=>["lon","lat"], :shape=>[128,64], :val=>{[0,0]=>2.4552479553e2}}
+  },
+  File.join(dir, 'regular_gaussian_surface.grib2') =>
+  {"2t"=>{:dims=>["lon","lat"], :shape=>[128,64], :val=>{[0,0]=>2.4552479553e2}}
+  },
+  File.join(dir, 'regular_gaussian_model_level.grib1') =>
+  {"t"=>{:dims=>["lon","lat"], :shape=>[128,64], :val=>{[0,0]=>1.9907820129e2}}
+  },
+  File.join(dir, 'regular_gaussian_model_level.grib2') =>
+  {"t"=>{:dims=>["lon","lat"], :shape=>[128,64], :val=>{[0,0]=>1.9907820129e2}}
+  },
+  File.join(dir, 'regular_gaussian_pressure_level.grib1') =>
+  {"t"=>{:dims=>["lon","lat"], :shape=>[128,64], :val=>{[0,0]=>2.4746011353e2}}
+  },
+  File.join(dir, 'regular_gaussian_pressure_level.grib2') =>
+  {"t"=>{:dims=>["lon","lat"], :shape=>[128,64], :val=>{[0,0]=>2.4746011353e2}}
+  },
+  File.join(dir, 'regular_gaussian_pressure_level_constant.grib1') =>
+  {"t"=>{:dims=>["lon","lat"], :shape=>[128,64], :val=>{[0,0]=>1.0}}
+  },
+  File.join(dir, 'regular_gaussian_pressure_level_constant.grib2') =>
+  {"t"=>{:dims=>["lon","lat"], :shape=>[128,64], :val=>{[0,0]=>1.0}}
+  },
+  File.join(dir, 'regular_latlon_surface.grib1') =>
+  {"2t"=>{:dims=>["lon","lat"], :shape=>[16,31], :val=>{[0,0]=>279.0}}
+  },
+  File.join(dir, 'regular_latlon_surface.grib2') =>
+  {"2t"=>{:dims=>["lon","lat"], :shape=>[16,31], :val=>{[0,0]=>279.0}}
+  },
+  File.join(dir, 'regular_latlon_surface_constant.grib1') =>
+  {"2t"=>{:dims=>["lon","lat"], :shape=>[16,31], :val=>{[0,0]=>1.0}}
+  },
+  File.join(dir, 'regular_latlon_surface_constant.grib2') =>
+  {"2t"=>{:dims=>["lon","lat"], :shape=>[16,31], :val=>{[0,0]=>1.0}}
+  },
+  File.join(dir, 'tp_ecmwf.grib') =>
+  {"tp"=>{:dims=>["lon","lat","step"], :shape=>[16,31,4], :val=>{[0,0,0]=>2.0546913147e-3}}
+  }
+}
 fname_not_exist = 'not_exist'
-fnames = [
-          fname_reg_gaus_surf_g1,
-          fname_reg_gaus_surf_g2,
-          fname_reg_gaus_ml_g1,
-          fname_reg_gaus_ml_g2,
-          fname_reg_gaus_pl_g1,
-          fname_reg_gaus_pl_g2,
-          fname_reg_gaus_pl_const_g1,
-          fname_reg_gaus_pl_const_g2,
-          fname_reg_ll_surf_g1,
-          fname_reg_ll_surf_g2,
-          fname_reg_ll_surf_const_g1
-         ]
 
 describe NumRu::Grib do
   describe 'Grib.is_a_Grib?' do
     it 'returns whether the specific file is a GRIB1/2 file' do
-      NumRu::Grib.is_a_Grib?(fname_reg_gaus_surf_g1).should be_true
-      NumRu::Grib.is_a_Grib?(fname_reg_gaus_surf_g2).should be_true
+      GribData.each do |fname, hash|
+        NumRu::Grib.is_a_Grib?(fname).should be_true
+      end
     end
   end
 
   describe 'Grib.open' do
-    context 'specified a file name' do
-      describe 'grib1' do
-        it 'returns Grib object' do
-          (@file = NumRu::Grib.open(fname_reg_gaus_surf_g1)).should be_a_kind_of(NumRu::Grib)
-        end
-        after { @file.close if @file }
-      end
-      describe 'grib2' do
-        it 'returns Grib object' do
-          (@file = NumRu::Grib.open(fname_reg_gaus_surf_g2)).should be_a_kind_of(NumRu::Grib)
-        end
-        after { @file.close if @file }
+    before { @files = [] }
+    after { @files.each{|file| file.close} }
+    it 'returns Grib object' do
+      GribData.each do |fname, hash|
+        file = NumRu::Grib.open(fname)
+        file.should be_a_kind_of(NumRu::Grib)
+        @files.push file
       end
     end
-    context 'specified a file name which does not exist' do
-      it { proc { NumRu::Grib.open(fname_not_exist)}.should raise_error }
+    it 'raise exception when file name which does not exist is given' do
+      proc { NumRu::Grib.open(fname_not_exist)}.should raise_error
     end
   end
 
   describe '#close' do
-    file = NumRu::Grib.open(fname_reg_gaus_surf_g1)
-    it { file.close }
+    before do
+      @files = []
+      GribData.each{|fname,hash| @files.push NumRu::Grib.open(fname) }
+    end
+    it do
+      @files.each do |file|
+        proc {file.close}.should_not raise_error
+      end
+    end
   end
 
   describe 'instance methods' do
     before do
-      @files = fnames.map{|fname| NumRu::Grib.open(fname)}
+      @files = {}
+      GribData.each{|fname,hash| @files[fname] = NumRu::Grib.open(fname) }
     end
     after do
-      @files.each{|file| file.close}
+      @files.each{|fname,file| file.close}
     end
     it '#path returns path' do
-      @files.each {|file| file.path.should be_a_kind_of(String)}
+      @files.each {|fname,file| file.path.should == fname}
     end
     it '#var_names returns Array of variable names' do
-      @files.each {|file|  file.var_names.should be_a_kind_of(Array)}
+      @files.each {|fname, file| file.var_names.should == GribData[fname].keys}
     end
     it '#var returns GribVar object' do
-      @files.each do |file|
-        vname = file.var_names[0]
-        file.var(vname).should be_a_kind_of NumRu::GribVar
+      @files.each do |fname, file|
+        GribData[fname].each do |vname, hash|
+          file.var(vname).should be_a_kind_of NumRu::GribVar
+        end
       end
     end
   end
@@ -85,63 +106,68 @@ end
 
 describe NumRu::GribVar do
   before do
-    @vars = []
     @files = []
-    fnames.each do |fname|
+    @vars = Hash.new
+    GribData.each do |fname,vars|
       file = NumRu::Grib.open(fname)
       @files.push file
-      file.var_names.each{|vname| @vars.push file.var(vname)}
+      vars.each{|vname, hash| @vars[file.var(vname)] = hash}
     end
   end
-  after do
-    @files.each{|file| file.close}
-  end
+  after { @files.each{|file| file.close}}
   it '#rank returns rank' do
-    @vars.each{|var| var.rank.should be_a_kind_of Fixnum}
+    @vars.each{|var,hash| var.rank.should == hash[:shape].length}
   end
   it '#total returns total length' do
-    @vars.each{|var| var.total.should be_a_kind_of Fixnum}
+    @vars.each{|var,hash| var.total.should == hash[:shape].inject(1,:*)}
   end
   it '#dim_names returns Array of dimension names' do
-    @vars.each{|var| var.dim_names.should be_a_kind_of Array}
+    @vars.each{|var,hash| var.dim_names.should == hash[:dims]}
   end
   it '#shape returns Array of shape' do
-    @vars.each{|var| var.shape.should be_a_kind_of Array}
+    @vars.each{|var,hash| var.shape.should == hash[:shape]}
   end
   it '#dim returns GribDim' do
-    @vars.each do |var|
-      var.dim_names.each{|dname| var.dim(dname).should be_a_kind_of NumRu::GribDim}
+    @vars.each do |var,hash|
+      hash[:dims].each{|dname| var.dim(dname).should be_a_kind_of NumRu::GribDim}
     end
   end
   it '#att_names returns Array of attribute names' do
-    @vars.each{|var| var.att_names.should be_a_kind_of Array}
+    @vars.each{|var,hash| var.att_names.should be_a_kind_of Array}
   end
   it '#att returns attribute value' do
-    @vars.each do |var|
+    @vars.each do |var,hash|
       var.att_names.each{|aname| var.att(aname).should_not == nil }
     end
   end
   it '#typecode returns type code' do
-    @vars.each{|var| var.typecode.should be_a_kind_of Fixnum}
+    @vars.each{|var,hash| var.typecode.should be_a_kind_of Fixnum}
   end
   it '#missing_value returns missing value' do
-    @vars.each{|var| var.missing_value.should be_a_kind_of Numeric}
+    @vars.each{|var,hash| var.missing_value.should be_a_kind_of Numeric}
   end
   it '#get returns value' do
-    @vars.each{|var| var.get.should be_a_kind_of NArrayMiss}
+    @vars.each do |var,hash|
+      hash[:val].each do |idx,val|
+        v = var.get[*idx]
+        ((v-val).abs/((v+val)/2)).should < 1e-10
+      end
+    end
   end
 end
 
 describe NumRu::GribDim do
   before do
-    @dims = []
     @files = []
-    fnames.each do |fname|
+    @dims = Hash.new
+    GribData.each do |fname, vars|
       file = NumRu::Grib.open(fname)
       @files.push file
-      file.var_names.each do |vname|
+      vars.each do |vname, hash|
         var = file.var(vname)
-        var.dim_names{|dname| @dims.push var.dim(dname)}
+        hash[:dims].each_with_index do |dname,i|
+          @dims[var.dim(dname)] = hash[:shape][i]
+        end
       end
     end
   end
@@ -149,19 +175,22 @@ describe NumRu::GribDim do
     @files.each{|file| file.close}
   end
   it '#length returns length' do
-    @dims.each{|dim| dim.length.should be_a_kind_of Integer}
+    @dims.each{|dim,len| dim.length.should == len}
   end
   it '#val returns value' do
-    @dims.each{|dim| dim.val.should be_a_kind_of NArray}
+    @dims.each do |dim,len|
+      dim.val.should be_a_kind_of NArray
+      dim.val.length.should == len
+    end
   end
   it '#typecode returns typecode' do
-    @dims.each{|dim| dim.typecode.should be_a_kind_of Fixnum}
+    @dims.each{|dim,len| dim.typecode.should be_a_kind_of Fixnum}
   end
   it '#att_names returns Array of attribute names' do
-    @dims.each{|dim| dim.att_names.should be_a_kind_of Array}
+    @dims.each{|dim,len| dim.att_names.should be_a_kind_of Array}
   end
   it '#att returns attributes value' do
-    @dims.each do |dim|
+    @dims.each do |dim,len|
       dim.att_names.each{|aname| dim.att(aname).should_not == nil }
     end
   end

metadata

@@ -1,13 +1,13 @@
 --- !ruby/object:Gem::Specification 
 name: rb-grib
 version: !ruby/object:Gem::Version 
-  hash: 19
-  prerelease: 
+  hash: 23
+  prerelease: false
   segments: 
   - 0
-  - 1
-  - 4
-  version: 0.1.4
+  - 2
+  - 0
+  version: 0.2.0
 platform: ruby
 authors: 
 - Seiya Nishizawa
@@ -15,7 +15,8 @@ autorequire:
 bindir: bin
 cert_chain: []
 
-date: 2011-11-29 00:00:00 Z
+date: 2011-12-15 00:00:00 +09:00
+default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
   name: rspec
@@ -31,6 +32,34 @@ dependencies:
         version: "0"
   type: :development
   version_requirements: *id001
+- !ruby/object:Gem::Dependency 
+  name: narray
+  prerelease: false
+  requirement: &id002 !ruby/object:Gem::Requirement 
+    none: false
+    requirements: 
+    - - ">="
+      - !ruby/object:Gem::Version 
+        hash: 3
+        segments: 
+        - 0
+        version: "0"
+  type: :runtime
+  version_requirements: *id002
+- !ruby/object:Gem::Dependency 
+  name: narray_miss
+  prerelease: false
+  requirement: &id003 !ruby/object:Gem::Requirement 
+    none: false
+    requirements: 
+    - - ">="
+      - !ruby/object:Gem::Version 
+        hash: 3
+        segments: 
+        - 0
+        version: "0"
+  type: :runtime
+  version_requirements: *id003
 description: This class library enable you to handle GRIB file.
 email: 
 - seiya@gfd-dennou.org
@@ -42,6 +71,7 @@ extra_rdoc_files: []
 
 files: 
 - .gitignore
+- BSDL
 - Gemfile
 - LICENSE.txt
 - README.rdoc
@@ -59,14 +89,25 @@ files:
 - data/regular_latlon_surface_constant.grib1
 - data/regular_latlon_surface_constant.grib2
 - data/regular_latlon_surface_prec.grib1
+- data/tp_ecmwf.grib
 - ext/extconf.rb
 - ext/grib.c
 - lib/numru/grib.rb
+- lib/numru/grib/definitions/grib1/localConcepts/kwbc/name.def
+- lib/numru/grib/definitions/grib1/localConcepts/kwbc/paramId.def
+- lib/numru/grib/definitions/grib1/localConcepts/kwbc/shortName.def
+- lib/numru/grib/definitions/grib1/localConcepts/kwbc/units.def
+- lib/numru/grib/definitions/grib2/localConcepts/rjtd/name.def
+- lib/numru/grib/definitions/grib2/localConcepts/rjtd/paramId.def
+- lib/numru/grib/definitions/grib2/localConcepts/rjtd/shortName.def
+- lib/numru/grib/definitions/grib2/localConcepts/rjtd/units.def
 - lib/numru/grib/grib.rb
+- lib/numru/grib/setenv.rb
 - lib/numru/grib/version.rb
 - rb-grib.gemspec
 - spec/grib_read_spec.rb
 - spec/spec_helper.rb
+has_rdoc: true
 homepage: http://ruby.gfd-dennou.org/products/rb-grib/
 licenses: []
 
@@ -96,7 +137,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
 requirements: []
 
 rubyforge_project: rb-grib
-rubygems_version: 1.8.10
+rubygems_version: 1.3.7
 signing_key: 
 specification_version: 3
 summary: Ruby class library to hanlde GRIB file