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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2009-02-24 16:23:19 -0800 (Tue, 24 Feb 2009)
New Revision: 14140

Added:
   seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/convert_skewness_to_angle.f90
Modified:
   seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90
Log:
slightly improved the output when checking a mesh;
also added a small program to convert mesh skewness to mesh angles


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 23:51:23 UTC (rev 14139)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90	2009-02-25 00:23:19 UTC (rev 14140)
@@ -201,7 +201,7 @@
 
   print *
   print *,'minimum length of an edge in the whole mesh (m) = ',distmin
-  print *,'minimum length divided in ',NGLLX-1,' average segments (m) = ',distmin/real(NGLLX-1)
+! print *,'minimum length divided in ',NGLLX-1,' average segments (m) = ',distmin/real(NGLLX-1)
   print *,'maximum length of an edge in the whole mesh (m) = ',distmax
   print *
 
@@ -255,7 +255,7 @@
   print *,'max equiangle skewness = ',equiangle_skewness_max
 ! print *,'min equiangle skewness = ',equiangle_skewness_min
   print *
-  print *,'skewness max deviation angle from 90 degrees = ',90.*equiangle_skewness_max
+  print *,'max deviation angle from a right angle (90 degrees) is therefore = ',90.*equiangle_skewness_max
   print *
   print *,'worst angle in the mesh is therefore ',90.*(1. - equiangle_skewness_max)
   print *,'or ',180. - 90.*(1. - equiangle_skewness_max),' degrees'

Added: seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/convert_skewness_to_angle.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/convert_skewness_to_angle.f90	                        (rev 0)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/convert_skewness_to_angle.f90	2009-02-25 00:23:19 UTC (rev 14140)
@@ -0,0 +1,53 @@
+!=====================================================================
+!
+!               S p e c f e m 3 D  V e r s i o n  1 . 4
+!               ---------------------------------------
+!
+!                 Dimitri Komatitsch and Jeroen Tromp
+!    Seismological Laboratory - California Institute of Technology
+!         (c) California Institute of Technology September 2006
+!
+! This program is free software; you can redistribute it and/or modify
+! it under the terms of the GNU General Public License as published by
+! the Free Software Foundation; either version 2 of the License, or
+! (at your option) any later version.
+!
+! This program is distributed in the hope that it will be useful,
+! but WITHOUT ANY WARRANTY; without even the implied warranty of
+! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+! GNU General Public License for more details.
+!
+! You should have received a copy of the GNU General Public License along
+! with this program; if not, write to the Free Software Foundation, Inc.,
+! 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
+!
+!=====================================================================
+
+! computes and displays the min or max angle in a mesh that give rise
+! to a certain value of skewness
+
+! Dimitri Komatitsch, University of Pau, France, February 2009.
+
+  program convert_skewness_to_angle
+
+  implicit none
+
+  include "constants.h"
+
+  double precision :: equiangle_skewness
+
+  print *
+  print *,'enter equiangle skewness to convert (0.=perfect, 1.=bad):'
+  read(5,*) equiangle_skewness
+
+  print *
+  print *,'equiangle skewness = ',equiangle_skewness
+  print *
+  print *,'deviation angle from a right angle (90 degrees) is therefore = ',90.*equiangle_skewness
+  print *
+  print *,'worst angle in the mesh is therefore ',90.*(1. - equiangle_skewness)
+  print *,'or ',180. - 90.*(1. - equiangle_skewness),' degrees'
+  print *
+
+  end program convert_skewness_to_angle
+