def calculate_new_surface_params(iL, temp, mass, dt):
    N = estimateParticlesBBN(mass*BDATMMFRAC) # assume all energy is distributed through the atmosphere
    Q = iL * dt
    dT = Q / ((3/2) * N * kb)
    return temp + dT