[cig-commits] r22461 - in seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES: . benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA global_PREM_kernels/amplitude global_PREM_kernels/amplitude/DATA global_PREM_kernels/traveltime global_PREM_kernels/traveltime/DATA global_s362ani global_s362ani/DATA noise_examples/global_long noise_examples/global_short noise_examples/regional noise_examples/test_global noise_examples/test_regional regional_Greece_small regional_Greece_small/DATA regional_MiddleEast regional_MiddleEast/DATA
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Sun Jun 30 07:28:52 PDT 2013
Author: dkomati1
Date: 2013-06-30 07:28:52 -0700 (Sun, 30 Jun 2013)
New Revision: 22461
Added:
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/small_benchmark_run_to_test_more_complex_Earth/
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/small_benchmark_run_to_test_very_simple_Earth/
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/vi_all_Par_files_to_add_parameters_to_them.bash
Removed:
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/small_benchmark_run_to_test_a_new_machine/
Modified:
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/process.kernel.sh
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/process.sh
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/process.kernel.sh
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/process.sh
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/process.kernel.sh
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/process.sh
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/process.kernel.sh
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/process.sh
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/process.kernel.sh
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/process.sh
Log:
updated "EXAMPLES" from the trunk; will need to re-add some GPU parameters to all the Par_files later
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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.
@@ -47,9 +54,19 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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 +75,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 +89,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 +100,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 +111,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 +125,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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.
@@ -47,9 +54,19 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 100.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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 +76,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 +89,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 +100,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 +125,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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.
@@ -47,9 +54,19 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 100.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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 +76,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 +89,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 +100,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 +125,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .true.
ROTATION = .true.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 30.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 1000
@@ -122,5 +125,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-# set to true to use GPUs
-GPU_MODE = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/process.kernel.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/process.kernel.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/process.kernel.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -37,7 +37,7 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# compiles for an adjoint simulation
cp $currentdir/DATA/Par_file DATA/Par_file
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/process.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/process.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/amplitude/process.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -33,7 +33,7 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# compiles for a forward simulation
cp $currentdir/DATA/Par_file DATA/Par_file
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file 2013-06-30 14:28:52 UTC (rev 22461)
@@ -1,6 +1,6 @@
# forward or adjoint simulation
-SIMULATION_TYPE = 1
+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
@@ -47,7 +47,6 @@
GRAVITY = .true.
ROTATION = .true.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 30.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 1000
@@ -122,5 +125,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-# set to true to use GPUs
-GPU_MODE = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/process.kernel.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/process.kernel.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/process.kernel.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -37,7 +37,7 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# compiles for an adjoint simulation
cp $currentdir/DATA/Par_file DATA/Par_file
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/process.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/process.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_PREM_kernels/traveltime/process.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -33,7 +33,7 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# compiles for a forward simulation
cp $currentdir/DATA/Par_file DATA/Par_file
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/DATA/Par_file 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/DATA/Par_file 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .true.
ROTATION = .true.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 15.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 1000
@@ -122,5 +125,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-# set to true to use GPUs
-GPU_MODE = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/process.kernel.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/process.kernel.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/process.kernel.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -37,7 +37,7 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# compiles for an adjoint simulation
cp $currentdir/DATA/Par_file DATA/Par_file
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/process.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/process.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/global_s362ani/process.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -33,7 +33,7 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# compiles for a forward simulation
cp $currentdir/DATA/Par_file DATA/Par_file
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 205.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_1_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 205.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 205.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_2_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 205.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_3_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_long/Par_file_NOISE_3_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 205.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_1_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_2_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_3_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/global_short/Par_file_NOISE_3_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_1_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_2_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_3_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/regional/Par_file_NOISE_3_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 50.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -47,7 +47,6 @@
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -55,8 +54,15 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 1.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -83,15 +89,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
-# mesh files
-LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 10000
@@ -121,3 +124,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/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/DATA/Par_file 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/DATA/Par_file 2013-06-30 14:28:52 UTC (rev 22461)
@@ -1,39 +1,44 @@
# 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 = 90.d0 # angular size of a chunk
-ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+ANGULAR_WIDTH_XI_IN_DEGREES = 20.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 20.d0
CENTER_LATITUDE_IN_DEGREES = 40.d0
CENTER_LONGITUDE_IN_DEGREES = 25.d0
GAMMA_ROTATION_AZIMUTH = 0.d0
# 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 = 48
-NEX_ETA = 48
+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
-MODEL = 1D_transversely_isotropic_prem
+#
+# 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 = .true.
@@ -42,15 +47,23 @@
GRAVITY = .true.
ROTATION = .true.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# record length in minutes
-RECORD_LENGTH_IN_MINUTES = 5.1d0
+RECORD_LENGTH_IN_MINUTES = 2.5d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -63,7 +76,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_COARSE saves movie only at corners of elements
MOVIE_VOLUME_TYPE = 2
MOVIE_TOP_KM = -100.0
MOVIE_BOTTOM_KM = 1000.0
@@ -77,18 +89,15 @@
# 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
# path to store the local database files on each node
-# mesh files
LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO = 500
+NTSTEP_BETWEEN_OUTPUT_INFO = 50
# interval in time steps for temporary writing of seismograms
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
@@ -116,13 +125,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-# set to true to use GPUs
-GPU_MODE = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
-# set to true to use ADIOS
-ADIOS_ENABLED = .true.
-ADIOS_FOR_FORWARD_ARRAYS = .true.
-ADIOS_FOR_MPI_ARRAYS = .true.
-ADIOS_FOR_ARRAYS_SOLVER = .true.
-ADIOS_FOR_SOLVER_MESHFILES = .true.
-ADIOS_FOR_AVS_DX = .true.
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/process.kernel.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/process.kernel.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/process.kernel.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -36,9 +36,9 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-#./configure F90=ifort MPIF90=/usr/local/openmpi-ifort/bin/mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=/usr/local/openmpi-ifort/bin/mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# configures package with gfortran compiler
-#./configure F90=gfortran MPIF90=mpif90 FLAGS_CHECK="-O3" FLAGS_NO_CHECK="-O3"
+#./configure F90=gfortran MPIF90=mpif90 FLAGS_CHECK="-O3"
# compiles for an adjoint simulation
cp $currentdir/DATA/Par_file DATA/Par_file
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/process.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/process.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_Greece_small/process.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -32,15 +32,15 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-#./configure F90=ifort MPIF90=/usr/local/openmpi-ifort/bin/mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=/usr/local/openmpi-ifort/bin/mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# configures package with gfortran compiler
-#./configure F90=gfortran MPIF90=mpif90 FLAGS_CHECK="-O3" FLAGS_NO_CHECK="-O3"
+#./configure F90=gfortran MPIF90=mpif90 FLAGS_CHECK="-O3"
# compiles for a forward simulation
cp $currentdir/DATA/Par_file DATA/Par_file
make clean
-make #>& $currentdir/tmp_make_output.log
-make xcombine_vol_data #>> $currentdir/tmp_make_output.log
+make >& $currentdir/tmp_make_output.log
+make xcombine_vol_data >> $currentdir/tmp_make_output.log
# backup of constants setup
cp setup/* $currentdir/OUTPUT_FILES/
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/DATA/Par_file 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/DATA/Par_file 2013-06-30 14:28:52 UTC (rev 22461)
@@ -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
@@ -42,7 +47,6 @@
GRAVITY = .true.
ROTATION = .true.
ATTENUATION = .true.
-ATTENUATION_NEW = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
@@ -50,7 +54,16 @@
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 15.0d0
+# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
+# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
+# but requires a significant amount of disk space for temporary storage.
+PARTIAL_PHYS_DISPERSION_ONLY = .true.
+UNDO_ATTENUATION = .false.
+NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
+
# 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.
@@ -63,7 +76,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_COARSE saves movie only at corners of elements
MOVIE_VOLUME_TYPE = 2
MOVIE_TOP_KM = -100.0
MOVIE_BOTTOM_KM = 1000.0
@@ -77,15 +89,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
-# mesh files
LOCAL_PATH = ./DATABASES_MPI
-# temporary wavefield/kernel/movie files
-LOCAL_TMP_PATH = ./DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 1000
@@ -116,5 +125,6 @@
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-# set to true to use GPUs
-GPU_MODE = .false.
+# output kernels on a regular grid instead of on the mesh points
+SAVE_REGULAR_KL = .false.
+
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/process.kernel.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/process.kernel.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/process.kernel.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -37,7 +37,7 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# compiles for an adjoint simulation
cp $currentdir/DATA/Par_file DATA/Par_file
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/process.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/process.sh 2013-06-30 14:11:14 UTC (rev 22460)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/regional_MiddleEast/process.sh 2013-06-30 14:28:52 UTC (rev 22461)
@@ -33,7 +33,7 @@
# using default configuration
cd ../../
# configures package with ifort compiler
-./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" FLAGS_NO_CHECK="-O3 -assume byterecl" > tmp.log
+./configure F90=ifort MPIF90=mpif90 FLAGS_CHECK="-O3 -assume byterecl" > tmp.log
# compiles for a forward simulation
cp $currentdir/DATA/Par_file DATA/Par_file
Copied: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/vi_all_Par_files_to_add_parameters_to_them.bash (from rev 22460, seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/vi_all_Par_files_to_add_parameters_to_them.bash)
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/vi_all_Par_files_to_add_parameters_to_them.bash (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/EXAMPLES/vi_all_Par_files_to_add_parameters_to_them.bash 2013-06-30 14:28:52 UTC (rev 22461)
@@ -0,0 +1,4 @@
+#!/bin/bash
+
+vi ../DATA/Par_file ./global_s362ani/DATA/Par_file ./regional_MiddleEast/DATA/Par_file ./small_benchmark_run_to_test_more_complex_Earth/Par_file_this_example ./noise_examples/test_regional/Par_file_NOISE_3_noattenuation ./noise_examples/test_regional/Par_file_NOISE_2_noattenuation ./noise_examples/test_regional/Par_file_NOISE_1_noattenuation ./noise_examples/test_regional/Par_file_NOISE_1_attenuation ./noise_examples/test_regional/Par_file_NOISE_2_attenuation ./noise_examples/test_global/Par_file_NOISE_3_noattenuation ./noise_examples/test_global/Par_file_NOISE_2_noattenuation ./noise_examples/test_global/Par_file_NOISE_1_noattenuation ./noise_examples/test_global/Par_file_NOISE_1_attenuation ./noise_examples/test_global/Par_file_NOISE_2_attenuation ./noise_examples/global_short/Par_file_NOISE_3_noattenuation ./noise_examples/global_short/Par_file_NOISE_2_noattenuation ./noise_examples/global_short/Par_file_NOISE_1_noattenuation ./noise_examples/global_short/Par_file_NOISE_1_attenuation ./noise_examples/global_short/Par_file_NOISE_2_attenuation ./noise_examples/regional/Par_file_NOISE_3_noattenuation ./noise_examples/regional/Par_file_NOISE_2_noattenuation ./noise_examples/regional/Par_file_NOISE_1_noattenuation ./noise_examples/regional/Par_file_NOISE_1_attenuation ./noise_examples/regional/Par_file_NOISE_2_attenuation ./noise_examples/global_long/Par_file_NOISE_3_noattenuation ./noise_examples/global_long/Par_file_NOISE_2_noattenuation ./noise_examples/global_long/Par_file_NOISE_1_noattenuation ./noise_examples/global_long/Par_file_NOISE_1_attenuation ./noise_examples/global_long/Par_file_NOISE_2_attenuation ./regional_Greece_small/DATA/Par_file ./small_benchmark_run_to_test_very_simple_Earth/Par_file_this_example ./global_PREM_kernels/amplitude/DATA/Par_file ./global_PREM_kernels/traveltime/DATA/Par_file ./benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file ./benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file ./benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
+
More information about the CIG-COMMITS
mailing list