[cig-commits] [commit] devel: 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 branch  : devel
Link       : https://github.com/geodynamics/specfem3d_globe/compare/754300e5a26e2a5d17069e661645fe819a7857a1...cd83b686a7c6014e132933f685157d43151a665e

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

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