[cig-commits] r21373 - in seismo/3D/SPECFEM3D_GLOBE/trunk: DATA EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA EXAMPLES/global_PREM_kernels/amplitude/DATA EXAMPLES/global_PREM_kernels/traveltime/DATA EXAMPLES/global_s362ani/DATA EXAMPLES/noise_examples/global_long EXAMPLES/noise_examples/global_short EXAMPLES/noise_examples/regional EXAMPLES/noise_examples/test_global EXAMPLES/noise_examples/test_regional EXAMPLES/regional_Greece_small/DATA EXAMPLES/regional_MiddleEast/DATA EXAMPLES/small_benchmark_run_to_test_a_new_machine
elliott.sales.de.andrade at geodynamics.org
elliott.sales.de.andrade at geodynamics.org
Sat Feb 16 22:50:26 PST 2013
Author: elliott.sales.de.andrade
Date: 2013-02-16 22:50:25 -0800 (Sat, 16 Feb 2013)
New Revision: 21373
Modified:
seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_s362ani/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_Greece_small/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_MiddleEast/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/small_benchmark_run_to_test_a_new_machine/Par_file_this_example
Log:
Fix example Par_files.
I only added the missing parameters so that things don't crash because
of it. I have no idea whether the *other* parameters are still correct.
I also made the comments the same so that the diff would be smaller.
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file 2013-02-17 06:50:25 UTC (rev 21373)
@@ -36,7 +36,7 @@
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = s362ani
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 15.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
@@ -121,4 +121,3 @@
# output kernels on a regular grid instead of on the mesh points
SAVE_REGULAR_KL = .false.
-
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,6 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 0 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -14,8 +15,8 @@
GAMMA_ROTATION_AZIMUTH = 20.d0
# number of elements at the surface along the two sides of the first chunk
-# (must be multiple of 32 and 16 * multiple of NPROC below)
-NEX_XI = 512
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 512
NEX_ETA = 512
# number of MPI processors along the two sides of the first chunk
@@ -23,15 +24,21 @@
NPROC_ETA = 16
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362wmani, s362ani_prem, s29ea
-MODEL = 1D_transversely_isotropic_prem
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_transversely_isotropic_prem
# parameters describing the Earth model
OCEANS = .false.
@@ -48,8 +55,11 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
+MOVIE_COARSE = .true.
NTSTEP_BETWEEN_FRAMES = 100
HDUR_MOVIE = 0.d0
@@ -58,10 +68,8 @@
# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
-# type 4 saves the trace and deviatoric stress in the whole volume
-#MOVIE_VOLUME_COARSE saves movie only at corners of elements
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
MOVIE_VOLUME_TYPE = 2
-MOVIE_VOLUME_COARSE = .true.
MOVIE_TOP_KM = -100.0
MOVIE_BOTTOM_KM = 1000.0
MOVIE_WEST_DEG = -90.0
@@ -74,7 +82,7 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
@@ -85,7 +93,7 @@
NTSTEP_BETWEEN_OUTPUT_INFO = 100
# interval in time steps for temporary writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
NTSTEP_BETWEEN_READ_ADJSRC = 1000
# output format for the seismograms (one can use either or all of the three formats)
@@ -96,7 +104,7 @@
# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
ROTATE_SEISMOGRAMS_RT = .false.
-# decide wheter master process writes all the seismograms or if all processes do it in parallel
+# decide if master process writes all the seismograms or if all processes do it in parallel
WRITE_SEISMOGRAMS_BY_MASTER = .false.
# save all seismograms in one large combined file instead of one file per seismogram
@@ -110,3 +118,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,6 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 0 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -15,7 +16,7 @@
# number of elements at the surface along the two sides of the first chunk
# (must be multiple of 16 and 8 * multiple of NPROC below)
-NEX_XI = 256
+NEX_XI = 256
NEX_ETA = 256
# number of MPI processors along the two sides of the first chunk
@@ -23,15 +24,21 @@
NPROC_ETA = 8
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362wmani, s362ani_prem, s29ea
-MODEL = 1D_transversely_isotropic_prem
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_transversely_isotropic_prem
# parameters describing the Earth model
OCEANS = .false.
@@ -48,8 +55,11 @@
RECORD_LENGTH_IN_MINUTES = 100.0d0
# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
+MOVIE_COARSE = .true.
NTSTEP_BETWEEN_FRAMES = 100
HDUR_MOVIE = 0.d0
@@ -59,9 +69,7 @@
# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
-#MOVIE_VOLUME_COARSE saves movie only at corners of elements
MOVIE_VOLUME_TYPE = 2
-MOVIE_VOLUME_COARSE = .true.
MOVIE_TOP_KM = -100.0
MOVIE_BOTTOM_KM = 1000.0
MOVIE_WEST_DEG = -90.0
@@ -74,7 +82,7 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
@@ -85,7 +93,7 @@
NTSTEP_BETWEEN_OUTPUT_INFO = 100
# interval in time steps for temporary writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
NTSTEP_BETWEEN_READ_ADJSRC = 1000
# output format for the seismograms (one can use either or all of the three formats)
@@ -110,3 +118,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,6 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 0 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -15,7 +16,7 @@
# number of elements at the surface along the two sides of the first chunk
# (must be multiple of 16 and 8 * multiple of NPROC below)
-NEX_XI = 256
+NEX_XI = 256
NEX_ETA = 256
# number of MPI processors along the two sides of the first chunk
@@ -23,15 +24,21 @@
NPROC_ETA = 8
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362wmani, s362ani_prem, s29ea
-MODEL = 1D_isotropic_prem
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
# parameters describing the Earth model
OCEANS = .false.
@@ -48,8 +55,11 @@
RECORD_LENGTH_IN_MINUTES = 100.0d0
# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
+MOVIE_COARSE = .true.
NTSTEP_BETWEEN_FRAMES = 100
HDUR_MOVIE = 0.d0
@@ -59,9 +69,7 @@
# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
-#MOVIE_VOLUME_COARSE saves movie only at corners of elements
MOVIE_VOLUME_TYPE = 2
-MOVIE_VOLUME_COARSE = .true.
MOVIE_TOP_KM = -100.0
MOVIE_BOTTOM_KM = 1000.0
MOVIE_WEST_DEG = -90.0
@@ -74,7 +82,7 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
@@ -85,7 +93,7 @@
NTSTEP_BETWEEN_OUTPUT_INFO = 100
# interval in time steps for temporary writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
NTSTEP_BETWEEN_READ_ADJSRC = 1000
# output format for the seismograms (one can use either or all of the three formats)
@@ -110,3 +118,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,6 +1,6 @@
# forward or adjoint simulation
-SIMULATION_TYPE = 3
+SIMULATION_TYPE = 3
NOISE_TOMOGRAPHY = 0 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
@@ -24,10 +24,10 @@
NPROC_ETA = 8
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
@@ -36,7 +36,7 @@
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 30.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
@@ -82,7 +82,7 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
@@ -118,3 +118,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,6 +1,6 @@
# forward or adjoint simulation
-SIMULATION_TYPE = 3
+SIMULATION_TYPE = 3
NOISE_TOMOGRAPHY = 0 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
@@ -24,10 +24,10 @@
NPROC_ETA = 8
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
@@ -36,7 +36,7 @@
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 30.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
@@ -82,7 +82,7 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
@@ -118,3 +118,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_s362ani/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_s362ani/DATA/Par_file 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_s362ani/DATA/Par_file 2013-02-17 06:50:25 UTC (rev 21373)
@@ -24,10 +24,10 @@
NPROC_ETA = 5
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
@@ -36,7 +36,7 @@
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = s362ani
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 15.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
@@ -82,7 +82,7 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
@@ -118,3 +118,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 205.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 205.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 205.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 205.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_3_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/Par_file_NOISE_3_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 3
-NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 205.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_3_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/Par_file_NOISE_3_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 3
-NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_3_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/Par_file_NOISE_3_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 3
-NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 50.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 3
-NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 10
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1
-NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .true. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,7 +1,7 @@
# forward or adjoint simulation
SIMULATION_TYPE = 3
-NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
@@ -24,19 +24,19 @@
NPROC_ETA = 20
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_isotropic_prem
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 1.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
@@ -82,13 +82,12 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
LOCAL_PATH = /scratch/yangl/DATABASES_MPI
-#/scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -118,3 +117,7 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
+
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_Greece_small/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_Greece_small/DATA/Par_file 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_Greece_small/DATA/Par_file 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,11 +1,11 @@
# forward or adjoint simulation
-SIMULATION_TYPE = 1
+SIMULATION_TYPE = 1
NOISE_TOMOGRAPHY = 0 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
-NCHUNKS = 1
+NCHUNKS = 1
# angular width of the first chunk (not used if full sphere with six chunks)
ANGULAR_WIDTH_XI_IN_DEGREES = 20.d0 # angular size of a chunk
@@ -16,23 +16,28 @@
# number of elements at the surface along the two sides of the first chunk
# (must be multiple of 16 and 8 * multiple of NPROC below)
-NEX_XI = 64
+NEX_XI = 64
NEX_ETA = 64
# number of MPI processors along the two sides of the first chunk
-NPROC_XI = 2
+NPROC_XI = 2
NPROC_ETA = 2
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = 1D_transversely_isotropic_prem
# parameters describing the Earth model
@@ -50,6 +55,8 @@
RECORD_LENGTH_IN_MINUTES = 2.5d0
# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
MOVIE_COARSE = .true.
@@ -62,7 +69,6 @@
# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
-#MOVIE_VOLUME_COARSE saves movie only at corners of elements
MOVIE_VOLUME_TYPE = 2
MOVIE_TOP_KM = -100.0
MOVIE_BOTTOM_KM = 1000.0
@@ -76,7 +82,7 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .true.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
@@ -112,3 +118,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_MiddleEast/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_MiddleEast/DATA/Par_file 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_MiddleEast/DATA/Par_file 2013-02-17 06:50:25 UTC (rev 21373)
@@ -1,11 +1,11 @@
# forward or adjoint simulation
-SIMULATION_TYPE = 1
+SIMULATION_TYPE = 1
NOISE_TOMOGRAPHY = 0 # flag of noise tomography, three steps (1,2,3). If earthquake simulation, set it to 0.
SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
-NCHUNKS = 1
+NCHUNKS = 1
# angular width of the first chunk (not used if full sphere with six chunks)
ANGULAR_WIDTH_XI_IN_DEGREES = 45.d0 # angular size of a chunk
@@ -16,23 +16,28 @@
# number of elements at the surface along the two sides of the first chunk
# (must be multiple of 16 and 8 * multiple of NPROC below)
-NEX_XI = 128
+NEX_XI = 128
NEX_ETA = 128
# number of MPI processors along the two sides of the first chunk
-NPROC_XI = 8
+NPROC_XI = 8
NPROC_ETA = 8
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s20rts, s40rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = s29ea
# parameters describing the Earth model
@@ -50,6 +55,8 @@
RECORD_LENGTH_IN_MINUTES = 15.0d0
# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
MOVIE_COARSE = .true.
@@ -62,7 +69,6 @@
# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
-#MOVIE_VOLUME_COARSE saves movie only at corners of elements
MOVIE_VOLUME_TYPE = 2
MOVIE_TOP_KM = -100.0
MOVIE_BOTTOM_KM = 1000.0
@@ -76,7 +82,7 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
@@ -112,3 +118,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/small_benchmark_run_to_test_a_new_machine/Par_file_this_example
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/small_benchmark_run_to_test_a_new_machine/Par_file_this_example 2013-02-15 23:02:49 UTC (rev 21372)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/small_benchmark_run_to_test_a_new_machine/Par_file_this_example 2013-02-17 06:50:25 UTC (rev 21373)
@@ -24,10 +24,10 @@
NPROC_ETA = 2
# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135f_no_mud, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
#
# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135f_no_mud_onecrust
#
# fully 3D models:
# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
@@ -36,7 +36,7 @@
#
# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
+# associated reference model rather than the default 3D crustal model
# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
MODEL = s362ani
@@ -55,7 +55,7 @@
RECORD_LENGTH_IN_MINUTES = 55.0d0
# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
#MOVIE_COARSE does not work with create_movie_AVS_DX
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
@@ -82,7 +82,7 @@
# save mesh files to check the mesh
SAVE_MESH_FILES = .false.
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+# restart files (number of runs can be 1 or higher, choose 1 for no restart files)
NUMBER_OF_RUNS = 1
NUMBER_OF_THIS_RUN = 1
@@ -118,3 +118,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
More information about the CIG-COMMITS
mailing list