Global (60N-60S) half hourly, 0.5 degree lat/lon, geostaionary satellite IR products (in hourly files) file description (data collection begins Feb 26, 1999) Hourly files are unix compressed and time tagged so that: globl_YYYYMMDDHH_0.5-deg_i2.Z (ie globl_1999030605_0.5-deg_i2.Z) YYYY=year, MM=month, DD=day, HH=hour (contains first :00 data, then :30 data) Note!! A satellite zenith angle correction (dependent upon 1. pixel temperature 2. satellite zenith angle 3. latitude 4. season) has been applied to IR pixels sza > 26 deg before product processing. Each hourly IR product file contains global (60N-60S) 0.5 deg lat/lon geostationary satellite IR products for both half hours (:00 and :30) in integer*2 big-endian binary (ie SGI, HP, IBM) for; 1. odd # satellites only [1=GMS, 3=GOES-8, 5=MET-5 (5W-to-70E only)] 2. even # satellites only [2=GOES-10, 4=MET-7, 5=MET-5 (70E-133E only)] 3. all satellites merged [data closest to geos' subsat used in overlap] Product arrays: (MISSING VAULES = -9999) SCALING FACTOR || 1. GPI fraction (ratio of pixels colder than 235K)/(all pixels) *10000 2. Standard deviation of pixels about 0.5 deg IR mean temperature *10 3. 0.5 deg IR mean temperature *10 4. Maximum IR pixel temp for 0.5 deg region *10 5. Maximum IR pixel temp for 0.5 deg region *10 6. Satellite ID for data (1=GMS, 2=GOES-10, 3=GOES-8, 4=MET-7, 5=MET-5) Each 0.5 degree lat/lon product array contains 720x240 cell-centered values from 60N to 60S, 0E to 360E such that the grid is filled: Lat Lon 0.0E 0.5E 1.0E . . . . . . 359.0E 359.5E 360.0E 60.0N (1,1) (2,1) . . . . . . (719,1) (720,1) 59.5N (1,2) (2,2) . . . . . . (719,2) (720,2) 59.0N . . . . . . . . . . . 59.0S (1,239) (2,239) . . . . . . (719,239) (720,239) 59.5S (1,240) (2,240) . . . . . . (719,240) (720,240) 60.0S The Product files can be read in FORTRAN with this code: c NOTE!! If SGI compile with option -old_rl (recl=bytes, not 4 byte words) c P5MAX = the 6 product arrays (above) c P5LEV = [1=odd # sats, 2=even # sats, 3=merged sats (above)] Parameter (I5MAX=720, J5MAX=240, P5MAX=6, P5LEV=3) integer*2 iprod(I5MAX,J5MAX,P5MAX,P5LEV) character*80 mergein do 15 iyrmndy = 19990323, 19990323 do 20 ihr = 0, 23 do 25 ihlfhr = 0, 30, 30 mergein = "/export/lnx57/wd52rj/mcidas/data/glbl_ir/mrg/globl &_YYYY0M0D0H_0.5-deg_i2" write(mergein(52:59),'(i8)') iyrmndy if (ihr.lt.10) then write(mergein(61:61),'(i1)') ihr else write(mergein(60:61),'(i2)') ihr endif write(*,*) mergein if (ihlfhr.eq.0) then open (unit=10, file = mergein, access = 'direct', 1 form='unformatted', recl=720*240*6*3*2) irrec = 1 else irrec = 2 endif write(*,*) irrec read (10, rec=irrec) iprod do 900 ilv=1,P5LEV write (6,'(//,'' lev '',/,i8)') ilv do 900 ipd=1,P5MAX write (6,'(/,'' prod '',/,i8)') ipd write (6,'(12(1x,i9))') ((iprod(i5,j5,ipd,ilv),i5=35,46), & j5=112,128) 900 continue 25 continue 20 continue 15 continue stop end If you use GrADS graphics, the products can be dispalyed with this control file. NOTE!! big-endian option (SGI binary) and template option DSET /export/lnx57/wd52rj/mcidas/data/glbl_ir/mrg/globl_%y4%m2%d2%h2_0.5-deg_i2 OPTIONS yrev big_endian template UNDEF -9999 TITLE IR product files XDEF 720 LINEAR 0.25 0.5 YDEF 240 LINEAR -59.75 0.5 ZDEF 01 LEVELS 1 TDEF 9999 LINEAR 0z26feb1999 30mn VARS 18 gpio 1 -1,40,2,-1 gpi fraction (*10000) odd sats stdo 1 -1,40,2,-1 standard deviation (*10) odd sats tmpo 1 -1,40,2,-1 mean temp (*10) odd sats maxo 1 -1,40,2,-1 max pixel temp (*10) odd sats mino 1 -1,40,2,-1 min pixel temp (*10) odd sats sido 1 -1,40,2,-1 satellite ID odd sats gpie 1 -1,40,2,-1 gpi fraction (*10000) even sats stde 1 -1,40,2,-1 standard deviation (*10) even sats tmpe 1 -1,40,2,-1 mean temp (*10) even sats maxe 1 -1,40,2,-1 max pixel temp (*10) even sats mine 1 -1,40,2,-1 min pixel temp (*10) even sats side 1 -1,40,2,-1 satellite ID (*10) even sats gpim 1 -1,40,2,-1 gpi fraction (*10000) merged sats stdm 1 -1,40,2,-1 standard deviation (*10) merged sats tmpm 1 -1,40,2,-1 mean temp (*10) merged sats maxm 1 -1,40,2,-1 max pixel temp (*10) merged sats minm 1 -1,40,2,-1 min pixel temp (*10) merged sats sidm 1 -1,40,2,-1 satellite ID merged sats ENDVARS