iconPlot 0.0.8 → 0.0.9

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 81daefe787ba86d077804f2afbb20f24bb8625e0
-  data.tar.gz: e2538b75ce60c906d8c90afd08654bc4a9f7a9f1
+  metadata.gz: df59cf8fcb883833c0355e39533c5da8e3075d25
+  data.tar.gz: 8346f47771d4ce6573732e64986e928a05eb9878
 SHA512:
-  metadata.gz: b2a7478b62d8e95335816b006caeb81471f749cf1d014875647a9504a2efe03bca72c5a7db1036650407ff9c5e0a6d59a04865cef19b869de7163bf1213ae5aa
-  data.tar.gz: c7da25dd5cd29dcc0232e049bf86da2192ab56cab2ff232a5b1a2f247d1158d8d1ace670f39c7a664607fe5382b6628a934867a4f0cb3ce9705b23558a2b5ef0
+  metadata.gz: 6b0086f6cb9c924a88b124b560298971805f810437da1bd3b0682f3339aca3bc7a36af54ae8d5f976df00710cc7aa951e50061b483de699b605dde68593187b9
+  data.tar.gz: 25de6154c6143f65914063ae23a1e4096e57fbae0e4c5636bed530f456c5f8e8206186e0b254de92a436d19feebcc16f30e6388232ae51a7f53cc0259a553c3f

data/gemspec

@@ -3,7 +3,7 @@ $:.unshift File.join(File.dirname(__FILE__),"..","lib")
 
 spec = Gem::Specification.new do |s|
   s.name              = "iconPlot"
-  s.version           = '0.0.8'
+  s.version           = '0.0.9'
   s.platform          = Gem::Platform::RUBY
   s.files             = ["lib/iconPlot.rb"] + ["gemspec"] + ["bin/nclsh"] + Dir.glob("lib/*ncl")
   s.executables       << 'nclsh'

data/lib/iconPlot.rb

