[cig-commits] r22692 - seismo/3D/SPECFEM3D/trunk/src/specfem3D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-07-31 09:49:59 -0700 (Wed, 31 Jul 2013)
New Revision: 22692

Modified:
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90
Log:
fixed the CPML array allocation bug mentioned by Timo L?\195?\164hivaara


Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90	2013-07-31 16:48:31 UTC (rev 22691)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90	2013-07-31 16:49:59 UTC (rev 22692)
@@ -370,8 +370,14 @@
     read(27) CPML_width_x
     read(27) CPML_width_y
     read(27) CPML_width_z
+
     allocate(is_CPML(NSPEC_AB),stat=ier) !need to be optimized
     if(ier /= 0) stop 'error allocating array is_CPML'
+
+! make sure there are no PMLs by default,
+! and then below if NSPEC_CPML > 0 we will need the real flags for this mesh from the disk
+    is_CPML(:) = .false.
+
     if( NSPEC_CPML > 0 ) then
       allocate(CPML_regions(NSPEC_CPML),stat=ier)
       if(ier /= 0) stop 'error allocating array CPML_regions'