[cig-commits] [commit] devel, master: fixed a bug in Matthias Meschede's crustal model smoothing routines: latitude was used instead of longitude (b173684)

[cig_noreply at geodynamics.org]

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

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f

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

commit b173684dc995f509d41db6f327de2f1c2d9e5654
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Fri May 2 15:00:24 2014 +0200

    fixed a bug in Matthias Meschede's crustal model smoothing routines: latitude was used instead of longitude


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

b173684dc995f509d41db6f327de2f1c2d9e5654
 src/meshfem3D/model_crustmaps.f90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/meshfem3D/model_crustmaps.f90 b/src/meshfem3D/model_crustmaps.f90
index cb81c5f..4d9dd2b 100644
--- a/src/meshfem3D/model_crustmaps.f90
+++ b/src/meshfem3D/model_crustmaps.f90
@@ -370,7 +370,7 @@
 
   ! checks if inside/outside of critical region for mesh stretching
   if( SMOOTH_CRUST ) then
-    dist = dsqrt( (lon-LAT_CRITICAL_EUROPE)**2 + (lat-LAT_CRITICAL_EUROPE )**2 )
+    dist = dsqrt( (lon-LON_CRITICAL_EUROPE)**2 + (lat-LAT_CRITICAL_EUROPE )**2 )
     if( dist < CRITICAL_RANGE_EUROPE ) then
       ! sets up smoothing points
       ! by default uses CAP smoothing with crustmap resolution, e.g. 1/4 degree