GeoClaw with WRF

Overview

This repository is developed for storm surge simulation of geoclaw with user-defined wind and SLP fields.
Here is an example of the simulation with uniform and constant wind field which blows from west to east (plotted by visclaw_julia ).

Requirements

• clawpack v5.4.1
  This repo doesn't support clawpack v5.5.0 yet.
• NetCDF4 libraries (netcdf-c and netcdf-fortran)

Usage

• install the required libraries and packages
• set the environmental variables
  You can check whether the variables are set appropriately by running the following commands in your shell.

echo $CLAW # top directory of your clawpack
echo $PYTHONPATH # top directory of your clawpack
echo $FC # Fortran compiler
echo $NETCDF_F_ROOT # directory of your netcdf-fortran

• modify data.py
  The variable storm_type is updated and the new type 4 is defined.
  In order to avoid errors caused by this update, run

./replace_datapy.sh

Then $CLAW/geoclaw/src/python/geoclaw/data.py was updated and original data.py was saved as data.py_org.

• run a test simulation
  Test simulations can be implemented by

make clobber
make .exe
make storm
make .data
make .output

Ongoing and Future Work

Wind and SLP should be dumped at open boundaries when the storm field covers the boundaries of simulation.
The sources on boundary condition are now being developed to adopt this dumping effect.

License

BSD 3-Clause

Acknowledgements

The author would like to thank Marc Kjerland Ph.D. for collaboration on the early stages of this work.

Author

• Takuya Miyashita miyashita (at) hydrocoast.jp