PRO HOV850vMERRA

;This program is a test to produce a hovmoller plot in IDL
;Here I ingest MERRA 850v and perform a series of data organization and then contour the data
;Created by Teddy Allen , Sept. 21, 2009

;longitude 76-80 has weird values! CHECK

RESTORE,'V850_coarse.79-05.sav'

set_plot,'ps'
device,/color, filename='c:\testdata\hovm850vMERRA.ps'
loadct,13


file_a='C:\Documents and Settings\Teddy\My Documents\TRMMmerra\2002TRMMmerra_9080W.txt' ;read in 365 daily data 
data=read_ascii(file_a)
f1=data.(0)

file_b='C:\Documents and Settings\Teddy\My Documents\TRMMmerra\2002TRMMmerra_8070W.txt' 
data=read_ascii(file_b)
f2=data.(0)

file_c='C:\Documents and Settings\Teddy\My Documents\TRMMmerra\2002TRMMmerra_7060W.txt' 
data=read_ascii(file_c)
f3=data.(0)

month = ['JAN','FEB','MAR','APR','MAY','JUN','JUL','AUG','SEP','OCT','NOV','DEC']
x=findgen(365)
avgLATv850=total(V850[*,20:22,*],2)/3.    ;averages within 12.75N-22.75N
avgLATv850=ROTATE(avgLATv850,4)

time=findgen(365)
stepmax= max(avglatv850(81:90,8037:8401))
stepmin= min(avglatv850(81:90,8037:8401))
skip=(stepmax-stepmin)/10
step1=skip+stepmin
step2=step1+skip
step3=step2+skip
step4=step3+skip
step5=step4+skip
step6=step5+skip
step7=step6+skip
step8=step7+skip
step9=step8+skip
step10=step9+skip
lev=[stepmin,step1,step2,step3,step4,step5,step6,step7,step8,step9,stepmax] ;used in CONTOUR

levs=findgen(10)                                                            ;used in MAKE_KEY
month=['JAN','FEB','MAR','APR','MAY','JUN','JUL','AUG','SEP','OCT','NOV','DEC',' ']

!P.MULTI = [0, 1, 2, 0, 1]
CONTOUR, avgLATv850(81:90,8037:8401),lon(81:90),time,color=0,$
thick=3,/follow,nlevels=10,levels=lev,/fill,xrange=[270,300],$
ystyle=1,yrange=[0,364],title='2002 MERRA v850 (averaged from 12.75N - 22.75N)',$
xtitle='longitude',charsize=1.0,charthick=2.0,yticks=12,ytickname=month,$
xticks=6

MAKE_KEY, COLORS=findgen(n_elements(levs)) * !D.N_colors/n_elements(levs),$
labels = round_any([(stepmin),(step1),(step2),(step3),(step4),(step5),(step6),(step7),(step8),(step9),(stepmax)],sig=3),$
units = 'm/s' ,/orientation, charsize = .6, charthick=1.5, /noborder,title='  m/s'

PLOT,x,f1,ystyle=8,xstyle=1,xrange=[0,365],title='2002 TRMM daily rainfall (averaged from 12.75N - 22.75N)',charthick=2.0,xticks=12,xtickname=month
OPLOT,x,f1,color=255
OPLOT,x,f2,color=150,thick=2
OPLOT,x,f3,color=50,thick=2
AXIS,yaxis=1,xstyle=1,xrange=[0,365],ytitle='accumulated rainfall (mm)',charthick=2.0,/save

xyouts,25,27,'270-280 lon',color=255, charthick=2
xyouts,25,24,'280-290 lon',color=150, charthick=2
xyouts,25,21,'290-300 lon',color=50, charthick=2

device,/close
stop
END