[cig-commits] [commit] devel: Modif coupling (99d8198)

[cig_noreply at geodynamics.org]

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

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

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

commit 99d8198dbc061608db76aa2247affa1081a8f81f
Author: Clément Durochat <c.durochat at gmail.com>
Date:   Fri Oct 17 17:57:44 2014 +0200

    Modif coupling


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

99d8198dbc061608db76aa2247affa1081a8f81f
 src/generate_databases/get_model.f90             | 23 -----------------------
 src/generate_databases/model_external_values.f90 | 21 +++++++++++++++++++++
 2 files changed, 21 insertions(+), 23 deletions(-)

diff --git a/src/generate_databases/get_model.f90 b/src/generate_databases/get_model.f90
index cbbeb1b..f5b1150 100644
--- a/src/generate_databases/get_model.f90
+++ b/src/generate_databases/get_model.f90
@@ -593,26 +593,3 @@
   end select
 
   end subroutine get_model_binaries
-
-!----------------------------------------------------------------
-
-  subroutine FindLayer(x,y,z)
-    use external_model
-    implicit none
-    integer il
-    double precision radius
-    double precision :: x,y,z
-    radius =  dsqrt(x**2 + y**2 + (z+zref)**2) / 1000.d0
-
-    !write(124,*) 'RADIUS ',radius,x,y,z,z+zref,zref
-    il = 1
-    do while (radius .gt. zlayer(il).and.il.lt.nlayer)
-       il = il + 1
-    end do
-    il = il - 1
-    ilayer = il 
-    
-    !write(124,*) 'r, i : ',z,zref,radius, ilayer
-  end subroutine FindLayer
-
-!----------------------------------------------------------------
diff --git a/src/generate_databases/model_external_values.f90 b/src/generate_databases/model_external_values.f90
index 5cf6482..493f8d1 100644
--- a/src/generate_databases/model_external_values.f90
+++ b/src/generate_databases/model_external_values.f90
@@ -284,6 +284,27 @@
 !! !! ================= VM VM CUSTOM SUBROUTINE FOR DSM COUPLING
 !----------------------------------------------------------------
 
+  subroutine FindLayer(x,y,z)
+
+    use external_model
+    implicit none
+    integer il
+    double precision radius
+    double precision :: x,y,z
+    radius =  dsqrt(x**2 + y**2 + (z+zref)**2) / 1000.d0
+
+    !write(124,*) 'RADIUS ',radius,x,y,z,z+zref,zref
+    il = 1
+    do while (radius .gt. zlayer(il).and.il.lt.nlayer)
+       il = il + 1
+    end do
+    il = il - 1
+    ilayer = il 
+    
+    !write(124,*) 'r, i : ',z,zref,radius, ilayer
+  end subroutine FindLayer
+
+!----------------------------------------------------------------
 
  subroutine model_1D(x_eval,y_eval,z_eval, &
                              rho_final,vp_final,vs_final,r1)