[ncl-talk] draw_significance_values

[Kunal Bali]

Dear NCL users,

please ignore my previous mail. I have sorted out.



[image: Mailtrack]
<https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality1&>
Sent
with Mailtrack
<https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality1&>

---
Kunal Bali
Research Scholar
Centre for Atmospheric Sciences (Room-322, Block 6)
Indian Institute of Technology Delhi
Hauz Khas, New Delhi - 110016
India




On Wed, May 30, 2018 at 2:31 PM, Kunal Bali <kunal.bali9 at gmail.com> wrote:

> Dear NCL users,
>
> I am trying to draw the significance values at 95% on the correlation map
> but not able to succeed yet.
>
> First I calculate and plot the correlation using
>
>  ts1 = var_a(lat|:,lon|:,time|:)
>   ts2 = var_b(lat|:,lon|:,time|:)
>   printVarSummary(ts1)
>   printVarSummary(ts2)
>
> ;************************************************
> ; calculate cross correlations
> ;************************************************
> ; maxlag = 0                               ; set lag
> ; note, the max lag should not be more than N/4
>
> ; ccr = esccr(ts1,ts2(lat|:,lon|:,time|:), maxlag) ; calc cross
> correlations
> *  ccr = escorc(ts1,ts2) *
>   printVarSummary(ccr)
>   printMinMax(ccr,0)
>
> ; copy meta data and coordinate variables using contributed functions
>   copy_VarAtts(ts2, ccr)
>   copy_VarCoords_1(ts1,ccr)
>
>
>   *plot* = gsn_csm_contour_map(wks,ccr,res)
>
>
> *Then, I tried to calculate significance and draw over plot using*
>
>
> ;**********************************************
> ;Calculating significance
> ;**********************************************
>
>
> * prob_Dec=rtest(ccr,11,0) sig_Dec =100*(1-prob_Dec)*
>
> printVarSummary(sig_Dec)
>
> rescn2 = res
> rescn2 at cnMinLevelValF = 95
> rescn2 at cnMaxLevelValF = 100
> rescn2 at cnLevelSpacingF = 5.
> rescn2 at cnFillMode = "RasterFill"
> plotB=gsn_csm_contour(wks,sig_Dec,rescn2) ; contours are at 95,100
>
> opt = True
> opt at gsnShadeFillType = "pattern" ; pattern fill
> opt at gsnShadeHigh = 2 ; use pattern #2
> *plotB*=gsn_contour_shade(plotB,-999.,95.,opt)
>
>  overlay(*plot,plotB*)
>  draw(plot)
>  frame(wks)
>
> end
>
> After running the script, the summary and warning displayed as
>
> Variable: ts1
> Type: float
> Total Size: 141504 bytes
>             35376 values
> Number of Dimensions: 3
> Dimensions and sizes:    [lat | 67] x [lon | 48] x [time | 11]
> Coordinates:
>             lat: [   5..  38]
>             lon: [68.125..97.5]
>             time: [3134040..8392920]
> Number Of Attributes: 8
>   vmin :    -1e+15
>   vmax :    1e+15
>   fmissing_value :    1e+15
>   missing_value :    1e+15
>   _FillValue :    1e+15
>   units :    kg m-2 s-1
>
>
> Variable: ts2
> Type: float
> Total Size: 141504 bytes
>             35376 values
> Number of Dimensions: 3
> Dimensions and sizes:    [lat | 67] x [lon | 48] x [time | 11]
> Coordinates:
>             lat: [ 5..38]
>             lon: [68.125..97.5]
>             time: [928806..1016454]
> Number Of Attributes: 5
>   missing_value :    -32767
>   _FillValue :    -32767
>   units :    K
>
>
> Variable: ccr
> Type: float
> Total Size: 12864 bytes
>             3216 values
> Number of Dimensions: 2
> Dimensions and sizes:    [67] x [48]
> Coordinates:
> Number Of Attributes: 1
>   _FillValue :    1e+15
> (0)    min=-0.855062   max=0.561651
>
> Variable: sig_Dec
> Type: float
> Total Size: 12864 bytes
>             3216 values
> Number of Dimensions: 2
> Dimensions and sizes:    [67] x [48]
> Coordinates:
> warning:mpMinLonF is not a valid resource in sign_Dec_contour at this time
> warning:mpMaxLonF is not a valid resource in sign_Dec_contour at this time
> warning:mpMinLatF is not a valid resource in sign_Dec_contour at this time
> warning:mpMaxLatF is not a valid resource in sign_Dec_contour at this time
> warning:mpShapeMode is not a valid resource in sign_Dec_contour at this
> time
> warning:mpGridLineThicknessF is not a valid resource in sign_Dec_contour
> at this time
> warning:mpGridSpacingF is not a valid resource in sign_Dec_contour at this
> time
>
>
>
>
> In case full script, I have attached it.
>
> please provide some thoughts on that if you have time.
>
> Thank You
> ---
> Kunal Bali
>
>
>
> ‌
> [image: Mailtrack]
> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality1&> Sent
> with Mailtrack
> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality1&>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180530/ceb3213d/attachment.html>