[cig-commits] r22479 - in seismo/3D/SPECFEM3D_GLOBE: branches/SPECFEM3D_GLOBE_SUNFLOWER/DATA branches/SPECFEM3D_GLOBE_SUNFLOWER/setup branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries trunk/DATA trunk/setup trunk/src/auxiliaries trunk/src/specfem3D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Mon Jul 1 08:41:04 PDT 2013
Author: dkomati1
Date: 2013-07-01 08:41:04 -0700 (Mon, 01 Jul 2013)
New Revision: 22479
Modified:
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/setup/constants.h.in
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/check_buffers_2D.f90
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/combine_AVS_DX.f90
seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/combine_vol_data.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file
seismo/3D/SPECFEM3D_GLOBE/trunk/setup/constants.h.in
seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/check_buffers_1D.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/check_buffers_2D.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/create_movie_AVS_DX.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/create_movie_GMT_global.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_boundary_kernel.f90
Log:
done merging all the easy changes in src/auxiliaries
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/DATA/Par_file 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/DATA/Par_file 2013-07-01 15:41:04 UTC (rev 22479)
@@ -1,27 +1,27 @@
# 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 = .true. # save last frame of forward simulation or not
+SAVE_FORWARD = .false. # save last frame of forward simulation or not
# number of chunks (1,2,3 or 6)
-NCHUNKS = 1
+NCHUNKS = 6
# 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
-CENTER_LATITUDE_IN_DEGREES = 40.d0
-CENTER_LONGITUDE_IN_DEGREES = 25.d0
-GAMMA_ROTATION_AZIMUTH = 0.d0
+CENTER_LATITUDE_IN_DEGREES = 90.d0
+CENTER_LONGITUDE_IN_DEGREES = 0.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 = 240
+NEX_ETA = 240
# number of MPI processors along the two sides of the first chunk
-NPROC_XI = 2
-NPROC_ETA = 2
+NPROC_XI = 5
+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
@@ -31,9 +31,14 @@
#
# 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 = s362ani
# parameters describing the Earth model
OCEANS = .true.
@@ -48,13 +53,13 @@
ABSORBING_CONDITIONS = .true.
# record length in minutes
-RECORD_LENGTH_IN_MINUTES = 0.1d0
+RECORD_LENGTH_IN_MINUTES = 15.0d0
# save AVS or OpenDX movies
MOVIE_SURFACE = .false.
MOVIE_VOLUME = .false.
-MOVIE_COARSE = .true.
-NTSTEP_BETWEEN_FRAMES = 50
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 100
HDUR_MOVIE = 0.d0
# save movie in volume. Will save element if center of element is in prescribed volume
@@ -63,7 +68,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
@@ -75,9 +79,9 @@
MOVIE_STOP = 40000
# save mesh files to check the mesh
-SAVE_MESH_FILES = .true.
+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
@@ -88,7 +92,7 @@
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 = 1000
# interval in time steps for temporary writing of seismograms
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
@@ -103,7 +107,7 @@
ROTATE_SEISMOGRAMS_RT = .false.
# decide if master process writes all the seismograms or if all processes do it in parallel
-WRITE_SEISMOGRAMS_BY_MASTER = .false.
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
# save all seismograms in one large combined file instead of one file per seismogram
# to avoid overloading shared non-local file systems such as GPFS for instance
@@ -119,7 +123,7 @@
# set to true to use GPUs
GPU_MODE = .false.
-# set to true to use GPUs
+# set to true to use ADIOS
ADIOS_ENABLED = .true.
ADIOS_FOR_FORWARD_ARRAYS = .true.
ADIOS_FOR_MPI_ARRAYS = .true.
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/setup/constants.h.in
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/setup/constants.h.in 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/setup/constants.h.in 2013-07-01 15:41:04 UTC (rev 22479)
@@ -68,6 +68,8 @@
! I/O unit for source and receiver vtk file
integer, parameter :: IOVTK = 98
+! I/O Bluegene serial access
+ logical, parameter :: IO_BLUEGENE_SERIAL = .false.
! R_EARTH is the radius of the bottom of the oceans (radius of Earth in m)
double precision, parameter :: R_EARTH = 6371000.d0
@@ -192,8 +194,9 @@
! was found by trial and error
double precision, parameter :: SOURCE_DECAY_MIMIC_TRIANGLE = 1.628d0
-! maximum number of sources to locate simultaneously
+! maximum number of sources and receivers to locate simultaneously
integer, parameter :: NSOURCES_SUBSET_MAX = 100
+ integer, parameter :: NREC_SUBSET_MAX = 200
! use a force source located exactly at a grid point instead of a CMTSOLUTION source
! this can be useful e.g. for asteroid impact simulations
@@ -271,10 +274,19 @@
! (default is to use transverse isotropy only between MOHO and 220 )
logical, parameter :: USE_FULL_TISO_MANTLE = .false.
+! regular kernel parameters
+ character(len=*), parameter :: PATHNAME_KL_REG = 'DATA/kl_reg_grid.txt'
+ integer, parameter :: NM_KL_REG_LAYER = 100
+ integer, parameter :: NM_KL_REG_PTS = 300000
+ real, parameter :: KL_REG_MIN_LON = 0.0
+ real, parameter :: KL_REG_MAX_LON = 360.0
+ real, parameter :: KL_REG_MIN_LAT = -90.0
+ real, parameter :: KL_REG_MAX_LAT = +90.0
+
! forces attenuation arrays to be fully 3D and defined on all GLL points
! (useful for more accurate attenuation profiles and adjoint inversions)
! (if set to .false., 3D attenuation is only used for models with 3D attenuation distributions)
- logical, parameter :: USE_3D_ATTENUATION_ARRAYS = .true.
+ logical, parameter :: USE_3D_ATTENUATION_ARRAYS = .false.
!!-----------------------------------------------------------
@@ -656,3 +668,29 @@
integer, parameter :: NSPEC_DOUBLING_BASICBRICK = 8
integer, parameter :: NGLOB_DOUBLING_BASICBRICK = 27
+! this for LDDRK
+ logical, parameter :: USE_LDDRK = .false. ! .true.
+
+! the maximum CFL of LDDRK is significantly higher than that of the Newmark scheme,
+! in a ratio that is theoretically 1.327 / 0.697 = 1.15 / 0.604 = 1.903 for a solid with Poisson's ratio = 0.25
+! and for a fluid (see the manual of the 2D code, SPECFEM2D, Tables 4.1 and 4.2, and that ratio does not
+! depend on whether we are in 2D or in 3D). However in practice a ratio of about 1.5 to 1.7 is often safer
+! (for instance for models with a large range of Poisson's ratio values).
+! Since the SPECFEM3D_GLOBE code computes the time step using the Newmark scheme, for LDDRK we will simply
+! multiply that time step by this ratio when LDDRK is on and when flag INCREASE_CFL_FOR_LDDRK is true.
+ logical, parameter :: INCREASE_CFL_FOR_LDDRK = .true.
+ double precision, parameter :: RATIO_BY_WHICH_TO_INCREASE_IT = 1.5d0
+
+ integer, parameter :: NSTAGE = 6
+ real(kind=CUSTOM_REAL), parameter, dimension(NSTAGE) :: ALPHA_LDDRK = &
+ (/0.0_CUSTOM_REAL,-0.737101392796_CUSTOM_REAL, -1.634740794341_CUSTOM_REAL,&
+ -0.744739003780_CUSTOM_REAL,-1.469897351522_CUSTOM_REAL,-2.813971388035_CUSTOM_REAL/)
+
+ real(kind=CUSTOM_REAL), parameter, dimension(NSTAGE) :: BETA_LDDRK = &
+ (/0.032918605146_CUSTOM_REAL,0.823256998200_CUSTOM_REAL,0.381530948900_CUSTOM_REAL,&
+ 0.200092213184_CUSTOM_REAL,1.718581042715_CUSTOM_REAL,0.27_CUSTOM_REAL/)
+
+ real(kind=CUSTOM_REAL), parameter, dimension(NSTAGE) :: C_LDDRK = &
+ (/0.0_CUSTOM_REAL,0.032918605146_CUSTOM_REAL,0.249351723343_CUSTOM_REAL,&
+ 0.466911705055_CUSTOM_REAL,0.582030414044_CUSTOM_REAL,0.847252983783_CUSTOM_REAL/)
+
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/check_buffers_2D.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/check_buffers_2D.f90 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/check_buffers_2D.f90 2013-07-01 15:41:04 UTC (rev 22479)
@@ -138,7 +138,7 @@
MOVIE_VOLUME,MOVIE_COARSE,ATTENUATION_3D,RECEIVERS_CAN_BE_BURIED, &
PRINT_SOURCE_TIME_FUNCTION,SAVE_MESH_FILES, &
ATTENUATION,REFERENCE_1D_MODEL,THREE_D_MODEL,ABSORBING_CONDITIONS, &
- INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE,&
+ INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE, &
LOCAL_PATH,LOCAL_TMP_PATH,MODEL, &
SIMULATION_TYPE,SAVE_FORWARD, &
NPROC,NPROCTOT,NEX_PER_PROC_XI,NEX_PER_PROC_ETA, &
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/combine_AVS_DX.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/combine_AVS_DX.f90 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/combine_AVS_DX.f90 2013-07-01 15:41:04 UTC (rev 22479)
@@ -1006,13 +1006,15 @@
if(ivalue == 3) then
print *
print *,'elevation min, max = ',minval(elevation_sphere),maxval(elevation_sphere)
- if(TOPOGRAPHY .or. CRUSTAL) print *,'elevation has been normalized for topography'
- print *
+ if(TOPOGRAPHY .or. CRUSTAL) then
+ print *,'elevation has been normalized for topography'
+ print *
! print percentage of oceans at surface of the globe
+ print *
+ print *,'the oceans represent ',100. * below_zero / (above_zero + below_zero),' % of the surface of the mesh'
+ endif
print *
- print *,'the oceans represent ',100. * below_zero / (above_zero + below_zero),' % of the surface of the mesh'
- print *
endif
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/combine_vol_data.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/combine_vol_data.f90 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/auxiliaries/combine_vol_data.f90 2013-07-01 15:41:04 UTC (rev 22479)
@@ -209,6 +209,7 @@
print *, 'nspec(it) = ', nspec(1:num_node)
print *, 'nglob(it) = ', nglob(1:num_node)
+ print *, 'nelement(it) = ', nelement(1:num_node)
call write_integer_fd(efd,sum(nelement(1:num_node)))
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file 2013-07-01 15:41:04 UTC (rev 22479)
@@ -49,7 +49,7 @@
ATTENUATION = .true.
# absorbing boundary conditions for a regional simulation
-ABSORBING_CONDITIONS = .false.
+ABSORBING_CONDITIONS = .true.
# record length in minutes
RECORD_LENGTH_IN_MINUTES = 15.0d0
@@ -62,8 +62,6 @@
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 = .false.
@@ -94,6 +92,7 @@
NUMBER_OF_THIS_RUN = 1
# path to store the local database files on each node
+# mesh files
LOCAL_PATH = ./DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/setup/constants.h.in
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/setup/constants.h.in 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/setup/constants.h.in 2013-07-01 15:41:04 UTC (rev 22479)
@@ -86,7 +86,6 @@
! I/O unit for source and receiver vtk file
integer, parameter :: IOVTK = 98
-
! I/O Bluegene serial access
logical, parameter :: IO_BLUEGENE_SERIAL = .false.
@@ -253,7 +252,21 @@
! the mesher) and use them for the computation of boundary kernel (in the solver)
logical, parameter :: SAVE_BOUNDARY_MESH = .false.
+! mesh coloring:
+! add mesh coloring for the GPU + MPI implementation
+ logical, parameter :: USE_MESH_COLORING_GPU = .false.
+ integer, parameter :: MAX_NUMBER_OF_COLORS = 1000
+! enhanced coloring:
+! using Droux algorithm
+! try several times with one more color before giving up
+ logical, parameter :: USE_DROUX_OPTIMIZATION = .false.
+ integer, parameter :: MAX_NB_TRIES_OF_DROUX_1993 = 15
+
+! using balancing algorithm
+! postprocess the colors to balance them if Droux (1993) algorithm is not used
+ logical, parameter :: BALANCE_COLORS_SIMPLE_ALGO = .true.
+
!!-----------------------------------------------------------
!!
!! adjoint kernel outputs
@@ -288,6 +301,12 @@
real, parameter :: KL_REG_MIN_LAT = -90.0
real, parameter :: KL_REG_MAX_LAT = +90.0
+! forces attenuation arrays to be fully 3D and defined on all GLL points
+! (useful for more accurate attenuation profiles and adjoint inversions)
+! (if set to .false., 3D attenuation is only used for models with 3D attenuation distributions)
+ logical, parameter :: USE_3D_ATTENUATION_ARRAYS = .false.
+
+
!!-----------------------------------------------------------
!!
!! time stamp information
@@ -342,9 +361,12 @@
double precision, parameter :: PI = 3.141592653589793d0
double precision, parameter :: TWO_PI = 2.d0 * PI
double precision, parameter :: PI_OVER_FOUR = PI / 4.d0
+ double precision, parameter :: PI_OVER_TWO = PI / 2.0d0
! to convert angles from degrees to radians
double precision, parameter :: DEGREES_TO_RADIANS = PI / 180.d0
+! to convert angles from radians to degrees
+ double precision, parameter :: RADIANS_TO_DEGREES = 180.d0 / PI
! 3-D simulation
integer, parameter :: NDIM = 3
@@ -468,6 +490,16 @@
integer, parameter :: THREE_D_MODEL_S40RTS = 11
integer, parameter :: THREE_D_MODEL_GAPP2 = 12
+! define flag for regions of the global Earth for attenuation
+ integer, parameter :: NUM_REGIONS_ATTENUATION = 5
+
+ integer, parameter :: IREGION_ATTENUATION_INNER_CORE = 1
+ integer, parameter :: IREGION_ATTENUATION_CMB_670 = 2
+ integer, parameter :: IREGION_ATTENUATION_670_220 = 3
+ integer, parameter :: IREGION_ATTENUATION_220_80 = 4
+ integer, parameter :: IREGION_ATTENUATION_80_SURFACE = 5
+ integer, parameter :: IREGION_ATTENUATION_UNDEFINED = 6
+
! number of standard linear solids for attenuation
integer, parameter :: N_SLS = 3
@@ -477,6 +509,10 @@
integer, parameter :: ATTENUATION_COMP_RESOLUTION = 1
integer, parameter :: ATTENUATION_COMP_MAXIMUM = 5000
+! for lookup table for attenuation every 100 m in radial direction of Earth model
+ integer, parameter :: NRAD_ATTENUATION = 70000
+ double precision, parameter :: TABLE_ATTENUATION = R_EARTH_KM * 10.0d0
+
! for determination of the attenuation period range
! if this is set to .true. then the hardcoded values will be used
! otherwise they are computed automatically from the Number of elements
@@ -624,8 +660,12 @@
! parameters for GLL model (used for iterative model inversions)
character(len=*), parameter :: PATHNAME_GLL_modeldir = 'DATA/GLL/'
+
+ integer, parameter :: GLL_REFERENCE_1D_MODEL = REFERENCE_MODEL_1DREF
integer, parameter :: GLL_REFERENCE_MODEL = THREE_D_MODEL_S29EA
+
!!-- uncomment for using S362ANI as reference model (Hejun Zhu)
+! integer, parameter :: GLL_REFERENCE_1D_MODEL = REFERENCE_MODEL_1DREF
! integer, parameter :: GLL_REFERENCE_MODEL = THREE_D_MODEL_S362ANI
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/check_buffers_1D.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/check_buffers_1D.f90 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/check_buffers_1D.f90 2013-07-01 15:41:04 UTC (rev 22479)
@@ -139,7 +139,9 @@
MOVIE_VOLUME,MOVIE_COARSE,ATTENUATION_3D,RECEIVERS_CAN_BE_BURIED, &
PRINT_SOURCE_TIME_FUNCTION,SAVE_MESH_FILES, &
ATTENUATION,REFERENCE_1D_MODEL,THREE_D_MODEL,ABSORBING_CONDITIONS, &
- INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE,LOCAL_PATH,MODEL,SIMULATION_TYPE,SAVE_FORWARD, &
+ INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE, &
+ LOCAL_PATH,MODEL, &
+ SIMULATION_TYPE,SAVE_FORWARD, &
NPROC,NPROCTOT,NEX_PER_PROC_XI,NEX_PER_PROC_ETA, &
NSPEC, &
NSPEC2D_XI, &
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/check_buffers_2D.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/check_buffers_2D.f90 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/check_buffers_2D.f90 2013-07-01 15:41:04 UTC (rev 22479)
@@ -141,7 +141,9 @@
MOVIE_VOLUME,MOVIE_COARSE,ATTENUATION_3D,RECEIVERS_CAN_BE_BURIED, &
PRINT_SOURCE_TIME_FUNCTION,SAVE_MESH_FILES, &
ATTENUATION,REFERENCE_1D_MODEL,THREE_D_MODEL,ABSORBING_CONDITIONS, &
- INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE,LOCAL_PATH,MODEL,SIMULATION_TYPE,SAVE_FORWARD, &
+ INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE, &
+ LOCAL_PATH,MODEL, &
+ SIMULATION_TYPE,SAVE_FORWARD, &
NPROC,NPROCTOT,NEX_PER_PROC_XI,NEX_PER_PROC_ETA, &
NSPEC, &
NSPEC2D_XI, &
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/create_movie_AVS_DX.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/create_movie_AVS_DX.f90 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/create_movie_AVS_DX.f90 2013-07-01 15:41:04 UTC (rev 22479)
@@ -842,7 +842,9 @@
MOVIE_VOLUME,MOVIE_COARSE,ATTENUATION_3D,RECEIVERS_CAN_BE_BURIED, &
PRINT_SOURCE_TIME_FUNCTION,SAVE_MESH_FILES, &
ATTENUATION,REFERENCE_1D_MODEL,THREE_D_MODEL,ABSORBING_CONDITIONS, &
- INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE,LOCAL_PATH,MODEL,SIMULATION_TYPE,SAVE_FORWARD, &
+ INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE, &
+ LOCAL_PATH,MODEL, &
+ SIMULATION_TYPE,SAVE_FORWARD, &
NPROC,NPROCTOT,NEX_PER_PROC_XI,NEX_PER_PROC_ETA, &
NSPEC,NSPEC2D_XI,NSPEC2D_ETA, &
NSPEC2DMAX_XMIN_XMAX,NSPEC2DMAX_YMIN_YMAX,NSPEC2D_BOTTOM,NSPEC2D_TOP, &
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/create_movie_GMT_global.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/create_movie_GMT_global.f90 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/create_movie_GMT_global.f90 2013-07-01 15:41:04 UTC (rev 22479)
@@ -191,7 +191,9 @@
MOVIE_VOLUME,MOVIE_COARSE,ATTENUATION_3D,RECEIVERS_CAN_BE_BURIED, &
PRINT_SOURCE_TIME_FUNCTION,SAVE_MESH_FILES, &
ATTENUATION,REFERENCE_1D_MODEL,THREE_D_MODEL,ABSORBING_CONDITIONS, &
- INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE,LOCAL_PATH,MODEL,SIMULATION_TYPE,SAVE_FORWARD, &
+ INCLUDE_CENTRAL_CUBE,INFLATE_CENTRAL_CUBE, &
+ LOCAL_PATH,MODEL, &
+ SIMULATION_TYPE,SAVE_FORWARD, &
NPROC,NPROCTOT,NEX_PER_PROC_XI,NEX_PER_PROC_ETA, &
NSPEC,NSPEC2D_XI,NSPEC2D_ETA, &
NSPEC2DMAX_XMIN_XMAX,NSPEC2DMAX_YMIN_YMAX,NSPEC2D_BOTTOM,NSPEC2D_TOP, &
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_boundary_kernel.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_boundary_kernel.f90 2013-07-01 12:46:50 UTC (rev 22478)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_boundary_kernel.f90 2013-07-01 15:41:04 UTC (rev 22479)
@@ -63,7 +63,7 @@
real(kind=CUSTOM_REAL) :: b_kl(NGLLX,NGLLY,NSPEC2D_DISC)
logical :: fluid_solid_boundary
-! --- local variables ---
+ ! --- local variables ---
integer ispec2D,i,j,k,iglob,ispec
real(kind=CUSTOM_REAL), dimension(NDIM,NGLLX,NGLLY,NGLLZ) :: displl, accell, b_displl, Kdvect
real(kind=CUSTOM_REAL), dimension(NDIM) :: normal, temp1, temp2, temp3
More information about the CIG-COMMITS
mailing list