[cig-commits] r12744 - in seismo/3D/SPECFEM3D_GLOBE/trunk: . version41_beta/src

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2008-08-28 07:58:33 -0700 (Thu, 28 Aug 2008)
New Revision: 12744

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/model_iasp91.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/model_iasp91.f90
Log:
added reference to original IASP91 paper; removed comment about model from Sebastien Chevrot


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/model_iasp91.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/model_iasp91.f90	2008-08-28 00:28:28 UTC (rev 12743)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/model_iasp91.f90	2008-08-28 14:58:33 UTC (rev 12744)
@@ -28,6 +28,9 @@
   subroutine model_iasp91(myrank,x,rho,vp,vs,Qkappa,Qmu,idoubling,ONE_CRUST,check_doubling_flag, &
                      RICB,RCMB,RTOPDDOUBLEPRIME,R771,R670,R400,R220,R120,RMOHO,RMIDDLE_CRUST)
 
+! from B. L. N. Kennett and E. R. Engdahl, Traveltimes for global earthquake location
+! and phase identification, Geophysical Journal International, vol. 105, p. 429-465 (1991)
+
   implicit none
 
   include "constants.h"
@@ -164,13 +167,6 @@
     Qmu=143.0d0
     Qkappa=57827.0d0
 
-! from Sebastien Chevrot: for the IASP91 model
-! Depth        R                Vp                    Vs
-! 0-20       6351-6371         5.80                  3.36
-! 20-35      6336-6351         6.50                  3.75
-! 35-120     6251-6336   8.78541-0.74953 x       6.706231-2.248585 x
-! with x = r / 6371
-
   else if(r > R220 .and. r <= R120) then
     rho=2.6910d0+0.6924d0*x
     vp=25.41389-17.69722*x

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/model_iasp91.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/model_iasp91.f90	2008-08-28 00:28:28 UTC (rev 12743)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/model_iasp91.f90	2008-08-28 14:58:33 UTC (rev 12744)
@@ -28,11 +28,14 @@
   subroutine model_iasp91(myrank,x,rho,vp,vs,Qkappa,Qmu,idoubling,ONE_CRUST,check_doubling_flag, &
                      RICB,RCMB,RTOPDDOUBLEPRIME,R771,R670,R400,R220,R120,RMOHO,RMIDDLE_CRUST)
 
+! from B. L. N. Kennett and E. R. Engdahl, Traveltimes for global earthquake location
+! and phase identification, Geophysical Journal International, vol. 105, p. 429-465 (1991)
+
   implicit none
 
   include "constants.h"
 
-! given a normalized radius x, gives the non-dimesionalized density rho,
+! given a normalized radius x, gives the non-dimensionalized density rho,
 ! speeds vp and vs, and the quality factors Qkappa and Qmu
 
   logical check_doubling_flag
@@ -164,13 +167,6 @@
     Qmu=143.0d0
     Qkappa=57827.0d0
 
-! from Sebastien Chevrot: for the IASP91 model
-! Depth        R                Vp                    Vs
-! 0-20       6351-6371         5.80                  3.36
-! 20-35      6336-6351         6.50                  3.75
-! 35-120     6251-6336   8.78541-0.74953 x       6.706231-2.248585 x
-! with x = r / 6371
-
   else if(r > R220 .and. r <= R120) then
     rho=2.6910d0+0.6924d0*x
     vp=25.41389-17.69722*x