[cig-commits] r21815 - seismo/3D/SPECFEM3D/trunk/src/generate_databases

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-04-11 09:01:18 -0700 (Thu, 11 Apr 2013)
New Revision: 21815

Modified:
   seismo/3D/SPECFEM3D/trunk/src/generate_databases/get_absorbing_boundary.f90
Log:
add one comment on case when set PML_CONDITIONS and PML_INSTEAD_OF_FREE_SURFACE to be .true.


Modified: seismo/3D/SPECFEM3D/trunk/src/generate_databases/get_absorbing_boundary.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/generate_databases/get_absorbing_boundary.f90	2013-04-11 15:54:04 UTC (rev 21814)
+++ seismo/3D/SPECFEM3D/trunk/src/generate_databases/get_absorbing_boundary.f90	2013-04-11 16:01:18 UTC (rev 21815)
@@ -513,6 +513,10 @@
     if( STACEY_INSTEAD_OF_FREE_SURFACE ) then
        write(IMAIN,*) '     absorbing boundary includes free surface'
     endif
+!in curent code when set PML_CONDITIONS and PML_INSTEAD_OF_FREE_SURFACE to be .true. we should also
+!provide the free_surface_file, since we need the free_surface_file to determin ibelm_top which is the
+!outer boundary of top CPML layer. Thus when set PML_CONDITIONS and PML_INSTEAD_OF_FREE_SURFACE to be .true.
+!we should suppress the following four lines.
     if( PML_INSTEAD_OF_FREE_SURFACE .and. itop /= 0 ) then
        print*,'please check Par_file/free_surface_file and recompile solver'
        stop 'error: number of free surface faces should be zero when PML_INSTEAD_OF_FREE_SURFACE is set to .true.'