[cig-commits] [commit] devel, master: cosmetic changes to be able to later apply an old patch from Thomas Gillet to significantly speed up the mesher for model s362ani (40002cb)

[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 40002cb370e888493a8a1b814191555766bf15ed
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Wed Apr 23 01:06:39 2014 +0200

    cosmetic changes to be able to later apply an old patch from Thomas Gillet to significantly speed up the mesher for model s362ani


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

40002cb370e888493a8a1b814191555766bf15ed
 src/meshfem3D/model_s362ani.f90 | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/src/meshfem3D/model_s362ani.f90 b/src/meshfem3D/model_s362ani.f90
index 06d000b..52fd052 100644
--- a/src/meshfem3D/model_s362ani.f90
+++ b/src/meshfem3D/model_s362ani.f90
@@ -32,7 +32,7 @@
 !
 ! In this model, radial anisotropy is confined to the uppermost mantle.
 ! The model (and the corresponding mesh) incorporate
-! tomography on the 650~km and 410~km discontinuities in the 1D reference model REF.
+! tomography on the 650-km and 410-km discontinuities in the 1D reference model REF.
 !
 ! s362wmani: A version of S362ANI with anisotropy allowed throughout the mantle.
 !
@@ -1091,8 +1091,8 @@
   ncon=0
 
   do iver=1,nver
-    if(xlat > verlat(iver)-2.*verrad(iver)) then
-      if(xlat < verlat(iver)+2.*verrad(iver)) then
+    if(xlat > verlat(iver) - 2.*verrad(iver)) then
+      if(xlat < verlat(iver) + 2.*verrad(iver)) then
         ver8=DEGREES_TO_RADIANS*(verlat(iver))
         xla8=DEGREES_TO_RADIANS*(xlat)
         dd=sin(ver8)*sin(xla8) + cos(ver8)*cos(xla8)* cos(DEGREES_TO_RADIANS*(xlon-verlon(iver)))
@@ -1755,6 +1755,7 @@
   implicit none
 
   integer :: i,M,MP1,k,l,LP1
+
   real(kind=4) :: THETA,DSFL3,COSEC,SFL3
 
   real(kind=4) :: X(M+1),XP(M+1),XCOSEC(M+1)