<div>Dear NCL people,</div><div><br></div><div>I'm struggling at computing the EOFs of a zonal wind stress field. I'm always getting this warning from the eofunc functions:</div><div><br></div><div><span style="font-family:monospace"><span style="color:#000000;background-color:#ffffff;"> Warning! LAPACK routine returned eigenvalue <= 0.
</span><br> Setting it to zero...<br>
<br></span>I've seen that the original LAPACK's dspevx function sets a 1.e-6 threshold for computing EOFs; if there is a value below that threshold, zero EOFs are returned (https://github.com/yyr/ncl/blob/master/ni/src/lib/nfpfort/eof_scripps.f). I'm using the where function to set all values below 1.e-6 as fill values (something like taux = where(taux.le.1.e-6,taux@_FillValue,taux)) before the eofunc function is called, yet I'm getting the same warning and the same empty EOFs. Is there anything else I'm missing?</div><div><br></div><div>Thanks,</div><div>Giorgio<br><span style="font-family:monospace"></span></div>