; -IDL script to create .GIF files of AMSU-B Channel 16 (89GHz) 
;- Created: 29 Jun 99 Kurt F. Brueske
;- Modified: 11 Jun 00 Open tbdata.dat, lat.dat, lon.dat, parameters.dat and nhcbulletin.txt  generated by
;-                    ty:/home2/kurtb/amsu/run/epac/limbfixer_epac'
;-           11 Jun 00 Copied from ch16_epac41.pro and modified for use with sectorized AMSU-B
;-                    data.  Numscans no longer provided in parameters.dat since value is always 60... 
;-	      5 Dec 00 Copied from *epac*.pro

;- WESTERN NORTH PACIFIC BASIN PROCESSING

;- Set up Z buffer in 2D mode (allows script to run in background)...

set_plot,'Z'
device,z_buffer=0,set_colors=256,set_resolution=[725,550],set_character_size = [ 6, 9 ]

;- Get storm information...

info=strarr(6)
openr,1,'/home2/kurtb/amsu/run/derrick/science2002/procfiles/parameters.dat'
readf,1,info
close,1

openr,1,'/home2/kurtb/amsu/run/derrick/science2002/procfiles/parameters.dat'
readf,1,latin,lonin
close,1

;- Define arrays and number of scan lines in sectorized AMSU-B file...

numscans=60
t=fltarr(5,90,numscans)
lat=fltarr(90,numscans)
lon=lat

;- Open data files...
openr,1,'/home2/kurtb/amsu/run/derrick/science2002/procfiles/tbdatab.dat'
readf,1,t
close,1
openr,2,'/home2/kurtb/amsu/run/derrick/science2002/procfiles/latb.dat'
readf,2,lat
close,2
openr,3,'/home2/kurtb/amsu/run/derrick/science2002/procfiles/lonb.dat'
readf,3,lon
close,3

;- Create plots...

sat1 = info[5]
sat = String(sat1, Format='(A7)')

erase
;-Channel 16 
plotamsub,reform(t[0,*,*]),numscans,lat,lon,latin,lonin,info[3],info[2],8,0,0,sat
restore,filename='/home2/kurtb/amsu/idl/amsubcolors.sav'
tvlct,r,g,b
write_gif,'/home2/kurtb/amsu/run/derrick/science2002/procfiles/amsub_ch16.gif',tvrd(),r,g,b


exit