[cig-commits] r21317 - seismo/2D/SPECFEM2D/trunk/src/specfem2D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-02-01 05:47:37 -0800 (Fri, 01 Feb 2013)
New Revision: 21317

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_color_image.F90
Log:
fixed the blue color detection for water bug detected by Dara Farrell


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_color_image.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_color_image.F90	2013-02-01 13:33:04 UTC (rev 21316)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_color_image.F90	2013-02-01 13:47:37 UTC (rev 21317)
@@ -323,6 +323,7 @@
   double precision, dimension(4,3,numat) :: poroelastcoef
   double precision, dimension(numat) :: porosity,tortuosity
   double precision, dimension(NGLLX,NGLLX,nspec) :: vpext
+  double precision :: vp_of_the_model
 
 ! display acoustic layers as constant blue, because they likely correspond to water in the case of ocean acoustics
 ! or in the case of offshore oil industry experiments.
@@ -427,7 +428,18 @@
     if(DRAW_WATER_IN_BLUE .and. .not. elastic(ispec) .and. .not. poroelastic(ispec)) then
       do j = 1,NGLLZ
         do i = 1,NGLLX
-          vp_display(ibool(i,j,ispec)) = -1
+
+          !--- if external medium, get elastic parameters of current grid point
+          if(assign_external_model) then
+            vp_of_the_model = vpext(i,j,ispec)
+          else
+            vp_of_the_model = sqrt((lambdal_relaxed + 2.d0*mul_relaxed) / rhol)
+          endif
+
+! test that water is indeed water and not an acoustic version of a sediment for instance
+! thus check that Vp is the typical Vp of water
+          if(abs(vp_of_the_model - 1500.d0) <= 45.d0) vp_display(ibool(i,j,ispec)) = -1
+
         enddo
       enddo
     endif