[cig-commits] r22538 - in seismo/3D/SPECFEM3D/trunk: doc/USER_MANUAL src/shared src/specfem3D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-07-08 05:13:21 -0700 (Mon, 08 Jul 2013)
New Revision: 22538

Modified:
   seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.pdf
   seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.tex
   seismo/3D/SPECFEM3D/trunk/src/shared/constants.h.in
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/locate_source.f90
Log:
moved flag USE_BEST_LOCATION to the "constants.h.in" file to make it easier to change it if/when needed, and added a line about this to the users manual. 


Modified: seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.pdf
===================================================================
(Binary files differ)

Modified: seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.tex
===================================================================
--- seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.tex	2013-07-08 00:29:40 UTC (rev 22537)
+++ seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.tex	2013-07-08 12:13:21 UTC (rev 22538)
@@ -1555,6 +1555,12 @@
 the solver and thus only its direction matters here; its norm is ignored
 and the norm of the force used is the factor force source times the
 source time function. 
+When using this option, by default the code can locate the force source
+anywhere between mesh points in order to honor its exact location; this
+is more precise than using the closest GLL mesh point, but it is also a bit slower.
+If needed, you can change that default behavior and force the code to use the closest
+GLL mesh point instead by setting flag \texttt{USE\_BEST\_LOCATION} to \texttt{.false.} instead of \texttt{.true.} in file
+\texttt{src/shared/constants.h.in} and running the \texttt{configure} script again and recompiling the code.
 \item [{\texttt{USE\_RICKER\_TIME\_FUNCTION}}] Turn this flag on to use
 a Ricker source time function instead of the source time functions
 set by default to represent a (tilted) \texttt{FORCESOLUTION} force

Modified: seismo/3D/SPECFEM3D/trunk/src/shared/constants.h.in
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/shared/constants.h.in	2013-07-08 00:29:40 UTC (rev 22537)
+++ seismo/3D/SPECFEM3D/trunk/src/shared/constants.h.in	2013-07-08 12:13:21 UTC (rev 22538)
@@ -125,6 +125,10 @@
 !!
 !!-----------------------------------------------------------
 
+! interpolates position to find best possible location closest to target position (default)
+! note: this can be turned off to speedup location search and locate the source at the closest GLL point instead
+  logical, parameter :: USE_BEST_LOCATION = .true.
+
 ! flag to print the details of source location
   logical, parameter :: SHOW_DETAILS_LOCATE_SOURCE = .false.
 

Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/locate_source.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/locate_source.f90	2013-07-08 00:29:40 UTC (rev 22537)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/locate_source.f90	2013-07-08 12:13:21 UTC (rev 22538)
@@ -163,13 +163,6 @@
 
   !-----------------------------------------------------------------------------------
 
-  ! interpolates position to find best possible location closest to target position (default)
-  ! note: this can be turned off to speedup location search 
-  !       and position source at closest GLL point
-  logical, parameter :: USE_BEST_LOCATION = .true.
-
-  !-----------------------------------------------------------------------------------
-
   ! get the base pathname for output files
   call get_value_string(OUTPUT_FILES, 'OUTPUT_FILES', OUTPUT_FILES_PATH(1:len_trim(OUTPUT_FILES_PATH)))