<div dir="ltr"><div><div><div><div>I think the only NCL examples of <b>reconstruction</b> are based on a simple text book:<br><br> <a href="https://www.ncl.ucar.edu/Applications/eof.shtml">https://www.ncl.ucar.edu/Applications/eof.shtml</a><br><br></div><div>Unweighted:<br></div> Example: eof_0.ncl and eof_0_640 [use 6.4.0 if available]<br><br></div>Weighted:<br> Example: eof_0a.ncl [modify for 6.4.0 if available]<br><br>---<br></div>FYI: NCL's EOF works on anomalies<br><br></div>Good luck<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Dec 15, 2017 at 2:52 AM, 蘇容瑩 <span dir="ltr"><<a href="mailto:kikyourobin@g.ncu.edu.tw" target="_blank">kikyourobin@g.ncu.edu.tw</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"><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">Hello all,</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">I input omega data(time, lev, lat, lon) and do EOF for every grid point </font><font face="microsoft jhenghei, sans-serif">to get characteristic vertical structure of </font><span style="font-family:arial,sans-serif">omega.</span> </pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">I used eofunc to get 3 modes of eigenvectors.</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">However, I met a problem, when I reconstructed mode 1(w1) and mode 2 (w2).</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">I tested 3 grid points and compared reconstructed omega and original omega.</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">grid 1 : 0N, 150E in western Pacific </font></pre><pre><font style="color:rgb(0,0,0)" face="microsoft jhenghei, sans-serif">grid 2 : 10N, 120W in eastern Pacific</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">grid 3 : 10S, 110W in eastern Pacific cold tongue</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">Result in grid 1 & 2 match my intuition. w1 >w2</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">However, result in grid 3 against my intuition. w2>w1</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">After examining data, I found ts_mean of mode 2 larger then mode 1 so that w2 larger than w1. </font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">I attached my code for EOF reconstruction below.</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">Is there anything wrong with my conceptions or coding process?</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif">Thanks for help!</font></pre><pre style="color:rgb(0,0,0)"><font face="microsoft jhenghei, sans-serif"><br></font></pre><pre style="color:rgb(0,0,0)"><br></pre><pre><font face="microsoft jhenghei, sans-serif"><font color="#000000"> do mode = 0, 2</font></font></pre><pre><font face="microsoft jhenghei, sans-serif"><font color="#000000"> w_re = new((/nSeason,klev,nlat,mlon/)<wbr>,float,msg_val)
w_re = 0.
dimName = (/"time", "plev", "lat", "lon"/)
do dim = 0, 3
w_re!dim = dimName(dim)
if (dim .ne. 0) then
w_re&$dimName(dim)$ = w_eof&$dimName(dim)$
end if
end do
do i = 0, nlat-1
do j = 0, mlon-1
do time = 0, nSeason-1
if (num(ismissing(w_eof(mode,:,i,<wbr>j))) .le. 5) then
w_re(time,:,i,j) = w_eof(mode,:,i,j)*(w_ts(mode,t<wbr>ime,i,j) + w_ts_mean(mode,i,j))
end if
end do
end do
end do
</font></font></pre><span class="HOEnZb"><font color="#888888"><div><font face="microsoft jhenghei, sans-serif"><br></font></div><div><br></div>-- <br><div class="m_-512435487528678876gmail-m_1838311307133038767m_-4311819739578672805gmail-m_6749479378614629118gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font face="Arial, Helvetica, sans-serif" color="#000000"><span style="font-size:13.3333px;background-color:rgb(255,255,255)"><br></span></font></div><div dir="ltr"><font face="Arial, Helvetica, sans-serif" color="#000000"><span style="font-size:13.3333px;background-color:rgb(255,255,255)">-------------------</span></font></div><div dir="ltr"><span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:13.3333px;background-color:rgb(255,255,255)">Rung-Ying, Su</span></div><div dir="ltr"><span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:13.3333px;background-color:rgb(238,238,238)"><div dir="ltr"><span style="font-size:13.3333px">Department of Atmospheric Sciences </span><br></div><div dir="ltr">National Central University (<a href="http://www.ncu.edu.tw/" target="_blank">http://www.ncu.edu.tw/</a>)</div><div dir="ltr"><br></div></span></div><div dir="ltr"><span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:13.3333px;background-color:rgb(238,238,238)"><br></span></div><div dir="ltr"><span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:13.3333px;background-color:rgb(238,238,238)"><br></span></div><div dir="ltr"><span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:13.3333px;background-color:rgb(238,238,238)"><br></span></div><div dir="ltr"><span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:13.3333px;background-color:rgb(238,238,238)"><br></span></div></div></div></div></div></div></div></div>
</font></span></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></div>