Hi, First to clarify, the EOF functions are internal Ferret functions, not scripts. The EOF demo script calls all three EOF functions. It first lists the result of the STAT function which will give you the number of eof's returned, percent variance each explains and the eigenvalues. That output is listed here: Next, after you hit RETURN, it plots the spatial variation, which is the variable defined as let exy = EOF_SPACE(sample_function, 0.1) Here the first eigenfunction is shown as the result at L=1, and the second eigenfunction is the result at L=2: SET VIEW vll; CONTOUR/L=1/TITLE="EOF 1" exy ! this is the first spatial eigenfunctionTo save the first eigenfunction, you would save the variable exy[L=1], and likewise the second is exy[L=2]. In the figure, the upper frames are the original function, for comparison, and the lower ones are the eigenfunctions. Finally the script computes and plots the time amplitude functions. The demo goes on to compute eof's using the example coads_climatology dataset. You could do something similar using a different dataset or a different region. Call EOF_TFUNC to get the time amplitude function and EOF_SPACE to get the spatial eigenfunctions. On 8/4/2010 2:53 AM, medepalli santha raju wrote:
|