[cig-commits] r20490 - seismo/2D/SPECFEM2D/trunk/src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Sat Jul 7 08:57:53 PDT 2012
Author: dkomati1
Date: 2012-07-07 08:57:52 -0700 (Sat, 07 Jul 2012)
New Revision: 20490
Modified:
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/setup_sources_receivers.F90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
suppressed warnings about variables that were too long
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90 2012-07-07 15:47:01 UTC (rev 20489)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90 2012-07-07 15:57:52 UTC (rev 20490)
@@ -57,7 +57,7 @@
b_absorb_acoustic_bottom,b_absorb_acoustic_top,IS_BACKWARD_FIELD,&
is_PML,nspec_PML,npoin_PML,ibool_PML,spec_to_PML,which_PML_elem,&
K_x_store,K_z_store,d_x_store,d_z_store,alpha_x_store,alpha_z_store,&
- rmemory_potential_acoustic,rmemory_potential_acoustic_corner,&
+ rmemory_potential_acoustic,rmemory_pot_acoustic_corner,&
rmemory_acoustic_dux_dx,rmemory_acoustic_dux_dz,&
rmemory_acoustic_dux_dx_corner,rmemory_acoustic_dux_dz_corner,deltat,&
PML_BOUNDARY_CONDITIONS)
@@ -134,14 +134,14 @@
integer, dimension(nspec) :: spec_to_PML
integer, dimension(NGLLX,NGLLZ,nspec) :: ibool_PML
- real(kind=CUSTOM_REAL), dimension(2,NGLLX,NGLLZ,nspec_PML) :: rmemory_potential_acoustic,rmemory_potential_acoustic_corner
+ real(kind=CUSTOM_REAL), dimension(2,NGLLX,NGLLZ,nspec_PML) :: rmemory_potential_acoustic,rmemory_pot_acoustic_corner
real(kind=CUSTOM_REAL), dimension(2,NGLLX,NGLLZ,nspec_PML) :: &
rmemory_acoustic_dux_dx,rmemory_acoustic_dux_dz,&
rmemory_acoustic_dux_dx_corner,rmemory_acoustic_dux_dz_corner
real(kind=CUSTOM_REAL), dimension(npoin_PML) :: K_x_store,K_z_store,d_x_store,d_z_store,alpha_x_store,alpha_z_store
real(kind=CUSTOM_REAL), dimension(NGLLX,NGLLZ,nspec_PML) :: potential_dot_dot_acoustic_PML, &
- potential_dot_dot_acoustic_PML_corner
+ pot_dotdot_acoustic_PML_corner
real(kind=CUSTOM_REAL), dimension(NGLLX,NGLLZ,nspec_PML) ::PML_dux_dxl,PML_dux_dzl,&
PML_dux_dxl_new,PML_dux_dzl_new
real(kind=CUSTOM_REAL) :: coef0_x, coef1_x, coef2_x, coef0_z, coef1_z, coef2_z,bb,aa,&
@@ -515,8 +515,8 @@
coef1_x = (1 - exp(- bb * deltat / 2)) * aa
coef2_x = (1 - exp(- bb * deltat / 2)) * exp(- bb * deltat / 2) * aa
- rmemory_potential_acoustic_corner(1,i,j,ispec_PML)=&
- coef0_x * rmemory_potential_acoustic_corner(1,i,j,ispec_PML) &
+ rmemory_pot_acoustic_corner(1,i,j,ispec_PML)=&
+ coef0_x * rmemory_pot_acoustic_corner(1,i,j,ispec_PML) &
+ (potential_acoustic(iglob)+deltat*potential_dot_acoustic(iglob)) * coef1_x &
+ potential_acoustic(iglob) * coef2_x
else
@@ -527,8 +527,8 @@
coef1_x = (1 - exp(- bb * deltat / 2)) * aa
coef2_x = (1 - exp(- bb * deltat / 2)) * exp(- bb * deltat / 2) * aa
- rmemory_potential_acoustic_corner(1,i,j,ispec_PML)=&
- coef0_x * rmemory_potential_acoustic_corner(1,i,j,ispec_PML) &
+ rmemory_pot_acoustic_corner(1,i,j,ispec_PML)=&
+ coef0_x * rmemory_pot_acoustic_corner(1,i,j,ispec_PML) &
+ (potential_acoustic(iglob)+deltat*potential_dot_acoustic(iglob)) * coef1_x &
+ potential_acoustic(iglob) * coef2_x
@@ -541,8 +541,8 @@
coef1_x = (1 - exp(- bb * deltat / 2)) * aa
coef2_x = (1 - exp(- bb * deltat / 2)) * exp(- bb * deltat / 2) * aa
- rmemory_potential_acoustic_corner(2,i,j,ispec_PML)=&
- coef0_x * rmemory_potential_acoustic_corner(2,i,j,ispec_PML) &
+ rmemory_pot_acoustic_corner(2,i,j,ispec_PML)=&
+ coef0_x * rmemory_pot_acoustic_corner(2,i,j,ispec_PML) &
+ (potential_acoustic(iglob)+deltat*potential_dot_acoustic(iglob)) * coef1_x &
+ potential_acoustic(iglob) * coef2_x
@@ -554,8 +554,8 @@
coef1_x = (1 - exp(- bb * deltat / 2)) * aa
coef2_x = (1 - exp(- bb * deltat / 2)) * exp(- bb * deltat / 2) * aa
- rmemory_potential_acoustic_corner(2,i,j,ispec_PML)=&
- coef0_x * rmemory_potential_acoustic_corner(2,i,j,ispec_PML) &
+ rmemory_pot_acoustic_corner(2,i,j,ispec_PML)=&
+ coef0_x * rmemory_pot_acoustic_corner(2,i,j,ispec_PML) &
+ (potential_acoustic(iglob)+deltat*potential_dot_acoustic(iglob)) * coef1_x &
+ potential_acoustic(iglob) * coef2_x
endif
@@ -597,9 +597,9 @@
ispec_PML=spec_to_PML(ispec)
iPML=ibool_PML(i,j,ispec)
- potential_dot_dot_acoustic_PML_corner(i,j,ispec_PML)= wxgll(i)*wzgll(j)/ kappal*jacobian(i,j,ispec) * &
+ pot_dotdot_acoustic_PML_corner(i,j,ispec_PML)= wxgll(i)*wzgll(j)/ kappal*jacobian(i,j,ispec) * &
((d_x_store(iPML)*k_z_store(iPML)+d_z_store(iPML)*k_x_store(iPML))*potential_dot_acoustic(iglob)+ &
- rmemory_potential_acoustic_corner(1,i,j,ispec_PML)+rmemory_potential_acoustic_corner(2,i,j,ispec_PML)&
+ rmemory_pot_acoustic_corner(1,i,j,ispec_PML)+rmemory_pot_acoustic_corner(2,i,j,ispec_PML)&
-k_z_store(iPML)*d_x_store(iPML)*alpha_x_store(iPML)*potential_acoustic(iglob)-&
k_x_store(iPML)*d_z_store(iPML)*alpha_z_store(iPML)*potential_acoustic(iglob)+&
d_z_store(iPML)*d_x_store(iPML)*potential_acoustic(iglob))
@@ -651,7 +651,7 @@
+ potential_dot_dot_acoustic_PML(i,j,ispec_PML)
potential_dot_dot_acoustic(iglob) = potential_dot_dot_acoustic(iglob) &
- - potential_dot_dot_acoustic_PML_corner(i,j,ispec_PML)
+ - pot_dotdot_acoustic_PML_corner(i,j,ispec_PML)
endif
else
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/setup_sources_receivers.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/setup_sources_receivers.F90 2012-07-07 15:47:01 UTC (rev 20489)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/setup_sources_receivers.F90 2012-07-07 15:57:52 UTC (rev 20490)
@@ -93,7 +93,10 @@
integer :: ixmin, ixmax, izmin, izmax
! Local variables
- integer irec,i_source,ispec,ispec_acoustic_surface
+ integer i_source,ispec,ispec_acoustic_surface
+#ifndef USE_MPI
+ integer irec
+#endif
do i_source=1,NSOURCES
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-07-07 15:47:01 UTC (rev 20489)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-07-07 15:57:52 UTC (rev 20490)
@@ -1004,7 +1004,7 @@
real(kind=CUSTOM_REAL), dimension(:,:,:,:,:), allocatable :: rmemory_displ_elastic_corner
real(kind=CUSTOM_REAL), dimension(:,:,:,:), allocatable :: rmemory_potential_acoustic,&
- rmemory_potential_acoustic_corner,&
+ rmemory_pot_acoustic_corner,&
rmemory_acoustic_dux_dx,rmemory_acoustic_dux_dz,&
rmemory_acoustic_dux_dx_corner,rmemory_acoustic_dux_dz_corner
@@ -2910,7 +2910,7 @@
if (any_acoustic .and. nspec_PML>0) then
allocate(rmemory_potential_acoustic(2,NGLLX,NGLLZ,nspec_PML))
- allocate(rmemory_potential_acoustic_corner(2,NGLLX,NGLLZ,nspec_PML))
+ allocate(rmemory_pot_acoustic_corner(2,NGLLX,NGLLZ,nspec_PML))
allocate(rmemory_acoustic_dux_dx(2,NGLLX,NGLLZ,nspec_PML))
allocate(rmemory_acoustic_dux_dz(2,NGLLX,NGLLZ,nspec_PML))
@@ -2921,7 +2921,7 @@
else
allocate(rmemory_potential_acoustic(1,1,1,1))
- allocate(rmemory_potential_acoustic_corner(1,1,1,1))
+ allocate(rmemory_pot_acoustic_corner(1,1,1,1))
allocate(rmemory_acoustic_dux_dx(1,1,1,1))
allocate(rmemory_acoustic_dux_dz(1,1,1,1))
@@ -2950,7 +2950,7 @@
allocate(rmemory_displ_elastic_corner(1,1,1,1,1))
allocate(rmemory_potential_acoustic(1,1,1,1))
- allocate(rmemory_potential_acoustic_corner(1,1,1,1))
+ allocate(rmemory_pot_acoustic_corner(1,1,1,1))
allocate(rmemory_acoustic_dux_dx(1,1,1,1))
allocate(rmemory_acoustic_dux_dz(1,1,1,1))
@@ -4917,7 +4917,7 @@
b_absorb_acoustic_bottom,b_absorb_acoustic_top,.false.,&
is_PML,nspec_PML,npoin_PML,ibool_PML,spec_to_PML,which_PML_elem,&
K_x_store,K_z_store,d_x_store,d_z_store,alpha_x_store,alpha_z_store,&
- rmemory_potential_acoustic,rmemory_potential_acoustic_corner,&
+ rmemory_potential_acoustic,rmemory_pot_acoustic_corner,&
rmemory_acoustic_dux_dx,rmemory_acoustic_dux_dz,&
rmemory_acoustic_dux_dx_corner,rmemory_acoustic_dux_dz_corner,deltat,&
PML_BOUNDARY_CONDITIONS)
@@ -4937,7 +4937,7 @@
b_absorb_acoustic_bottom,b_absorb_acoustic_top,.true.,&
is_PML,nspec_PML,npoin_PML,ibool_PML,spec_to_PML,which_PML_elem,&
K_x_store,K_z_store,d_x_store,d_z_store,alpha_x_store,alpha_z_store,&
- rmemory_potential_acoustic,rmemory_potential_acoustic_corner,&
+ rmemory_potential_acoustic,rmemory_pot_acoustic_corner,&
rmemory_acoustic_dux_dx,rmemory_acoustic_dux_dz,&
rmemory_acoustic_dux_dx_corner,rmemory_acoustic_dux_dz_corner,deltat,&
PML_BOUNDARY_CONDITIONS)
More information about the CIG-COMMITS
mailing list