|
In my code below, IDL is referring the COLORBAR function as undefined and I do not understand why the function is not being called when its proper syntax. Could someone point me in the right direction please? Thanks for the help!
trmm_used *=3.
for ilat=0, num_lat-1 do begin
for ilon=0, num_lon-1 do begin
values = trmm_used[ilon,ilat,*,0,*,*]
values_list= where(values GE 0, amt)
if amt ge 0 then begin
rainfall_mean = mean(values[values_list])
monthly_rainfall[ilon,ilat] = rainfall_mean
endif
endfor
endfor
WINDOW, 3, XSIZE=900, YSIZE=900
LOADCT, 39, NCOLORS=20
DEVICE, DECOMPOSED=0
!P.COLOR=2
!P.BACKGROUND=255
ERASE, 255
MAP_SET, 0,0, /CYLINDRICAL, LIMIT=[LatMin, LonMin, LatMax, LonMax], $
POSITION=[.1,.1,.95,.95], /ISOTROPIC, /NOBORDER, /NOERASE, LONLAB=Longitude, $
LATLAB=Latitude, TITLE='January Rainfall'
CONTOUR, LEVELS=findgen(20)*0.1, C_COLORS=indgen(20), monthly_rainfall, $
longitude[indices_lon], latitude[indices_lat], $
/OVERPLOT, /CELL_FILL, POSITION=[0.15,0.15,0.95,0.80]
XCHAR=1
YCHAR=1
COLORBAR, RGB_TABLE=39, POSITION=[0.15,0.85,0.95,0.90], $
RANGE=[Min(monthly_rainfall), MAX(monthly_rainfall)], FORMAT='(F0.2)'
MAP_CONTINENTS, /COASTS, /COUNTRIES, COLOR=0, THICK=1, /HIRES
AXIS, XAXIS=0, XRANGE=[LonMin, LonMax], XSTYLE=1, COLOR=0, CHARSIZE=XCHAR, XTHICK=2, CHARTHICK=2, $
XTICKINTERVAL=1, XMINOR=4
AXIS, YAXIS=0, YRANGE=[LonMin, LonMax], YSTYLE=1, COLOR=0, CHARSIZE=YCHAR, YTHICK=2, CHARTHICK=2, $
YTICKINTERVAL=1, YMINOR=4
AXIS, XAXIS=1, XRANGE=[LonMin, LonMax], XSTYLE=1, COLOR=0, CHARSIZE=0.01, XTHICK=2, CHARTHICK=1, $
XTICKINTERVAL=1, XMINOR=4
AXIS, YAXIS=1, YRANGE=[LonMin, LonMax], YSTYLE=1, COLOR=0, CHARSIZE=0.01, YTHICK=2, CHARTHICK=1, $
YTICKINTERVAL=1, YMINOR=4
END
|