# Re: Warning! LAPACK routine returned eigenvalue <= 0

From: Erika Folova <e.folova_at_nyahnyahspammersnyahnyah>
Date: Mon Jun 10 2013 - 13:35:36 MDT

Hallo,

Here some information:

ncl 43> print("X: min="+min(u)+" max="+max(u))
(0) X: min=-82.4277 max=156.545

then, you are correct that
utrend&lat( 0) = 90.
.
.
.
utrend&lat(72) = -90.

So what should I do, do I have to reverse the pole points?

Thank you

On Tue, Jun 11, 2013 at 4:22 AM, Dennis Shea <shea@ucar.edu> wrote:

> print( utrend&lat )
>
> Looks like the pole points are wrong. Maybe ..
>
> utrend&lat( 0) = 90.
> utrend&lat(72) = -90.
>
>
> clat = utrend&lat
> clat = sqrt( cos(rad*clat) )
>
> Also, if X is the original (unweighted) array
>
> print("X: min="+min(X)+" max="+max)X))
>
> If there are NaNs, the user must eliminate.
>
>
>
> On 6/10/13 1:13 PM, Erika Folova wrote:
>
>> I think the problem is in these lines:
>>
>> 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*
>>>
>>>
>>>
>>> Erika,
>>>
>>>
>>>
>>>
>>>
>>
>>
>> ______________________________**_________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<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:36:57 2013

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