[cig-commits] [commit] devel: Change in mesh-chunk subroutines (c497ce0)

[cig_noreply at geodynamics.org]

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

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d/compare/ae876feb0b07b860f769cd822d71da0c57556416...c497ce046eaca64c0b32596123eee513e0ac361c

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

commit c497ce046eaca64c0b32596123eee513e0ac361c
Author: Clément Durochat <c.durochat at gmail.com>
Date:   Tue Oct 28 17:07:48 2014 +0100

    Change in mesh-chunk subroutines


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

c497ce046eaca64c0b32596123eee513e0ac361c
 src/meshfem3D/earth_chunk_HEX8_ReadIasp91.f90 | 4 ++--
 src/meshfem3D/meshfem3D.f90                   | 5 +++--
 2 files changed, 5 insertions(+), 4 deletions(-)

diff --git a/src/meshfem3D/earth_chunk_HEX8_ReadIasp91.f90 b/src/meshfem3D/earth_chunk_HEX8_ReadIasp91.f90
index a1fbd37..68d8d7d 100644
--- a/src/meshfem3D/earth_chunk_HEX8_ReadIasp91.f90
+++ b/src/meshfem3D/earth_chunk_HEX8_ReadIasp91.f90
@@ -238,7 +238,7 @@ end subroutine StorePointZ
    zpoint(1)=zlayer(nlayer) - Z_DEPTH_BLOCK
    write(*,*) zlayer(nlayer) ,  Z_DEPTH_BLOCK
    !! niveau de depart
-   call FindLayer_in_mesh_chunk(ilayer,zlayer,zpoint(1),nlayer)
+   call FindLayer_for_earth_chunk_mesh(ilayer,zlayer,zpoint(1),nlayer)
    write(*,*) '              INITIALISATION calcul du niveau de depart : '
    write(*,*)
    write(*,*) 'zlayer : ', zlayer
@@ -342,7 +342,7 @@ end subroutine StorePointZ
       ProfForGemini(ilay-1,2)  =  zz(ilay+1)
       ProfForGemini(ilay-1,3)  = 0.5d0 * (zz(ilay) + zz(ilay+1))
 
-      call FindLayer_in_mesh_chunk(niveau,zlayer, ProfForGemini(ilay-1,3),nlayer)
+      call FindLayer_for_earth_chunk_mesh(niveau,zlayer, ProfForGemini(ilay-1,3),nlayer)
       Niveau_elm(ilay-1)=niveau
       write(*,'(i5,2f15.3,i10)') ilay,zz(ilay),zz(ilay+1),niveau
    enddo
diff --git a/src/meshfem3D/meshfem3D.f90 b/src/meshfem3D/meshfem3D.f90
index b16669d..ec2d9f0 100644
--- a/src/meshfem3D/meshfem3D.f90
+++ b/src/meshfem3D/meshfem3D.f90
@@ -420,9 +420,10 @@
 
     elseif (NGNOD == 27) then
       ! creates mesh in MESH/
-      call earth_chunk_HEX27_Mesher()
+!! to uncomment when finished ==> call earth_chunk_HEX27_Mesher()
       ! done with mesher
-      stop 'Done creating a chunk of the earth Mesh (HEX27 elements), see directory MESH/'
+      stop 'Creating a chunk of the earth Mesh not implmented yet with HEX27 elements'
+!! to uncomment when finished ==> stop 'Done creating a chunk of the earth Mesh (HEX27 elements), see directory MESH/'
 
     else
       stop 'Bad number of nodes per hexahedron : NGNOD must be equal to 8 or 27'