January 29 2007   7:38:13.471 AM
 
TRIANGULATION_REFINE
  FORTRAN90 version
  Read a "linear" or "quadratic" triangulation
  and write out a refined triangulation.
 
  In particular:
 
  Read a dataset of NODE_NUM1 points in 2 dimensions.
 
  Read a triangulation dataset of TRIANGLE_NUM1
  triangles using 3 or 6 nodes per triangle.
 
  Subdivide each triangle into 4 triangles,
  generate new nodes as midpoints of current nodes.
 
  Write out the new node and triangulation data.
 
  If the input triangulation was Delaunay, then
  the output triangulation will be Delaunay.
 
  Read the header of "order3_nodes.txt".
 
  Spatial dimension DIM_NUM =        2
  Number of points NODE_NUM1  =       20
 
  Read the data in "order3_nodes.txt".
 
  First 5 nodes:
 
  Row         1             2       
  Col
    1    0.00000       0.00000    
    2    1.00000       0.00000    
    3    2.00000       0.00000    
    4    3.00000       0.00000    
    5    4.00000       0.00000    
 
 
  Read the header of "order3_triangles.txt".
 
  Triangle order =        3
  Number of triangles TRIANGLE_NUM1  =       24
 
  Read the data in "order3_triangles.txt".
 
  First 5 triangles:
 
  Row       1      2      3
  Col
 
    1       6      1      2
    2       7      6      2
    3       7     11      6
    4      12     11      7
    5      16     11     12
 
 
  Number of nodes in refined mesh =            63
  Number of triangles in refined mesh =        96
 
  First 5 output nodes:
 
  Row         1             2       
  Col
    1    0.00000       0.00000    
    2    1.00000       0.00000    
    3    2.00000       0.00000    
    4    3.00000       0.00000    
    5    4.00000       0.00000    
 
 
  First 5 output triangles
 
  Row       1      2      3
  Col
 
    1       6     22     24
    2      22      1     21
    3      24     21      2
    4      21     24     22
    5       7     34     25
 
 
TRIANGULATION_REFINE
  Normal end of execution.
 
January 29 2007   7:38:13.500 AM