[cig-commits] [commit] devel: fixed reading external models for SH simulations (3e6e20c)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem2d

On branch  : devel
Link       : https://github.com/geodynamics/specfem2d/compare/7a7ca08afb0578095ac833a886583b5fdee17ee3...caf697fd40a47a7cef1cbc80e13c600c3ec41167

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

commit 3e6e20cd0d510ce9a38ace525a0b0911c9cd3497
Author: rmodrak <rmodrak at princeton.edu>
Date:   Thu Nov 13 16:51:24 2014 -0500

    fixed reading external models for SH simulations


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

3e6e20cd0d510ce9a38ace525a0b0911c9cd3497
 src/specfem2D/read_external_model.f90 | 8 +++++---
 src/specfem2D/specfem2D.F90           | 2 +-
 2 files changed, 6 insertions(+), 4 deletions(-)

diff --git a/src/specfem2D/read_external_model.f90 b/src/specfem2D/read_external_model.f90
index f89d0dc..f479c14 100644
--- a/src/specfem2D/read_external_model.f90
+++ b/src/specfem2D/read_external_model.f90
@@ -50,7 +50,7 @@
                 coord,kmato,rhoext,vpext,vsext,gravityext,Nsqext, &
                 QKappa_attenuationext,Qmu_attenuationext, &
                 c11ext,c13ext,c15ext,c33ext,c35ext,c55ext,c12ext,c23ext,c25ext, &
-                READ_EXTERNAL_SEP_FILE,ATTENUATION_VISCOELASTIC_SOLID)
+                READ_EXTERNAL_SEP_FILE,ATTENUATION_VISCOELASTIC_SOLID,p_sv)
 
   implicit none
   include "constants.h"
@@ -64,7 +64,7 @@
   double precision, dimension(NDIM,nglob) :: coord
 
   ! Material properties
-  logical :: any_acoustic,any_gravitoacoustic,any_elastic,any_poroelastic,READ_EXTERNAL_SEP_FILE,ATTENUATION_VISCOELASTIC_SOLID
+  logical :: any_acoustic,any_gravitoacoustic,any_elastic,any_poroelastic,READ_EXTERNAL_SEP_FILE,ATTENUATION_VISCOELASTIC_SOLID,p_sv
   integer, dimension(nspec) :: kmato
   logical, dimension(nspec) :: acoustic,gravitoacoustic,elastic,poroelastic
   double precision, dimension(NGLLX,NGLLZ,nspec) :: rhoext,vpext,vsext,gravityext,Nsqext
@@ -176,7 +176,9 @@
     do j = 1,NGLLZ
       do i = 1,NGLLX
         iglob = ibool(i,j,ispec)
-        if(.not. (i == 1 .and. j == 1) .and. &
+
+
+        if(p_sv .and. (.not. (i == 1 .and. j == 1)) .and. &
           ((vsext(i,j,ispec) >= TINYVAL .and. previous_vsext < TINYVAL) .or. &
            (vsext(i,j,ispec) < TINYVAL  .and. previous_vsext >= TINYVAL)))  &
           call exit_MPI('external velocity model cannot be both fluid and solid inside the same spectral element')
diff --git a/src/specfem2D/specfem2D.F90 b/src/specfem2D/specfem2D.F90
index d152f72..0081d13 100644
--- a/src/specfem2D/specfem2D.F90
+++ b/src/specfem2D/specfem2D.F90
@@ -2271,7 +2271,7 @@
                 coord,kmato,rhoext,vpext,vsext,gravityext,Nsqext, &
                 QKappa_attenuationext,Qmu_attenuationext, &
                 c11ext,c13ext,c15ext,c33ext,c35ext,c55ext,c12ext,c23ext,c25ext, &
-                READ_EXTERNAL_SEP_FILE,ATTENUATION_VISCOELASTIC_SOLID)
+                READ_EXTERNAL_SEP_FILE,ATTENUATION_VISCOELASTIC_SOLID,p_sv)
   endif
 
 !