[cig-commits] [commit] devel: fixed an important bug in the call to compute_pressure_whole_medium() (cf40549)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/specfem2d

On branch  : devel
Link       : https://github.com/geodynamics/specfem2d/compare/065f6d28ff5e8a1ab64a09a912c2510e9a10c616...cf40549f4c1d258bf2cfb1af1cce7f53f4661e58

>---------------------------------------------------------------

commit cf40549f4c1d258bf2cfb1af1cce7f53f4661e58
Author: Dimitri Komatitsch <komatits at users.noreply.github.com>
Date:   Wed Feb 26 02:09:15 2014 +0100

    fixed an important bug in the call to compute_pressure_whole_medium()
    
    potential_acoustic was used instead of potential_dot_dot_acoustic when dumping the pressure field to disk (but JPEG pictures were fine).


>---------------------------------------------------------------

cf40549f4c1d258bf2cfb1af1cce7f53f4661e58
 src/specfem2D/specfem2D.F90 | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/src/specfem2D/specfem2D.F90 b/src/specfem2D/specfem2D.F90
index 19ee03b..0b1e2f4 100644
--- a/src/specfem2D/specfem2D.F90
+++ b/src/specfem2D/specfem2D.F90
@@ -1,4 +1,3 @@
-
   program specfem2D
 
 !========================================================================
@@ -8563,7 +8562,7 @@ if(coupled_elastic_poro) then
         else if(imagetype_wavefield_dumps == 4 .and. p_sv) then
 
           if (myrank == 0) write(IOUT,*) 'dumping the pressure field...'
-          call compute_pressure_whole_medium(potential_acoustic,displ_elastic,&
+          call compute_pressure_whole_medium(potential_dot_dot_acoustic,displ_elastic,&
                      displs_poroelastic,displw_poroelastic,elastic,poroelastic,vector_field_display, &
                      xix,xiz,gammax,gammaz,ibool,hprime_xx,hprime_zz,nspec, &
                      nglob,nglob_acoustic,nglob_elastic,nglob_poroelastic,assign_external_model, &