January 16 2007 9:48:16.618 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 an associated triangulation dataset of TRIANGLE_NUM1 triangles which uses 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 "order6_nodes.txt". Spatial dimension DIM_NUM = 2 Number of points NODE_NUM1 = 37 Read the data in "order6_nodes.txt". First 5 nodes: Row 1 2 Col 1 0.00000 0.00000 2 0.00000 1.00000 3 0.00000 2.00000 4 0.00000 3.00000 5 0.00000 4.00000 Read the header of "order6_triangles.txt". Triangle order = 6 Number of triangles TRIANGLE_NUM1 = 10 Read the data in "order6_triangles.txt". First 5 triangles: Row 1 2 3 4 5 6 Col 1 1 12 3 7 8 2 2 1 9 12 6 11 7 3 3 15 5 12 13 4 4 18 5 15 14 13 17 5 9 33 11 19 20 10 Number of boundary nodes = 20 Number of interior nodes = 5 Number of edges = 25 New number of nodes = 117 New number of triangles = 40 Triangle_neighbors Row 1 2 3 Col 1 2 -1 -1 2 -1 -1 1 3 -1 4 -1 4 -1 3 -1 5 -1 6 -1 6 -1 5 -1 7 -1 -1 8 8 7 -1 -1 9 -1 10 -1 10 -1 9 -1 38 0.250000 0.500000 39 0.750000 1.50000 40 0.750000 2.00000 41 0.250000 2.00000 42 0.00000 1.50000 43 0.00000 0.500000 44 0.500000 1.50000 45 0.250000 1.50000 46 0.250000 1.00000 47 0.250000 0.00000 48 0.750000 0.00000 49 1.00000 0.500000 50 1.00000 1.50000 51 0.750000 0.500000 52 0.750000 1.00000 53 0.500000 0.500000 54 0.500000 2.00000 55 1.50000 2.00000 56 1.50000 2.50000 57 0.500000 3.50000 58 0.00000 3.50000 59 0.00000 2.50000 60 1.00000 2.50000 61 0.500000 3.00000 62 0.500000 2.50000 63 1.50000 4.00000 64 0.500000 4.00000 65 2.00000 2.50000 66 2.00000 3.50000 67 1.00000 3.50000 68 1.50000 3.00000 69 1.50000 3.50000 70 1.75000 0.00000 71 3.25000 0.00000 72 3.25000 0.250000 73 1.75000 0.750000 74 1.00000 0.750000 75 1.00000 0.250000 76 2.50000 0.250000 77 1.75000 0.500000 78 1.75000 0.250000 79 3.25000 1.00000 80 1.75000 1.00000 81 4.00000 0.250000 82 4.00000 0.750000 83 2.50000 0.750000 84 3.25000 0.500000 85 3.25000 0.750000 86 2.25000 2.00000 87 2.75000 2.00000 88 3.00000 2.25000 89 3.00000 2.75000 90 2.75000 2.75000 91 2.25000 2.25000 92 2.75000 2.25000 93 2.75000 2.50000 94 2.50000 2.25000 95 2.75000 3.00000 96 2.25000 3.00000 97 2.00000 2.75000 98 2.00000 2.25000 99 2.50000 2.75000 100 2.25000 2.75000 101 2.25000 2.50000 102 3.25000 1.00000 103 3.75000 1.00000 104 3.75000 1.75000 105 3.25000 3.25000 106 3.00000 3.25000 107 3.00000 1.75000 108 3.50000 1.75000 109 3.25000 2.50000 110 3.25000 1.75000 111 3.75000 4.00000 112 3.25000 4.00000 113 4.00000 1.75000 114 4.00000 3.25000 115 3.50000 3.25000 116 3.75000 2.50000 117 3.75000 3.25000 TRIANGLE_NODE2 Row 1 2 3 4 5 6 Col 1 1 7 2 38 46 43 2 12 8 7 40 44 39 3 3 2 8 42 45 41 4 7 8 2 44 45 46 5 1 6 7 47 53 38 6 9 11 6 49 51 48 7 12 7 11 39 52 50 8 6 11 7 51 52 53 9 3 12 4 54 62 59 10 15 13 12 56 60 55 11 5 4 13 58 61 57 12 12 13 4 60 61 62 13 18 14 17 63 69 66 14 5 13 14 57 67 64 15 15 17 13 65 68 56 16 14 13 17 67 68 69 17 9 19 10 70 78 75 18 33 20 19 72 76 71 19 11 10 20 74 77 73 20 19 20 10 76 77 78 21 35 21 34 79 85 82 22 11 20 21 73 83 80 23 33 34 20 81 84 72 24 21 20 34 83 84 85 25 15 22 23 86 94 91 26 26 27 22 88 92 87 27 28 23 27 90 93 89 28 22 27 23 92 93 94 29 15 23 16 91 101 98 30 28 24 23 95 99 90 31 17 16 24 97 100 96 32 23 24 16 99 100 101 33 25 30 27 102 110 107 34 35 31 30 104 108 103 35 29 27 31 106 109 105 36 30 31 27 108 109 110 37 37 32 36 111 117 114 38 29 31 32 105 115 112 39 35 36 31 113 116 104 40 32 31 36 115 116 117 NODE_XY2: Row 1 2 Col 1 0.00000 0.00000 2 0.00000 1.00000 3 0.00000 2.00000 4 0.00000 3.00000 5 0.00000 4.00000 6 0.500000 0.00000 7 0.500000 1.00000 8 0.500000 2.00000 9 1.00000 0.00000 10 1.00000 0.500000 11 1.00000 1.00000 12 1.00000 2.00000 13 1.00000 3.00000 14 1.00000 4.00000 15 2.00000 2.00000 16 2.00000 2.50000 17 2.00000 3.00000 18 2.00000 4.00000 19 2.50000 0.00000 20 2.50000 0.500000 21 2.50000 1.00000 22 2.50000 2.00000 23 2.50000 2.50000 24 2.50000 3.00000 25 3.00000 1.00000 26 3.00000 2.00000 27 3.00000 2.50000 28 3.00000 3.00000 29 3.00000 4.00000 30 3.50000 1.00000 31 3.50000 2.50000 32 3.50000 4.00000 33 4.00000 0.00000 34 4.00000 0.500000 35 4.00000 1.00000 36 4.00000 2.50000 37 4.00000 4.00000 38 0.250000 0.500000 39 0.750000 1.50000 40 0.750000 2.00000 41 0.250000 2.00000 42 0.00000 1.50000 43 0.00000 0.500000 44 0.500000 1.50000 45 0.250000 1.50000 46 0.250000 1.00000 47 0.250000 0.00000 48 0.750000 0.00000 49 1.00000 0.500000 50 1.00000 1.50000 51 0.750000 0.500000 52 0.750000 1.00000 53 0.500000 0.500000 54 0.500000 2.00000 55 1.50000 2.00000 56 1.50000 2.50000 57 0.500000 3.50000 58 0.00000 3.50000 59 0.00000 2.50000 60 1.00000 2.50000 61 0.500000 3.00000 62 0.500000 2.50000 63 1.50000 4.00000 64 0.500000 4.00000 65 2.00000 2.50000 66 2.00000 3.50000 67 1.00000 3.50000 68 1.50000 3.00000 69 1.50000 3.50000 70 1.75000 0.00000 71 3.25000 0.00000 72 3.25000 0.250000 73 1.75000 0.750000 74 1.00000 0.750000 75 1.00000 0.250000 76 2.50000 0.250000 77 1.75000 0.500000 78 1.75000 0.250000 79 3.25000 1.00000 80 1.75000 1.00000 81 4.00000 0.250000 82 4.00000 0.750000 83 2.50000 0.750000 84 3.25000 0.500000 85 3.25000 0.750000 86 2.25000 2.00000 87 2.75000 2.00000 88 3.00000 2.25000 89 3.00000 2.75000 90 2.75000 2.75000 91 2.25000 2.25000 92 2.75000 2.25000 93 2.75000 2.50000 94 2.50000 2.25000 95 2.75000 3.00000 96 2.25000 3.00000 97 2.00000 2.75000 98 2.00000 2.25000 99 2.50000 2.75000 100 2.25000 2.75000 101 2.25000 2.50000 102 3.25000 1.00000 103 3.75000 1.00000 104 3.75000 1.75000 105 3.25000 3.25000 106 3.00000 3.25000 107 3.00000 1.75000 108 3.50000 1.75000 109 3.25000 2.50000 110 3.25000 1.75000 111 3.75000 4.00000 112 3.25000 4.00000 113 4.00000 1.75000 114 4.00000 3.25000 115 3.50000 3.25000 116 3.75000 2.50000 117 3.75000 3.25000 TRIANGULATION_REFINE Normal end of execution. January 16 2007 9:48:16.849 AM