~meshing/meshing/urop : revision 21

21

//

22

//////////////////////////////////////////////////////////////////////////

23

24

// Arguments: -f filchner-ronne_iceshelf/shorelines.geo -n ../RTopo105b_50S.nc -no -b -85.0:-20.0,-89.0:-75.0 -64.0:-30.0,-89.0:-70.0 -30.0:-20.0,-89.0:-75.0 -p 1 -r 'latitude <= -83'

25

// Source netCDF located at ../RTopo105b_50S.nc

26

// Output to filchner-ronne_iceshelf/shorelines.geo

27

// Boundaries restricted to [1]

28

// Region defined by ((latitude <= -83) and (((longitude >= -85.0) and (longitude <= -20.0) and (latitude >= -89.0) and (latitude <= -75.0)) or ((longitude >= -64.0) and (longitude <= -30.0) and (latitude >= -89.0) and (latitude <= -70.0)) or ((longitude >= -30.0) and (longitude <= -20.0) and (latitude >= -89.0) and (latitude <= -75.0))))

29

//

30

// Including iceshelf ocean cavities

31

// Paths found: 348

32

33

IP = newp;

34

IL = newl;

35

ILL = newll;

36

IS = news;

37

IFI = newf;

38

Point ( IP + 0 ) = { 0, 0, 0 };

39

Point ( IP + 1 ) = { 0, 0, 6.37101e+06 };

40

PolarSphere ( IS + 0 ) = { IP, IP + 1 };

41

42

// Ice-Land mass number 1

43

// Path 1 points 317/59599 area 23814.8 (required closing in 2 parts of the path)

44

Point ( IP + 2 ) = { -8.69030439, 13.84905714, 0.00000000 };

45

Point ( IP + 3 ) = { -8.69867315, 13.86687400, 0.00000000 };

46

Point ( IP + 4 ) = { -8.69464005, 13.86940314, 0.00000000 };

47

Point ( IP + 5 ) = { -8.69060438, 13.87193226, 0.00000000 };

48

Point ( IP + 6 ) = { -8.68656982, 13.87445905, 0.00000000 };

49

Point ( IP + 7 ) = { -8.68253267, 13.87698582, 0.00000000 };

50

Point ( IP + 8 ) = { -8.67849664, 13.87951027, 0.00000000 };

51

Point ( IP + 9 ) = { -8.67445802, 13.88203469, 0.00000000 };

52

Point ( IP + 10 ) = { -8.67041867, 13.88455794, 0.00000000 };

53

Point ( IP + 11 ) = { -8.66638044, 13.88707886, 0.00000000 };

54

Point ( IP + 12 ) = { -8.67472736, 13.90495739, 0.00000000 };

55

Point ( IP + 13 ) = { -8.69552662, 13.93829710, 0.00000000 };

56

Point ( IP + 14 ) = { -8.70393815, 13.95630078, 0.00000000 };

57

Point ( IP + 15 ) = { -8.69987905, 13.95883144, 0.00000000 };

58

Point ( IP + 16 ) = { -8.69581735, 13.96136209, 0.00000000 };

59

Point ( IP + 17 ) = { -8.69175678, 13.96389039, 0.00000000 };

60

Point ( IP + 18 ) = { -8.68769361, 13.96641867, 0.00000000 };

61

Point ( IP + 19 ) = { -8.68363157, 13.96894461, 0.00000000 };

62

Point ( IP + 20 ) = { -8.69203679, 13.98700338, 0.00000000 };

63

Point ( IP + 21 ) = { -8.70045717, 14.00509746, 0.00000000 };

64

Point ( IP + 22 ) = { -8.70889849, 14.02323412, 0.00000000 };

65

Point ( IP + 23 ) = { -8.72993953, 14.05711480, 0.00000000 };

66

Point ( IP + 24 ) = { -8.73844666, 14.07537939, 0.00000000 };

67

Point ( IP + 25 ) = { -8.74697317, 14.09368832, 0.00000000 };

68

Point ( IP + 26 ) = { -8.75552491, 14.11204902, 0.00000000 };

69

Point ( IP + 27 ) = { -8.75142052, 14.11459468, 0.00000000 };

70

Point ( IP + 28 ) = { -8.74731350, 14.11714031, 0.00000000 };

71

Point ( IP + 29 ) = { -8.74320762, 14.11968358, 0.00000000 };

72

Point ( IP + 30 ) = { -8.73909913, 14.12222682, 0.00000000 };

73

Point ( IP + 31 ) = { -8.73498989, 14.12476887, 0.00000000 };

74

Point ( IP + 32 ) = { -8.74353629, 14.14318555, 0.00000000 };

75

Point ( IP + 33 ) = { -8.76480433, 14.17758788, 0.00000000 };

76

Point ( IP + 34 ) = { -8.78618554, 14.21217323, 0.00000000 };

77

