[cig-commits] r20687 - in seismo/2D/SPECFEM2D/trunk: EXAMPLES/semi_infinite_homo src/specfem2D

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2012-09-05 02:27:25 -0700 (Wed, 05 Sep 2012)
New Revision: 20687

Modified:
   seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/assemble_MPI.F90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
fix the wrong output the file name of PML_ELEMENT_FILE in specfem2D.F90 and two small bugs in other two files


Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D	2012-09-05 04:00:23 UTC (rev 20686)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D	2012-09-05 09:27:25 UTC (rev 20687)
@@ -101,7 +101,7 @@
 # III: (model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu).
 # For istropic elastic/acoustic material use I and set Vs to zero to make a given model acoustic, for anisotropic elastic use II,
 # and for isotropic poroelastic material use III. The mesh can contain acoustic, elastic, and poroelastic models simultaneously.
-1 1 2700.d0 3000.d0 0.d0 0 0 20. 10. 0 0 0 0 0 0
+1 1 2700.d0 3000.d0 1732.05d0 0 0 20. 10. 0 0 0 0 0 0
 
 # external mesh or not
 read_external_mesh              = .false.

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/assemble_MPI.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/assemble_MPI.F90	2012-09-05 04:00:23 UTC (rev 20686)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/assemble_MPI.F90	2012-09-05 09:27:25 UTC (rev 20687)
@@ -445,7 +445,7 @@
 
 
 !-----------------------------------------------
-! Assembling accel_elastic for elastic elements :
+! Assembling accel_elastic for poroelastic elements :
 ! the buffers are filled, the ISEND and IRECV are started here, then
 ! contributions are added.
 ! The previous version included communication overlap using persistent

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2012-09-05 04:00:23 UTC (rev 20686)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2012-09-05 09:27:25 UTC (rev 20687)
@@ -2828,10 +2828,6 @@
       is_PML(:) = .false.
       which_PML_elem(:,:) = .false.
 
-
-      write(*,*)read_external_mesh,'read_external_mesh'
-      write(*,*)CPML_element_file,'CPML_element_file'
-
       call pml_init(nspec,nglob,anyabs,ibool,nelemabs,codeabs,numabs,&
                   nspec_PML,is_PML,which_PML_elem,spec_to_PML, &
                   icorner_iglob,NELEM_PML_THICKNESS,&