 
     lat = (j) -> Math.PI / 256 * (128 - j)
 
-    #avgt.init = (263.15 for i in [0...256])
-    #avgt.init = (253.15 + 40 * Math.cos(lat(i)) for i in [0...256])
-    avgt.init = (303.15 for i in [0...256])
+    rand = -> 1 + 2 * (Math.random() - 0.5) / 100
+
+    avgt.init = (273.15 for i in [0...256])
+    #avgt.init = (263.15 + 20 * Math.cos(lat(i)) for i in [0...256])
+    #avgt.init = (303.15 for i in [0...256])
 
     avgt.evolve =  (data, energyIn) ->
         [gland, land, gair, air] = data

            0.7
        else
            0.3
    shtAbsorbLand = 150 / 342
    shtAbsorbAir  =  50 / 342
    lngAbsorbAir  =  370 / 390
    landOutputRatio = (390 + 24 + 78) / 390

    radius = p.radius

    lcapacity = (t) ->
        if t < 263.15
            2106
        else if t < (263.15 + 2)
            335000 / 2
        else if t < 373.15
            4218
        else if t < (383.15 - 2)
            2250000 / 2
    acapacity = (t) -> 1005

    shtAbsorb = (t) ->
