[cig-commits] r21751 - seismo/3D/SPECFEM3D/trunk/src/decompose_mesh

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-04-07 16:35:37 -0700 (Sun, 07 Apr 2013)
New Revision: 21751

Modified:
   seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/decompose_mesh.F90
Log:
now check that the number of CPML elements is not zero when PML_CONDITIONS is set to true


Modified: seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/decompose_mesh.F90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/decompose_mesh.F90	2013-04-07 23:31:34 UTC (rev 21750)
+++ seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/decompose_mesh.F90	2013-04-07 23:35:37 UTC (rev 21751)
@@ -662,6 +662,10 @@
        read(98,*) nspec_cpml, CPML_width
     endif
 
+! sanity check
+    if( PML_CONDITIONS .and. nspec_cpml <= 0) &
+        stop 'error: PML_CONDITIONS is set to true but nspec_cpml <= 0 in file absorbing_cpml_file'
+
     ! C-PML spectral elements global indexing
     allocate(CPML_to_spec(nspec_cpml),stat=ier)
     if(ier /= 0) stop 'error allocating array CPML_to_spec'