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

[danielpeter at geodynamics.org]

Author: danielpeter
Date: 2011-04-04 08:57:44 -0700 (Mon, 04 Apr 2011)
New Revision: 18168

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D.f90
Log:
fixes a check for anisotropic elements in meshfem3D.f90

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D.f90	2011-04-04 04:19:03 UTC (rev 18167)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D.f90	2011-04-04 15:57:44 UTC (rev 18168)
@@ -678,13 +678,15 @@
 
     enddo
 
-    ! store number of anisotropic elements found in the mantle
-    if(nspec_tiso /= 0 .and. iregion_code /= IREGION_CRUST_MANTLE) &
+    ! checks number of anisotropic elements found in the mantle
+    if(iregion_code /= IREGION_CRUST_MANTLE .and. nspec_tiso /= 0 ) &
       call exit_MPI(myrank,'found anisotropic elements outside of the mantle')
 
-    if(iregion_code == IREGION_CRUST_MANTLE .and. nspec_tiso == 0) &
-      call exit_MPI(myrank,'found no anisotropic elements in the mantle')
-
+    if( TRANSVERSE_ISOTROPY ) then
+      if(iregion_code == IREGION_CRUST_MANTLE .and. nspec_tiso == 0) &
+        call exit_MPI(myrank,'found no anisotropic elements in the mantle')
+    endif
+    
     ! computes total area and volume
     call meshfem3D_compute_area(myrank,NCHUNKS,iregion_code, &
                               area_local_bottom,area_local_top,&