[ncl-talk] computing SOI index for 100 ensembles

[Adam Phillips]

Hi Sri,
When contacting ncl-talk, please include the error message and point out
the exact line where the error is occurring. There is not enough
information here for folks on ncl-talk to assist you.
Adam

On Thu, Apr 15, 2021 at 2:29 AM Sri nandini via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hello dear ncl-users,
>
> I have some errors while computing the SOI model index over a
> 100ensembles. Could someone help me out below in my script (i only
> attached the part where the error is); in the looping over the ensemble
> part for computing the climatologies and than the standard deviation? I
> would like compute the SOI index for each ensemble and than plot out the
> timeseries.
>
>
>    T      = (fslp->slp(tStrt:tLast,:,{latT},{lonT}))
>    D      = (fslp->slp(tStrt:tLast,:,{latD},{lonD}))
>    printVarSummary(T)               ;[time | 661] x [ens | 100]
>    printVarSummary(D)               ;[time | 661] x [ens | 100]
>
>    TClm   = new ((/nmos,nens/), typeof(T), T at _FillValue)     ; monthly
> climatologies
>    DClm   = new ((/nmos,nens/), typeof(D), D at _FillValue)
>
>    do nmo=0,nmos-1                                ; loop over each month
> for current member
>     do n=0,nens-1                                      ; loop over each
> ensemble member
>       TClm(nmo,n) = avg(T(tClmStrt+nmo:tClmLast:nmos,n))  ;; convenience
> and efficiency
>       DClm(nmo,n) = avg(D(tClmStrt+nmo:tClmLast:nmos,n))
>    end do
>   end do
>    printVarSummary(TClm)
>    printVarSummary(DClm)
>    TAnom    = T
>    DAnom    = D
>
>    do nmo=0,nmos-1                                   ; anomalies
> reference clim
>       do ne=0,nens-1
>       TAnom(nmo::nmos,n) = T(nmo::nmos,n) - TClm(nmo)
>       DAnom(nmo::nmos,n) = D(nmo::nmos,n) - DClm(nmo)
>       end do
>    end do
>    printVarSummary(TAnom)
>    printVarSummary(DAnom)
>
>    TAnomStd   = stddev(TAnom(tClmStrt:tClmLast))     ; stddev of
> anomalies over clStrt & clLast
>    DAnomStd   = stddev(DAnom(tClmStrt:tClmLast))
>    printVarSummary(TAnomStd)
> printVarSummary(DAnomStd) ; signal and noise
>    soi_signal = (TAnom/TAnomStd) - (DAnom/DAnomStd)  ; (ntim)
>    copy_VarCoords(TAnom, soi_signal)
>    soi_signal at long_name = "SOI: 20th Century Reanalysis:
> "+yrStrt+"-"+yrLast
>    printVarSummary(soi_signal)
>
> --
> Dr. Sri, Nandini-Weiss
> Research Associate
>
> Universität Hamburg
> Center for Earth System Research and Sustainability (CEN)
> Cluster of Excellence 'Climate, Climatic Change, and Society' (CLICCS)
>
> Bundesstrasse 53, 20146 Hamburg
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk



-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210415/696925bb/attachment.html>