[cig-commits] [commit] devel: Remove unused constants. (8c28984)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Thu May 1 00:50:18 PDT 2014


Repository : ssh://geoshell/specfem3d

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d/compare/cb32c88d6155d7974561a6f72fc17aea596e2c4d...50aa953c1db3f565d76415f5305410a529996b75

>---------------------------------------------------------------

commit 8c28984862717ad1276aae68f61f2aeebf992bd4
Author: Elliott Sales de Andrade <esalesde at physics.utoronto.ca>
Date:   Fri Apr 25 00:47:02 2014 -0400

    Remove unused constants.


>---------------------------------------------------------------

8c28984862717ad1276aae68f61f2aeebf992bd4
 setup/constants.h.in                                  | 12 +-----------
 src/generate_databases/pml_set_local_dampingcoeff.f90 |  4 ++--
 src/meshfem3D/meshfem3D.f90                           |  2 --
 src/specfem3D/compute_forces_poro_fluid_part.f90      |  3 +--
 src/specfem3D/pml_compute_accel_contribution.f90      |  6 ++----
 src/specfem3D/pml_compute_memory_variables.f90        | 12 ++++--------
 6 files changed, 10 insertions(+), 29 deletions(-)

diff --git a/setup/constants.h.in b/setup/constants.h.in
index 09c9756..0540c2b 100644
--- a/setup/constants.h.in
+++ b/setup/constants.h.in
@@ -51,9 +51,6 @@
 ! if running on a Beowulf, also modify name of nodes in filter_machine_file.f90
   logical, parameter :: LOCAL_PATH_IS_ALSO_GLOBAL = . at LOCAL_PATH_IS_ALSO_GLOBAL@.
 
-! apply heuristic rule to modify doubling regions to balance angles
-  logical, parameter :: APPLY_HEURISTIC_RULE = .true.
-
 ! use inlined products of Deville et al. (2002) to speedup the calculations to compute internal forces
   logical, parameter :: USE_DEVILLE_PRODUCTS = .false.
 
@@ -97,10 +94,8 @@
   integer, parameter :: IOABS_AC = 1
 ! I/O unit for plotting source time function
   integer, Parameter :: IOSTF = 71
-! I/O unit for interface file
-  integer, parameter :: IIN_INTERFACES = 43
 ! I/O unit for noise simulations
-  integer, parameter :: IIN_NOISE = 44,IOUT_NOISE = 45
+  integer, parameter :: IIN_NOISE = 44, IOUT_NOISE = 45
 ! I/O unit for SU formatted seismograms
   integer, parameter :: IOUT_SU = 46
 ! I/O unit for SU formatted adjoint sources
