fix for acoustic-elastic PML coupling on GPUs; adds parameters PRINT_SOURCE_TIME_FUNCTION and SAVE_MESH_FILES to Par_file for turning on/off optional mesh outputs

DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -336,6 +339,9 @@ nbregions                       = 5              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 100
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/Rayleigh_wave_no_crack/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -321,6 +324,9 @@ nbregions                       = 1              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 200
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/Rayleigh_wave_no_crack_simpler_way_to_generate_the_wave/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -321,6 +324,9 @@ nbregions                       = 1              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 100
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/Rayleigh_wave_with_crack/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -327,6 +330,9 @@ nbregions                       = 3              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 200
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/anisotropy/anisotropic_isotropic_model/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -325,6 +328,9 @@ nbregions                       = 3              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 100
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/anisotropy/anisotropic_isotropic_model/DATA/Par_file_homogeneous_with_attenuation

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -330,6 +333,9 @@ nbregions                       = 3              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 100
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/anisotropy/anisotropic_zinc_crystal/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -321,6 +324,9 @@ nbregions                       = 1              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 100
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/axisymmetric_examples/axisymmetric_case_AXISYM_option/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -339,6 +342,9 @@ nbregions                       = 2              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 500
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/axisymmetric_examples/axisymmetric_case_AXISYM_option/REF_SEIS/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -339,6 +342,9 @@ nbregions                       = 2              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 500
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/axisymmetric_examples/global_Earth_ak135f/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -325,6 +328,9 @@ nbregions                       = 2              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 500
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/fluid_solid_kernel/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -326,6 +329,9 @@ nbregions                       = 3              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 100
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/initial_mode_LDDRK/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -311,6 +314,9 @@ nbregions                       = 1              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 1000
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/initial_plane_wave_with_free_surface/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -310,6 +313,9 @@ nbregions                       = 1              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 200
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/meshing/Gmsh_example_CPML_MPI/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -333,6 +336,9 @@ nbregions                       = 2              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 1000
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/meshing/Gmsh_example_Stacey_MPI/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -323,6 +326,9 @@ nbregions                       = 2              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 300
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/meshing/Gmsh_example_Stacey_MPI/DATA/Par_file.serial

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -323,6 +326,9 @@ nbregions                       = 2              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 300
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/meshing/Marmousi_mesh_of_the_model/CUBIT_meshing/4_export_mesh.py

@@ -212,9 +212,9 @@
     #       this will help when creating a new nummaterial_velocity_file for Marmousi2 with different properties for each surface.
     print("# each surface with its own block material definition")
     print("#")
-    print("# using elastic material definition")    
+    print("# using elastic material definition")
     print("#")
-    
+
     for i in surfaces:
         isurf = i  # surface ids start at 1
         print("# surface: ",isurf)

EXAMPLES/applications/meshing/Marmousi_mesh_of_the_model/DATA/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -322,6 +325,9 @@ nbregions                       = 1              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 2000
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/meshing/Marmousi_mesh_of_the_model/DATA/Par_file.mesh-improved

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -322,6 +325,9 @@ nbregions                       = 1              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 2000
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not

EXAMPLES/applications/meshing/Marmousi_mesh_of_the_model/REF_SEIS/Par_file

@@ -142,6 +142,9 @@ noise_source_time_function_type = 4
 # This is only useful for GPU version (for now)
 write_moving_sources_database   = .false.
 
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .true.
+
 #-----------------------------------------------------------
 #
 # Receivers
@@ -322,6 +325,9 @@ nbregions                       = 1              # then set below the different
 # (every how many time steps we display information about the simulation. costly, do not use a very small value)
 NTSTEP_BETWEEN_OUTPUT_INFO      = 2000
 
+# save mesh files (VTK, ASCII or Gnuplot) to check the mesh
+SAVE_MESH_FILES                 = .true.
+
 # meshing output
 output_grid_Gnuplot             = .false.        # generate a GNUPLOT file containing the grid, and a script to plot it
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not