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 12:56:40 2013
This archive was generated by hypermail 2.1.8 : Tue Jun 11 2013 - 12:03:58 MDT