pro multi
    
    ; Bo Milvang-Jensen (milvang@astro.ku.dk)
    ; For the RelViz Web exhibition, June 1996
    ; Last rev.: 07-Jun-1996

    ; Black hole properties
    M = 1.0
    a_over_M = 0.998

    ; Initial position of photon
    r = 5.0
    th = 0.0
    ph = 0.0
    t = 0.0

    ; Initial direction of photon
    dr = -0.22
    dth = 0.0
    dph = -0.2

    ; Affine parameter control
    dl = 0.2
    nl = 300

    ; 3D scene format
    scene_format = 'Inventor'

    ; Drawing style of geodesic
    pointset = 0 ; 1 means Line

    ; Color of geodesic
    cr = 0.0
    cg = 1.0
    cb = 1.0

    ; Coordinate system widget
    coo = 'off'
    coo_length = 5.0

    ; Horizon widget
    horizon = 'off'
    auto_radius = 'on'
    manual_radius = 1.0 ; Only used if auto_radius = 'off'
    hcr = 1.0
    hcg = 1.0
    hcb = 0.0
    transp = 0.50

    ; Static Limit / Ergosphere widget
    make_r_stat = 'off'
    scr = 0.0
    scg = 1.0
    scb = 1.0

    ; The rest
    out_file = '/m/RelViz/WWW/milvang/test/multi.iv'
    append = 0 ; 0 means do not append
    verbose = 0 ; 0 means not verbose

    nrays = 20

    p1 = 0.0
    p2 = 1.3

    cr1 = 0.0 & cr2 = 0.0
    cg1 = 1.0 & cg2 = 0.0
    cb1 = 1.0 & cb2 = 1.0

    for i=0,nrays-1 do begin
    if (i gt 0      ) then append = 1 ; We are after the first round
    if (i eq nrays-1) then begin      ; We are at the last round
      coo = 'on'
      horizon = 'off'
      make_r_stat = 'off'
    end
    th = p1  + i*p2 /(nrays-1)
    cr = cr1 + i*(cr2-cr1)/(nrays-1)
    cg = cg1 + i*(cg2-cg1)/(nrays-1)
    cb = cb1 + i*(cb2-cb1)/(nrays-1)
    print, 'th = ', th
    script1, M, a_over_M, r, th, ph, t, dr, dth, dph, dl, nl, $
	     scene_format, pointset, cr, cg, cb, coo, coo_length, $
	     horizon, auto_radius, manual_radius, $
	     hcr, hcg, hcb, transp, make_r_stat, scr, scg, scb, $
	     out_file, append=append, verbose=verbose
    end
end