How much distance that tsunami wave can travel from the gauge point?

[Teja sri]

Hello,

I am working on tsunami inundation hazard modeling using Clawpack. I understand how to obtain tsunami wave heights at the shoreline from the Chile 2010 example, but I am struggling to determine the inundation values.

Could you please help me with the following?

1. How can I determine how much dry land is flooded due to the tsunami from  the gauge point?

2. How do I read fort.q00 files?

3. How do I read gauge files from the _output directory?

4. How to i get tsunami velocities?

5. What are the parameters that have to be included extra in setrun.py to get the flooded length from gauge point?

Thank you in advance!

[Kyle Mandli]

I am not entirely sure what you are asking for in some of these, but here is some guidance.  For all of these questions, the Python code provided [1, 2] will be able to read in the output files and interpret them.  The plotting code uses the Python code of interest and would be a place to look.  For instance,

from clawpack.pyclaw.solution import Solution

sol = Solution(0, path="_output")

print(sol)

sol.states[0].q[0, :, :]

reads in the first time output (probably the initial condition) solution and accesses the first grid available's depth field (note that there is a lot of complexity here because of the adaptive grids).  Gauges have a similar mechanism with

from clawpack.pyclaw.gauges import GaugeSolution

gauge = GaugeSolution(1, path="_output")

print(gauge)

gauge.q[0, :]

which will read in gauge 1 and access the q[0] column of output (the depth h).

The interpretation of the columns of data come from the q vector, which in GeoClaw is usually h (depth), hu (x-momentum), hv (y-momentum), and eta (surface).  The bathymetry (b), depth and surface are related by

h = eta - b

so with two of them you can figure out everything else.  The variable eta is always defined this way, so if h = 0 (there is no water), eta = b.

There is a lot more documentation and pointers to papers at [3].  Also take a look at the Jupyter notebooks provided [4] as they contain a lot of useful workflows for working with input and output data in GeoClaw.  There are also some streaming videos on [5] that may be helpful.  Feel free to ask additional follow-up questions, but I would recommend detailing what you have tried and what is going wrong so we can better help you.  Also, someone else may have some pointers about specific aspects of your questions.

1. https://www.clawpack.org/pyclaw/solution.html
2. http://www.clawpack.org/gauges.html
3. http://www.clawpack.org/
4. http://www.clawpack.org/gallery/notebooks.html
5. http://www.clawpack.org/community.html

Kyle

--
You received this message because you are subscribed to the Google Groups "claw-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to claw-users+...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/claw-users/311ac9a0-6e44-45d3-ba1a-a73962cafd75n%40googlegroups.com.

[Randall J LeVeque]

Kyle sent some good references for getting started, the Jupyter notebooks on GeoClaw and Tsunami modeling in particular may be useful, and you might also look at some of the webinars and tutorials on the Clawpack Community page.

I should note that the chile2010 example in GeoClaw is not set up to model coastal inundation properly -- the grids are far too coarse for that, and finer coastal topography DEMs would also need to be included as input for inundation runs.  Some of the other examples and notebooks give a better idea of how to set things up for coastal modeling, including using gauges and also fgmax grids to capture maximum inundation and flow speeds.

 - Randy

To view this discussion visit https://groups.google.com/d/msgid/claw-users/ba1e7f7b-2904-4dbf-bf26-db4fc5e53d52%40Spark.

[Teja sri]

Thankyou for the references.

To view this discussion visit https://groups.google.com/d/msgid/claw-users/CADwYU3%2B_yrOtznKr-fPhxTF14MnL3ieujL8EZuR%3DFOOXdndksA%40mail.gmail.com.