[cig-commits] r20388 - in seismo/3D/SPECFEM3D_GLOBE/trunk/src: meshfem3D specfem3D

[danielpeter at geodynamics.org]

Author: danielpeter
Date: 2012-06-19 15:12:59 -0700 (Tue, 19 Jun 2012)
New Revision: 20388

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/create_mass_matrices.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/prepare_timerun.f90
Log:
fixes test for negative mass matrix

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/create_mass_matrices.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/create_mass_matrices.f90	2012-06-19 22:12:35 UTC (rev 20387)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/create_mass_matrices.f90	2012-06-19 22:12:59 UTC (rev 20388)
@@ -675,6 +675,7 @@
   endif
 
   ! check that mass matrix is positive
-  if(minval(rmassz(:)) <= 0.) call exit_MPI(myrank,'negative rmassz matrix term')
+  ! note: in fictitious elements it is still zero
+  if(minval(rmassz(:)) < 0.) call exit_MPI(myrank,'negative rmassz matrix term')
 
   end subroutine create_mass_matrices

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/prepare_timerun.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/prepare_timerun.f90	2012-06-19 22:12:35 UTC (rev 20387)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/prepare_timerun.f90	2012-06-19 22:12:59 UTC (rev 20388)
@@ -300,7 +300,7 @@
 
   ! suppress fictitious mass matrix elements in central cube
   ! because the slices do not compute all their spectral elements in the cube
-  where(rmass_inner_core(:) <= 0.) rmass_inner_core = 1.
+  where(rmass_inner_core(:) <= 0.) rmass_inner_core = 1._CUSTOM_REAL
 
   end subroutine prepare_timerun_centralcube