August 31 2010 11:31:39.863 AM 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 ********** 89.560726 2 -91.129774 12.028447 3 -20.381664 83.864768 4 ********** 44.626983 5 24.119403 -23.869638 I X Y 1 -2.680773 1.563130 2 -1.590515 0.209936 3 -0.355727 1.463716 4 -1.832381 0.778888 5 0.420963 -0.416604 I X Y Z 1 -0.006867 -0.003409 0.999971 2 -0.019284 -0.977854 0.208397 3 0.100184 -0.037222 0.994272 4 -0.184052 -0.687484 0.702488 5 0.834632 0.373688 -0.404657 STRIPACK triangulation data structure, n = 100 Node Longitude Latitude neighbors of node 1 -0.153597E+03 0.895607E+02 38 37 93 62 2 -0.911298E+02 0.120284E+02 77 78 88 61 3 -0.203817E+02 0.838648E+02 63 56 65 43 4 -0.104988E+03 0.446270E+02 55 92 18 77 81 33 5 0.241194E+02 -0.238696E+02 83 59 52 23 76 10 15 58 6 -0.171248E+03 -0.348536E+01 74 91 45 72 9 11 7 -0.286954E+02 -0.767242E+02 29 35 49 85 8 89 8 -0.367729E+02 -0.890361E+02 85 66 97 17 89 7 9 0.171571E+03 -0.275254E+02 72 22 14 20 99 11 6 10 0.693378E+02 -0.283961E+02 76 28 15 5 11 -0.178220E+03 -0.507687E+02 99 73 66 12 80 74 6 9 12 -0.133228E+03 -0.660311E+02 100 98 53 39 80 11 66 13 -0.163160E+03 0.720944E+02 48 79 92 55 33 93 14 0.122320E+03 -0.203821E+02 22 31 28 20 9 15 0.642560E+02 -0.981319E+01 28 25 58 5 10 16 0.295023E+02 0.291478E+02 57 90 51 70 83 58 17 0.840693E+02 -0.871005E+02 97 73 76 94 89 8 18 -0.138225E+03 0.271539E+02 92 46 91 74 78 77 4 19 0.122508E+03 0.263536E+02 68 82 36 40 75 50 20 0.120599E+03 -0.318623E+02 28 84 99 9 14 21 0.887531E+02 0.640246E+02 69 44 71 32 67 68 48 22 0.123552E+03 -0.177684E+02 9 72 40 36 31 14 23 0.103820E+02 -0.527626E+02 35 29 94 76 5 52 24 0.595745E+02 0.843371E+02 43 65 69 25 0.830652E+02 0.177119E+02 28 31 41 32 57 58 15 26 -0.321825E+02 0.311365E+02 70 30 81 54 59 27 -0.519402E+02 -0.776118E+01 59 54 87 60 47 95 28 0.847357E+02 -0.305973E+02 76 84 20 14 31 25 15 10 29 -0.103256E+02 -0.719311E+02 89 94 23 35 7 30 -0.134548E+02 0.459816E+02 70 51 90 65 56 33 81 26 31 0.934890E+02 0.190248E+02 28 14 22 36 82 41 25 32 0.728854E+02 0.394286E+02 41 67 21 71 57 25 33 -0.871323E+02 0.715202E+02 93 13 55 4 81 30 56 34 0.157574E+03 0.284812E+02 96 72 45 50 35 -0.158027E+02 -0.628710E+02 23 52 95 49 7 29 36 0.111056E+03 0.202167E+02 82 31 22 40 19 37 0.147185E+03 0.861588E+02 69 48 93 1 38 43 38 0.701556E+02 0.898456E+02 62 63 43 37 1 39 -0.100984E+03 -0.437764E+02 53 49 88 80 12 40 0.127784E+03 0.915577E+01 36 22 72 75 19 41 0.879828E+02 0.286286E+02 82 67 32 25 31 42 -0.715993E+02 0.972092E+01 61 60 87 81 43 0.619087E+02 0.859968E+02 69 37 38 63 3 65 24 44 0.427371E+02 0.723462E+02 65 90 71 21 69 45 0.168327E+03 0.284264E+02 72 6 91 46 64 50 34 46 0.176486E+03 0.411945E+02 91 18 92 79 64 45 47 -0.583177E+02 -0.175581E+02 60 88 49 95 27 48 0.151476E+03 0.771530E+02 93 37 69 21 68 64 79 13 49 -0.579336E+02 -0.633897E+02 35 95 47 88 39 53 85 7 50 0.155140E+03 0.314153E+02 45 64 68 19 75 96 34 51 0.167503E+02 0.485306E+02 70 16 90 30 52 -0.124155E+02 -0.289219E+02 5 59 95 35 23 53 -0.116311E+03 -0.691526E+02 98 86 85 49 39 12 54 -0.439221E+02 0.205865E+02 81 87 27 59 26 55 -0.116921E+03 0.577111E+02 33 13 92 4 56 -0.515310E+02 0.804770E+02 63 93 33 30 65 3 57 0.596435E+02 0.416032E+02 71 90 16 58 25 32 58 0.426777E+02 -0.431297E+00 25 57 16 83 5 15 59 -0.111790E+02 -0.225357E+02 5 83 70 26 54 27 95 52 60 -0.661493E+02 -0.141289E+02 27 87 42 61 88 47 61 -0.832032E+02 0.952255E+01 88 60 42 81 77 2 62 -0.108028E+03 0.896255E+02 93 63 38 1 63 -0.304459E+02 0.882711E+02 93 56 3 43 38 62 64 0.172549E+03 0.443357E+02 79 48 68 50 45 46 65 0.155562E+02 0.816766E+02 56 30 90 44 69 24 43 3 66 -0.133373E+03 -0.732106E+02 85 86 100 12 11 73 97 8 67 0.109707E+03 0.421243E+02 68 21 32 41 82 68 0.124496E+03 0.453171E+02 82 19 50 64 48 21 67 69 0.752518E+02 0.804327E+02 65 44 21 48 37 43 24 70 0.104304E+02 0.371117E+02 26 59 83 16 51 30 71 0.589827E+02 0.564857E+02 32 21 44 90 57 72 0.144266E+03 0.105470E+02 22 9 6 45 34 96 75 40 73 0.158819E+03 -0.788930E+02 97 66 11 99 84 76 17 74 -0.138809E+03 -0.353146E+01 91 6 11 80 78 18 75 0.132000E+03 0.175842E+02 72 96 50 19 40 76 0.476525E+02 -0.652443E+02 94 17 73 84 28 10 5 23 77 -0.925964E+02 0.157311E+02 61 81 4 18 78 2 78 -0.122258E+03 0.359429E+01 80 88 2 77 18 74 79 -0.172355E+03 0.694581E+02 48 64 46 92 13 80 -0.111509E+03 -0.353142E+02 88 78 74 11 12 39 81 -0.675576E+02 0.305383E+02 87 54 26 30 33 4 77 61 42 82 0.111473E+03 0.296892E+02 19 68 67 41 31 36 83 0.265387E+02 0.661842E+00 58 16 70 59 5 84 0.117454E+03 -0.429165E+02 28 76 73 99 20 85 -0.768762E+02 -0.762193E+02 49 53 86 66 8 7 86 -0.122215E+03 -0.717751E+02 98 100 66 85 53 87 -0.608673E+02 0.886783E+01 60 27 54 81 42 88 -0.842928E+02 -0.223947E+02 61 2 78 80 39 49 47 60 89 0.161212E+02 -0.872731E+02 17 94 29 7 8 90 0.192673E+02 0.527248E+02 16 57 71 44 65 30 51 91 -0.172219E+03 0.217580E+02 6 74 18 46 45 92 -0.126451E+03 0.492486E+02 13 79 46 18 4 55 93 -0.142235E+03 0.816445E+02 37 48 13 33 56 63 62 1 94 0.411170E+02 -0.694360E+02 23 29 89 17 76 95 -0.337352E+02 -0.326767E+02 59 27 47 49 35 52 96 0.144550E+03 0.174276E+02 34 50 75 72 97 0.159422E+03 -0.813325E+02 17 8 66 73 98 -0.123617E+03 -0.694430E+02 53 12 100 86 99 0.157305E+03 -0.511263E+02 9 20 84 73 11 100 -0.131230E+03 -0.718968E+02 86 98 12 66 NB = 0 boundary arcs. NA = 294 arcs. NT = 196 triangles. STRIPACK (TRLIST) output, n = 100 Node Longitude Latitude 1 -0.153597E+03 0.895607E+02 2 -0.911298E+02 0.120284E+02 3 -0.203817E+02 0.838648E+02 4 -0.104988E+03 0.446270E+02 5 0.241194E+02 -0.238696E+02 6 -0.171248E+03 -0.348536E+01 7 -0.286954E+02 -0.767242E+02 8 -0.367729E+02 -0.890361E+02 9 0.171571E+03 -0.275254E+02 10 0.693378E+02 -0.283961E+02 11 -0.178220E+03 -0.507687E+02 12 -0.133228E+03 -0.660311E+02 13 -0.163160E+03 0.720944E+02 14 0.122320E+03 -0.203821E+02 15 0.642560E+02 -0.981319E+01 16 0.295023E+02 0.291478E+02 17 0.840693E+02 -0.871005E+02 18 -0.138225E+03 0.271539E+02 19 0.122508E+03 0.263536E+02 20 0.120599E+03 -0.318623E+02 21 0.887531E+02 0.640246E+02 22 0.123552E+03 -0.177684E+02 23 0.103820E+02 -0.527626E+02 24 0.595745E+02 0.843371E+02 25 0.830652E+02 0.177119E+02 26 -0.321825E+02 0.311365E+02 27 -0.519402E+02 -0.776118E+01 28 0.847357E+02 -0.305973E+02 29 -0.103256E+02 -0.719311E+02 30 -0.134548E+02 0.459816E+02 31 0.934890E+02 0.190248E+02 32 0.728854E+02 0.394286E+02 33 -0.871323E+02 0.715202E+02 34 0.157574E+03 0.284812E+02 35 -0.158027E+02 -0.628710E+02 36 0.111056E+03 0.202167E+02 37 0.147185E+03 0.861588E+02 38 0.701556E+02 0.898456E+02 39 -0.100984E+03 -0.437764E+02 40 0.127784E+03 0.915577E+01 41 0.879828E+02 0.286286E+02 42 -0.715993E+02 0.972092E+01 43 0.619087E+02 0.859968E+02 44 0.427371E+02 0.723462E+02 45 0.168327E+03 0.284264E+02 46 0.176486E+03 0.411945E+02 47 -0.583177E+02 -0.175581E+02 48 0.151476E+03 0.771530E+02 49 -0.579336E+02 -0.633897E+02 50 0.155140E+03 0.314153E+02 51 0.167503E+02 0.485306E+02 52 -0.124155E+02 -0.289219E+02 53 -0.116311E+03 -0.691526E+02 54 -0.439221E+02 0.205865E+02 55 -0.116921E+03 0.577111E+02 56 -0.515310E+02 0.804770E+02 57 0.596435E+02 0.416032E+02 58 0.426777E+02 -0.431297E+00 59 -0.111790E+02 -0.225357E+02 60 -0.661493E+02 -0.141289E+02 61 -0.832032E+02 0.952255E+01 62 -0.108028E+03 0.896255E+02 63 -0.304459E+02 0.882711E+02 64 0.172549E+03 0.443357E+02 65 0.155562E+02 0.816766E+02 66 -0.133373E+03 -0.732106E+02 67 0.109707E+03 0.421243E+02 68 0.124496E+03 0.453171E+02 69 0.752518E+02 0.804327E+02 70 0.104304E+02 0.371117E+02 71 0.589827E+02 0.564857E+02 72 0.144266E+03 0.105470E+02 73 0.158819E+03 -0.788930E+02 74 -0.138809E+03 -0.353146E+01 75 0.132000E+03 0.175842E+02 76 0.476525E+02 -0.652443E+02 77 -0.925964E+02 0.157311E+02 78 -0.122258E+03 0.359429E+01 79 -0.172355E+03 0.694581E+02 80 -0.111509E+03 -0.353142E+02 81 -0.675576E+02 0.305383E+02 82 0.111473E+03 0.296892E+02 83 0.265387E+02 0.661842E+00 84 0.117454E+03 -0.429165E+02 85 -0.768762E+02 -0.762193E+02 86 -0.122215E+03 -0.717751E+02 87 -0.608673E+02 0.886783E+01 88 -0.842928E+02 -0.223947E+02 89 0.161212E+02 -0.872731E+02 90 0.192673E+02 0.527248E+02 91 -0.172219E+03 0.217580E+02 92 -0.126451E+03 0.492486E+02 93 -0.142235E+03 0.816445E+02 94 0.411170E+02 -0.694360E+02 95 -0.337352E+02 -0.326767E+02 96 0.144550E+03 0.174276E+02 97 0.159422E+03 -0.813325E+02 98 -0.123617E+03 -0.694430E+02 99 0.157305E+03 -0.511263E+02 100 -0.131230E+03 -0.718968E+02 triangle vertices neighbors arcs kt n1 n2 n3 kt1 kt2 kt3 ka1 ka2 ka3 1 1 38 37 148 2 4 190 1 3 2 1 37 93 147 3 1 188 2 1 3 1 93 62 186 4 2 268 4 2 4 1 62 38 150 1 3 194 3 4 5 2 77 78 85 6 8 97 5 7 6 2 78 88 195 7 5 290 6 5 7 2 88 61 184 8 6 260 8 6 8 2 61 77 185 5 7 264 7 8 9 3 63 56 181 10 12 253 9 11 10 3 56 65 132 11 9 162 10 9 11 3 65 43 108 12 10 128 12 10 12 3 43 63 151 9 11 195 11 12 13 4 55 92 62 14 18 66 13 17 14 4 92 18 81 15 13 92 14 13 15 4 18 77 85 16 14 98 15 14 16 4 77 81 185 17 15 263 16 15 17 4 81 33 134 18 16 165 18 16 18 4 33 55 63 13 17 68 17 18 19 5 83 59 183 20 26 259 19 25 20 5 59 52 177 21 19 245 20 19 21 5 52 23 107 22 20 127 21 20 22 5 23 76 106 23 21 124 22 21 23 5 76 10 48 24 22 50 23 22 24 5 10 15 49 25 23 51 24 23 25 5 15 58 70 26 24 78 26 24 26 5 58 83 75 19 25 83 25 26 27 6 74 91 83 28 32 94 27 31 28 6 91 45 164 29 27 218 28 27 29 6 45 72 141 30 28 178 29 28 30 6 72 9 43 31 29 44 30 29 31 6 9 11 47 32 30 48 32 30 32 6 11 74 54 27 31 57 31 32 33 7 29 35 104 34 38 121 33 37 34 7 35 49 145 35 33 184 34 33 35 7 49 85 174 36 34 236 35 34 36 7 85 8 39 37 35 39 36 35 37 7 8 89 42 38 36 42 38 36 38 7 89 29 128 33 37 156 37 38 39 8 85 66 187 40 36 270 40 39 40 8 66 97 189 41 39 274 41 40 41 8 97 17 77 42 40 87 43 41 42 8 17 89 80 37 41 90 42 43 43 9 72 22 101 44 30 116 45 44 44 9 22 14 66 45 43 73 46 45 45 9 14 20 68 46 44 75 47 46 46 9 20 99 93 47 45 107 49 47 47 9 99 11 50 31 46 53 48 49 48 10 76 28 127 49 23 153 52 50 49 10 28 15 69 24 48 77 51 52 50 11 99 73 192 51 47 283 54 53 51 11 73 66 189 52 50 275 55 54 52 11 66 12 59 53 51 64 56 55 53 11 12 80 58 54 52 61 58 56 54 11 80 74 194 32 53 288 57 58 55 12 100 98 196 56 59 292 59 63 56 12 98 53 178 57 55 246 60 59 57 12 53 39 152 58 56 198 62 60 58 12 39 80 154 53 57 200 61 62 59 12 66 100 188 55 52 271 63 64 60 13 48 79 173 61 65 234 65 71 61 13 79 92 167 62 60 222 67 65 62 13 92 55 13 63 61 66 69 67 63 13 55 33 18 64 62 68 70 69 64 13 33 93 139 65 63 176 72 70 65 13 93 48 147 60 64 187 71 72 66 14 22 31 103 67 44 119 74 73 67 14 31 28 111 68 66 132 76 74 68 14 28 20 92 45 67 106 75 76 69 15 28 25 111 70 49 133 79 77 70 15 25 58 115 25 69 138 78 79 71 16 57 90 182 72 76 255 80 85 72 16 90 51 130 73 71 159 81 80 73 16 51 70 129 74 72 157 82 81 74 16 70 83 183 75 73 257 84 82 75 16 83 58 26 76 74 83 86 84 76 16 58 57 115 71 75 137 85 86 77 17 97 73 189 78 41 273 88 87 78 17 73 76 193 79 77 285 89 88 79 17 76 94 106 80 78 123 91 89 80 17 94 89 128 42 79 154 90 91 81 18 92 46 167 82 14 223 93 92 82 18 46 91 164 83 81 217 95 93 83 18 91 74 27 84 82 94 96 95 84 18 74 78 194 85 83 287 99 96 85 18 78 77 5 15 84 97 98 99 86 19 68 82 190 87 91 276 100 104 87 19 82 36 135 88 86 168 101 100 88 19 36 40 102 89 87 117 102 101 89 19 40 75 155 90 88 202 103 102 90 19 75 50 176 91 89 242 105 103 91 19 50 68 175 86 90 239 104 105 92 20 28 84 127 93 68 152 108 106 93 20 84 99 192 46 92 282 107 108 94 21 69 44 163 95 100 216 109 114 95 21 44 71 162 96 94 212 110 109 96 21 71 32 138 97 95 175 111 110 97 21 32 67 137 98 96 172 112 111 98 21 67 68 190 99 97 277 113 112 99 21 68 48 172 100 98 232 115 113 100 21 48 69 146 94 99 186 114 115 101 22 72 40 155 102 43 203 118 116 102 22 40 36 88 103 101 117 120 118 103 22 36 31 135 66 102 169 119 120 104 23 35 29 33 105 107 121 122 126 105 23 29 94 128 106 104 155 125 122 106 23 94 76 79 22 105 123 124 125 107 23 52 35 144 104 21 183 126 127 108 24 43 65 11 109 110 128 129 130 109 24 65 69 163 110 108 214 131 129 110 24 69 43 149 108 109 191 130 131 111 25 28 31 67 112 69 132 134 133 112 25 31 41 136 113 111 170 135 134 113 25 41 32 137 114 112 173 136 135 114 25 32 57 138 115 113 174 139 136 115 25 57 58 76 70 114 137 138 139 116 26 70 30 129 117 120 158 140 143 117 26 30 81 134 118 116 166 141 140 118 26 81 54 180 119 117 251 142 141 119 26 54 59 121 120 118 145 144 142 120 26 59 70 183 116 119 258 143 144 121 27 59 54 119 122 126 145 146 150 122 27 54 87 180 123 121 250 147 146 123 27 87 60 158 124 122 206 148 147 124 27 60 47 169 125 123 226 149 148 125 27 47 95 171 126 124 230 151 149 126 27 95 59 177 121 125 243 150 151 127 28 76 84 193 92 48 284 152 153 128 29 89 94 80 105 38 154 155 156 129 30 70 51 73 130 116 157 160 158 130 30 51 90 72 131 129 159 161 160 131 30 90 65 161 132 130 211 163 161 132 30 65 56 10 133 131 162 164 163 133 30 56 33 139 134 132 177 167 164 134 30 33 81 17 117 133 165 166 167 135 31 36 82 87 136 103 168 171 169 136 31 82 41 156 112 135 205 170 171 137 32 41 67 156 97 113 204 172 173 138 32 71 57 182 114 96 256 174 175 139 33 56 93 181 64 133 252 176 177 140 34 96 72 191 141 143 281 179 181 141 34 72 45 29 142 140 178 180 179 142 34 45 50 166 143 141 220 182 180 143 34 50 96 176 140 142 241 181 182 144 35 52 95 177 145 107 244 185 183 145 35 95 49 171 34 144 229 184 185 146 37 69 48 100 147 149 186 189 192 147 37 48 93 65 2 146 187 188 189 148 37 38 43 151 149 1 196 193 190 149 37 43 69 110 146 148 191 192 193 150 38 62 63 186 151 4 267 197 194 151 38 63 43 12 148 150 195 196 197 152 39 53 49 174 153 57 237 199 198 153 39 49 88 170 154 152 227 201 199 154 39 88 80 195 58 153 289 200 201 155 40 72 75 191 89 101 280 202 203 156 41 82 67 190 137 136 278 204 205 157 42 61 60 184 158 160 262 207 209 158 42 60 87 123 159 157 206 208 207 159 42 87 81 180 160 158 249 210 208 160 42 81 61 185 157 159 265 209 210 161 44 65 90 131 162 163 211 213 215 162 44 90 71 182 95 161 254 212 213 163 44 69 65 109 161 94 214 215 216 164 45 91 46 82 165 28 217 219 218 165 45 46 64 168 166 164 224 221 219 166 45 64 50 175 142 165 240 220 221 167 46 92 79 61 168 81 222 225 223 168 46 79 64 173 165 167 233 224 225 169 47 60 88 184 170 124 261 228 226 170 47 88 49 153 171 169 227 231 228 171 47 49 95 145 125 170 229 230 231 172 48 68 64 175 173 99 238 235 232 173 48 64 79 168 60 172 233 234 235 174 49 53 85 179 35 152 247 236 237 175 50 64 68 172 91 166 238 239 240 176 50 75 96 191 143 90 279 241 242 177 52 59 95 126 144 20 243 244 245 178 53 98 86 196 179 56 294 248 246 179 53 86 85 187 174 178 269 247 248 180 54 81 87 159 122 118 249 250 251 181 56 63 93 186 139 9 266 252 253 182 57 71 90 162 71 138 254 255 256 183 59 83 70 74 120 19 257 258 259 184 60 61 88 7 169 157 260 261 262 185 61 81 77 16 8 160 263 264 265 186 62 93 63 181 150 3 266 267 268 187 66 85 86 179 188 39 269 272 270 188 66 86 100 196 59 187 293 271 272 189 66 73 97 77 40 51 273 274 275 190 67 82 68 86 98 156 276 277 278 191 72 96 75 176 155 140 279 280 281 192 73 99 84 93 193 50 282 286 283 193 73 84 76 127 78 192 284 285 286 194 74 80 78 195 84 54 291 287 288 195 78 80 88 154 6 194 289 290 291 196 86 98 100 55 188 178 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.969195 131.628039 0.034239 3 85.484374 -155.053471 0.071149 4 85.311456 -121.089274 0.075477 1 89.846988 -152.551489 0.004997 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. August 31 2010 11:31:39.886 AM