iconPlot 0.0.6 → 0.0.8

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 81daefe787ba86d077804f2afbb20f24bb8625e0
+  data.tar.gz: e2538b75ce60c906d8c90afd08654bc4a9f7a9f1
+SHA512:
+  metadata.gz: b2a7478b62d8e95335816b006caeb81471f749cf1d014875647a9504a2efe03bca72c5a7db1036650407ff9c5e0a6d59a04865cef19b869de7163bf1213ae5aa
+  data.tar.gz: c7da25dd5cd29dcc0232e049bf86da2192ab56cab2ff232a5b1a2f247d1158d8d1ace670f39c7a664607fe5382b6628a934867a4f0cb3ce9705b23558a2b5ef0

data/gemspec

@@ -3,9 +3,10 @@ $:.unshift File.join(File.dirname(__FILE__),"..","lib")
 
 spec = Gem::Specification.new do |s|
   s.name              = "iconPlot"
-  s.version           = '0.0.6'
+  s.version           = '0.0.8'
   s.platform          = Gem::Platform::RUBY
-  s.files             = ["lib/iconPlot.rb"] + ["gemspec"] + ["contrib/nclsh"] + Dir.glob("lib/*ncl")
+  s.files             = ["lib/iconPlot.rb"] + ["gemspec"] + ["bin/nclsh"] + Dir.glob("lib/*ncl")
+  s.executables       << 'nclsh'
   s.test_file         = "test/test_iconPlot.rb"
   s.description       = "Plot with ncl via Ruby: requires NCL 6.* and CDO 1.5.*"
   s.summary           = "Plot ICON output with ncl via Ruby"

data/lib/iconPlot.rb

@@ -6,19 +6,13 @@ require 'shellwords'
 class IconPlot < Struct.new(:caller,:plotter,:libdir,:otype,:display,:cdo,:debug,:isIcon)
   VERSION = '0.0.4'
 
-  def IconPlot.gemPath
-    gemSearcher = Gem::GemPathSearcher.new
-    gemspec     = gemSearcher.find('iconPlot')
-    gemspec.gem_dir
-  end
   def initialize(*args)
     super(*args)
