<div dir="ltr">Hi Sri,<div>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.</div><div>Adam</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 15, 2021 at 2:29 AM Sri nandini via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello dear ncl-users,<br>
<br>
I have some errors while computing the SOI model index over a <br>
100ensembles. Could someone help me out below in my script (i only <br>
attached the part where the error is); in the looping over the ensemble <br>
part for computing the climatologies and than the standard deviation? I <br>
would like compute the SOI index for each ensemble and than plot out the <br>
timeseries.<br>
<br>
<br>
   T      = (fslp->slp(tStrt:tLast,:,{latT},{lonT}))<br>
   D      = (fslp->slp(tStrt:tLast,:,{latD},{lonD}))<br>
   printVarSummary(T)               ;[time | 661] x [ens | 100]<br>
   printVarSummary(D)               ;[time | 661] x [ens | 100]<br>
<br>
   TClm   = new ((/nmos,nens/), typeof(T), T@_FillValue)     ; monthly <br>
climatologies<br>
   DClm   = new ((/nmos,nens/), typeof(D), D@_FillValue)<br>
<br>
   do nmo=0,nmos-1                                ; loop over each month <br>
for current member<br>
    do n=0,nens-1                                      ; loop over each <br>
ensemble member<br>
      TClm(nmo,n) = avg(T(tClmStrt+nmo:tClmLast:nmos,n))  ;; convenience <br>
and efficiency<br>
      DClm(nmo,n) = avg(D(tClmStrt+nmo:tClmLast:nmos,n))<br>
   end do<br>
  end do<br>
   printVarSummary(TClm)<br>
   printVarSummary(DClm)<br>
   TAnom    = T<br>
   DAnom    = D<br>
<br>
   do nmo=0,nmos-1                                   ; anomalies <br>
reference clim<br>
      do ne=0,nens-1<br>
      TAnom(nmo::nmos,n) = T(nmo::nmos,n) - TClm(nmo)<br>
      DAnom(nmo::nmos,n) = D(nmo::nmos,n) - DClm(nmo)<br>
      end do<br>
   end do<br>
   printVarSummary(TAnom)<br>
   printVarSummary(DAnom)<br>
<br>
   TAnomStd   = stddev(TAnom(tClmStrt:tClmLast))     ; stddev of <br>
anomalies over clStrt & clLast<br>
   DAnomStd   = stddev(DAnom(tClmStrt:tClmLast))<br>
   printVarSummary(TAnomStd)<br>
printVarSummary(DAnomStd) ; signal and noise<br>
   soi_signal = (TAnom/TAnomStd) - (DAnom/DAnomStd)  ; (ntim)<br>
   copy_VarCoords(TAnom, soi_signal)<br>
   soi_signal@long_name = "SOI: 20th Century Reanalysis: "+yrStrt+"-"+yrLast<br>
   printVarSummary(soi_signal)<br>
<br>
-- <br>
Dr. Sri, Nandini-Weiss<br>
Research Associate<br>
<br>
Universität Hamburg<br>
Center for Earth System Research and Sustainability (CEN)<br>
Cluster of Excellence 'Climate, Climatic Change, and Society' (CLICCS)<br>
<br>
Bundesstrasse 53, 20146 Hamburg<br>
<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist,  </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a>   </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>