[ncl-talk] fatal:Subscript out of range, error in subscript #0 fatal:An error occurred reading unknown fatal:["Execute.c":8635]:Execute: Error occurred at or near line 201
Rahpeni Fajarianti
rahpenifajarianti at gmail.com
Mon Dec 16 04:26:14 MST 2019
Hello NCL user. I wanna make mjo clivar 14, but I got some error.
This is my script:
ncl 0> ;******************************************************
ncl 1> ;
ncl 2> ; mjoclivar_14.ncl
ncl 3> ;
ncl 4> ;***********************************************************
ncl 5> ; Combined EOFs
ncl 6> ; Latest Update: July, 2016: Eun-Pa Lim; Bureau of Meteorology,
Australia
ncl 7> ;***********************************************************
ncl 8> ;;
ncl 9> ;; The following are automatically loaded from 6.2.0 onward
ncl 10> ;;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
ncl 11> ;;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
ncl 12> ;;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
ncl 13>
ncl 14> undef("read_rename")
ncl 15> function read_rename(f[1]:file, varName[1]:string \
ncl 15> ,iStrt[1]:integer, iLast[1]:integer \
ncl 15> ,latS[1]:numeric , latN[1]:numeric )
ncl 16> ; Utility to force specific named dimensions
ncl 17> ; This is done for historical reasons (convenience)
ncl 18> begin
ncl 19> work = f->$varName$(iStrt:iLast,{latS:latN},:) ;
(time,lat,lon)
ncl 20> work!0 = "time" ; CAM model
names
ncl 21> work!1 = "lat"
ncl 22> work!2 = "lon"
ncl 23> return(work)
ncl 24> end
ncl 25> ; =========================> MAIN <==============================
ncl 26> begin
ncl 27> neof = 2
ncl 28>
/ncl 29> latS = -15
ncl 30> latN = 15
ncl 31>
ncl 32> ymdStrt = 20180101 ; start yyyymmdd
ncl 33> ymdLast = 20181231 ; last
ncl 34>
ncl 35> yrStrt = ymdStrt/10000
ncl 36> yrLast = ymdLast/10000
ncl 37>
ncl 38> pltDir = "/home/peni/" ; plot
directory
ncl 39> pltType = "png"
ncl 40> pltName = "mjoclivar"
ncl 41>
ncl 42> diri = "/home/peni/" ; input
directory
ncl 43>
ncl 44> filolr = "anomolr.nc"
ncl 45> filu200 = "uanom2.nc"
ncl 46> filu850 = "uanom.nc"
ncl 47>
ncl 48> ;************************************************
ncl 49> ; create BandPass Filter
ncl 50> ;************************************************
ncl 51> ihp = 2 ; bpf=>band pass filter
rncl 52> nWgt = 201
ncl 53> sigma = 1.0 ; Lanczos sigma
r ncl 54> fca = 1./100.
ncl 55> fcb = 1./20.
"olr", "u850", "u200"
ceof = newncl 56> wgt = filwgts_lanczos (nWgt, ihp, fca, fcb,
sigma )
ncl 57>
ncl 58> ;***********************************************************
ncl 59> ; Find the indices corresponding to the start/end times
ncl 60> ;***********************************************************
ncl 61> f = addfile (diri+filolr , "r")
ncl 62> TIME = f->time ; days since ...
ncl 63>
ncl 64> YMD = cd_calendar(TIME, -2) ; entire (time,6)
ncl 65>
tncl 66> iStrt = ind(YMD.eq.ymdStrt) ; index start
ncl 67> iLast = ind(YMD.eq.ymdLast) ; index last
ncl 68> delete([/ TIME, YMD /])
ncl 69>
nncl 70> ;***********************************************************
ncl 71> ; Read anomalies
ncl 72> ;***********************************************************
ncl 73>
ncl 74> work = read_rename(f,"anomolr",iStrt,iLast,latS,latN) ;
(time,lat,lon)
ncl 75> OLR = dim_avg_n_Wrap(work, 1) ;
(time,lon)
ncl 76>
ncl 77> f = addfile (diri+filu850 , "r")
ncl 78> work = read_rename(f,"uanom",iStrt,iLast,latS,latN) ;
(time,lat,lon)
ncl 79> U850 = dim_avg_n_Wrap(work, 1) ; (time,lon)
ncl 80>
ncl 81> f = addfile (diri+filu200 , "r")
ncl 82> work = read_rename(f,"uanom2",iStrt,iLast,latS,latN) ;
(time,lat,lon)
ncl 83> U200 = dim_avg_n_Wrap(work, 1) ; (time,lon)
ncl 84>
^ncl 85> dimw = dimsizes( work )
ncl 86> ntim = dimw(0)
ncl 87> nlat = dimw(1)
ncl 88> mlon = dimw(2)
ncl 89> delete(work)
ncl 90>
ncl 91> lon = OLR&lon
ncl 92> time = OLR&time
ncl 93> date = cd_calendar(time, -2) ; yyyymmdd
ncl 94>
ncl 95> ;************************************************
ncl 96> ; Apply the band pass filter to the original anomalies
ncl 97> ;************************************************
ncl 98> olr = wgt_runave_n_Wrap ( OLR, wgt, 0, 0) ; (time,lon)
ncl 99> u850 = wgt_runave_n_Wrap (U850, wgt, 0, 0)
ncl 100> u200 = wgt_runave_n_Wrap (U200, wgt, 0, 0)
ncl 101>
ncl 102> ;************************************************
ncl 103> ; remove temporal means of band pass series: *not* necessary
ncl 104> ;************************************************
ncl 105> olr = dim_rmvmean_n( olr, 0) ; (time,lon)
ncl 106> u850 = dim_rmvmean_n(u850, 0)
ncl 107> u200 = dim_rmvmean_n(u200, 0)
ncl 108>
ncl 109> ;************************************************
ncl 110> ; Compute the temporal variance at each lon
ncl 111> ;************************************************
ncl 112> var_olr = dim_variance_n_Wrap( olr, 0) ; (lon)
ncl 113> var_u850 = dim_variance_n_Wrap(u850, 0)
ncl 114> var_u200 = dim_variance_n_Wrap(u200, 0)
ncl 115>
ncl 116> ;************************************************
ncl 117> ; Compute the zonal mean of the temporal variance
ncl 118> ;************************************************
ncl 119> zavg_var_olr = dim_avg_n_Wrap( var_olr , 0)
ncl 120> zavg_var_u850 = dim_avg_n_Wrap( var_u850, 0)
ncl 121> zavg_var_u200 = dim_avg_n_Wrap( var_u200, 0)
ncl 122>
ncl 123> ;************************************************
ncl 124> ; Normalize by sqrt(avg_var*)
ncl 125> ;************************************************
ncl 126> olr = olr/sqrt(zavg_var_olr ) ; (time,lon)
ncl 127> u850 = u850/sqrt(zavg_var_u850)
ncl 128> u200 = u200/sqrt(zavg_var_u200)
ncl 129>
ncl 130> ;************************************************
ncl 131> ; Combine the normalized data into one variable
ncl 132> ;************************************************
ncl 133> cdata = new ( (/3*mlon,ntim/), typeof(olr),
getFillValue(olr))
ncl 134> do ml=0,mlon-1
ncl 135> cdata(ml ,:) = (/ olr(:,ml) /)
ncl 136> cdata(ml+ mlon,:) = (/ u850(:,ml) /)
ncl 137> cdata(ml+2*mlon,:) = (/ u200(:,ml) /)
ncl 138> end do
ncl 139>
=ncl 140> ;************************************************
ncl 141> ; Compute **combined** EOF; Sign of EOF is arbitrary
ncl 142> ;************************************************
ncl 143> eof_cdata = eofunc(cdata , neof, False) ; (neof,3*mlon)
ncl 144> print("==============")
ncl 145> printVarSummary(eof_cdata)
ncl 146> printMinMax(eof_cdata, True)
ncl 147>
ncl 148> eof_ts_cdata = eofunc_ts(cdata,eof_cdata,False) ; (neof,3*ntim)
ncl 149> print("==============")
ncl 150> printVarSummary(eof_ts_cdata)
ncl 151> printMinMax(eof_ts_cdata, True)
ncl 152>
incl 153> ;************************************************
ncl 154> ; For clarity, explicitly extract each variable. Create time
series
ncl 155> ;************************************************
ncl 156>
ncl 157> nvar = 3 ; "olr", "u850", "u200"
ncl 158> ceof = new( (/nvar,neof,mlon/), typeof(cdata),
getFillValue(cdata))
ncl 159>
ncl 160> do n=0,neof-1
ncl 161> ceof(0,n,:) = eof_cdata(n,0:mlon-1) ; olr
ncl 162> ceof(1,n,:) = eof_cdata(n,mlon:2*mlon-1) ; u850
ncl 163> ceof(2,n,:) = eof_cdata(n,2*mlon:) ; u200
ncl 164> end do
ncl 165>
(ncl 166> ceof!0 = "var"
ncl 167> ceof!1 = "eof"
ncl 168> ceof!2 = "lon"
ncl 169> ceof&lon = olr&lon
ncl 170>
ncl 171> ceof_ts = new( (/nvar,neof,ntim/), typeof(cdata),
getFillValue(cdata))
ncl 172> ceof_ts(0,:,:) = eofunc_ts_Wrap(
olr(lon|:,time|:),ceof(0,:,:),False) ; (0,neof,ntim)
ncl 173> ceof_ts(1,:,:) =
eofunc_ts_Wrap(u850(lon|:,time|:),ceof(1,:,:),False) ; (1,neof,ntim)
ncl 174> ceof_ts(2,:,:) =
eofunc_ts_Wrap(u200(lon|:,time|:),ceof(2,:,:),False) ; (2,neof,ntim)
ncl 175>
ncl 176> ;**********************************************t*
ncl 177> ; Add code contributed by Marcus N. Morgan, Florida Institute of
Technology; Feb 2015
ncl 178> ; Calculate % variance (pcv_ )accounted for by OLR, U850 and U200
ncl 179> ;************************************************
ncl 180>
ncl 181> pcv_eof_olr = new(neof,typeof(ceof))
ncl 182> pcv_eof_u850 = new(neof,typeof(ceof))
ncl 183> pcv_eof_u200 = new(neof,typeof(ceof))
ncl 184>
ncl 185> do n=0,neof-1
ncl 186> pcv_eof_olr(n) = avg((ceof(0,n,:)*sqrt(ceof at eval
(n)))^2)*100
ncl 187> pcv_eof_u850(n) = avg((ceof(1,n,:)*sqrt(ceof at eval
(n)))^2)*100
ncl 188> pcv_eof_u200(n) = avg((ceof(2,n,:)*sqrt(ceof at eval
(n)))^2)*100
ncl 189> ;;print("pcv: neof="+(n+1)+": "+pcv_eof_olr(n)+"
"+pcv_eof_u850(n)+" "+pcv_eof_u200(n))
ncl 190> end do
ncl 191>
ncl 192> ;************************************************
ncl 193> ; Change sign of EOFs for spatial structures of PC1 and PC2
ncl 194> ; to represent convection over the tropical Indian Ocean and the
tropical western Pacific Ocean, respectively
ncl 195> ; (Ad hoc approach)
ncl 196> ;************************************************
ncl 197>
ncl 198> imax_olr_eof1 = maxind(ceof(0,0,:))
ncl 199> imax_olr_eof2 = maxind(ceof(0,1,:))
ncl 200>
gncl 201> lonmax_eof1 = ceof&lon(imax_olr_eof1) ; longitude of max
value (i.e. suppressed convection)
ncl 202> lonmax_eof2 = ceof&lon(imax_olr_eof2)
ncl 203>
Pncl 204> if (lonmax_eof1.ge.60 .and. lonmax_eof1.lt.180) then ; Change
the sign of EOF1
ncl 205> ceof(:,0,:) = -ceof(:,0,:) ; if OLR
is positive
ncl 206> ceof_ts(:,0,:) = -ceof_ts(:,0,:) ; over
the tropical Indian Ocean
ncl 207> eof_cdata(0,:) = -eof_cdata(0,:)
ncl 208> eof_ts_cdata(0,:) = -eof_ts_cdata(0,:)
ncl 209> end if
ncl 210>
ncl 211> if (lonmax_eof2.ge.120 .and. lonmax_eof2.lt.180) then ; Change
the sign of EOF2
ncl 212> ceof(:,1,:) = -ceof(:,1,:) ; if OLR
is positive
ncl 213> ceof_ts(:,1,:) = -ceof_ts(:,1,:) ; over
the tropical western Pacific Ocean
ncl 214> eof_cdata(1,:) = -eof_cdata(1,:)
ncl 215> eof_ts_cdata(1,:) = -eof_ts_cdata(1,:)
ncl 216> end if
ncl 217>
ncl 218> print("==============")
= 0. ; reference line ncl 219>
printVarSummary(eof_cdata)
ncl 220> printMinMax(eof_cdata, True)
ncl 221>
=ncl 222> ;************************************************
ncl 223> ; Compute cross correlation of each variable's EOF time series at
zero-lag
ncl 224> ;************************************************
ncl 225> r_olr_u850 = escorc(ceof_ts(0,:,:) , ceof_ts(1,:,:) ) ; (neof)
ncl 226> r_olr_u200 = escorc(ceof_ts(0,:,:) , ceof_ts(2,:,:) )
ncl 227> r_u850_u200 = escorc(ceof_ts(1,:,:) , ceof_ts(2,:,:) )
ncl 228>
ncl 229> print("==============")
ncl 230> do n=0,neof-1
ncl 231> print("neof="+n \
oncl 231> +" r_olr_u850=" +sprintf("%4.3f",r_olr_u850(n)) \
ncl 231> +" r_olr_u200=" +sprintf("%4.3f",r_olr_u200(n)) \
; modify the panel plotncl 231> +"
r_u850_u200="+sprintf("%4.3f",r_u850_u200(n)) )
ncl 232> end do
ncl 233> print("==============")
ncl 234>
ncl 235> ;************************************************
ncl 236> ; Compute cross correlation of the multivariate EOF; EOF 1 vs EOF 2
ncl 237> ;************************************************
ncl 238>
ncl 239> mxlag = 25
ncl 240> rlag_01 = esccr(eof_ts_cdata(0,:),eof_ts_cdata(1,:), mxlag)
; (N,mxlag+1)
ncl 241> rlag_10 = esccr(eof_ts_cdata(1,:),eof_ts_cdata(0,:), mxlag)
; (N,mxlag+1)
ncl 242> ccr_12 = new ( (/2*mxlag+1/), float)
ncl 243>
(ncl 244> ccr_12(mxlag:) = rlag_10(0:mxlag)
ncl 245> ccr_12(0:mxlag) = rlag_01(::-1) ; reverse order
ncl 246> ;;print(ccr_12)
ncl 247>
ncl 248>
ncl 249> ;************************************************
ncl 250> ; Normalize the multivariate EOF 1&2 component time series
ncl 251> ; Compute (PC1^2+PC2^2): values > 1 indicate "strong" periods
ncl 252> ;************************************************
ncl 253> eof_ts_cdata(0,:) = eof_ts_cdata(0,:)/stddev(eof_ts_cdata(0,:))
ncl 254> eof_ts_cdata(1,:) = eof_ts_cdata(1,:)/stddev(eof_ts_cdata(1,:))
ncl 255>
ncl 256> mjo_ts_index = eof_ts_cdata(0,:)^2 + eof_ts_cdata(1,:)^2
ncl 257> mjo_ts_index_smt = runave(mjo_ts_index, 91, 0) ; 91-day running
mean
ncl 258>
rncl 259> nGood = num(.not.ismissing(mjo_ts_index)) ; # non-missing
ncl 260> nStrong = num(mjo_ts_index .ge. 1.0)
ncl 261> print("nGood="+nGood+" nStrong="+nStrong+"
nOther="+(nGood-nStrong))
ncl 262>
2ncl 263> ;************************************************
ncl 264> ; Write PC results to netCDF for use in another example.
ncl 265> ;************************************************
ncl 266> mjo_ts_index!0 = "time"
ncl 267> mjo_ts_index&time = time
ncl 268> mjo_ts_index at long_name = "MJO PC INDEX"
ncl 269> mjo_ts_index at info = "(PC1^2 + PC2^2)"
ncl 270>
*ncl 271> PC1 = eof_ts_cdata(0,:)
ncl 272> PC1!0 = "time"
*******ncl 273> PC1&time = time
ncl 274> PC1 at long_name = "PC1"
ncl 275> PC1 at info = "PC1/stddev(PC1)"
ncl 276>
ncl 277> PC2 = eof_ts_cdata(1,:)
ncl 278> PC2!0 = "time"
ncl 279> PC2&time = time
ncl 280> PC2 at long_name = "PC2"
ncl 281> PC2 at info = "PC2/stddev(PC2)"
ncl 282>
"ncl 283> diro = "/home/peni/"
ncl 284> filo = "MJO_PC_INDEX.nc"
ncl 285> system("/bin/rm -f "+diro+filo) ; remove any pre-existing file
ncl 286> ncdf = addfile(diro+filo,"c") ; open output netCDF file
ncl 287> ; make time an UNLIMITED
dimension
ncl 288> filedimdef(ncdf,"time",-1,True) ; recommended for most
applications
ncl 289> ; output variables directly
ncl 290> ncdf->MJO_INDEX = mjo_ts_index
ncl 291> ncdf->PC1 = PC1
ncl 292> ncdf->PC2 = PC2
ncl 293>
ncl 294> ;------------------------------------------------------------
ncl 295> ; PLOTS
ncl 296> ;------------------------------------------------------------
Pncl 297>
Cncl 298> yyyymmdd = cd_calendar(time, -2)
ncl 299> yrfrac = yyyymmdd_to_yyyyfrac(yyyymmdd, 0.0)
ncl 300> delete([/ yrfrac at long_name, lon at long_name /])
ncl 301>
ncl 302> day = ispan(-mxlag, mxlag, 1)
ncl 303> ;day at long_name = "lag (day)"
ncl 304>
ncl 305> pltPath = pltDir+pltName
ncl 306>
ncl 307> wks = gsn_open_wks(pltType,pltPath)
ncl 308> gsn_define_colormap(wks,"default")
ncl 309> plot = new(3,graphic)
ncl 310>
ncl 311> ;************************************************
ncl 312> ; Multivariate EOF plots
ncl 313> ;************************************************
ncl 314> rts = True
ncl 315> rts at gsnDraw = False ; don't draw yet
ncl 316> rts at gsnFrame = False ; don't advance frame yet
ncl 317> rts at gsnScale = True ; force text scaling
ncl 318>
ncl 319> rts at vpHeightF = 0.40 ; Changes the aspect ratio
ncl 320> rts at vpWidthF = 0.85
ncl 321> rts at vpXF = 0.10 ; change start locations
ncl 322> rts at vpYF = 0.75 ; the plot
ncl 323> rts at xyLineThicknesses = (/2, 2, 2/)
ncl 324> rts at xyLineColors = (/"black","red","blue"/)
ncl 325> rts at gsnYRefLine = 0. ; reference line
ncl 326> rts at trXMaxF = max(lon)
ncl 327> rts at trXMinF = min(lon)
ncl 328>
ncl 329> rts at pmLegendDisplayMode = "Always" ; turn on legend
ncl 330> rts at pmLegendSide = "Top" ; Change
location of
ncl 331> rts at pmLegendParallelPosF = 1.16 ; move units
right
ncl 332> rts at pmLegendOrthogonalPosF = -0.50 ; move units
down
ncl 333> rts at pmLegendWidthF = 0.15 ; Change width
and
ncl 334> rts at pmLegendHeightF = 0.15 ; height of
legend.
ncl 335> rts at lgLabelFontHeightF = 0.0175
ncl 336>
ncl 337>
ncl 338> rtsP = True ; modify the
panel plot
ncl 339> ; rtsP at gsnMaximize = True ; large format
ncl 340> rtsP at gsnPanelMainString = "Multivariate EOF: 15S-15N:
"+yrStrt+"-"+yrLast
ncl 341>
ncl 342> do n=0,neof-1
ncl 343> rts at xyExplicitLegendLabels = (/"OLR: "+sprintf("%4.1f",
pcv_eof_u200(n)) +"%" \
ncl 343> ,"U850: "+sprintf("%4.1f",
pcv_eof_u850(n))+"%" \
ncl 343> ,"U200: "+sprintf("%4.1f",
pcv_eof_olr(n))+"%" /)
ncl 344> rts at gsnLeftString = "EOF "+(n+1)
ncl 345> rts at gsnRightString = sprintf("%3.1f",ceof at pcvar(n)) +"%"
ncl 346> plot(n) = gsn_csm_xy (wks,lon,ceof(:,n,:),rts)
ncl 347> end do
ncl 348> gsn_panel(wks,plot(0:1),(/2,1/),rtsP) ; now draw as one plot
ncl 349>
ncl 350> ;-----------------------------------------
ncl 351> ; The following doesn't work with some older versions of NCL
ncl 352> ; With old versions, the user must delete each individually.
ncl 353> ;-----------------------------------------
ncl 354> delete([/ rts at xyExplicitLegendLabels, rts at pmLegendDisplayMode \
ncl 354> , rts at xyLineThicknesses , rts at gsnLeftString \
ncl 354> , rts at gsnRightString , rts at xyLineColors \
ncl 354> , rts at trXMaxF , rts at trXMinF /]
)
ncl 355>
ncl 356> lag = ispan(-mxlag,mxlag,1)
ncl 357> lag at long_name = "lag (days)"
ncl 358>
ncl 359> plot(0) = gsn_csm_xy (wks, lag ,ccr_12,rts)
ncl 360> rtsP at gsnPanelMainString = "Cross Correlation: Multivariate
EOF: 15S-15N: " \
ncl 360> + yrStrt+"-"+yrLast
ncl 361> rtsP at gsnPaperOrientation = "portrait" ; force portrait
ncl 362> gsn_panel(wks,plot(0),(/1,1/),rtsP) ; now draw as one plot
ncl 363>
ncl 364> ;************************************************
ncl 365> ; MJO "strong" index
ncl 366> ;************************************************
ncl 367> rts at gsnYRefLine = 1.0
ncl 368> rts at gsnYRefLineColor = "black"
ncl 369> rts at xyMonoDashPattern = True
ncl 370> rts at xyLineColors = (/"black", "blue"/)
ncl 371> rts at xyLineThicknesses = (/1, 2/)
ncl 372> rts at pmLegendDisplayMode = "Always" ; turn on legend
ncl 373> rts at pmLegendWidthF = 0.12 ; Change width
and
ncl 374> rts at pmLegendHeightF = 0.10 ; height of
legend.
ncl 375> rts at pmLegendParallelPosF = 0.86 ; move units
right
ncl 376> rts at pmLegendOrthogonalPosF = -0.40 ; move units
down
ncl 377> rts at xyExplicitLegendLabels = (/"daily", "91-day runavg" /)
ncl 378>
ncl 379> mjo_ind_plt = new ( (/2,ntim/), typeof(mjo_ts_index))
ncl 380> mjo_ind_plt(0,:) = mjo_ts_index
ncl 381> mjo_ind_plt(1,:) = (/ mjo_ts_index_smt /)
ncl 382> plot(0) = gsn_csm_xy(wks, yrfrac,mjo_ind_plt,rts)
ncl 383>
ncl 384> rtsP at gsnPanelMainString = "MJO Index: (PC1^2+ PC2^2) :
15S-15N: "+yrStrt+"-"+yrLast
ncl 385> gsn_panel(wks,plot(0),(/1,1/),rtsP) ; now draw as one plot
ncl 386>
ncl 387> end
sub drveof: ier,jopt= 10 0 returned from vcmssm/crmssm
warning:eofunc: 10 eigenvectors failed to converge
(0) ==============
Variable: eof_cdata
Type: double
Total Size: 6960 bytes
870 values
Number of Dimensions: 2
Dimensions and sizes: [2] x [435]
Coordinates:
Number Of Attributes: 5
_FillValue : -999000000
method : no transpose
matrix : covariance
pcvar : ( -9.99e+08, -9.99e+08 )
eval : ( -999000000, -999000000 )
(0)
(0) min=-999000000 max=-999000000
warning:eofunc_ts: 2 eigenvectors failed to converge
(0) ==============
Variable: eof_ts_cdata
Type: double
Total Size: 5840 bytes
730 values
Number of Dimensions: 2
Dimensions and sizes: [2] x [365]
Coordinates:
Number Of Attributes: 3
_FillValue : -999000000
matrix : covariance
ts_mean : ( -999000000, -999000000 )
(0)
(0) min=-999000000 max=-999000000
warning:eofunc_ts: 2 eigenvectors failed to converge
warning:eofunc_ts: 2 eigenvectors failed to converge
warning:eofunc_ts: 2 eigenvectors failed to converge
warning:avg: Entire input array contains missing values, can't compute
average
warning:avg: Entire input array contains missing values, can't compute
average
warning:avg: Entire input array contains missing values, can't compute
average
warning:avg: Entire input array contains missing values, can't compute
average
warning:avg: Entire input array contains missing values, can't compute
average
warning:avg: Entire input array contains missing values, can't compute
average
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading unknown
fatal:["Execute.c":8635]:Execute: Error occurred at or near line 201
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191216/242a08f0/attachment.html>
More information about the ncl-talk
mailing list