[cig-commits] r19887 - seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-03-27 15:31:55 -0700 (Tue, 27 Mar 2012)
New Revision: 19887

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D.f90
Log:
put calculation of the exact area back for the top surface


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D.f90	2012-03-27 22:29:16 UTC (rev 19886)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D.f90	2012-03-27 22:31:55 UTC (rev 19887)
@@ -1213,7 +1213,8 @@
     write(IMAIN,*) '   calculated top area: ',area_total_top
 
     ! compare to exact theoretical value
-    if(NCHUNKS == 6 .and. .not. TOPOGRAPHY) then
+    if((NCHUNKS == 6 .or. (abs(ANGULAR_WIDTH_XI_IN_DEGREES - 90.d0) < TINYVAL .and. &
+                           abs(ANGULAR_WIDTH_ETA_IN_DEGREES - 90.d0) < TINYVAL)) .and. .not. TOPOGRAPHY) then
       select case(iregion_code)
         case(IREGION_CRUST_MANTLE)
           write(IMAIN,*) '            exact area: ',dble(NCHUNKS)*(4.0d0/6.0d0)*PI*R_UNIT_SPHERE**2