Possible Projects for GeoClaw Hackathon

Feel free to add things below for the GeoClaw Hackathon 2025

Or even better: Briefly summarize here, and then link to an issue in this repository or in some other Clawpack repository, creating a new one if necessary (e.g. in geoclaw/issues or another repository such as visclaw if that's more relevant). Having it as an issue makes it easier for other to join the discussion about it, and makes sure it doesn't get lost after the worshop ends.

Possible topics for the hackathon:

• Add sediment transport to D-Claw (and GeoClaw?)

• Question about conservation of mass for a synthetic depositional problem (Katy will explain)

• More general approach for defining arrival time in Geoclaw (e.g., based on velocity rather than altitude, or similar). Dave has a branch already. There are some tricky details for lahars into lakes.

• Add precipitation as a standard feature in GeoClaw and D-Claw (Katy)

  • Create a general 'dq' feature, similar to dtopo that allows within-simulation changes to q. This could be used for adding precipitation, among other things.
  • Update the standard D-Claw variables to include a storage layer for water thinner than drytol (see photo of notes on Katy's whiteboard).

• Consider putting somethings Dave has developed into geoclaw python libraries (or digger, or elsewhere) (Katy, Dave)

  • regions to GIS files
  • FGMax to geotif or netcdf

• PyVista/GeoVista visualization tools in 3D, see visclaw/issues/311 and some more recent examples with structures and debris tracking. https://www.scipy2025.scipy.org/

• 2D and 3D visualizations using Julia VisClaw.jl (Takuya)

• Add more examples to geoclaw/examples, e.g.

  • Periodic in longitude on sphere, with GeoVista plotting as in this example
  • Illustrating restart capabilities
  • More complete example of realistic tsunami from source to inundation
  • A volcanic tsunami event, such as in Sandanbata et al. (2024), using BoussClaw (Takuya)
  • Tsunami run-up tests with dispersive effects on a simple sloping beach (e.g., Benchmark Methods for Tsunami Model Validation and Verification) (Takuya)
  • Storm surges or meteotsunamis under pre-computed atmospheric conditions using the OWI_integration branch (Takuya)

• Adopt xarray and rioxarray more heavily, and/or adopt pooch for fetching remote data files.

• Shallow water on the cubed sphere with source terms - problems with ghost cell interpolation at cubed sphere patch boundaries (Donna)

• Discussion of coupling features in ForestClaw. Possible avenue to coupling GeoClaw with other applications? (Donna)

• Serre-Green-Nagdi (SGN) solver (Donna)

• Discuss eventual SGN+D-Claw code

• Python tools for running large (>50K) GeoClaw sweeps, focusing on tsunamis generated by random earthquakes, to expand into storm surge / asteroid impact applications.

  • Monitoring, diagnostic, bug checking, etc. (Donsub)
  • Basic scripts for accessing the generated datasets using PyTorch ML models (Donsub)

• Create brief addition to documentation that describes the typical file structure of an application directory (Katy).

• Creating a new Docker image(s) for Clawpack/GeoClaw, in particular one with PETSc installed (work on Friday with Jed). See https://github.com/clawpack/docker-files/pull/15.

• The GeoClaw BMI developed here for the National Water Model, and possible applications for connecting GeoClaw to other models.

  See also this paper: E. Hamidi et al., Coupling Coastal and Hydrologic Models through Next Generation National Water Model Framework. J. Hydrol. Eng. 30, (2024). https://doi.org/10.1061/JHYEFF.HEENG-6343

• Discuss switching to MyST for examples and/or documentation.

• Merge Keh-Ming's Mfluid code into current AMRclaw framework for adaptivity and parallelization

• Help Carlos with filpatch problem in new hyperbolic formulation with dispersion

General housekeeping things we could tackle:

• Fix regression testing using pytest, add more regression tests

• Are gauges properly working (output at initial and final times? info for surges?)

• Are flagregions checked based on cell centers or edges? (Causes problems in regression testing)

• Test if GeoClaw works with more than 2 ghost cells?

• Think about ways to make num_eqn aka nvar aka meqn a constant value instead of variable using iso-c bindings. Test if this speeds up the code.

• Speed up limiters using a "merge" command.

• Add priority ordering option when specifying topofiles to make sure the right one is used where overlapping.

Requests from users or developers who can't attend...

• ???