Thanks to visit codestin.com
Credit goes to github.com

Skip to content

Commit 96724f5

Browse files
author
Daniel Peter
committed
customized shakemovie version v5.1.1-shakemovie to run synthetics for website http://global.shakemovie.princeton.edu 
1 parent 6280047 commit 96724f5

File tree

180 files changed

+1873
-775682
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

180 files changed

+1873
-775682
lines changed

‎DATA/CMTSOLUTION.default

Lines changed: 0 additions & 13 deletions
This file was deleted.

‎DATA/Par_file.in

Lines changed: 19 additions & 19 deletions
Original file line numberDiff line numberDiff line change
@@ -15,13 +15,13 @@ CENTER_LONGITUDE_IN_DEGREES = 10.d0
1515
GAMMA_ROTATION_AZIMUTH = 20.d0
1616

1717
# number of elements at the surface along the two sides of the first chunk
18-
# (must be multiple of 16 and 8 * multiple of NPROC below)
19-
NEX_XI = 240
20-
NEX_ETA = 240
18+
# (must be multiple of 32 and 16 * multiple of NPROC below)
19+
NEX_XI = 256
20+
NEX_ETA = 256
2121

2222
# number of MPI processors along the two sides of the first chunk
23-
NPROC_XI = 5
24-
NPROC_ETA = 5
23+
NPROC_XI = 8
24+
NPROC_ETA = 8
2525

2626
# 1D models with real structure:
2727
# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
@@ -36,7 +36,7 @@ NPROC_ETA = 5
3636
#
3737
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
3838
# to take the 1D crustal model from the
39-
# associated reference model rather than the default 3D crustal model
39+
# associated reference model rather than the default 3D crustal model
4040
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
4141
MODEL = s362ani
4242

@@ -52,15 +52,15 @@ ATTENUATION = .true.
5252
ABSORBING_CONDITIONS = .false.
5353

5454
# record length in minutes
55-
RECORD_LENGTH_IN_MINUTES = 15.0d0
55+
RECORD_LENGTH_IN_MINUTES = 100.0d0
5656

5757
# save AVS or OpenDX movies
58-
#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
58+
#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
5959
#MOVIE_COARSE does not work with create_movie_AVS_DX
6060
MOVIE_SURFACE = .false.
6161
MOVIE_VOLUME = .false.
62-
MOVIE_COARSE = .false.
63-
NTSTEP_BETWEEN_FRAMES = 100
62+
MOVIE_COARSE = .true.
63+
NTSTEP_BETWEEN_FRAMES = 500
6464
HDUR_MOVIE = 0.d0
6565

6666
# save movie in volume. Will save element if center of element is in prescribed volume
@@ -69,13 +69,13 @@ HDUR_MOVIE = 0.d0
6969
# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
7070
# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
7171
# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
72-
MOVIE_VOLUME_TYPE = 2
72+
MOVIE_VOLUME_TYPE = 6
7373
MOVIE_TOP_KM = -100.0
7474
MOVIE_BOTTOM_KM = 1000.0
75-
MOVIE_WEST_DEG = -90.0
76-
MOVIE_EAST_DEG = 90.0
75+
MOVIE_WEST_DEG = 0.0
76+
MOVIE_EAST_DEG = 180.0
7777
MOVIE_NORTH_DEG = 90.0
78-
MOVIE_SOUTH_DEG = -90.0
78+
MOVIE_SOUTH_DEG = -30.0
7979
MOVIE_START = 0
8080
MOVIE_STOP = 40000
8181

@@ -87,25 +87,25 @@ NUMBER_OF_RUNS = 1
8787
NUMBER_OF_THIS_RUN = 1
8888

8989
# path to store the local database files on each node
90-
LOCAL_PATH = ./DATABASES_MPI
90+
LOCAL_PATH = SYN_DATABASES_MPI
9191

9292
# interval at which we output time step info and max of norm of displacement
9393
NTSTEP_BETWEEN_OUTPUT_INFO = 1000
9494

9595
# interval in time steps for temporary writing of seismograms
9696
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
97-
NTSTEP_BETWEEN_READ_ADJSRC = 1000
97+
NTSTEP_BETWEEN_READ_ADJSRC = 100000
9898

9999
# output format for the seismograms (one can use either or all of the three formats)
100-
OUTPUT_SEISMOS_ASCII_TEXT = .true.
100+
OUTPUT_SEISMOS_ASCII_TEXT = .false.
101101
OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
102-
OUTPUT_SEISMOS_SAC_BINARY = .false.
102+
OUTPUT_SEISMOS_SAC_BINARY = .true.
103103

104104
# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
105105
ROTATE_SEISMOGRAMS_RT = .false.
106106

