<div dir="ltr"><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear NCL experts,</div><div><br></div><div>I have a script that computes for the mean spectrum similar to example 5 here:</div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/specx_anal.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/specx_anal.shtml</a><br></div><div><br></div><div>I am applying this to 40 summers. Here's the part that I am confused</div><div><br></div><div>spcavg = new (ntim/2, typeof(x))<br>spcavg = 0.0<br>dm = 1<br>sm = 12<br>pct = 0.10<br>r1zsum = 0.0<br>do n=0,nseg-1<br>  dof    = specx_anal(x(n,:),dm,sm,pct)<br>  spcavg = spcavg + dof@spcx<br>  r1     = dof@xlag1<br>  r1zsum = r1zsum  + 0.5*(log((1+r1)/(1-r1)))<br>end do<br></div><div><br></div><div>This part only works if I set:</div><div> <b>spcavg = new (ntim/80, typeof(x)) </b> <br></div><div><br></div><div>Am I correct to assume that spcavg is always ntim/df?</div><div>df is 2*nseg = 80.</div><div><br></div><div><br></div><div>Sincerely,</div><div>Lyndz</div></div></div></div></div></div></div></div></div></div></div></div>