in a user-oriented code called SOLA. Although SOLA does not treat free surfaces, an extended version, SOLA-SURF, is also available [10] that uses the surface height function method (see Section ILA). The basic simplicity and flexibility of the SOLA codes make them excellent foundations for the development of more sophisticated codes.