[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