[ncl-talk] ContourPlotDraw and PlotManagerDraw error while ploting SST

Amit Kumar amitkumar.777 at hotmail.com
Tue Aug 30 00:57:21 MDT 2016


Hi All,

While ploting south hemisphere co-ordinates i'm getting following error.

warning:ContourPlotDraw: out of range coordinates encountered; standard RasterFill rendering method may be unreliable;
 consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values

I'm using

  res at mpMinLonF     = 40.0
  res at mpMaxLonF     = 110.0
  res at mpMinLatF     = -50.0
  res at mpMaxLatF     = 25.0

to zoom into southern hemisphere, but it is not plotting southern latitudes is giving above error.

Is there any other way to write southern latitudes?


Regards

Amit

________________________________
From: Mary Haley <haley at ucar.edu>
Sent: Tuesday, August 30, 2016 1:58:13 AM
To: Amit Kumar
Cc: ncl-talk at ucar.edu; Guido Cioni
Subject: Re: [ncl-talk] ContourPlotDraw and PlotManagerDraw error while ploting SST

Amit,

Guido's answer is correct, and I wanted to point out that you can also set this resource inside your NCL script. See:

http://www.ncl.ucar.edu/Document/Language/error_messages.shtml#WorkspaceReallocation

The above also recommends trying raster contouring, if you have a large array that you are trying to plot:


  res at cnFillMode<http://www.ncl.ucar.edu/Document/Graphics/Resources/cn.shtml#cnFillMode> = "RasterFill"
  res at cnRasterSmoothingOn<http://www.ncl.ucar.edu/Document/Graphics/Resources/cn.shtml#cnRasterSmoothingOn> = True   ; optional, but may produce better results

Try it first without the "cnRasterSmoothingOn" set to True. If it looks too blocky, then add that resource


--Mary


On Mon, Aug 29, 2016 at 3:55 AM, Guido Cioni <guidocioni at gmail.com<mailto:guidocioni at gmail.com>> wrote:
As stated here http://www.ncl.ucar.edu/Document/Graphics/hlures.shtml [at the bottom] try to increase the default contour memory in the .hluresfile.
Regards,

Guido Cioni
http://guidocioni.altervista.org

On 29 Aug 2016, at 10:17, Amit Kumar <amitkumar.777 at hotmail.com<mailto:amitkumar.777 at hotmail.com>> wrote:


Hi all,

I'm trying to plot daily file of SST using NCL. But I'm getting following error:

fatal:ContourPlotDraw: Workspace reallocation would exceed maximum size 100000000
fatal:ContourPlotDraw: draw error
fatal:ContourPlotDraw: draw error
fatal:PlotManagerDraw: error in plot draw
fatal:_NhlPlotManagerDraw: Draw error

I'm able to plot HEM, OLR which are in same format and with same NCL script but with SST I'm encountering these errors.


My NC file description is as shown below:

dimensions:
      x = 1744
      y = 657
      time = 1  // unlimited
   variables:
      double longitude ( y, x )
         standard_name :    longitude
         long_name :    longitude
         units :    degrees_east
         _CoordinateAxisType :    Lon

      double latitude ( y, x )
         standard_name :    latitude
         long_name :    latitude
         units :    degrees_north
         _CoordinateAxisType :    Lat

      double time ( time )
         standard_name :    time
         units :    hours since 1-1-1 00:00:00
         calendar :    standard
         axis :    T

      double SST_DLY ( time, y, x )
         coordinates :    latitude longitude
         _FillValue :    -9e+33
         missing_value :    -9e+33
         fillvalue :    -9.e+33


My script is as below:

load "/usr/share/ncarg/nclscripts/csm/gsn_code.ncl"
load "/usr/share/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/share/ncarg/nclscripts/csm/contributed.ncl"

begin

 a = addfile("/media/amit/E/Research/SWM/SST/AUG/output/india/07AUG.nc","r")
 ;filename = "/media/amit/E/India_shapefile/Districtbound(Satmet).shp"
 filename = "/media/amit/E/India_shapefile/State_bound(Satmet).shp"

  var     = doubletofloat(a->SST_DLY(0,:,:))
  lat     =  a->latitude
  lon     =  a->longitude

wks = gsn_open_wks("png","07AUG")                ; send graphics to PNG file

  res                      = True                  ; plot mods desired
  res at gsnDraw              = False
  res at gsnFrame             = False
  res at gsnAddCyclic         = False
  res at gsnStringFontHeightF = 0.015
  res at cnFillOn             = True                  ; turn on color fill
  res at gsnMaximize          = True     ; maximize plot in frame
  res at mpOutlineOn   = False   ; Use outlines from shapefile
  res at mpFillOn      = False
  res at cnLinesOn         = False    ; turn off contour lines
  res at cnLineLabelsOn    = False    ; turn off line labels
  res at cnLevelSelectionMode = "ExplicitLevels"
  res at cnLevels             = (/273,283,293,294,295,296,297,298,299,300,301,302,303,304,305/) ;   15 contour values
  res at cnFillPalette            = "rainbow"
  res at lbLabelBarOn     = False
  res at lbOrientation   = "Vertical"
  res at sfXArray             = lon
  res at sfYArray        =  lat
 ; res at mpMinLonF     = 1.0
 ; res at mpMaxLonF     = 163.0
 ; res at mpMinLatF     = -81.0
 ; res at mpMaxLatF     = 25.0
  res at pmTickMarkDisplayMode     = "Always"
  ;res at tiXAxisString            = "Longitude"
  ;res at tiYAxisString            = "Lattitude"
  res at tiXAxisFontHeightF     = 0.02
  res at tiYAxisFontHeightF     = 0.02
plot = gsn_csm_contour_map(wks,var, res)        ; create plot

;---Resource list for shapefile outlines
 plres                     = True
 plres at gsLineColor       = "black"
 plres at gsLineThicknessF  = "1.0"
 dumstr                = unique_string("poly")
 id = gsn_add_shapefile_polylines(wks,plot,filename,plres)
 draw(plot)
 frame(wks)
end



_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk


_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160830/8c808f07/attachment.html 


More information about the ncl-talk mailing list