DEM plot: memory allocation error

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
Received on Mon Apr 8 13:18:11 2013