28 March 2023 6:29:44.411 PM stripack_test(): FORTRAN90 version Test stripack(). TEST01 TRANS converts Cartesian to spherical coordinates. TRMESH creates a triangulation. TRPRNT prints out a triangulation. TRLIST creates a triangle list. TRLPRT prints a triangle list. TRPLOT plots a triangulation. AREAS computes areas. BNODES computes boundary nodes. GETNP gets the next nearest node to a given node. NEARND returns the nearest node to a given point. DELARC removes a boundary arc if possible. CRLIST constructs the Voronoi diagram. VRPLOT plots the Voronoi diagram. SCOORD prints the Voronoi region boundary associated with a point. INSIDE determines if a point is inside a Voronoi region. I RLON RLAT 1 -43.386615 -4.089131 2 ********** 42.377641 3 175.179041 -5.711186 4 176.790641 -34.226284 5 -52.440530 81.438552 I X Y 1 -0.757239 -0.071369 2 -2.222367 0.739629 3 3.057451 -0.099679 4 3.085579 -0.597361 5 -0.915260 1.421371 I X Y Z 1 0.724885 -0.685169 -0.071308 2 -0.447986 -0.587379 0.674014 3 -0.991516 0.083625 -0.099514 4 -0.825526 0.046289 -0.562463 5 0.090749 -0.118012 0.988857 STRIPACK triangulation data structure, n = 100 Node Longitude Latitude neighbors of node 1 -0.433866E+02 -0.408913E+01 19 99 61 27 24 2 -0.127332E+03 0.423776E+02 56 15 73 11 3 0.175179E+03 -0.571119E+01 85 91 81 46 58 4 0.176791E+03 -0.342263E+02 81 76 67 48 46 5 -0.524405E+02 0.814386E+02 72 79 59 37 77 66 6 0.137011E+03 0.683895E+02 45 47 82 49 97 31 7 -0.136060E+03 0.747595E+02 69 66 68 18 45 31 8 -0.731043E+02 -0.635081E+02 54 94 70 87 24 83 16 9 0.124891E+03 -0.829969E+02 34 71 36 13 40 10 10 0.176635E+03 -0.816425E+02 40 54 48 34 9 11 -0.115358E+03 0.367038E+02 73 53 20 52 56 2 12 -0.143034E+01 0.248271E+02 19 65 32 60 25 44 13 0.620730E+02 -0.852115E+02 36 70 43 94 40 9 14 0.131860E+03 0.136293E+01 97 86 74 17 95 91 15 -0.145589E+03 0.230576E+02 2 56 51 85 92 73 16 -0.107751E+03 -0.650564E+02 83 96 48 54 8 17 0.102113E+03 -0.108687E+02 74 42 78 63 80 95 14 18 0.153403E+03 0.812319E+02 7 68 50 47 45 19 -0.243381E+02 -0.560500E+01 24 89 65 12 44 99 1 20 -0.105802E+03 0.335619E+02 53 28 61 64 72 52 11 21 0.143437E+03 -0.238771E+02 95 80 29 81 22 -0.108712E+03 0.259671E+02 53 73 23 28 23 -0.977807E+02 0.232156E+01 73 38 75 27 61 28 22 24 -0.507845E+02 -0.271745E+02 27 75 96 83 8 87 89 19 1 25 0.290529E+02 0.432011E+02 60 41 35 62 57 44 12 26 0.449208E+02 0.802740E+02 59 62 100 50 37 27 -0.585453E+02 -0.111325E+02 61 23 75 24 1 28 -0.104759E+03 0.315485E+02 53 22 23 61 20 29 0.141761E+03 -0.277277E+02 80 76 81 21 30 0.500986E+01 -0.399881E+02 84 65 89 87 55 39 31 0.171695E+03 0.520593E+02 97 51 56 52 69 7 45 6 32 0.300540E+02 -0.124385E+02 84 39 42 98 41 60 12 65 33 0.712078E+02 0.272081E+02 93 74 86 35 34 0.156955E+03 -0.793237E+02 48 67 71 9 10 35 0.650617E+02 0.487396E+02 86 49 82 100 62 25 41 93 33 36 0.944912E+02 -0.604609E+02 9 71 78 42 39 55 70 13 37 -0.537675E+02 0.833130E+02 59 26 50 77 5 38 -0.119687E+03 -0.163459E+02 73 92 90 75 23 39 0.246397E+02 -0.349206E+02 55 36 42 32 84 30 40 -0.177281E+03 -0.890681E+02 94 54 10 9 13 41 0.486300E+02 0.256903E+02 35 25 60 32 98 93 42 0.733772E+02 -0.133811E+02 74 98 32 39 36 78 17 43 0.173244E+01 -0.791272E+02 70 94 13 44 -0.406140E+01 0.318596E+02 25 57 99 19 12 45 0.140445E+03 0.741470E+02 31 7 18 47 6 46 -0.179144E+03 -0.311665E+02 48 90 92 58 3 81 4 47 0.135962E+03 0.743656E+02 18 50 100 82 6 45 48 -0.155651E+03 -0.645219E+02 54 16 96 75 90 46 4 67 34 10 49 0.109877E+03 0.540149E+02 82 35 86 97 6 50 0.795949E+02 0.846935E+02 100 47 18 68 77 37 26 51 0.174147E+03 0.420675E+02 97 91 85 15 56 31 52 -0.124682E+03 0.654808E+02 31 56 11 20 72 66 69 53 -0.110591E+03 0.295325E+02 28 20 11 73 22 54 -0.111360E+03 -0.780964E+02 94 8 16 48 10 40 55 0.131296E+02 -0.482870E+02 87 70 36 39 30 56 -0.130452E+03 0.504373E+02 11 52 31 51 15 2 57 -0.548867E+01 0.558985E+02 25 62 79 72 99 44 58 -0.161331E+03 -0.102693E+02 46 92 85 3 59 -0.145281E+02 0.766300E+02 79 62 26 37 5 60 0.415050E+02 0.254981E+02 25 12 32 41 61 -0.707562E+02 0.336386E+02 23 27 1 99 88 64 20 28 62 0.185394E+01 0.642706E+02 25 35 100 26 59 79 57 63 0.119309E+03 -0.531741E+02 71 67 76 80 17 78 64 -0.672837E+02 0.415322E+02 88 72 20 61 65 0.280801E+01 -0.120214E+02 12 19 89 30 84 32 66 -0.837802E+02 0.824276E+02 77 68 7 69 52 72 5 67 0.157284E+03 -0.527688E+02 71 34 48 4 76 63 68 0.165457E+03 0.825487E+02 7 66 77 50 18 69 -0.128878E+03 0.705323E+02 31 52 66 7 70 0.454945E+01 -0.749018E+02 13 36 55 87 8 94 43 71 0.113084E+03 -0.620449E+02 9 34 67 63 78 36 72 -0.602036E+02 0.576893E+02 88 99 57 79 5 66 52 20 64 73 -0.134980E+03 0.160191E+02 92 38 23 22 53 11 2 15 74 0.873713E+02 -0.384519E+01 14 86 33 93 98 42 17 75 -0.111572E+03 -0.433350E+02 96 24 27 23 38 90 48 76 0.139712E+03 -0.376096E+02 80 63 67 4 81 29 77 -0.803563E+02 0.830426E+02 5 37 50 68 66 78 0.104552E+03 -0.584122E+02 71 63 17 42 36 79 -0.370814E+02 0.676424E+02 5 72 57 62 59 80 0.134871E+03 -0.285000E+02 21 95 17 63 76 29 81 0.160829E+03 -0.224053E+02 29 76 4 46 3 91 95 21 82 0.107574E+03 0.590554E+02 6 47 100 35 49 83 -0.100605E+03 -0.552484E+02 8 24 96 16 84 0.173872E+02 -0.170729E+02 65 30 39 32 85 -0.162501E+03 0.204309E+02 51 91 3 58 92 15 86 0.787442E+02 0.347516E+02 74 14 97 49 35 33 87 -0.101151E+02 -0.547193E+02 55 30 89 24 8 70 88 -0.602886E+02 0.386926E+02 61 99 72 64 89 -0.109405E+02 -0.338914E+02 65 19 24 87 30 90 -0.122575E+03 -0.220995E+02 92 46 48 75 38 91 0.166467E+03 0.567734E+01 85 51 97 14 95 81 3 92 -0.154504E+03 -0.834936E+01 46 90 38 73 15 85 58 93 0.683831E+02 0.262216E+02 35 41 98 74 33 94 -0.463472E+02 -0.825261E+02 70 8 54 40 13 43 95 0.136195E+03 -0.134407E+02 81 91 14 17 80 21 96 -0.104803E+03 -0.523235E+02 24 75 48 16 83 97 0.135546E+03 0.453353E+02 86 14 91 51 31 6 49 98 0.595731E+02 0.228565E+02 41 32 42 74 93 99 -0.506777E+02 0.366477E+02 61 1 19 44 57 72 88 100 0.581384E+02 0.784635E+02 62 35 82 47 50 26 NB = 0 boundary arcs. NA = 294 arcs. NT = 196 triangles. STRIPACK (TRLIST) output, n = 100 Node Longitude Latitude 1 -0.433866E+02 -0.408913E+01 2 -0.127332E+03 0.423776E+02 3 0.175179E+03 -0.571119E+01 4 0.176791E+03 -0.342263E+02 5 -0.524405E+02 0.814386E+02 6 0.137011E+03 0.683895E+02 7 -0.136060E+03 0.747595E+02 8 -0.731043E+02 -0.635081E+02 9 0.124891E+03 -0.829969E+02 10 0.176635E+03 -0.816425E+02 11 -0.115358E+03 0.367038E+02 12 -0.143034E+01 0.248271E+02 13 0.620730E+02 -0.852115E+02 14 0.131860E+03 0.136293E+01 15 -0.145589E+03 0.230576E+02 16 -0.107751E+03 -0.650564E+02 17 0.102113E+03 -0.108687E+02 18 0.153403E+03 0.812319E+02 19 -0.243381E+02 -0.560500E+01 20 -0.105802E+03 0.335619E+02 21 0.143437E+03 -0.238771E+02 22 -0.108712E+03 0.259671E+02 23 -0.977807E+02 0.232156E+01 24 -0.507845E+02 -0.271745E+02 25 0.290529E+02 0.432011E+02 26 0.449208E+02 0.802740E+02 27 -0.585453E+02 -0.111325E+02 28 -0.104759E+03 0.315485E+02 29 0.141761E+03 -0.277277E+02 30 0.500986E+01 -0.399881E+02 31 0.171695E+03 0.520593E+02 32 0.300540E+02 -0.124385E+02 33 0.712078E+02 0.272081E+02 34 0.156955E+03 -0.793237E+02 35 0.650617E+02 0.487396E+02 36 0.944912E+02 -0.604609E+02 37 -0.537675E+02 0.833130E+02 38 -0.119687E+03 -0.163459E+02 39 0.246397E+02 -0.349206E+02 40 -0.177281E+03 -0.890681E+02 41 0.486300E+02 0.256903E+02 42 0.733772E+02 -0.133811E+02 43 0.173244E+01 -0.791272E+02 44 -0.406140E+01 0.318596E+02 45 0.140445E+03 0.741470E+02 46 -0.179144E+03 -0.311665E+02 47 0.135962E+03 0.743656E+02 48 -0.155651E+03 -0.645219E+02 49 0.109877E+03 0.540149E+02 50 0.795949E+02 0.846935E+02 51 0.174147E+03 0.420675E+02 52 -0.124682E+03 0.654808E+02 53 -0.110591E+03 0.295325E+02 54 -0.111360E+03 -0.780964E+02 55 0.131296E+02 -0.482870E+02 56 -0.130452E+03 0.504373E+02 57 -0.548867E+01 0.558985E+02 58 -0.161331E+03 -0.102693E+02 59 -0.145281E+02 0.766300E+02 60 0.415050E+02 0.254981E+02 61 -0.707562E+02 0.336386E+02 62 0.185394E+01 0.642706E+02 63 0.119309E+03 -0.531741E+02 64 -0.672837E+02 0.415322E+02 65 0.280801E+01 -0.120214E+02 66 -0.837802E+02 0.824276E+02 67 0.157284E+03 -0.527688E+02 68 0.165457E+03 0.825487E+02 69 -0.128878E+03 0.705323E+02 70 0.454945E+01 -0.749018E+02 71 0.113084E+03 -0.620449E+02 72 -0.602036E+02 0.576893E+02 73 -0.134980E+03 0.160191E+02 74 0.873713E+02 -0.384519E+01 75 -0.111572E+03 -0.433350E+02 76 0.139712E+03 -0.376096E+02 77 -0.803563E+02 0.830426E+02 78 0.104552E+03 -0.584122E+02 79 -0.370814E+02 0.676424E+02 80 0.134871E+03 -0.285000E+02 81 0.160829E+03 -0.224053E+02 82 0.107574E+03 0.590554E+02 83 -0.100605E+03 -0.552484E+02 84 0.173872E+02 -0.170729E+02 85 -0.162501E+03 0.204309E+02 86 0.787442E+02 0.347516E+02 87 -0.101151E+02 -0.547193E+02 88 -0.602886E+02 0.386926E+02 89 -0.109405E+02 -0.338914E+02 90 -0.122575E+03 -0.220995E+02 91 0.166467E+03 0.567734E+01 92 -0.154504E+03 -0.834936E+01 93 0.683831E+02 0.262216E+02 94 -0.463472E+02 -0.825261E+02 95 0.136195E+03 -0.134407E+02 96 -0.104803E+03 -0.523235E+02 97 0.135546E+03 0.453353E+02 98 0.595731E+02 0.228565E+02 99 -0.506777E+02 0.366477E+02 100 0.581384E+02 0.784635E+02 triangle vertices neighbors arcs kt n1 n2 n3 kt1 kt2 kt3 ka1 ka2 ka3 1 1 19 99 91 2 5 103 1 4 2 1 99 61 186 3 1 269 2 1 3 1 61 27 108 4 2 125 3 2 4 1 27 24 110 5 3 131 5 3 5 1 24 19 89 1 4 100 4 5 6 2 56 15 74 7 9 79 6 8 7 2 15 73 77 8 6 82 7 6 8 2 73 11 54 9 7 57 9 7 9 2 11 56 57 6 8 60 8 9 10 3 85 91 177 11 14 247 10 13 11 3 91 81 196 12 10 294 11 10 12 3 81 46 19 13 11 18 12 11 13 3 46 58 169 14 12 228 14 12 14 3 58 85 184 10 13 264 13 14 15 4 81 76 125 16 19 151 15 19 16 4 76 67 189 17 15 273 16 15 17 4 67 48 146 18 16 184 17 16 18 4 48 46 167 19 17 225 20 17 19 4 46 81 12 15 18 18 19 20 20 5 72 79 182 21 25 259 21 25 21 5 79 59 185 22 20 268 22 21 22 5 59 37 123 23 21 147 23 22 23 5 37 77 158 24 22 206 24 23 24 5 77 66 193 25 23 285 26 24 25 5 66 72 178 20 24 250 25 26 26 6 45 47 88 27 31 97 27 31 27 6 47 82 171 28 26 233 28 27 28 6 82 49 149 29 27 189 29 28 29 6 49 97 174 30 28 241 30 29 30 6 97 31 132 31 29 162 32 30 31 6 31 45 36 26 30 36 31 32 32 7 69 66 179 33 37 252 33 38 33 7 66 68 193 34 32 284 34 33 34 7 68 18 86 35 33 95 35 34 35 7 18 45 88 36 34 98 37 35 36 7 45 31 31 37 35 36 39 37 37 7 31 69 135 32 36 167 38 39 38 8 54 94 162 39 44 213 40 45 39 8 94 70 165 40 38 220 41 40 40 8 70 87 180 41 39 255 42 41 41 8 87 24 113 42 40 137 43 42 42 8 24 83 112 43 41 134 44 43 43 8 83 16 78 44 42 84 46 44 44 8 16 54 80 38 43 87 45 46 45 9 34 71 147 46 50 186 47 51 46 9 71 36 153 47 45 197 48 47 47 9 36 13 64 48 46 68 49 48 48 9 13 40 67 49 47 71 50 49 49 9 40 10 51 50 48 53 52 50 50 9 10 34 53 45 49 55 51 52 51 10 40 54 162 52 49 214 54 53 52 10 54 48 80 53 51 86 56 54 53 10 48 34 146 50 52 185 55 56 54 11 73 53 101 55 8 116 58 57 55 11 53 20 92 56 54 105 59 58 56 11 20 52 96 57 55 109 61 59 57 11 52 56 134 9 56 165 60 61 58 12 19 65 90 59 63 101 62 66 59 12 65 32 141 60 58 178 63 62 60 12 32 60 140 61 59 174 64 63 61 12 60 25 114 62 60 138 65 64 62 12 25 44 118 63 61 142 67 65 63 12 44 19 91 58 62 104 66 67 64 13 36 70 157 65 47 204 69 68 65 13 70 43 165 66 64 222 70 69 66 13 43 94 165 67 65 221 72 70 67 13 94 40 162 48 66 215 71 72 68 14 97 86 174 69 73 240 73 77 69 14 86 74 143 70 68 179 74 73 70 14 74 17 81 71 69 89 75 74 71 14 17 95 85 72 70 93 76 75 72 14 95 91 196 73 71 292 78 76 73 14 91 97 176 68 72 245 77 78 74 15 56 51 133 75 6 163 80 79 75 15 51 85 177 76 74 248 81 80 76 15 85 92 184 77 75 263 83 81 77 15 92 73 159 7 76 208 82 83 78 16 83 96 112 79 43 133 85 84 79 16 96 48 172 80 78 236 88 85 80 16 48 54 52 44 79 86 87 88 81 17 74 42 164 82 70 219 90 89 82 17 42 78 154 83 81 198 91 90 83 17 78 63 191 84 82 280 92 91 84 17 63 80 190 85 83 276 94 92 85 17 80 95 97 71 84 111 93 94 86 18 68 50 175 87 34 244 96 95 87 18 50 47 170 88 86 231 99 96 88 18 47 45 26 35 87 97 98 99 89 19 24 89 113 90 5 136 102 100 90 19 89 65 127 58 89 154 101 102 91 19 44 99 166 1 63 223 103 104 92 20 53 28 104 93 55 119 106 105 93 20 28 61 109 94 92 127 107 106 94 20 61 64 187 95 93 270 108 107 95 20 64 72 192 96 94 281 110 108 96 20 72 52 178 56 95 251 109 110 97 21 95 80 85 98 100 111 112 114 98 21 80 29 124 99 97 150 113 112 99 21 29 81 125 100 98 152 115 113 100 21 81 95 196 97 99 293 114 115 101 22 53 73 54 102 104 116 117 120 102 22 73 23 105 103 101 122 118 117 103 22 23 28 109 104 102 128 121 118 104 22 28 53 92 101 103 119 120 121 105 23 73 38 159 106 102 209 123 122 106 23 38 75 161 107 105 211 124 123 107 23 75 27 110 108 106 130 126 124 108 23 27 61 3 109 107 125 129 126 109 23 61 28 93 103 108 127 128 129 110 24 27 75 107 111 4 130 132 131 111 24 75 96 172 112 110 235 135 132 112 24 96 83 78 42 111 133 134 135 113 24 87 89 128 89 41 156 136 137 114 25 60 41 140 115 61 173 139 138 115 25 41 35 152 116 114 196 140 139 116 25 35 62 151 117 115 193 141 140 117 25 62 57 181 118 116 258 143 141 118 25 57 44 166 62 117 224 142 143 119 26 59 62 185 120 123 267 144 148 120 26 62 100 151 121 119 192 145 144 121 26 100 50 170 122 120 230 146 145 122 26 50 37 158 123 121 207 149 146 123 26 37 59 22 119 122 147 148 149 124 29 80 76 190 125 98 275 153 150 125 29 76 81 15 99 124 151 152 153 126 30 84 65 141 127 131 176 155 160 127 30 65 89 90 128 126 154 157 155 128 30 89 87 113 129 127 156 158 157 129 30 87 55 180 130 128 257 159 158 130 30 55 39 156 131 129 202 161 159 131 30 39 84 136 126 130 169 160 161 132 31 97 51 176 133 30 246 164 162 133 31 51 56 74 134 132 163 166 164 134 31 56 52 57 135 133 165 168 166 135 31 52 69 179 37 134 253 167 168 136 32 84 39 131 137 141 169 170 177 137 32 39 42 155 138 136 200 171 170 138 32 42 98 164 139 137 218 172 171 139 32 98 41 163 140 138 217 175 172 140 32 41 60 114 60 139 173 174 175 141 32 65 84 126 136 59 176 177 178 142 33 93 74 195 143 145 291 180 182 143 33 74 86 69 144 142 179 181 180 144 33 86 35 148 145 143 188 183 181 145 33 35 93 152 142 144 195 182 183 146 34 48 67 17 147 53 184 187 185 147 34 67 71 188 45 146 272 186 187 148 35 86 49 174 149 144 242 190 188 149 35 49 82 28 150 148 189 191 190 150 35 82 100 171 151 149 232 194 191 151 35 100 62 120 116 150 192 193 194 152 35 41 93 163 145 115 216 195 196 153 36 71 78 191 154 46 278 199 197 154 36 78 42 82 155 153 198 201 199 155 36 42 39 137 156 154 200 203 201 156 36 39 55 130 157 155 202 205 203 157 36 55 70 180 64 156 256 204 205 158 37 50 77 175 23 122 243 206 207 159 38 73 92 77 160 105 208 210 209 160 38 92 90 168 161 159 226 212 210 161 38 90 75 173 106 160 237 211 212 162 40 94 54 38 51 67 213 214 215 163 41 98 93 195 152 139 289 216 217 164 42 74 98 195 138 81 290 218 219 165 43 70 94 39 66 65 220 221 222 166 44 57 99 183 91 118 261 223 224 167 46 48 90 173 168 18 238 227 225 168 46 90 92 160 169 167 226 229 227 169 46 92 58 184 13 168 265 228 229 170 47 50 100 121 171 87 230 234 231 171 47 100 82 150 27 170 232 233 234 172 48 96 75 111 173 79 235 239 236 173 48 75 90 161 167 172 237 238 239 174 49 86 97 68 29 148 240 241 242 175 50 68 77 193 158 86 283 243 244 176 51 97 91 73 177 132 245 249 246 177 51 91 85 10 75 176 247 248 249 178 52 72 66 25 179 96 250 254 251 179 52 66 69 32 135 178 252 253 254 180 55 87 70 40 157 129 255 256 257 181 57 62 79 185 182 117 266 260 258 182 57 79 72 20 183 181 259 262 260 183 57 72 99 194 166 182 287 261 262 184 58 92 85 76 14 169 263 264 265 185 59 79 62 181 119 21 266 267 268 186 61 99 88 194 187 2 286 271 269 187 61 88 64 192 94 186 282 270 271 188 63 71 67 147 189 191 272 274 279 189 63 67 76 16 190 188 273 277 274 190 63 76 80 124 84 189 275 276 277 191 63 78 71 153 188 83 278 279 280 192 64 88 72 194 95 187 288 281 282 193 66 77 68 175 33 24 283 284 285 194 72 88 99 186 183 192 286 287 288 195 74 93 98 163 164 142 289 290 291 196 81 91 95 72 100 11 292 293 294 Number of boundary nodes NB = 0 Number of arcs NA = 294 Number of triangles NT = 196 TRPLOT created the triangulation plot file: "stripack_test_del.eps". Relative area of convex hull = 1.00 Output from BNODES: Number of boundary nodes = 0 Number of arcs = 294 Number of triangles = 196 Subroutine DELARC was not tested. Nodes 99 and 100 do not form a removable boundary arc. VRPLOT created the Voronoi plot file: "stripack_test_vor.eps". Voronoi region for node 1 Triangle Latitude Longitude Circumradius 2 18.051045 -31.986187 0.433376 3 14.441190 -56.966266 0.399766 4 12.235949 -60.298739 0.408986 5 -15.619603 -47.021386 0.210693 1 -18.728105 -35.032476 0.292542 INSIDE correctly reports that node N0 is inside its Voronoi region! TEST02 EDGE can be used to force an edge between two nodes. DEL can be used to delete a node. EDGE has forced an edge between two nodes. Call DELNOD to delete node 9 Call DELNOD to delete node 8 Call DELNOD to delete node 7 Call DELNOD to delete node 6 Call DELNOD to delete node 5 Call DELNOD to delete node 4 stripack_test(): Normal end of execution. 28 March 2023 6:29:44.414 PM