Point ( IP + 35 ) = { -8.79484106, 14.23080008, 0.00000000 };

78

Point ( IP + 36 ) = { -8.79070023, 14.23335834, 0.00000000 };

79

Point ( IP + 37 ) = { -8.79937120, 14.25203279, 0.00000000 };

80

Point ( IP + 38 ) = { -8.80806804, 14.27076048, 0.00000000 };

81

Point ( IP + 39 ) = { -8.80391559, 14.27332259, 0.00000000 };

82

Point ( IP + 40 ) = { -8.79976429, 14.27588231, 0.00000000 };

83

Point ( IP + 41 ) = { -8.79561034, 14.27844200, 0.00000000 };

84

Point ( IP + 42 ) = { -8.79145755, 14.28099931, 0.00000000 };

85

Point ( IP + 43 ) = { -8.80015263, 14.29978318, 0.00000000 };

86

Point ( IP + 44 ) = { -8.80886382, 14.31860464, 0.00000000 };

87

Point ( IP + 45 ) = { -8.80469935, 14.32116580, 0.00000000 };

88

Point ( IP + 46 ) = { -8.81342617, 14.34003558, 0.00000000 };

89

Point ( IP + 47 ) = { -8.82217914, 14.35895939, 0.00000000 };

90

Point ( IP + 48 ) = { -8.81800294, 14.36152442, 0.00000000 };

91

Point ( IP + 49 ) = { -8.81382407, 14.36408942, 0.00000000 };

92

Point ( IP + 50 ) = { -8.80964638, 14.36665202, 0.00000000 };

93

Point ( IP + 51 ) = { -8.80546602, 14.36921458, 0.00000000 };

94

Point ( IP + 52 ) = { -8.80128492, 14.37177593, 0.00000000 };

95

Point ( IP + 53 ) = { -8.79710499, 14.37433488, 0.00000000 };

96

Point ( IP + 54 ) = { -8.79292240, 14.37689380, 0.00000000 };

97

Point ( IP + 55 ) = { -8.78874098, 14.37945032, 0.00000000 };

98

Point ( IP + 56 ) = { -8.78455690, 14.38200680, 0.00000000 };

99

Point ( IP + 57 ) = { -8.78037399, 14.38456089, 0.00000000 };

100

Point ( IP + 58 ) = { -8.77618843, 14.38711494, 0.00000000 };

101

Point ( IP + 59 ) = { -8.78490289, 14.40611230, 0.00000000 };

102

Point ( IP + 60 ) = { -8.79363361, 14.42514791, 0.00000000 };

103

Point ( IP + 61 ) = { -8.78943623, 14.42770581, 0.00000000 };

104

Point ( IP + 62 ) = { -8.78524004, 14.43026132, 0.00000000 };

105

Point ( IP + 63 ) = { -8.78104117, 14.43281678, 0.00000000 };

106

Point ( IP + 64 ) = { -8.77684349, 14.43536985, 0.00000000 };

107

Point ( IP + 65 ) = { -8.78557230, 14.45446274, 0.00000000 };

108

Point ( IP + 66 ) = { -8.80732017, 14.49024341, 0.00000000 };

109

Point ( IP + 67 ) = { -8.81612371, 14.50948271, 0.00000000 };

110

Point ( IP + 68 ) = { -8.81190373, 14.51204598, 0.00000000 };

111

Point ( IP + 69 ) = { -8.80768107, 14.51460919, 0.00000000 };

112

Point ( IP + 70 ) = { -8.80345767, 14.51717117, 0.00000000 };

113

Point ( IP + 71 ) = { -8.79923545, 14.51973076, 0.00000000 };

114

Point ( IP + 72 ) = { -8.80803811, 14.53902759, 0.00000000 };

115

Point ( IP + 73 ) = { -8.81685735, 14.55836362, 0.00000000 };

116

Point ( IP + 74 ) = { -8.82569921, 14.57774669, 0.00000000 };

117

Point ( IP + 75 ) = { -8.83456783, 14.59718587, 0.00000000 };

118

Point ( IP + 76 ) = { -8.84345732, 14.61667365, 0.00000000 };

119

Point ( IP + 77 ) = { -8.85236366, 14.63620126, 0.00000000 };

120

Point ( IP + 78 ) = { -8.84810683, 14.63877506, 0.00000000 };

121

Point ( IP + 79 ) = { -8.85702945, 14.65835359, 0.00000000 };

122

Point ( IP + 80 ) = { -8.86597919, 14.67798905, 0.00000000 };

123

Point ( IP + 81 ) = { -8.86171021, 14.68056681, 0.00000000 };

124

Point ( IP + 82 ) = { -8.85743852, 14.68314451, 0.00000000 };

125

Point ( IP + 83 ) = { -8.85316803, 14.68571978, 0.00000000 };

126

Point ( IP + 84 ) = { -8.84889484, 14.68829499, 0.00000000 };

