[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