[cig-commits] r22526 - seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-07-06 16:04:30 -0700 (Sat, 06 Jul 2013)
New Revision: 22526

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/save_header_file.f90
Log:
first step of saving memory when 3D attenuation is off


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/save_header_file.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/save_header_file.f90	2013-07-06 22:56:26 UTC (rev 22525)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/save_header_file.f90	2013-07-06 23:04:30 UTC (rev 22526)
@@ -478,9 +478,17 @@
   write(IOUT,*) 'integer, parameter :: NCORNERSCHUNKS_VAL = ',NCORNERSCHUNKS
 
   if(ATTENUATION) then
-    att1     = NGLLX
-    att2     = NGLLY
-    att3     = NGLLZ
+!! DK DK July 2013: to save a huge amount of memory, when 3D attenuation is off it is sufficient to save a single point
+!! DK DK July 2013: per spectral element because the Q attenuation factor is then constant per layer of the geological model
+    if(ATTENUATION_3D) then
+      att1     = NGLLX
+      att2     = NGLLY
+      att3     = NGLLZ
+    else
+      att1     = 1
+      att2     = 1
+      att3     = 1
+    endif
     att4     = NSPEC(IREGION_CRUST_MANTLE)
     att5     = NSPEC(IREGION_INNER_CORE)
   else