127

Point ( IP + 85 ) = { -8.85784480, 14.70798904, 0.00000000 };

128

Point ( IP + 86 ) = { -8.88008959, 14.74492534, 0.00000000 };

129

Point ( IP + 87 ) = { -8.88911208, 14.76476412, 0.00000000 };

130

Point ( IP + 88 ) = { -8.89815602, 14.78465299, 0.00000000 };

131

Point ( IP + 89 ) = { -8.90721736, 14.80458301, 0.00000000 };

132

Point ( IP + 90 ) = { -8.90291155, 14.80717275, 0.00000000 };

133

Point ( IP + 91 ) = { -8.89860302, 14.80976243, 0.00000000 };

134

Point ( IP + 92 ) = { -8.90767532, 14.82974756, 0.00000000 };

135

Point ( IP + 93 ) = { -8.91677541, 14.84979131, 0.00000000 };

136

Point ( IP + 94 ) = { -8.91245646, 14.85238384, 0.00000000 };

137

Point ( IP + 95 ) = { -8.90813478, 14.85497629, 0.00000000 };

138

Point ( IP + 96 ) = { -8.90381432, 14.85756630, 0.00000000 };

139

Point ( IP + 97 ) = { -8.91291963, 14.87766827, 0.00000000 };

140

Point ( IP + 98 ) = { -8.92204256, 14.89781208, 0.00000000 };

141

Point ( IP + 99 ) = { -8.91770964, 14.90040613, 0.00000000 };

142

Point ( IP + 100 ) = { -8.92684956, 14.92060322, 0.00000000 };

143

Point ( IP + 101 ) = { -8.93601760, 14.94085982, 0.00000000 };

144

Point ( IP + 102 ) = { -8.93167216, 14.94345793, 0.00000000 };

145

Point ( IP + 103 ) = { -8.92732398, 14.94605597, 0.00000000 };

146

Point ( IP + 104 ) = { -8.92297703, 14.94865156, 0.00000000 };

147

Point ( IP + 105 ) = { -8.91862734, 14.95124707, 0.00000000 };

148

Point ( IP + 106 ) = { -8.91427888, 14.95384013, 0.00000000 };

149

Point ( IP + 107 ) = { -8.90992768, 14.95643311, 0.00000000 };

150

Point ( IP + 108 ) = { -8.90557771, 14.95902364, 0.00000000 };

151

Point ( IP + 109 ) = { -8.90122500, 14.96161409, 0.00000000 };

152

Point ( IP + 110 ) = { -8.89687154, 14.96420327, 0.00000000 };

153

Point ( IP + 111 ) = { -8.89251932, 14.96679000, 0.00000000 };

154

Point ( IP + 112 ) = { -8.88816435, 14.96937664, 0.00000000 };

155

Point ( IP + 113 ) = { -8.88381062, 14.97196084, 0.00000000 };

156

Point ( IP + 114 ) = { -8.87945415, 14.97454496, 0.00000000 };

157

Point ( IP + 115 ) = { -8.87509892, 14.97712662, 0.00000000 };

158

Point ( IP + 116 ) = { -8.87074094, 14.97970820, 0.00000000 };

159

Point ( IP + 117 ) = { -8.86638421, 14.98228733, 0.00000000 };

160

Point ( IP + 118 ) = { -8.86202473, 14.98486637, 0.00000000 };

161

Point ( IP + 119 ) = { -8.85766650, 14.98744297, 0.00000000 };

162

Point ( IP + 120 ) = { -8.85330552, 14.99001948, 0.00000000 };

163

Point ( IP + 121 ) = { -8.84894579, 14.99259354, 0.00000000 };

164

Point ( IP + 122 ) = { -8.84458332, 14.99516751, 0.00000000 };

165

Point ( IP + 123 ) = { -8.84022209, 14.99773903, 0.00000000 };

166

Point ( IP + 124 ) = { -8.83585812, 15.00031046, 0.00000000 };

167

Point ( IP + 125 ) = { -8.83149340, 15.00288062, 0.00000000 };

168

Point ( IP + 126 ) = { -8.82712993, 15.00544834, 0.00000000 };

169

Point ( IP + 127 ) = { -8.82276372, 15.00801596, 0.00000000 };

170

Point ( IP + 128 ) = { -8.81839876, 15.01058113, 0.00000000 };

171

Point ( IP + 129 ) = { -8.81403105, 15.01314621, 0.00000000 };

172

Point ( IP + 130 ) = { -8.80966460, 15.01570885, 0.00000000 };

173

Point ( IP + 131 ) = { -8.80529540, 15.01827139, 0.00000000 };

174

Point ( IP + 132 ) = { -8.80092745, 15.02083149, 0.00000000 };

175

Point ( IP + 133 ) = { -8.79655676, 15.02339148, 0.00000000 };

176

