https://github.com/johnrudge/rivers2d

River inversion
https://github.com/johnrudge/rivers2d

Last synced: 23 days ago
JSON representation

River inversion

• Host: GitHub
• URL: https://github.com/johnrudge/rivers2d
• Owner: johnrudge
• Created: 2022-02-22T16:31:30.000Z (over 2 years ago)
• Default Branch: master
• Last Pushed: 2022-02-22T16:33:06.000Z (over 2 years ago)
• Last Synced: 2023-12-20T19:42:31.142Z (11 months ago)
• Language: Python
• Size: 145 KB
• Stars: 0
• Watchers: 1
• Forks: 0
• Open Issues: 0
• Metadata Files:
  • Readme: README.md

Awesome Lists containing this project

README

        
River Inversion 2D

==================

The files in this repository are a set of python scripts which implement the 

river profile to uplift history inversion alogrithm described in Rudge, Roberts,

White and Richardson (2015) [1]. The code uses the non-negative least square 

(NNLS) algorithm of Mathieu Blondel (https://gist.github.com/mblondel/4421380) 

to perform the inversion.

Installation

============

In order to run these scripts, you need to have the following installed:

* python 3.6

* numpy (1.7 or later)

* scipy (0.12 or later)

* FEniCS 2019.1.0 (www.fenicsproject.org)

* scikit-learn (0.13 or later)

For VTK output of river data, you also need the pyevtk library

installed.

Usage

=====

The main routine is rivers_2d.py, which can be run with

python3 rivers_2d.py

which will invert the example data for Madagascar, contained in the 

madagascar_data folder. All obs_river* files have format:

x (lon, m), y (lat, m), z (elevation, m), d (distance, m), A (area, m^2)

References

==========

[1] Rudge J.F., Roberts G.G., White N., Richardson C.N. Uplift histories of 

Africa and Australia from linear inverse modeling of drainage inventories (2015) 

J. Geophys. Res. Earth Surf. 120:1-21 http://dx.doi.org/10.1002/2014JF003297