画不出图来，提示fatal:ContourPlotDraw: Workspace reallocation would exceed ma...

lisixiuhe

;*******************************************************
; lcnative_2.ncl
;*******************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
begin
;************************************************
; open file and read in data
;************************************************
  f     = addfile ("CoordinatesInfo.nc", "r")
  nf     = addfile ("nc.nc", "r")
;print(f)
  t = nf->ps(0,0,:,:)
;print(t)   
  lat2d     = f->lat
  lon2d     = f->lon
  dimll  = dimsizes(lat2d)                ; get size of dimensions
;print(dimll)
  nlat   = dimll(0)
  mlon   = dimll(1)
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("png","aaa_out")         ; open a workstation
gsn_define_colormap(wks,"gui_default")    ; choose colormap
  ;gsn_define_colormap(wks,"test111") ; select color map
  res                   = True              ; plot mods desired
  res@tiMainString      = "test" ;Native Lambert Conformal"
;************************************************
; lambert conformal projections are limited using
; the corners method rather than the latlon method
; seen for cylindrical equidistant projections
;************************************************
  res@mpLimitMode       = "Corners"         ; choose range of map
  res@mpLeftCornerLatF       = lat2d(0,0)
  res@mpLeftCornerLonF       = lon2d(0,0)
  res@mpRightCornerLatF      = lat2d(nlat-1,mlon-1)
  res@mpRightCornerLonF      = lon2d(nlat-1,mlon-1)
;************************************************
; native lambert grids read in from grib files are
; different than those read in from netCDF files.
; NCL automatically looks for the parallel and
; meridian information on the grib file and attaches
; this information as attributes to the lat2d array.
;************************************************
  res@mpProjection          = "LambertConformal"
  res@mpLambertParallel1F   = nf->Latin1(0)
  res@mpLambertParallel2F   = nf->Latin2(0)
  res@mpLambertMeridianF    = nf->LoV(0)       ; ncl adds from grib file

  res@pmTickMarkDisplayMode = "Always"
  res@mpFillOn              = True           ; turn on map fill
  res@mpOutlineDrawOrder    = "PostDraw"      ; draw continental outline last
  res@mpOutlineBoundarySets = "AllBoundaries" ; all boundaries
;****************************************************************************
; usually, when data is placed onto a map, it is TRANSFORMED to the specified
; projection. Since this model is already on a native lambert conformal grid,
; we want to turn OFF the tranformation.
;****************************************************************************
  res@tfDoNDCOverlay = True
;****************************************************************************
  res@cnFillOn              = True            ; color plot desired
  res@cnLinesOn             = False           ; turn on contour lines
  res@cnLineLabelsOn    = False     ;turn on label lines

  res@gsnSpreadColors       = True            ; use full range of colormap
  res@gsnAddCyclic          = False           ; regional data
  ;res@tmXTLabelsOn          = False
  ;res@tmYRLabelsOn          = False
  res@mpGridAndLimbOn = True                ;turn on grid line
  
plot = gsn_csm_contour_map(wks,t,res)
end

lisixiuhe

自己动手，丰衣足食啊http://www.ncl.ucar.edu/Document/Functions/Built-in/NhlGetWorkspaceObjectId.shtml

lisixiuhe

Variable: t
Type: float
Total Size: 13694804 bytes
            3423701 values
Number of Dimensions: 2
Dimensions and sizes:        [y | 1801] x [x | 1901]
Coordinates:
Number Of Attributes: 10
  z :         0
  navigation_dim :        nav
  record :        valtime, reftime
  _FillValue :        1e+37
  long_name :        surface pressure
  units :        pascal
  valid_range :        (  0, 120000 )
  LAPS_var :        PS
  lvl_coord :        AGL
  LAPS_units :        PA
这是ps变量的概要信息