[cig-commits] [commit] devel: Quentin Brissaud fixed two minor bugs in his gravity routines (bd53bae)

[cig_noreply at geodynamics.org]

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