<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 face="microsoft jhenghei, sans-serif" style="color:rgb(0,0,0)">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><div><font face="microsoft jhenghei, sans-serif"><br></font></div><div><br></div>-- <br><div class="gmail-m_1838311307133038767m_-4311819739578672805gmail-m_6749479378614629118gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font color="#000000" face="Arial, Helvetica, sans-serif"><span style="font-size:13.3333px;background-color:rgb(255,255,255)"><br></span></font></div><div dir="ltr"><font color="#000000" face="Arial, Helvetica, sans-serif"><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>
</div>