compute_recharge function

Simulation using a recharge model