<div dir="ltr"><div><div><div>Hi Allison,<br><br></div>Can you please try the attached? <br><br></div><div>I am not sure why you want to rotate.<br></div><div>The 1st two modes look very distinct.<br></div><div><br></div><div>I did not plot anything. I'll let that to you<br><br></div><div><br></div>THX<br></div>D<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 28, 2016 at 8:07 AM, Allison Serakos <span dir="ltr"><<a href="mailto:sera0021@umn.edu" target="_blank">sera0021@umn.edu</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">script name: ncl-talk_reof.ncl<div>data file name: HadISST_sst.nc.gz</div><div><br></div><div>Thank you,</div><div>Allison</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 27, 2016 at 1:49 PM, Allison Serakos <span dir="ltr"><<a href="mailto:sera0021@umn.edu" target="_blank">sera0021@umn.edu</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">
<p class="MsoNormal">Hello,</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I am using the eofunc_varimax_Wrap function to analyze global
rotated EOFs produced from the eofunc_Wrap function. However, when I run my script,
it hangs when it reaches the eofunc_varimax_Wrap function. The script does not
finish and I get no error message. </p>
<p class="MsoNormal"><span style="color:black"> </span></p>
<p class="MsoNormal"><span style="color:black">The dataset that I am using is the </span>HadISST
1.1 monthly average sea surface temperature<span style="color:black"> dataset.
It is a 3-dimensional (time, latitude, longitude) global 1</span><span style="font-family:Symbol;color:black">°</span><span style="color:black"> x 1</span><span style="font-family:Symbol;color:black">°</span><span style="color:black"> grid dataset with
1755 time values, 180 latitude, and 360 longitude. I produce annually averaged values
and spatially weighted it prior to the EOF analysis. </span></p>
<p class="MsoNormal"><span style="color:black"> </span></p>
<p class="MsoNormal">I have been using a Mac machine with NCL version 6.3.0 and a
Linux machine with NCL version 6.2.1 and the program hangs when conducting the
rotation. </p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I tried running the script with global data for different
time periods (all between 113 and 102 years of global data) on both the Mac and
Linux machines <span style="background-image:initial;background-repeat:initial">and found that the function works only for some time periods and crashes
with others. These periods differ between the two Mac and Linux machines that I
used.</span><span style="font-size:10pt;font-family:Times"> </span>I also
tried reshaping the data from a 3D matrix to a 2D matrix using the reshape
function and combining the spatial components. In this situation, the rotation
output was still the same as above for the respective time periods and
machines. Additionally, I tried reducing the spatial domain and found that when
I reduce the latitude range (anything smaller than 69S<a name="m_2278485676934046223_m_-7288369395758337668__GoBack"></a>-90N)
for any time period the script runs on both machines and the function does not
crash. <span style="font-size:10pt;font-family:Times"></span></p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I used ‘ncl –x’ to see if there is an error in my script,
but the script only crashes before the eofunc_varimax_Warp function with no
messages. When I run the script, NCL uses about 592.5MB of the 12GB of memory
on the Mac machine, so I do not believe this an issue of lack of memory space. </p>
<p class="MsoNormal"><span style="color:black"> </span></p>
<p class="MsoNormal"><span style="color:black">I am wondering if there is a known problem with
the eofunc_varimax_Wrap function? I am also wondering if you have any thoughts
as to why the function may be crashing under certain situations and not others?</span></p>
<p class="MsoNormal"><span style="color:black"> </span></p>
<p class="MsoNormal"><span style="color:black">Here is a copy of my EOF and Varimax Rotation
script:</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">;----------------------------------------------------------------
</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">;***** EOF *****</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">neval = 6</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">optEOF = True</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">optEOF@method =
"transpose"</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">optEOF@jopt = 0</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black"> </span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">eof =
eofunc_Wrap(av,neval,optEOF)</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black"> </span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">;***** VARIMAX ROTATION *****</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">eof_rot =
eofunc_varimax_Wrap(eof,-1)</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">eof_varimax_reorder(eof_rot)</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Menlo Regular";color:black">;----------------------------------------------------------------</span><span style="color:black"></span></p>
<p class="MsoNormal"><span style="color:black"> </span></p>
<p class="MsoNormal"><span style="color:black">Thanks for your help in advance. </span></p>
<p class="MsoNormal"><span style="color:black"> </span></p>
<p class="MsoNormal"><span style="color:black">Kind regards,</span></p>
<p class="MsoNormal"><span style="color:black">Allison</span></p><span><font color="#888888">
<p class="MsoNormal"> </p>-- <br><div data-smartmail="gmail_signature"><div dir="ltr">Allison Serakos<div>Graduate Research Assistant</div><div>Department of Soil, Water, and Climate</div><div>University of Minnesota</div><div>1991 Upper Buford Circle, Rm S535</div><div>St. Paul, MN 55108</div></div></div>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div data-smartmail="gmail_signature"><div dir="ltr">Allison Serakos<div>Graduate Research Assistant</div><div>Department of Soil, Water, and Climate</div><div>University of Minnesota</div><div>1991 Upper Buford Circle, Rm S535</div><div>St. Paul, MN 55108</div></div></div>
</div>
</div></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></div>