[cig-commits] [commit] : added a comment about the way I interpolate vp_final in the transition layer between the HR and MR blocks (7f92a29)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/specfem3d

On branch  : 
Link       : https://github.com/geodynamics/specfem2d/compare/1e201257d91c794056b990a43329e05d04f77454...0000000000000000000000000000000000000000

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

commit 7f92a29a94e77a032934b9130c2c51423ee6ea9c
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Tue Dec 29 23:03:41 2009 +0000

    added a comment about the way I interpolate vp_final in the transition layer between the HR and MR blocks


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

7f92a29a94e77a032934b9130c2c51423ee6ea9c
 interpolate_gocad_block_HR.f90 | 12 ++++++++++++
 1 file changed, 12 insertions(+)

diff --git a/interpolate_gocad_block_HR.f90 b/interpolate_gocad_block_HR.f90
index 871ed09..e4eaaf4 100644
--- a/interpolate_gocad_block_HR.f90
+++ b/interpolate_gocad_block_HR.f90
@@ -106,6 +106,9 @@
               VP_MIN_GOCAD,VP_VS_RATIO_GOCAD_TOP,VP_VS_RATIO_GOCAD_BOTTOM, &
               IMPOSE_MINIMUM_VP_GOCAD,THICKNESS_TAPER_BLOCK_HR, &
               vp_hauksson,vs_hauksson,doubling_index,HAUKSSON_REGIONAL_MODEL,MOHO_MAP_LUPEI)
+! the vp_final we start from is the one we computed by interpolation in the HR block,
+! and as we get closer to the transition with the MR block we mix it with the value on the edge of MR;
+! it is not a linear interpolation between two fixed values but rather a linear variation of the amount of mixing
     vp_final = vp_ref_MR * (1. - gamma_interp_x) + vp_final * gamma_interp_x
 
 ! x = xmax
@@ -116,6 +119,9 @@
               VP_MIN_GOCAD,VP_VS_RATIO_GOCAD_TOP,VP_VS_RATIO_GOCAD_BOTTOM, &
               IMPOSE_MINIMUM_VP_GOCAD,THICKNESS_TAPER_BLOCK_HR, &
               vp_hauksson,vs_hauksson,doubling_index,HAUKSSON_REGIONAL_MODEL, MOHO_MAP_LUPEI)
+! the vp_final we start from is the one we computed by interpolation in the HR block,
+! and as we get closer to the transition with the MR block we mix it with the value on the edge of MR;
+! it is not a linear interpolation between two fixed values but rather a linear variation of the amount of mixing
     vp_final = vp_ref_MR * gamma_interp_x + vp_final * (1. - gamma_interp_x)
 
 ! y = ymin
@@ -126,6 +132,9 @@
               VP_MIN_GOCAD,VP_VS_RATIO_GOCAD_TOP,VP_VS_RATIO_GOCAD_BOTTOM, &
               IMPOSE_MINIMUM_VP_GOCAD,THICKNESS_TAPER_BLOCK_HR, &
               vp_hauksson,vs_hauksson,doubling_index,HAUKSSON_REGIONAL_MODEL, MOHO_MAP_LUPEI)
+! the vp_final we start from is the one we computed by interpolation in the HR block,
+! and as we get closer to the transition with the MR block we mix it with the value on the edge of MR;
+! it is not a linear interpolation between two fixed values but rather a linear variation of the amount of mixing
     vp_final = vp_ref_MR * (1. - gamma_interp_y) + vp_final * gamma_interp_y
 
 ! y = ymax
@@ -136,6 +145,9 @@
               VP_MIN_GOCAD,VP_VS_RATIO_GOCAD_TOP,VP_VS_RATIO_GOCAD_BOTTOM, &
               IMPOSE_MINIMUM_VP_GOCAD,THICKNESS_TAPER_BLOCK_HR, &
               vp_hauksson,vs_hauksson,doubling_index,HAUKSSON_REGIONAL_MODEL, MOHO_MAP_LUPEI)
+! the vp_final we start from is the one we computed by interpolation in the HR block,
+! and as we get closer to the transition with the MR block we mix it with the value on the edge of MR;
+! it is not a linear interpolation between two fixed values but rather a linear variation of the amount of mixing
     vp_final = vp_ref_MR * gamma_interp_y + vp_final * (1. - gamma_interp_y)
 
   endif