[ncl-talk] Trend_plot_with_90%_CI
Kunal Bali
kunal.bali9 at gmail.com
Mon Aug 19 05:02:12 MDT 2019
I think I have sorted out this issue by changing
copy_VarAtts(var_a, sig_Dec)
to
copy_VarMeta(var_a, sig_Dec)
---
[image: plot03.png]
Kunal Bali
On Mon, Aug 19, 2019 at 4:09 PM Kunal Bali <kunal.bali9 at gmail.com> wrote:
> 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/08fa8a46/attachment-0001.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/08fa8a46/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot03.png
Type: image/png
Size: 175257 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190819/08fa8a46/attachment-0003.png>
More information about the ncl-talk
mailing list