[ncl-talk] Data file for example bar_2.ncl

[Dennis Shea]

The 'soi.nc'  file is no longer available. However, you can use the *SOI.nc
*file:
[1] use the *Trenberth 11-pt filter* and
[2] Use a decadal filter

*See attached*
====

 f = addfile ("SOI.nc", "r")      ; add file
 date  = f->time                     ; YYYYMM
 dsoi  = f->SOI_DARWIN     ; Darwin SOI

;---Trenberth's 11-point (month) smoother effectively removes fluctuations
;   with periods of less than 8 months but includes all others.
;   At 24 months 80% of the variance is retained.
;   Use reflective boundaries (kopt=1) to fill out plot

  kopt  = 1
  wgtk  = (/ 0.0270, 0.05856, 0.09030, 0.11742, 0.13567, \
             0.1421, 0.13567, 0.11742, 0.09030, 0.05856, 0.027 /)
  dsoik = *wgt_runave_n_Wrap *
<http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_runave_n_Wrap.shtml>(dsoi,
wgtk, kopt, 0)
  *printVarSummary*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml>
(dsoik)
  *printMinMax*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/printMinMax.shtml>
(dsoik,0)

;---Decadal filter

  ihp   = 0                                ; low_pass
  sigma = 1.0                           ; Lanczos sigma
  nWgt  = 121                          ; loose 60 (nWgt/2) months each end
  fca   = 1./120.                       ; decadal
  wgtd  = *filwgts_lanczos
<http://www.ncl.ucar.edu/Document/Functions/Built-in/filwgts_lanczos.shtml>
*(nWgt, ihp, fca, -999., sigma )   ; generate eights
  dsoid = *wgt_runave_n_Wrap*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_runave_n_Wrap.shtml>
( dsoi, wgtd, 0, 0 )             ; apply weights
*  printVarSummary*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml>
(dsoid)
  *printMinMax*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/printMinMax.shtml>
(dsoid,0)

On Fri, Mar 22, 2019 at 11:09 PM Atul Saini <atulsainimail at gmail.com> wrote:

> Hello all,
>              Data file for example bar_2.ncl is required (soi.nc). We
> have SOI.nc at https://www.ncl.ucar.edu/Applications/Data/ but it doesn't
> have the variables used in the example script.
>
>
> Regards,
>
>
>
> _______________________________________________
> 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/20190324/3639aced/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: bar_2.ncl_talk
Type: application/octet-stream
Size: 3325 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190324/3639aced/attachment.obj>