@@ -27,10 +27,10 @@ To fully define this case we **don't need all** of the possible input files, onl
2.`sw.input`
3.`extbc.input`
## First setup
# 3. First setup
Let's first setup `parameters.input` and `sw.input`, which are mandatory. Remember to put these files inside the `input` directory.
### Simulation control `parameters.input`.
## Simulation control `parameters.input`.
The simulation duration should be 20s, with a CFL value of 0.5, spatial output with a frequency of 1s and observation output frequency of 0.25s.
...
...
@@ -38,16 +38,16 @@ We will run this test on a single compute node, therefore ***no domain decomposi
We will use `NETCDF` output, default boundaries should be `REFLECTIVE`, and a reasonable screen output is between 200-1000.
### Shallow water model control `sw.input`
## Shallow water model control `sw.input`
The initial conditions for the experiment are those of a quiescent (i.e., zero velocity) water surface at sea level. In the typical shallow water variable naming `h` is water depth, `z` is the elevation of the bed relative to some arbitrary datum, and therefore `h+z` is the elevation of the water surface relative to the same datum. This datum must be consistent with that of the DEM. In the DEM `z=0` is the mean sea level. With this information, set `initialMode` and `initialValue`.
This benchmark typically uses Manning's friction law, with a roughness value of 0.01. Set this up in the file. Finally, because of the small scale of the experiment, and because in a tsunami we are particularly interested on the advancement of the wave over dry land, we will use a dry depth tolerance of 0.0001 m.
### Running the simulation
# 4. Running the simulation
Copy the example sbatch script that was created in a previous tutorial and modify it to include the path to this tsunami case. Then run the sbatch script.
### Visualising results
# 5. Visualising results
Let's check out the results. Although the output of this exercise is rather small, we will nonetheless visualise it in a lightweight way in the remote system.
We first start by forwarding graphics, often referred to as forwarding X from the remote server into our local machine, with the same SSH command with which you login, but including `-X`.
...
...
@@ -83,7 +83,7 @@ This is because
Nonetheless, it is important to highlight that still water equilibrium is a very important property of the numerical scheme signalling well-balancing between the fluxes and source terms in the numerical solver.
### Run the full simulation
# 6. Run the full simulation
Now we want to run the full simulation, with boundaries which model the incoming tsunami wave. For this we need to create the `extbc.input` file (`extbc` stands for external boundary conditions).
Details on how to setup boundary conditions are in the [SERGHEI BC wiki page](https://gitlab.com/serghei-model/serghei/-/wikis/User-Guide/Boundary-conditions).