@@ -35,7 +35,7 @@ class IconPlot < Struct.new(:caller,:plotter,:libdir,:otype,:display,:cdo,:debug
                  raise ArgumentError
                end
 
-    opts[:tStrg] =ofile if opts[:tStrg].nil?
+    opts[:tStrg] = opts[:tStrg].nil? ? ofile : opts[:tStrg]
 
     cmd   = [self.caller,self.plotter].join(' ')
     cmd << " -altLibDir=#{self.libdir} #{varIdent} -iFile=#{ifile} -oFile=#{ofile} -oType=#{self.otype}"
@@ -68,16 +68,16 @@ class IconPlot < Struct.new(:caller,:plotter,:libdir,:otype,:display,:cdo,:debug
       exit -1
     end
     Cdo.debug = true
-    unit = Cdo.showunit(:in => "-selname,#{varname} #{ifile}").first
+    unit = Cdo.showunit(:input => "-selname,#{varname} #{ifile}").first
     ExtCsvDiagram.plot_xy(icon,"datetime",varname,
                           "ICON: #{operation} on #{varname} (file:#{ifile})", # Change title here
-    :label_position => 'below',:skipColumnCheck => true,
-      :type => 'lines',:groupBy => ["depth"], :onlyGroupTitle => true,
+                        :label_position => 'below',:skipColumnCheck => true,
+                        :type => 'lines',:groupBy => ["depth"], :onlyGroupTitle => true,
       #                     :addSettings => ["logscale y"],     # Commend theses out for large scale values like Vert_Mixing_V
       #                     :yrange => '[0.0001:10]',           # Otherwise you'll see nothing reasonable
       :terminal => true ? 'png' : 'x11',
       :ylabel => "#{varname} [#{Shellwords.escape(unit)}]",     # Correct the label if necessary
-      :input_time_format => "'%Y%m%d %H:%M:%S'",
+      :input_time_format => "'%Y-%m-%d %H:%M:%S'",
         :filename => ofile,
         :output_time_format => '"%d.%m.%y \n %H:%M"',:size => "1600,600")
       return "#{ofile}.png"
@@ -85,7 +85,7 @@ class IconPlot < Struct.new(:caller,:plotter,:libdir,:otype,:display,:cdo,:debug
   def scatterPlot(ifile,ofile,xVarname,yVarname,opts={})
     # is there a variable which discribes different regions in the ocean
     regionVar = opts[:regionVar].nil? ? 'rregio_c' : opts[:regionVar]
-    hasRegion = Cdo.showname(:in => ifile).join.split.include?(regionName)
+    hasRegion = Cdo.showname(:input => ifile).join.split.include?(regionName)
     unless hasRegion
       warn "Variable '#{regionName}' for showing regions is NOT found in the input '#{ifile}'!"
       warn "Going on without plotting regions!"
@@ -100,7 +100,7 @@ class IconPlot < Struct.new(:caller,:plotter,:libdir,:otype,:display,:cdo,:debug
     FileUtils.rm(file) if File.exists?(file)
   end
   def show(*files)
-    files.flatten.each {|file| IO.popen("#{self.display} #{file} &") }
+    files.flatten.each {|file| out = IO.popen("#{self.display} #{file} &").read; puts out if self.debug }
   end
   def defaultPlot(ifile,ofile,opts={})
     show(scalarPlot(ifile,ofile,'T',opts))
@@ -117,13 +117,13 @@ class IconPlot < Struct.new(:caller,:plotter,:libdir,:otype,:display,:cdo,:debug
     IO.popen("echo 'date|time|depth|#{varname}' > #{dataFile}")
     Cdo.debug = true
     Cdo.outputkey('date,time,level,value', 
-                  :in => "-#{operation} -selname,#{varname} #{ifile} >>#{dataFile} 2>&1")
+                  :input => "-#{operation} -selname,#{varname} #{ifile} >>#{dataFile} 2>/dev/null")
 
     # postprocessing for correct time values
     data = []
     File.open(dataFile).each_with_index {|line,lineIndex|
       next if line.chomp.empty?
-      _t = line.chomp.gsub(/ +/,'|').split('|')
+      _t = line.chomp.sub(/^ /,'').gsub(/ +/,'|').split('|')
       if 0 == lineIndex then
         data << _t
         next
@@ -131,14 +131,14 @@ class IconPlot < Struct.new(:caller,:plotter,:libdir,:otype,:display,:cdo,:debug
       if "0" == _t[1] then
         _t[1] = '00:00:00'
       else
-        time = _t[1].reverse
-        timeStr = ''
-        while time.size > 2 do
-          timeStr << time[0,2] << ':'
-          time = time[2..-1]
-        end
-        timeStr << time.ljust(2,'0') unless time.size == 0
-        _t[1] = timeStr.reverse
+       #time = _t[1].reverse
+       #timeStr = ''
+       #while time.size > 2 do
+       #  timeStr << time[0,2] << ':'
+       #  time = time[2..-1]
+       #end
+       #timeStr << time.ljust(2,'0') unless time.size == 0
+       #_t[1] = timeStr.reverse
       end
       data << _t
     }

data/lib/icon_plot.ncl

@@ -55,7 +55,8 @@
   end if
 ;------------------------------------------------------------------------------
 ;==============================================================================
-  VERSION = "0.0.6"
+  VERSION = "0.0.9"
+;  setfileoption("nc", "filestructure", "advanced")
 ;----- MAIN PROGRAM -----------------------------------------------------------
   if (isvar("help")) then
     dq = str_get_dq()
@@ -63,7 +64,7 @@
     print("; Basic Usage:")
     print(";   ncl icon_plot.ncl 'iFile=" + dq + "path/file.nc" + dq + "' 'varName=" + dq + "ELEV" + dq + "' 'oFile=" + dq + "test" + dq + "' timeStep=1 levIndex=1")
     print(";")
-    print("; Select an special area of the globe with mapLLC (LowerLeftCorner) and mapURC (UpperRightCorner)   ")
+    print("; Select a special area of the globe with mapLLC (LowerLeftCorner) and mapURC (UpperRightCorner)   ")
     print(";   ncl icon_plot.ncl 'iFile=" + dq + "path/file.nc" + dq + "' 'varName=" + dq + "W" + dq + "'    'mapLLC=(/35.0, -8/)' 'mapURC=(/55, 8/)'")
     print(";")
     print("; Masking:")
@@ -98,7 +99,9 @@
     print(";")
     print("; Optional Parameter:")
     print("; oType         : output graphics format (ps, eps, png, default: eps)")
-    print("; resolution    : resolution string whish is used for remapping the icon data")
+    print("; xsize         : width of plot (PNG only)")
+    print("; ysize         : height of plot (PNG only)")
+    print("; resolution    : resolution string which is used for remapping the icon data")
     print(";                 to regular grid (default: r90x45)")
     print("; vecRefLength  : reference vector length (default: 8.0)")
     print("; vecColByLen   : should vectors get coloured according to their lengths (default:False)")
@@ -112,15 +115,15 @@
     print("; mapLLC        : (lon,lat) value array of the Lower Left Corner of the map")
     print("; mapURC        : (lon,lat) ------- || ------- Upper Right Corner of the map")
     print(";                  use -180 to 180 for longitude !")
-    print("; secLC, secRC  : start and and point of a vertical section (no default value)")
-    print(";                 section is plotted it secLC and secRC are present")
-    print("; secMode       : perform transection along a great crice or on a straight line")
+    print("; secLC, secRC  : start and end point of a vertical section (no default value)")
+    print(";                 section is plotted if secLC and secRC are present")
+    print("; secMode       : perform transection along a great cricle or on a straight line")
     print(";                 possible values: circle, straight (default)")
     print(";               : ATTENTION: using 'circle' for sections with a small angle wrt the equator is NOT RECOMMENDED")
     print("; showSecMap    : display the location of the vertical section in an additional map (default:true)")
     print("; centerLon     : center longitute for satelite view (default:30.)")
     print("; centerLat     : center latitude for satelite view (default:20.)")
-    print("; satDist       : satelite distance (default:20.0)")
+    print("; satDist       : satellite distance (default:20.0)")
     print("; scaleFactor   : optional scale factor")
     print("; fillMode      : how to color contour plots,i.e. sets cnFillMode: CellFill or RasterFill (default:RasterFill)")
     print("; selMode       : mode for level spacing of the plot:")
@@ -128,14 +131,17 @@
     print(";                 manual (automatic linear spacing with usage of minVar and maxVar)")
     print(";                 auto (default: let do ncl the spacing)")
     print("; scaleLimit    : Limits the number of decades for levels in selMode-halflog")
-    print("; minVar/maxVar : min/maximal Value to plot for selMode = 'manual' or 'halflog'")
+    print("; minVar/maxVar : min/maximal value to plot for selMode = 'manual' or 'halflog'")
     print("; numLevs       : set a number of labels for manual selMode only. NumCol=numLevs+2")
     print(";                 (Not set directly, but used for computing the contour level spacing. default=10)")
     print("; plotLevs      : set individual contour levels, selMode is set to 'manual' automatically")
     print("; mapType       : projection type (default: lonlat), other: ortho (not compatible with showGrid), NHps")
     print(";                 and SHps (polar stereographic projections centered at the north and south pole respectively),")
     print(";                 sat (satelite view, use option centerLat, centerLon and satDist)")
-    print("; hov           : plot hovmoeller diagram, expects variable with two 'active' dimensions, i.e. var(:,:,0,0)")
+    print("; hov[Depth]    : plot hovmoeller diagram, expects variable with only two 'active' dimensions, i.e. time/depth")
+    print("; hovCut        : plot hovmoeller diagram, input is 'levIndex' and 'hovLC' 'hovRC'")
+    print("; hovLC, hovRC  : start and end point of a section for hovCut - currently same longitude allowed only,")
+    print(";                 the Northern and Southern values are used as plotting boundaries; switches to hovCut")
     print("; mapLine       : (logical) draws continent lines on the plot (foreground/transparent)")
     print("; maskName      : variable to mask with. maskName is expected NOT to have time dimension")
     print(";                  'maskName=" + dq + "none" + dq + "' is same as default (no mask variable)")
@@ -146,17 +152,17 @@
     print("; lStrg         : left string")
     print("; rStrg         : right string")
     print("; tStrg         : title string")
-    print("; bStrg         : base string - default is prgr name and time stamp only")
+    print("; bStrg         : base string - default is program name and time stamp only")
     print("; maxView       : (logical) maximize plot area on paper (not for buildbot -> convert to png)")
     print("; colormap      : string for predefined colormaps of the plot (e.g. 'colormap=" + dq + "BlAqGrYeOrReVi200" + dq + "')")
     print("; withLines     : draw lines between each contour color (default:True)")
-    print("; withLineLabels: label the countour lines with theri values (auto-enable withLines, default:False)")
+    print("; withLineLabels: label the countour lines with their values (auto-enable withLines, default:False)")
     print("; showGrid      : display polygon lines with filled colors instead of contour plot")
     print("; markCells     : mark cell centers with dots")
     print("; showNcd       : display NDC Grid to find Normalized Device Coordinates on the plot")
     print("; k2c           : if True, perform Kelvin2Celsius shift (default:True)")
     print(";")
-    print("; ATMOSPEHRE STUFF:")
+    print("; ATMOSPHERE STUFF:")
     print("; atmLev        : chooses vertical plot levels for atmosphere input (h: height,p:pressure,m:modellevel,default:p)")
     print(";---------------------------------------------------------------")
     print("; altLibDir     : Alternative directory for loading icon_plot_lib.ncl")
@@ -374,30 +380,37 @@
     end if
   end if
 
-  if ( (isvar("secLC").or.isvar("secRC")).ne.(isvar("secLC").and.isvar("secRC"))) then
-    print("Please provide secLC AND secRC for generatin a section plot!")
+  if ((isvar("secLC").or.isvar("secRC")).ne.(isvar("secLC").and.isvar("secRC"))) then
+    print("Please provide secLC AND secRC for generating a section plot!")
     exit
   end if
   if (isvar("secLC").and.isvar("secRC")) plotMode = "section" end if
-  if (isvar("hov")) plotMode = "hovmoeller" end if
+
+  if ((isvar("hovLC").or.isvar("hovRC")).ne.(isvar("hovLC").and.isvar("hovRC"))) then
+    print("Please provide hovLC AND hovRC for generating a Hovmoeller plot!")
+    exit
+  end if
+  if (isvar("hovLC").and.isvar("hovRC")) plotMode = "hovmoeller" end if
+  if (isvar("hovDepth").or.isvar("hovCut").or.(isvar("hov"))) plotMode = "hovmoeller" end if
 
   if (plotMode.eq."vector" .or. plotMode.eq."overlay" .or. plotMode.eq."section") then
     if (horizonalGridType .ne. "unstructured") then ; regular grid expected, i.e. remappgin is NOT required
       rFile = addfile(iFile+".nc","r")
     else
-      ; performe some remapping to a regular grid because ncl cannot draw vector
+      ; perform some remapping to a regular grid because ncl cannot draw vector
       ; from unstructured grids
       ; TODO addVars=(/ "PS","PHIS"/)
-      remapFilename = setRemapFilename(iFile,iType,resolution,atmLev)
+      operator = "remapnn"
+      remapFilename = setRemapFilename(iFile,iType,resolution,atmLev,operator)
       print("remapFilename:"+remapFilename)
       if (.not. checkRemappedFile(iFile,remapFilename,ivarname,iType,atmLev,atmPLevs,atmHLevs) ) then
         addVars=(/""/)
         print("PERFORM remapnn again!")
-        remapForVecPlot(iFile,remapFilename,resolution,useMask,plotMode,DEBUG,addVars)
+        remapForVecPlot(iFile,remapFilename,resolution,useMask,plotMode,DEBUG,addVars,operator)
       end if
       rFile = addfile( remapFilename+".nc", "r" )
     end if
-   ;if (plotMode .eq. "hoffmueller") then ; perform zonmean
+   ;if (plotMode .eq. "hovmoeller") then ; perform zonmean
    ;  zonmeanFilename = setZonmeanFilename(remapFilename,ivarname,atmLev)
    ;  zonmean4HoffmuellerPlot(remapFilename,ivarname,zonmeanFilename)
    ;end if
@@ -437,6 +450,7 @@
     print("oFile       = "+oFile)
     print("Graphics format is " +oType)
     print("plotMode    = "+plotMode)
+    print("showGrid    = "+showGrid)
 
     if (plotMode.eq."scalar" .or. plotMode.eq."section" .or. plotMode.eq."hovmoeller") then
       print("varName     = "+varName)
@@ -465,24 +479,6 @@
 
     var = selField(varName,File,timeStep,levIndex,horizonalGridType)
 
-    if (horizonalGridType .eq. "unstructured") then
-      if (isvar("gridFile")) then
-        x      = getCoordinateFromFile(lonCo,gridFile)
-        y      = getCoordinateFromFile(latCo,gridFile)
-        if (showGrid) then
-          bounds = getBoundsFromFile(lonCo,latCo,gridFile)
-        end if
-      else
-        lonlatSizes = getCoordSizes(var,File)
-        x           = new(lonlatSizes(0),double)
-        y           = new(lonlatSizes(1),double)
-        getLonLats(var,File,x,y)
-        if (showGrid) then
-          bounds  = getBoundsOfCoordinates(var,File)
-        end if
-      end if
-    end if
-
     scaleVar(var,scaleFactor)
 
     ; set minVar_maxVar for plotting
@@ -562,6 +558,10 @@
   ; make the plot
   ;---------------------------------------------------------------
   ; preparations
+  if (oType .eq. "png") then
+  if(isvar("xsize")) oType@wkWidth  = xsize end if
+  if(isvar("ysize")) oType@wkHeight = ysize end if
+  end if
   wks = gsn_open_wks(oType,oFile)
 
   gsn_define_colormap(wks,colormap)
@@ -570,7 +570,7 @@
     drawNDCGrid(wks)
   end if
 
-  ResC = setDefaultResource(True,withLines,withLineLabels,fillMode)
+  ResC = setDefaultResource(True,withLines,withLineLabels,fillMode,mapType)
 
   if (.not. userColors) setDefaultColors(ResC,minVar,maxVar,DEBUG) end if
 
@@ -578,10 +578,28 @@
 
   if (useMask) setMaskColor(wks,ResC) end if
 
-
-  if (getHorizGridType(iFile,ivarname,isIcon) .eq. "unstructured") then
+  if (horizonalGridType .eq. "unstructured" .or. horizonalGridType .eq. "curvilinear") then
     if (plotMode.eq."scalar" .or. plotMode.eq."overlay")
-      setCoordinates(ResC,x,y)
+      if (horizonalGridType .eq. "unstructured") then
+        if (isvar("gridFile")) then
+          ;TODO: fix x,y creation for external gridfile (plot restart data)
+          coords = setCoordinatesFromExtraFile(lonCo,latCo,gridFile,var,ResC)
+          if (showGrid) then
+            bounds = getBoundsFromFile(lonCo,latCo,gridFile)
+          end if
+        else
+          if (showGrid) then
+            bounds  = getBoundsOfCoordinates(var,File)
+          end if
+        end if
+      end if
+    end if
+    if (.not.isvar("coords")) then
+      if (isvar("rFile")) then
+        coords = setCoordinates(ResC,rFile,var)
+      else
+        coords = setCoordinates(ResC,File,var)
+      end if
     end if
   end if
 
@@ -608,7 +626,11 @@
 
   setMapType(ResC,mapType,centerLon,centerLat,satDist)
 
-  selMapCut(ResC,mapLLC,mapURC)
+  if (isvar("limitMap")) then
+    limitMap2Coords(ResC)
+  else
+    selMapCut(ResC,mapLLC,mapURC)
+  end if
 
   setMapVisibility(ResC,mapLine)
 
@@ -654,29 +676,67 @@
     draw(plot)
   end if
   if (plotMode.eq."scalar") then
-    if (DEBUG .and. horizonalGridType.eq."unstructured") print("Gridtype is "+getGridType(var)) end if
+    ;if (DEBUG .and. horizonalGridType.eq."unstructured") print("Gridtype is "+getGridType(var)) end if
+      if ( horizonalGridType .eq. "curvilinear") then
+        ResC@cnFillMode           = "CellFill"
+      end if
+    if (showGrid ) then
+      if ( horizonalGridType .eq. "curvilinear") then
+        ResC@cnLinesOn            = False
+        ResC@cnFillMode           = "CellFill"
+        ResC@cnCellFillEdgeColor   = "black"          ;-- set edge color (default: -1 transparent)
+        ResC@cnCellFillMissingValEdgeColor = "grey"   ;-- set miss val edge color (default: -1 transparent)
+      end if
+    end if
     plot = gsn_csm_contour_map(wks,var,ResC)
-
     if (showGrid) then
-      if ( horizonalGridType .ne. "unstructured") then
+      if (horizonalGridType .eq. "unstructured") then
+        ;=======================================================================================
+        ;plotGrid(wks,plot,var,coords(1,:),bounds,File,ResC,DEBUG)
+        boundslon = bounds(0,:,:)
+        boundslat = bounds(1,:,:)
+        pres = plotGrid_62(wks,plot,var,boundslon,boundslat)
+        gsid = gsn_add_polygon(wks,plot,ndtooned(boundslon),ndtooned(boundslat),pres)
+        ; }}}
+        ;=======================================================================================
+      end if
+      if ( horizonalGridType .ne. "curvilinear" .and. horizonalGridType .ne. "unstructured") then
         print("#=====================================================================")
         print("Grid representation is only implemented for the unstructured ICON grid")
         exit
       end if
-      plotGrid(wks,plot,var,x,bounds,File,ResC,DEBUG)
     end if
     if ( markCells ) then
+      x = coords(0,:)
+      y = coords(1,:)
       plotCellMarkers(x,y,wks,plot,True)
     end if
+    draw(plot)
   end if
   if (plotMode.eq."scatter") then
     plot = gsn_csm_xy(wks,xvar,yvar,ResC)
   end if
   if (plotMode .eq. "hovmoeller") then ; create hoffmuelle diagram
-    ; select the whole time series
-    ; and plot it with
-    ResC@trYReverse    = True
-    plot = gsn_csm_hov(wks, var(depth|:,time|:,lat|0,lon|0), ResC)
+    if (isvar("hovDepth").or.(isvar("hov"))) then
+      ; select the whole time series
+      ; and plot it with
+      ResC@trYReverse    = True
+      plot = gsn_csm_hov(wks, var(depth|:,time|:,lat|0,lon|0), ResC)
+    end if
+    if (isvar("hovRC")) then
+      startLon = hovRC(0)
+      northLat = hovRC(1)
+      endLon   = hovLC(0)
+      southLat = hovLC(1)
+      hovCut = True
+    end if
+    if (isvar("hovCut")) then
+      ;plot = gsn_csm_hov(wks, var(depth|0,time|:,lat|:,lon|0), ResC)
+      ;plot = gsn_csm_lat_time(wks, var(depth|0,time|:,lat|:,lon|0), ResC)
+      ;plot = gsn_csm_lat_time(wks, var(depth|0,lat|:,lon|165,time|:), ResC)
+      ;plot = gsn_csm_hov(wks, var({depth|150},{lat|70:-80},{lon|330},time|:), ResC)
+      plot = gsn_csm_hov(wks, var(depth|levIndex,{lat|northLat:southLat},{lon|startLon},time|:), ResC)
+    end if
   end if
   if (plotMode.eq."section") then
 
@@ -687,7 +747,7 @@
     end if
     points   = ispan(0,secPoints-1,1)*1.0
 
-    res  = setDefaultSectionResource(points,secPoints,secLC,secRC)
+    res  = setDefaultSectionResource(points,secPoints,secLC,secRC,withLines,withLineLabels)
     setAutomaticPlotCaptions(res,plotMode,varName,File,iFile,timeStep,levIndex,iType,atmLev,k2c)
     setLevels(selMode,res,minVar,maxVar,scaleLimit,numLevs,DEBUG)
     setSectionVertLabel(res,iType,atmLev)

data/lib/icon_plot_lib.ncl

@@ -1,22 +1,22 @@
 ;---------------------------------------------------------------
 ; ICON PLOT LIB
 ; This library aimes to provide an abstract way to access and plot ICON data
-; fields from NetCDF files. It should obay some general requirements:
-; * Don-Repeat-Yourself-principle: Every method has one responsebility, which no
+; fields from NetCDF files. It should obey some general requirements:
+; * Dont-Repeat-Yourself-Principle: Every method has one responsibility, which no
 ;   other method has. Or in other words: Don't just CopyAndPaste.
 ; * Use as much information from the input files as possible. If something is
-;   missing, create a better file or let the user input this peace of
-;   information. Examples are coordinates of data variabels and bounds of them.
+;   missing, create a better file or let the user input this piece of
+;   information. Examples are coordinates of data variables and bounds of them.
 ;   There are attributes for this.
 ; * naming conventions:
 ;   # check* methods are procedure which can exit the whole program
-;   # get* methods are functions with return, what they descibe with their
+;   # get* methods are functions with return, what they describe with their
 ;     name. This is similar to NCL itself, e.g. getfilevardims()
 ;   # set* methods are procedures
 ;---------------------------------------------------------------
 ; Authors       : Ralf Mueller (ralf.mueller@zmaw.de)
 ;                 Stephan Lorenz (stephan.lorenz@zmaw.de)
-; VERSION = "0.0.6"
+; VERSION = "0.0.9"
 ;-------------------------------------------------------------------------------
 ; some global defaults
   DEFAULTLON   = "clon"
@@ -29,9 +29,10 @@
   DEFAULTCOLOR = "BlueDarkRed18"      ; ncolors = 256
   DEFAULTCOLOR = "BlueWhiteOrangeRed" ; ncolors = 256
 ;-------------------------------------------------------------------------------
-procedure setCDO(path)
+undef("setCDO")
+procedure setCDO(path2Cdo)
 begin
-  CDO=path
+  CDO=path2Cdo
 end
 ;-------------------------------------------------------------------------------
 ; signum function
@@ -44,16 +45,62 @@ begin
     return(val/abs(val))
   end if
 end
+;---------------------------------------------------------------
+; return ncl version string
+undef("getNclVersion")
+function getNclVersion()
+begin
+  version_string = systemfunc("ncl -V")
+  version_array  = str_split(version_string, ".")
+  ; use integers instead
+  return stringtointeger(version_array)
+end
+;---------------------------------------------------------------
+; return true if ncl uses groups for netcdf4 input (release 6.1.1 with netcdf4 input)
+undef("versionWithNC4GroupNotationSupport")
+function versionWithNC4GroupNotationSupport()
+begin
+  return (/6,1,1/)
+end
+undef("versionIsGE")
+function versionIsGE(otherVersion)
+begin
+  return True
+  diff = otherVersion .gt. getNclVersion()
+  print(""+diff)
+  do i=0,dimsizes(diff)
+    if (.not. diff(i)) then
+      return False
+    end if
+  end do
+  return True
+end
+;-------------------------------------------------------------------------------
+; return a filehandle that can access variables names without leading '/' even
+; if the input it nc4 and ncl has version 6.1.1 or above
+undef("nc4Aware_addfile")
+function nc4Aware_addfile(filename,debug)
+begin
+  filehandle = addfile( filename, "r" )
+  filetype   = get_file_version(filehandle)
+  if debug then
+    print("[nc4Aware_addfile]Read in file:'" + filename +"'")
+    print("[nc4Aware_addfile]    Filetype:'" + filetype +"'")
+    print("[nc4Aware_addfile]NCL  version:'" + str_join(getNclVersion(),".") +"'")
+    print("[nc4Aware_addfile]version > 611'" + versionIsGE((/6,1,1/)) )
+  end if
+  if ((get_file_version(filehandle) .eq. "NETCDF4") .and. versionIsGE((/6,1,1/))) then
+    print("[nc4Aware_addfile]switch to root group '/'")
+   ;filehandle = filehandle=>/
+  end if
+  return filehandle
+end
 ;-------------------------------------------------------------------------------
 ; Is the 'time' variable a dimesions in the input dimensions array?
 undef("has_time")
 function has_time(dims)
 begin
-  if (any(dims.eq."time")) then
-    return(True)
-  else
-    return(False)
-  end if
+  return any(dims.eq."time")
 end
 ;-------------------------------------------------------------------------------
 ; Has the input variable the dimension 'time'?
@@ -74,7 +121,7 @@ function has_var(filehandle,varname)
 begin
   vNames = getfilevarnames(filehandle)
   n      = dimsizes (vNames)
-  do i=0,n-1              ; loop thru each variable
+  do i=0,n-1              ; loop through each variable
     if (vNames(i).eq.varname) return True end if
   end do
   return False
@@ -103,7 +150,7 @@ begin
 end
 ;-------------------------------------------------------------------------------
 ; Read a horizontal field of the variable 'varname' from the given filehandle.
-; If it is time dependend or has a vertival axis, use the given values for
+; If it is time dependant or has a vertical axis, use the given values for
 ; timstep and levelindex for selecting the field.
 ;
 ; Currently this limited to 3 dimensions, so that 4D tracer variables are ignored
@@ -368,15 +415,17 @@ end
 ;-------------------------------------------------------------------------------
 ; Return the coordinates of a given variable in degrees
 undef("getLonLats")
-procedure getLonLats(variable,filehandle,x,y)
+procedure getLonLats(variable,filehandle,x,y,isicon)
 begin
   lonlat = str_split(variable@coordinates," ")
   lon    = lonlat(0)
   lat    = lonlat(1)
   x      = filehandle->$lon$
   y      = filehandle->$lat$
-  x      = x * RAD2DEG
-  y      = y * RAD2DEG
+  if (isicon) then
+    x      = x * RAD2DEG
+    y      = y * RAD2DEG
+  end if
 end
 ;-------------------------------------------------------------------------------
 ; get variable from a certain file
@@ -401,7 +450,7 @@ end
 ; Examples:
 ;   (/1,2,5,10,20,50,100,200,500/) or
 ;   (/-1,-5e-1,-2e-1,-1e-1,-5e-2,-2e-2,-1e-2,0,1e-2,2e-2,5e-2,1e-1,2e-1,5e-1,1/)
-; scaleLimit determines the numer of decades on the positive or negtive axis of the output array
+; scaleLimit determines the number of decades on the positive or negative axis of the output array
 undef("createLevels")
 function createLevels(minVar, maxVar, scaleLimit)
 begin
@@ -486,6 +535,7 @@ begin
 end
 ;---------------------------------------------------------------
 ; Print information about dimensions and attributes of a given varaiable
+undef("printVar")
 procedure printVar(varname, filehandle)
 begin
   dims  = getfilevardims(filehandle,varname)
@@ -590,9 +640,9 @@ end
 ;---------------------------------------------------------------
 ; set filename for the automatically remapped file
 undef("setRemapFilename")
-function setRemapFilename(filename,itype,resolution,atmlev)
+function setRemapFilename(filename,itype,resolution,atmlev,operator)
 begin
-  return setNewFilename(filename,"remapnn_"+resolution,itype,atmlev)
+  return setNewFilename(filename,operator+"_"+resolution,itype,atmlev)
 end
 undef("setZonmeanFilename")
 function setZonmeanFilename(filename,itype,atmlev)
@@ -600,7 +650,7 @@ begin
   return setNewFilename(filename,"zonmean",itype,atmlev)
 end
 ;---------------------------------------------------------------
-; Check, if two files have the same number of timestemps. Exit otherwise
+; Check, if two files have the same number of time steps. Exit otherwise
 undef("checkRemappedFile")
 function checkRemappedFile(infilename,remapfilename,varname,itype,atmlev,atmplevs,atmhlevs)
 begin
@@ -676,13 +726,13 @@ end
 ; Perform a remapping (wich CDO) to a regular grid with a given resolution and return the
 ; filename of the remapped file. If the file is already present, the name is returned only
 undef("remapForVecPlot")
-procedure remapForVecPlot(iFile,remapFilename,resolution,useMask,plotMode,debug,addvars)
+procedure remapForVecPlot(iFile,remapFilename,resolution,useMask,plotMode,debug,addvars,operator)
 begin
   if (plotMode.eq."section") then vecVars=(/varName/) end if
   print("#=====================================================================================")
   print("Looking for remapped file: "+remapFilename)
   print("Use CDO to perform remapping: Create intermediate file: "+remapFilename)
-  print("Remove this intermediate file, if it was NOT automatically created from your inpur file "+iFile+"!")
+  print("Remove this intermediate file, if it was NOT automatically created from your input file "+iFile+"!")
   variables = str_join(vecVars,",")
   if ( addvars(0) .ne. "")
     variables = str_concat((/variables,",",str_join(addvars,",")/))
@@ -694,7 +744,7 @@ begin
   end if
   if (plotMode.eq."overlay") variables = variables+","+varName end if
 
-  cmd = CDO+" -P 8 -remapnn,"+resolution+" -selname,"+variables+" "+iFile+" "+remapFilename
+  cmd = CDO+" -P 8 -"+operator+","+resolution+" -selname,"+variables+" "+iFile+" "+remapFilename
   if debug then
     print(cmd)
   end if
@@ -712,7 +762,7 @@ begin
   system(cmd)
 end
 ;---------------------------------------------------------------
-; Return the bounds of the coordinates of agiven variable
+; Return the bounds of the coordinates of a given variable
 undef("getBoundsOfCoordinates")
 function getBoundsOfCoordinates(variable,filehandle)
 begin
@@ -729,7 +779,7 @@ begin
     return (/boundslon, boundslat/)
 end
 ;---------------------------------------------------------------
-; Return the bounds of the coordinates of agiven variable
+; Return the bounds of the coordinates of a given variable
 undef("getBoundsFromFile")
 function getBoundsFromFile(lonname,latname,filename)
 begin
@@ -839,23 +889,12 @@ begin
   end if
 end
 ;---------------------------------------------------------------
-; Create a default NCL resource
-undef("setDefaultResource")
-function setDefaultResource(verticallabelbar,withLines,withLineLabels,fillmode)
+; set the lines and their labels
+undef("setLineAndLineLabels")
+procedure setLineAndLineLabels(resource,withLines,withLineLabels)
 begin
-  resource                          = True
-  resource@gsnMaximize              = False
-  resource@gsnFrame                 = False
-  resource@gsnDraw                  = True
-  resource@gsnSpreadColors          = True
-  resource@cnFillOn                 = True
-
 ; lines/labels inside the plot ================================================
-  if (withLines) then
-    resource@cnLinesOn              = True
-  else
-    resource@cnLinesOn              = False
-  end if
+  resource@cnLinesOn                = False
   resource@cnLineLabelsOn           = False
   if (withLineLabels ) then
     withLines                       = True
@@ -864,7 +903,23 @@ begin
     resource@cnLineLabelFontThicknessF = 0.006
   end if
   if (withLines) resource@cnLinesOn = True end if
+end
+;---------------------------------------------------------------
+; Create a default NCL resource
+undef("setDefaultResource")
+function setDefaultResource(verticallabelbar,withLines,withLineLabels,fillmode,maptype)
+begin
+  resource                          = True
+  resource@gsnMaximize              = False
+  resource@gsnFrame                 = False
+  resource@gsnDraw                  = True
+  resource@gsnSpreadColors          = True
+  resource@cnFillOn                 = True
+
+  resource@mapType                  = maptype
 
+; lines/labels inside the plot ================================================
+  setLineAndLineLabels(resource,withLines,withLineLabels)
   ; representation of missing values ===========================================
   delete(resource@cnMissingValFillPattern);     no longer filling missing value areas
   delete(resource@cnMissingValFillColor);       no longer filling missing value areas
@@ -909,7 +964,7 @@ begin
 end
 ;---------------------------------------------------------------
 undef("setDefaultSectionResource")
-function setDefaultSectionResource(points,secpoints,seclc,secrc)
+function setDefaultSectionResource(points,secpoints,seclc,secrc,withLines,withLineLabels)
 begin
   resource                              = True;           plot mods desired
   resource@gsnFrame                     = False;                don't turn page yet
@@ -930,10 +985,6 @@ begin
   resource@cnFillOn                     = True;          turn on color
   resource@lbLabelAutoStride            = True;          nice label bar label stride
   resource@gsnSpreadColors              = True;          use full range of colormap
-  resource@cnLinesOn                    = False;         turn off countour lines
- ; if (useContourLines)
-  resource@cnLinesOn             = False
-  resource@cnLineLabelsOn        = False
   resource@lbOrientation                = "vertical";    vertical label bar
   resource@pmLabelBarOrthogonalPosF     = -0.05;         move label bar closer to plot
   resource@lbTitlePosition              = "Bottom"
@@ -953,6 +1004,9 @@ begin
   ;resource@gsnYAxisIrregular2Log        = True
   resource@cnFillMode           = "RasterFill"
   resource@cnRasterSmoothingOn  = False
+
+; lines/labels inside the plot ================================================
+  setLineAndLineLabels(resource,withLines,withLineLabels)
   return resource
 end
 ;---------------------------------------------------------------
@@ -1026,7 +1080,8 @@ end
 undef("setBaseString")
 procedure setBaseString(workstation,plot,resource,basestring)
 begin
-  resource@txFontHeightF = 0.01 ; if this is changed, please adjust the implementation of getLshiftForNDCString()
+  resource@txFontHeightF = 0.01 ; if this is changed, please adjust 
+                                ; the implementation of getLshiftForNDCString()
   gsn_text_ndc(workstation,basestring,getLshiftForNDCString(basestring),.01,resource)
 end
 ;---------------------------------------------------------------
@@ -1268,10 +1323,12 @@ end
 undef("checkMinMaxVar")
 procedure checkMinMaxVar(minvar,maxvar)
 begin
+  if ( (.not. ismissing(minvar)) .and. (.not. ismissing(maxvar)) ) then
   if ( minvar .gt. maxvar ) then
     print("minVar has to be larger than maxVar")
     exit
   end if
+  end if
 end
 ;---------------------------------------------------------------
 ; set the color of masked locations
@@ -1282,13 +1339,53 @@ begin
   resource@cnMissingValPerimOn         = True;          turn on missing value outlines
   resource@cnMissingValPerimThicknessF = 2.5;   set the missing value outline thickness
 end
+undef("setCoordsForResouce")
+procedure setCoordsForResouce(x,y,variable,resource)
+begin
+  if (dimsizes(dimsizes(x)).ne.1) then
+     variable@lon2d = x
+     variable@lat2d = y
+  end if
+     resource@sfXArray = x
+     resource@sfYArray = y
+end
 ;---------------------------------------------------------------
 ; set the coordinates for an ICON data field
 undef("setCoordinates")
-procedure setCoordinates(resource,lons,lats)
+function setCoordinates(resource,filehandle,variable)
 begin
-  resource@sfXArray = lons
-  resource@sfYArray = lats
+  ret = True
+  if (isatt(variable,"coordinates")) then
+    lonlat = str_split(variable@coordinates," ")
+    lon    = lonlat(0)
+    lat    = lonlat(1)
+  else
+    lon = "lon"
+    lat = "lat"
+    ret = False
+  end if
+  x      = filehandle->$lon$
+  y      = filehandle->$lat$
+  if ("radian" .eq. x@units) x = x * RAD2DEG end if
+  if ("radian" .eq. y@units) y = y * RAD2DEG end if
+
+  setCoordsForResouce(x,y,variable,resource)
+
+  if (ret) then
+  return (/x,y/)
+  else
+  return (/1,1/)
+  end if
+end
+undef("setCoordinatesFromExtraFile")
+function setCoordinatesFromExtraFile(lon_name,lat_name,filename,variable,resource)
+begin
+  x = getCoordinateFromFile(lon_name,filename)
+  y = getCoordinateFromFile(lat_name,filename)
+
+  setCoordsForResouce(x,y,variable,resource)
+
+  return (/x,y/)
 end
 ;---------------------------------------------------------------
 ; set the coordinates bounds for an ICON data field
@@ -1361,7 +1458,7 @@ end
 undef("getVarLevelType")
 function getVarLevelType(filename,varname)
 begin
-  f       = addfile(filename+".nc","r")
+  f       = addfile(filename,"r")
   vertdim = getVertDim(f,f->$varname$)
   if (ismissing(vertdim)) then
     delete(f)
@@ -1385,7 +1482,7 @@ begin
     gridtype = "unstructured"
   else
     print("CDO:"+CDO)
-    print("Call cdo griddes for determine the horizontal gridtype")
+    print("Call cdo griddes to determine the horizontal gridtype")
     gridtype = systemfunc(CDO+" -griddes -selname,"+varname+" -seltimestep,1 "+filename+" | grep gridtype | cut -d ' ' -f 4")
   end if
   print("Found horizontal grid of type: "+gridtype)
@@ -1408,7 +1505,7 @@ begin
 
   lats = fspan(startLat,endLat,npoints)
   lons = fspan(startLon,endLon,npoints)
-  ; dont let the lons be greater than 360, because theses values will not be in
+  ; don't let the lons be greater than 360, because these values will not be in
   ; the data. Otherwise an interpolation error will occur
   lons = where(lons.gt.360, lons-360.0, lons)
 
@@ -1531,7 +1628,7 @@ begin
   var_at_p@units     = var@units
   var_at_p@long_name = var@long_name
 
-  extrapolate = True ; switch for doing extrapolation below the ground:
+  extrapolate = True ; switch to do extrapolation below the ground:
   intmethod   = 1 ; 1: method for temperature, -1: method for geopotential, 0: other vars
   tlow =  var(nlevels-1,:,:)  ; temperature at lowest model level
                                        ; (use ground temperature if available)
@@ -1641,7 +1738,7 @@ begin
   var_at_p@units     = theta3d@units
   var_at_p@long_name = theta3d@long_name
 
-  extrapolate = True ; switch for doing extrapolation below the ground:
+  extrapolate = True ; switch to do extrapolation below the ground:
   intmethod   = 1 ; 1: method for temperature, -1: method for geopotential, 0: other vars
   tlow =  theta3dll(nlevels-1,:,:)  ; temperature at lowest model level
                                        ; (use ground temperature if available)
@@ -1661,6 +1758,11 @@ undef("getFlags")
 function getFlags(var,boundslon,boundslat,resource)
 begin
   flags = new(dimsizes(var),logical,"No_FillValue")
+
+  if (resource@mapType .EQ. "NHps" .OR. resource@mapType .EQ. "SHps") then
+    flags = True
+    return flags
+  end if
   do i = 0,dimsizes(var) - 1
     flags(i) = where(all(boundslon(i,:) .gt. resource@mpMaxLonF) .or. \
                      all(boundslon(i,:) .lt. resource@mpMinLonF) .or. \
@@ -1764,7 +1866,7 @@ begin
   if (.not. ismissing(vhig(0))) then
     do j = 0, dimsizes(vhig) -1
       pres@gsFillColor = colors(i+1)
-      if ( flags(vhig(i))) then
+      if ( flags(vhig(j))) then
       gsn_polygon(wks,plot,boundslon(vhig(j),:),boundslat(vhig(j),:),pres)
       end if
     end do
@@ -1811,6 +1913,50 @@ begin
   performGridPlot(variable,levels,colors,boundslon,boundslat,workstation,plot,flags,debug)
 end
 ;---------------------------------------------------------------
+; Plot the ICON grid: New version, req. ncl 6.2
+undef("plotGrid_62")
+function plotGrid_62(wks,plot,var,boundslon,boundslat)
+begin
+  pres             = True
+  pres@gsEdgesOn   = True   ; Turn on edges
+  pres@gsFillIndex = 0      ; Solid fill, the default
+  ;---Create color array for triangles
+  ntri     = dimsizes(var)                  ;-- Number of triangles
+  gscolors = new(ntri,integer)
+  gscolors = -1                           ;-- Initialize to transparent
+
+  getvalues plot@contour
+   "cnLevels"     : levels
+   "cnFillColors" : colors
+  end getvalues
+  ;---All triangles less than lowest level
+  vind = ind(var .lt. levels(0))
+  if(.not.all(ismissing(vind))) then
+    gscolors(vind) = colors(0)
+  end if
+
+  ;---All triangles inbetween levels
+  do i = 1, dimsizes(levels) - 1
+    vind := ind(var .ge. levels(i-1) .and. var .lt. levels(i))
+    if(.not.all(ismissing(vind))) then
+      gscolors(vind) = colors(i)
+    end if
+  end do
+
+  ;---All triangles higher than highest level
+  delete(vind)
+  vind = ind(var .ge. levels(dimsizes(levels)-1))
+  if(.not.all(ismissing(vind))) then
+    gscolors(vind) = colors(dimsizes(levels)-1)
+  end if
+
+  pres@gsColors = gscolors
+  pres@gsSegments = ispan(0,dimsizes(var) * 3,3)
+
+  return pres
+  ; }}}
+end
+;---------------------------------------------------------------
 ; read in mask variable
 undef("readMaskVar")
 function readMaskVar()
