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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-03-01 10:02:54 -0800 (Fri, 01 Mar 2013)
New Revision: 21427

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
Log:
fix one bug in pml when turn on READ_EXTERNAL_SEP_FILE in par_file


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90	2013-03-01 16:42:15 UTC (rev 21426)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90	2013-03-01 18:02:54 UTC (rev 21427)
@@ -462,6 +462,10 @@
               else
                 rhol = rhoext(i,j,ispec)
               endif
+                cpl = vpext(i,j,ispec)
+                !assuming that in fluid(acoustic) part input cpl is defined by sqrt(kappal/rhol), &
+                !which is not the same as in cpl input in elastic part
+                kappal = rhol*cpl*cpl 
             else
               if(CUSTOM_REAL == SIZE_REAL) then
                 lambdal_relaxed = sngl(poroelastcoef(1,1,kmato(ispec)))