[cig-commits] r16880 - seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN

[carltape at geodynamics.org]

Author: carltape
Date: 2010-06-03 07:54:27 -0700 (Thu, 03 Jun 2010)
New Revision: 16880

Modified:
   seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/compute_rho_estimate.f90
   seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_regions_mesh.f90
Log:
Added a call in create_regions_mesh.f90 to overwrite ALL rho values in the model using the Brocher (2005) Vp-rho scaling specified by compute_rho_estimate.f90, which has an added comment to indicate this. The default synthetics are now those of the m00 model of Tape et. al (2009,2010).


Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/compute_rho_estimate.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/compute_rho_estimate.f90	2010-06-03 14:41:26 UTC (rev 16879)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/compute_rho_estimate.f90	2010-06-03 14:54:27 UTC (rev 16880)
@@ -26,6 +26,8 @@
   subroutine compute_rho_estimate(rho,vp)
 
 ! compute rho estimate in Gocad block and in Hauksson model from Vp
+! -- this may be called in create_regions_mesh.f90 to apply the
+!    Vp-rho scaling throughout the entire volume
 
   implicit none
 

Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_regions_mesh.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_regions_mesh.f90	2010-06-03 14:41:26 UTC (rev 16879)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_regions_mesh.f90	2010-06-03 14:54:27 UTC (rev 16880)
@@ -745,6 +745,9 @@
   flag_sediments(i,j,k,ispec) = point_is_in_sediments
   if(point_is_in_sediments) not_fully_in_bedrock(ispec) = .true.
 
+! CHT: overwrite ALL density values using the Brocher2005 vp-rho scaling
+  call compute_rho_estimate(rho,vp)
+
 ! define elastic parameters in the model
 ! distinguish between single and double precision for reals
   if(ANISOTROPY) then