         land = avgt.init
         air = avgt.init
         evolve = ->
-            [time, x, v] = step(time, x, v, au.fromSI_T(p.period))
+            [time, x, v] = step(time, x, v, au.fromSI_T(planet.period))
             [x1, y1, x2, y2, x3, y3] = x
             [vx1, vy1, vx2, vy2, vx3, vy3] = v
             phase = [x1, y1, vx1, vy1, x2, y2, vx2, vy2, x3, y3, vx3, vy3]
             self.postMessage({orb: phase})
 
-            lumTotal = r.total(time, x)
+            lumTotal = rad.total(time, x)
             self.postMessage({lum: lumTotal})
 
-            lumA = r.a(time, x)
-            lumB = r.b(time, x)
+            lumA = rad.a(time, x)
+            lumB = rad.b(time, x)
             twilight = (
                 (
                     [

    'cs!/wahlque/universe/wahlque/starB'
    'cs!/wahlque/universe/wahlque/planet/planet'
    'cs!/wahlque/universe/wahlque/planet/radiation'
    'cs!avgtfield'
], (exports, solver, vec3, b3, au, a, b, p, r, avgt) ->
    handle = 0
    time = 0
    start = () ->
        m1 = a.mass
        m2 = b.mass
        acceleration = b3.acceleration(au, m1, m2)
        step = solver.solve(acceleration)

        [x1, y1] = a.initPosition
        [vx1, vy1] = a.initVelocity
        [x2, y2] = b.initPosition
        [vx2, vy2] = b.initVelocity
        [x3, y3] = p.initPosition
        [vx3, vy3] = p.initVelocity
        x = [x1, y1, x2, y2, x3, y3]
        v = [vx1, vy1, vx2, vy2, vx3, vy3]

        l1 = a.luminosity
        l2 = b.luminosity

        lng = (i) -> 2 * Math.PI / 256 * i
        lat = (j) -> Math.PI / 256 * (128 - j)

            self.postMessage({twlt: twilight})

            data = [global, land, global, air]
            [gland, land, gair, air] = avgt.evolve(data, rad.averageIn(time, x))
            self.postMessage({tmp: [gland, land, gair, air]})
            self.postMessage({msg:
                '<ul>' +
