[ncl-talk] (no subject)

[Dr. Muhammad Afzaal]

Dear NCL community!

I am trying to plot significance stipples above trend plot. I took help
from examples, and threads from NCL-talk. Trend is displayed correctly, but
stipples(>=0.95) do not appear. I checked the values, those are inbetween 0
and 1. Following is my code for plotting. Am I missing some resource?
Please give suggestions to correct it.
Thank you in anticipation

Please not p(0,:,:) is significance and p(1,:,:) is trend
;===============================================================

  wks          = gsn_open_wks ("png","spi_shp_trnd") ; send graphics to PNG
file
  cmap = read_colormap_file("BlueWhiteOrangeRed")  ; read color data

  res1                  = True
  res1 at gsnDraw          = False              ; don't draw
  res1 at gsnFrame         = False              ; don't advance frame
  res1 at gsnAddCyclic     = False            ; This is very important for
subset data
  res1 at gsnMaximize      = True

  lat_min = 12
  lat_max = 33
  lon_min = 30
  lon_max = 60

  res1 at mpMinLatF           = min(prc&lat)
  res1 at mpMaxLatF          = max(prc&lat)
  res1 at mpMinLonF          = min(prc&lon)
  res1 at mpMaxLonF         = max(prc&lon)
  res1 at mpCenterLonF     = (res1 at mpMinLonF + res1 at mpMaxLonF) / 2.

;=============plot trend===============================
  res1 at cnFillPalette    = cmap(::-1,:)       ; set and reverse color map
  res1 at cnLinesOn        =  False                 ; Turn off contour lines
  res1 at cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
  res1 at cnMinLevelValF   = -0.0025            ; set min contour level
  res1 at cnMaxLevelValF   =  0.0025            ; set max contour level
  res1 at cnLevelSpacingF  =  0.0005            ; set contour spacing
  res1 at cnLabelMasking   = True               ; Turn off countour label

  plot1=gsn_csm_contour_map(wks,p(1,:,:),res1)

;==========plot significance==========================
  res2 = True
  res2 at gsnDraw       = False
  res2 at gsnFrame       = False
  res2 at cnFillOn       = True
  res2 at cnMonoFillColor =True
  res2 at cnLinesOn      = False
  res2 at cnLineLabelsOn = False
  res2 at cnLevelSelectionMode = "Explicitlevels"
  res2 at cnLevels     =(/ 0.95 /)
  res2 at cnFillPattern     =(/-1,17/)
  res2 at cnInfoLabelOn    = False
  res2 at cnLabelMasking   = False
  res2 at gsnAddCyclic     = False
  res2 at cnFillScaleF = 0.7
  res2 at cnFillDotSizeF       = 0.005

  plot2=gsn_csm_contour(wks,p(0,:,:),res2)

  overlay(plot1,plot2)

  draw(plot1)
  frame(wks)

end

;===================

Regards
Dr Afzaal
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180118/2af0ceaf/attachment.html>