@@ -1822,10 +1968,11 @@ undef("plotVecOrStream")
 function plotVecOrStream(useScalar,useStreamlines,uvar,vvar,scalarvar,resource,workstation,showMap)
 begin
   ; streamline setup
-  resource@stArrowStride    = 3;           arrows start every third
-  resource@stArrowLengthF   = 0.002;       default is dynamic
-  resource@stMinDistanceF   = 0.025;       streamline dist
-  resource@stLineThicknessF = 2.0;         streamline thickness
+  resource@stArrowStride    = 4;           arrows start every third
+  resource@stArrowLengthF   = 0.0015;      default is dynamic
+  resource@stMinDistanceF   = 0.008;       streamline dist
+  resource@stLineThicknessF = 1.0;         streamline thickness
+ ; resource@stCrossoverCheckCount = 100;
   if (useScalar) then
     if (useStreamlines) then
       if (showMap)
@@ -1857,6 +2004,25 @@ begin
   end if
   return vc
 end
+;---------------------------------------------------------------
+; limit the map to given coords
+undef("limitMap2Coords")
+procedure limitMap2Coords(resource)
+begin
+  lons = resource@sfXArray
+  lats = resource@sfYArray
+  resource@mpMinLatF         = floor(min(lats))
+  resource@mpMaxLatF         = ceil(max(lats))
+  resource@mpMinLonF         = floor(min(lons))
+  resource@mpMaxLonF         = ceil(max(lons))
+
+
+;   bounds of plotting area for stereographic plots
+  resource@mpLeftCornerLonF  = floor(min(lons))
+  resource@mpLeftCornerLatF  = floor(min(lats))
+  resource@mpRightCornerLonF = ceil(max(lons))
+  resource@mpRightCornerLatF = ceil(max(lats))
+end
 
 ;
 ; vim:ft=ncl

