<div dir="ltr">That's a great catch. The way I handle it is:<div><br></div><div><div><font face="monospace, monospace"> ; calculate the PCs</font></div><div><font face="monospace, monospace"> retVal = eofunc_ts_Wrap( data(var|:,lon|:,time|:), eof, False )</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> do e = 0, nEof-1</font></div><div><font face="monospace, monospace"> retVal(e,:) = ( retVal(e,:) + retVal@ts_mean(e) ) / pcStddev(e)</font></div><div><font face="monospace, monospace"> end do</font></div></div><div><br></div><div><font face="arial, helvetica, sans-serif">Where pcStddev is the standard deviation calculated from the original PCs over a 30+ year period.</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">It would be great if eofunc_ts had an option to not subtract the mean automatically.</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Carl</font></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jul 25, 2017 at 12:34 AM, Eun-Pa Lim <span dir="ltr"><<a href="mailto:eun-pa.lim@bom.gov.au" target="_blank">eun-pa.lim@bom.gov.au</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" style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<div dir="ltr" style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<div>
<div dir="ltr" style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="font-size:12pt">Dear NCL developers</span>
<div>
<p><br>
</p>
<p>My colleagues Harry Hendon and Abhik Santra at the Bureau of Meteorology want to report an issue regarding the computation of the RMM1 and RMM2 time series (eof_ts_cdata) and the expansion coefficients of olr, u200 and u850 (ceof_ts) in mjoclivar_14.ncl. </p>
<pre style="word-wrap:break-word;white-space:pre-wrap">> eof_ts_cdata = eofunc_ts(cdata,eof_cdata,<wbr>False) <br></pre>
<p>Harry and Abhik were trying to replace cdata with forecast data of a particular year to see if forecast amplitudes of RMM1 and RMM2 are realistic or not. </p>
<p>Because eofunc_ts_Wrap removes the time mean of the input data by default (i.e. 1-year time mean of forecast data for their case), their forecast MJO amplitudes were systematically biased. <br>
</p>
<p><br>
</p>
<p>This problem will occur when people try to compute the MJO amplitues with a subset of a dataset.<br>
</p>
<p><br>
</p>
<p>Could you please add a comment where eof_ts_cdata and ceof_ts appear in mjoclivar_14.ncl that users should compute the expansion coefficients manually rather than using eofunc_ts_Wrap if the removal of the time mean of the input data is not intended/desirable?<br>
</p>
<p><br>
</p>
<p>Please let me know if what I've explained is not clear. <br>
</p>
<p><br>
</p>
<p>Many thanks!<br>
</p>
<p><br>
</p>
<p>Eun-Pa<br>
</p>
<p><br>
</p>
<p>Research Scientist<br>
</p>
<p>The Australian Bureau of Meteorology<br>
</p>
<p>GPO Box 1289, Melbourne, VIC 3000, AUS<br>
</p>
<p><br>
</p>
</div>
</div>
</div>
</div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><br><table border="0" cellpadding="2" cellspacing="2" style="max-width:100%;border-collapse:collapse;border-spacing:0px;color:rgb(51,51,51);font-size:14px;border:3px solid rgb(170,170,170);font-family:Times;line-height:12px;background-color:transparent"><tbody><tr><td align="center" height="71" width="71"><span style="font-size:11px"><span style="font-family:arial,helvetica,sans-serif"><a href="https://ncics.org/" target="_blank"><img src="https://ncics.org/wp-content/themes/ncics/img/NCICS_logo-131x131.png" alt=" "><br></a></span></span><br><font size="4"><b style="text-align:start;color:rgb(38,58,143);font-family:arial,helvetica,sans-serif"> <a href="http://www.cyclonecenter.org/" style="color:rgb(38,58,143);text-decoration:none" target="_blank">Cyclone</a></b><a href="http://www.cyclonecenter.org/" style="text-align:start;color:rgb(38,58,143);font-family:arial,helvetica,sans-serif;text-decoration:none" target="_blank">Center.org</a> </font></td><td valign="top"><span style="font-size:11px"><span style="font-family:arial,helvetica,sans-serif"><b>Carl J. Schreck III, PhD</b><br></span></span><span style="font-size:11px"><span style="font-family:arial,helvetica,sans-serif"><b>Research Scholar / Research </b></span></span><b style="font-family:arial,helvetica,sans-serif;font-size:11px">Assistant Professor</b><span style="font-size:11px"><span style="font-family:arial,helvetica,sans-serif"><br></span></span><a href="http://www.cicsnc.org/" style="font-family:arial,helvetica,sans-serif;font-size:11px;color:rgb(38,58,143);text-decoration:none;font-weight:bold" target="_blank">Cooperative Institute for Climate and Satellites NC (CICS</a><b><a href="http://www.cicsnc.org/" style="font-family:arial,helvetica,sans-serif;font-size:11px;color:rgb(38,58,143);text-decoration:none" target="_blank">-NC)</a> /</b><br><span style="font-size:11px"><span style="font-family:arial,helvetica,sans-serif"><a href="https://meas.sciences.ncsu.edu/" style="color:rgb(38,58,143);text-decoration:none;font-weight:bold" target="_blank">Department of Marine, Earth, and Atmospheric Sciences (MEAS)</a> <br><a href="http://ncsu.edu/" style="color:rgb(38,58,143);text-decoration:none;font-weight:bold" target="_blank">North Carolina State University</a><br><a href="http://ncei.noaa.gov/" style="color:rgb(38,58,143);text-decoration:none;font-weight:bold" target="_blank">NOAA National Centers for Environmental Information (NCEI)</a><br>151 Patton Ave, Asheville, NC 28801<br>e: </span></span><a href="mailto:carl_schreck@ncsu.edu" style="color:rgb(38,58,143);font-family:arial,helvetica,sans-serif;font-size:11px;text-decoration:none;font-weight:bold" target="_blank">carl_schreck@ncsu.edu</a><br style="font-family:arial,helvetica,sans-serif;font-size:11px"><span style="font-size:11px"><span style="font-family:arial,helvetica,sans-serif">o: +1 828 257 3140<br></span></span><span style="font-family:arial,helvetica,sans-serif;font-size:11px">c: +1 828 484 1702<br></span><a href="http://scholar.google.com/citations?hl=en&user=th8ONEcAAAAJ&view_op=list_works&sortby=pubdate" style="color:rgb(38,58,143);font-family:arial,helvetica,sans-serif;font-size:11px;text-decoration:none;font-weight:bold" target="_blank">Publications</a><br><a href="https://ncics.org/mjo/" style="color:rgb(38,58,143);font-family:arial,helvetica,sans-serif;font-size:11px;text-decoration:none;font-weight:bold" target="_blank">ncics.org/mjo</a><br style="font-family:arial,helvetica,sans-serif;font-size:11px"></td></tr></tbody></table></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div>