Parametric steady state solar pond analyses—a numerical approach