Skip to content

Commit 77065ed

Browse files
committed
updates regional_Berkeley/ example (adds Berkeley source time function to Par_file; adds reference solution)
1 parent 5541742 commit 77065ed

22 files changed

+43947
-17
lines changed

EXAMPLES/regional_Berkeley/DATA/Par_file

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -269,6 +269,16 @@ USE_MONOCHROMATIC_CMT_SOURCE = .false.
269269
# print source time function
270270
PRINT_SOURCE_TIME_FUNCTION = .true.
271271

272+
## Berkeley source time function
273+
STF_IS_UCB_HEAVISIDE = .true.
274+
# UCB source frequency content (i.e., heaviside function)
275+
SOURCE_T1 = 500.d0
276+
SOURCE_T2 = 400.d0
277+
SOURCE_T3 = 50.d0
278+
SOURCE_T4 = 40.d0
279+
# UCB source time shift
280+
TAU = 500.d0
281+
272282
#-----------------------------------------------------------
273283
#
274284
# Seismograms

EXAMPLES/regional_Berkeley/README

Lines changed: 15 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -2,9 +2,22 @@
22
README
33
----------------------------------------------------------
44

5-
This example creates a regional mesh using transversely isotropic PREM
5+
This example creates a regional mesh using the 3D Berkeley model SEMUCB_A3d
66
and runs a spectral-element simulation.
77

8+
In the Par_file, we also add the following parameters to use a specific Berkeley source time function:
9+
```
10+
## Berkeley source time function
11+
STF_IS_UCB_HEAVISIDE = .true.
12+
# UCB source frequency content (i.e., heaviside function)
13+
SOURCE_T1 = 500.d0
14+
SOURCE_T2 = 400.d0
15+
SOURCE_T3 = 50.d0
16+
SOURCE_T4 = 40.d0
17+
# UCB source time shift
18+
TAU = 500.d0
19+
```
20+
821
The simulation is setup to use 4 CPUs only, thus can be run locally on a desktop machine
922
(needs to have MPI installed).
1023

@@ -18,24 +31,9 @@ type:
1831
to setup this directory for the simulation and to run a corresponding
1932
script "run_mesher_solver.bash".
2033

21-
The simulation is accurate down to a shortest period of ~ 15 s.
34+
The simulation is accurate down to a shortest period of ~ 21 s.
2235
(minimum period resolved ~ (256/NEX_XI) * (ANGULAR_WIDTH_XI_IN_DEGREES/90) * 17)
2336

2437
Reference seismograms for comparisons are provided in the directory: REF_SEIS/
25-
26-
27-
2. kernel simulation:
28-
---------------------
29-
30-
type:
31-
./run_this_example.kernel.sh
32-
33-
to setup this directory for the simulation and to submit a corresponding
34-
script "run_mesher_solver.kernel.bash".
35-
36-
The simulation will use the adjoint source given in directory: SEM/
37-
to compute a travel-time sensitivity kernel.
38-
39-
4038
Please modify these scripts according to your specifics.
4139

0 commit comments

Comments
 (0)