[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