Re: DEM plot: memory allocation error

Hi Noel,

Memory issues aside, it doesn't seem like you'll be able to graphically depict both the full areal extent of that data (1500miles wide?) and resolve its full detail. What if you were to resample it at a lower resolution; perhaps something as simple as:

 plot = gsn_csm_contour_map(wks, dem90m(::2,::2), res1)

Just a thought…

Rick

On Apr 11, 2013, at 9:01 AM, Noel Aloysius <noel.aloysius@gmail.com> wrote:

> Hi Mary,
>
> The variable type is float. Should I convert the type to integer and work.
>
> Noel
>
> Noel Aloysius
>
>
> On Thu, Apr 11, 2013 at 10:53 AM, Mary Haley <haley@ucar.edu> wrote:
> Hi Noel,
>
> Sorry, I've been away from ncl-talk for a bit.
>
> This is not something you can fix with wsMaximumSize, unfortunately. I think internally, NCL might be making a copy of your array, and simply running out of memory.
>
> What type is the variable?
>
> I'm hoping Dave Brown will have an idea for a work-around, that might prevent a copy from taking place.
>
> --Mary
>
>
> On Apr 8, 2013, at 1:17 PM, Noel Aloysius wrote:
>
> > Dear NCL-talk,
> >
> > I have a 90m DEM with the following dimensions (/27260,26610/), ~ 2.9GB in size. The script below creates a spatial plot of the DEM.
> >
> > The plot is successfully created if I restrict the region as follows,
> >
> > plot = gsn_csm_contour_map(wks, dem90m(5000:27259,5000:26609), res1)
> >
> > I get an error otherwise.
> >
> > fatal:NhlRealloc Failed:[errno=12]
> > fatal:ContourPlotDraw: dynamic memory allocation error
> > fatal:ContourPlotDraw: draw error
> > warning:WorkstationDeactivate: workstation not active or not opened
> >
> > Is there a way around?
> >
> > Thanks,
> > Noel
> >
> >
> > ;*****************************************************************
> > ffig = diro + "dem90m"
> > wks = gsn_open_wks(type, ffig)
> > plot = new(1, graphic)
> > gsn_define_colormap(wks, "ncl_default")
> > maps = NewList("fifo")
> >
> > ;--------------------------------------------------------
> > setvalues NhlGetWorkspaceObjectId()
> > "wsMaximumSize": 2000000000
> > end setvalues
> > ;--------------------------------------------------------
> >
> > res1 = True ; plot mods desired
> > res1@gsnDraw = False
> > res1@gsnFrame = False
> >
> > res1@tiMainString = "90m DEM"
> > res1@gsnPaperOrientation = "auto"
> > res1@gsnMaximize = True
> > res1@gsnSpreadColors = True
> > res1@gsnAddCyclic = False
> >
> > res1@lbLabelBarOn = True
> >
> > res1@cnFillMode = "RasterFill"
> > res1@cnFillOn = True
> > res1@cnLinesOn = False
> >
> > res1@mpCenterLonF = (regBnds(2)+regBnds(3))/2 ; center the map arounf this lon.
> > res1@mpMinLatF = regBnds(0) ; select a subregion
> > res1@mpMaxLatF = regBnds(1)
> > res1@mpMinLonF = regBnds(2)
> > res1@mpMaxLonF = regBnds(3)
> >
> > plot = gsn_csm_contour_map(wks, dem90m(5000:27259,5000:26609), res1)
> > ;plot = gsn_csm_contour_map(wks, dem90m, res1)
> >
> > ; draw polygon
> > lnres1 = True
> > lnres1@gsLineColor = "black"
> > lnres1@gsLineThicknessF = 2.0
> > ListPush(maps, gsn_add_shapefile_polylines(wks,plot,fshp1,lnres1))
> >
> > draw(wks)
> > frame(wks)
> > ;***************************************************************
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 11 09:16:33 2013