 
             aoutput = aoutputBtm + aoutputUpr
 
-            ltransfer = 1000 * lcapacity(1) * 0.3 * transfer(land, i) * dA(lat(i), 10) * 100
-            atransfer = 1 * lcapacity(1) * 1 * transfer(air, i) * dA(lat(i), 5000) * 100
+            ltransfer = 1000 * lcapacity(1) * 0.3 * transfer(land, i) * dA(lat(i), 7) * 10
+            atransfer = 1 * lcapacity(1) * 1 * transfer(air, i) * dA(lat(i), 4000) * 10
 
-            lradius = linput - loutput + ltransfer
-            aradius = ainput - aoutput + atransfer
+            lradius = (linput - loutput + ltransfer) * rand()
+            aradius = (ainput - aoutput + atransfer) * rand()
 
-            nland[i] = land[i] + lradius / lcapacity(land[i]) / dV(10) / 1000 * period
-            nair[i] = air[i] + aradius / acapacity(air[i]) / dV(5000) / 1 * period
+            nland[i] = land[i] + lradius / lcapacity(land[i]) / dV(7) / 1000 * period
+            nair[i] = air[i] + aradius / acapacity(air[i]) / dV(4000) / 1 * period
         [gland, nland, gair, nair]
 
     avgt

            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) ->
