print, 'Ausgabemedium:                     '
print, '                                   '
print, '               output.ps        (0)'
print, '               Tektronix-Schirm (1)'
print, '               X-Schirm         (2)'
print, '               Windows          (3)'
print, '               Programmende     (4)'
read, medium

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

           /portrait, $
;            xoffset=2., yoffset=3., $
;            xsize=16., ysize=24.
            xoffset=0., yoffset=5., $
;            xsize=17.5, ysize=16.
            xsize=19.25, ysize=17.6
     end
  1: begin
       set_plot, 'TEK'
     end
  2: begin
       set_plot, 'X'
;       window, 0, xsize=650.*210/297, ysize=650.
;       !p.color=255
;       !p.background=0
     end
  3: begin
       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'

dummy = ' '

dat = dblarr(3, 180, 91)

cont = dblarr(181, 91)
theta = dblarr(181, 91)
phi = dblarr(181, 91)
x = dblarr(181, 91)
y = dblarr(181, 91)
  
print, 'Welches File soll ich darstellen?'
read, dummy

openr, 1, dummy
readf, 1, dat
close, 1

for i=0, 179 do begin
   theta(i,*) = dat(0,i,*)
   phi(i,*) = dat(1,i,*)
   cont(i,*) = dat(2,i,*)
endfor

cont(180,*) = dat(2,0,*)
theta(180,*) = dat(0,0,*)
phi(180,*) = dat(1,0,*)

x = theta * cos(phi*!DPI/180.)
y = theta * sin(phi*!DPI/180.)

contour, cont, x, y, $
   /follow, $
   charsize = 2.0, $
   c_charsize = 1.0, $
   charthick = 2.0, $
;   title = dummy, $
   title = 'Probe Transmitter Antenna FM!CLCP 2040 MHz!C!CAntenna Pattern [dB]', $
;   title = 'Probe Transmitter Antenna FM!CRCP 2098 MHz!C!CAntenna Pattern [dB]', $
;   xtitle = 'theta * cos(phi)', $
;   ytitle = 'theta * sin(phi)', $
   xtitle = '!7h!5 * cos !7u!5 [deg]', $
   ytitle = '!7h!5 * sin !7u!5 [deg]', $
   xrange = [-90., 90.], $
   xstyle = 1, $
   yrange = [-90., 90.], $
   ystyle = 1, $
;   levels = [-40., -35., -30., -25., -20., -15., -10., -5., 0., 5.]
   levels = [-8., -4., -1., 0., 1., 2., 3., 4., 5.], $
   c_labels = [1, 1, 1, 1, 1, 1, 1, 1, 1]

contour, theta, x, y, $
   /overplot, $
   levels = [10., 20., 30., 40., 50., 60., 70., 80., 90.], $
   c_linestyle = 1

xyouts, 0., -150., /data, '!7h!5 = Elevation Angle [0!9%!5, 90!9%!5]!C!B!7u!5 = Azimuth Angle [0!9%!5, 360!9%!5]', size = 1.5, alignment = 0.5
   
if medium eq 0 then device, /close
end