cdo 1.3.4 → 1.3.6

checksums.yaml

@@ -1,7 +1,7 @@
 ---
-SHA1:
-  metadata.gz: 839669f481246fb63e5ad6b7871127ef6b1c710d
-  data.tar.gz: fae0e2fa60efa590697c9cf73f0e1fc3db103cd8
+SHA256:
+  metadata.gz: ae07b4039ff9ecbaea1bdea8a4d50c83b63bdd5d65c6b87f5d0c8f67632dc937
+  data.tar.gz: 23d76998b8cd25d1b0ba6619651c255c33920e104181912ccdbe588bfb580e08
 SHA512:
-  metadata.gz: 12c44ad5bb29e96457a623d07060a87557f79b2ee9246faf1b0e477865ff5603219118260fa2db0623a149dda45f133db132980ddd18c41c44949dfb86c06fa0
-  data.tar.gz: '02851b39cb5c9fd19941e0e85c594d01fdcd51349bf0df6255477fb92ae1ede816a42d438fe19ba077ed7dc1958bcfc38d789031186faf09361b76c43c97ddd3'
+  metadata.gz: d870e95cb67677c2f00aa707ce458ee314c7767c857c9085c65a7e838eba6b5e00e40c2c529b1c603b9444562de3e85f2cae812386516261460ff842660aa50e
+  data.tar.gz: 70d7277e46d3c12175fe9168a86aa60e60f0d8837a3d36a94906ce9c3ac725d7c2f8428fc8e841f8d5fbba1c92f7fd4f5bf2e2ba6bbbcda6935815507ac21f42

data/gemspec

@@ -3,9 +3,9 @@ $:.unshift File.join(File.dirname(__FILE__),"..","lib")
 
 spec = Gem::Specification.new do |s|
   s.name              = "cdo"
-  s.version           = '1.3.4'
+  s.version           = '1.3.6'
   s.platform          = Gem::Platform::RUBY
-  s.files             = ["lib/cdo.rb","lib/cdo_lib.rb"] + ["gemspec","LICENSE"]
+  s.files             = ["lib/cdo.rb","gemspec","LICENSE"]
   s.test_file         = "test/test_cdo.rb"
   s.description       = "Easy access to the Climate Data operators"
   s.summary           = "Easy access to the Climate Data operators"
@@ -16,6 +16,8 @@ spec = Gem::Specification.new do |s|
   s.required_ruby_version = ">= 1.9"
   s.add_development_dependency('unifiedPlot')
   s.add_development_dependency('minitest')
+  s.add_development_dependency('facets')
+  s.add_development_dependency('colorize')
   s.add_development_dependency('rake')
 end
 

data/lib/cdo.rb

@@ -3,6 +3,11 @@ require 'open3'
 require 'logger'
 require 'stringio'
 
+class Hash
+  alias :include? :has_key?
+end
+
+
 class Cdo
 
   # hardcoded fallback list of output operators - from 1.8.0 there is an
@@ -69,7 +74,8 @@ class Cdo
 
   # collect the complete list of possible operators
   def getOperators(path2cdo)
-    if version <= '1.7.0' then
+    case
+    when version <= '1.7.0' then
       cmd       = path2cdo + ' 2>&1'
       help      = IO.popen(cmd).readlines.map {|l| l.chomp.lstrip}
       if 5 >= help.size
@@ -78,12 +84,21 @@ class Cdo
         exit
       end
 
-      @operators = help[(help.index("Operators:")+1)..help.index(help.find {|v| v =~ /CDO version/ }) - 2].join(' ').split
+      operators = help[(help.index("Operators:")+1)..help.index(help.find {|v| v =~ /CDO version/ }) - 2].join(' ').split
+    when version <= '1.9.1' then
+      cmd       = "#{path2cdo} --operators"
+      operators = IO.popen(cmd).readlines.map {|l| l.split(' ').first }
     else
