Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/brews/erebusfall
Ice-volume correction to marine-isotope proxy records in Python.
https://github.com/brews/erebusfall
d18o deuterium ice-volume-correction marine-isotope paleoceanography paleoclimate python
Last synced: about 1 month ago
JSON representation
Ice-volume correction to marine-isotope proxy records in Python.
- Host: GitHub
- URL: https://github.com/brews/erebusfall
- Owner: brews
- License: gpl-3.0
- Created: 2017-11-07T17:15:34.000Z (about 7 years ago)
- Default Branch: master
- Last Pushed: 2019-01-19T03:11:32.000Z (almost 6 years ago)
- Last Synced: 2024-10-09T15:42:31.330Z (about 1 month ago)
- Topics: d18o, deuterium, ice-volume-correction, marine-isotope, paleoceanography, paleoclimate, python
- Language: Python
- Homepage:
- Size: 60.5 KB
- Stars: 2
- Watchers: 3
- Forks: 0
- Open Issues: 2
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
# erebusfall
[](https://travis-ci.org/brews/erebusfall)
Apply a simple ice-volume correction to marine-isotope
[δ18O](https://en.wikipedia.org/wiki/%CE%9418O) and [δD](https://en.wikipedia.org/wiki/Deuterium)
proxy records, in Python.
The correction uses the [LR04](https://doi.org/10.1594/PANGAEA.701576)
benthic stack from [Lisiecki and Raymo, 2005](https://doi.org/10.1029/2004PA001071). The stack is scaled so that the [LGM](https://en.wikipedia.org/wiki/Last_Glacial_Maximum)-to-present
change is assumed to be 1.0 ‰ in accordance with the pore-water estimate of
[Schrag et al. 1996](https://doi.org/10.1126/science.272.5270.1930). The
package is adapted from the ice-volume correction in [Tierney et al. 2017](https://doi.org/10.1130/G39457.1).
## Example
Start by importing `erebusfall`, `numpy`:
```python
import erebusfall as ef
import numpy as np
```
...and creating a hypothetical proxy time
series...
```python
age_ka = np.arange(0, 20, 1)
proxy = np.random.normal(loc=-2.0, size=len(age_ka))
```
Now we can plug this into `ef.icevol_correction()` with a few key options:
```python
proxy_adjusted = ef.icevol_correction(age_ka, proxy,
proxytype='d18o',
timeunit='ka')
```
We first plug in the proxy age and proxy values. The `proxytype='d18o'`
indicates that we're dealing with δ18O. We can set `proxytype='dd'` for a δD
record. The `timeunit` argument indicates that `age_ka` is in thousands of
years before present. There are also options for for "years BP" and
"million years BP". The output from the function, `proxy_adjusted`, is the
corrected isotope proxy as a numpy array.
See `help(ef.icevol_correction)` for more documentation.
## Installation
You can install the package [from PyPI](https://pypi.python.org/pypi/erebusfall) with
```
pip install erebusfall
```
If you want to use `conda`:
```
conda install erebusfall -c sbmalev
```
## Development and Support
Source code is [hosted online](https://github.com/brews/erebusfall) under an Open
Source license. Please feel free to file any
[bugs and issues](https://github.com/brews/erebusfall/issues) you find.