Skip to content

A python package for calculating a numerical solution to the plume equations, designed primarily for buoyant plumes of subglacial discharge at ocean-terminating glaciers.

License

Notifications You must be signed in to change notification settings

soroushr/pyplume

 
 

Repository files navigation

pyplume

DOI

Pyplume calculates a numerical solution to the equations formulated for a point source plume released at the base of a vertical wall. The code includes capabilities for calculating the melt rate in the case where the vertical wall is made of ice.

See the examples and docstrings for more details on how to use pyplume. The links below will launch interactive jupyter notebooks with Binder so you can see how the code works and experiment with the examples. Here's a link explaining how to interact with jupyter notebooks if you haven't used them before.

The first release (v0.1.0) contains the code used in Halbach et al. (submitted). Example 3 is a simplified, interactive verion of the code used in the paper, with a synthetic discharge timeseries. The full code used in the paper is available in the examples folder as kronebreen_upwelling.py. The full datasets as used in the paper are available upon request from the authors.

Examples:

1. Simple plume

A minimal example demonstrating how to use pyplume with linear ambient temperature and salinity profiles and plot the results.

Binder

2. Comparison to analytical solution

A slightly more detailed look at using pyplume and comparing the results to an analytical solution assuming uniform ambient temperature and salinity, with the inlet Richardson number set to $\Gamma_0=1$.

Binder

3. Discharge upwelling timeseries

An example demonstrating how pyplume can be used with a discharge timeseries to calculate how the plume behaves on seasonal timescales. The example also demonstrates a method which can be used to estimate plume driven upwelling, as described in Halbach et al. (submitted).

Binder

More examples and documentation will be added in future.

Acknowledgements

The code was developed by Alistair Everett with inspiration from Tom Cowton's iceplume Fortran model (Cowton et al., 2015). Some of this development work was done as part of the TIGRIF (RCN project number: 243808/E40) and TW-ICE research projects at the Norwegian Polar Institute.

This code relies upon many previous studies which derived the solutions and equations built into this model, references have been included throughout the code. If you use pyplume, in addition to citing this code, please also consider citing some of the important papers on which this code is based.

About

A python package for calculating a numerical solution to the plume equations, designed primarily for buoyant plumes of subglacial discharge at ocean-terminating glaciers.

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Languages

  • Python 100.0%