107107
# decide if master process writes all the seismograms or if all processes do it in parallel
108-
WRITE_SEISMOGRAMS_BY_MASTER = .true.
108+
WRITE_SEISMOGRAMS_BY_MASTER = .false.
109109

110110
# save all seismograms in one large combined file instead of one file per seismogram
111111
# to avoid overloading shared non-local file systems such as GPFS for instance

‎DATA/Par_file.default renamed to ‎DATA/Par_file.movie

Lines changed: 19 additions & 19 deletions
Original file line numberDiff line numberDiff line change
@@ -15,13 +15,13 @@ CENTER_LONGITUDE_IN_DEGREES = 10.d0
1515
GAMMA_ROTATION_AZIMUTH = 20.d0
1616

1717
# number of elements at the surface along the two sides of the first chunk
18-
# (must be multiple of 16 and 8 * multiple of NPROC below)
19-
NEX_XI = 240
20-
NEX_ETA = 240
18+
# (must be multiple of 32 and 16 * multiple of NPROC below)
19+
NEX_XI = 128
20+
NEX_ETA = 128
2121

2222
# number of MPI processors along the two sides of the first chunk
23-
NPROC_XI = 5
24-
NPROC_ETA = 5
23+
NPROC_XI = 8
24+
NPROC_ETA = 8
2525

2626
# 1D models with real structure:
2727
# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
@@ -31,12 +31,12 @@ NPROC_ETA = 5
3131
#
3232
# fully 3D models:
3333
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
34-
# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
34+
# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
3535
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
3636
#
3737
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
3838
# to take the 1D crustal model from the
39-
# associated reference model rather than the default 3D crustal model
39+
# associated reference model rather than the default 3D crustal model
4040
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
4141
MODEL = s362ani
4242

@@ -52,12 +52,12 @@ ATTENUATION = .true.
5252
ABSORBING_CONDITIONS = .false.
5353

5454
# record length in minutes
55-
RECORD_LENGTH_IN_MINUTES = 15.0d0
55+
RECORD_LENGTH_IN_MINUTES = 100.0d0
5656

5757
# save AVS or OpenDX movies
58-
#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
58+
#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
5959
#MOVIE_COARSE does not work with create_movie_AVS_DX
60-
MOVIE_SURFACE = .false.
60+
MOVIE_SURFACE = .true.
6161
MOVIE_VOLUME = .false.
6262
MOVIE_COARSE = .false.
6363
NTSTEP_BETWEEN_FRAMES = 100
@@ -69,13 +69,13 @@ HDUR_MOVIE = 0.d0
6969
# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
7070
# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
7171
# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
72-
MOVIE_VOLUME_TYPE = 2
72+
MOVIE_VOLUME_TYPE = 6
7373
MOVIE_TOP_KM = -100.0
7474
MOVIE_BOTTOM_KM = 1000.0
75-
MOVIE_WEST_DEG = -90.0
76-
MOVIE_EAST_DEG = 90.0
75+
MOVIE_WEST_DEG = 0.0
76+
MOVIE_EAST_DEG = 180.0
7777
MOVIE_NORTH_DEG = 90.0
78-
MOVIE_SOUTH_DEG = -90.0
78+
MOVIE_SOUTH_DEG = -30.0
7979
MOVIE_START = 0
8080
MOVIE_STOP = 40000
8181

@@ -87,25 +87,25 @@ NUMBER_OF_RUNS = 1
8787
NUMBER_OF_THIS_RUN = 1
8888

8989
# path to store the local database files on each node
90-
LOCAL_PATH = ./DATABASES_MPI
90+
LOCAL_PATH = MOVIE_DATABASES_MPI
9191

9292
# interval at which we output time step info and max of norm of displacement
9393
NTSTEP_BETWEEN_OUTPUT_INFO = 1000
9494

9595
# interval in time steps for temporary writing of seismograms
9696
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
97-
NTSTEP_BETWEEN_READ_ADJSRC = 1000
97+
NTSTEP_BETWEEN_READ_ADJSRC = 100000
9898

9999
# output format for the seismograms (one can use either or all of the three formats)
100-
OUTPUT_SEISMOS_ASCII_TEXT = .true.
100+
OUTPUT_SEISMOS_ASCII_TEXT = .false.
101101
OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
102-
OUTPUT_SEISMOS_SAC_BINARY = .false.
102+
OUTPUT_SEISMOS_SAC_BINARY = .true.
103103

104104
# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
105105
ROTATE_SEISMOGRAMS_RT = .false.
106106

107107
# decide if master process writes all the seismograms or if all processes do it in parallel
108-
WRITE_SEISMOGRAMS_BY_MASTER = .true.
108+
WRITE_SEISMOGRAMS_BY_MASTER = .false.
109109

