adcircpy 0.9.1

Python package for working with ADCIRC input and output files.

AdcircPy

A Python interface for handling inputs and outputs for the ADCIRC hydrodynamic model.

Basic Installation:

You can install this software through conda or through pip.
The dependecies of Python packages installed through pip must be satisfied by the operating system (see Method 2 for details). On a normal user's computer, the recommended method of installation is through conda. Whether you are using the pip method or the conda method, it is generally good practice to always make use of a virtual environment for your project. If you are a Python newcomer, make sure to read upon Python environments and their usage as you dive more into Python development.

If you are doing development on a Windows machine it is highly recommended to use the Windows Subsystem for Linux, and run this software from a Widows-Subsystem-for-Linux shell. MingW64, Cygwin, the git-shell-for-windows, or any other off-the-shelf POSIX layer on Windows should work.

Method 1: Conda

The advantage of using conda is that it provides the necessary dependencies as precompiled binaries, therefore no local compilations are required, and cuts the installation time substantially. This is the recommended method to use.

Install Miniconda3 (recommended) or Anaconda3 on your system. Rememeber to make this installation on the POSIX layer when using Windows.

Create a new conda environment for your project, and inside the environment run:

conda update -n base -c defaults conda
conda install -c conda-forge gdal
pip install AdcircPy

Method 2: Pip-only (requires compilation of system libraries).

Pip requires that some system libraries are precompiled and installed before running pip install command. It is highly recommended that you use the system packge manager to satisfy these dependencies, or, they may be compiled from source. The .docker directory contains a full piecewise build on an Alpine distribution, which may be used for reference in case further guidance is required with respect to the full compilation procedure. Note that "dev" versions mean that the headers are required along with the compiled library. Also note that some of these libraries you may already have in your system.

Full system dependency list:

• g++
• gcc
• m4
• make
• autoconf
• zlib-dev
• libc-dev
• curl-dev
• linux-headers
• hdf5-dev
• netcdf-c-dev
• python3-dev
• sqlite-dev
• Proj4-dev
• GDAL-dev
• openblas-dev (for scipy)
• freetype-dev (for matplotlib)
• py3-pip

Once these dependencies have been met, you may simply do

pip install AdcircPy

NOTE: If you are developer and would like to debug the package, you may install this package on developer mode by first cloning the package, and then doing:

pip install -e .

Usage

CLI interface

After installation a few new commands are available on the command line interface.

A few of the available commands are:

• PlotMesh
• GenerateTidalRun
• GenerateBestTrackFile
• PlotMaxele
• HighWaterMarkValidation
• PlotTidalStationsOutput

Since this package is in development, the list above may not be up to date, therefore you may check setup.py, where the entrypoints for the present release are listed.

As an API

This package is an API to handle input and output files for the ADCIRC hydrodynamic model. You can load an Adcirc mesh and plot it by doing:

from AdcircPy import read_mesh
mesh = read_mesh('/path/to/fort.14')
# The most basic you can do with this software is create plots:
mesh.make_plot(show=True)

Example of fort.15 generation:

Tidal only run:

from AdcircPy import read_mesh
from datetime import datetime, timedelta
mesh = read_mesh('/path/to/fort.14', epsg=4326, vertical_datum='LSML')
start_date = datetime.now()
end_date = start_date + timedelta(days=5)
tidal_run = mesh.TidalRun(start_date, end_date, spinup_days=7)
tidal_run.dump()  # will write to stdout if no output path is given.

Example where global outputs are requested:

from AdcircPy import read_mesh
from AdcircPy import ElevationGlobalOutput as EGO
from datetime import datetime, timedelta
mesh = read_mesh('/path/to/fort.14')
start_date = datetime.now()
end_date = start_date + timedelta(days=5)
tidal_run = mesh.TidalRun(
            start_date, end_date,
            ElevationGlobalOutput=EGO(sampling_frequency=15))
tidal_run.dump()

* Examples will be updated as development continues. *

Please report bugs to [e-mail]