@@ -168,9 +163,6 @@
 ! (INRIA research report section 6.1:  http://hal.inria.fr/docs/00/07/32/19/PDF/RR-3471.pdf)
   real(kind=CUSTOM_REAL), parameter :: CPML_Rcoef = 0.001_CUSTOM_REAL
 
-! empirical damping modified factor
-  real(kind=CUSTOM_REAL), parameter :: damping_factor = 1.2_CUSTOM_REAL
-
 ! flags for the seven CPML regions
   integer, parameter :: CPML_X_ONLY = 1
   integer, parameter :: CPML_Y_ONLY = 2
@@ -416,8 +408,6 @@
   double precision, parameter :: GRAV = 6.6723d-11
 ! number of layers in PREM
   integer, parameter :: NR = 640
-! for lookup table for gravity every 100 m in radial direction of Earth model
-  integer, parameter :: NRAD_GRAVITY = 70000
 ! R_EARTH is the radius of the bottom of the oceans (radius of Earth in m)
   double precision, parameter :: R_EARTH = 6371000.d0
 ! same radius in km
diff --git a/src/generate_databases/pml_set_local_dampingcoeff.f90 b/src/generate_databases/pml_set_local_dampingcoeff.f90
index 3819ae3..7062155 100644
--- a/src/generate_databases/pml_set_local_dampingcoeff.f90
+++ b/src/generate_databases/pml_set_local_dampingcoeff.f90
@@ -35,7 +35,7 @@ subroutine pml_set_local_dampingcoeff(myrank,xstore,ystore,zstore)
                                     K_store_x,K_store_y,K_store_z,alpha_store_x,alpha_store_y,alpha_store_z,CPML_to_spec, &
                                     CPML_width_x,CPML_width_y,CPML_width_z,NPOWER,&
                                     CUSTOM_REAL,SIZE_REAL,NGLLX,NGLLY,NGLLZ,nspec_cpml,PML_INSTEAD_OF_FREE_SURFACE, &
-                                    IMAIN,FOUR_THIRDS,CPML_REGIONS,f0_FOR_PML,PI, &
+                                    IMAIN,CPML_REGIONS,f0_FOR_PML,PI, &
                                     CPML_X_ONLY,CPML_Y_ONLY,CPML_Z_ONLY,CPML_XY_ONLY,CPML_XZ_ONLY,CPML_YZ_ONLY,CPML_XYZ,&
                                     SIMULATION_TYPE,SAVE_FORWARD,nspec => NSPEC_AB,is_CPML,&
                                     mask_ibool_interior_domain,nglob_interface_PML_acoustic,points_interface_PML_acoustic,&
@@ -1471,7 +1471,7 @@ function pml_damping_profile_l(myrank,iglob,dist,vp,delta)
   !   vp:    P-velocity
   !   delta: thickness of the C-PML layer
 
-  use generate_databases_par, only: CUSTOM_REAL,NPOWER,CPML_Rcoef,damping_factor
+  use generate_databases_par, only: CUSTOM_REAL,NPOWER,CPML_Rcoef
 
   implicit none
 
diff --git a/src/meshfem3D/meshfem3D.f90 b/src/meshfem3D/meshfem3D.f90
index 772e7af..c7e98d9 100644
--- a/src/meshfem3D/meshfem3D.f90
+++ b/src/meshfem3D/meshfem3D.f90
@@ -785,8 +785,6 @@
 
   enddo
 
-  close(IIN_INTERFACES)
-
   if(myrank == 0) then
     write(IMAIN,*)
     write(IMAIN,*) '**************************'
diff --git a/src/specfem3D/compute_forces_poro_fluid_part.f90 b/src/specfem3D/compute_forces_poro_fluid_part.f90
index cb2fd32..1558852 100644
--- a/src/specfem3D/compute_forces_poro_fluid_part.f90
+++ b/src/specfem3D/compute_forces_poro_fluid_part.f90
@@ -44,8 +44,7 @@
 ! compute forces for the fluid poroelastic part
 
   use constants,only: CUSTOM_REAL,NGLLX,NGLLY,NGLLZ,NDIM, &
-                      N_SLS, &
-                      ONE_THIRD,FOUR_THIRDS
+                      ONE_THIRD
 
   implicit none
 
diff --git a/src/specfem3D/pml_compute_accel_contribution.f90 b/src/specfem3D/pml_compute_accel_contribution.f90
index 77ffb58..9c0a6af 100644
--- a/src/specfem3D/pml_compute_accel_contribution.f90
+++ b/src/specfem3D/pml_compute_accel_contribution.f90
@@ -40,8 +40,7 @@ subroutine pml_compute_accel_contribution_elastic(ispec,ispec_CPML,displ,veloc,r
   use pml_par, only: CPML_regions,d_store_x,d_store_y,d_store_z,K_store_x,K_store_y,K_store_z,&
                      alpha_store_x, alpha_store_y, alpha_store_z, &
                      NSPEC_CPML,accel_elastic_CPML,displ_old
-  use constants, only: CUSTOM_REAL,NDIM,NGLLX,NGLLY,NGLLZ,CPML_X_ONLY,CPML_Y_ONLY,CPML_Z_ONLY, &
-                       CPML_XY_ONLY,CPML_XZ_ONLY,CPML_YZ_ONLY,CPML_XYZ
+  use constants, only: CUSTOM_REAL,NDIM,NGLLX,NGLLY,NGLLZ
 
   implicit none
 
@@ -200,8 +199,7 @@ subroutine pml_compute_accel_contribution_acoustic(ispec,ispec_CPML,potential_ac
   use pml_par, only: CPML_regions,NSPEC_CPML,d_store_x,d_store_y,d_store_z,K_store_x,K_store_y,K_store_z,&
                      alpha_store_x, alpha_store_y, alpha_store_z, &
                      NSPEC_CPML,potential_dot_dot_acoustic_CPML,potential_acoustic_old
-  use constants, only: CUSTOM_REAL,NGLLX,NGLLY,NGLLZ,CPML_X_ONLY,CPML_Y_ONLY,CPML_Z_ONLY, &
-                       CPML_XY_ONLY,CPML_XZ_ONLY,CPML_YZ_ONLY,CPML_XYZ
+  use constants, only: CUSTOM_REAL,NGLLX,NGLLY,NGLLZ
 
   implicit none
 
diff --git a/src/specfem3D/pml_compute_memory_variables.f90 b/src/specfem3D/pml_compute_memory_variables.f90
index 78e0cde..f41b187 100644
--- a/src/specfem3D/pml_compute_memory_variables.f90
+++ b/src/specfem3D/pml_compute_memory_variables.f90
@@ -52,8 +52,7 @@ subroutine pml_compute_memory_variables_elastic(ispec,ispec_CPML,tempx1,tempy1,t
                      PML_dux_dxl_old, PML_dux_dyl_old, PML_dux_dzl_old, &
                      PML_duy_dxl_old, PML_duy_dyl_old, PML_duy_dzl_old, &
                      PML_duz_dxl_old, PML_duz_dyl_old, PML_duz_dzl_old
-  use constants, only: CUSTOM_REAL,NGLLX,NGLLY,NGLLZ,FOUR_THIRDS, &
-                       CPML_X_ONLY,CPML_Y_ONLY,CPML_Z_ONLY,CPML_XY_ONLY,CPML_XZ_ONLY,CPML_YZ_ONLY,CPML_XYZ
+  use constants, only: CUSTOM_REAL,NGLLX,NGLLY,NGLLZ,FOUR_THIRDS
 
   implicit none
 
@@ -430,8 +429,7 @@ subroutine pml_compute_memory_variables_acoustic(ispec,ispec_CPML,temp1,temp2,te
                      PML_dpotential_dxl,PML_dpotential_dyl,PML_dpotential_dzl,&
                      PML_dpotential_dxl_old,PML_dpotential_dyl_old,PML_dpotential_dzl_old
 
-  use constants, only: CUSTOM_REAL,NGLLX,NGLLY,NGLLZ,FOUR_THIRDS,USE_DEVILLE_PRODUCTS, &
-                       CPML_X_ONLY,CPML_Y_ONLY,CPML_Z_ONLY,CPML_XY_ONLY,CPML_XZ_ONLY,CPML_YZ_ONLY,CPML_XYZ
+  use constants, only: CUSTOM_REAL,NGLLX,NGLLY,NGLLZ,USE_DEVILLE_PRODUCTS
 
   implicit none
 
@@ -643,8 +641,7 @@ subroutine pml_compute_memory_variables_acoustic_elastic(ispec_CPML,iface,iglob,
   use specfem_par, only: NGLOB_AB,it,deltat
   use pml_par,only : CPML_regions,k_store_x,k_store_y,k_store_z,d_store_x,d_store_y,d_store_z,&
                      alpha_store_x,alpha_store_y,alpha_store_z,displ_old
-  use constants, only: CUSTOM_REAL,NDIM,NGLLX,NGLLY,NGLLZ,&
-                       CPML_X_ONLY,CPML_Y_ONLY,CPML_Z_ONLY,CPML_XY_ONLY,CPML_XZ_ONLY,CPML_YZ_ONLY,CPML_XYZ
+  use constants, only: CUSTOM_REAL,NDIM,NGLLX,NGLLY,NGLLZ
 
   implicit none
 
@@ -746,8 +743,7 @@ subroutine pml_compute_memory_variables_elastic_acoustic(ispec_CPML,iface,iglob,
   use specfem_par, only: NGLOB_AB,it,deltat
   use pml_par,only : CPML_regions,k_store_x,k_store_y,k_store_z,d_store_x,d_store_y,d_store_z,&
                      alpha_store_x,alpha_store_y,alpha_store_z
-  use constants, only: CUSTOM_REAL,NDIM,NGLLX,NGLLY,NGLLZ,&
-                       CPML_X_ONLY,CPML_Y_ONLY,CPML_Z_ONLY,CPML_XY_ONLY,CPML_XZ_ONLY,CPML_YZ_ONLY,CPML_XYZ
+  use constants, only: CUSTOM_REAL,NDIM,NGLLX,NGLLY,NGLLZ
 
   implicit none
 



More information about the CIG-COMMITS mailing list