110110
# save all seismograms in one large combined file instead of one file per seismogram
111111
# to avoid overloading shared non-local file systems such as GPFS for instance

‎DATA/Par_file.syn

Lines changed: 120 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,120 @@
1+
2+
# forward or adjoint simulation
3+
SIMULATION_TYPE = 1
4+
NOISE_TOMOGRAPHY = 0 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
5+
SAVE_FORWARD = .false. # save last frame of forward simulation or not
6+
7+
# number of chunks (1,2,3 or 6)
8+
NCHUNKS = 6
9+
10+
# angular width of the first chunk (not used if full sphere with six chunks)
11+
ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
12+
ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
13+
CENTER_LATITUDE_IN_DEGREES = 40.d0
14+
CENTER_LONGITUDE_IN_DEGREES = 10.d0
15+
GAMMA_ROTATION_AZIMUTH = 20.d0
16+
17+
# number of elements at the surface along the two sides of the first chunk
18+
# (must be multiple of 32 and 16 * multiple of NPROC below)
19+
NEX_XI = 256
20+
NEX_ETA = 256
21+
22+
# number of MPI processors along the two sides of the first chunk
23+
NPROC_XI = 8
24+
NPROC_ETA = 8
25+
26+
# 1D models with real structure:
27+
# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
28+
#
29+
# 1D models with only one fictitious averaged crustal layer:
30+
# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
31+
#
32+
# fully 3D models:
33+
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
34+
# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
35+
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
36+
#
37+
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
38+
# to take the 1D crustal model from the
39+
# associated reference model rather than the default 3D crustal model
40+
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
41+
MODEL = s362ani
42+
43+
# parameters describing the Earth model
44+
OCEANS = .true.
45+
ELLIPTICITY = .true.
46+
TOPOGRAPHY = .true.
47+
GRAVITY = .true.
48+
ROTATION = .true.
49+
ATTENUATION = .true.
50+
51+
# absorbing boundary conditions for a regional simulation
52+
ABSORBING_CONDITIONS = .false.
53+
54+
# record length in minutes
55+
RECORD_LENGTH_IN_MINUTES = 100.0d0
56+
57+
# save AVS or OpenDX movies
58+
#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
59+
#MOVIE_COARSE does not work with create_movie_AVS_DX
60+
MOVIE_SURFACE = .false.
61+
MOVIE_VOLUME = .false.
62+
MOVIE_COARSE = .true.
63+
NTSTEP_BETWEEN_FRAMES = 500
64+
HDUR_MOVIE = 0.d0
65+
66+
# save movie in volume. Will save element if center of element is in prescribed volume
67+
# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
68+
# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
69+
# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
70+
# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
71+
# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
72+
MOVIE_VOLUME_TYPE = 6
73+
MOVIE_TOP_KM = -100.0
74+
MOVIE_BOTTOM_KM = 1000.0
75+
MOVIE_WEST_DEG = 0.0
76+
MOVIE_EAST_DEG = 180.0
77+
MOVIE_NORTH_DEG = 90.0
78+
MOVIE_SOUTH_DEG = -30.0
79+
MOVIE_START = 0
80+
MOVIE_STOP = 40000
81+
82+
# save mesh files to check the mesh
83+
SAVE_MESH_FILES = .false.
84+
85+
# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
86+
NUMBER_OF_RUNS = 1
87+
NUMBER_OF_THIS_RUN = 1
88+
89+
# path to store the local database files on each node
90+
LOCAL_PATH = SYN_DATABASES_MPI
91+
92+
# interval at which we output time step info and max of norm of displacement
93+
NTSTEP_BETWEEN_OUTPUT_INFO = 1000
94+
95+
# interval in time steps for temporary writing of seismograms
96+
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
97+
NTSTEP_BETWEEN_READ_ADJSRC = 100000
98+
99+
# output format for the seismograms (one can use either or all of the three formats)
100+
OUTPUT_SEISMOS_ASCII_TEXT = .false.
101+
OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
102+
OUTPUT_SEISMOS_SAC_BINARY = .true.
103+
104+
# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
105+
ROTATE_SEISMOGRAMS_RT = .false.
106+
107+
# decide if master process writes all the seismograms or if all processes do it in parallel
108+
WRITE_SEISMOGRAMS_BY_MASTER = .false.
109+
110+
# save all seismograms in one large combined file instead of one file per seismogram
111+
# to avoid overloading shared non-local file systems such as GPFS for instance
112+
SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
113+
USE_BINARY_FOR_LARGE_FILE = .false.
114+
115+
# flag to impose receivers at the surface or allow them to be buried
116+
RECEIVERS_CAN_BE_BURIED = .true.
117+
118+
# print source time function
119+
PRINT_SOURCE_TIME_FUNCTION = .false.
120+

0 commit comments

Comments
 (0)