[cig-commits] r11918 - seismo/3D/SPECFEM3D_GLOBE/trunk

[dmichea at geodynamics.org]

Author: dmichea
Date: 2008-05-06 07:40:30 -0700 (Tue, 06 May 2008)
New Revision: 11918

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/read_compute_parameters.f90
Log:
fixed a bug : integer limit overflow for huge meshes.



Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/read_compute_parameters.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/read_compute_parameters.f90	2008-05-05 19:17:31 UTC (rev 11917)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/read_compute_parameters.f90	2008-05-06 14:40:30 UTC (rev 11918)
@@ -2248,12 +2248,12 @@
             nb_lay_sb = 0
             nspec_sb = 0
         endif
-        tmp_sum = tmp_sum + ((NEX_XI / ratio_sampling_array(iter_layer)) * (NEX_ETA / ratio_sampling_array(iter_layer)) * &
+        tmp_sum = tmp_sum + (((NEX_XI / ratio_sampling_array(iter_layer)) * (NEX_ETA / ratio_sampling_array(iter_layer)) * &
                 (ner(iter_layer) - doubling*nb_lay_sb)) + &
                 doubling * ((NEX_XI / ratio_sampling_array(iter_layer)) * (NEX_ETA / ratio_sampling_array(iter_layer)) * &
-                (nspec_sb/4))
+                (nspec_sb/4))) / NPROC
     enddo
-    NSPEC(iter_region) = tmp_sum / NPROC
+    NSPEC(iter_region) = tmp_sum
 enddo
 
   if(INCLUDE_CENTRAL_CUBE) NSPEC(IREGION_INNER_CORE) = NSPEC(IREGION_INNER_CORE) + &