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

dkomati1 at geodynamics.org dkomati1 at geodynamics.org
Wed Mar 20 19:46:45 PDT 2013


Author: dkomati1
Date: 2013-03-20 19:46:44 -0700 (Wed, 20 Mar 2013)
New Revision: 21595

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
Log:
better print statements for PML


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90	2013-03-21 02:22:59 UTC (rev 21594)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90	2013-03-21 02:46:44 UTC (rev 21595)
@@ -187,7 +187,7 @@
 #else
      nspec_PML_tot = nspec_PML
 #endif
-     if(myrank == 0) write(IOUT,*) "number of PML spectral elements on side ",ibound,":",nspec_PML_tot
+     if(myrank == 0) write(IOUT,*) "Number of PML spectral elements on side ",ibound,":",nspec_PML_tot
 
      enddo ! end loop on the 4 boundaries
 
@@ -318,7 +318,7 @@
 #else
      nspec_PML_tot = nspec_PML
 #endif
-     if(myrank == 0) write(IOUT,*) "number of PML spectral elements :", nspec_PML_tot
+     if(myrank == 0) write(IOUT,*) "Total number of PML spectral elements :", nspec_PML_tot
 
      endif
 
@@ -340,7 +340,7 @@
 #else
      nspec_PML_tot = nspec_PML
 #endif
-     if(myrank == 0) write(IOUT,*) "number of PML spectral elements :", nspec_PML_tot
+     if(myrank == 0) write(IOUT,*) "Total number of PML spectral elements :", nspec_PML_tot
 
   endif
 
@@ -679,18 +679,18 @@
   d0_z_top = - (NPOWER + 1) * vpmax * log(Rcoef) / (2.d0 * thickness_PML_z_top)
 
 !   if (myrank == 0) then
-      write(IOUT,*)
-      write(IOUT,*) 'PML properties -------',myrank,'myrank'
-      write(IOUT,*) '     Vpmax=', vpmax
-      write(IOUT,*) '     log(Rcoef)=',log(Rcoef)
-      write(IOUT,*) '     thickness_PML_z_bottom =',thickness_PML_z_bottom
-      write(IOUT,*) '     thickness_PML_x_right  =',thickness_PML_x_right
-      write(IOUT,*) '     thickness_PML_z_top    =',thickness_PML_z_top
-      write(IOUT,*) '     thickness_PML_x_left   =',thickness_PML_x_left
-      write(IOUT,*) '     d0_bottom       =', d0_z_bottom
-      write(IOUT,*) '     d0_right        =', d0_x_right
-      write(IOUT,*) '     d0_top          =', d0_z_top
-      write(IOUT,*) '     d0_left         =', d0_x_left
+!     write(IOUT,*)
+!     write(IOUT,*) 'PML properties -------',myrank,'myrank'
+!     write(IOUT,*) '     Vpmax=', vpmax
+!     write(IOUT,*) '     log(Rcoef)=',log(Rcoef)
+!     write(IOUT,*) '     thickness_PML_z_bottom =',thickness_PML_z_bottom
+!     write(IOUT,*) '     thickness_PML_x_right  =',thickness_PML_x_right
+!     write(IOUT,*) '     thickness_PML_z_top    =',thickness_PML_z_top
+!     write(IOUT,*) '     thickness_PML_x_left   =',thickness_PML_x_left
+!     write(IOUT,*) '     d0_bottom       =', d0_z_bottom
+!     write(IOUT,*) '     d0_right        =', d0_x_right
+!     write(IOUT,*) '     d0_top          =', d0_z_top
+!     write(IOUT,*) '     d0_left         =', d0_x_left
 !   endif
 
    d_x = 0.d0



More information about the CIG-COMMITS mailing list