[cig-commits] r22917 - seismo/2D/SPECFEM2D/trunk/src/specfem2D
xie.zhinan at geodynamics.org
xie.zhinan at geodynamics.org
Wed Oct 2 17:20:11 PDT 2013
Author: xie.zhinan
Date: 2013-10-02 17:20:11 -0700 (Wed, 02 Oct 2013)
New Revision: 22917
Modified:
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
Log:
update several comments
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90 2013-10-02 22:20:57 UTC (rev 22916)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90 2013-10-03 00:20:11 UTC (rev 22917)
@@ -235,7 +235,7 @@
! derivatives of potential
PML_dux_dxl_old(i,j) = dux_dxi*xixl + dux_dgamma*gammaxl
PML_dux_dzl_old(i,j) = dux_dxi*xizl + dux_dgamma*gammazl
-
+ ! the subroutine of lik_parameter_computation is presented at the end of compute_forces_viscoelastic.F90
call lik_parameter_computation(time_n,deltat,kappa_z,beta_z,alpha_z,kappa_x,beta_x,alpha_x,&
CPML_region_local,31,A5,A6,A7,singularity_type_zx,bb_zx_1,bb_zx_2,&
coef0_zx_1,coef1_zx_1,coef2_zx_1,coef0_zx_2,coef1_zx_2,coef2_zx_2)
@@ -367,6 +367,7 @@
alpha_z = alpha_z_store(i,j,ispec_PML)
beta_x = alpha_x + d_x / kappa_x
beta_z = alpha_z + d_z / kappa_z
+ ! the subroutine of l_parameter_computation is presented at the end of compute_forces_viscoelastic.F90
call l_parameter_computation(time_n,deltat,kappa_x,beta_x,alpha_x,kappa_z,beta_z,alpha_z, &
CPML_region_local,A0,A1,A2,A3,A4,singularity_type,&
bb_1,coef0_1,coef1_1,coef2_1,bb_2,coef0_2,coef1_2,coef2_2)
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90 2013-10-02 22:20:57 UTC (rev 22916)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90 2013-10-03 00:20:11 UTC (rev 22917)
@@ -718,9 +718,9 @@
! 2004 paper and his 2007 book. See also file doc/problem_attenuation_reference_Specfem2D_fixed_by_Xie_Zhinan.pdf
e1_sum = 0._CUSTOM_REAL; e11_sum = 0._CUSTOM_REAL; e13_sum = 0._CUSTOM_REAL
do i_sls = 1,N_SLS
- e1_sum = e1_sum + e1(i,j,ispec,i_sls)
- e11_sum = e11_sum + e11(i,j,ispec,i_sls)
- e13_sum = e13_sum + e13(i,j,ispec,i_sls)
+ e1_sum = e1_sum + e1(i,j,ispec,i_sls)
+ e11_sum = e11_sum + e11(i,j,ispec,i_sls)
+ e13_sum = e13_sum + e13(i,j,ispec,i_sls)
enddo
sigma_xx = sigma_xx + lambdalplusmul_relaxed_viscoel * e1_sum + TWO * mul_relaxed_viscoelastic * e11_sum
@@ -729,10 +729,10 @@
sigma_zx = sigma_xz
if(PML_BOUNDARY_CONDITIONS .and. is_PML(ispec)) then
- sigma_xx = lambdaplus2mu_unrelaxed_elastic*dux_dxl + lambdal_unrelaxed_elastic*PML_duz_dzl(i,j)
- sigma_zz = lambdaplus2mu_unrelaxed_elastic*duz_dzl + lambdal_unrelaxed_elastic*PML_dux_dxl(i,j)
- sigma_zx = mul_unrelaxed_elastic * (PML_duz_dxl(i,j) + dux_dzl)
- sigma_xz = mul_unrelaxed_elastic * (PML_dux_dzl(i,j) + duz_dxl)
+ sigma_xx = lambdaplus2mu_unrelaxed_elastic*dux_dxl + lambdal_unrelaxed_elastic*PML_duz_dzl(i,j)
+ sigma_zz = lambdaplus2mu_unrelaxed_elastic*duz_dzl + lambdal_unrelaxed_elastic*PML_dux_dxl(i,j)
+ sigma_zx = mul_unrelaxed_elastic * (PML_duz_dxl(i,j) + dux_dzl)
+ sigma_xz = mul_unrelaxed_elastic * (PML_dux_dzl(i,j) + duz_dxl)
endif
else
! no attenuation
@@ -1283,7 +1283,7 @@
if(p_sv)then !P-SV waves
b_absorb_elastic_bottom(1,i,ib_bottom(ispecabs),it) = (tx + traction_x_t0)*weight
b_absorb_elastic_bottom(3,i,ib_bottom(ispecabs),it) = (tz + traction_z_t0)*weight
- else!SH (membrane) waves
+ else !SH (membrane) waves
b_absorb_elastic_bottom(2,i,ib_bottom(ispecabs),it) = ty*weight
endif
endif
@@ -1394,7 +1394,7 @@
if(p_sv)then !P-SV waves
b_absorb_elastic_top(1,i,ib_top(ispecabs),it) = (tx- traction_x_t0)*weight
b_absorb_elastic_top(3,i,ib_top(ispecabs),it) = (tz- traction_z_t0)*weight
- else!SH (membrane) waves
+ else !SH (membrane) waves
b_absorb_elastic_top(2,i,ib_top(ispecabs),it) = ty*weight
endif
endif
@@ -1403,7 +1403,7 @@
if(p_sv)then !P-SV waves
accel_elastic(1,iglob) = accel_elastic(1,iglob) - b_absorb_elastic_top(1,i,ib_top(ispecabs),NSTEP-it+1)
accel_elastic(3,iglob) = accel_elastic(3,iglob) - b_absorb_elastic_top(3,i,ib_top(ispecabs),NSTEP-it+1)
- else!SH (membrane) waves
+ else !SH (membrane) waves
accel_elastic(2,iglob) = accel_elastic(2,iglob) - b_absorb_elastic_top(2,i,ib_top(ispecabs),NSTEP-it+1)
endif
endif
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90 2013-10-02 22:20:57 UTC (rev 22916)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90 2013-10-03 00:20:11 UTC (rev 22917)
@@ -460,7 +460,7 @@
! PML fixed parameters to compute parameter in PML
double precision, parameter :: NPOWER = 2.d0
double precision, parameter :: Rcoef = 0.001d0
- double precision, parameter :: damping_modified_factor = 1.2d0
+ double precision, parameter :: damping_modified_factor = 1.0d0
double precision, parameter :: K_MAX_PML = 1.d0 ! from Gedney page 8.11
double precision :: ALPHA_MAX_PML
More information about the CIG-COMMITS
mailing list