Numerical model of the Messinian Mediterranean combining hydrological water balance, river erosion, and flexural isostasy: TISC code and input dataset for the Lago-Mare

We hereby share the code and input/output datasets for the numerical model
of the Messinian Salinity Crisis published in the Science Advances paper
linked to this submission. The Messinian salinity crisis (MSC) was a short
period of isolation of the Mediterranean Sea causing the precipitation of
a million cubic kilometers of salt. The puzzling sedimentary record that
formed after salt deposition yields conflicting levels of desiccation,
ranging from a full exposure of most of the Mediterranean seafloor based
on shallow sediment found in the abyss to a nearly full Mediterranean
scenario as suggested by similar, fresher-water deposits ubiquitous along
the coastline: the so-called Lago-Mare formation. We use the TISC
landscape evolution model (integrating fluvial drainage, hydrological
balance in endorheic lakes, fluvial erosion, and flexural isostasy) to
explore the expected evolution for the drawdown stage of the Messinian
salinity crisis, known as the 'Lago-Mare'. The model is
constrained by published climate and sediment budget estimates. 
The results show that the propagation of an erosional wave into the
surrounding continents added a gradual sea level rise to the climatic
oscillations of the Mediterranean. This retrogressive fluvial incision of
the spillways of the Paratethys and the Pannonian basins also explains the
Mediterranean transition to fresher water conditions during the MSC.