[ncl-talk] Sub: Set lonL/lonR error
Adv
advita6 at gmail.com
Tue May 17 13:09:34 MDT 2016
Hi,
I am getting this error. Could someone help me to fix this bug?
Thanks
fatal:["NclFile.c":2103]:Subscript out of range, error in subscript #2
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 51 in file
regress_dennis.ncl
load "./contributed.ncl_beta_640"
begin
;************************************************
; Specify geographical region and time span (year-month start and end
;************************************************
latS = 40
; latS = 0
latN = 49
lonL = -116
lonR = -90
ymStrt = 195101
ymLast = 201012
pltType = "x11"
pltName = "regress"
pltTitle = "Globe: "+(ymStrt/100)+"-"+(ymLast/100)
;************************************************
; Read from netCDF file: variable is type short...unpack
;************************************************
version= "v2" ; "v2", "v2c"
diri = "./"
; fili = "air.sig995.mon.mean."+version+".nc"
fili = "air.sig995.mon.mean.nc"
f = addfile(diri+fili,"r")
YYYYMM = cd_calendar( f->time, -1)
yyy = yyyymm/100
dimx = dimsizes(x)
ntim = dimx(0) ; all years and months
nlat = dimx(1)
mlon = dimx(2)
year = ispan(yyyy(0), yyyy(ntim-1), 1)
nyrs = dimsizes(year)
;************************************************
; Areal averages: cos(lat) is good enough
;************************************************
wgt = cos(0.01745329*x&lat)
xann = month_to_annual(x , 1) ; [year| 60]x[lat|
91]x[lon| 180]
xavg = wgt_areaave_Wrap(xann , wgt, 1.0, 1) ; [year| 110]
xann&year = year
xann at long_name = "Annual Air Temperature ("+version+": sigma=0.995)"
printVarSummary(xann)
printMinMax(xann,0)
;return
;************************************************
; Calculate the regression coefficients (slopes)
;************************************************
rc = regCoef(year,xann(lat|:,lon|:,year|:))
rc at long_name = "Trend"
rc at units = "degC/"+nyrs+"_year"
printMinMax(rc,0)
;************************************************
; plotting parameters
;************************************************
wks = gsn_open_wks("x11","regress_4_"+version) ; specifies a ps
plot
res = True
res at gsnMaximize = True ; make large
res at cnFillOn = True ; turn on color
res at cnLinesOn = False ; turn off contour lines
res at cnLineLabelsOn = False ; turn off contour line
labels
;;res at cnFillMode = "RasterFill"
res at cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
; res at cnMinLevelValF = -2.0 ; set min contour level
res at cnMinLevelValF = -3.0 ; set min contour level
; res at cnMaxLevelValF = 2.0 ; set max contour level
res at cnMaxLevelValF = 5.0 ; set max contour level
; res at cnLevelSpacingF = 0.2 ; set contour interval
res at cnLevelSpacingF = 0.5 ; set contour interval
res at mpFillOn = False ; turn off default
background gray
;res at mpCenterLonF = 180
res at gsnCenterString = year(0)+"-"+year(nyrs-1)
res at tiMainString = "20th Century Reanalysis ("+version+"):
sig995" ; fili
plot = gsn_csm_contour_map_ce(wks,rc,res)
;************************************************
; Perform linear regression on selected annual mean
;************************************************
rcTest = regline_stats(year, xavg) ; degC/year
; rcTest at long_name = "trend: ("+toint(LAT)+","+toint(LON)+")"
rcTest at long_name = "trend:"
rcTest at units = "degC/year"
rcTest = rcTest*nyrs ; (C/year)*(nyrs)
rcTest at units = "degC/"+nyrs+"_year"
pltarry = new ( (/2,nyrs/), typeof(xann))
pltarry(0,:) = xavg ; use markers
pltarry(1,:) = rcTest at Yest ; regression values
resxy = True ; plot mods desired
resxy at xyMarkLineModes = (/"Markers","Lines"/); choose which have
markers
resxy at xyMarkerColors = (/"Red", "Blue" /); choose which have
markers
resxy at xyMarkers = 16 ; choose type of
marker
resxy at xyMarkerSizeF = 0.0075 ; Marker size
(default 0.01)
resxy at xyDashPatterns = 0 ; solid line
resxy at xyLineThicknesses = (/2,3/)
resxy at xyLineColors = (/ "blue" , "black" /)
resxy at tmYLFormat = "f" ; not necessary but
nicer labels
;;resxy at trXMinF = min(year)
resxy at trXMaxF = max(year)
resxy at tiMainString = "regline: 20th Renalysis ("+version+"):
trend="+sprintf("%5.2f", rcTest)
;---Make legend smaller and move into plot
resxy at pmLegendDisplayMode = "Always" ; turn on legend
resxy at pmLegendSide = "Top" ; Change location of
resxy at pmLegendParallelPosF = 0.225 ; move units Leftt
resxy at pmLegendOrthogonalPosF = -0.30 ; move units down
resxy at pmLegendWidthF = 0.10 ; Change width and
resxy at pmLegendHeightF = 0.100 ; height of legend.
resxy at lgPerimOn = True ; turn off/on box
around
resxy at lgLabelFontHeightF = .015 ; label font height
resxy at xyExplicitLegendLabels = (/"data" , "regline" /)
plot = gsn_csm_xy
(wks,year,pltarry,resxy) ; create plot
end
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160517/29d0342d/attachment.html
More information about the ncl-talk
mailing list