[ncl-talk] No valid values in scalar field

[Ciao Kai Liang]

Hi,

I have one variable "dox" with (360x720x24). It is ozone global map and
contains 24 different perturbation results. I would like to plot them on
multiple panels in a single sheet. However, I don't know why I can't plot
the first 3 (see pic). The output of "dox"from printVarSummary and plotting
script shows below.

Thx
CiaoKai

dox-----------------------------------------------------------------------------
Variable: dox
Type: float
Total Size: 24883200 bytes
            6220800 values
Number of Dimensions: 3
Dimensions and sizes:   [lat | 360] x [lon | 720] x [24]
Coordinates:
            lat: [-89.75..89.75]
            lon: [-179.75..179.75]
Number Of Attributes: 1
  _FillValue :  9.96921e+36

plot
script----------------------------------------------------------------------------------------------

plots = new(24,graphic)
  wks_type = "png"
  wks_name = "dox_panel"
  wks = gsn_open_wks(wks_type,wks_name)         ; send graphics to PDF file

  res                              = True
  res at gsnDraw              = False              ; don't draw
  res at gsnFrame             = False              ; don't advance frame

  res at cnLevelSelectionMode = "ManualLevels"     ; set manual contour levels
  res at cnMinLevelValF       = -2.                ; set min contour level
  res at cnMaxLevelValF       =  2.                ; set max contour level
  res at cnLevelSpacingF      =  1.                ; set contour spacing

  res at cnFillOn          = True
  res at cnLinesOn         = False
  res at cnLineLabelsOn    = False
  res at lbLabelBarOn      = False
  res at cnFillPalette     = "BlWhRe"              ; Blue-Red colormap
  res at gsnRightString    = "ppb"                 ; "~" is txFuncCode

  leftstring = (/"a","g","m","b","h","n",\
                 "c","i","o","d","j","p",\
                 "e","k","q","f","l","r",\
                 "s","t","u","v","w","X"/)

  do i=0,23
    res at gsnLeftString  =  "("+leftstring(i)+")"
    plots(i) = gsn_csm_contour_map_ce(wks,dox(:,:,i),res)
  end do

  pnlres = True
  pnlres at txString          = "Dox"
  pnlres at gsnPanelLabelBar  = True
  pnlres at lbLabelAutoStride = True

  gsn_panel(wks,plots,(/6,4/),pnlres)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160618/af709ff1/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: htap2_dox_panel.png
Type: image/png
Size: 370189 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160618/af709ff1/attachment-0001.png