-      cmd = "#{path2cdo} --operators"
-
-      @operators = IO.popen(cmd).readlines.map {|l| l.split(' ').first }
+      cmd       = "#{path2cdo} --operators"
+      operators = {}
+      IO.popen(cmd).readlines.map {|line|
+        lineContent = line.chomp.split(' ')
+        name = lineContent[0]
+        iCounter, oCounter = lineContent[-1].tr(')','').tr('(','').split('|')
+        operators[name] = { iStreams: iCounter.to_i, oStreams: oCounter.to_i }
+      }
     end
+    return operators
   end
 
   def getNoOuputOperators(path2cdo)
@@ -285,9 +300,15 @@ class Cdo
     puts collectLogs
   end
 
+  def hasCdo(path=@cdo)
+    executable = system("#{path} -V >/dev/null 2>&1")
+    fullpath   = File.exists?(path) and File.executable?(path)
+
+    return (executable or fullpath)
+  end
   # check if cdo backend is working
   def check
-    return false unless system("#@cdo -h 1>/dev/null 2>&1")
+    return false unless hasCdo
 
     retval = _call("#{@cdo} -V")
     pp retval if @debug

data/test/test_cdo.rb

@@ -1,8 +1,7 @@
 $:.unshift File.join(File.dirname(__FILE__),"..","lib")
+require 'cdo'
 
 require 'minitest/autorun'
-require 'cdo'
-require 'pp'
 
 
 #===============================================================================
@@ -15,6 +14,7 @@ class TestCdo < Minitest::Test
 
   @@show           = ENV.has_key?('SHOW')
   @@maintainermode = ENV.has_key?('MAINTAINERMODE')
+  @@debug          = ENV.has_key?('DEBUG')
 
   def setup
     @cdo = Cdo.new
@@ -25,6 +25,8 @@ class TestCdo < Minitest::Test
     if ENV['CDO']
       assert_equal(ENV['CDO'],@cdo.cdo)
     else
+      pp DEFAULT_CDO_PATH
+      pp @cdo.cdo
       assert_equal(DEFAULT_CDO_PATH,@cdo.cdo)
     end
     newCDO="#{ENV['HOME']}/local/bin/cdo-dev"
@@ -41,6 +43,14 @@ class TestCdo < Minitest::Test
   def test_V
     puts @cdo.version(verbose=true)
   end
