print, 'Output Medium:                      '
print, '                                    '
print, '               output.ps         (0)'
print, '               output.eps        (1)'
print, '               Tektronix-Screen  (2)'
print, '               X-Screen          (3)'
print, '               Windows           (4)'
print, '               terminate program (5)'
read, medium

case medium of
  0: begin
       set_plot, 'PS'
       device, filename='output.ps', $
         /landscape,  $
           xoffset=2.5, yoffset=26., $
           xsize=24., ysize=16.

;           /portrait, $
;            xoffset=2., yoffset=3., $
;            xsize=16., ysize=24.
     end
  1: begin
        set_plot, 'PS'
	device, filename='output.eps', $
	/encapsulated, $
	xsize=24., ysize=16.
     end
  2: begin
       set_plot, 'TEK'
     end
  3: begin
       set_plot, 'X'
;       window, 0, xsize=650.*210/297, ysize=650.
;       !p.color=255
;       !p.background=0
     end
  4: begin
       set_plot, 'WIN'
;       window, 0, xsize=600., ysize=400.
;       !p.color=255
;       !p.background=0
     end
  else: begin
          print, 'Programmende'
          stop
        end
endcase

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

temp = dindgen(101)
tau  = dindgen(101)
p    = dindgen(101)

temp = temp - 30.0

;tau = 0.0035 * (temp - 85.) * (temp - 85.) + 4.0
;p = 9.93 - 0.047 * temp

tmt = [55.2, 39.6, 24.8, 11.2, -10.0, -24.0]
pct = [6.81, 6.94, 7.5, 7.92, 8.61, 9.27]

tmq = [58.0, 39.0, 24.0, 10.0, -11.0, -30.0]
pcq = [6.26, 6.64, 7.1, 7.5, 8.0, 8.75]

polyp = svdfit((tmt+tmq)/2., (pct+pcq)/2., 3)

print, polyp(0)
print, polyp(1)
print, polyp(2)

p = polyp(0) + polyp(1) * temp + polyp(2) * temp * temp

;oplot, temp, p, linestyle = 1
;xyouts, temp(20), p(20), 'P2', /data

plot, temp, p, $
  linestyle = 0, $
  xrange = [-30.0, 60.0], $
  xstyle = 1, $
  yrange = [6.0, 13.0], $
  ystyle = 9, $
  xtitle = 'TEMPERATURE [!9%!5C]', $
  ytitle = 'POWER [W]'

xyouts, temp(10), p(10), 'P', /data

oplot, tmt, pct, $
  psym = 1
oplot, tmq, pcq, $
  psym = 1

axis, yaxis = 1, $
  yrange = [0.0, 40.0], /save, $
  ystyle = 1, $
  ytitle = '!7s!5 [sec]'

;oplot, temp, tau, $
;  linestyle = 2
;
;xyouts, temp(10), tau(10), $
;  '!7s!5', /data

tmt  = [55.2, 39.6, 24.8, 11.2, -10.0, -24.0]
wtmt = [8., 12., 15.+55./60., 19.+48./60., 25.+40./60., 38.+43./60.]

tmr  = [55., 39.3, 24., 11., -10.7, -24.0]
wtmr = [8.+15./60., 11.5, 15.+51./60., 23.+46./60., 26.+52./60., 34.+35./60.]

polytau = svdfit((tmt+tmr)/2., (wtmt+wtmr)/2., 3)

print, polytau(0)
print, polytau(1)
print, polytau(2)

tau = polytau(0) + polytau(1) * temp + polytau(2) * temp * temp
oplot, temp, tau, linestyle = 2
xyouts, temp(10), tau(10), '!7s!5', /data

oplot, tmt, wtmt, $
  psym = 2
oplot, tmr, wtmr, $
  psym = 2

oplot, [20.], [35.5], psym = 1
oplot, [20.], [32.5], psym = 2
xyouts, 22., 35., /data, $
  'Steady State Power Consumption'
xyouts, 22., 32., /data, $
  'Warm-Up Time until Lock'

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