Point ( IP + 134 ) = { -8.79218733, 15.02594904, 0.00000000 };

177

Point ( IP + 135 ) = { -8.78781515, 15.02850649, 0.00000000 };

178

Point ( IP + 136 ) = { -8.78344424, 15.03106150, 0.00000000 };

179

Point ( IP + 137 ) = { -8.77907057, 15.03361641, 0.00000000 };

180

Point ( IP + 138 ) = { -8.77469817, 15.03616888, 0.00000000 };

181

Point ( IP + 139 ) = { -8.77032302, 15.03872125, 0.00000000 };

182

Point ( IP + 140 ) = { -8.76594712, 15.04127234, 0.00000000 };

183

Point ( IP + 141 ) = { -8.76157249, 15.04382099, 0.00000000 };

184

Point ( IP + 142 ) = { -8.75719512, 15.04636954, 0.00000000 };

185

Point ( IP + 143 ) = { -8.75281900, 15.04891564, 0.00000000 };

186

Point ( IP + 144 ) = { -8.74844015, 15.05146164, 0.00000000 };

187

Point ( IP + 145 ) = { -8.74406255, 15.05400520, 0.00000000 };

188

Point ( IP + 146 ) = { -8.73968221, 15.05654865, 0.00000000 };

189

Point ( IP + 147 ) = { -8.73530314, 15.05908967, 0.00000000 };

190

Point ( IP + 148 ) = { -8.73092132, 15.06163057, 0.00000000 };

191

Point ( IP + 149 ) = { -8.73983492, 15.08206140, 0.00000000 };

192

Point ( IP + 150 ) = { -8.74876584, 15.10253500, 0.00000000 };

193

Point ( IP + 151 ) = { -8.74437138, 15.10507982, 0.00000000 };

194

Point ( IP + 152 ) = { -8.73997820, 15.10762219, 0.00000000 };

195

Point ( IP + 153 ) = { -8.73558226, 15.11016445, 0.00000000 };

196

Point ( IP + 154 ) = { -8.73118759, 15.11270427, 0.00000000 };

197

Point ( IP + 155 ) = { -8.72679018, 15.11524398, 0.00000000 };

198

Point ( IP + 156 ) = { -8.72239202, 15.11778240, 0.00000000 };

199

Point ( IP + 157 ) = { -8.71799514, 15.12031838, 0.00000000 };

200

Point ( IP + 158 ) = { -8.71359551, 15.12285424, 0.00000000 };

201

Point ( IP + 159 ) = { -8.70919716, 15.12538767, 0.00000000 };

202

Point ( IP + 160 ) = { -8.70479605, 15.12792097, 0.00000000 };

203

Point ( IP + 161 ) = { -8.71368972, 15.14847224, 0.00000000 };

204

Point ( IP + 162 ) = { -8.72261091, 15.16908436, 0.00000000 };

205

Point ( IP + 163 ) = { -8.71819709, 15.17162157, 0.00000000 };

206

Point ( IP + 164 ) = { -8.71378456, 15.17415633, 0.00000000 };

207

Point ( IP + 165 ) = { -8.70936926, 15.17669096, 0.00000000 };

208

Point ( IP + 166 ) = { -8.70495525, 15.17922315, 0.00000000 };

209

Point ( IP + 167 ) = { -8.70053848, 15.18175522, 0.00000000 };

210

Point ( IP + 168 ) = { -8.69612300, 15.18428485, 0.00000000 };

211

Point ( IP + 169 ) = { -8.69170476, 15.18681434, 0.00000000 };

212

Point ( IP + 170 ) = { -8.68728780, 15.18934140, 0.00000000 };

213

Point ( IP + 171 ) = { -8.68286809, 15.19186832, 0.00000000 };

214

Point ( IP + 172 ) = { -8.67844765, 15.19439397, 0.00000000 };

215

Point ( IP + 173 ) = { -8.67402849, 15.19691716, 0.00000000 };

216

Point ( IP + 174 ) = { -8.66960657, 15.19944023, 0.00000000 };

217

Point ( IP + 175 ) = { -8.67847811, 15.22013567, 0.00000000 };

218

Point ( IP + 176 ) = { -8.68736700, 15.24087468, 0.00000000 };

219

Point ( IP + 177 ) = { -8.68293229, 15.24340163, 0.00000000 };

220

Point ( IP + 178 ) = { -8.67849888, 15.24592613, 0.00000000 };

221

Point ( IP + 179 ) = { -8.67406271, 15.24845050, 0.00000000 };

222

Point ( IP + 180 ) = { -8.66962783, 15.25097242, 0.00000000 };

223

Point ( IP + 181 ) = { -8.66519019, 15.25349420, 0.00000000 };

224

Point ( IP + 182 ) = { -8.66075384, 15.25601355, 0.00000000 };

225

Point ( IP + 183 ) = { -8.65631473, 15.25853275, 0.00000000 };

