pro read_inundation

path='/model/john/data/FW_FRAC_LAT_LON_GRID_0_25/'
files=file_search(path,'*.DAT',count=nfiles)

;files are 0.25 x 0.25 deg. binary (little endian)
dx=0.25
tmparr=fltarr(360./dx,180./dx)

;ps_start
region='globe'
for i=0,nfiles-1 do begin
;for i=0,10 do begin
  openr,unit,files[i],/get_lun
  readu,unit,tmparr
  free_lun,unit
  
  ;filter out the following
  ;1. ocean, large inland water bodies, cities etc -9999
  ;2. snow -999
  ;3. retrievals during rainfall -99
  ;also screen out zero values
  gd=where(tmparr gt 0, complement=bd, countgd)
  tmparr[bd]=!values.f_nan
  tmparr=rotate(tmparr,7)
  
  ;title/date
  tmpstr=file_basename(files[i],'.DAT')
  yrmo=strmid(tmpstr,0,7)
  dy1=strmid(tmpstr,1,2,/reverse_offset)
  dy0=string(fix(dy1)-10,format='(I02)')
  datestr=yrmo+'_'+dy0+'-'+dy1
  
  ps_start,filename=path+tmpstr+'_'+region+'.ps'
  plot_tm5,globgrid=tmparr,dx=dx,coldiv=findgen(11)/20.,region=region,cbartitle='Fractional Inundation',plottitle=datestr,/noplot
  ps_end,/gif

endfor

ps_end

end