[cig-commits] r14133 - seismo/3D/SPECFEM3D_SESAME/trunk

dkomati1 at geodynamics.org dkomati1 at geodynamics.org
Mon Feb 23 18:04:09 PST 2009


Author: dkomati1
Date: 2009-02-23 18:04:08 -0800 (Mon, 23 Feb 2009)
New Revision: 14133

Modified:
   seismo/3D/SPECFEM3D_SESAME/trunk/write_AVS_DX_mesh_quality_data.f90
Log:
fixed a bug in the calculation of the diagonal aspect ratio of a mesh element


Modified: seismo/3D/SPECFEM3D_SESAME/trunk/write_AVS_DX_mesh_quality_data.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/write_AVS_DX_mesh_quality_data.f90	2009-02-24 01:22:42 UTC (rev 14132)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/write_AVS_DX_mesh_quality_data.f90	2009-02-24 02:04:08 UTC (rev 14133)
@@ -246,8 +246,8 @@
      dist2 = dsqrt((xelm(2) - xelm(8))**2 + (yelm(2) - yelm(8))**2 + (zelm(2) - zelm(8))**2)
      dist3 = dsqrt((xelm(3) - xelm(5))**2 + (yelm(3) - yelm(5))**2 + (zelm(3) - zelm(5))**2)
      dist4 = dsqrt((xelm(4) - xelm(6))**2 + (yelm(4) - yelm(6))**2 + (zelm(4) - zelm(6))**2)
-     distmin = dmin1(distmin,dist1,dist2,dist3,dist4)
-     distmax = dmax1(distmax,dist1,dist2,dist3,dist4)
+     distmin = dmin1(dist1,dist2,dist3,dist4)
+     distmax = dmax1(dist1,dist2,dist3,dist4)
      diagonal_aspect_ratio = distmax / distmin
 
 ! write mesh quality information for each element



More information about the CIG-COMMITS mailing list