226

Point ( IP + 184 ) = { -8.65187693, 15.26104951, 0.00000000 };

227

Point ( IP + 185 ) = { -8.63860280, 15.24280334, 0.00000000 };

228

Point ( IP + 186 ) = { -8.62535725, 15.22459437, 0.00000000 };

229

Point ( IP + 187 ) = { -8.62092931, 15.22710213, 0.00000000 };

230

Point ( IP + 188 ) = { -8.61649862, 15.22960975, 0.00000000 };

231

Point ( IP + 189 ) = { -8.61206719, 15.23211608, 0.00000000 };

232

Point ( IP + 190 ) = { -8.60763707, 15.23461997, 0.00000000 };

233

Point ( IP + 191 ) = { -8.60320419, 15.23712372, 0.00000000 };

234

Point ( IP + 192 ) = { -8.59877260, 15.23962504, 0.00000000 };

235

Point ( IP + 193 ) = { -8.59433827, 15.24212621, 0.00000000 };

236

Point ( IP + 194 ) = { -8.58990523, 15.24462494, 0.00000000 };

237

Point ( IP + 195 ) = { -8.58546944, 15.24712353, 0.00000000 };

238

Point ( IP + 196 ) = { -8.58103495, 15.24961969, 0.00000000 };

239

Point ( IP + 197 ) = { -8.57659770, 15.25211570, 0.00000000 };

240

Point ( IP + 198 ) = { -8.57216176, 15.25460928, 0.00000000 };

241

Point ( IP + 199 ) = { -8.56772307, 15.25710271, 0.00000000 };

242

Point ( IP + 200 ) = { -8.56328568, 15.25959370, 0.00000000 };

243

Point ( IP + 201 ) = { -8.55014771, 15.24137704, 0.00000000 };

244

Point ( IP + 202 ) = { -8.53704792, 15.22321517, 0.00000000 };

245

Point ( IP + 203 ) = { -8.53262039, 15.22569725, 0.00000000 };

246

Point ( IP + 204 ) = { -8.52819010, 15.22817917, 0.00000000 };

247

Point ( IP + 205 ) = { -8.52375910, 15.23065980, 0.00000000 };

248

Point ( IP + 206 ) = { -8.51932940, 15.23313802, 0.00000000 };

249

Point ( IP + 207 ) = { -8.51489695, 15.23561607, 0.00000000 };

250

Point ( IP + 208 ) = { -8.51046581, 15.23809171, 0.00000000 };

251

Point ( IP + 209 ) = { -8.50603192, 15.24056719, 0.00000000 };

252

Point ( IP + 210 ) = { -8.50159935, 15.24304024, 0.00000000 };

253

Point ( IP + 211 ) = { -8.48857479, 15.22489372, 0.00000000 };

254

Point ( IP + 212 ) = { -8.47557811, 15.20678396, 0.00000000 };

255

Point ( IP + 213 ) = { -8.47115536, 15.20924816, 0.00000000 };

256

Point ( IP + 214 ) = { -8.46672986, 15.21171220, 0.00000000 };

257

Point ( IP + 215 ) = { -8.46230568, 15.21417383, 0.00000000 };

258

Point ( IP + 216 ) = { -8.45787875, 15.21663529, 0.00000000 };

259

Point ( IP + 217 ) = { -8.45345313, 15.21909435, 0.00000000 };

260

Point ( IP + 218 ) = { -8.44902478, 15.22155324, 0.00000000 };

261

Point ( IP + 219 ) = { -8.44459773, 15.22400971, 0.00000000 };

262

Point ( IP + 220 ) = { -8.43167057, 15.20591907, 0.00000000 };

263

Point ( IP + 221 ) = { -8.41878080, 15.18788260, 0.00000000 };

264

Point ( IP + 222 ) = { -8.41436152, 15.19033141, 0.00000000 };

265

Point ( IP + 223 ) = { -8.40994356, 15.19277781, 0.00000000 };

266

Point ( IP + 224 ) = { -8.40552286, 15.19522404, 0.00000000 };

267

Point ( IP + 225 ) = { -8.40110347, 15.19766787, 0.00000000 };

268

Point ( IP + 226 ) = { -8.39668135, 15.20011153, 0.00000000 };

269

Point ( IP + 227 ) = { -8.38385111, 15.18210097, 0.00000000 };

270

Point ( IP + 228 ) = { -8.37104819, 15.16412668, 0.00000000 };

271

Point ( IP + 229 ) = { -8.36663583, 15.16656160, 0.00000000 };

272

Point ( IP + 230 ) = { -8.36222478, 15.16899412, 0.00000000 };

273

Point ( IP + 231 ) = { -8.35781100, 15.17142647, 0.00000000 };

274

Point ( IP + 232 ) = { -8.35339853, 15.17385642, 0.00000000 };

275

