[cig-commits] r15572 - in seismo/2D/SPECFEM2D/trunk: . DATA

[cmorency at geodynamics.org]

Author: cmorency
Date: 2009-08-21 14:11:59 -0700 (Fri, 21 Aug 2009)
New Revision: 15572

Modified:
   seismo/2D/SPECFEM2D/trunk/DATA/Par_file
   seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
Another fix


Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file	2009-08-21 20:32:56 UTC (rev 15571)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file	2009-08-21 21:11:59 UTC (rev 15572)
@@ -113,7 +113,7 @@
 19 1 1946.8d0 2100.d0   1212.27d0 0 0 10.d0 10.d0 0 0 0 0 0 0
 20 1 1857.9d0 1900.d0   1096.80d0 0 0 10.d0 10.d0 0 0 0 0 0 0
 21 1 2087.6d0 2500.d0   1443.35d0 0 0 10.d0 10.d0 0 0 0 0 0 0
-22 3 1986.3 1040.3 0.2 2.0 1d-11 0.0 1d-11 5.341d9 5.341d8 5.341d8 0.0d-4 3.204d9 10.d0
+22 3 1986.3 1040.3 0.4 2.0 1d-11 0.0 1d-11 5.341d9 5.341d8 5.341d8 0.0d-4 3.204d9 10.d0
 #22 3 2200.d0 786.3d0  0.4 2.0 1d-11 0.0 1d-11 5.341d9 2d9 3d9 0.0d-4 3.204d9 10.d0
 #22 1 1634.52d0 2272.89d0   1472.71d0 0 0 10.d0 10.d0 0 0 0 0 0 0
 # define the different regions of the model in the (nx,nz) spectral element mesh

Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-08-21 20:32:56 UTC (rev 15571)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-08-21 21:11:59 UTC (rev 15572)
@@ -5078,19 +5078,15 @@
      if(iglob /= iglob2) &
             call exit_MPI( 'error in solid/porous iglob detection')
 
-           displ(1,iglob)=(valence_poroelastic(iglob)*displs_poroelastic(1,iglob)&
-                              +valence_elastic(iglob)*displ_elastic(1,iglob))/ &
-                          (valence_poroelastic(iglob)+valence_elastic(iglob))
-           displ(2,iglob)=(valence_poroelastic(iglob)*displs_poroelastic(2,iglob) &
-                              +valence_elastic(iglob)*displ_elastic(2,iglob))/ &
-                          (valence_poroelastic(iglob)+valence_elastic(iglob))
+           displ(1,iglob)=(2.d0*displs_poroelastic(1,iglob) + &
+                          displ_elastic(1,iglob2))/3.d0
+           displ(2,iglob)=(displs_poroelastic(2,iglob)  + &
+                          displ_elastic(2,iglob2))/3.d0
 
-           veloc(1,iglob)=(valence_poroelastic(iglob)*velocs_poroelastic(1,iglob) &
-                              +valence_elastic(iglob)*veloc_elastic(1,iglob))/ &
-                          (valence_poroelastic(iglob)+valence_elastic(iglob))
-           veloc(2,iglob)=(valence_poroelastic(iglob)*velocs_poroelastic(2,iglob) &
-                              +valence_elastic(iglob)*veloc_elastic(2,iglob))/ &
-                          (valence_poroelastic(iglob)+valence_elastic(iglob))
+           veloc(1,iglob)=(2.d0*velocs_poroelastic(1,iglob)  +&
+                          veloc_elastic(1,iglob2))/3.d0
+           veloc(2,iglob)=(2.d0*velocs_poroelastic(2,iglob)  +&
+                          veloc_elastic(2,iglob2))/3.d0
 
         enddo
       enddo