GeoClaw input files of the 2004 Sumatra-Andaman tsunami scenarios published in Ulrich et al. (2021)

This dataset contains the input files of the GeoClaw scenarios of the 2004 Sumatra-Andaman tsunami presented in:

Ulrich, T., Gabriel, A. A., Madden, E. H. (2021). Stress, rigidity and sediment strength control megathrust earthquake and tsunami dynamics. doi: 10.31223/osf.io/s9263.

runconverterLMU_WGS84.sh: contains all the steps to transform a SeisSol surface output to a Geoclaw tt3 file.
It uses the displacement converter of Samoa to rasterize a SeisSol surface output to NetCDF.
See README_build_displacement-converter.txt for the procedure to download and build the displacement converter.
Note that the SAMPLER (https://github.com/SeisSol/SAMPLER) will replace the displacement-converter in the future.

convert_geographic_SeisSol_geom.py: to transform the geometry array of a SeisSol surface output file to the geocentric coordinate system (latitude, longitude).
tapperNetcdf.py: to apply a Hanning window on a NetCDF file. This prevents sharp displacement discontinuities at the limits of the region of imposed displacements, which could generate spurious waves.
convert_netcdf_tt3.py: to convert a NetCDF displacement file to the tt3 format (GeoClaw).

The GeoClaw simulations require the following files:

displacement_tt3_files.tar.gz : rasterized input files in tt3 format for 4 earthquake scenarios.
gebco_2019_n25.0_s-21.0_w55.0_e110.0.nc: input bathymetry and topography data in NetCDF format downloaded from https://www.gebco.net/.
setrun.py which defines the simulation parameters.
a Makefile, plateform specific see e.g. https://github.com/clawpack/geoclaw/blob/master/examples/tsunami/chile2010/Makefile
setplot_fig4.py: configures GeoClaw for generating outputs for figure 4.                  
setplot_animation.py: configures GeoClaw for generating outputs for the supplementary animations.      
GeoClaw simulations are run with  `make .plots`.