[cig-commits] [commit] devel: Last modif to validate new coupling with DSM (b4e3c80)

[cig_noreply at geodynamics.org]

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

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d/compare/0479069f470e39347345ab3d4aff453320ac453e...57f2934cebe2e3d82d758afd2da6bf77d96c901d

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

commit b4e3c8083382e1ce972098449f98f64a916642de
Author: Clément Durochat <c.durochat at gmail.com>
Date:   Tue Oct 21 15:03:08 2014 +0200

    Last modif to validate new coupling with DSM


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

b4e3c8083382e1ce972098449f98f64a916642de
 src/generate_databases/get_model.f90 | 19 ++++++++++---------
 1 file changed, 10 insertions(+), 9 deletions(-)

diff --git a/src/generate_databases/get_model.f90 b/src/generate_databases/get_model.f90
index f5b1150..d727059 100644
--- a/src/generate_databases/get_model.f90
+++ b/src/generate_databases/get_model.f90
@@ -111,6 +111,15 @@
   ! material properties on all GLL points: taken from material values defined for
   ! each spectral element in input mesh
   do ispec = 1, nspec
+
+    if (COUPLE_WITH_EXTERNAL_CODE) then
+      iglob = ibool(3,3,3,ispec)
+      xmesh = xstore_dummy(iglob)
+      ymesh = ystore_dummy(iglob)
+      zmesh = zstore_dummy(iglob)
+      call FindLayer(xmesh,ymesh,zmesh)
+    endif
+
     ! loops over all gll points in element
     do k = 1, NGLLZ
       do j = 1, NGLLY
@@ -170,15 +179,7 @@
  
           !! VM VM for coupling with DSM
           !! find the # layer where the middle of the element is located
-          if (COUPLE_WITH_EXTERNAL_CODE) then 
-
-            if( (NGLLX == 5) .and. (NGLLY == 5) .and. (NGLLZ == 5) ) then
-              if (i==3 .and. j==3 .and. k==3) call FindLayer(xmesh,ymesh,zmesh)
-            else
-              stop 'bad number of GLL points for coupling with DSM'
-            endif
-
-          end if
+          if ( COUPLE_WITH_EXTERNAL_CODE .and. (i==3 .and. j==3 .and. k==3) ) call FindLayer(xmesh,ymesh,zmesh)
 
           ! material index 1: associated material number
           ! 1 = acoustic, 2 = elastic, 3 = poroelastic, -1 = undefined tomographic