[ncl-talk] Plot significant difference only (shaded color map)

[Rashed Mahmood]

Change the following to True to get color filled areas:
res at cnFillOn            = False    ; This is in your resources.

If this does not solve the issue then please attach a figure  and full
script, and if possible some data.
Good luck!



On Wed, Jul 31, 2019 at 11:45 AM S Br <sbr.climate at gmail.com> wrote:

> Hi Rashed,
> I already added but I have sent only the main part of the script.
>
> On Wed, Jul 31, 2019 at 6:50 PM Rashed Mahmood <rashidcomsis at gmail.com>
> wrote:
>
>> Try adding following at the end, after overlay:
>>
>> draw(plot)
>> frame(wks)
>>
>> On Wed, Jul 31, 2019 at 7:57 AM S Br via ncl-talk <ncl-talk at ucar.edu>
>> wrote:
>>
>>> Hello NCL,
>>> I am plotting the difference between two files as per the script written
>>> below, but I want to display (contour shade) only the value that are
>>> significance at 95% level. I am able to stripe the significance difference
>>> but can't display the difference that are significant (I want to retain the
>>> original color shade of the difference)
>>>
>>> diffposglob1 = avgposglobst1-avgstctrl
>>> probtposglob1 = 100.*(1.
>>> -ttest(avgposglobst1,varposglobst1,50,avgstctrl,varctrl,50,True,False)) ;
>>> all areas whose value > 95 are significant
>>>
>>> wks = gsn_open_wks("x11","surtemp_sig_diff")
>>>  res                     = True
>>>   res at gsnDraw             = False            ; don't draw
>>>   res at gsnFrame            = False            ; don't advance frame
>>>   res at cnFillOn            = False            ; turn on color
>>>   res at cnLinesOn           = True            ; turns off contour lines.
>>>   res at cnLineLabelsOn      = False            ; turn off line labels
>>>   res at cnFillPalette       = "temp_diff_18lev"   ; set color map
>>>   res at cnInfoLabelOn       = False           ; turn off cn info label
>>> explanation
>>>   res at lbLabelBarOn        = True           ; turn off individual cb's
>>>   res at cnLevelSelectionMode =  "ManualLevels"
>>>    res at cnMinLevelValF       = -1.5
>>>    res at cnMaxLevelValF       =  1.5
>>>    res at cnLevelSpacingF      =  0.3
>>>    res at tmXBLabelFontHeightF = 0.030                 ; resize tick labels
>>>    res at tmYLLabelFontHeightF = 0.030
>>>    res at gsnStringFontHeightF = 0.030
>>> ;***********************************************************
>>>   res2 = True                            ; res2 probability plots
>>>   res2 at gsnDraw             = False       ; Do not draw plot
>>>   res2 at gsnFrame            = False       ; Do not advance frome
>>>   res2 at cnInfoLabelOn       = False
>>>   res2 at cnLinesOn           = False       ; do not draw contour lines
>>>   res2 at cnLineLabelsOn      = False       ; do not draw contour labels
>>>   res2 at cnFillScaleF        = 0.3         ; add extra density
>>> ;  sgres at cnFillDotSizeF       = 0.003
>>>   res2 at cnLevelSelectionMode = "ExplicitLevels"       ; set explicit
>>> cnlev
>>>   res2 at cnLevels   = (/95./)              ; only have 1 contour level
>>>   res2 at gsnRightString            = ""                ; Turn off
>>> subtitles
>>>   res2 at gsnLeftString             = ""
>>>   res at gsnAddCyclic                = True            ; add cyclic point
>>> ;***********************************************************
>>>   opt = True                                         ; set up parameters
>>> for pattern fill
>>>   opt at gsnShadeFillType = "pattern"                   ; specify pattern
>>> fill
>>>   opt at gsnShadeHigh      = 15                        ; stipple pattern
>>> ;  opt at gsnShadeFillScaleF = 5
>>> ;  opt at gsnShadeDotSizeF = 1                           ; make dots larger
>>> ;***********************************************************
>>> plot = gsn_csm_contour_map(wks,diffposglob1(0,:,:),res)
>>>   plot2   = gsn_csm_contour(wks,probtposglob1(0,:,:),res2)
>>>   plot2   = gsn_contour_shade(plot2,-999,95., opt)      ; stipple all
>>> areas >= 95  contour
>>>   overlay (plot, plot2)
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190731/a959768a/attachment.html>