; -IDL script to create .GIF files of AMSU-A Channel 5 
;- Created: 29 Jun 99 Kurt F. Brueske
;- Modified: 5 Aug 99 Open tbdata.dat, lat.dat, lon.dat, parameters.dat and nhcbulletin.txt  generated by
;-                    ty:/home2/kurtb/amsu/run/epac/limbfixer_epac'
;-          24 May 00 Copied from gifscript_ch5_epac41.pro and modified for use with sectorized AMSU-A
;-                    data.  Numscans no longer provided in parameters.dat since value is always 22... 
;-           7 Jun 00 Found through experimentation that by adding 360. to 'lon' we can universally use
;-                    'plotamsu_ch5/6/7/8.pro'...
;-	     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-A file...

numscans=22
t=fltarr(15,30,numscans)
lat=fltarr(30,numscans)
lon=lat

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

;- Create plots...

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

erase
;-Channel 5
plotamsu_ch5,reform(t[4,*,*]),numscans,lat,lon,latin,lonin,info[3],info[2],8,1,1,sat
loadct,39
tvlct,r,g,b,/get
write_gif,'/home2/kurtb/amsu/run/derrick/science2002/procfiles/amsua_ch5.gif',tvrd(),r,g,b


exit