Point ( IP + 233 ) = { -8.34898334, 15.17628620, 0.00000000 };

276

Point ( IP + 234 ) = { -8.34456946, 15.17871359, 0.00000000 };

277

Point ( IP + 235 ) = { -8.33182890, 15.16076249, 0.00000000 };

278

Point ( IP + 236 ) = { -8.31912501, 15.14286492, 0.00000000 };

279

Point ( IP + 237 ) = { -8.31472086, 15.14528363, 0.00000000 };

280

Point ( IP + 238 ) = { -8.31031398, 15.14770216, 0.00000000 };

281

Point ( IP + 239 ) = { -8.30590640, 15.15011941, 0.00000000 };

282

Point ( IP + 240 ) = { -8.30150013, 15.15253427, 0.00000000 };

283

Point ( IP + 241 ) = { -8.29709114, 15.15494895, 0.00000000 };

284

Point ( IP + 242 ) = { -8.29268347, 15.15736124, 0.00000000 };

285

Point ( IP + 243 ) = { -8.28004248, 15.13948616, 0.00000000 };

286

Point ( IP + 244 ) = { -8.26742825, 15.12164684, 0.00000000 };

287

Point ( IP + 245 ) = { -8.26303026, 15.12405051, 0.00000000 };

288

Point ( IP + 246 ) = { -8.25862956, 15.12645400, 0.00000000 };

289

Point ( IP + 247 ) = { -8.25423018, 15.12885512, 0.00000000 };

290

Point ( IP + 248 ) = { -8.24982809, 15.13125605, 0.00000000 };

291

Point ( IP + 249 ) = { -8.24542731, 15.13365460, 0.00000000 };

292

Point ( IP + 250 ) = { -8.23287040, 15.11584134, 0.00000000 };

293

Point ( IP + 251 ) = { -8.22034946, 15.09808096, 0.00000000 };

294

Point ( IP + 252 ) = { -8.21595833, 15.10047094, 0.00000000 };

295

Point ( IP + 253 ) = { -8.21156449, 15.10286074, 0.00000000 };

296

Point ( IP + 254 ) = { -8.20717198, 15.10524816, 0.00000000 };

297

Point ( IP + 255 ) = { -8.19469760, 15.08752091, 0.00000000 };

298

Point ( IP + 256 ) = { -8.18224949, 15.06982890, 0.00000000 };

299

Point ( IP + 257 ) = { -8.16983224, 15.05218063, 0.00000000 };

300

Point ( IP + 258 ) = { -8.16182437, 15.03221050, 0.00000000 };

301

Point ( IP + 259 ) = { -8.15383620, 15.01228967, 0.00000000 };

302

Point ( IP + 260 ) = { -8.14586301, 14.99240940, 0.00000000 };

303

Point ( IP + 261 ) = { -8.15022272, 14.99003981, 0.00000000 };

304

Point ( IP + 262 ) = { -8.14226334, 14.97021084, 0.00000000 };

305

Point ( IP + 263 ) = { -8.13432810, 14.95043916, 0.00000000 };

306

Point ( IP + 264 ) = { -8.13867561, 14.94807293, 0.00000000 };

307

Point ( IP + 265 ) = { -8.14302441, 14.94570435, 0.00000000 };

308

Point ( IP + 266 ) = { -8.14737053, 14.94333558, 0.00000000 };

309

Point ( IP + 267 ) = { -8.13943751, 14.92361938, 0.00000000 };

310

Point ( IP + 268 ) = { -8.13151928, 14.90394310, 0.00000000 };

311

Point ( IP + 269 ) = { -8.13585325, 14.90157769, 0.00000000 };

312

Point ( IP + 270 ) = { -8.12794865, 14.88195194, 0.00000000 };

313

Point ( IP + 271 ) = { -8.12006334, 14.86237423, 0.00000000 };

314

Point ( IP + 272 ) = { -8.12438523, 14.86001215, 0.00000000 };

315

Point ( IP + 273 ) = { -8.11651344, 14.84048461, 0.00000000 };

316

Point ( IP + 274 ) = { -8.10866539, 14.82101312, 0.00000000 };

317

Point ( IP + 275 ) = { -8.11297525, 14.81865435, 0.00000000 };

318

Point ( IP + 276 ) = { -8.10514062, 14.79923261, 0.00000000 };

319

Point ( IP + 277 ) = { -8.09732044, 14.77984984, 0.00000000 };

320

Point ( IP + 278 ) = { -8.10161833, 14.77749438, 0.00000000 };

321

Point ( IP + 279 ) = { -8.10591750, 14.77513658, 0.00000000 };

322

Point ( IP + 280 ) = { -8.11021599, 14.77277754, 0.00000000 };

323

Point ( IP + 281 ) = { -8.10239691, 14.75344928, 0.00000000 };

324

Point ( IP + 282 ) = { -8.09460131, 14.73417627, 0.00000000 };

325

