July 21 2007 12:08:21.938 PM STRIPACK_PRB FORTRAN90 version Test the routines in the STRIPACK library. 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 21.122342 87.102062 2 139.613946 35.991132 3 143.552850 -40.443870 4 66.324136 28.998345 5 121.540238 55.771544 I X Y 1 0.368654 1.520218 2 2.436723 0.628164 3 2.505470 -0.705879 4 1.157575 0.506117 5 2.121277 0.973397 I X Y Z 1 0.047160 0.018219 0.998721 2 -0.616294 0.524249 0.587660 3 -0.612186 0.452121 -0.648703 4 0.351220 0.801017 0.484784 5 -0.294239 0.479399 0.826801 STRIPACK triangulation data structure, n = 100 Node Longitude Latitude neighbors of node 1 0.211223E+02 0.871021E+02 56 83 68 88 2 0.139614E+03 0.359911E+02 55 12 48 23 66 65 42 3 0.143553E+03 -0.404439E+02 54 75 71 60 45 46 4 0.663241E+02 0.289983E+02 82 21 53 96 77 17 5 0.121540E+03 0.557715E+02 14 6 100 80 17 77 66 6 0.111906E+03 0.738008E+02 100 5 14 94 28 7 -0.324334E+02 -0.351967E+02 58 37 91 34 35 70 98 26 27 20 8 -0.142740E+03 -0.283336E+01 87 59 22 78 11 71 86 9 0.127606E+03 0.873706E+02 81 76 56 99 13 10 0.864860E+02 -0.275172E+02 45 38 33 73 43 11 -0.163744E+03 0.864229E+01 78 12 55 60 71 8 12 -0.165187E+03 0.204824E+02 78 31 64 48 2 55 11 13 0.147159E+03 0.848577E+02 99 28 94 39 81 9 14 0.137440E+03 0.687496E+02 66 50 94 6 5 15 0.449177E+02 -0.667815E+02 73 69 52 98 19 46 43 16 0.112020E+02 0.772993E+02 57 72 62 80 88 17 0.748834E+02 0.474197E+02 77 5 80 82 4 18 -0.743588E+02 0.607811E+01 91 37 47 95 30 19 0.341184E+02 -0.772043E+02 98 32 67 51 75 46 15 20 0.232588E+02 -0.172647E+02 84 49 53 58 7 27 21 0.968858E+01 0.415201E+02 72 74 24 58 53 4 82 62 22 -0.146812E+03 0.244241E+02 78 8 59 31 23 0.173334E+03 0.682516E+02 94 61 66 2 48 39 24 -0.403673E+02 0.340182E+02 74 85 47 29 58 21 25 -0.132832E+03 0.514770E+02 93 95 85 90 48 64 31 26 0.643207E+00 -0.483685E+02 98 52 27 7 27 0.432169E+01 -0.471638E+02 52 84 20 7 26 28 0.125549E+03 0.811014E+02 99 88 80 100 6 94 13 29 -0.436391E+02 0.218079E+02 47 37 58 24 30 -0.800176E+02 0.507818E+01 95 87 86 36 97 34 91 18 31 -0.143136E+03 0.420014E+02 64 12 78 22 59 93 25 32 -0.338979E+02 -0.794178E+02 98 70 79 67 19 33 0.721551E+02 -0.613709E+01 38 96 49 69 73 10 34 -0.645614E+02 -0.355907E+02 91 30 97 35 7 35 -0.746421E+02 -0.523082E+02 97 92 70 7 34 36 -0.124987E+03 -0.370872E+02 86 40 71 75 41 92 97 30 37 -0.494288E+02 0.176399E+02 47 18 91 7 58 29 38 0.867377E+02 -0.699280E+01 10 45 65 77 96 33 39 -0.164204E+03 0.739755E+02 48 90 81 13 94 23 40 -0.133503E+03 -0.282445E+02 86 71 36 41 -0.140171E+03 -0.775778E+02 92 36 75 51 63 44 42 0.139266E+03 0.214907E+02 65 45 60 55 2 43 0.860593E+02 -0.351287E+02 73 15 46 45 10 44 -0.963578E+02 -0.841319E+02 63 67 79 70 92 41 45 0.118057E+03 -0.110889E+02 10 43 46 3 60 42 65 38 46 0.111942E+03 -0.536276E+02 45 43 15 19 75 54 3 47 -0.499320E+02 0.210595E+02 29 24 85 89 95 18 37 48 -0.160732E+03 0.632352E+02 90 39 23 2 12 64 25 49 0.459542E+02 -0.695919E+01 84 69 33 96 53 20 50 0.144852E+03 0.694805E+02 66 61 94 14 51 -0.173889E+03 -0.855112E+02 63 41 75 19 67 52 0.109928E+02 -0.501828E+02 26 98 15 69 84 27 53 0.169756E+02 0.663479E+01 49 96 4 21 58 20 54 0.139697E+03 -0.554016E+02 46 75 3 55 0.143452E+03 0.278084E+02 60 11 12 2 42 56 0.113432E+03 0.877142E+02 76 83 1 88 99 9 57 -0.267235E+02 0.753294E+02 88 68 81 90 74 72 16 58 -0.438373E+00 0.623304E+01 21 24 29 37 7 20 53 59 -0.129822E+03 0.333480E+02 22 8 87 93 31 60 0.155528E+03 -0.189092E+02 71 11 55 42 45 3 61 0.149375E+03 0.695154E+02 66 23 94 50 62 0.298943E+02 0.657977E+02 72 21 82 80 16 63 -0.122286E+03 -0.842433E+02 41 51 67 44 64 -0.147499E+03 0.516615E+02 48 12 31 25 65 0.107345E+03 0.555208E+01 45 42 2 66 77 38 66 0.130414E+03 0.473541E+02 77 65 2 23 61 50 14 5 67 -0.936835E+02 -0.860002E+02 63 51 19 32 79 44 68 -0.320601E+02 0.874180E+02 83 81 57 88 1 69 0.521500E+02 -0.299845E+02 84 52 15 73 33 49 70 -0.782874E+02 -0.661765E+02 92 44 79 32 98 7 35 71 -0.175709E+03 -0.199219E+02 3 75 36 40 86 8 11 60 72 0.234922E+02 0.667855E+02 16 57 74 21 62 73 0.658562E+02 -0.323476E+02 15 43 10 33 69 74 -0.309081E+02 0.602420E+02 90 85 24 21 72 57 75 0.157367E+03 -0.685929E+02 46 19 51 41 36 71 3 54 76 0.120732E+03 0.879047E+02 9 81 83 56 77 0.850219E+02 0.270780E+02 96 38 65 66 5 17 4 78 -0.147635E+03 0.237313E+02 31 12 11 8 22 79 -0.505542E+02 -0.788055E+02 70 44 67 32 80 0.665294E+02 0.699968E+02 82 17 5 100 28 88 16 62 81 -0.546884E+02 0.868637E+02 83 76 9 13 39 90 57 68 82 0.426467E+02 0.570777E+02 21 4 17 80 62 83 0.139195E+02 0.892081E+02 1 56 76 81 68 84 0.269261E+02 -0.204594E+02 27 52 69 49 20 85 -0.675344E+02 0.355200E+02 24 74 90 25 95 89 47 86 -0.136401E+03 -0.125803E+02 71 40 36 30 87 8 87 -0.953961E+02 0.187796E+02 95 93 59 8 86 30 88 0.239845E+02 0.821346E+02 68 57 16 80 28 99 56 1 89 -0.710667E+02 0.313203E+02 95 47 85 90 -0.948251E+02 0.734196E+02 81 39 48 25 85 74 57 91 -0.592951E+02 -0.745010E+01 30 34 7 37 18 92 -0.911198E+02 -0.655060E+02 35 97 36 41 44 70 93 -0.130069E+03 0.408697E+02 87 95 25 31 59 94 0.144377E+03 0.787549E+02 6 14 50 61 23 39 13 28 95 -0.750781E+02 0.283030E+02 85 25 93 87 30 18 47 89 96 0.638432E+02 0.120454E+02 53 49 33 38 77 4 97 -0.794561E+02 -0.483212E+02 34 30 36 92 35 98 0.917342E+01 -0.665706E+02 15 52 26 7 70 32 19 99 0.120844E+03 0.870531E+02 56 88 28 13 9 100 0.111087E+03 0.738062E+02 28 80 5 6 NB = 0 boundary arcs. NA = 294 arcs. NT = 196 triangles. stripack (trlist) output, n = 100 Node Longitude Latitude 1 0.211223E+02 0.871021E+02 2 0.139614E+03 0.359911E+02 3 0.143553E+03 -0.404439E+02 4 0.663241E+02 0.289983E+02 5 0.121540E+03 0.557715E+02 6 0.111906E+03 0.738008E+02 7 -0.324334E+02 -0.351967E+02 8 -0.142740E+03 -0.283336E+01 9 0.127606E+03 0.873706E+02 10 0.864860E+02 -0.275172E+02 11 -0.163744E+03 0.864229E+01 12 -0.165187E+03 0.204824E+02 13 0.147159E+03 0.848577E+02 14 0.137440E+03 0.687496E+02 15 0.449177E+02 -0.667815E+02 16 0.112020E+02 0.772993E+02 17 0.748834E+02 0.474197E+02 18 -0.743588E+02 0.607811E+01 19 0.341184E+02 -0.772043E+02 20 0.232588E+02 -0.172647E+02 21 0.968858E+01 0.415201E+02 22 -0.146812E+03 0.244241E+02 23 0.173334E+03 0.682516E+02 24 -0.403673E+02 0.340182E+02 25 -0.132832E+03 0.514770E+02 26 0.643207E+00 -0.483685E+02 27 0.432169E+01 -0.471638E+02 28 0.125549E+03 0.811014E+02 29 -0.436391E+02 0.218079E+02 30 -0.800176E+02 0.507818E+01 31 -0.143136E+03 0.420014E+02 32 -0.338979E+02 -0.794178E+02 33 0.721551E+02 -0.613709E+01 34 -0.645614E+02 -0.355907E+02 35 -0.746421E+02 -0.523082E+02 36 -0.124987E+03 -0.370872E+02 37 -0.494288E+02 0.176399E+02 38 0.867377E+02 -0.699280E+01 39 -0.164204E+03 0.739755E+02 40 -0.133503E+03 -0.282445E+02 41 -0.140171E+03 -0.775778E+02 42 0.139266E+03 0.214907E+02 43 0.860593E+02 -0.351287E+02 44 -0.963578E+02 -0.841319E+02 45 0.118057E+03 -0.110889E+02 46 0.111942E+03 -0.536276E+02 47 -0.499320E+02 0.210595E+02 48 -0.160732E+03 0.632352E+02 49 0.459542E+02 -0.695919E+01 50 0.144852E+03 0.694805E+02 51 -0.173889E+03 -0.855112E+02 52 0.109928E+02 -0.501828E+02 53 0.169756E+02 0.663479E+01 54 0.139697E+03 -0.554016E+02 55 0.143452E+03 0.278084E+02 56 0.113432E+03 0.877142E+02 57 -0.267235E+02 0.753294E+02 58 -0.438373E+00 0.623304E+01 59 -0.129822E+03 0.333480E+02 60 0.155528E+03 -0.189092E+02 61 0.149375E+03 0.695154E+02 62 0.298943E+02 0.657977E+02 63 -0.122286E+03 -0.842433E+02 64 -0.147499E+03 0.516615E+02 65 0.107345E+03 0.555208E+01 66 0.130414E+03 0.473541E+02 67 -0.936835E+02 -0.860002E+02 68 -0.320601E+02 0.874180E+02 69 0.521500E+02 -0.299845E+02 70 -0.782874E+02 -0.661765E+02 71 -0.175709E+03 -0.199219E+02 72 0.234922E+02 0.667855E+02 73 0.658562E+02 -0.323476E+02 74 -0.309081E+02 0.602420E+02 75 0.157367E+03 -0.685929E+02 76 0.120732E+03 0.879047E+02 77 0.850219E+02 0.270780E+02 78 -0.147635E+03 0.237313E+02 79 -0.505542E+02 -0.788055E+02 80 0.665294E+02 0.699968E+02 81 -0.546884E+02 0.868637E+02 82 0.426467E+02 0.570777E+02 83 0.139195E+02 0.892081E+02 84 0.269261E+02 -0.204594E+02 85 -0.675344E+02 0.355200E+02 86 -0.136401E+03 -0.125803E+02 87 -0.953961E+02 0.187796E+02 88 0.239845E+02 0.821346E+02 89 -0.710667E+02 0.313203E+02 90 -0.948251E+02 0.734196E+02 91 -0.592951E+02 -0.745010E+01 92 -0.911198E+02 -0.655060E+02 93 -0.130069E+03 0.408697E+02 94 0.144377E+03 0.787549E+02 95 -0.750781E+02 0.283030E+02 96 0.638432E+02 0.120454E+02 97 -0.794561E+02 -0.483212E+02 98 0.917342E+01 -0.665706E+02 99 0.120844E+03 0.870531E+02 100 0.111087E+03 0.738062E+02 triangle vertices neighbors arcs kt n1 n2 n3 kt1 kt2 kt3 ka1 ka2 ka3 1 1 56 83 182 2 4 259 1 3 2 1 83 68 192 3 1 282 2 1 3 1 68 88 184 4 2 264 4 2 4 1 88 56 183 1 3 263 3 4 5 2 55 12 62 6 11 64 5 10 6 2 12 48 67 7 5 72 6 5 7 2 48 23 111 8 6 130 7 6 8 2 23 66 110 9 7 128 8 7 9 2 66 65 191 10 8 280 9 8 10 2 65 42 165 11 9 220 11 9 11 2 42 55 167 5 10 224 10 11 12 3 54 75 173 13 17 237 12 16 13 3 75 71 151 14 12 194 13 12 14 3 71 60 64 15 13 67 14 13 15 3 60 45 166 16 14 221 15 14 16 3 45 46 168 17 15 226 17 15 17 3 46 54 173 12 16 238 16 17 18 4 82 21 105 19 23 120 18 22 19 4 21 53 104 20 18 118 19 18 20 4 53 96 177 21 19 245 20 19 21 4 96 77 157 22 20 203 21 20 22 4 77 17 28 23 21 27 23 21 23 4 17 82 86 18 22 96 22 23 24 5 14 6 31 25 30 32 24 30 25 5 6 100 33 26 24 34 25 24 26 5 100 80 129 27 25 157 26 25 27 5 80 17 86 28 26 97 28 26 28 5 17 77 22 29 27 27 29 28 29 5 77 66 191 30 28 278 31 29 30 5 66 14 72 24 29 80 30 31 31 6 14 94 73 32 24 81 33 32 32 6 94 28 69 33 31 75 35 33 33 6 28 100 129 25 32 158 34 35 34 7 58 37 131 35 43 162 36 44 35 7 37 91 87 36 34 98 37 36 36 7 91 34 137 37 35 171 38 37 37 7 34 35 146 38 36 187 39 38 38 7 35 70 148 39 37 190 40 39 39 7 70 98 139 40 38 176 41 40 40 7 98 26 124 41 39 149 42 41 41 7 26 27 125 42 40 150 43 42 42 7 27 20 100 43 41 114 45 43 43 7 20 58 99 34 42 111 44 45 44 8 87 59 189 45 50 274 46 51 45 8 59 22 107 46 44 124 47 46 46 8 22 78 108 47 45 126 48 47 47 8 78 11 61 48 46 63 49 48 48 8 11 71 64 49 47 68 50 49 49 8 71 86 160 50 48 210 52 50 50 8 86 87 133 44 49 166 51 52 51 9 81 76 194 52 55 288 53 56 52 9 76 56 182 53 51 260 54 53 53 9 56 99 183 54 52 262 55 54 54 9 99 13 68 55 53 74 57 55 55 9 13 81 71 51 54 78 56 57 56 10 45 38 155 57 60 201 58 61 57 10 38 33 142 58 56 181 59 58 58 10 33 73 145 59 57 185 60 59 59 10 73 43 80 60 58 88 62 60 60 10 43 45 168 56 59 227 61 62 61 11 78 12 65 62 47 70 65 63 62 11 12 55 5 63 61 64 66 65 63 11 55 60 167 64 62 223 69 66 64 11 60 71 14 48 63 67 68 69 65 12 78 31 108 66 61 125 71 70 66 12 31 64 122 67 65 144 73 71 67 12 64 48 121 6 66 142 72 73 68 13 99 28 127 69 54 154 76 74 69 13 28 94 32 70 68 75 77 76 70 13 94 39 112 71 69 131 79 77 71 13 39 81 159 55 70 208 78 79 72 14 66 50 178 73 30 249 82 80 73 14 50 94 179 31 72 251 81 82 74 15 73 69 145 75 80 184 83 89 75 15 69 52 181 76 74 258 84 83 76 15 52 98 124 77 75 148 85 84 77 15 98 19 92 78 76 104 86 85 78 15 19 46 96 79 77 108 87 86 79 15 46 43 168 80 78 228 90 87 80 15 43 73 59 74 79 88 89 90 81 16 57 72 188 82 85 271 91 94 82 16 72 62 106 83 81 121 92 91 83 16 62 80 190 84 82 276 93 92 84 16 80 88 128 85 83 155 95 93 85 16 88 57 184 81 84 265 94 95 86 17 80 82 190 23 27 275 96 97 87 18 91 37 35 88 91 98 99 102 88 18 37 47 130 89 87 160 100 99 89 18 47 95 175 90 88 241 101 100 90 18 95 30 132 91 89 165 103 101 91 18 30 91 137 87 90 172 102 103 92 19 98 32 139 93 77 177 105 104 93 19 32 67 141 94 92 179 106 105 94 19 67 51 180 95 93 255 107 106 95 19 51 75 161 96 94 213 109 107 96 19 75 46 173 78 95 239 108 109 97 20 84 49 176 98 100 244 110 113 98 20 49 53 177 99 97 246 112 110 99 20 53 58 104 43 98 117 111 112 100 20 27 84 126 97 42 152 113 114 101 21 72 74 188 102 106 270 115 122 102 21 74 24 113 103 101 134 116 115 103 21 24 58 116 104 102 137 119 116 104 21 58 53 99 19 103 117 118 119 105 21 82 62 190 106 18 277 123 120 106 21 62 72 82 101 105 121 122 123 107 22 59 31 138 108 45 175 127 124 108 22 31 78 65 46 107 125 126 127 109 23 94 61 179 110 112 250 129 132 110 23 61 66 178 8 109 248 128 129 111 23 48 39 158 112 7 207 133 130 112 23 39 94 70 109 111 131 132 133 113 24 74 85 193 114 102 284 135 134 114 24 85 47 174 115 113 240 136 135 115 24 47 29 130 116 114 161 138 136 116 24 29 58 131 103 115 163 137 138 117 25 93 95 196 118 123 292 139 146 118 25 95 85 195 119 117 291 140 139 119 25 85 90 193 120 118 283 141 140 120 25 90 48 158 121 119 206 143 141 121 25 48 64 67 122 120 142 145 143 122 25 64 31 66 123 121 144 147 145 123 25 31 93 138 117 122 174 146 147 124 26 98 52 76 125 40 148 151 149 125 26 52 27 126 41 124 153 150 151 126 27 52 84 181 100 125 257 152 153 127 28 99 88 183 128 68 261 156 154 128 28 88 80 84 129 127 155 159 156 129 28 80 100 26 33 128 157 158 159 130 29 47 37 88 131 115 160 164 161 131 29 37 58 34 116 130 162 163 164 132 30 95 87 196 133 90 294 167 165 133 30 87 86 50 134 132 166 168 167 134 30 86 36 149 135 133 192 169 168 135 30 36 97 154 136 134 199 170 169 136 30 97 34 146 137 135 188 173 170 137 30 34 91 36 91 136 171 172 173 138 31 59 93 189 123 107 273 174 175 139 32 98 70 39 140 92 176 178 177 140 32 70 79 171 141 139 232 180 178 141 32 79 67 170 93 140 230 179 180 142 33 38 96 157 143 57 204 182 181 143 33 96 49 177 144 142 247 183 182 144 33 49 69 176 145 143 243 186 183 145 33 69 73 74 58 144 184 185 186 146 34 97 35 147 37 136 189 187 188 147 35 97 92 154 148 146 198 191 189 148 35 92 70 172 38 147 234 190 191 149 36 86 40 160 150 134 212 193 192 150 36 40 71 160 151 149 211 195 193 151 36 71 75 13 152 150 194 196 195 152 36 75 41 161 153 151 214 197 196 153 36 41 92 164 154 152 217 200 197 154 36 92 97 147 135 153 198 199 200 155 38 45 65 165 156 56 219 202 201 156 38 65 77 191 157 155 279 205 202 157 38 77 96 21 142 156 203 204 205 158 39 48 90 120 159 111 206 209 207 159 39 90 81 186 71 158 267 208 209 160 40 86 71 49 150 149 210 211 212 161 41 75 51 95 162 152 213 215 214 162 41 51 63 180 163 161 254 216 215 163 41 63 44 169 164 162 229 218 216 164 41 44 92 172 153 163 235 217 218 165 42 65 45 155 166 10 219 222 220 166 42 45 60 15 167 165 221 225 222 167 42 60 55 63 11 166 223 224 225 168 43 46 45 16 60 79 226 227 228 169 44 63 67 180 170 163 253 231 229 170 44 67 79 141 171 169 230 233 231 171 44 79 70 140 172 170 232 236 233 172 44 70 92 148 164 171 234 235 236 173 46 75 54 12 17 96 237 238 239 174 47 85 89 195 175 114 290 242 240 175 47 89 95 195 89 174 289 241 242 176 49 84 69 181 144 97 256 243 244 177 49 96 53 20 98 143 245 246 247 178 50 66 61 110 179 72 248 252 249 179 50 61 94 109 73 178 250 251 252 180 51 67 63 169 162 94 253 254 255 181 52 69 84 176 126 75 256 257 258 182 56 76 83 194 1 52 287 259 260 183 56 88 99 127 53 4 261 262 263 184 57 88 68 3 185 85 264 266 265 185 57 68 81 192 186 184 281 268 266 186 57 81 90 159 187 185 267 269 268 187 57 90 74 193 188 186 285 272 269 188 57 74 72 101 81 187 270 271 272 189 59 87 93 196 138 44 293 273 274 190 62 82 80 86 83 105 275 276 277 191 65 66 77 29 156 9 278 279 280 192 68 83 81 194 185 2 286 281 282 193 74 90 85 119 113 187 283 284 285 194 76 81 83 192 182 51 286 287 288 195 85 95 89 175 174 118 289 290 291 196 87 95 93 117 189 132 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_prb_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_prb_vor.eps". Voronoi region for node 1 Triangle Latitude Longitude Circumradius 2 87.609916 63.587615 0.034421 3 87.956253 -2.161075 0.022714 4 83.476666 -8.966004 0.074517 1 82.431990 70.568015 0.106345 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_PRB Normal end of execution. July 21 2007 12:08:21.967 PM