[cig-commits] r8577 - seismo/2D/SPECFEM2D/trunk

[walter at geodynamics.org]

Author: walter
Date: 2007-12-07 15:57:56 -0800 (Fri, 07 Dec 2007)
New Revision: 8577

Modified:
   seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
fixed issue due to line length in specfem2D.F90. fixed bug encountered with attenuation (arrays e1, e11 and e13 were allocated but not initialized); depending on the platform/compiler, the simulation could blow up.

Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2007-09-11 14:42:31 UTC (rev 8576)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2007-12-07 23:57:56 UTC (rev 8577)
@@ -549,6 +549,9 @@
   allocate(e1(NGLLX,NGLLZ,nspec_allocate,N_SLS))
   allocate(e11(NGLLX,NGLLZ,nspec_allocate,N_SLS))
   allocate(e13(NGLLX,NGLLZ,nspec_allocate,N_SLS))
+  e1(:,:,:,:) = 0._CUSTOM_REAL
+  e11(:,:,:,:) = 0._CUSTOM_REAL
+  e13(:,:,:,:) = 0._CUSTOM_REAL
 
   allocate(dux_dxl_n(NGLLX,NGLLZ,nspec_allocate))
   allocate(duz_dzl_n(NGLLX,NGLLZ,nspec_allocate))
@@ -864,9 +867,11 @@
 #endif
 
 ! for acoustic
-  if(TURN_ANISOTROPY_ON .and. .not. any_elastic_glob) call exit_MPI('cannot have anisotropy if acoustic simulation only')
+  if(TURN_ANISOTROPY_ON .and. .not. any_elastic_glob) &
+    call exit_MPI('cannot have anisotropy if acoustic simulation only')
 
-  if(TURN_ATTENUATION_ON .and. .not. any_elastic_glob) call exit_MPI('currently cannot have attenuation if acoustic simulation only')
+  if(TURN_ATTENUATION_ON .and. .not. any_elastic_glob) &
+    call exit_MPI('currently cannot have attenuation if acoustic simulation only')
 
 ! for attenuation
   if(TURN_ANISOTROPY_ON .and. TURN_ATTENUATION_ON) then