Re: "Eigenvectors failed to converge" with eofunc

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon, 02 Mar 2009 08:47:41 -0700

My speculation is that you have some grid points will all
missing values.

Aideen Foley wrote:
> Hi everyone,
> I am just wondering if anyone has seen this error message when using eofunc?
> I am working with a lat x lon x time temperature data set, where each
> entry on the t axis is a January from a 30 year series.
> (So basically I extracted 30 of each month from the main dataset
> seperately and saved them as new files).
>
> Here is my script, using the instructions from example 2 in the
> documentation to reorder axes.
>
> /eof_jan = addfile("meantemp_jan_EOF.nc","r")
> jan1 = eof_jan->V5(:,:,:)/
>
> / jan1!0 = "time" ; name dimensions if not already
> done
> jan1!1 = "lat" ; must be named to reorder
> jan1!2 = "lon" /
> / neval = 3 ; calculate all EOFs
> ev = eofunc(jan1(lat|:,lon|:,time|:),neval,False)
> /
>
> /jan_output = addfile("jan_out.nc <http://jan_out.nc>","c") /
> //
> /jan_output->ev = jan1/
> //
> //
> The error message says:
>
> /sub drveof: ier,jopt= 10 0 returned from vcmssm/crmssm/
> /warning:eofunc: 10 eigenvectors failed to converge/
>
> Any idea how to fix this?
>
> Thanks,
>
> Aideen Foley
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Mar 02 2009 - 08:47:41 MST

This archive was generated by hypermail 2.2.0 : Mon Mar 02 2009 - 16:45:42 MST