[cig-commits] [commit] devel,master: fixed the plane wave code and made sure the "EXAMPLES/canyon" example now works fine; also removed useless white spaces in source files. (44fbb26)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:19:21 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 44fbb26226c63e770fc994134d9a3ae62406482b
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Thu Sep 1 13:42:50 2011 +0000
fixed the plane wave code and made sure the "EXAMPLES/canyon" example now works fine;
also removed useless white spaces in source files.
>---------------------------------------------------------------
44fbb26226c63e770fc994134d9a3ae62406482b
canyon/Par_file_canyon | 31 ++++++++++++++++---------------
canyon/SOURCE_canyon | 6 +++---
2 files changed, 19 insertions(+), 18 deletions(-)
diff --git a/canyon/Par_file_canyon b/canyon/Par_file_canyon
index 6c299cb..a03cf07 100644
--- a/canyon/Par_file_canyon
+++ b/canyon/Par_file_canyon
@@ -1,9 +1,10 @@
# title of job
-title = Calcul Mexique Alejandro
+title = Canyon Mexico Paco P-SV
# forward or adjoint simulation
SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint + kernels
-SAVE_FORWARD = .true. # save the last frame, needed for adjoint simulation
+NOISE_TOMOGRAPHY = 0 # 0 = earthquake simulation, 1/2/3 = noise simulation
+SAVE_FORWARD = .false. # save the last frame, needed for adjoint simulation
# parameters concerning partitioning
nproc = 1 # number of processes
@@ -40,17 +41,17 @@ anglerec = 0.d0 # angle to rotate components at
rec_normal_to_surface = .false. # base anglerec normal to surface (external mesh and curve file needed)
# first receiver line
-nrec = 20 # number of receivers
-xdeb = 1. # first receiver x in meters
+nrec = 5 # number of receivers
+xdeb = 14.5 # first receiver x in meters
zdeb = 9. # first receiver z in meters
-xfin = 18. # last receiver x in meters (ignored if onlyone receiver)
+xfin = 17.5 # last receiver x in meters (ignored if onlyone receiver)
zfin = 9. # last receiver z in meters (ignored if onlyone receiver)
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
+output_color_image = .false. # output color image of the results
imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
cutsnaps = 1. # minimum amplitude in % for snapshots
meshvect = .true. # display mesh on vector plots or not
@@ -63,17 +64,17 @@ sizemax_arrows = 1.d0 # maximum size of arrows on vec
gnuplot = .false. # generate a GNUPLOT file for the grid
output_grid = .false. # save the grid in a text file or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
+output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
# define models as
-# I: (model_number 1 rho Vp Vs 0 0 QKappa Qmu 0 0 0 0 0 0) or
+# I: (model_number 1 rho Vp Vs 0 0 Qp Qs 0 0 0 0 0 0) or
# II: (model_number 2 rho c11 c13 c15 c33 c35 c55 0 0 0 0 0 0) or
-# III: (model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu).
+# III: (model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qs).
# For istropic elastic/acoustic material use I and set Vs to zero to make a given model acoustic, for anisotropic elastic use II,
# and for isotropic poroelastic material use III. The mesh can contain acoustic, elastic, and poroelastic models simultaneously.
-1 1 1.d0 2.d0 1.d0 0 0 10.d0 10.d0 0 0 0 0 0 0
+1 1 1.d0 2.d0 1.d0 0 0 9999.d0 9999.d0 0 0 0 0 0 0
# external mesh or not
read_external_mesh = .true.
@@ -86,11 +87,11 @@ absorbing_conditions = .true.
# data concerning mesh, when generated using third-party app (more info in README)
# (see also absorbing_conditions above)
-mesh_file = ./mesh/canyon_mesh_file # file containing the mesh
-nodes_coords_file = ./mesh/canyon_nodes_coords_file # file containing the nodes coordinates
-materials_file = ./mesh/canyon_materials_file # file containing the material number for each element
-free_surface_file = ./mesh/canyon_free_surface_file # file containing the free surface
-absorbing_surface_file = ./mesh/canyon_absorbing_surface_file # file containing the absorbing surface
+mesh_file = ./DATA/mesh/canyon_mesh_file # file containing the mesh
+nodes_coords_file = ./DATA/mesh/canyon_nodes_coords_file # file containing the nodes coordinates
+materials_file = ./DATA/mesh/canyon_materials_file # file containing the material number for each element
+free_surface_file = ./DATA/mesh/canyon_free_surface_file # file containing the free surface
+absorbing_surface_file = ./DATA/mesh/canyon_absorbing_surface_file # file containing the absorbing surface
tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
#-----------------------------------------------------------------------------
diff --git a/canyon/SOURCE_canyon b/canyon/SOURCE_canyon
index 787006e..e5ade1b 100644
--- a/canyon/SOURCE_canyon
+++ b/canyon/SOURCE_canyon
@@ -1,8 +1,8 @@
# source parameters
source_surf = .false. # source inside the medium or at the surface
-xs = 2. # source location x in meters
-zs = 3. # source location z in meters
-source_type = 1 # 1 for plane P waves, 2 for plane SV waves, 3 for Rayleigh wave
+xs = 8. # source location x in meters
+zs = 9. # source location z in meters
+source_type = 2 # 1 for plane P waves, 2 for plane SV waves, 3 for Rayleigh wave
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 1.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
More information about the CIG-COMMITS
mailing list