; Procedure to compare SODAR and MASS-DIMM integrals
; Jan 20, 2005
  ; Ignore Nov. 24 
;---------------------------------------------------------
pro mdsod


dds = [25, 26, 27, 28, 30, 1]
mms = [11, 11, 11, 11, 11, 12]
ndays = 6

sodall = fltarr(ndays,40)
mdall = fltarr(ndays,40)

 restore, 'massdimm.idl'

for id = 0,ndays-1 do begin
  print, 'Date: ', dds[id],'/', mms[id]

;--------MASS-DIMM data --------------------------------
; mddat structure:
; prof = { UTyear: 0, UTmonth: 0, UTday: 0, UTmass:0.0, $
;	  fsee:0.0, Intcn2:FLTARR(6), UTdimm:0.0, see:0.0 }

dd1 = dds[id]+1   & mm1 = mms[id]
if (dd1 gt 30) then begin dd1=1 & mm1=mm1+1 & endif
morn = where((mddat.UTmonth eq mm1) and (mddat.UTday eq dd1) and (mddat.UTmass lt 12) )
nmd = n_elements(morn)
md = mddat(morn)

utmd = fltarr(nmd) & glint = fltarr(nmd)
for i=0,nmd-1 do begin 
  utmd[i] = md[i].UTmass
  glint[i] = (md[i].see^(5./3.)*6.8 - total(md[i].intcn2) )  > 0.
endfor

; ------ SODAR data
sod = selsod(mms[id], dds[id])
utsod = sod.uthour -1.
altsod = sod[0].alt

; SODAR seeing
nsod = n_elements(sod)   
sodsee = fltarr(nsod)
wt =  (1. - (altsod-250.)/250. ) > 0.
wt(where (wt gt 1)) = 1
for i=0,nsod-1 do sodsee[i] = (total(sod[i].cn2*wt)/6.8e-13)^0.6

; Temporal binning of MASS-DIMM
for i=0,nsod-1 do begin
  sel = where( (utmd ge utsod[i] - 0.33333) and   (utmd le utsod[i]) )
  nsel = n_elements(sel)
  if (nsel ge 4 ) then begin
    sodall[id,i] = sodsee[i]
    mdall[id,i] = total(glint(sel))/nsel
  endif
endfor


  plot, sodall[id,*]
  oplot, (mdall[id,*]/6.8)^0.6, linestyle=2

  stop


endfor
print, 'All nights are processed!'

sel = where(sodall gt 0.)

  sodsee = sodall(sel)
  mdsee = (mdall(sel)/6.8)^0.6

!p.charsize = 1.5
  plot, sodsee, mdsee, psym=1, xtitle='SODAR, arcsec', ytitle='MASS-DIMM, arcsec', xr=[0,1.5],yr=[0,1]
  oplot, [0,2], [0,2]


stop



end