09 May 2025 9:46:09.488 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 66.783542 -48.489571 2 129.166149 -7.860773 3 -44.177537 -58.857949 4 13.357664 -19.667317 5 174.068995 -49.827337 I X Y 1 1.165593 -0.846303 2 2.254375 -0.137196 3 -0.771043 -1.027265 4 0.233135 -0.343259 5 3.038077 -0.869651 I X Y Z 1 0.261262 0.609088 -0.748835 2 -0.625637 0.768032 -0.136766 3 0.370900 -0.360402 -0.855888 4 0.916188 0.217551 -0.336558 5 -0.641640 0.066658 -0.764104 STRIPACK triangulation data structure, n = 100 Node Longitude Latitude neighbors of node 1 0.667835E+02 -0.484896E+02 50 11 65 22 26 67 2 0.129166E+03 -0.786077E+01 60 54 55 75 23 3 -0.441775E+02 -0.588579E+02 62 27 45 38 4 0.133577E+02 -0.196673E+02 62 92 96 11 6 44 5 0.174069E+03 -0.498273E+02 81 67 51 71 90 6 0.309913E+02 0.114446E+02 14 13 44 4 11 7 -0.135919E+03 -0.512025E+02 51 79 42 43 29 71 8 0.867905E+02 0.728875E+02 85 74 86 9 -0.174193E+03 -0.163631E+02 90 71 29 69 82 10 -0.115322E+03 0.313872E+02 34 32 33 19 88 29 11 0.296786E+02 -0.235304E+02 96 65 1 50 14 6 4 12 0.362239E+02 -0.666181E+02 52 26 21 64 13 0.302907E+02 0.248295E+02 14 84 53 44 6 14 0.498420E+02 0.156744E+02 11 50 94 84 13 6 15 0.155553E+03 0.434495E+02 68 41 39 36 16 35 46 16 0.143308E+03 0.528206E+02 36 24 57 35 15 17 -0.866100E+02 0.647460E+02 33 47 25 97 63 56 19 18 0.163487E+03 -0.236635E+02 82 31 54 81 90 19 -0.133659E+03 0.567453E+02 56 28 39 41 88 10 33 17 20 -0.790998E+02 0.281930E+02 33 32 34 89 66 47 21 0.331934E+02 -0.624351E+02 26 22 92 38 64 12 22 0.344250E+02 -0.475313E+02 26 1 65 96 92 21 23 0.111573E+03 0.182395E+01 75 87 60 2 24 0.125433E+03 0.660724E+02 86 74 57 16 36 100 25 -0.230040E+02 0.564630E+02 47 66 44 78 97 17 26 0.553252E+02 -0.567703E+02 22 21 12 52 67 1 27 -0.601155E+02 -0.594432E+02 62 43 45 3 28 -0.156951E+03 0.718806E+02 56 49 100 39 19 29 -0.154389E+03 -0.478753E+01 71 7 43 34 10 88 69 9 30 -0.632827E+01 -0.869600E+02 64 91 93 77 95 31 0.158557E+03 -0.215194E+01 82 61 73 55 54 18 32 -0.803685E+02 0.282215E+02 33 10 34 20 33 -0.976161E+02 0.378871E+02 32 20 47 17 19 10 34 -0.839761E+02 -0.934118E+01 10 29 43 80 98 89 20 32 35 0.130896E+03 0.419419E+02 57 83 55 46 15 16 36 0.149370E+03 0.552947E+02 15 39 100 24 16 37 -0.105784E+03 -0.712472E+02 99 48 79 38 -0.381394E+02 -0.595902E+02 45 91 64 21 92 62 3 39 -0.177632E+03 0.658725E+02 100 36 15 41 19 28 40 0.153458E+03 -0.856277E+02 93 91 99 51 70 95 77 41 -0.163145E+03 0.482346E+02 68 88 19 39 15 42 -0.108624E+03 -0.534356E+02 79 48 43 7 43 -0.906439E+02 -0.336601E+02 48 45 27 62 80 34 29 7 42 44 -0.246949E+02 0.204353E+02 66 98 62 4 6 13 53 78 25 45 -0.673557E+02 -0.620868E+02 43 48 99 91 38 3 27 46 0.152037E+03 0.301519E+02 35 55 73 68 15 47 -0.514714E+02 0.496846E+02 17 33 20 66 25 48 -0.916602E+02 -0.618604E+02 99 45 43 42 79 37 49 -0.167572E+03 0.828902E+02 56 63 59 58 100 28 50 0.820929E+02 -0.299667E+02 67 81 60 87 94 14 11 1 51 -0.172414E+03 -0.685339E+02 99 79 7 71 5 67 70 40 52 0.559947E+02 -0.790248E+02 64 95 67 26 12 53 0.291075E+02 0.296747E+02 84 85 78 44 13 54 0.144298E+03 -0.135395E+02 60 81 18 31 55 2 55 0.142832E+03 0.153813E+02 35 83 75 2 54 31 73 46 56 -0.129633E+03 0.694934E+02 19 17 63 49 28 57 0.127682E+03 0.443228E+02 16 24 74 72 83 35 58 0.121705E+03 0.812241E+02 59 76 86 100 49 59 0.135893E+03 0.888661E+02 63 97 76 58 49 60 0.125329E+03 -0.171509E+02 23 87 50 81 54 2 61 0.161645E+03 -0.192034E+00 82 69 68 73 31 62 -0.484724E+02 -0.539756E+02 38 92 4 44 98 80 43 27 3 63 -0.135495E+03 0.851232E+02 56 17 97 59 49 64 -0.609883E+01 -0.769799E+02 95 52 12 21 38 91 30 65 0.295140E+02 -0.373895E+02 1 11 96 22 66 -0.523745E+02 0.366279E+02 20 89 98 44 25 47 67 0.143811E+03 -0.660909E+02 95 70 51 5 81 50 1 26 52 68 0.163929E+03 0.161950E+02 69 88 41 15 46 73 61 69 -0.165893E+03 0.631046E+01 88 68 61 82 9 29 70 0.152056E+03 -0.781002E+02 51 67 95 40 71 -0.163284E+03 -0.396871E+02 51 7 29 9 90 5 72 0.114497E+03 0.236696E+02 75 83 57 74 84 94 87 73 0.151498E+03 0.154724E+02 31 61 68 46 55 74 0.922727E+02 0.696095E+02 85 84 72 57 24 86 8 75 0.118493E+03 0.147957E+02 87 23 2 55 83 72 76 0.102074E+03 0.838695E+02 59 97 85 86 58 77 0.667420E+02 -0.879755E+02 93 40 95 30 78 -0.223275E+00 0.566115E+02 44 53 85 97 25 79 -0.128811E+03 -0.655965E+02 48 42 7 51 99 37 80 -0.799306E+02 -0.135021E+02 43 62 98 34 81 0.152750E+03 -0.450369E+02 90 18 54 60 50 67 5 82 0.171813E+03 -0.157201E+02 90 9 69 61 31 18 83 0.119083E+03 0.232300E+02 75 55 35 57 72 84 0.536273E+02 0.294998E+02 94 72 74 85 53 13 14 85 0.437528E+02 0.707211E+02 86 76 97 78 53 84 74 8 86 0.964560E+02 0.760676E+02 74 24 100 58 76 85 8 87 0.102389E+03 0.333907E+01 75 72 94 50 60 23 88 -0.152444E+03 0.284887E+02 68 69 29 10 19 41 89 -0.585651E+02 0.186871E+02 20 34 98 66 90 0.179390E+03 -0.356253E+02 9 82 18 81 5 71 91 -0.792181E+02 -0.809316E+02 64 38 45 99 40 93 30 92 0.154133E+02 -0.437359E+02 96 4 62 38 21 22 93 -0.861061E+02 -0.892455E+02 91 40 77 30 94 0.908418E+02 0.465063E+01 72 84 14 50 87 95 0.598406E+02 -0.804275E+02 64 30 77 40 70 67 52 96 0.251220E+02 -0.375177E+02 11 4 92 22 65 97 -0.155385E+02 0.826795E+02 63 17 25 78 85 76 59 98 -0.566256E+02 0.163150E+02 34 80 62 44 66 89 99 -0.972249E+02 -0.758879E+02 79 51 40 91 45 48 37 100 0.148684E+03 0.765223E+02 86 24 36 39 28 49 58 NB = 0 boundary arcs. NA = 294 arcs. NT = 196 triangles. STRIPACK (TRLIST) output, n = 100 Node Longitude Latitude 1 0.667835E+02 -0.484896E+02 2 0.129166E+03 -0.786077E+01 3 -0.441775E+02 -0.588579E+02 4 0.133577E+02 -0.196673E+02 5 0.174069E+03 -0.498273E+02 6 0.309913E+02 0.114446E+02 7 -0.135919E+03 -0.512025E+02 8 0.867905E+02 0.728875E+02 9 -0.174193E+03 -0.163631E+02 10 -0.115322E+03 0.313872E+02 11 0.296786E+02 -0.235304E+02 12 0.362239E+02 -0.666181E+02 13 0.302907E+02 0.248295E+02 14 0.498420E+02 0.156744E+02 15 0.155553E+03 0.434495E+02 16 0.143308E+03 0.528206E+02 17 -0.866100E+02 0.647460E+02 18 0.163487E+03 -0.236635E+02 19 -0.133659E+03 0.567453E+02 20 -0.790998E+02 0.281930E+02 21 0.331934E+02 -0.624351E+02 22 0.344250E+02 -0.475313E+02 23 0.111573E+03 0.182395E+01 24 0.125433E+03 0.660724E+02 25 -0.230040E+02 0.564630E+02 26 0.553252E+02 -0.567703E+02 27 -0.601155E+02 -0.594432E+02 28 -0.156951E+03 0.718806E+02 29 -0.154389E+03 -0.478753E+01 30 -0.632827E+01 -0.869600E+02 31 0.158557E+03 -0.215194E+01 32 -0.803685E+02 0.282215E+02 33 -0.976161E+02 0.378871E+02 34 -0.839761E+02 -0.934118E+01 35 0.130896E+03 0.419419E+02 36 0.149370E+03 0.552947E+02 37 -0.105784E+03 -0.712472E+02 38 -0.381394E+02 -0.595902E+02 39 -0.177632E+03 0.658725E+02 40 0.153458E+03 -0.856277E+02 41 -0.163145E+03 0.482346E+02 42 -0.108624E+03 -0.534356E+02 43 -0.906439E+02 -0.336601E+02 44 -0.246949E+02 0.204353E+02 45 -0.673557E+02 -0.620868E+02 46 0.152037E+03 0.301519E+02 47 -0.514714E+02 0.496846E+02 48 -0.916602E+02 -0.618604E+02 49 -0.167572E+03 0.828902E+02 50 0.820929E+02 -0.299667E+02 51 -0.172414E+03 -0.685339E+02 52 0.559947E+02 -0.790248E+02 53 0.291075E+02 0.296747E+02 54 0.144298E+03 -0.135395E+02 55 0.142832E+03 0.153813E+02 56 -0.129633E+03 0.694934E+02 57 0.127682E+03 0.443228E+02 58 0.121705E+03 0.812241E+02 59 0.135893E+03 0.888661E+02 60 0.125329E+03 -0.171509E+02 61 0.161645E+03 -0.192034E+00 62 -0.484724E+02 -0.539756E+02 63 -0.135495E+03 0.851232E+02 64 -0.609883E+01 -0.769799E+02 65 0.295140E+02 -0.373895E+02 66 -0.523745E+02 0.366279E+02 67 0.143811E+03 -0.660909E+02 68 0.163929E+03 0.161950E+02 69 -0.165893E+03 0.631046E+01 70 0.152056E+03 -0.781002E+02 71 -0.163284E+03 -0.396871E+02 72 0.114497E+03 0.236696E+02 73 0.151498E+03 0.154724E+02 74 0.922727E+02 0.696095E+02 75 0.118493E+03 0.147957E+02 76 0.102074E+03 0.838695E+02 77 0.667420E+02 -0.879755E+02 78 -0.223275E+00 0.566115E+02 79 -0.128811E+03 -0.655965E+02 80 -0.799306E+02 -0.135021E+02 81 0.152750E+03 -0.450369E+02 82 0.171813E+03 -0.157201E+02 83 0.119083E+03 0.232300E+02 84 0.536273E+02 0.294998E+02 85 0.437528E+02 0.707211E+02 86 0.964560E+02 0.760676E+02 87 0.102389E+03 0.333907E+01 88 -0.152444E+03 0.284887E+02 89 -0.585651E+02 0.186871E+02 90 0.179390E+03 -0.356253E+02 91 -0.792181E+02 -0.809316E+02 92 0.154133E+02 -0.437359E+02 93 -0.861061E+02 -0.892455E+02 94 0.908418E+02 0.465063E+01 95 0.598406E+02 -0.804275E+02 96 0.251220E+02 -0.375177E+02 97 -0.155385E+02 0.826795E+02 98 -0.566256E+02 0.163150E+02 99 -0.972249E+02 -0.758879E+02 100 0.148684E+03 0.765223E+02 triangle vertices neighbors arcs kt n1 n2 n3 kt1 kt2 kt3 ka1 ka2 ka3 1 1 50 11 51 2 6 57 1 5 2 1 11 65 50 3 1 54 2 1 3 1 65 22 97 4 2 117 3 2 4 1 22 26 93 5 3 110 4 3 5 1 26 67 109 6 4 136 6 4 6 1 67 50 162 1 5 222 5 6 7 2 60 54 172 8 11 242 7 10 8 2 54 55 125 9 7 157 8 7 9 2 55 75 173 10 8 243 9 8 10 2 75 23 99 11 9 121 11 9 11 2 23 60 100 7 10 122 10 11 12 3 62 27 110 13 15 138 12 14 13 3 27 45 111 14 12 139 13 12 14 3 45 38 136 15 13 174 15 13 15 3 38 62 138 12 14 179 14 15 16 4 62 92 138 17 21 178 16 20 17 4 92 96 98 18 16 118 17 16 18 4 96 11 50 19 17 55 18 17 19 4 11 6 29 20 18 30 19 18 20 4 6 44 28 21 19 27 21 19 21 4 44 62 152 16 20 202 20 21 22 5 81 67 162 23 26 221 22 25 23 5 67 51 167 24 22 232 23 22 24 5 51 71 35 25 23 35 24 23 25 5 71 90 39 26 24 41 26 24 26 5 90 81 81 22 25 90 25 26 27 6 14 13 56 28 29 62 28 29 28 6 13 44 58 20 27 64 27 28 29 6 11 14 51 27 19 56 29 30 30 7 51 79 166 31 35 229 31 36 31 7 79 42 147 32 30 194 32 31 32 7 42 43 148 33 31 195 33 32 33 7 43 29 115 34 32 146 34 33 34 7 29 71 40 35 33 42 37 34 35 7 71 51 24 30 34 35 36 37 36 8 85 74 193 37 38 287 38 39 37 8 74 86 101 38 36 124 40 38 38 8 86 85 195 36 37 289 39 40 39 9 90 71 25 40 43 41 43 46 40 9 71 29 34 41 39 42 44 43 41 9 29 69 116 42 40 147 45 44 42 9 69 82 181 43 41 257 47 45 43 9 82 90 82 39 42 92 46 47 44 10 34 32 88 45 49 102 48 52 45 10 32 33 87 46 44 101 49 48 46 10 33 19 77 47 45 85 50 49 47 10 19 88 86 48 46 99 51 50 48 10 88 29 116 49 47 148 53 51 49 10 29 34 115 44 48 145 52 53 50 11 96 65 97 2 18 116 54 55 51 11 50 14 59 29 1 66 56 57 52 12 52 26 109 53 55 137 58 60 53 12 26 21 93 54 52 111 59 58 54 12 21 64 96 55 53 114 61 59 55 12 64 52 168 52 54 234 60 61 56 13 14 84 60 57 27 67 63 62 57 13 84 53 170 58 56 237 65 63 58 13 53 44 153 28 57 205 64 65 59 14 50 94 165 60 51 226 68 66 60 14 94 84 190 56 59 278 67 68 61 15 68 41 146 62 67 192 69 74 62 15 41 39 85 63 61 97 70 69 63 15 39 36 132 64 62 170 71 70 64 15 36 16 68 65 63 76 72 71 65 15 16 35 70 66 64 78 73 72 66 15 35 46 131 67 65 166 75 73 67 15 46 68 157 61 66 213 74 75 68 16 36 24 103 69 64 127 77 76 69 16 24 57 102 70 68 125 79 77 70 16 57 35 129 65 69 164 78 79 71 17 33 47 92 72 77 107 80 86 72 17 47 25 105 73 71 131 81 80 73 17 25 97 108 74 72 134 82 81 74 17 97 63 179 75 73 253 83 82 75 17 63 56 158 76 74 215 84 83 76 17 56 19 83 77 75 95 87 84 77 17 19 33 46 71 76 85 86 87 78 18 82 31 122 79 82 154 88 93 79 18 31 54 125 80 78 158 89 88 80 18 54 81 172 81 79 241 91 89 81 18 81 90 26 82 80 90 94 91 82 18 90 82 43 78 81 92 93 94 83 19 56 28 112 84 76 141 96 95 84 19 28 39 114 85 83 143 98 96 85 19 39 41 62 86 84 97 100 98 86 19 41 88 146 47 85 191 99 100 87 20 33 32 45 88 92 101 103 108 88 20 32 34 44 89 87 102 104 103 89 20 34 89 128 90 88 162 105 104 90 20 89 66 185 91 89 268 106 105 91 20 66 47 105 92 90 130 109 106 92 20 47 33 71 87 91 107 108 109 93 21 26 22 4 94 53 110 112 111 94 21 22 92 98 95 93 119 113 112 95 21 92 38 138 96 94 180 115 113 96 21 38 64 137 54 95 176 114 115 97 22 65 96 50 98 3 116 120 117 98 22 96 92 17 94 97 118 119 120 99 23 75 87 192 100 10 282 123 121 100 23 87 60 164 11 99 224 122 123 101 24 86 74 37 102 104 124 126 128 102 24 74 57 174 69 101 245 125 126 103 24 36 100 132 104 68 169 129 127 104 24 100 86 178 101 103 250 128 129 105 25 47 66 91 106 72 130 132 131 106 25 66 44 151 107 105 201 133 132 107 25 44 78 153 108 106 204 135 133 108 25 78 97 196 73 107 293 134 135 109 26 52 67 169 5 52 235 136 137 110 27 62 43 150 111 12 200 140 138 111 27 43 45 149 13 110 197 139 140 112 28 56 49 158 113 83 216 142 141 113 28 49 100 161 114 112 219 144 142 114 28 100 39 132 84 113 168 143 144 115 29 43 34 126 49 33 160 145 146 116 29 88 69 187 41 48 272 147 148 117 30 64 91 137 118 121 175 149 152 118 30 91 93 139 119 117 181 150 149 119 30 93 77 145 120 118 188 151 150 120 30 77 95 144 121 119 186 153 151 121 30 95 64 168 117 120 233 152 153 122 31 82 61 181 123 78 258 155 154 123 31 61 73 183 124 122 261 156 155 124 31 73 55 156 125 123 211 159 156 125 31 55 54 8 79 124 157 158 159 126 34 43 80 150 127 115 199 161 160 127 34 80 98 184 128 126 263 163 161 128 34 98 89 185 89 127 266 162 163 129 35 57 83 175 130 70 246 165 164 130 35 83 55 173 131 129 244 167 165 131 35 55 46 156 66 130 212 166 167 132 36 39 100 114 103 63 168 169 170 133 37 99 48 154 134 135 206 171 172 134 37 48 79 147 135 133 193 173 171 135 37 79 99 166 133 134 228 172 173 136 38 45 91 155 137 14 209 177 174 137 38 91 64 117 96 136 175 176 177 138 38 92 62 16 15 95 178 179 180 139 40 93 91 118 140 145 181 182 189 140 40 91 99 155 141 139 208 183 182 141 40 99 51 166 142 140 230 184 183 142 40 51 70 167 143 141 231 185 184 143 40 70 95 186 144 142 269 187 185 144 40 95 77 120 145 143 186 190 187 145 40 77 93 119 139 144 188 189 190 146 41 68 88 187 86 61 273 191 192 147 42 79 48 134 148 31 193 196 194 148 42 48 43 149 32 147 198 195 196 149 43 48 45 154 111 148 207 197 198 150 43 62 80 184 126 110 264 199 200 151 44 66 98 185 152 106 267 203 201 152 44 98 62 184 21 151 265 202 203 153 44 53 78 171 107 58 238 204 205 154 45 48 99 133 155 149 206 210 207 155 45 99 91 140 136 154 208 209 210 156 46 55 73 124 157 131 211 214 212 157 46 73 68 183 67 156 260 213 214 158 49 56 63 75 159 112 215 217 216 159 49 63 59 179 160 158 254 218 217 160 49 59 58 176 161 159 248 220 218 161 49 58 100 178 113 160 251 219 220 162 50 67 81 22 163 6 221 223 222 163 50 81 60 172 164 162 240 225 223 164 50 60 87 100 165 163 224 227 225 165 50 87 94 191 59 164 280 226 227 166 51 99 79 135 30 141 228 229 230 167 51 67 70 186 142 23 270 231 232 168 52 64 95 121 169 55 233 236 234 169 52 95 67 186 109 168 271 235 236 170 53 84 85 193 171 57 285 239 237 171 53 85 78 196 153 170 294 238 239 172 54 60 81 163 80 7 240 241 242 173 55 83 75 188 9 130 275 243 244 174 57 74 72 189 175 102 277 247 245 175 57 72 83 188 129 174 276 246 247 176 58 59 76 180 177 160 255 249 248 177 58 76 86 195 178 176 290 252 249 178 58 86 100 104 161 177 250 251 252 179 59 63 97 74 180 159 253 256 254 180 59 97 76 194 176 179 288 255 256 181 61 82 69 42 182 122 257 259 258 182 61 69 68 187 183 181 274 262 259 183 61 68 73 157 123 182 260 261 262 184 62 98 80 127 150 152 263 264 265 185 66 89 98 128 151 90 266 267 268 186 67 95 70 143 167 169 269 270 271 187 68 69 88 116 146 182 272 273 274 188 72 75 83 173 175 192 275 276 283 189 72 74 84 193 190 174 286 279 277 190 72 84 94 60 191 189 278 281 279 191 72 94 87 165 192 190 280 284 281 192 72 87 75 99 188 191 282 283 284 193 74 85 84 170 189 36 285 286 287 194 76 97 85 196 195 180 292 291 288 195 76 85 86 38 177 194 289 290 291 196 78 85 97 194 108 171 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 -26.600690 55.531286 0.411511 3 -28.383795 54.027186 0.390634 4 -35.137862 51.228764 0.307410 5 -45.922325 50.663234 0.195921 6 -61.701213 95.921512 0.365986 1 -48.244698 106.690352 0.457425 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. 09 May 2025 9:46:09.491 PM