Simulating Precipitation with ADCIRC-Hydrology

This article walks through a workflow for acquiring rainfall data and using it with a fork of ADCIRC that includes a source term. This was last updated for version X. 

Acquisition of Data

The code can accept data from any source so long as it is formatted 'OWI-style'. In particular, there are two modes of operation: 

1. Global precipitation intensity - reads from a fort.425 ascii file with data for every node. 
2. Local precipitation intensity - reads from a fort.27 metadata file and a fort.425 data file. This applies precipitation over a specified box. 

The following script and required module from Maria Contreras downloads NCEP rainfall data, reprojects the data, and rewrites the files to proper format for us. It requires the gdal package which can be installed through conda.

Pp_prep_storm.py

OM2D_functions.py

There is an experimental version in that works that will hopefully run faster as the python version is pretty slow right now. Check here in the future for updates. 

Running 

To include precipitation, all that is needed as a control line at the bottom of your fort.15. For example: 

&rainfallcontrol NCRAIN=2, CRAIN_TIMINC=3600

Here, the NCRAIN option controls the mode of operation: 2 for local data, 12 for global data, and negative for indicating that the rainfall begins at hotstart time. CRAIN_TIMINC controls the time delta in between successive rain data. 

Related articles

• Page:
  Simulating Precipitation with ADCIRC-Hydrology