Re: Warning! LAPACK routine returned eigenvalue <= 0

I think the problem is in these lines:

  rad = 4.*atan(1.)/180.
  clat = utrend&lat
  clat = sqrt( cos(rad*clat) )

Variable: clat
Type: float
Total Size: 292 bytes
            73 values
Number of Dimensions: 1
Dimensions and sizes: [lat | 73]
Coordinates:
Number Of Attributes: 2
  long_name : latitude
  units : degrees_north
*(0) -nan*
(1) 0.2088527
(2) 0.2952216
(3) 0.3612842
(4) 0.4167112
(5) 0.4652308
(6) 0.5087426
(7) 0.5483665
(8) 0.5848249
(9) 0.6186141
(10) 0.6500909
(11) 0.6795208
(12) 0.7071068
(13) 0.7330072
(14) 0.7573483
(15) 0.7802317
(16) 0.8017403
(17) 0.8219429
(18) 0.8408964
(19) 0.8586485
(20) 0.8752397
(21) 0.8907039
(22) 0.9050702
(23) 0.9183635
(24) 0.9306049
(25) 0.9418126
(26) 0.952002
(27) 0.9611865
(28) 0.9693775
(29) 0.9765844
(30) 0.9828153
(31) 0.9880769
(32) 0.9923748
(33) 0.9957132
(34) 0.9980956
(35) 0.999524
(36) 1
(37) 0.999524
(38) 0.9980956
(39) 0.9957132
(40) 0.9923748
(41) 0.9880769
(42) 0.9828153
(43) 0.9765844
(44) 0.9693775
(45) 0.9611865
(46) 0.952002
(47) 0.9418126
(48) 0.9306049
(49) 0.9183635
(50) 0.9050702
(51) 0.8907039
(52) 0.8752397
(53) 0.8586485
(54) 0.8408964
(55) 0.8219429
(56) 0.8017403
(57) 0.7802317
(58) 0.7573483
(59) 0.7330072
(60) 0.7071068
(61) 0.6795208
(62) 0.6500909
(63) 0.6186141
(64) 0.5848249
(65) 0.5483665
(66) 0.5087426
(67) 0.4652308
(68) 0.4167112
(69) 0.3612842
(70) 0.2952216
(71) 0.2088527
*(72) -nan*

Variable: utrend
Type: float
Total Size: 510630912 bytes
            127657728 values
Number of Dimensions: 4
Dimensions and sizes: [time | 528] x [level | 23] x [lat | 73] x [lon | 144]
Coordinates:
            time: [508416..893376]
            level: [1..1000]
            lat: [90..-90]
            lon: [ 0..357.5]

But I dont know why I got "NaN".

Thank you,

On Tue, Jun 11, 2013 at 3:56 AM, Erika Folova <e.folova@gmail.com> wrote:

> Hallo everyone,
>
> would you please to help me, I got some errors when I tried to calculate
> eof
> for variable (*level, lat*, time):
>
> Warning! LAPACK routine returned eigenvalue <= 0.
> Setting it to zero...
> warning:ContourPlotInitialize: no valid values in scalar field;
> ContourPlot not possible:[errno=1101]
> warning:LlSetValues: invalid range for Y Axis log mode; setting trYLog off
>
>
> Part of my codes:
> [SNIP]
>
> wx = xw (level|:,lat|:,lon|:,time|:)
> delete(xw)
> x= wx (:,:,{0},:)
> printVarSummary(x)
>
>
> eof = eofunc_Wrap(x, neof, optEOF)
> eof = -1*eof
> eof_ts = eofunc_ts_Wrap (x, eof, optETS)
>
> printVarSummary( eof ) ; examine EOF
> variables
> printVarSummary( eof_ts )
>
> eof_ts = dim_standardize_n( eof_ts, 0, 1) ; normalize
>
> eof_regres = eof ; regress
> do ne=0,neof-1
> eof_regres(ne,:,:) = (/ regCoef(eof_ts(ne,:),
> x(level|:,lat|:,time|:)) /)
> end do
>
> ----
>
> plot = gsn_csm_pres_hgt(wks, eof_regres(0,:,:), res ) ; place holder
>
> end
>
>
> ==============================
>
> Variable: x
> Type: float
> Total Size: 3546048 bytes
> 886512 values
> Number of Dimensions: 3
> Dimensions and sizes: [level | 23] x [lat | 73] x [time | 528]
> Coordinates:
> level: [1..1000]
> lat: [90..-90]
> time: [508416..893376]
> Number Of Attributes: 7
> lon : 0
> missing_value : -32767
> long_name : Wgt: U velocity
> units : m s**-1
> missing_value_original : -32767
> _FillValue : -32767
> _FillValue_original : -32767
> Warning! LAPACK routine returned eigenvalue <= 0.
> Setting it to zero...
>
> Variable: eof
> Type: float
> Total Size: 6716 bytes
> 1679 values
> Number of Dimensions: 3
> Dimensions and sizes: [evn | 1] x [level | 23] x [lat | 73]
> Coordinates:
> evn: [1..1]
> level: [1..1000]
> lat: [90..-90]
> Number Of Attributes: 7
> eval_transpose : 0
> eval : 0
> pcvar : -nan
> matrix : covariance
> method : transpose
> _FillValue : -32767
> long_name : EOF: Wgt: U velocity
>
> Variable: eof_ts
> Type: float
> Total Size: 2112 bytes
> 528 values
> Number of Dimensions: 2
> Dimensions and sizes: [evn | 1] x [time | 528]
> Coordinates:
> evn: [1..1]
> time: [508416..893376]
> Number Of Attributes: 4
> ts_mean : -nan
> matrix : covariance
> _FillValue : -32767
> long_name : EOF: Amplitude: Wgt: U velocity
> *ncl 77> print(max(eof_regres))*
> *(0) 9.96921e+36*
> *ncl 78> print(max(eof)) *
> *(0) -nan*
> *ncl 79> print(max(eof_ts))*
> *(0) -nan*
> *ncl 80> print(min(eof)) *
> *(0) -nan*
>
>
> Thanks for your help.
>
> Erika,
>
>
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 10 13:13:37 2013