;program to compare and 'integrated' depol ratio for the whole beam
; and the microwave radiometer total liquid water path.
;
;IN workspace should be the variable "file" which contains the
; filenames to all the lidar files.

tintdep=0.
twater=0.
ttime=0.
for ii=0,n_elements(file)-1 do begin

  restore,file[ii]
  time=time10
  dep=dep10
  water=water10
  timeax=time10   ;time05ax

  intdep=fltarr(n_elements(time))*0

  for i=0,n_elements(time)-1 do begin
    iwhere=where((dep[i,*] gt 0) and (dep[i,*] le .12),numlq)

    if numlq ne 0 then intdep[i]=total(dep[i,iwhere])*numlq*0.030
  endfor

  tintdep=[tintdep,intdep]
  twater=[twater,water]
  ttime=[ttime,time+ii*24]

  !p.multi=[0,1,2]
  title='MWR, Integrated Liquid Water Path, '+strmid(file[ii],0,6)
  plot,timeax,water,xstyle=1,ystyle=1,$
    title=title,$
    ytitle='LWP (cm)',min_value=0,color=1
  plot,time,intdep,xstyle=1,ystyle=1,$
    title='Integrated Liquid Depol Ratios',$
    xtitle='Time (hrs)',ytitle='Int. Dep. Ratio',color=1
  !p.multi=0

  wait,1
endfor
end