Finite element model for transient two-layer cooling pond behavior