[ncl-talk] Trend_plot_with_90%_CI

Kunal Bali kunal.bali9 at gmail.com
Mon Aug 19 04:39:33 MDT 2019


unfortunately, it's not.


---
Kunal Bali






On Mon, Aug 19, 2019 at 2:11 PM S Br <sbr.climate at gmail.com> wrote:

> Please try adding
> res2 at gsnAddCyclic     = False
>
> it may work fine.
>
> Cheers
> S
>
> On Sun, Aug 18, 2019 at 7:05 PM Kunal Bali via ncl-talk <ncl-talk at ucar.edu>
> wrote:
>
>> Dear NCL users,
>>
>> I also tried the another function (given below) to calculate 90% CI, but
>> the CI values are not overlaying correctly.
>> please let me know the correct way to plot it.
>> (Here, I calculated trend values seperately using CDO), and then using
>> resultant file to calculate 90% significant values.
>>
>> [image: plot03.png]
>>
>>      var_a    = a->BCSMASS(0,:,:).
>>
>>
>> ;**********************************************
>>>
>>> ;Calculating significance
>>>
>>> ;**********************************************
>>>
>>>
>>>>  prob_Dec=rtest(var_a,5,0)
>>>
>>>  sig_Dec =100*(1-prob_Dec)
>>>
>>>  sig_Dec = where(ismissing(var_a),sig_Dec at _FillValue,sig_Dec)
>>>
>>>
>>>>  printVarSummary(sig_Dec)
>>>
>>>  printMinMax(sig_Dec, True)
>>>
>>>
>>>>  ;copy_VarAtts(var_a, sig_Dec)
>>>
>>>  ;copy_VarCoords_1(var_a,sig_Dec)
>>>
>>>
>>>
>>>   sig_Dec = mask(sig_Dec, sig_Dec.ge.99,False)
>>>
>>>   sig_Dec = mask(sig_Dec, sig_Dec.lt.90,False)
>>>
>>>
>>>
>>> printVarSummary(sig_Dec)
>>>
>>> printMinMax(sig_Dec, True)
>>>
>>>
>>
>>   res2=True
>>>
>>>   res2 at gsnLeftString     = ""
>>>
>>>   res2 at gsnRightString     = ""
>>>
>>>   res2 at gsnDraw              = False           ; Do not draw plot
>>>
>>>   res2 at gsnFrame             = False           ; Do not advance frome
>>>
>>>   ;res2 at cnMonoFillColor      = False
>>>
>>>   res2 at cnLinesOn            = False    ; turn off contour lines
>>>
>>>   res2 at lbLabelBarOn         = False                            ; turn
>>>> off label bar
>>>
>>>   res2 at cnFillColor          = "black"
>>>
>>>   res2 at cnLevelSelectionMode = "ManualLevels"   ; manually specify
>>>> contour levels
>>>
>>>   res2 at cnMinLevelValF       = 90        ; min level
>>>
>>>   res2 at cnMaxLevelValF       = 99           ; max level
>>>
>>>   res2 at cnLevelSpacingF      =  0.1             ; contour interval
>>>
>>>   res2 at gsnSpreadColors      = False
>>>
>>>   ;res2 at cnConstFEnableFill = True
>>>
>>>
>>>
>>>   res2 at trYReverse             =  True
>>>
>>>   res2 at tfDoNDCOverlay         =  False                ;-- transform to
>>>> standard lat/lon
>>>
>>>   res2 at trGridType             = "TriangularMesh"
>>>
>>>
>>>>
>>>>   plotb=gsn_csm_contour(wks,sig_Dec, res2) ; contours are at 90,100
>>>
>>>
>>>> opt = True
>>>
>>> opt at gsnShadeFillType = "pattern" ; pattern fill
>>>
>>> opt at gsnShadeHigh =17    ; use pattern #2
>>>
>>> opt at gsnShadeFillDotSizeF = 0.005
>>>
>>> opt at gsnShadeFillScaleF = 1.5
>>>
>>>
>>>> plotc = gsn_contour_shade(plotb,-999.,90.,opt)
>>>
>>>
>>>>
>>>
>>> ;************************************************
>>>
>>> ; Attach text to plot using plot coordinates.
>>>
>>> ;************************************************
>>>
>>>   txres               = True                     ; text mods desired
>>>
>>>   txres at txFontHeightF = 0.02                     ; font smaller.
>>>> default big
>>>
>>>   txres at txFontThicknessF = 1
>>>
>>>   txres at txBackgroundFillColor = "white"
>>>
>>>
>>>>
>>>>   ;---Resource list for shapefile outlines
>>>
>>>  plres                 = True
>>>
>>>  plres at gsLineColor       = "black"
>>>
>>>  plres at gsLineThicknessF  = "6"
>>>
>>>  dumstr           = unique_string("poly")
>>>
>>>
>>>>
>>>>
>>>
>>>  id = gsn_add_shapefile_polylines(wks,plotc,filename,plres)
>>>
>>>
>>>>
>>>
>>>
>>>>  overlay(plot,plotc)
>>>
>>>  draw(plot)
>>>
>>>  frame(wks)
>>>
>>> end
>>>
>>>
>>>
>>
>>
>>
>>
>>
>>
>> ---
>> Kunal Bali
>>
>>
>>
>>
>>
>>
>> On Sun, Aug 18, 2019 at 2:08 PM Kunal Bali <kunal.bali9 at gmail.com> wrote:
>>
>>> Dear NCL users,
>>>
>>> I want to highlight the 90% CI over trend plot.
>>>
>>> I used the given code for climatology, but not working properly during
>>> trend calculation.
>>>
>>>
>>>> arr95 = var_a(0,:,:)
>>>> arr95 = arr95 at _FillValue
>>>> arr5 = arr95
>>>> ntime = dimsizes(var_a&time)
>>>> i95 = toint(ntime*.95)-1  ; approximate 95% index
>>>> ip    = dim_pqsort_n(var_a, 2, 0)   ; ascending order, sort var_a in
>>>> ascending order at the 0th dimension.
>>>> var_a95 = var_a(i95,:,:)
>>>>
>>>
>>> But I am not getting how to incorporate trend function
>>>
>>>
>>>>   var_b    = b->BCSMASS
>>>>
>>>    ts         = var_b(lat|:,lon|:,time|:)
>>>
>>>    time     = ts&time                       ; year since 2001-01-01
>>>
>>>    rc         = regCoef(time,ts)
>>>
>>>
>>>
>>>
>>>
>>> please provide some information on that. For more clarification I have
>>> shared my file with dropbox.
>>>
>>> Thank you
>>>
>>>
>>> regards
>>> Kunal
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> https://www.dropbox.com/s/15z3b06liu8z8a2/YM1.nc?dl=0
>>>
>> _______________________________________________
>> 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/20190819/474d7b5b/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot03.png
Type: image/png
Size: 194610 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190819/474d7b5b/attachment.png>


More information about the ncl-talk mailing list