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

[S Br]

Hi Rashed,
I could solve the issue by using the mask function.

Thanks.


On Wed, Jul 31, 2019 at 9:18 PM Rashed Mahmood <rashidcomsis at gmail.com>
wrote:

> 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/20190802/3573178f/attachment.html>