[cig-commits] [commit] : 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). (d1fd9ba)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/specfem3d

On branch  : 
Link       : https://github.com/geodynamics/specfem2d/compare/1e201257d91c794056b990a43329e05d04f77454...0000000000000000000000000000000000000000

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

commit d1fd9ba38147dcdc4d8da4d7ceea403d00ab6087
Author: Carl Tape <carltape at gi.alaska.edu>
Date:   Thu Jun 3 14:54:27 2010 +0000

    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).


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

d1fd9ba38147dcdc4d8da4d7ceea403d00ab6087
 compute_rho_estimate.f90 | 2 ++
 create_regions_mesh.f90  | 3 +++
 2 files changed, 5 insertions(+)

diff --git a/compute_rho_estimate.f90 b/compute_rho_estimate.f90
index 50ce76a..19dcb84 100644
--- a/compute_rho_estimate.f90
+++ b/compute_rho_estimate.f90
@@ -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
 
diff --git a/create_regions_mesh.f90 b/create_regions_mesh.f90
index 99dd020..71c6b83 100644
--- a/create_regions_mesh.f90
+++ b/create_regions_mesh.f90
@@ -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