+  def test_hasCdo
+    assert(@cdo.hasCdo)
+    @cdo.cdo = 'cccccccc'
+    assert_equal( false, @cdo.hasCdo)
+
+    @cdo.cdo='/bin/cdo'
+    assert(@cdo.hasCdo) if File.exist?(@cdo.cdo)
+  end
   def test_getOperators
     %w[for random stdatm info showlevel sinfo remap geopotheight mask topo thicknessOfLevels].each {|op|
       if ["thicknessOfLevels"].include?(op)
@@ -56,7 +66,7 @@ class TestCdo < Minitest::Test
   end
 
   def test_outputOperators
-    @cdo.debug = true
+    @cdo.debug = @@debug
     levels = @cdo.showlevel(:input => "-stdatm,0")
     assert_equal([0,0].map(&:to_s),levels)
 
@@ -72,7 +82,7 @@ class TestCdo < Minitest::Test
   end
   def test_CDO_version
     assert("1.4.3.1" < @cdo.version,"Version too low: #{@cdo.version}")
-    assert("1.8.1" < @cdo.version,"Version too low: #{@cdo.version}")
+    assert("1.7.1" < @cdo.version,"Version too low: #{@cdo.version}")
     assert("3.0" > @cdo.version,"Version too high: #{@cdo.version}")
   end
   def test_args
@@ -84,7 +94,7 @@ class TestCdo < Minitest::Test
     (1...info.size).each {|i| assert_equal(0.0,info[i].split[-1].to_f)}
   end
   def test_operator_options
-    @cdo.debug=true
+    @cdo.debug=@@debug
     targetLevels = [0,10,50,100,200,400,1000]
     levels = @cdo.showlevel(:input => " -stdatm,#{targetLevels.join(',')}")
     [0,1].each {|i| assert_equal(targetLevels.join(' '),levels[i])}
@@ -92,7 +102,7 @@ class TestCdo < Minitest::Test
     assert_equal(["P T"],names)
   end
   def test_chain
-    @cdo.debug = true
+    @cdo.debug = @@debug
     ofile = @cdo.setname('veloc',:input => " -copy -random,r1x1",:options => "-f nc")
     assert_equal(["veloc"],@cdo.showname(:input => ofile))
   end
@@ -127,6 +137,7 @@ class TestCdo < Minitest::Test
     outs.clear
     outs << @cdo.stdatm(0,10,20,:output => 'test_force')
     mtime0 = File.stat(outs[-1]).mtime
+    sleep(1)
     outs << @cdo.stdatm(0,10,20,:output => 'test_force')
     mtime1 = File.stat(outs[-1]).mtime
     assert(mtime0 != mtime1)
@@ -194,7 +205,7 @@ class TestCdo < Minitest::Test
   end
 
   def test_verticalLevels
-    @cdo.debug = true
+    @cdo.debug = @@debug
     targetThicknesses = [50.0,  100.0,  200.0,  300.0,  450.0,  600.0,  800.0, 1000.0, 1000.0, 1000.0]
     sourceLevels = %W{25 100 250 500 875 1400 2100 3000 4000 5000}
     thicknesses = @cdo.thicknessOfLevels(:input => "-selname,T #{@cdo.stdatm(*sourceLevels,:options => '-f nc')}")
@@ -212,7 +223,7 @@ class TestCdo < Minitest::Test
   end
 
   def test_errorException
-    @cdo.debug = true
+    @cdo.debug = @@debug
     # stdout operators get get wrong input
     assert_raises ArgumentError do
       @cdo.showname(:input => '-for,d')
@@ -303,6 +314,11 @@ class TestCdo < Minitest::Test
   if @@maintainermode  then
     require 'unifiedPlot'
 
+    def test_longChain
+      ifile = "-enlarge,global_0.3 -settaxis,2000-01-01 -expr,'t=sin(for*3.141529/180.0)' -for,1,10"
+      t = @cdo.fldmax(input: "-div -sub -timmean -seltimestep,2,3 #{ifile} -seltimestep,1 #{ifile}  -gridarea #{ifile}",returnArray: "t")
+      assert_equal(8.981299259858133e-09,t[0])
+    end
     def test_returnArray
       temperature = @cdo.stdatm(0,:options => '-f nc',:returnCdf => true).var('T').get.flatten[0]
       assert_raises ArgumentError do
@@ -314,7 +330,7 @@ class TestCdo < Minitest::Test
       assert_equal("1013.25 898.543456035875",pressure.flatten.to_a.join(' '))
     end
     def test_returnMaArray
-      @cdo.debug = true
+      @cdo.debug = @@debug
       topo = @cdo.topo(:options => '-f nc',:returnMaArray => 'topo')
       assert_equal(-1890.0,topo.mean.round)
       bathy = @cdo.setrtomiss(0,10000,
@@ -382,17 +398,17 @@ class TestCdo < Minitest::Test
       input = "-settunits,days  -setyear,2000 -for,1,4"
       cdfFile = @cdo.copy(:options =>"-f nc",:input=>input)
       cdf     = @cdo.readCdf(cdfFile)
-      assert_equal(['lon','lat','time','for'],cdf.var_names)
+      assert_empty(['lon','lat','for','time'] - cdf.var_names)
     end
     def test_selIndexListFromIcon
       input = "~/data/icon/oce.nc"
     end
     def test_readArray
-      @cdo.debug = true
+      @cdo.debug = @@debug
       assert_equal([40,80],@cdo.readArray(@cdo.sellonlatbox(-10,10,-20,20,:input => '-topo',:options => '-f nc'), 'topo').shape)
     end
     def test_doc
-      @cdo.debug = true
+      @cdo.debug = @@debug
       @cdo.help(:remap)
       @cdo.help(:infov)
       @cdo.help(:topo)
@@ -400,7 +416,7 @@ class TestCdo < Minitest::Test
       @cdo.help
     end
     def test_fillmiss
-      @cdo.debug = true
+      @cdo.debug = @@debug
       # check up-down replacement
       rand = @cdo.setname('v',:input => '-random,r1x10 ', :options => ' -f nc',:output => '/tmp/rand.nc')
       cdf  = @cdo.openCdf(rand)

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: cdo
 version: !ruby/object:Gem::Version
-  version: 1.3.4
+  version: 1.3.6
 platform: ruby
 authors:
 - Ralf Mueller
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2017-09-06 00:00:00.000000000 Z
+date: 2018-04-03 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: unifiedPlot
@@ -38,6 +38,34 @@ dependencies:
     - - ">="
       - !ruby/object:Gem::Version
         version: '0'
+- !ruby/object:Gem::Dependency
+  name: facets
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: colorize
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
 - !ruby/object:Gem::Dependency
   name: rake
   requirement: !ruby/object:Gem::Requirement
@@ -61,7 +89,6 @@ files:
 - LICENSE
 - gemspec
 - lib/cdo.rb
-- lib/cdo_lib.rb
 - test/test_cdo.rb
 homepage: https://code.zmaw.de/projects/cdo/wiki/Cdo%7Brbpy%7D
 licenses:
@@ -83,7 +110,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 2.6.11
+rubygems_version: 2.7.4
 signing_key: 
 specification_version: 4
 summary: Easy access to the Climate Data operators

data/lib/cdo_lib.rb

@@ -1,415 +0,0 @@
-require 'pp'
-require 'open3'
-require 'logger'
-require 'stringio'
-
-# Copyright (C) 2011-2013 Ralf Mueller, ralf.mueller@zmaw.de
-# See COPYING file for copying and redistribution conditions.
-# 
-# This program is free software; you can redistribute it and/or modify
-# it under the terms of the GNU General Public License as published by
-# the Free Software Foundation; version 2 of the License.
-# 
-# This program is distributed in the hope that it will be useful,
-# but WITHOUT ANY WARRANTY; without even the implied warranty of
-# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
-# GNU General Public License for more details.
-
-# ==============================================================================
-# CDO calling mechnism
-module Cdo
-
-  VERSION = "1.2.5"
-  @@file  = StringIO.new
-
-  State = {
-    :debug       => false,
-    :returnCdf   => false,
-    :operators   => [],
-    :forceOutput => true,
-    :env         => {},
-    :log         => false,
-    :logger      => Logger.new(@@file),
-  }
-  State[:debug] = ENV.has_key?('DEBUG')
-  State[:logger].formatter = proc do |serverity, time, progname, msg|
-    msg
-  end
-
-  @@CDO = ENV['CDO'].nil? ? 'cdo' : ENV['CDO']
-
-  # Since cdo-1.5.4 undocumented operators are given with the -h option. For
-  # earlier version, they have to be provided manually
-  @@undocumentedOperators = %w[anomaly beta boxavg change_e5lsm change_e5mask
-    change_e5slm chisquare chvar cloudlayer cmd com command complextorect
-    covar0 covar0r daycount daylogs del29feb delday delete deltap deltap_fl
-    delvar diffv divcoslat dumplogo dumplogs duplicate eca_r1mm enlargegrid
-    ensrkhistspace ensrkhisttime eof3d eof3dspatial eof3dtime export_e5ml
-    export_e5res fc2gp fc2sp fillmiss fisher fldcovar fldrms fourier fpressure
-    gather gengrid geopotheight ggstat ggstats globavg gp2fc gradsdes
-    gridverify harmonic hourcount hpressure import_e5ml import_e5res
-    import_obs imtocomplex infos infov interpolate intgrid intgridbil
-    intgridtraj intpoint isosurface lmavg lmean lmmean lmstd log lsmean
-    meandiff2test mergegrid mod moncount monlogs mrotuv mrotuvb mulcoslat ncode
-    ncopy nmltest normal nvar outputbounds outputboundscpt outputcenter
-    outputcenter2 outputcentercpt outputkey outputtri outputvector outputvrml
-    pardup parmul pinfo pinfov pressure_fl pressure_hl read_e5ml remapcon1
-    remapdis1 retocomplex scalllogo scatter seascount select selgridname
-    seloperator selvar selzaxisname setrcaname setvar showvar sinfov smemlogo
-    snamelogo sort sortcode sortlevel sortname sorttaxis sorttimestamp sortvar
-    sp2fc specinfo spectrum sperclogo splitvar stimelogo studentt template1
-    template2 test test2 testdata thinout timcount timcovar tinfo transxy trms
-    tstepcount vardes vardup varmul varquot2test varrms vertwind write_e5ml
-    writegrid writerandom yearcount]
-
-  @@outputOperatorsPattern = /(diff|info|output|griddes|zaxisdes|show|ncode|ndate|nlevel|nmon|nvar|nyear|ntime|npar|gradsdes|pardes)/
-
-  private
-  def Cdo.getOperators(force=false)
-    # Do NOT compute anything, if it is not required
-    return State[:operators] unless (State[:operators].empty? or force)
-    cmd       = @@CDO + ' 2>&1'
-    help      = IO.popen(cmd).readlines.map {|l| l.chomp.lstrip}
-    if 5 >= help.size
-      warn "Operators could not get listed by running the CDO binary (#{@@CDO})"
-      pp help if Cdo.debug
-      exit
-    end
-    # in version 1.5.6 the output of '-h' has changed
-    State[:operators] = case 
-                        when Cdo.version < "1.5.6"
-                          (help[help.index("Operators:")+1].split + @@undocumentedOperators).uniq
-                        when Cdo.version <= "1.7.0"
-                          help[(help.index("Operators:")+1)..help.index(help.find {|v| v =~ /CDO version/}) - 2].join(' ').split
-                        else
-                          IO.popen(@@CDO + ' --operators').readlines.map {|l| l.split(' ').first}
-                        end
-  end
-
-  def Cdo.hasError(cmd,retvals)
-    if (State[:debug])
-      puts("RETURNCODE: #{retvals[:returncode]}")
-    end
-    if ( 0 != retvals[:returncode] )
-      puts("Error in calling:")
-      puts(">>> "+cmd+"<<<")
-      puts(retvals[:stderr])
-      return true
-    else
-      return false
-    end
-  end
-
-  def Cdo.env=(envHash)
-    State[:env] = envHash
-  end
-  def Cdo.env; State[:env]; end
-
-  def Cdo.call(cmd)
-    if (State[:debug])
-      puts '# DEBUG ====================================================================='
-      pp Cdo.env unless Cdo.env.empty?
-      puts 'CMD: '
-      puts cmd
-      puts '# DEBUG ====================================================================='
-    end
-    stdin, stdout, stderr, wait_thr = Open3.popen3(Cdo.env,cmd)
-
-    {
-      :stdout     => stdout.read,
-      :stderr     => stderr.read,
-      :returncode => wait_thr.value.exitstatus
-    }
-  end
-
-  def Cdo.run(cmd,ofile='',options='',returnCdf=false,force=nil,returnArray=nil,returnMaArray=nil)
-    cmd = "#{@@CDO} -O #{options} #{cmd} "
-    case ofile
-    when $stdout
-      retvals = Cdo.call(cmd)
-      State[:logger].info(cmd+"\n") if State[:log]
-      unless hasError(cmd,retvals)
-        return retvals[:stdout].split($/).map {|l| l.chomp.strip}
-      else
-        raise ArgumentError,"CDO did NOT run successfully!"
-      end
-    else
-      force = State[:forceOutput] if force.nil?
-      if force or not File.exists?(ofile.to_s)
-        ofile = MyTempfile.path if ofile.nil?
-        cmd << "#{ofile}"
-        retvals = call(cmd)
-        State[:logger].info(cmd+"\n") if State[:log]
-        if hasError(cmd,retvals)
-          raise ArgumentError,"CDO did NOT run successfully!"
-        end
-      else
-        warn "Use existing file '#{ofile}'" if Cdo.debug
-      end
-    end
-    if not returnArray.nil?
-      Cdo.readArray(ofile,returnArray)
-    elsif not returnMaArray.nil?
-      Cdo.readMaArray(ofile,returnMaArray)
-    elsif returnCdf or State[:returnCdf]
-      Cdo.readCdf(ofile)
-    else
-      return ofile
-    end
-  end
-
-  def Cdo.parseArgs(args)
-    # splitinto hash-like args and the rest
-    operatorArgs = args.reject {|a| a.class == Hash}
-    opts = operatorArgs.empty? ? '' : ',' + operatorArgs.join(',')
-    io   = args.find {|a| a.class == Hash}
-    io   = {} if io.nil?
-    args.delete_if   {|a| a.class == Hash}
-    # join input streams together if possible
-    io[:input] = io[:input].join(' ') if io[:input].respond_to?(:join)
-    return [io,opts]
-  end
-
-  def Cdo.method_missing(sym, *args, &block)
-    ## args is expected to look like [opt1,...,optN,:input => iStream,:output => oStream] where
-    # iStream could be another CDO call (timmax(selname(Temp,U,V,ifile.nc))
-    puts "Operator #{sym.to_s} is called" if State[:debug]
-    if getOperators.include?(sym.to_s)
-      io, opts = Cdo.parseArgs(args)
-      if @@outputOperatorsPattern.match(sym)
-        run(" -#{sym.to_s}#{opts} #{io[:input]} ",$stdout)
-      else
-        run(" -#{sym.to_s}#{opts} #{io[:input]} ",io[:output],io[:options],io[:returnCdf],io[:force],io[:returnArray],io[:returnMaArray])
-      end
-    else
-      raise ArgumentError,"Operator #{sym.to_s} not found"
-    end
-  end
-
-  def Cdo.loadCdf
-    begin
-      require "numru/netcdf_miss"
-      include NumRu
-    rescue LoadError
-      warn "Could not load ruby's netcdf bindings. Please install it."
-      raise
-    end
-  end
-
-  def Cdo.getSupportedLibs(force=false)
-    return unless (State[:libs].nil? or force)
-    _, _, stderr, _ = Open3.popen3(@@CDO + " -V")
-    supported       = stderr.readlines.map(&:chomp)
-    with            = supported.grep(/(with|Features)/)[0].split(':')[1].split.map(&:downcase)
-    libs            = supported.grep(/library version/).map {|l| 
-      l.strip.split(':').map {|l| 
-        l.split.first.downcase
-      }[0,2]
-    }
-    State[:libs] = {}
-    with.flatten.each {|k| State[:libs][k]=true}
-    libs.each {|lib,version| State[:libs][lib] = version}
-  end
-
-  public
-  def Cdo.debug=(value)
-    State[:debug] = value
-  end
-  def Cdo.debug
-    State[:debug]
-  end
-  def Cdo.forceOutput=(value)
-    State[:forceOutput] = value
-  end
-  def Cdo.forceOutput
-    State[:forceOutput]
-  end
-  def Cdo.log=(value)
-    State[:log] = value
-  end
-
-  def Cdo.log
-    State[:log]
-  end
-
-  def Cdo.version
-    cmd     = @@CDO + ' 2>&1'
-    help    = IO.popen(cmd).readlines.map {|l| l.chomp.lstrip}
-    regexp  = %r{CDO version (\d.*), Copyright}
-    line    = help.find {|v| v =~ regexp}
-    version = regexp.match(line)[1]
-  end
-
-  def Cdo.setReturnCdf(value=true)
-    if value
-      Cdo.loadCdf
-    end
-    State[:returnCdf] = value
-  end
-
-  def Cdo.unsetReturnCdf
-    setReturnCdf(false)
-  end
-
-  def Cdo.returnCdf
-    State[:returnCdf]
-  end
-
-  def Cdo.hasCdo?(bin=@@CDO)
-    return true if File.exists?(@@CDO) and File.executable?(@@CDO)
-    ENV['PATH'].split(File::PATH_SEPARATOR).each {|path| 
-      return true if File.exists?([path,bin].join(File::SEPARATOR))
-    }
-  end
-
-  # test if @@CDO can be used
-  def Cdo.checkCdo
-    unless hasCdo?(@@CDO)
-      warn "Testing application #@@CDO is not available!"
-      exit 1
-    else
-      puts "Using CDO: #@@CDO"
-      puts IO.popen(@@CDO + " -V").readlines
-    end
-    return true
-  end
-
-  def Cdo.setCdo(cdo)
-    puts "Will use #{cdo} instead of #@@CDO" if Cdo.debug
-    @@CDO = cdo
-    Cdo.getOperators(true)
-    Cdo.getSupportedLibs(true)
-  end
-
-  def Cdo.getCdo
-    @@CDO
-  end
-
-  def Cdo.operators
-    Cdo.getOperators if State[:operators].empty?
-    State[:operators]
-  end
-
-  def Cdo.libs
-    getSupportedLibs
-    State[:libs]
-  end
-
-  def Cdo.hasLib?(lib)
-    return Cdo.libs.has_key?(lib)
-    return false
-  end
-
-  def Cdo.libsVersion(lib)
-    unless Cdo.hasLib?(lib)
-      raise ArgumentError, "Cdo does NOT have support for '#{lib}'"
-    else
-      if State[:libs][lib].kind_of? String
-        return State[:libs][lib]
-      else
-        warn "No version information available about '#{lib}'"
-        return false
-      end
-    end
-  end
-
-  def Cdo.showlog
-    @@file.rewind
-    puts @@file.read
-  end
-
-  #==================================================================
-  # Addional operotors:
-  #------------------------------------------------------------------
-  def Cdo.boundaryLevels(args)
-    ilevels         = Cdo.showlevel(:input => args[:input])[0].split.map(&:to_f)
-    bound_levels    = Array.new(ilevels.size+1)
-    bound_levels[0] = 0
-    (1..ilevels.size).each {|i| 
-      bound_levels[i] =bound_levels[i-1] + 2*(ilevels[i-1]-bound_levels[i-1])
-    }
-    bound_levels
-  end
-
-  def Cdo.thicknessOfLevels(args)
-    bound_levels = Cdo.boundaryLevels(args)
-    delta_levels    = []
-    bound_levels.each_with_index {|v,i| 
-      next if 0 == i
-      delta_levels << v - bound_levels[i-1]
-    }
-    delta_levels
-  end
-
-  def Cdo.readCdf(iFile)
-    Cdo.loadCdf unless State[:returnCdf] 
-    NetCDF.open(iFile)
-  end
-
-  def Cdo.openCdf(iFile)
-    Cdo.loadCdf unless State[:returnCdf] 
-    NetCDF.open(iFile,'r+')
-  end
-
-  def Cdo.readArray(iFile,varname)
-    filehandle = Cdo.readCdf(iFile)
-    if filehandle.var_names.include?(varname)
-      # return the data array
-      filehandle.var(varname).get
-    else
-      raise ArgumentError, "Cannot find variable '#{varname}'"
-    end
-  end
-
-  def Cdo.readMaArray(iFile,varname)
-    filehandle = Cdo.readCdf(iFile)
-    if filehandle.var_names.include?(varname)
-      # return the data array
-      filehandle.var(varname).get_with_miss
-    else
-      raise ArgumentError,"Cannot find variable '#{varname}'"
-    end
-  end
-
-  def Cdo.help(operator=nil)
-    if operator.nil?
-      puts Cdo.call([@@CDO,'-h'].join(' '))[:stderr]
-    else
-      operator = operator.to_s unless String == operator.class
-      if Cdo.operators.include?(operator)
-        puts Cdo.call([@@CDO,'-h',operator].join(' '))[:stdout]
-      else
-        puts "Unknown operator #{operator}"
-      end
-    end
-  end
-end
-
-# Helper module for easy temp file handling
-module MyTempfile
-  require 'tempfile'
-  @@_tempfiles           = []
-  @@persistent_tempfiles = false
-  @@N                    = 10000000
-
-  def MyTempfile.setPersist(value)
-    @@persistent_tempfiles = value
-  end
-
-  def MyTempfile.path
-    unless @@persistent_tempfiles
-      t = Tempfile.new(self.class.to_s)
-      @@_tempfiles << t
-      @@_tempfiles << t.path
-      t.path
-    else
-      t = "_"+rand(@@N).to_s
-      @@_tempfiles << t
-      t
-    end
-  end
-
-  def MyTempfile.showFiles
-    @@_tempfiles.each {|f| print(f+" ") if f.kind_of? String}
-  end
-end