[cig-commits] [commit] devel: 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
cig_noreply at geodynamics.org
Tue Apr 22 16:09:34 PDT 2014
Repository : ssh://geoshell/specfem3d_globe
On branch : devel
Link : https://github.com/geodynamics/specfem3d_globe/compare/1a5b4a182b85c2dbf7da31d567c33c1985eb474c...9f9faf71ddb9dd81afedf5968b94ef1b47757f75
>---------------------------------------------------------------
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)
More information about the CIG-COMMITS
mailing list