Point ( IP + 283 ) = { -8.09888789, 14.73182052, 0.00000000 };

326

Point ( IP + 284 ) = { -8.09110554, 14.71259652, 0.00000000 };

327

Point ( IP + 285 ) = { -8.08333745, 14.69341087, 0.00000000 };

328

Point ( IP + 286 ) = { -8.07558907, 14.67427113, 0.00000000 };

329

Point ( IP + 287 ) = { -8.06786388, 14.65518586, 0.00000000 };

330

Point ( IP + 288 ) = { -8.07212552, 14.65283897, 0.00000000 };

331

Point ( IP + 289 ) = { -8.07638842, 14.65048976, 0.00000000 };

332

Point ( IP + 290 ) = { -8.08064870, 14.64814039, 0.00000000 };

333

Point ( IP + 291 ) = { -8.07292547, 14.62910793, 0.00000000 };

334

Point ( IP + 292 ) = { -8.06521630, 14.61011320, 0.00000000 };

335

Point ( IP + 293 ) = { -8.06946483, 14.60776708, 0.00000000 };

336

Point ( IP + 294 ) = { -8.07371462, 14.60541865, 0.00000000 };

337

Point ( IP + 295 ) = { -8.06601196, 14.58647470, 0.00000000 };

338

Point ( IP + 296 ) = { -8.05833225, 14.56758443, 0.00000000 };

339

Point ( IP + 297 ) = { -8.05067189, 14.54873900, 0.00000000 };

340

Point ( IP + 298 ) = { -8.04302537, 14.52993073, 0.00000000 };

341

Point ( IP + 299 ) = { -8.04725058, 14.52759106, 0.00000000 };

342

Point ( IP + 300 ) = { -8.03961687, 14.50883009, 0.00000000 };

343

Point ( IP + 301 ) = { -8.03200582, 14.49012207, 0.00000000 };

344

Point ( IP + 302 ) = { -8.03622138, 14.48778454, 0.00000000 };

345

Point ( IP + 303 ) = { -8.04043434, 14.48544685, 0.00000000 };

346

Point ( IP + 304 ) = { -8.04464855, 14.48310687, 0.00000000 };

347

Point ( IP + 305 ) = { -8.03703837, 14.46445073, 0.00000000 };

348

Point ( IP + 306 ) = { -8.02944186, 14.44583115, 0.00000000 };

349

Point ( IP + 307 ) = { -8.03364454, 14.44349437, 0.00000000 };

350

Point ( IP + 308 ) = { -8.02606068, 14.42492142, 0.00000000 };

351

Point ( IP + 309 ) = { -8.01849924, 14.40640066, 0.00000000 };

352

Point ( IP + 310 ) = { -8.01095667, 14.38792342, 0.00000000 };

353

Point ( IP + 311 ) = { -8.00342756, 14.36948217, 0.00000000 };

354

Point ( IP + 312 ) = { -7.99591530, 14.35108525, 0.00000000 };

355

Point ( IP + 313 ) = { -7.98425662, 14.33506236, 0.00000000 };

356

Point ( IP + 314 ) = { -7.97262462, 14.31907788, 0.00000000 };

357

Point ( IP + 315 ) = { -7.96101581, 14.30312331, 0.00000000 };

358

Point ( IP + 316 ) = { -7.95685588, 14.30543790, 0.00000000 };

359

Point ( IP + 317 ) = { -7.95269337, 14.30775235, 0.00000000 };

360

// **** END 316/316True

361

LoopStart0 = IP + 2;

362

LoopEnd0 = IP + 317;

363

BSpline ( IL + 0 ) = { IP + 2 : IP + 317 };

364

Physical Line( 3 ) = { IL + 0 };

365

366

// Closing path with parallels and merdians, from (-60.94166565, -83.00000000) to (-57.89166641, -83.00000000)

367

// Drawing meridian to max latitude index 317 at -60.941666.2, -83.000000.2 (to match -83.000000.2)

368

// Drawing parallel index 317 at -60.941666.2 (to match -57.891666.2), -83.000000.2

369

Point ( IP + 318 ) = { -7.96605406, 14.27794652, 0.00000000 };

370

Point ( IP + 319 ) = { -7.99096164, 14.26402139, 0.00000000 };

371

Point ( IP + 320 ) = { -8.01584487, 14.25005281, 0.00000000 };

372

Point ( IP + 321 ) = { -8.04070368, 14.23604083, 0.00000000 };

373

Point ( IP + 322 ) = { -8.06553800, 14.22198548, 0.00000000 };

374

Point ( IP + 323 ) = { -8.09034775, 14.20788680, 0.00000000 };

375

Point ( IP + 324 ) = { -8.11513285, 14.19374485, 0.00000000 };

376

Point ( IP + 325 ) = { -8.13989324, 14.17955966, 0.00000000 };

377

