print, '               df.ps         (0)'
print, '               df.eps        (1)'
read, medium

case medium of
  0: begin
       set_plot, 'PS'
       device, filename='df.ps', $
          /landscape,  $
            xoffset=2.5, yoffset=26., $
            xsize=20., ysize=15.
     end
  1: begin
        set_plot, 'PS'
	device, filename='df.eps', $
        /encapsul, $
	xsize=20., ysize=15.
     end
endcase

xyouts, 0.5, 0.5, /normal, '!5'

; all sizes = 1.2, all thicks = 2
!x.charsize=1.2
!x.thick=2
!y.charsize=1.2
!y.thick=2
!p.charsize=1.2
!p.charthick=2.0
!p.thick=2

;tit = 'Nominal Descent'
dummy = ' '
fn = ' '
ln = long(0.)
fn = 'fr_contributions.out'
openr, 1, fn
while not eof(1) do begin
  readf, 1, dummy
  ln = ln + 1
endwhile
close, 1

dat = dblarr(19, ln)
openr, 1, fn
readf, 1, dat
close, 1

!p.multi = [0, 1, 2, 0, 0]

!x.range = [0., 160.]
;!x.range = [0., 180.]

f0overc = 20.4/3.

dat[3,*] = -dat[3,*] * f0overc / 1000.
dat[4,*] = -dat[4,*] * f0overc / 1000.
dat[5,*] = -dat[5,*] * f0overc / 1000.
dat[6,*] = -dat[6,*] * f0overc / 1000.
dat[7,*] = -dat[7,*] * f0overc / 1000.
dat[8,*] = -dat[8,*] * f0overc / 1000.
dat[9,*] = -dat[9,*] * f0overc / 1000.
dat[10,*] = dat[10,*] * f0overc / 1000.
dat[11,*] = -dat[11,*] * f0overc / 1000.

!p.region = [0.0, 0.5, 1.0, 1.0]
plot, dat[0,*], dat[3,*], $
  position = [0.1, 0.5, 1.0, 0.9], $
  linestyle = 0, $
  yrange = [-10., 20.], $
  ystyle = 1, $
  ytitle = '  ', $
  xtitle = ' ', $
  xtickname = replicate(' ',30), $
  xstyle = 9

oplot, dat[0,*], dat[10,*] - dat[9,*], $
  linestyle = 5

for alt = 160, 0, -10 do begin
  if (alt le dat[1,0]) and (alt ge dat[1,ln-1]) then begin
    j = long(0.)
    while (not ((alt le dat[1,j]) and (alt ge dat[1,j+1]))) do j = j + 1
    tval = dat[0,j] + (dat[0,j+1] - dat[0,j]) / (dat[1,j+1] - dat[1,j]) * (dat[1,j+1] - alt)
    ttick = [!x.crange[0], tval, !x.crange[1]]
    if alt mod 50.0 then begin
      tlen = 0.01
      if alt ge 30 then begin
        tstring = ' '
      endif else begin
        tstring = strcompress(string(alt), /remove_all)
      endelse
    endif else begin
      tlen = 0.02
      tstring = strcompress(string(alt), /remove_all)
    endelse
    if (tval gt !x.crange[0]) and (tval lt !x.crange[1]) then $
      axis, xaxis = 1, $
        xtitle = 'ALTITUDE [km]', $
	xstyle = 1, $
	xrange = !x.crange, $
	xticks = 2, $
        xtickv = ttick, $
	xticklen = tlen, $
	xtickname = [' ', tstring, ' ']
  endif
endfor

oplot, !x.crange, [0., 0.], thick = 0.05

xyouts, 100.0, 15., 'Doppler Shift from:', charsize = 1.2

oplot, [100.0, 120.0], [12., 12.], $
  linestyle = 0
xyouts, 125., 12., 'Total', charsize = 1.2

oplot, [100.0, 120.0], [9., 9.], $
  linestyle = 5
xyouts, 125.0, 9., 'Orbiter', charsize = 1.2

;xyouts, 90., 15., tit, charsize = 1.6

!p.region = [0.0, 0.0, 1.0, 0.5]
plot, dat[0,*], dat[7,*], $                ; zonal wind LOS
  position = [0.1, 0.1, 1.0, 0.5], $
  linestyle = 1, $
  yrange = [-0.5, 0.05], $
  ystyle = 1, $
  ytitle = ' ', $
  xtitle = 'DESCENT TIME [min]', $
  xstyle = 9
oplot, dat[0,*], dat[9,*], $          ; Titan rotation
  linestyle = 2
oplot, dat[0,*], dat[8,*], $       ; meridional winds
  linestyle = 3
oplot, dat[0,*], dat[5,*], $        ; descent vel.
  linestyle = 4
for alt = 160, 0, -10 do begin
  if (alt le dat[1,0]) and (alt ge dat[1,ln-1]) then begin
    j = long(0.)
    while (not ((alt le dat[1,j]) and (alt ge dat[1,j+1]))) do j = j + 1
    tval = dat[0,j] + (dat[0,j+1] - dat[0,j]) / (dat[1,j+1] - dat[1,j]) * (dat[1,j+1] - alt)
    ttick = [!x.crange[0], tval, !x.crange[1]]
    if alt mod 50.0 then begin
      tlen = 0.01
      if alt ge 30 then begin
        tstring = ' '
      endif else begin
        tstring = strcompress(string(alt), /remove_all)
      endelse
    endif else begin
      tlen = 0.02
      tstring = strcompress(string(alt), /remove_all)
    endelse
    if (tval gt !x.crange[0]) and (tval lt !x.crange[1]) then $
      axis, xaxis = 1, $
        xtitle = ' ', $
	xstyle = 1, $
	xrange = !x.crange, $
	xticks = 2, $
        xtickv = ttick, $
	xticklen = tlen, $
	xtickname = [' ', ' ', ' ']
  endif
endfor

xyouts, 100.0, -0.2, 'Doppler shift from:', charsize = 1.2

oplot, [100.0, 120.0], [-0.25, -0.25], $
  linestyle = 1
xyouts, 125.0, -0.25, 'zonal wind', charsize = 1.2

oplot, [100.0, 120.0], [-0.3, -0.3], $
  linestyle = 2
xyouts, 125.0, -0.3, 'Titan rotation', charsize = 1.2

oplot, [100.0, 120.0], [-0.35, -0.35], $
  linestyle = 3
xyouts, 125.0, -0.35, 'meridional wind', charsize = 1.2

oplot, [100.0, 120.0], [-0.4, -0.4], $
  linestyle = 4
xyouts, 125.0, -0.4, 'Probe descent', charsize = 1.2

xyouts, -20.0, 0.05, 'DOPPLER SHIFT !7d!Xf [kHz]', alignment = 0.5, $
  charsize=1.6, orientation=90.0

!p.multi = 0
!p.charsize=1.
!p.charthick=1.
!p.thick=1.
!x.thick = 1.
!y.thick = 1.
!x.range = 0
!x.style = 0

if (medium eq 0) or (medium eq 1) then device, /close
end