Tsunami simulation - Revision history

Yjzhang at 17:20, 10 May 2017

2017-05-10T17:20:00Z

Yjzhang at 17:35, 18 December 2014

2014-12-18T17:35:03Z

Jan at 11:18, 24 July 2013

2013-07-24T11:18:40Z

Jan at 08:33, 24 July 2013

2013-07-24T08:33:33Z

Yjzhang at 18:07, 23 July 2013

2013-07-23T18:07:22Z

Yjzhang at 18:06, 23 July 2013

2013-07-23T18:06:49Z

Yjzhang at 18:05, 23 July 2013

2013-07-23T18:05:38Z

Yjzhang at 17:59, 23 July 2013

2013-07-23T17:59:24Z

Yjzhang at 17:58, 23 July 2013

2013-07-23T17:58:13Z

Yjzhang: Created page with "A sample tsunami run (impact of Alaska tsunami waves on Cannon Beach, OR) can be downloaded [http://www.ccrm.vims.edu/yinglong/wiki_files/ here]. Typically you need to fol..."

2013-07-23T17:57:55Z

Created page with "A sample tsunami run (impact of Alaska tsunami waves on Cannon Beach, OR) can be downloaded [http://www.ccrm.vims.edu/yinglong/wiki_files/ here]. Typically you need to fol..."

New page

