<div dir="ltr"><div><div><div><div>Thanks very much. <br></div>I did not realize the loop overwrites the iteration because was in a vector so I thought it just keeps every value separately. It works now.<br></div>Thanks for pointing out the huge tapering value - I was wondering why the spectra looked really weird. <br></div>Cheers, <br></div>Ioana <br><br><div><div><div><div><div><div class="gmail_extra"><div class="gmail_quote">On 31 October 2017 at 14:42, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>The 'specx_anal' returns a <b>scalar</b> [degrees of freedom] and multiple attributes including the spectrum ( spc[*] ). <br><br> <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/specx_anal.shtml" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/Built-in/<wbr>specx_anal.shtml</a><br><br>Each iteration of your loop *overwrites* the previous iteration's attributes. Please do a 'printVarSummary' after the loop<br><br>spec=new((/nn,64/),"double") <wbr> ; wrong<span class=""><br>do ii=0,nn(0)-1<br> spec(ii,:) = specx_anal(indextotal(ii,:),d,<wbr>sm,pct)<br>end do<br><br></span>----<br></div><div>A taper of 0.7 means 70% of the values will be tapered. That is extremely large. Usually, taper ~ 0.1<br></div>Maybe<br><br>maxconf=0.95<br>minconf=0.05<span class=""><br>d = 0<br>sm = 3<br></span>pct = 0.10<br>do ii=0,nn-1<br> sdof= specx_anal(indextotal(ii,:),d,<wbr>sm,pct)<br><br> <span style="color:rgb(255,0,0)">xx=specx_ci(sdof,minconf,maxco<wbr>nf)</span><span class="m_2760879022576473304gmail-HOEnZb"><font color="#888888"><br></font></span><br></div><div> ;plot or whatever<br></div><div>end do<br><br>===<br><a href="https://www.ncl.ucar.edu/Document/Functions/Shea_util/specx_ci.shtml" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/Shea_util/<wbr>specx_ci.shtml</a><br><br></div>is prototyped as:<br><br><pre>function specx_ci (
sdof : numeric,
lowval : numeric,
highval : numeric
)<br><br><br></pre><pre>That should be<br><br>function specx_ci (
sdof[1] : numeric, <=== scalar
lowval : numeric,
highval : numeric
)</pre><div><br></div><div> Good luck<br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Tue, Oct 31, 2017 at 8:00 AM, Ioana Colfescu <span dir="ltr"><<a href="mailto:ioana.colfescu@ncas.ac.uk" target="_blank">ioana.colfescu@ncas.ac.uk</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div><div><div><div><div>Hi, <br></div>I'm trying to use specx_ci for a vector of values that are the output of specx_anal. <br></div>If I do it individually ( for every vector component) it works but when I do try to apply it for the entire vector I get the following error. <br><br>fatal:chiinv: The two input arrays must have the same dimension sizes<br><br></div>Can anyone tell me how to fix this and what's going on ? <br></div>My code is below and the part where I get the error from is in red.<br></div>Thanks. <br><div><br>begin<br>amo = addfile("<a href="http://int_ext_AMO_anom.1871-1998.nc" target="_blank">int_ext_AMO_anom.1871<wbr>-1998.nc</a>","r")<br>index0=amo->AMO_Cext<br>index1=amo->AMO_C1int<br>index2=amo->AMO_C2int<br>index3=amo->AMO_C3int<br>index4=amo->AMO_C4int<br>indextotal=(/index0,index1,ind<wbr>ex2,index3,index4/)<br>nn=dimsizes(indextotal)<br>;*************<br>d = 0<br>sm = 3<br>pct = 0.70<br><br>spec=new((/nn(0),64/),"double"<wbr>)<br>do ii=0,nn(0)-1<br>spec(ii,:) = specx_anal(indextotal(ii,:),d,<wbr>sm,pct)<br>end do<br><br>maxconf=0.95<br>minconf=0.05<br>printVarSummary(spec)<br><span style="color:rgb(255,0,0)">xx=specx_ci(spec(0,:),minconf,<wbr>maxconf)</span><span class="m_2760879022576473304HOEnZb"><font color="#888888"><br><br><br><div><div><div><div><div><div><div><div><br clear="all"><div><div class="m_2760879022576473304m_-9220217469070176495gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><span style="background-color:rgb(255,255,255)"><span style="color:rgb(61,133,198)"><span style="font-family:georgia,serif"><span><b>Ioana Colfescu </b><br></span></span></span></span></div><div><span style="background-color:rgb(255,255,255)"><span style="color:rgb(61,133,198)"><span style="font-family:georgia,serif"><span>Research Scientist <br></span></span></span></span></div><div><span style="background-color:rgb(255,255,255)"><span style="color:rgb(61,133,198)"><span style="font-family:georgia,serif"><span>National Centre for Atmospheric Science <br></span></span></span></span></div><div><span style="background-color:rgb(255,255,255)"><span style="color:rgb(61,133,198)"><span style="font-family:georgia,serif"><span>School of Earth and Environment <br></span></span></span></span></div><div><span style="background-color:rgb(243,243,243)"><span style="color:rgb(159,197,232)"><span><span style="font-family:georgia,serif"><span><span style="background-color:rgb(255,255,255)"><span style="color:rgb(61,133,198)">University of Leeds </span></span><br></span></span></span></span></span></div><div dir="ltr"><span style="color:rgb(159,197,232)"><span><span style="font-family:georgia,serif"><span><br></span></span></span></span></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div></div></div></div></div></div></div></font></span></div></div>
<br></div></div>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div></div></div></div></div></div></div>