Tsunami Simulation: TsunAWI Testcase "Coquimbo 2015", Setups and Example Results

This dataset contains three setups with example results for TsunAWI simulations for the Chilean coast with a focus on
the Coquimbo region. They are used in the LEXIS project, WP6 "Earthquake and Tsunami Pilot".

coquimbo_150m.tgz - coarse setup with down to 150m resolution in Coquimbo and other areas of interest along the Chilean coast

coquimbo_20m.tgz - fine setup with down to 20m resolution in Coquimbo, also high resolution along other areas of interest along the Chilean coast

coquimbo_20m_new.tgz - setup tuned for the LEXIS workflow aiming at a real time estimate of the inundation on the City of Coquimbo, thus a focus with 20m resolution only in Coquimbo region.

Each dataset includes a README explaining the included files. The following description is taken as pars pro toto from coquimbo_20m_new.tgz.

Computational Mesh

mesh/Coquimbo_20m/nod2D.out 
Location of mesh vertices in geographical coordinates, vertices on the open boundary are flagged with "1"
Resolution in Coquimbo area: 20m

mesh/Coquimbo_20m/nodhn.out
Water depth at the respective vertices. Negative values denote elevation on land. 
In the ocean: ‘The GEBCO_08 Grid, version 20090202 (http://www.gebco.net)’ 
On land: SRTM, Shuttle radar topography mission, 1 arcsec resolution, https://www2.jpl.nasa.gov/srtm/ Bilinear interpolation to the mesh vertices.

mesh/Coquimbo_20m/elem2D.out
Mesh geometry: List of the vertices forming the triangular elements.

Input Parameters

namelist.tsunami
Parameters read by TsunAWI at run time, e.g.,

path to the mesh and the output

the time step dt and the time frame of the simulation T_end

physical parameters like Manning friction, viscosity

toggle the output 

Additional command line parameters
The results below are for an idealized source (Cosine bell at the location of the 2015 Coquimbo event, started with 
ompTsuna.x -xy -71.674 -31.573 -mw 8.3 -id coquimbo_20m

Result

coquimbo_20m_new.nc 
Example output in netcdf format, raw data on the triangular mesh, see
ncdump -h coquimbo_20m.nc    

coquimbo_20m_new_inundation.nc
The same raw data but restricted to the inundated area, for faster postprocessing.

The TsunAWI version used is 
https://gitlab.awi.de/tsunawi/tsunawi/-/tags/LEXIS
with one minor bugfix: Add 
    use mpi_parallel
to line 881 of tsunawi/src/parameters.F90    

Compiled with tsunawi/src/Makefile.in_karolina
and run on karolina.it4i.cz in MPI-parallel mode with 128 tasks.

This work was supported by the LEXIS project funded by the EU’s Horizon 2020 research and innovation programme (2014-2020) under grant agreement No 825532.