A sample tsunami run (impact of Alaska tsunami waves on Cannon Beach, OR) can be downloaded
[http://www.ccrm.vims.edu/yinglong/wiki_files/ here].

Typically you need to follow these 2 steps in tsunami simulations:

* Deformation run (EX15/Def/ in the sample run): this simulates the earthquake and the set-up of the initial surface
waves. For this you need bdef.gr3 (refer to the user manual) which
specifies the total seafloor deformation. We typically run this stage for 10 sec duration and
at the end of the run you'll find an output files called 10_hotstart which is then used as
hotstart.in for the next stage. In addition, use mod_depth.f, which takes info in bdef.gr3 and
hgrid.gr3 (pre-earthquake depths) to generate hgrid.new (post-earthquake depths). The latter
is used in the next stage. For completeness you need to attach the boundary condition
part (b.tmp) of hgrid.gr3 to the end of hgrid.new;
* Propagation and inundation stage (EX15/ in the sample run): this run continues from the deformation
run above with no further seafloor movement. You'll notice that many input files are identical
to the Def/ run, but be careful of differences in hgrid.gr3 (linked to Def/hgrid.new) and param.in
(imm=0 etc). For this run we also turn off the advection in the open ocean (up to say 20m isobath) because it's negligible
there (adv.gr3). Also note that hotstart.in is from the previous stage. A relatively small
time step of 1 sec (by SELFE standard) is used because of the constraint from inundation stage.
The operational time step for tsunami applications is generally in the range of a few seconds.
The inundation results may also be sensitive to the min. depth used in the run (1cm in this example).
Also since only 2 S layers are used in the vertical (quasi-2D mode) we set the bottom friction to 0.
If you have to use non-zero friction you'll need to use more layers and make sure the bottom boundary
layer is resolved.
* After the run is done you can look at global outputs (elevation, depth-averaged velocity etc).
The maximum elevation (maxelev.gr3) and depth-averaged velocity (maxdahv.gr3) are also part of the
outputs. The maximum inundation can be computed from maxelev.gr3.

@@ Line 14: / Line 14: @@
        to the Def/ run, but be careful of differences in hgrid.gr3 (linked to Def/hgrid.new) and param.in
        (imm=0 etc).
-*For this run we also turn off the advection in the open ocean (up to say 20m isobath) because it's negligible there (adv.gr3) - see for info on CFL number.
+*For this run we also turn off the advection in the open ocean (up to say 60m isobath) because it's negligible there (adv.gr3) - see for info on CFL number.
 *Also note that hotstart.in is from the previous stage. A relatively small time step of 1 sec (by SCHISM standard) is used because of the constraint from inundation stage.
        The operational time step for tsunami applications is generally in the range of a few seconds.

@@ Line 1: / Line 1: @@
-A sample tsunami run (impact of Alaska tsunami waves on Cannon Beach, OR) can be downloaded [http://www.ccrm.vims.edu/yinglong/wiki_files/tsunami_ex15.tgz here]. Note that the files are compatible with the serial version of SELFE, but the idea for parallel version is similar.
+A sample tsunami run (impact of Alaska tsunami waves on Cannon Beach, OR) can be downloaded [http://www.ccrm.vims.edu/yinglong/wiki_files/tsunami_ex15.tgz here]. Note that the files are compatible with the serial version of SCHISM, but the idea for parallel version is similar.
 Typically you need to follow these 2 steps in tsunami simulations:
@@ Line 15: / Line 15: @@
        (imm=0 etc).
 *For this run we also turn off the advection in the open ocean (up to say 20m isobath) because it's negligible there (adv.gr3) - see for info on CFL number.
-*Also note that hotstart.in is from the previous stage. A relatively small time step of 1 sec (by SELFE standard) is used because of the constraint from inundation stage.
+*Also note that hotstart.in is from the previous stage. A relatively small time step of 1 sec (by SCHISM standard) is used because of the constraint from inundation stage.
        The operational time step for tsunami applications is generally in the range of a few seconds.
        The inundation results may also be sensitive to the min. depth used in the run (1cm in this example).

@@ Line 4: / Line 4: @@
 <OL>
    <LI> Deformation run (EX15/Def/ in the sample run): this simulates the earthquake and the set-up of the initial surface
-        waves. For this you need bdef.gr3 (refer to the user manual) which
+        waves. For this you need:
-       specifies the total seafloor deformation. We typically run this stage for 10 sec duration and
+*bdef.gr3 (refer to the user manual) which specifies the total seafloor deformation.
-       at the end of the run you'll find an output files called 10_hotstart which is then used as
+*turn on hotout, dt=1s, 10 steps, rnday=1.158e-4
-       hotstart.in for the next (propagation and inundation) stage. In addition, use mod_depth.f (inside the bundle), which takes info in bdef.gr3 and
+*We typically run this stage for 10 sec duration and at the end of the run you'll find an output files called 10_hotstart which is then used as hotstart.in for the next (propagation and inundation) stage.
-       hgrid.gr3 (pre-earthquake depths) to generate hgrid.new (post-earthquake depths). The latter
+*In addition, use mod_depth.f (inside the bundle), which takes info in bdef.gr3 and hgrid.gr3 (pre-earthquake depths) to generate hgrid.new (post-earthquake depths). The latter is used in the next stage. For completeness you need to attach the boundary condition part (b.tmp in the bundle) of hgrid.gr3 to the end of hgrid.new;
-       is used in the next stage. For completeness you need to attach the boundary condition
    <LI> Propagation and inundation stage (EX15/ in the sample run): this run continues from the deformation
        run above with no further seafloor movement. You'll notice that many input files are identical

@@ Line 21: / Line 21: @@
        Also since only 2 S layers are used in the vertical (quasi-2D mode) we set the bottom friction to 0.
       If you have to use non-zero friction you'll need to use more layers and make sure the bottom boundary
-      layer is resolved. If you use the recent parallel version, you can simply use a 2D model.
+      layer is resolved. If you use the recent parallel version, you can simply use a 2D model with a proper Manning coefficient.
    <LI> After the run is done you can look at global outputs (elevation, depth-averaged velocity etc).
         The maximum elevation (maxelev.gr3) and depth-averaged velocity (maxdahv.gr3) are also part of the
         outputs (for parallel versions, use trunk/src/Utility/Combining_Scripts/combine_gr3.f90). The maximum inundation can be computed from maxelev.gr3.
 </OL>

← Older revision		Revision as of 17:59, 23 July 2013
Line 1:		Line 1:
−	A sample tsunami run (impact of Alaska tsunami waves on Cannon Beach, OR) can be downloaded	+	A sample tsunami run (impact of Alaska tsunami waves on Cannon Beach, OR) can be downloaded [http://www.ccrm.vims.edu/yinglong/wiki_files/tsunami_ex15.tgz here].
−	[http://www.ccrm.vims.edu/yinglong/wiki_files/tsunami_ex15.tgz here].

	Typically you need to follow these 2 steps in tsunami simulations:		Typically you need to follow these 2 steps in tsunami simulations: