You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: EXAMPLES/real_world/Myanmar_Mw7.7/README.md
+44-45Lines changed: 44 additions & 45 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -44,37 +44,35 @@ Furthermore, we want to use a (realistic) velocity model for this region.
44
44
We will take one of the EarthScope (IRIS) [EMC models](https://ds.iris.edu/ds/products/emc-earthmodels/).
45
45
For this region, a nice model seems to be [FWEA23](https://ds.iris.edu/ds/products/emc-fwea23/) by Liu et al. (2024).
46
46
47
-
48
-
49
47
1.**Setup model**:
50
48
51
-
We will first setup the topography surface and tomographic model for meshing our region.
52
-
In this example folder, we provide a bash script `setup_model.sh` to setup these model files.
53
-
54
-
Just run the setup script with the default target region:
55
-
```
56
-
> ./setup_model.sh
57
-
```
58
-
59
-
This might take a while to complete and download the rather large EMC model file.
60
-
Note that we added additional interfaces at 100km and 6km depth to facilitate the meshing with doubling layers.
61
-
62
-
Furthermore, USGS provides a global [Vs30 model and data](https://earthquake.usgs.gov/data/vs30/) set.
63
-
The script will download and extract a corresponding Vs30-interface for our region.
64
-
65
-
49
+
We will first setup the topography surface and tomographic model for meshing our region.
50
+
In this example folder, we provide a bash script `setup_model.sh` to setup these model files.
51
+
52
+
Just run the setup script with the default target region:
53
+
```
54
+
> ./setup_model.sh
55
+
```
56
+
57
+
This might take a while to complete and download the rather large EMC model file.
58
+
Note that we added additional interfaces at 100km and 6km depth to facilitate the meshing with doubling layers.
59
+
60
+
Furthermore, USGS provides a global [Vs30 model and data](https://earthquake.usgs.gov/data/vs30/) set.
61
+
The script will download and extract a corresponding Vs30-interface for our region.
62
+
63
+
66
64
2. **Wave simulation**:
67
65
68
-
After you have completed the mesh setup, you can run the in-house mesher `xmeshfem3D` and `xgenerate_databases` to create the spectral-element mesh of our region.
69
-
The seismic wave propagation solver `xspecfem3D` then creates the needed movie data files for our visualization.
70
-
71
-
To run the simulation, just type:
72
-
```
73
-
> ./run_this_example.sh
74
-
```
66
+
After you have completed the mesh setup, you can run the in-house mesher `xmeshfem3D` and `xgenerate_databases` to create the spectral-element mesh of our region.
67
+
The seismic wave propagation solver `xspecfem3D` then creates the needed movie data files for our visualization.
68
+
69
+
To run the simulation, just type:
70
+
```
71
+
> ./run_this_example.sh
72
+
```
75
73
76
74
77
-
3.**Visualization**:
75
+
3. **Visualization**:
78
76
79
77
For our simulation here, we turned on the surface movie and shake map outputs in `DATA/Par_file` like:
80
78
```
@@ -91,21 +89,21 @@ For this region, a nice model seems to be [FWEA23](https://ds.iris.edu/ds/produc
91
89
NTSTEP_BETWEEN_FRAMES = 100
92
90
HDUR_MOVIE = 0.0
93
91
```
94
-
92
+
95
93
To visualize the corresponding output data (`OUTPUT_FILES/moviedata***`), we can create movie snapshot files as `OUTPUT_FILES/AVS_*.inp` files:
96
94
```
97
95
> ./xcreate_movie_files.sh
98
96
```
99
-
97
+
100
98
Similar for the shakemap, we can plot the peak-ground velocity (PGV) values by
101
99
```
102
100
> ./xcreate_shakemap.sh 2
103
101
```
104
-
102
+
105
103
You can use for example [Paraview](https://www.paraview.org) to look at the created `OUTPUT_FILEs/AVS_movie*.inp` files.
106
104
107
105
108
-
That's it, try it out for yourself...
106
+
That's it, try it out for yourself...
109
107
110
108
111
109
## Visualization with Blender (optional)
@@ -117,14 +115,12 @@ A single wavefield snapshot like above can be plotted using [Blender](https://ww
117
115
> ./xcreate_snapshot.sh 10000
118
116
```
119
117
120
-
121
118
Regarding shakemap image coloring, for conversion to Modified Mercalli Intensity (MMI), Wald et al. suggest a table like:
151
147
152
148
149
+
153
150
## Spectral Acceleration (SA) - response spectrum at Bangkok Chatuchak building collapse
154
151
155
152
Our simulation outputs the ground motions at the Chatuchak building site in Bangkok, Thailand. In particular, we are interested in how the [Spectral Acceleration (SA)](https://en.wikipedia.org/wiki/Spectral_acceleration) for a building at this location looked like.
@@ -163,22 +160,24 @@ A simple script to produce this figure is provided in this example directory:
163
160
> ./xcreate_spectral_response.sh
164
161
```
165
162
163
+
166
164
## Sonification of traces (optional)
167
165
168
-
To hear how the ground shaking sounded, we can "sonify" the output traces, i.e., map the traces to an audible spectrum.
169
-
First, get the script from the [shakemovie](https://github.com/SPECFEM/shakemovie) repository:
For comparison, we provide a reference solution in folder `REF_SEIS/` with corresponding output files.
184
-
The simulation was run in parallel using 16 MPI processes, each using a single Nvidia A100 GPU. The total simulated time is 1000 s with a time-to-solution of ~ 24 min 15 s.
182
+
For comparison, we provide a reference solution in folder `REF_SEIS/` with corresponding output files.
183
+
The simulation was run in parallel using 16 MPI processes, each using a single Nvidia A100 GPU. The total simulated time is 1000 s with a time-to-solution of ~ 24 min 15 s.
0 commit comments