[cig-commits] r20828 - seismo/2D/SPECFEM2D/trunk/src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Fri Oct 12 08:41:08 PDT 2012
Author: dkomati1
Date: 2012-10-12 08:41:07 -0700 (Fri, 12 Oct 2012)
New Revision: 20828
Modified:
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
reduced the length of variable name rmemory_potential_acoustic_LDDRK because it was longer than the standard
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90 2012-10-12 15:04:20 UTC (rev 20827)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90 2012-10-12 15:41:07 UTC (rev 20828)
@@ -59,7 +59,7 @@
K_x_store,K_z_store,d_x_store,d_z_store,alpha_x_store,alpha_z_store,&
rmemory_potential_acoustic,&
rmemory_acoustic_dux_dx,rmemory_acoustic_dux_dz,&
- rmemory_potential_acoustic_LDDRK,alpha_LDDRK,beta_LDDRK, &
+ rmemory_potential_acoust_LDDRK,alpha_LDDRK,beta_LDDRK, &
rmemory_acoustic_dux_dx_LDDRK,rmemory_acoustic_dux_dz_LDDRK,&
deltat,PML_BOUNDARY_CONDITIONS)
@@ -151,7 +151,7 @@
!coefficients and memory variables when using CPML with LDDRK
integer :: stage_time_scheme,i_stage
real(kind=CUSTOM_REAL), dimension(Nstages) :: alpha_LDDRK,beta_LDDRK
- real(kind=CUSTOM_REAL), dimension(2,NGLLX,NGLLZ,nspec_PML) :: rmemory_potential_acoustic_LDDRK
+ real(kind=CUSTOM_REAL), dimension(2,NGLLX,NGLLZ,nspec_PML) :: rmemory_potential_acoust_LDDRK
real(kind=CUSTOM_REAL), dimension(NGLLX,NGLLZ,nspec_PML) :: &
rmemory_acoustic_dux_dx_LDDRK,rmemory_acoustic_dux_dz_LDDRK
@@ -571,12 +571,12 @@
if(stage_time_scheme == 6) then
bb = alpha_x_store(i,j,ispec_PML)
- rmemory_potential_acoustic_LDDRK(1,i,j,ispec_PML) = &
- alpha_LDDRK(i_stage) * rmemory_potential_acoustic_LDDRK(1,i,j,ispec_PML) &
+ rmemory_potential_acoust_LDDRK(1,i,j,ispec_PML) = &
+ alpha_LDDRK(i_stage) * rmemory_potential_acoust_LDDRK(1,i,j,ispec_PML) &
+ deltat * (-bb * rmemory_potential_acoustic(1,i,j,ispec_PML) + potential_acoustic(iglob))
rmemory_potential_acoustic(1,i,j,ispec_PML) = rmemory_potential_acoustic(1,i,j,ispec_PML) + &
- beta_LDDRK(i_stage) * rmemory_potential_acoustic_LDDRK(1,i,j,ispec_PML)
+ beta_LDDRK(i_stage) * rmemory_potential_acoust_LDDRK(1,i,j,ispec_PML)
rmemory_potential_acoustic(2,i,j,ispec_PML) =0.d0
end if
@@ -623,18 +623,18 @@
stop "something wrong in alpha definition"
endif
- rmemory_potential_acoustic_LDDRK(1,i,j,ispec_PML) = &
- alpha_LDDRK(i_stage) * rmemory_potential_acoustic_LDDRK(1,i,j,ispec_PML) &
+ rmemory_potential_acoust_LDDRK(1,i,j,ispec_PML) = &
+ alpha_LDDRK(i_stage) * rmemory_potential_acoust_LDDRK(1,i,j,ispec_PML) &
+ deltat * (-bb * rmemory_potential_acoustic(1,i,j,ispec_PML) + potential_acoustic(iglob))
- rmemory_potential_acoustic_LDDRK(2,i,j,ispec_PML) = &
- alpha_LDDRK(i_stage) * rmemory_potential_acoustic_LDDRK(2,i,j,ispec_PML) &
+ rmemory_potential_acoust_LDDRK(2,i,j,ispec_PML) = &
+ alpha_LDDRK(i_stage) * rmemory_potential_acoust_LDDRK(2,i,j,ispec_PML) &
+ deltat * (-bb * rmemory_potential_acoustic(2,i,j,ispec_PML) &
+ rmemory_potential_acoustic(1,i,j,ispec_PML))
rmemory_potential_acoustic(1,i,j,ispec_PML) = rmemory_potential_acoustic(1,i,j,ispec_PML) + &
- beta_LDDRK(i_stage) * rmemory_potential_acoustic_LDDRK(1,i,j,ispec_PML)
+ beta_LDDRK(i_stage) * rmemory_potential_acoust_LDDRK(1,i,j,ispec_PML)
rmemory_potential_acoustic(2,i,j,ispec_PML) = rmemory_potential_acoustic(2,i,j,ispec_PML) + &
- beta_LDDRK(i_stage) * rmemory_potential_acoustic_LDDRK(2,i,j,ispec_PML)
+ beta_LDDRK(i_stage) * rmemory_potential_acoust_LDDRK(2,i,j,ispec_PML)
end if
@@ -660,12 +660,12 @@
rmemory_potential_acoustic(1,i,j,ispec_PML) =0.d0
- rmemory_potential_acoustic_LDDRK(2,i,j,ispec_PML) = &
- alpha_LDDRK(i_stage) * rmemory_potential_acoustic_LDDRK(2,i,j,ispec_PML) &
+ rmemory_potential_acoust_LDDRK(2,i,j,ispec_PML) = &
+ alpha_LDDRK(i_stage) * rmemory_potential_acoust_LDDRK(2,i,j,ispec_PML) &
+ deltat * (-bb * rmemory_potential_acoustic(2,i,j,ispec_PML) + potential_acoustic(iglob))
rmemory_potential_acoustic(2,i,j,ispec_PML) = rmemory_potential_acoustic(2,i,j,ispec_PML) + &
- beta_LDDRK(i_stage) * rmemory_potential_acoustic_LDDRK(2,i,j,ispec_PML)
+ beta_LDDRK(i_stage) * rmemory_potential_acoust_LDDRK(2,i,j,ispec_PML)
end if
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-10-12 15:04:20 UTC (rev 20827)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-10-12 15:41:07 UTC (rev 20828)
@@ -1020,7 +1020,7 @@
real(kind=CUSTOM_REAL), dimension(:,:,:), allocatable :: &
rmemory_acoustic_dux_dx,rmemory_acoustic_dux_dz
- real(kind=CUSTOM_REAL), dimension(:,:,:,:), allocatable :: rmemory_potential_acoustic_LDDRK
+ real(kind=CUSTOM_REAL), dimension(:,:,:,:), allocatable :: rmemory_potential_acoust_LDDRK
real(kind=CUSTOM_REAL), dimension(:,:,:), allocatable :: &
rmemory_acoustic_dux_dx_LDDRK,rmemory_acoustic_dux_dz_LDDRK
@@ -2993,14 +2993,14 @@
if(time_stepping_scheme == 2)then
- allocate(rmemory_potential_acoustic_LDDRK(2,NGLLX,NGLLZ,nspec_PML),stat=ier)
+ allocate(rmemory_potential_acoust_LDDRK(2,NGLLX,NGLLZ,nspec_PML),stat=ier)
if(ier /= 0) stop 'error: not enough memory to allocate array rmemory_potential_acoustic'
allocate(rmemory_acoustic_dux_dx_LDDRK(NGLLX,NGLLZ,nspec_PML),stat=ier)
if(ier /= 0) stop 'error: not enough memory to allocate array rmemory_acoustic_dux_dx'
allocate(rmemory_acoustic_dux_dz_LDDRK(NGLLX,NGLLZ,nspec_PML),stat=ier)
if(ier /= 0) stop 'error: not enough memory to allocate array rmemory_acoustic_dux_dz'
- rmemory_potential_acoustic_LDDRK = ZERO
+ rmemory_potential_acoust_LDDRK = ZERO
rmemory_acoustic_dux_dx_LDDRK = ZERO
rmemory_acoustic_dux_dz_LDDRK = ZERO
endif
@@ -3034,7 +3034,7 @@
allocate(rmemory_acoustic_dux_dx(1,1,1))
allocate(rmemory_acoustic_dux_dz(1,1,1))
- allocate(rmemory_potential_acoustic_LDDRK(1,1,1,1))
+ allocate(rmemory_potential_acoust_LDDRK(1,1,1,1))
allocate(rmemory_acoustic_dux_dx_LDDRK(1,1,1))
allocate(rmemory_acoustic_dux_dz_LDDRK(1,1,1))
@@ -4998,7 +4998,7 @@
K_x_store,K_z_store,d_x_store,d_z_store,alpha_x_store,alpha_z_store,&
rmemory_potential_acoustic,&
rmemory_acoustic_dux_dx,rmemory_acoustic_dux_dz,&
- rmemory_potential_acoustic_LDDRK,alpha_LDDRK,beta_LDDRK, &
+ rmemory_potential_acoust_LDDRK,alpha_LDDRK,beta_LDDRK, &
rmemory_acoustic_dux_dx_LDDRK,rmemory_acoustic_dux_dz_LDDRK,&
deltat,PML_BOUNDARY_CONDITIONS)
if( SIMULATION_TYPE == 2 ) then
@@ -5019,7 +5019,7 @@
K_x_store,K_z_store,d_x_store,d_z_store,alpha_x_store,alpha_z_store,&
rmemory_potential_acoustic,&
rmemory_acoustic_dux_dx,rmemory_acoustic_dux_dz,&
- rmemory_potential_acoustic_LDDRK,alpha_LDDRK,beta_LDDRK, &
+ rmemory_potential_acoust_LDDRK,alpha_LDDRK,beta_LDDRK, &
rmemory_acoustic_dux_dx_LDDRK,rmemory_acoustic_dux_dz_LDDRK,&
deltat,PML_BOUNDARY_CONDITIONS)
endif
More information about the CIG-COMMITS
mailing list