[cig-commits] [commit] devel: Quentin Brissaud fixed two minor bugs in his gravity routines (bd53bae)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Apr 9 10:28:22 PDT 2014
Repository : ssh://geoshell/specfem2d
On branch : devel
Link : https://github.com/geodynamics/specfem2d/compare/35e0e3747dd39bd7418991a98ceb247c58167a78...80fabb71a99eadbaada38a043c9510f318d83015
>---------------------------------------------------------------
commit bd53baebf4b626fcba7c9a652d3e5898433df4d8
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Wed Apr 9 19:27:10 2014 +0200
Quentin Brissaud fixed two minor bugs in his gravity routines
>---------------------------------------------------------------
bd53baebf4b626fcba7c9a652d3e5898433df4d8
src/specfem2D/acoustic_forcing_boundary.f90 | 4 ++--
src/specfem2D/define_external_model.f90 | 2 +-
2 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/specfem2D/acoustic_forcing_boundary.f90 b/src/specfem2D/acoustic_forcing_boundary.f90
index a0181fb..2389e9e 100644
--- a/src/specfem2D/acoustic_forcing_boundary.f90
+++ b/src/specfem2D/acoustic_forcing_boundary.f90
@@ -58,6 +58,8 @@
x = coord(1,iglob)
delayed = 0
+! speed of light
+ c = 300000000.
if(forcing_type == 1) then !! First test function : same forcing for the whole boundary
! print *, ispec_acoustic
@@ -69,7 +71,6 @@
! infrasounds / seismic
accel_x = 0 !* Apo
- ! accel_z = 1000000000000
accel_z = A * (exp(-(alpha*(deltat*it-40-t0)/tho)**2) &
- exp(-(alpha*(deltat*it-70-t0)/tho)**2)) !* Apo
@@ -80,7 +81,6 @@
! (exp(-(alpha*(deltat*it+1000-t0)/tho)**2) &
! - exp(-(alpha*(deltat*it-1300-t0)/tho)**2)) !* Apo
-!WRITE(6,*) "acoustic_forcing.90 --> ACELL_Z : ", accel_z
! gravity wave /tsunami
! accel_x = 0 !* Apo
! accel_z = A * (exp(-(alpha*(deltat*it-1000-t0)/tho)**2) &
diff --git a/src/specfem2D/define_external_model.f90 b/src/specfem2D/define_external_model.f90
index 0b01ff5..d94f120 100644
--- a/src/specfem2D/define_external_model.f90
+++ b/src/specfem2D/define_external_model.f90
@@ -1221,7 +1221,7 @@
Nsq(i,j,ispec) = Nsq_atmos(ii-1) + frac * (Nsq_atmos(ii)-Nsq_atmos(ii-1))
if (Nsq(i,j,ispec) <= 0.0) then
write(*,*) 'STOP Negative Nsquare !!! :', &
- i,j,coord(1,iglob1),coord(2,iglob1),Nsq(i,j,ispec),gravity(i,j,ispec),vp(i,j,ispec)
+ i,j,coord(1,iglob),coord(2,iglob),Nsq(i,j,ispec),gravity(i,j,ispec),vp(i,j,ispec)
stop
endif
vs(i,j,ispec) = 0.d0
More information about the CIG-COMMITS
mailing list