data/test/test_iconPlot.rb

@@ -1,8 +1,8 @@
 $:.unshift File.join(File.dirname(__FILE__),"..","lib")
-require "test/unit"
+require "minitest/autorun"
 require "iconPlot"
 
-class TestIconPlot < Test::Unit::TestCase
+class TestIconPlot < Minitest::Test
 
   CALLER                 = "/home/ram/src/git/icon/scripts/postprocessing/tools/contrib/nclsh"
   PLOTTER                = "/home/ram/src/git/icon/scripts/postprocessing/tools/icon_plot.ncl"

metadata

@@ -1,55 +1,55 @@
 --- !ruby/object:Gem::Specification
 name: iconPlot
 version: !ruby/object:Gem::Version
-  version: 0.0.8
+  version: 0.0.9
 platform: ruby
 authors:
 - Ralf Mueller
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2013-03-26 00:00:00.000000000 Z
+date: 2014-06-17 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: cdo
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: extcsv
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: gnuplot
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
         version: '0'
 description: 'Plot with ncl via Ruby: requires NCL 6.* and CDO 1.5.*'
@@ -59,9 +59,9 @@ executables:
 extensions: []
 extra_rdoc_files: []
 files:
-- lib/iconPlot.rb
-- gemspec
 - bin/nclsh
+- gemspec
+- lib/iconPlot.rb
 - lib/icon_plot.ncl
 - lib/icon_plot_lib.ncl
 - test/test_iconPlot.rb
@@ -75,17 +75,17 @@ require_paths:
 - lib
 required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
-  - - '>='
+  - - ">="
     - !ruby/object:Gem::Version
       version: '1.8'
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
-  - - '>='
+  - - ">="
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 2.0.0
+rubygems_version: 2.2.2
 signing_key: 
 specification_version: 4
 summary: Plot ICON output with ncl via Ruby