<div dir="ltr">Hello,<div>I am not sure what you mean by your question. The values coming out of eofunc are unitless, while the values coming out of eofunc_ts are the original units. A common approach to get the original units on a pattern is to normalize the PC, and regress the original data onto the normalized PC. The resulting regression plot will be in the original units per standard deviation change in the PC.</div><div><br></div><div>Hope that makes sense. If not, please respond to the ncl-talk email list only and not to me personally.</div><div>Adam</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 25, 2016 at 9:23 AM, Adv <span dir="ltr">&lt;<a href="mailto:advita6@gmail.com" target="_blank">advita6@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_quote"><div dir="ltr"><div><div>Hi ,<br></div>Could someone suggest me a way to normalise the PC values and retain the same unit?<br></div><div><br><br></div><div>Thanks <br><br><br></div>;time=41,sta=392<br><div>xx1 = dats(time | :,sta|:)     ; reorder ... eofunc want &#39;time&#39; as rightmost dimension<br>printVarSummary(xx1)<br>printMinMax(xx1,0)<br>o=dats(sta|:,time|:)<br>printVarSummary(o)<br>evecv    = eofunc_Wrap    (o,neval,False)<br>evecv_ts = eofunc_ts_Wrap (o,evecv,False)<br>print(evecv_ts(0,:))<br>printVarSummary(evecv)<br>printMinMax(evecv(0,:),0);Normalised<br>; ================================&gt;   ; SUM OF THE SQUARES<br>                                      ; IF NORMALIZED, THEY SHOULD BE 1<br>;  sumsqr = dim_sum(evecv^2)<br>;  print(&quot;sum of squares: &quot; + sumsqr)<br><br>;To retain same unit;Denormalise<br>  do ne=0,neval-1<br>     evecv(ne,:) = evecv(ne,:)*sqrt( evecv@eval(ne) )<br>  end do<br>printMinMax(evecv(0,:),0);Normalised<br>evecv_ts = eofunc_ts_Wrap (o,evecv,False)<br><br>Output of min/max<br>    min=-1693.35   max=1080.62<br></div><div><br></div></div>
</div><br></div>
<br>_______________________________________________<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/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div 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>
</div>