Merge pull request #69 from komatits/devel

fixed a bug in Matthias Meschede's crustal model smoothing routines: lat...
SPECFEM · May 2, 2014 · bed4deb · bed4deb
2 parents 754300e + b173684
commit bed4deb
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/meshfem3D/model_crustmaps.f90 b/src/meshfem3D/model_crustmaps.f90
@@ -370,7 +370,7 @@ subroutine read_crustmaps(lat,lon,velp,vels,rhos,thicks)
 
   ! 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