PRO precip30

;This program loads annual precipitation data from the Minocqua
;NCDC site and plots 30 year averages

loadcomma,'/eddy/s2/bcook/wdepth/woodruff_annprecip',precip

;Define plot colors
;   0=black
;   1=red
;   2=green
;   3=blue
;   4=cyan
;   5=magenta
;   6=yellow
;   7=white
red=[0,1,0,0,0,1,1,1]*255
green=[0,0,1,0,1,0,1,1]*255
blue=[0,0,0,1,1,1,0,1]*255
tvlct,red,green,blue

set_plot,'ps'
device,file='/eddy/s2/bcook/wdepth/precip30.ps',/inches,xoffset=0.5,yoffset=5,xsize=7,ysize=5,/helvetica,/color

plot,precip(0,*),ts_smooth(reform(precip(1,*)),30,/backward),yrange=[0,120],/nodata,xtitle='Year',ytitle='Annual precipitation (cm)',charsize=1.5,charthick=2

FOR i=0,n_elements(precip(1,*))-1 DO BEGIN
   oplot,[precip(0,i),precip(0,i)],[0,precip(1,i)],thick=2,color=0
endfor

oplot,precip(0,*),ts_smooth(reform(precip(1,*)),30,/backward),color=1,thick=4

device,/close
set_plot,'x'

return
end