[cig-commits] r13947 - seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2009-01-24 15:32:42 -0800 (Sat, 24 Jan 2009)
New Revision: 13947

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90
Log:
fixed a bug in the calculation of the CRC for restart files in the outer core


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90	2009-01-24 23:08:02 UTC (rev 13946)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90	2009-01-24 23:32:42 UTC (rev 13947)
@@ -1753,19 +1753,19 @@
 
       read(55) displ_outer_core
       read(55) icrc_read
-      total_size_array = CUSTOM_REAL*NDIM*NGLOB_OUTER_CORE
+      total_size_array = CUSTOM_REAL*NGLOB_OUTER_CORE
       call compute_icrc(icrc,0_I2B,displ_outer_core,total_size_array,0_I2B,1)
       if(icrc /= icrc_read) stop 'CRC error in restart file'
 
       read(55) veloc_outer_core
       read(55) icrc_read
-      total_size_array = CUSTOM_REAL*NDIM*NGLOB_OUTER_CORE
+      total_size_array = CUSTOM_REAL*NGLOB_OUTER_CORE
       call compute_icrc(icrc,0_I2B,veloc_outer_core,total_size_array,0_I2B,1)
       if(icrc /= icrc_read) stop 'CRC error in restart file'
 
       read(55) accel_outer_core
       read(55) icrc_read
-      total_size_array = CUSTOM_REAL*NDIM*NGLOB_OUTER_CORE
+      total_size_array = CUSTOM_REAL*NGLOB_OUTER_CORE
       call compute_icrc(icrc,0_I2B,accel_outer_core,total_size_array,0_I2B,1)
       if(icrc /= icrc_read) stop 'CRC error in restart file'
 
@@ -2918,17 +2918,17 @@
     write(55) icrc
 
     write(55) displ_outer_core
-    total_size_array = CUSTOM_REAL*NDIM*NGLOB_OUTER_CORE
+    total_size_array = CUSTOM_REAL*NGLOB_OUTER_CORE
     call compute_icrc(icrc,0_I2B,displ_outer_core,total_size_array,0_I2B,1)
     write(55) icrc
 
     write(55) veloc_outer_core
-    total_size_array = CUSTOM_REAL*NDIM*NGLOB_OUTER_CORE
+    total_size_array = CUSTOM_REAL*NGLOB_OUTER_CORE
     call compute_icrc(icrc,0_I2B,veloc_outer_core,total_size_array,0_I2B,1)
     write(55) icrc
 
     write(55) accel_outer_core
-    total_size_array = CUSTOM_REAL*NDIM*NGLOB_OUTER_CORE
+    total_size_array = CUSTOM_REAL*NGLOB_OUTER_CORE
     call compute_icrc(icrc,0_I2B,accel_outer_core,total_size_array,0_I2B,1)
     write(55) icrc