Editing Module:Sandbox/Wnt

local p = {}

increment = 10 -- in meters  Basis is meters and kilograms throughout.
G = 6.674E-11 -- 6.674×10−11 N⋅m2/kg2
M = 1E24

function localg (r)
	-- black hole has gravity GM / r2
	return M * G / (r*r)
end

function deltapressure (pressure, gravity)
	-- for a given global increment, return the increase in pressure in atm, pressure in gs
    -- density of "air" per pressure in atm is (1.225 kg/m3) / atm
    -- pressure of air is 101325 newton / m2 at 1 atm
    -- newton weight is kg * gravity
    -- so x atm of pressure weigh x * gravity * 1225 kg per kilometer per square meter
    return pressure * gravity * 0.001225 * increment
end

function p.main (frame)
	-- we SHOULD start at the "exobase", which on Earth is roughly 7000 km from the center
	-- I am not finding decent figures for exosphere pressure
	-- for now let's calculate from an arbitrary (too high) 1 pascal and see from there what the dependence is
	-- exobase equivalent is where GM/r = same as on Earth, i.e. r = 7000 km * M / mass of Earth
	local r = 7000000 * M / 5.97237E24 -- kg, mass of Earth
	local p = 1 -- arbitrary 1 pascal
	local it = 0
	output = "radius gravity pressure"
	repeat
		it = it + 1
	    p = p + deltapressure (p, localg(r))
	    r = r - increment
	    output = output .. r .. "," .. localg(r) .. "," .. p .. "<br />"
	until (it > 100 or r<0)
	return output
end
return p