Examples

This page provides a short example and a list of example applications in Jupyter Notebooks. Examples in the form of Python scripts can also be found on the examples directory on GitHub.

Below is a list of Jupyter Notebooks with code, comments and figures:

• 1. A Basic Model
• 2. Pastas TimeSeries
• 3. Response functions
• 4. Model Diagnostic Checking
• 5. Adding river levels
• 6. Time Series Analysis with Pastas
• 7. Adding groundwater pumping
• 8. Synthetic Data
• 9. Response function comparison
• 10. Pastas Project
• 11. Adding Multiple Wells
• 12. Testing MCMC on Pastas Models
• 13. Non-linear recharge models
• 14. Threshold Non-linearities
• 15. ARMA(1,1) Noise Model for Pastas
• 16. Reading Dutch Datasets
• 17. Autocorrelation for irregular time series

Tip

The latest versions of the Jupyter Notebooks can be found in the examples folder on GitHub!

Short Example

Below is an example of a short script to simulate groundwater levels (the csv-files with the data can be found in the examples-directory on GitHub):

import pandas as pd
import pastas as ps

oseries = pd.read_csv('data/head_nb1.csv', parse_dates=['date'],
                      index_col='date', squeeze=True)
rain = pd.read_csv('data/rain_nb1.csv', parse_dates=['date'],
                   index_col='date', squeeze=True)
evap = pd.read_csv('data/evap_nb1.csv', parse_dates=['date'],
                   index_col='date', squeeze=True)

ml = ps.Model(oseries)
sm = ps.RechargeModel(rain, evap, ps.Gamma, name='recharge')
ml.add_stressmodel(sm)

ml.solve()
ml.plots.decomposition()

Visualisation 1. A Basic Model