Point ( IP + 326 ) = { -8.16462883, 14.16533128, 0.00000000 };

378

Point ( IP + 327 ) = { -8.18933955, 14.15105974, 0.00000000 };

379

Point ( IP + 328 ) = { -8.21402532, 14.13674510, 0.00000000 };

380

Point ( IP + 329 ) = { -8.23868607, 14.12238740, 0.00000000 };

381

Point ( IP + 330 ) = { -8.26332173, 14.10798668, 0.00000000 };

382

Point ( IP + 331 ) = { -8.28793221, 14.09354298, 0.00000000 };

383

Point ( IP + 332 ) = { -8.31251745, 14.07905635, 0.00000000 };

384

Point ( IP + 333 ) = { -8.33707736, 14.06452684, 0.00000000 };

385

Point ( IP + 334 ) = { -8.36161188, 14.04995448, 0.00000000 };

386

Point ( IP + 335 ) = { -8.38612093, 14.03533932, 0.00000000 };

387

Point ( IP + 336 ) = { -8.41060444, 14.02068141, 0.00000000 };

388

Point ( IP + 337 ) = { -8.43506232, 14.00598079, 0.00000000 };

389

Point ( IP + 338 ) = { -8.45949451, 13.99123750, 0.00000000 };

390

Point ( IP + 339 ) = { -8.48390093, 13.97645159, 0.00000000 };

391

Point ( IP + 340 ) = { -8.50828150, 13.96162311, 0.00000000 };

392

Point ( IP + 341 ) = { -8.53263616, 13.94675210, 0.00000000 };

393

Point ( IP + 342 ) = { -8.55696483, 13.93183861, 0.00000000 };

394

Point ( IP + 343 ) = { -8.58126743, 13.91688268, 0.00000000 };

395

Point ( IP + 344 ) = { -8.60554389, 13.90188435, 0.00000000 };

396

Point ( IP + 345 ) = { -8.62979413, 13.88684368, 0.00000000 };

397

Point ( IP + 346 ) = { -8.65401809, 13.87176070, 0.00000000 };

398

LoopStart1 = IP + 317;

399

LoopEnd1 = IP + 346;

400

BSpline ( IL + 1 ) = { IP + 317 : IP + 346, IP + 2 };

401

Physical Line( 4 ) = { IL + 1 };

402

403

Line Loop( ILL + 0 ) = { IL + 0, IL + 1 };

404

// **** END end of loop False316/316True

405

406

Delete { Point{1}; }

407

Delete { Point{2}; }

408

// Open boundaries (id 4): 2

409

// Closed boundaries (id 3): 1

410

411

Plane Surface( 10 ) = { ILL + 0 };

412

Physical Surface( 10 ) = { 10 };

413

414

Printf("Assigning characteristic mesh sizes...");

415

// Field[ IFI + 1] = Attractor;

416

// Field[ IFI + 1].EdgesList = { 999999, IL + 0 };

417

// Field [ IFI + 1 ].NNodesByEdge = 5e4;

418

//

419

// Field[ IFI + 2] = Threshold;

420

// Field[ IFI + 2].DistMax = 2e6;

421

// Field[ IFI + 2].DistMin = 3e4;

422

// Field[ IFI + 2].IField = IFI + 1;

423

// Field[ IFI + 2].LcMin = 5e4;

424

// Field[ IFI + 2].LcMax = 2e5;

425

//

426

// Background Field = IFI + 2;

427

// Dont extent the elements sizes from the boundary inside the domain

428

//Mesh.CharacteristicLengthExtendFromBoundary = 0;

429

430

Field[ IFI + 1] = MathEval;

431

Field[ IFI + 1].F = "1.0E5";

432

Background Field = IFI + 1;

433

434

//Set some options for better png output

435

General.Color.Background = {255,255,255};

436

General.Color.BackgroundGradient = {255,255,255};

437

General.Color.Foreground = Black;

438

Mesh.Color.Lines = {0,0,0};

439

440

General.Trackball = 0 ;

441

General.RotationX = 180;

442

General.RotationY = 0;

443

General.RotationZ = 270;

24

//////////////////////////////////////////////////////////////////////////

25

//

26

// Generation of boundary representation from arbitrary geophysical

27

// fields and initialisation for anisotropic, unstructured meshing.

28

//

29

30

//

31

// This program is free software: you can redistribute it and/or modify

32

// it under the terms of the GNU General Public License as published by

33

// the Free Software Foundation, either version 3 of the License, or

34

// (at your option) any later version.

35

//

36

// This program is distributed in the hope that it will be useful,

37

// but WITHOUT ANY WARRANTY; without even the implied warranty of

38

// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the

39

// GNU General Public License for more details.

40

//

41

// You should have received a copy of the GNU General Public License

42

// along with this program. If not, see <http://www.gnu.org/licenses/>.

43

//

44

//////////////////////////////////////////////////////////////////////////

45