<div dir="ltr">Ok, I remembered something about the correlation matrices. In version 6.1.2 there was a bug in <div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/covcorm.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/covcorm.shtml</a><br></div><div>where the options behaved the opposite of what they were supposed to do. Maybe that's the issue here.</div><div>I would recommend downloading the POP driver script</div><div><a href="http://ncl.ucar.edu/Applications/Scripts/PrnOscPat_driver.ncl">http://ncl.ucar.edu/Applications/Scripts/PrnOscPat_driver.ncl</a><br></div><div>save it under a different name and load that script at the beginning of your code. Instead of calling the built-in version, call the new version in your code. </div><div><br></div><div>Before you do that you can change the options for covcorm in line 25 in the new POP driver script. I would also include a few printVarSummarys throughout the code so that you can see if your array is suddenly 1D instead of 2D.</div><div><br></div><div>Hope this helps,</div><div>Maria</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 8, 2015 at 3:33 PM, Maria Gehne - NOAA Affiliate <span dir="ltr"><<a href="mailto:maria.gehne@noaa.gov" target="_blank">maria.gehne@noaa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Lena,<div>I can look into it tonight or tomorrow morning. I'm traveling right now so it's hard to look at your code and the pop routines. I wrote most of the original code so hopefully I can help you.</div><span class="HOEnZb"><font color="#888888"><div><br></div></font></span><div><span class="HOEnZb"><font color="#888888">Maria</font></span><div><div class="h5"><br><br>On Thursday, January 8, 2015, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Technically,the 'Principal Oscillation Pattern' is not a supported function. However, it will be looked at in th near future: today or tomorrow.<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 8, 2015 at 11:22 AM, Mary Haley <span dir="ltr"><<a>haley@ucar.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"><div class="gmail_default" style="font-size:small">Lena,</div><div class="gmail_default" style="font-size:small"><br>A friendly tip: if you reply to an email thread, it's helpful to keep all the original emails in the body of the message. This particular email didn't include the original email with the error message, so I had to go back a few emails to find it.<br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The error is referring to an "inverse_matrix" function, but your script is not calling this function directly, so it's difficult for us to debug your script without knowing what other function is referencing inverse_matrix.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Usually you can debug these things on your own by using "printVarSummary", which will tell you what the dimensionality and size of your data is.<br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Find the location in your script where you think inverse_matrix is being called by another function, and then insert some "printVarSummary" calls to look at the data before you call the offending function. This will hopefully shed some light on why your data is not the correct dimensionality.</div><span><font color="#888888"><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div></font></span><span><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 8, 2015 at 6:43 AM, Lena Frey <span dir="ltr"><<a>Frey_L@web.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div><div style="font-family:Verdana;font-size:12px"><div>Hi again!</div>
<div> </div>
<div>Thanks to Karin. The problem is, that I don't know why I get this error message.</div>
<div>I have used the ncl script from the ncl website to calculate the POPs for the mean sea level pressure.</div>
<div>I am using NCEP reanalysis data. The mean sea level pressure has the dimensions time, lat, lon.</div>
<div>So, following the ncl script it should be 2D before calculating the inverse matrix.</div>
<div> </div>
<div>You can find attached my ncl script.</div>
<div> </div>
<div>Thanks,</div><span><font color="#888888">
<div> </div>
<div>Lena</div></font></span></div></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></span></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div></div></div>
</blockquote></div><br></div>