- This is a Python code to numerically compute the solution of the linear stability analysis problem for the detachment-limited landscape evolution model (LEM).
- Using the developed code, the role of m and n (exponents of the specific drainage area and slope, respectively in the fluvial erosion term of the model) on the onset and selection of regular valley spacing can be examined.
The submitted research work can be found here: https://www.essoar.org/doi/10.1002/essoar.10511126.1.
-
numpy and scipy are the required python libraries.
-
This repository contains the code as a jupyter notebook, which can be installed using pip as
pip install notebook
or using Condaconda install -c conda-forge notebook
. -
Files s200.mat and w200.mat in the repository contain 200 Quadrature Points and Weight Functions, respectively, for the numerical solver.
Each number written in the list below explains the operations and functions for the particular Cell of the Jupyter Notebook.
- Import relevant packages from numpy and scipy libraries.
- Import quadrature points and weight functions arrays (s200.mat and w200.mat), define Python functions for trial and test functions used in the weak formulation, and define Python functions to numerically obtain the slope and its derivatives for generic values of m and n.
- Select values of parameters (CI,m,n,L,U,k), global variables and arrays used in the solver.
- For every Channelization Index (CI) and wavenumber (k), solve the eigenvalue problem using the spectral Galerkin technique with numerical quadrature.
- Extract the critical Channelization Index CIcr and the corresponding fastest growing (positive) spatial frequency (kcr).
- Save relevant arrays for plotting and further investigations.
For more information about this research, you can reach out to me: Shashank Kumar Anand.
Well-commented Python code for the numerical simulations of the Landscape Evolution Model (LEM) by the same author is available at github.com/ShashankAnand1996/LEM. The developed code has been tested to provide efficient and accurate solutions for the LEM in detachment-limited conditions. The published manuscript is also available: Linear layout of multiple flow-direction networks for landscape-evolution simulations.