Re: Warning! LAPACK routine returned eigenvalue <= 0

No. The indexing is correct.
============
Ho

Variable: utrend
[snip]
Coordinates:
[snip]
              lat: [90..-90]
[snip]

I have no idea how your code went from valid lat values
to NaNs

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