Re: Warning! LAPACK routine returned eigenvalue <= 0

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Jun 10 2013 - 13:15:57 MDT

On 6/10/13 1:13 PM, Erika Folova wrote:
> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 10 13:15:44 2013

This archive was generated by hypermail 2.1.8 : Tue Jun 11 2013 - 12:03:57 MDT