[cig-commits] r19201 - in seismo/2D/SPECFEM2D/trunk: EXAMPLES/canyon src/specfem2D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2011-11-14 16:25:52 -0800 (Mon, 14 Nov 2011)
New Revision: 19201

Modified:
   seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
done making the periodic boundary conditions more flexible


Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon	2011-11-15 00:06:36 UTC (rev 19200)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon	2011-11-15 00:25:52 UTC (rev 19201)
@@ -90,6 +90,15 @@
 # absorbing boundary active or not
 absorbing_conditions            = .true.
 
+# for horizontal periodic conditions: detect common points between left and right edges
+ADD_PERIODIC_CONDITIONS         = .false.
+
+# horizontal periodicity distance for periodic conditions
+PERIODIC_horiz_dist             = 0.3597d0
+
+# grid point detection tolerance for periodic conditions
+PERIODIC_DETECT_TOL             = 3.3334d-6
+
 #-----------------------------------------------------------------------------
 # PARAMETERS FOR EXTERNAL MESHING
 

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90	2011-11-15 00:06:36 UTC (rev 19200)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90	2011-11-15 00:25:52 UTC (rev 19201)
@@ -55,7 +55,8 @@
                   TURN_VISCATTENUATION_ON,Q0,freq0,p_sv, &
                   NSTEP,deltat,NTSTEP_BETWEEN_OUTPUT_SEISMO,NSOURCES, &
                   factor_subsample_image,USE_SNAPSHOT_NUMBER_IN_FILENAME,DRAW_WATER_CONSTANT_BLUE_IN_JPG,US_LETTER, &
-                  POWER_DISPLAY_COLOR,PERFORM_CUTHILL_MCKEE,SU_FORMAT,USER_T0)
+                  POWER_DISPLAY_COLOR,PERFORM_CUTHILL_MCKEE,SU_FORMAT,USER_T0, &
+                  ADD_PERIODIC_CONDITIONS,PERIODIC_horiz_dist,PERIODIC_DETECT_TOL)
 
 ! starts reading in parameters from input Database file
 
@@ -108,6 +109,15 @@
 !  simulation will start at t = - t0)
   double precision :: USER_T0
 
+!! DK DK for horizontal periodic conditions: detect common points between left and right edges
+  logical :: ADD_PERIODIC_CONDITIONS
+
+!! DK DK horizontal periodicity distance for periodic conditions
+  double precision :: PERIODIC_horiz_dist
+
+!! DK DK grid point detection tolerance for periodic conditions
+  double precision :: PERIODIC_DETECT_TOL
+
   ! local parameters
   integer :: ier
   character(len=80) :: datlin

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2011-11-15 00:06:36 UTC (rev 19200)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2011-11-15 00:25:52 UTC (rev 19201)
@@ -795,7 +795,10 @@
 ! to help locate elements with a negative Jacobian using OpenDX
   logical :: found_a_negative_jacobian
 
-!! DK DK the horizontal periodicity distance for periodic conditions
+!! DK DK for horizontal periodic conditions: detect common points between left and right edges
+  logical :: ADD_PERIODIC_CONDITIONS
+
+!! DK DK horizontal periodicity distance for periodic conditions
   double precision :: PERIODIC_horiz_dist
 
 !! DK DK grid point detection tolerance for periodic conditions
@@ -905,7 +908,8 @@
                   TURN_VISCATTENUATION_ON,Q0,freq0,p_sv, &
                   NSTEP,deltat,NTSTEP_BETWEEN_OUTPUT_SEISMO,NSOURCES, &
                   factor_subsample_image,USE_SNAPSHOT_NUMBER_IN_FILENAME,DRAW_WATER_CONSTANT_BLUE_IN_JPG,US_LETTER, &
-                  POWER_DISPLAY_COLOR,PERFORM_CUTHILL_MCKEE,SU_FORMAT,USER_T0)
+                  POWER_DISPLAY_COLOR,PERFORM_CUTHILL_MCKEE,SU_FORMAT,USER_T0, &
+                  ADD_PERIODIC_CONDITIONS,PERIODIC_horiz_dist,PERIODIC_DETECT_TOL)
 
   !
   !--- source information