SELEN: A Fortran 90 program for solving the “sea-level equation”

We present SELEN, a Fortran 90 computer program for solving the “sea-level equation”, the law that governs the sea-level variations associated with the melting of the Pleistocene ice sheets. As shown here by a number of examples, SELEN also allows to study a suite of geophysical processes accompanying the mass redistributions associated with glacial-isostatic adjustment, such as vertical deformations of the solid surface of the Earth, and variations of the shape of the geoid. SELEN is built in a simple and easily testable manner, and requires modest computer resources. The source code, which is freely available on request, can be easily modified to implement new features of the sea-level equation or to focus on specific geodynamical aspects of glacial-isostatic adjustment.