-
-    gempath = IconPlot.gemPath
+    
     defaults = {
-      :caller  => [gempath,'contrib','nclsh'].join(File::SEPARATOR),
-      :plotter => [gempath,'lib','icon_plot.ncl'].join(File::SEPARATOR),
-      :libdir  => [gempath,'lib'].join(File::SEPARATOR),
+      :caller  => Gem.bin_path('iconPlot','nclsh'),
+      :plotter => Gem.find_files("icon_plot.ncl")[0],
+      :libdir  => File.dirname(Gem.find_files("icon_plot.ncl")[0]),
       :otype   => 'png',
       :display => 'sxiv',
       :cdo     => ENV['CDO'].nil? ? 'cdo' : ENV['CDO'],
@@ -41,7 +35,7 @@ class IconPlot < Struct.new(:caller,:plotter,:libdir,:otype,:display,:cdo,:debug
                  raise ArgumentError
                end
 
-    opts[:tStrg] =ofile
+    opts[:tStrg] =ofile if opts[:tStrg].nil?
 
     cmd   = [self.caller,self.plotter].join(' ')
     cmd << " -altLibDir=#{self.libdir} #{varIdent} -iFile=#{ifile} -oFile=#{ofile} -oType=#{self.otype}"
@@ -52,7 +46,8 @@ class IconPlot < Struct.new(:caller,:plotter,:libdir,:otype,:display,:cdo,:debug
     out = IO.popen(cmd).read
     puts out if self.debug
 
-    return [FileUtils.pwd,"#{ofile}.#{self.otype}"].join(File::SEPARATOR)
+    #return [FileUtils.pwd,"#{ofile}.#{self.otype}"].join(File::SEPARATOR)
+    return "#{ofile}.#{self.otype}"
   end
   def scalarPlot(ifile,ofile,varname,opts={})
     plot(ifile,ofile,varname,'scalar',opts)

data/lib/icon_plot.ncl

@@ -9,23 +9,53 @@
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
 
+;------- lib location from the command line ------------------------------------------------------------------
+ libLoaded = False
+ if ( isvar("altLibDir") ) then
+   print("Load library: "+altLibDir+"/icon_plot_lib.ncl")
+   loadscript(altLibDir +"/icon_plot_lib.ncl")
+   libLoaded = True
+ end if
+;------------------------------------------------------------------------------
 ;------- config ---------------------------------------------------------------
   ; default settings from user file
-  CONFIGFILE = getenv("HOME")+"/.icon_plot.rc"
-  if (isfilepresent(CONFIGFILE)) then
-    print("Found configfile:"+CONFIGFILE)
-    loadscript(CONFIGFILE)
+  ; --- Skip any user configuration through a file
+  if (isvar("noConfig")) then
+    if ( isvar("DEBUG") .and. DEBUG) print("Subpress user configuration") end if
+  else
+  ; --- Use the configuration file from the commandline option
+    if ( isvar("configFile") ) then
+      if (isfilepresent(configFile)) then
+        print("Found configfile:"+configFile)
+        loadscript(configFile)
+      else
+        print("Could NOT find configfile:"+configFile)
+        exit()
+      end if
+    else
+  ; --- Use the default config file
+      CONFIGFILE = getenv("HOME")+"/.icon_plot.rc"
+      if (isfilepresent(CONFIGFILE)) then
+        print("Found configfile:"+CONFIGFILE)
+        loadscript(CONFIGFILE)
+      end if
+    end if
   end if
-
+;--- end config ---------------------------------------------------------------
+;--- default lib location  ----------------------------------------------------
   if ( .not. isvar("altLibDir") )  then
     print("load default library from /pool/data/ICON/tools")
     loadscript("/pool/data/ICON/tools/icon_plot_lib.ncl")
   else
-    print("Load library: "+altLibDir+"/icon_plot_lib.ncl")
-    loadscript(altLibDir +"/icon_plot_lib.ncl")
+    if ( .not. libLoaded ) then
+;--- lib location from the config file ----------------------------------------
+      print("Load library: "+altLibDir+"/icon_plot_lib.ncl")
+      loadscript(altLibDir +"/icon_plot_lib.ncl")
+    end if
   end if
-;--- end config ---------------------------------------------------------------
+;------------------------------------------------------------------------------
 ;==============================================================================
+  VERSION = "0.0.6"
 ;----- MAIN PROGRAM -----------------------------------------------------------
   if (isvar("help")) then
     dq = str_get_dq()
@@ -84,23 +114,28 @@
     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(";                 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("; 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:")
     print(";                 halflog (uses half logarythmic levels,i.e. 1,2,5 per decade), ")
     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("; contLevs      : set individual contour levels (not yet)")
+    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("; 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)")
@@ -114,7 +149,10 @@
     print("; bStrg         : base string - default is prgr 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("; 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(";")
@@ -122,11 +160,14 @@
     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")
+    print("; noConfig      : subpress reading of any config file, i.e. use commandline options only (default:False)")
+    print(";                 Use +noConfig or -noConfig=True")
     print("; cdo           : Path to an alternative cdo version")
     print(";---------------------------------------------------------------")
     print("; Authors       : Ralf Mueller (ralf.mueller@zmaw.de)")
     print(";                 Stephan Lorenz (stephan.lorenz@zmaw.de)")
-    print("; VERSION       : r8012 (2012-01-26)")
+    print(";                 Einar Olason (einar.olason@zmaw.de)")
+    print("; VERSION       : " + VERSION +" (2012-08-13)")
     print(";---------------------------------------------------------------")
     exit()
   end if
@@ -157,29 +198,34 @@
   end if
   print("Outputfile '"+str_concat((/oFile,".",oType/))+"' will be created in "+systemfunc("dirname "+iFile))
 
-  if(.not. isvar("numLevs")) then
-    numLevs  = NUMLEVS
+  if (.not. isvar("DEBUG"))          DEBUG          = False           end if
+  if (.not. isvar("showGrid"))       showGrid       = False           end if
+  if (.not. isvar("markCells"))      markCells      = False           end if
+  if (.not. isvar("showNdc"))        showNdc        = False           end if
+  if (.not. isvar("maxView"))        maxView        = False           end if
+  if (.not. isvar("scaleLimit"))     scaleLimit     = 2               end if
+  if (.not. isvar("selMode"))        selMode        = "auto"          end if
+  if (.not. isvar("plotMode"))       plotMode       = "scalar"        end if
+  if (.not. isvar("timeStep")) then
+    timeStep       = 0
   else
-    if (numLevs .lt. 1) then
-      print("WARNING: numLevs must be >= 1")
-      print(ABORTMSG)
-      exit
-    end if
+    t = floattointeger(timeStep)
+    delete(timeStep)
+    timeStep = t
   end if
-
-  if (.not. isvar("DEBUG"))       DEBUG       = False           end if
-  if (.not. isvar("showGrid"))    showGrid    = False           end if
-  if (.not. isvar("showNdc"))     showNdc     = False           end if
-  if (.not. isvar("maxView"))     maxView     = False           end if
-  if (.not. isvar("scaleLimit"))  scaleLimit  = 2               end if
-  if (.not. isvar("selMode"))     selMode     = "auto"          end if
-  if (.not. isvar("plotMode"))    plotMode   = "scalar"         end if
-  if (.not. isvar("timeStep"))    timeStep    = 0               end if
-  if (.not. isvar("levIndex"))    levIndex    = 0               end if
-  if (.not. isvar("scaleFactor")) scaleFactor = 1               end if
-  if (.not. isvar("colormap"))    colormap    = "BlueDarkRed18" end if
-  if (.not. isvar("k2c"))         k2c         = True            end if
-  if (.not. isvar("isIcon"))      isIcon      = False           end if
+  if (.not. isvar("levIndex"))       levIndex       = 0               end if
+  if (.not. isvar("scaleFactor"))    scaleFactor    = 1               end if
+  if (.not. isvar("fillMode"))       fillMode       = "RasterFill"    end if
+  if (.not. isvar("colormap")) then
+    colormap   = DEFAULTCOLOR
+    userColors = False
+  else
+    userColors = True
+  end if
+  if (.not. isvar("withLines"))      withLines      = True            end if
+  if (.not. isvar("withLineLabels")) withLineLabels = False           end if
+  if (.not. isvar("k2c"))            k2c            = True            end if
+  if (.not. isvar("isIcon"))         isIcon         = False           end if
 
   if (.not.isvar("mapType"))      mapType      = "lonlat"        end if
   if ( mapType .eq. "NHps" .and. .not.isvar("mapLLC")) mapLLC = (/-45., 45./)  end if
@@ -198,14 +244,15 @@
   if (.not.isvar("vecColByLen"))  vecColByLen  = False           end if
   if (.not.isvar("vecMinDist"))   vecMinDist   = 0.017           end if
   if (.not.isvar("streamLine"))   streamLine   = False           end if
-  if (.not.isvar("secPoints"))    secPoints    = 20              end if
-  ; default + optional coordinates
-  if (.not.isvar("lonCo"))         lonCo       =  DEFAULTLON     end if
-  if (.not.isvar("latCo"))         latCo       =  DEFAULTLAT     end if
-
-  if (.not.isvar("showSecMap")) showSecMap     =  True           end if
-  ; ATMOSPEHRE SETUP ----------------------------------------------------------
-  if (.not.isvar("atmLev"))     atmLev         =  "p"           end if
+  if (.not.isvar("secPoints"))    secPoints    = 161             end if
+;   default + optional coordinates
+  if (.not.isvar("lonCo"))         lonCo       = DEFAULTLON     end if
+  if (.not.isvar("latCo"))         latCo       = DEFAULTLAT     end if
+
+  if (.not.isvar("secMode"))    secMode        = "straight"     end if
+  if (.not.isvar("showSecMap")) showSecMap     = True           end if
+;   ATMOSPEHRE SETUP ----------------------------------------------------------
+  if (.not.isvar("atmLev"))     atmLev         = "p"           end if
   if (.not.isvar("atmPLevs")) then
     atmPLevSetup = (/10000,100000,2500/); unit: Pa
     atmPLevs     = ispan(atmPLevSetup(0),atmPLevSetup(1),atmPLevSetup(2))*1.0
@@ -280,6 +327,28 @@
   horizonalGridType = getHorizGridType(iFile,ivarname,isIcon)
   ;----------------------------------------------------------------------------
 
+  ; DETERMINE selMode automatically -------------------------------------------
+  if (isvar("minVar") .and. isvar("maxVar"))  then
+    if (selMode .eq. "auto") then
+      selMode = "manual"
+    end if
+  end if
+  if (isvar("plotLevs") )  then
+    if (selMode .eq. "auto") then
+      selMode = "manual"
+    end if
+  end if
+  if(.not. isvar("numLevs")) then
+    numLevs  = NUMLEVS
+  else
+    if (numLevs .lt. 1) then
+      print("WARNING: numLevs must be >= 1")
+      print(ABORTMSG)
+      exit
+    end if
+  end if
+  ;----------------------------------------------------------------------------
+
   ; DETERMINE THE PLOTMODE ----------------------------------------------------
   if (isvar("vecVars")) then
     ; expectes is a string or an array of size 2: "u-veloc v-veloc" or (/"u-veloc","v-veloc"/)
@@ -310,16 +379,16 @@
     exit
   end if
   if (isvar("secLC").and.isvar("secRC")) plotMode = "section" end if
-  if (isvar("hoff")) plotMode = "hoffmueller" end if
+  if (isvar("hov")) plotMode = "hovmoeller" end if
 
-  if (plotMode.eq."vector" .or. plotMode.eq."overlay" .or. plotMode.eq."section" .or. plotMode .eq. "hoffmueller") then
+  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
       ; from unstructured grids
       ; TODO addVars=(/ "PS","PHIS"/)
-      remapFilename = setRemapFilename(iFile,iType,atmLev)
+      remapFilename = setRemapFilename(iFile,iType,resolution,atmLev)
       print("remapFilename:"+remapFilename)
       if (.not. checkRemappedFile(iFile,remapFilename,ivarname,iType,atmLev,atmPLevs,atmHLevs) ) then
         addVars=(/""/)
@@ -328,10 +397,10 @@
       end if
       rFile = addfile( remapFilename+".nc", "r" )
     end if
-    if (plotMode .eq. "hoffmueller") then ; perform zonmean
-      zonmeanFilename = setZonmeanFilename(remapFilename,ivarname,atmLev)
-      zonmean4HoffmuellerPlot(remapFilename,ivarname,zonmeanFilename)
-    end if
+   ;if (plotMode .eq. "hoffmueller") then ; perform zonmean
+   ;  zonmeanFilename = setZonmeanFilename(remapFilename,ivarname,atmLev)
+   ;  zonmean4HoffmuellerPlot(remapFilename,ivarname,zonmeanFilename)
+   ;end if
   end if
   File  = addfile( iFile+".nc", "r" )
 
@@ -369,7 +438,7 @@
     print("Graphics format is " +oType)
     print("plotMode    = "+plotMode)
 
-    if (plotMode.eq."scalar" .or. plotMode.eq."section") then
+    if (plotMode.eq."scalar" .or. plotMode.eq."section" .or. plotMode.eq."hovmoeller") then
       print("varName     = "+varName)
     else
       print("vecVars     = "+vecVars)
@@ -416,12 +485,6 @@
 
     scaleVar(var,scaleFactor)
 
-    ; set selMode automatically
-    if (isvar("minVar") .and. isvar("maxVar"))  then
-      if (selMode .eq. "auto") then
-        selMode = "manual"
-      end if
-    end if
     ; set minVar_maxVar for plotting
     if(.not. isvar("minVar")) minVar = min(var) end if
     if(.not. isvar("maxVar")) maxVar = max(var) end if
@@ -435,6 +498,11 @@
       ;var = mask ( var, slm, 0.5)
     end if
   end if
+  if (plotMode.eq."hovmoeller") then
+    var = selRegularVar(varName,File,-1)
+    if(.not. isvar("minVar")) minVar = min(var) end if
+    if(.not. isvar("maxVar")) maxVar = max(var) end if
+  end if
   if (plotMode.eq."vector" .or. plotMode.eq."overlay") then; vector mode =======
     uvarname = vecVars(0)
     vvarname = vecVars(1)
@@ -488,11 +556,6 @@
       var = mask(var,maskVar,1)
     end if
   end if ; Reading data variables =============================================
-  if (plotMode .eq. "hoffmueller") then ; create hoffmuelle diagram
-    ; select the whole time series
-    ; and plot it with
-    ; plot = gsn_csm_lat_time(wks, shov(0,{-80:80},:), res)
-  end if
   if (iType .eq. "atm" .and. ivarname .eq. "T" .and. k2c) var = var - 273.15 end if ;atmosphere used to write out Kelvin
 
   ;---------------------------------------------------------------
@@ -507,20 +570,18 @@
     drawNDCGrid(wks)
   end if
 
-  ResC = True
+  ResC = setDefaultResource(True,withLines,withLineLabels,fillMode)
 
-  setDefaultResource(ResC,True)
+  if (.not. userColors) setDefaultColors(ResC,minVar,maxVar,DEBUG) end if
 
   if (maxView) ResC@gsnMaximize = True end if
 
   if (useMask) setMaskColor(wks,ResC) end if
 
 
-
   if (getHorizGridType(iFile,ivarname,isIcon) .eq. "unstructured") then
     if (plotMode.eq."scalar" .or. plotMode.eq."overlay")
       setCoordinates(ResC,x,y)
-      ;setBounds(ResC,File,x,y,DEBUG)
     end if
   end if
 
@@ -535,7 +596,7 @@
       print("Cannot plot vertical section of a 2D variable: "+varName)
       exit
     end if
-    trans           = setSection(secLC,secRC,secPoints,var)
+    trans           = setSection(secLC,secRC,secPoints,var,secMode)
     vertdim         = getVertDim(File,var)
     trans!0         = vertdim
     trans&$vertdim$ = var&$vertdim$
@@ -564,7 +625,11 @@
     else
       setDefaultVectorPlot(ResC,5.0,vecRefLength,"CurlyVector",vecMinDist)
 
-      vc = plotVecOrStream(vecColByLen,streamLine,uvar,vvar,velocity,ResC,wks,True)
+      vc = plotVecOrStream(vecColByLen,streamLine,\
+                          uvar,vvar, \
+                          sqrt(uvar*uvar + vvar*vvar), \
+                          ResC,wks,True)
+      draw(vc)
     end if
   end if
   if (plotMode.eq."overlay") then
@@ -577,7 +642,7 @@
 
     ResC@gsnDraw             = False
     ResC@gsnFrame            = False
-    vc   = plotVecOrStream(vecColByLen,streamLine,uvar,vvar,velocity,ResC2,wks,False)
+    vc   = plotVecOrStream(vecColByLen,streamLine,uvar,vvar,var,ResC2,wks,False)
 
     plot = gsn_csm_contour_map(wks,var,ResC)
     if (showGrid) then
@@ -600,10 +665,19 @@
       end if
       plotGrid(wks,plot,var,x,bounds,File,ResC,DEBUG)
     end if
+    if ( markCells ) then
+      plotCellMarkers(x,y,wks,plot,True)
+    end if
   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)
+  end if
   if (plotMode.eq."section") then
 
     if (showSecMap) then
@@ -627,20 +701,33 @@
 
    if (showSecMap) then
      ; map with section polygon
-     mres                  = True
-     mres@gsnFrame         = False;                don't turn page yet
-     mres@gsnDraw          = False;                don't draw yet
-     mres@vpWidthF         = 0.8;          set width of plot
-     mres@vpHeightF        = 0.3;          set height of plot
-     mres@vpXF             = 0.1
-     mres@vpYF             = 0.4
-     plot(1)               = gsn_csm_map_ce(wks,mres)
+     mres              = True
+     mres@gsnFrame     = False;                don't turn page yet
+     mres@gsnDraw      = False;                don't draw yet
+     mres@vpWidthF     = 0.8;          set width of plot
+     mres@vpHeightF    = 0.3;          set height of plot
+     mres@vpXF         = 0.1
+     mres@vpYF         = 0.4
+
+     newcolor          = NhlNewColor(wks,0.85,0.85,0.85);    add gray to colormap
+     mres@mpOutlineOn  = True;  (turn off continental outlines)
+     mres@mpFillOn     = False;  (turn on map fill)
+     mres@mpFillColors = (/"background","transparent","LightGray","transparent"/);  (fill land in gray and ocean in transparent)
+     mres@mpGrid       = False
+
+ ;    plot(1)               = gsn_map(wks,"CylindricalEquidistant",mres)
+     plot(1)               = gsn_csm_map(wks,mres)
      pres                  = True;            polyline mods desired
      pres@gsnFrame         = False;                don't turn page yet
      pres@gsnDraw          = False;                don't draw yet
      pres@gsLineColor      = "black";           color of lines
      pres@gsLineThicknessF = 2.0;             line thickness
-     dummy = gsn_add_polyline(wks,plot(1),(/ secLC(0),secRC(0) /),(/ secLC(1),secRC(1) /),pres)
+     pathRes               = getSectionPathForMap(secLC,secRC,secPoints,secMode)
+     if (DEBUG) then
+       print("PATHRES lon - lat")
+       print(pathRes@lon + "      " + pathRes@lat)
+     end if
+     pathLine              = gsn_add_polyline(wks,plot(1),pathRes@lon, pathRes@lat,pres)
 
      draw(plot(1))
    end if

data/lib/icon_plot_lib.ncl

@@ -16,14 +16,18 @@
 ;---------------------------------------------------------------
 ; Authors       : Ralf Mueller (ralf.mueller@zmaw.de)
 ;                 Stephan Lorenz (stephan.lorenz@zmaw.de)
+; VERSION = "0.0.6"
 ;-------------------------------------------------------------------------------
 ; some global defaults
-  DEFAULTLON = "clon"
-  DEFAULTLAT = "clat"
-  RAD2DEG    = 45./atan(1.)
-  ABORTMSG   = "Abort Script."
-  NUMLEVS    = 10
-  CDO        = "cdo"
+  DEFAULTLON   = "clon"
+  DEFAULTLAT   = "clat"
+  RAD2DEG      = 45./atan(1.)
+  ABORTMSG     = "Abort Script."
+  NUMLEVS      = 10
+  CDO          = "cdo"
+;  DEFAULTCOLOR = "BlWhRe"             ; ncolors = 103
+  DEFAULTCOLOR = "BlueDarkRed18"      ; ncolors = 256
+  DEFAULTCOLOR = "BlueWhiteOrangeRed" ; ncolors = 256
 ;-------------------------------------------------------------------------------
 procedure setCDO(path)
 begin
@@ -233,17 +237,19 @@ end
 undef("selRegularVar")
 function selRegularVar(varname,filehandle,timestep)
 begin
-    dims      = getfilevardims(filehandle,varname)
-    noOfDims  = dimsizes(dims)
-    noOfDims  = noOfDims(0)
-    usableVar = False
-    if ( has_time(dims) ) then
+  dims      = getfilevardims(filehandle,varname)
+  noOfDims  = dimsizes(dims)
+  noOfDims  = noOfDims(0)
+  print(""+noOfDims)
+  usableVar = False
+  if ( has_time(dims) ) then
+    if (timestep.eq.-1) then
       if (noOfDims .eq. 4) then        ; dims: (time,lev,x,y)
-        var = filehandle->$varname$(timestep,:,:,:)
+        var = filehandle->$varname$
         usableVar = True
       end if
       if (noOfDims .eq. 3) then        ; dims: (time,x,y)
-        var = filehandle->$varname$(timestep,:,:)
+        var = filehandle->$varname$(:,:,:)
         usableVar = True
       end if
       if (noOfDims .eq. 2) then        ; dims: (x,y)
@@ -255,24 +261,42 @@ begin
         usableVar = False
       end if
     else
-      if (noOfDims .eq. 3) then        ; dims: (lev,x,y)
-        var = filehandle->$varname$(:,:,:)
+      if (noOfDims .eq. 4) then        ; dims: (time,lev,x,y)
+        var = filehandle->$varname$(timestep,:,:,:)
+        usableVar = True
+      end if
+      if (noOfDims .eq. 3) then        ; dims: (time,x,y)
+        var = filehandle->$varname$(timestep,:,:)
         usableVar = True
       end if
       if (noOfDims .eq. 2) then        ; dims: (x,y)
         var = filehandle->$varname$(:,:)
         usableVar = True
       end if
-      if (noOfDims .eq. 1) then        ; dims: (x)
+      if (noOfDims .eq. 1) then        ; dims: (time)
         var = filehandle->$varname$(:)
-        usableVar = True
+        usableVar = False
       end if
     end if
-    if (.not. usableVar) then
-      print("Warning! Could not read " +  varname + " from input file")
-      exit
+  else
+    if (noOfDims .eq. 3) then        ; dims: (lev,x,y)
+      var = filehandle->$varname$(:,:,:)
+      usableVar = True
     end if
-    return var
+    if (noOfDims .eq. 2) then        ; dims: (x,y)
+      var = filehandle->$varname$(:,:)
+      usableVar = True
+    end if
+    if (noOfDims .eq. 1) then        ; dims: (x)
+      var = filehandle->$varname$(:)
+      usableVar = True
+    end if
+  end if
+  if (.not. usableVar) then
+    print("Warning! Could not read " +  varname + " from lonlat input file")
+    exit
+  end if
+  return var
 end
 ;-------------------------------------------------------------------------------
 ; abstract Field selection
@@ -566,9 +590,9 @@ end
 ;---------------------------------------------------------------
 ; set filename for the automatically remapped file
 undef("setRemapFilename")
-function setRemapFilename(filename,itype,atmlev)
+function setRemapFilename(filename,itype,resolution,atmlev)
 begin
-  return setNewFilename(filename,"remapnn",itype,atmlev)
+  return setNewFilename(filename,"remapnn_"+resolution,itype,atmlev)
 end
 undef("setZonmeanFilename")
 function setZonmeanFilename(filename,itype,atmlev)
@@ -778,17 +802,24 @@ undef("setLevels")
 procedure setLevels(selmode,resource,minvar,maxvar,scalelimit,numlevs,debug)
 begin
   if (selmode .eq. "manual") then
-    resource@cnLevelSelectionMode = "ManualLevels"
-    resource@cnMinLevelValF       = minvar
-    resource@cnMaxLevelValF       = maxvar
-;    diffLog10                 = log10(abs(maxVar-minVar))
-;;   if (diffLog10 .lt. 0) then
-;      resource@cnLevelSpacingF    = 10^(floor(diffLog10))/numLevs
-;;   else
-;;     resource@cnLevelSpacingF    = 10^(floor(diffLog10))/numLevs
-;;   end if
-    diffspacing               = abs(maxvar-minvar)/(int2flt(numlevs))
-    resource@cnLevelSpacingF      = diffspacing
+    if ( .not. isvar("plotLevs") ) then
+      resource@cnLevelSelectionMode = "ManualLevels"
+      resource@cnMinLevelValF       = minvar
+      resource@cnMaxLevelValF       = maxvar
+  ;    diffLog10                 = log10(abs(maxVar-minVar))
+  ;;   if (diffLog10 .lt. 0) then
+  ;      resource@cnLevelSpacingF    = 10^(floor(diffLog10))/numLevs
+  ;;   else
+  ;;     resource@cnLevelSpacingF    = 10^(floor(diffLog10))/numLevs
+  ;;   end if
+      diffspacing               = abs(maxvar-minvar)/(int2flt(numlevs))
+      resource@cnLevelSpacingF      = diffspacing
+    else
+      resource@cnLevelSelectionMode = "ExplicitLevels"
+      plotLevels = plotLevs
+      if (debug) print("plotlevels = "+plotLevels) end if
+      resource@cnLevels             = plotLevels
+    end if
   end if
 
   if (selmode .eq. "halflog") then
@@ -810,17 +841,39 @@ end
 ;---------------------------------------------------------------
 ; Create a default NCL resource
 undef("setDefaultResource")
-procedure setDefaultResource(resource,verticallabelbar)
-begin
-  resource@gsnMaximize          = False
-; resource@gsnPaperOrientation  = "landscape"  ; turns plot on paper
-; resource@wkPaperSize          = "A4"
-  resource@gsnFrame             = False
-  resource@gsnSpreadColors      = True
+function setDefaultResource(verticallabelbar,withLines,withLineLabels,fillmode)
+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@cnLineLabelsOn           = False
+  if (withLineLabels ) then
+    withLines                       = True
+    resource@cnLineLabelsOn         = True
+    resource@cnLineLabelFontHeightF = 0.008
+    resource@cnLineLabelFontThicknessF = 0.006
+  end if
+  if (withLines) resource@cnLinesOn = True end if
 
-  resource@cnFillOn             = True   ;  writes labelbar, no contour labels
-  resource@cnLinesOn            = False
-  resource@cnInfoLabelOn        = False
+  ; representation of missing values ===========================================
+  delete(resource@cnMissingValFillPattern);     no longer filling missing value areas
+  delete(resource@cnMissingValFillColor);       no longer filling missing value areas
+  resource@cnMissingValPerimOn          = True
+  resource@cnMissingValFillPattern      = -1;                set the missing value fill pattern to 5
+  resource@cnMissingValFillScaleF       = 0.9;              increase the density of the fill pattern (default   = 1.0)
+  resource@cnMissingValPerimColor       = "black";          change the missing value perimeter to black
+  resource@cnMissingValPerimDashPattern = 1;           set the dash pattern of the missing value perimeter to 1
+  resource@cnMissingValPerimThicknessF  = 1.0;         increase the thickness of the missing value perimeter 3X
 
   FontHeight0               = 0.012
   FontHeight1               = 0.015
@@ -835,8 +888,7 @@ begin
   resource@gsnStringFontHeightF = FontHeight0
 
   resource@mpFillOn               = True
-  resource@cnFillOn               = True
-
+; label settings ==============================================================
   resource@lbLabelBarOn             = True
   if (verticallabelbar .eq. True) then
     resource@lbOrientation            = "vertical"
@@ -848,11 +900,12 @@ begin
   resource@lbLabelAutoStride        = True
   resource@lbLabelFontHeightF       = FontHeight1     ;   color bar labels
 
-  resource@cnFillMode           = "RasterFill"
-  ;resource@cnFillMode           = "CellFill"
+  resource@cnFillMode           = fillmode
   resource@cnRasterSmoothingOn  = True
   resource@mpGreatCircleLinesOn = True
   resource@stMinArrowSpacingF   = 0.001
+
+  return resource
 end
 ;---------------------------------------------------------------
 undef("setDefaultSectionResource")
@@ -862,20 +915,32 @@ begin
   resource@gsnFrame                     = False;                don't turn page yet
   resource@gsnDraw                      = False;                don't draw yet
   resource@tmXBMode                     = "Explicit";     explicitly label x-axis
+  latDiff = stringtodouble(secrc(1)) - stringtodouble(seclc(1))
+  tic = latDiff/(stringtodouble(secpoints)-1)
+  p = (points * tic) + stringtodouble(seclc(1))
+ ; resource@tmXBValues                   = points
   resource@tmXBValues                   = (/points(0),points(secpoints-1)/);  points to label values
+  labels = new(dimsizes(p),string)
+  labels = p
+;print(""+labels)
+ ;resource@tmXBLabels                   = labels
   resource@tmXBLabels                   = (/ seclc(1) +"N, "+ seclc(0)+"E" , secrc(1)+"N, "+secrc(0)+"E" /)
 
+  resource@tmXBLabelFontHeightF  = 0.01
   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"
 ;   resource@gsnYAxisIrregular2Log      = True
 
-  delete(resource@cnMissingValFillPattern);     no longer filling missing value areas
-  delete(resource@cnMissingValFillColor);       no longer filling missing value areas
+;  delete(resource@cnMissingValFillPattern);     no longer filling missing value areas
+;  delete(resource@cnMissingValFillColor);       no longer filling missing value areas
   resource@cnMissingValFillColor        = "gray30"
   resource@cnMissingValPerimOn          = True
 
@@ -886,8 +951,11 @@ begin
   resource@cnMissingValPerimDashPattern = 1;           set the dash pattern of the missing value perimeter to 1
   resource@cnMissingValPerimThicknessF  = 3.0;         increase the thickness of the missing value perimeter 3X
   ;resource@gsnYAxisIrregular2Log        = True
+  resource@cnFillMode           = "RasterFill"
+  resource@cnRasterSmoothingOn  = False
   return resource
 end
+;---------------------------------------------------------------
 undef("setDefaultOverlayResource")
 procedure setDefaultOverlayResource(resource)
 begin
@@ -905,6 +973,30 @@ begin
   resource@lbLabelBarOn       = False;  switch of the vector label bar, because the speed is show as vector lenghts
 end
 ;---------------------------------------------------------------
+undef("setDefaultColors")
+procedure setDefaultColors(resource,minval,maxval,debug)
+begin
+; white is color 135
+  if ((minval*maxval) .LT. 0.0) then ; data goes fro positiv to negative or vs.
+    if (abs(maxval) .GT. abs(minval)) then
+      resource@gsnSpreadColorStart = 135  - toint(126*(abs(minval)/abs(maxval)))
+      resource@gsnSpreadColorEnd   = 255
+    else
+      resource@gsnSpreadColorStart = 2
+      resource@gsnSpreadColorEnd   = 135 + toint(126*(abs(maxval)/abs(minval)))
+    end if
+    resource@gsnContourZeroLineThicknessF = 2.0
+    if (debug) then
+      print("===================== COMPUTE the COLORS to make zero WHITE:")
+      print("minval: "+minval)
+      print("maxval: "+maxval)
+      print("100*(toint(abs(minval)/abs(maxval))): "+toint(100*(abs(minval)/abs(maxval))))
+      print("resource@gsnSpreadColorStart:" + resource@gsnSpreadColorStart)
+      print("resource@gsnSpreadColorEnd  :" + resource@gsnSpreadColorEnd)
+    end if
+  end if
+end
+;---------------------------------------------------------------
 undef("shift4SecMap")
 procedure shift4SecMap(resource)
 begin
@@ -1031,6 +1123,7 @@ begin
   else
     resource@lbTitleString      = ""
   end if
+  resource@tiMainFontHeightF  = 0.02
   if (itype .eq. "atm" .and. varname .eq. "T" .and. k2c) resource@lbTitleString = "C" end if
 
   resource@lbTitlePosition    = "Bottom"
@@ -1038,6 +1131,8 @@ begin
   resource@lbTitleFontHeightF = 0.02
   resource@lbLabelFontHeightF = 0.015
   resource@lbLeftMarginF      = 0.01
+  resource@cnLabelMasking   = True
+  resource@cnLabelMasking   = True
 end
 ;---------------------------------------------------------------
 ; vertical axes label for section plot
@@ -1183,10 +1278,9 @@ end
 undef("setMaskColor")
 procedure setMaskColor(wks,resource)
 begin
-  newcolor = NhlNewColor(wks,0.85,0.85,0.85)   ; add gray to colormap
-  ;resource@gsnSpreadColorStart  = 1
-  resource@gsnSpreadColorEnd     = 19  ; last color number not known?
-  resource@cnMissingValFillColor = 20  ; last color is used for missing values
+  resource@cnMissingValFillColor       = "gray90"
+  resource@cnMissingValPerimOn         = True;          turn on missing value outlines
+  resource@cnMissingValPerimThicknessF = 2.5;   set the missing value outline thickness
 end
 ;---------------------------------------------------------------
 ; set the coordinates for an ICON data field
@@ -1298,31 +1392,92 @@ begin
   return gridtype
 end
 ;---------------------------------------------------------------
+; compute the straight line between two corners
+undef("getStraightPath")
+function getStraightPath(startCorner,endCorner,npoints)
+begin
+  startLat = tofloat(startCorner(1))
+  endLat   = tofloat(endCorner(1))
+  startLon = tofloat(startCorner(0))
+  endLon   = tofloat(endCorner(0))
+
+  if (startLon .lt. 0.0) then
+    startLon = startLon + 360.0
+    endLon   = endLon   + 360.0
+  end if
+
+  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
+  ; the data. Otherwise an interpolation error will occur
+  lons = where(lons.gt.360, lons-360.0, lons)
+
+  distRes     = True
+  distRes@lat = lats
+  distRes@lon = lons
+
+  return distRes
+end
+;---------------------------------------------------------------
+; compute the distance between two points of great circle
+undef("getCirclePath")
+function getCirclePath(startCorner,endCorner,npoints)
+begin
+  startLat = startCorner(1)
+  endLat   = endCorner(1)
+
+  startLon = startCorner(0)
+  endLon   = endCorner(0)
+
+  distRes     = True
+  dist        = gc_latlon(startLat,startLon,endLat,endLon,npoints,2)
+  distRes@lat = dist@gclat
+  distRes@lon = dist@gclon
+
+  return distRes
+end
+;---------------------------------------------------------------
+; compute the path of the cross-section
+undef("getSectionPath")
+function getSectionPath(seclc,secrc,secpoints,secmode)
+begin
+  if (secmode .eq. "straight") then
+    pathRes = getStraightPath(seclc,secrc,secpoints)
+  else
+    pathRes = getCirclePath(seclc,secrc,secpoints)
+  end if
+  return pathRes
+end
+;---------------------------------------------------------------
+; compute the path of the cross-section but suitable for plotting in a separate Map
+undef("getSectionPathForMap")
+function getSectionPathForMap(seclc,secrc,secpoints,secmode)
+begin
+  pathRes     = getSectionPath(seclc,secrc,secpoints,secmode)
+  pathRes@lon = where(pathRes@lon.lt.pathRes@lon(0),pathRes@lon + 360.0,pathRes@lon)
+  return pathRes
+end
+;---------------------------------------------------------------
 ; setting for vertical cross sections based on remapped icon data
 undef("setSection")
-function setSection(secLC,secRC,npoints,var)
+function setSection(seclc,secrc,secpoints,var,secmode)
 begin
-  leftlat  = secLC(1)
-  rightlat = secRC(1)
+  pathRes = getSectionPath(seclc,secrc,secpoints,secmode)
 
-  leftlon  = secLC(0)
-  rightlon = secRC(0)
-
-  dist     = gc_latlon(leftlat,leftlon,rightlat,rightlon,npoints,2)
-  trans    = linint2_points(var&lon,var&lat,var,True,dist@gclon,dist@gclat,2)
+  trans   = linint2_points(var&lon,var&lat,var,True,pathRes@lon,pathRes@lat,2)
 
   return trans
 end
 ;---------------------------------------------------------------
 ; setting for vertical cross sections based on remapped icon data
 undef("setSectionFromHybridPress")
-function setSectionFromHybridPress(secLC,secRC,npoints,var,filehandle,timestep)
+function setSectionFromHybridPress(seclc,secrc,npoints,var,filehandle,timestep)
 begin
-  leftlat  = secLC(1)
-  rightlat = secRC(1)
+  leftlat  = seclc(1)
+  rightlat = secrc(1)
 
-  leftlon  = secLC(0)
-  rightlon = secRC(0)
+  leftlon  = seclc(0)
+  rightlon = secrc(0)
 
   ;config
   pressName  = "PS"
@@ -1557,8 +1712,8 @@ begin
 end
 ;---------------------------------------------------------------
 ; Create a resource for the grid plot
-undef("setupGridResource")
-function setupGridResource(var,levels,colors,boundslon,boundslat,wks,plot,debug)
+undef("performGridPlot")
+procedure performGridPlot(var,levels,colors,boundslon,boundslat,wks,plot,flags,debug)
 begin
   pres             = True
   pres@gsEdgesOn   = True   ; Turn on edges
@@ -1571,7 +1726,9 @@ begin
   if (.not. ismissing(vlow(0))) then
     do j = 0, dimsizes(vlow)-1
       pres@gsFillColor = colors(i)   ; first color
+      if ( flags(vlow(i))) then
       gsn_polygon(wks,plot,boundslon(vlow(j),:),boundslat(vlow(j),:),pres)
+      end if
     end do
   end if
   if (debug) then
@@ -1587,7 +1744,9 @@ begin
     if (.not. ismissing(vind(0))) then
     do j = 0, dimsizes(vind)-1
       pres@gsFillColor = colors(i+1)
-                       gsn_polygon( wks,plot,   boundslon(vind(j),:),boundslat(vind(j),:),pres)
+      ;if ( flags(vind(i))) then
+      gsn_polygon( wks,plot,   boundslon(vind(j),:),boundslat(vind(j),:),pres)
+      ;end if
     end do
     end if
     if (debug) then
@@ -1603,18 +1762,30 @@ begin
 
   ; if no index vhig is found with values larger levels(max-1), vhig is missing value:
   if (.not. ismissing(vhig(0))) then
-  do j = 0, dimsizes(vhig) -1
-    pres@gsFillColor = colors(i+1)
-                     gsn_polygon( wks,plot,   boundslon(vhig(j),:),boundslat(vhig(j),:),pres)
-  end do
+    do j = 0, dimsizes(vhig) -1
+      pres@gsFillColor = colors(i+1)
+      if ( flags(vhig(i))) then
+      gsn_polygon(wks,plot,boundslon(vhig(j),:),boundslat(vhig(j),:),pres)
+      end if
+    end do
   end if
   if (debug) then
     print ("i=" + i + ", Values   larger than      "+levels(i) +": " \
            + dimsizes(vhig) + " triangles considered, color= " + colors(i+1))
   end if
   delete(vhig)
-
-  return pres
+end
+;---------------------------------------------------------------
+; mark the cell location
+undef("plotCellMarkers")
+procedure plotCellMarkers(cell_lons,cell_lats,wks,plot,debug)
+begin
+  pres             = True
+  pres@gsMarkerIndex = 1
+  pres@gsMarkerSizeF = 0.003
+  do i = 0,dimsizes(cell_lons)-1
+    gsn_polymarker(wks,plot,cell_lons(i),cell_lats(i),pres)
+  end do
 end
 ;---------------------------------------------------------------
 ; Plot the ICON grid of a variable to the plot.
@@ -1637,8 +1808,7 @@ begin
   levels = getLevelsFromPlot(plot)
 
   if (debug) printGridPlotInfo(colors,levels,flags) end if
-
-  pres = setupGridResource(variable,levels,colors,boundslon,boundslat,workstation,plot,debug)
+  performGridPlot(variable,levels,colors,boundslon,boundslat,workstation,plot,flags,debug)
 end
 ;---------------------------------------------------------------
 ; read in mask variable
@@ -1651,6 +1821,11 @@ end
 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
   if (useScalar) then
     if (useStreamlines) then
       if (showMap)
@@ -1684,4 +1859,4 @@ begin
 end
 
 ;
-; vim:list
+; vim:ft=ncl

data/test/test_iconPlot.rb

@@ -20,15 +20,17 @@ class TestIconPlot < Test::Unit::TestCase
   if 'thingol' == `hostname`.chomp
     def test_simple
       ip = IconPlot.new(CALLER,PLOTTER,PLOTLIB,OFMT,PLOT_CMD,CDO,true)
+      ip.debug = true
       ofile          = 'test_icon_plot'
-      ip.show(ip.scalarPlot(OCE_PLOT_TEST_FILE,   ofile,"T",:levIndex => 0))
-      ip.show(ip.scalarPlot(OCE_PLOT_TEST_FILE,   ofile,"T",:levIndex => 2))
-      ip.show(ip.vectorPlot(OCE_PLOT_TEST_FILE,   ofile,"u-veloc v-veloc",:levIndex => 2))
-      ip.show(ip.scalarPlot(OCE_PLOT_TEST_FILE,   ofile,"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:mapType => "ortho"))
-      ip.show(ip.scalarPlot(OCE_PLOT_TEST_FILE,   ofile,"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:secLC => "-20,-60", :secRC => "-20,60"))
-      ip.show(ip.scalarPlot(OCELSM_PLOT_TEST_FILE,ofile,"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:secLC => "-20,-60", :secRC => "-20,60",:maskName => 'wet_c'))
-      ip.show(ip.scalarPlot(OCELSM_PLOT_TEST_FILE,ofile,"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:maskName => 'wet_c'))
-      ip.show(ip.scalarPlot(OCELSM_PLOT_TEST_FILE,ofile,"T",:levIndex => 2,:maskName => 'wet_c'))
+      ip.show(ip.scalarPlot(OCE_PLOT_TEST_FILE,   ofile+'_00',"T",:levIndex => 0))
+      ip.show(ip.scalarPlot(OCE_PLOT_TEST_FILE,   ofile+'_01',"T",:levIndex => 2))
+      ip.show(ip.vectorPlot(OCE_PLOT_TEST_FILE,   ofile+'_02',"u-veloc v-veloc",:levIndex => 2))
+      ip.show(ip.scalarPlot(OCE_PLOT_TEST_FILE,   ofile+'_03',"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:mapType => "ortho"))
+      ip.show(ip.scalarPlot(OCE_PLOT_TEST_FILE,   ofile+'_04',"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:secLC => "-20,-60", :secRC => "-20,60"))
+      ip.show(ip.scalarPlot(OCELSM_PLOT_TEST_FILE,ofile+'_05',"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:secLC => "-20,-60", :secRC => "-20,60",:maskName => 'wet_c'))
+      ip.show(ip.scalarPlot(OCELSM_PLOT_TEST_FILE,ofile+'_05',"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:secLC => "-50,-60", :secRC => "0,60",:maskName => 'wet_c',:secMode => 'circle'))
+      ip.show(ip.scalarPlot(OCELSM_PLOT_TEST_FILE,ofile+'_06',"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:maskName => 'wet_c'))
+      ip.show(ip.scalarPlot(OCELSM_PLOT_TEST_FILE,ofile+'_07',"T",:levIndex => 2,:maskName => 'wet_c'))
       ip.isIcon = false
       ip.show(ip.scalarPlot("remapnn_r90x45_oce.nc","reg_"+ofile,"T",:levIndex => 2,:mapType => "ortho"))
       ip.show(ip.scalarPlot("remapnn_r90x45_oce.nc","reg_"+ofile,"T",:vecVars => "u-veloc,v-veloc",:levIndex => 2,:mapType => "ortho"))

metadata

@@ -1,100 +1,93 @@
 --- !ruby/object:Gem::Specification
 name: iconPlot
 version: !ruby/object:Gem::Version
-  version: 0.0.6
-  prerelease: 
+  version: 0.0.8
 platform: ruby
 authors:
 - Ralf Mueller
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2012-09-21 00:00:00.000000000 Z
+date: 2013-03-26 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: cdo
   requirement: !ruby/object:Gem::Requirement
-    none: false
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
-    none: false
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: extcsv
   requirement: !ruby/object:Gem::Requirement
-    none: false
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
-    none: false
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: gnuplot
   requirement: !ruby/object:Gem::Requirement
-    none: false
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
-    none: false
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
-description: ! 'Plot with ncl via Ruby: requires NCL 6.* and CDO 1.5.*'
+description: 'Plot with ncl via Ruby: requires NCL 6.* and CDO 1.5.*'
 email: stark.dreamdetective@gmail.com
-executables: []
+executables:
+- nclsh
 extensions: []
 extra_rdoc_files: []
 files:
 - lib/iconPlot.rb
 - gemspec
-- contrib/nclsh
+- bin/nclsh
 - lib/icon_plot.ncl
 - lib/icon_plot_lib.ncl
 - test/test_iconPlot.rb
 homepage: https://github.com/Try2Code/iconPlot
 licenses:
 - GPLv2
+metadata: {}
 post_install_message: 
 rdoc_options: []
 require_paths:
 - lib
 required_ruby_version: !ruby/object:Gem::Requirement
-  none: false
   requirements:
-  - - ! '>='
+  - - '>='
     - !ruby/object:Gem::Version
       version: '1.8'
 required_rubygems_version: !ruby/object:Gem::Requirement
-  none: false
   requirements:
-  - - ! '>='
+  - - '>='
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 1.8.23
+rubygems_version: 2.0.0
 signing_key: 
-specification_version: 3
+specification_version: 4
 summary: Plot ICON output with ncl via Ruby
 test_files:
 - test/test_iconPlot.rb