
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" style="font-size:10pt;color:#2672EC;background-color:#FFFFFF;font-family:Arial,Helvetica,sans-serif;">

<p></p>

<div>Hi all,<br>

<br>

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

<br>

<span style="color: rgb(255, 0, 0);">fatal:ContourPlotDraw: Workspace reallocation would exceed maximum size 100000000</span><br>

<span style="color: rgb(255, 0, 0);">fatal:ContourPlotDraw: draw error</span><br>

<span style="color: rgb(255, 0, 0);">fatal:ContourPlotDraw: draw error</span><br>

<span style="color: rgb(255, 0, 0);">fatal:PlotManagerDraw: error in plot draw</span><br>

<span style="color: rgb(255, 0, 0);">fatal:_NhlPlotManagerDraw: Draw error</span><br>

</div>

<p></p>

<p><br>

</p>

<p>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.</p>

<p><br>

</p>

<p></p>

<div>My NC file description is as shown below:<br>

<br>

dimensions:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; x = 1744<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; y = 657<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time = 1&nbsp; // unlimited<br>

&nbsp;&nbsp; variables:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double longitude ( y, x )<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; standard_name :&nbsp;&nbsp; &nbsp;longitude<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; long_name :&nbsp;&nbsp; &nbsp;longitude<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; units :&nbsp;&nbsp; &nbsp;degrees_east<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; _CoordinateAxisType :&nbsp;&nbsp; &nbsp;Lon<br>

<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double latitude ( y, x )<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; standard_name :&nbsp;&nbsp; &nbsp;latitude<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; long_name :&nbsp;&nbsp; &nbsp;latitude<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; units :&nbsp;&nbsp; &nbsp;degrees_north<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; _CoordinateAxisType :&nbsp;&nbsp; &nbsp;Lat<br>

<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double time ( time )<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; standard_name :&nbsp;&nbsp; &nbsp;time<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; units :&nbsp;&nbsp; &nbsp;hours since 1-1-1 00:00:00<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calendar :&nbsp;&nbsp; &nbsp;standard<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; axis :&nbsp;&nbsp; &nbsp;T<br>

<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double SST_DLY ( time, y, x )<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; coordinates :&nbsp;&nbsp; &nbsp;latitude longitude<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; _FillValue :&nbsp;&nbsp; &nbsp;-9e&#43;33<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; missing_value :&nbsp;&nbsp; &nbsp;-9e&#43;33<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fillvalue :&nbsp;&nbsp; &nbsp;-9.e&#43;33<br>

<br>

<br>

My script is as below:<br>

<br>

<div>load &quot;/usr/share/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>

load &quot;/usr/share/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br>

load &quot;/usr/share/ncarg/nclscripts/csm/contributed.ncl&quot;<br>

<br>

begin<br>

<br>

&nbsp;a = addfile(&quot;/media/amit/E/Research/SWM/SST/AUG/output/india/07AUG.nc&quot;,&quot;r&quot;)<br>

&nbsp;;filename = &quot;/media/amit/E/India_shapefile/Districtbound(Satmet).shp&quot;<br>

&nbsp;filename = &quot;/media/amit/E/India_shapefile/State_bound(Satmet).shp&quot;<br>

<br>

&nbsp; var &nbsp;&nbsp; &nbsp;= doubletofloat(a-&gt;SST_DLY(0,:,:))<br>

&nbsp; lat &nbsp;&nbsp; &nbsp;=&nbsp; a-&gt;latitude<br>

&nbsp; lon &nbsp;&nbsp; &nbsp;=&nbsp; a-&gt;longitude<br>

<br>

wks = gsn_open_wks(&quot;png&quot;,&quot;07AUG&quot;)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; send graphics to PNG file<br>

<br>

&nbsp; res&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; plot mods desired<br>

&nbsp; res@gsnDraw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

&nbsp; res@gsnFrame&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

&nbsp; res@gsnAddCyclic&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

&nbsp; res@gsnStringFontHeightF = 0.015<br>

&nbsp; res@cnFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; turn on color fill<br>

&nbsp; res@gsnMaximize&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp; ; maximize plot in frame<br>

&nbsp; res@mpOutlineOn&nbsp;&nbsp; = False&nbsp;&nbsp; ; Use outlines from shapefile<br>

&nbsp; res@mpFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

&nbsp; res@cnLinesOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp; ; turn off contour lines<br>

&nbsp; res@cnLineLabelsOn&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp; ; turn off line labels<br>

&nbsp; res@cnLevelSelectionMode = &quot;ExplicitLevels&quot;<br>

&nbsp; res@cnLevels&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = (/273,283,293,294,295,296,297,298,299,300,301,302,303,304,305/) ;&nbsp;&nbsp; 15 contour values<br>

&nbsp; res@cnFillPalette&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;rainbow&quot;<br>

&nbsp; res@lbLabelBarOn &nbsp;&nbsp; &nbsp;= False<br>

&nbsp; res@lbOrientation&nbsp;&nbsp; = &quot;Vertical&quot;<br>

&nbsp; res@sfXArray &nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = lon<br>

&nbsp; res@sfYArray&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;=&nbsp; lat<br>

&nbsp;; res@mpMinLonF &nbsp;&nbsp; &nbsp;= 1.0<br>

&nbsp;; res@mpMaxLonF &nbsp;&nbsp; &nbsp;= 163.0<br>

&nbsp;; res@mpMinLatF &nbsp;&nbsp; &nbsp;= -81.0<br>

&nbsp;; res@mpMaxLatF &nbsp;&nbsp; &nbsp;= 25.0<br>

&nbsp; res@pmTickMarkDisplayMode &nbsp;&nbsp; &nbsp;= &quot;Always&quot;<br>

&nbsp; ;res@tiXAxisString&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= &quot;Longitude&quot;<br>

&nbsp; ;res@tiYAxisString&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= &quot;Lattitude&quot;<br>

&nbsp; res@tiXAxisFontHeightF &nbsp;&nbsp; &nbsp;= 0.02<br>

&nbsp; res@tiYAxisFontHeightF &nbsp;&nbsp; &nbsp;= 0.02 <br>

plot = gsn_csm_contour_map(wks,var, res)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; create plot<br>

<br>

;---Resource list for shapefile outlines<br>

&nbsp;plres&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>

&nbsp;plres@gsLineColor&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;black&quot;<br>

&nbsp;plres@gsLineThicknessF&nbsp; = &quot;1.0&quot;<br>

&nbsp;dumstr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; = unique_string(&quot;poly&quot;)<br>

&nbsp;id = gsn_add_shapefile_polylines(wks,plot,filename,plres)<br>

&nbsp;draw(plot)<br>

&nbsp;frame(wks)<br>

end</div>

<br>

</div>

<br>

<p></p>

</div>

</body>

</html>