; Procedure to plot SODAR results
; Input: month and date of the evening
pro plotprof, mm, dd
restore, 'sodar1.idl'

out = selsod(mm,dd)
; prof = { UTyear: 0, UTmonth: 0, UTday: 0, UThour:0.0, UTmin:0.0, $
;	    nlevels:40, alt:FLTARR(40), cn2:FLTARR(40) }
nprof = n_elements(out)

;  nalt = n_elements(out[0].cn2)
nalt = 40

ut = out.uthour -1.

print, 'Total profiles found: ', nprof
print, 'UT range: ',ut[0], ut[nprof-1]


tmp = fltarr(nprof,nalt) 
for i=0,nprof-1 do tmp(i,1:nalt-1) = out[i].cn2

nut = 450
utmin = 0.
utmax = 9.00
utstep = (utmax-utmin)/nut
utgrid = findgen(nut)*utstep + utmin
ntick = floor(utmax-utmin)
jtick = fix( (findgen(ntick))/utstep)

index = fltarr(nut)
for j=0,nut-1 do begin
  dt = abs(utgrid[j]-ut[0])
  index[j]=0
  for i=1,nprof-1 do if (abs(utgrid[j]-ut[i]) lt dt) then begin 
       dt = abs(utgrid[j]-ut[i])
       index[j] = i
     endif  
endfor



tmp = fltarr(nut,nalt)
for j=0,nut-1 do tmp[j,1:nalt-1]= out[index[j]].cn2
scale = max(tmp)
tmp = tmp/scale
scale = 1.


for j=0,ntick-1 do tmp(jtick[j],nalt-3:nalt-1)= scale
for j=0,ntick-1 do tmp(jtick[j],0:1)= scale
tmp2 = congrid(tmp, nut,200)
tmp2(*,100) = scale  ; & tmp2(*,200) = scale ; & tmp2(*,300) = scale 


;device, decomposed=0  ; to get the colors on X-display!
; loadct, 3  ; heat
; loadct, 4  ; BLUE/GREEN/RED/YELLOW
  tvscl, 1. - sqrt(tmp2)
; tvscl, congrid(tmp,400, 400, /interp) 
;  tvscl, sqrt(congrid(tmp,nut, 400)) 
 
; device, filename='---',/color, bits=8 ; to get color plots.

;stop

end