[cig-commits] r14136 - seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2009-02-24 08:35:47 -0800 (Tue, 24 Feb 2009)
New Revision: 14136

Modified:
   seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90
Log:
now display the CUBIT element number of elements with bad skewness


Modified: seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90	2009-02-24 14:51:55 UTC (rev 14135)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90	2009-02-24 16:35:47 UTC (rev 14136)
@@ -44,10 +44,10 @@
 ! number of points and of hex or quad elements
 ! number of points of a hex or quad element
 
-  character(len=100), parameter :: cubit_mesh_file = 'HOMOGENE_3D_lisse_300_in_meters.inp'
-  integer, parameter :: NPOIN = 98692, NSPEC = 90585, NGNOD = 8
-  logical, parameter :: IGNORE_OTHER_HEADERS = .false.
-  double precision, parameter :: delta_t = 1.d-3
+! character(len=100), parameter :: cubit_mesh_file = 'HOMOGENE_3D_lisse_300_in_meters.inp'
+! integer, parameter :: NPOIN = 98692, NSPEC = 90585, NGNOD = 8
+! logical, parameter :: IGNORE_OTHER_HEADERS = .false.
+! double precision, parameter :: delta_t = 1.d-3
 
 ! character(len=100), parameter :: cubit_mesh_file = 'regolite_3D_rego3d_70m_in_meters.inp'
 ! integer, parameter :: NPOIN = 4050696, NSPEC = 3410265, NGNOD = 8
@@ -59,10 +59,10 @@
 ! logical, parameter :: IGNORE_OTHER_HEADERS = .false.
 ! double precision, parameter :: delta_t = 5.d-3
 
-! character(len=100), parameter :: cubit_mesh_file = 'Eros_complet_regolite_fractures_2D_in_meters.inp'
-! integer, parameter :: NPOIN = 57630, NSPEC = 107904, NGNOD = 4
-! logical, parameter :: IGNORE_OTHER_HEADERS = .true.
-! double precision, parameter :: delta_t = 1.5d-4
+  character(len=100), parameter :: cubit_mesh_file = 'eros_complexe_2d_regolite_fractures_in_meters.inp'
+  integer, parameter :: NPOIN = 57807, NSPEC = 56983, NGNOD = 4
+  logical, parameter :: IGNORE_OTHER_HEADERS = .true.
+  double precision, parameter :: delta_t = 1.5d-4
 
   double precision, parameter :: VP_MAX = 3000.d0
 
@@ -97,7 +97,8 @@
 
   if(NGNOD /= 4 .and. NGNOD /= 8) stop 'NGNOD should be 4 or 8'
 
-  print *,'1 = output bad elements in OpenDX format'
+  print *
+  print *,'1 = output bad elements above a certain skewness threshold in OpenDX format'
   print *,'2 = do not output bad elements'
   print *
   print *,'enter value:'
@@ -154,7 +155,7 @@
     if(NGNOD == 4) then
 
 !! DK DK ignore other headers for 2D mesh of Eros with fractures, which has multiple material sets
-      if(IGNORE_OTHER_HEADERS .and. (i == 51099 .or. i == 105264)) read(10,*)
+      if(IGNORE_OTHER_HEADERS .and. (i == 51276 .or. i == 53917)) read(10,*)
 
       read(10,*) iread,ibool(1,i),ibool(2,i),ibool(3,i),ibool(4,i)
     else if(NGNOD == 8) then
@@ -236,7 +237,6 @@
 
   enddo
 
-  print *
   print *,'------------'
   print *,'mesh quality parameter definitions'
   print *
@@ -370,7 +370,13 @@
 
   enddo
 
-  if(ntotspecAVS_DX == 0) stop 'no elements in skewness range, no file created'
+  if(ntotspecAVS_DX == 0) then
+    stop 'no elements in skewness range, no file created'
+  else
+    print *
+    print *,'there are ',ntotspecAVS_DX,' elements in AVS or DX skewness range ',skewness_AVS_DX_min,skewness_AVS_DX_max
+    print *
+  endif
 
   open(unit=11,file='DX_mesh_quality.dx',status='unknown')
 
@@ -413,6 +419,7 @@
             ibool(4,ispec)-1, ibool(1,ispec)-1, ibool(8,ispec)-1, ibool(5,ispec)-1, &
             ibool(3,ispec)-1, ibool(2,ispec)-1, ibool(7,ispec)-1, ibool(6,ispec)-1
       endif
+      print *,'element ',ispec,' belongs to the range and has skewness = ',sngl(equiangle_skewness)
     endif
 
   enddo
@@ -457,10 +464,6 @@
 ! close OpenDX file
   close(11)
 
-  print *
-  print *,'there are ',ntotspecAVS_DX,' elements in AVS or DX skewness range ',skewness_AVS_DX_min,skewness_AVS_DX_max
-  print *
-
   endif
 
   end program check_mesh_quality_CUBIT_Abaqus