function Arrays_LT


global IEN LM ID B Node BNDY1 BNDY2 BNDY3 BNDY4
 

IEN=[
   1, 148,  88,  89
   2, 101,  87,  88
   3, 170, 144, 143
   4, 254, 141, 200
   5, 363, 278, 362
   6, 138,  82,  83
   7, 293, 200, 138
   8, 135,  85,  86
   9, 146,  65,  66
  10, 114, 112,  40
  11, 104,  60,  61
  12, 106,  18,  19
  13, 111,  19,  20
  14, 218, 110, 109
  15, 174,  39, 112
  16, 215,  38, 174
  17, 217, 176,  35
  18, 177,  34, 176
  19, 178,  24,  25
  20, 284, 216, 107
  21, 285, 183, 177
  22, 183, 179,  32
  23, 180,  31, 179
  24, 181,  29,  30
  25, 107,  28,  29
  26,  28,   2,  29
  27, 182, 178,  25
  28, 177,  33,  34
  29, 184,  23,  24
  30, 302, 215, 174
  31, 106, 105,  18
  32, 110,  21, 109
  33, 112,  39,  40
  34, 259, 114, 188
  35, 214, 106, 186
  36, 188,  41,  42
  37, 115, 104,  61
  38, 356, 259, 265
  39, 129,  47, 121
  40, 118,  60, 104
  41, 118,  59,  60
  42, 119,  58,  59
  43, 230, 130, 129
  44, 307, 190, 225
  45, 128, 123,  56
  46, 124,  55, 123
  47, 224, 192, 125
  48, 126,  54, 124
  49, 192, 191, 125
  50,  52,   3,  53
  51, 193, 192, 124
  52, 120,  57,  58
  53, 194, 118, 104
  54, 117, 116,  44
  55, 223, 195, 116
  56, 132, 115,  62
  57, 133,  64, 103
  58, 196, 101, 148
  59, 333, 213, 263
  60, 135, 102,  85
  61, 233, 134, 142
  62, 137,  84, 102
  63, 137,  83,  84
  64, 199, 137, 102
  65, 334, 254, 293
  66, 200,  81,  82
  67, 201,  69,  70
  68, 143,  80, 141
  69, 144,  79, 143
  70, 145,  87, 101
  71, 142,  67,  68
  72, 147,  78, 144
  73, 202, 144, 170
  74,  77,  76,   4
  75, 150,  89,  90
  76, 113,  16,  17
  77, 152,  90,  91
  78, 149,  16, 113
  79, 287, 218, 185
  80, 151,  15, 149
  81, 214, 187, 105
  82, 153,  14, 151
  83, 204, 113, 187
  84, 155,  13, 153
  85, 205, 149, 204
  86, 157,  12, 155
  87, 206, 151, 205
  88, 159,  10,  11
  89, 245, 155, 207
  90, 161,   9,  10
  91, 209, 159, 208
  92, 163,   8,   9
  93, 209, 161, 159
  94, 167,   7, 165
  95, 211, 167, 165
  96, 169,   6, 167
  97, 212, 167, 211
  98, 100,   1,   5
  99, 113,  17, 105
 100, 134,  66,  67
 101, 263, 136, 201
 102, 102,  84,  85
 103, 131,  43, 116
 104, 103,  64,  65
 105, 130,  45,  46
 106, 115,  61,  62
 107, 105,  17,  18
 108, 186, 106, 111
 109, 109,  22, 108
 110, 174,  38,  39
 111, 216,  27, 107
 112, 107,  27,  28
 113, 108,  22,  23
 114, 175,  36,  37
 115, 365, 301, 305
 116, 109,  21,  22
 117, 185, 109, 108
 118, 110,  20,  21
 119, 111, 106,  19
 120, 111,  20, 110
 121, 114,  40,  41
 122, 336, 266, 311
 123, 242, 150, 152
 124, 149,  15,  16
 125, 116,  43,  44
 126, 198, 103, 146
 127, 132,  62,  63
 128, 265, 131, 195
 129, 117,  44,  45
 130, 232, 133, 103
 131, 229, 118, 194
 132, 189, 117, 130
 133, 122,  48,  49
 134, 119,  59, 118
 135, 127,  49,  50
 136, 120,  58, 119
 137, 122, 121,  48
 138, 128,  57, 120
 139, 129,  46,  47
 140, 121,  47,  48
 141, 224,  52, 126
 142, 227, 127, 191
 143, 123,  55,  56
 144, 228, 226, 128
 145, 224, 126, 192
 146, 124,  54,  55
 147, 127,  50, 125
 148, 125,  50,  51
 149, 126,  52,  53
 150, 126,  53,  54
 151, 266, 193, 226
 152, 127, 122,  49
 153, 128,  56,  57
 154, 337, 311, 310
 155, 229, 222, 119
 156, 190, 129, 121
 157, 130, 117,  45
 158, 130,  46, 129
 159, 188,  42, 131
 160, 131,  42,  43
 161, 133, 132,  63
 162, 308, 223, 189
 163, 133,  63,  64
 164, 271, 197, 233
 165, 136,  68,  69
 166, 213, 142, 136
 167, 145, 135,  86
 168, 241, 101, 196
 169, 139,  70,  71
 170, 241, 145, 101
 171, 236,  72,  73
 172, 138,  83, 137
 173, 140, 139,  71
 174, 200,  82, 138
 175, 201,  70, 139
 176, 140,  71,  72
 177, 237, 236,  73
 178, 270, 140, 236
 179, 141,  80,  81
 180, 171, 170, 143
 181, 142, 134,  67
 182, 142,  68, 136
 183, 240, 202, 239
 184, 143,  79,  80
 185, 240,  76, 147
 186, 144,  78,  79
 187, 145,  86,  87
 188, 319, 196, 294
 189, 146, 103,  65
 190, 146,  66, 134
 191, 147,  76,  77
 192, 147,  77,  78
 193, 148, 101,  88
 194, 154,  91,  92
 195, 203, 196, 148
 196, 151,  14,  15
 197, 150, 148,  89
 198, 156,  92,  93
 199, 272, 205, 204
 200, 153,  13,  14
 201, 152, 150,  90
 202, 158,  93,  94
 203, 243, 242, 152
 204, 155,  12,  13
 205, 154, 152,  91
 206, 160,  94,  95
 207, 207, 153, 206
 208, 157,  11,  12
 209, 156, 154,  92
 210, 321, 244, 275
 211, 274, 205, 272
 212, 159,  11, 157
 213, 158, 156,  93
 214, 162,  95,  96
 215, 248, 247, 160
 216, 161,  10, 159
 217, 160, 158,  94
 218, 164,  96,  97
 219, 343, 272, 295
 220, 163,   9, 161
 221, 162, 160,  95
 222, 166,  97,  98
 223, 210, 163, 161
 224, 165,   8, 163
 225, 164, 162,  96
 226, 276, 209, 208
 227, 330, 277, 252
 228, 165,   7,   8
 229, 166, 164,  97
 230, 168,  98,  99
 231, 253, 169, 212
 232, 167,   6,   7
 233, 168, 166,  98
 234, 253, 100, 169
 235, 169, 100,   5
 236, 169,   5,   6
 237, 171, 143, 141
 238, 318, 278, 237
 239, 172, 138, 137
 240, 278, 270, 236
 241, 235, 102, 135
 242, 279, 139, 238
 243, 186, 111, 173
 244, 173, 111, 110
 245, 215, 175,  37
 246, 280, 112, 219
 247, 304, 281, 217
 248, 217,  36, 175
 249, 384, 365, 305
 250, 176,  34,  35
 251, 382, 381, 340
 252, 183,  33, 177
 253, 184,  24, 178
 254, 182,  25,  26
 255, 256, 182, 216
 256, 179,  31,  32
 257, 284, 181, 180
 258, 180,  30,  31
 259, 181, 107,  29
 260, 181,  30, 180
 261, 216,  26,  27
 262, 216, 182,  26
 263, 183,  32,  33
 264, 303, 284, 180
 265, 283, 178, 182
 266, 184, 108,  23
 267, 341, 185, 288
 268, 288, 108, 184
 269, 315, 198, 269
 270, 187, 113, 105
 271, 243, 152, 154
 272, 204, 149, 113
 273, 188, 114,  41
 274, 195, 131, 116
 275, 260, 230, 190
 276, 223, 117, 189
 277, 225, 121, 122
 278, 229, 119, 118
 279, 191, 127, 125
 280, 290, 192, 193
 281, 192, 126, 124
 282, 193, 124, 123
 283, 226, 123, 128
 284, 228, 120, 222
 285, 221, 104, 115
 286, 220, 115, 132
 287, 438, 356, 357
 288, 265, 188, 131
 289, 442, 375, 377
 290, 203, 148, 150
 291, 198, 146, 197
 292, 197, 146, 134
 293, 231, 220, 132
 294, 232, 103, 198
 295, 199, 172, 137
 296, 292, 235, 234
 297, 200, 141,  81
 298, 293, 138, 172
 299, 201, 136,  69
 300, 279, 201, 139
 301, 239, 237,  74
 302, 202, 147, 144
 303, 244, 154, 156
 304, 343, 295, 342
 305, 205, 151, 149
 306, 291, 187, 214
 307, 274, 206, 205
 308, 206, 153, 151
 309, 246, 244, 156
 310, 207, 155, 153
 311, 275, 244, 246
 312, 245, 208, 157
 313, 247, 246, 158
 314, 208, 159, 157
 315, 251, 249, 164
 316, 210, 161, 209
 317, 327, 251, 277
 318, 250, 163, 210
 319, 252, 251, 166
 320, 250, 211, 165
 321, 253, 212, 168
 322, 212, 169, 167
 323, 300, 299, 199
 324, 300, 235, 292
 325, 214, 105, 106
 326, 301, 291, 214
 327, 215,  37,  38
 328, 280, 174, 112
 329, 284, 256, 216
 330, 303, 179, 257
 331, 217,  35,  36
 332, 304, 175, 258
 333, 218, 173, 110
 334, 218, 109, 185
 335, 380, 302, 280
 336, 219, 112, 114
 337, 370, 289, 260
 338, 231, 132, 133
 339, 221, 194, 104
 340, 221, 115, 220
 341, 309, 191, 290
 342, 222, 120, 119
 343, 223, 116, 117
 344, 289, 189, 230
 345, 224, 125,  51
 346, 224,  51,  52
 347, 225, 190, 121
 348, 227, 122, 127
 349, 226, 193, 123
 350, 228, 128, 120
 351, 227, 225, 122
 352, 336, 309, 266
 353, 403, 346, 267
 354, 310, 228, 222
 355, 264, 194, 221
 356, 312, 229, 194
 357, 230, 189, 130
 358, 230, 129, 190
 359, 232, 231, 133
 360, 316, 232, 315
 361, 269, 198, 197
 362, 377, 273, 297
 363, 233, 197, 134
 364, 233, 142, 213
 365, 234, 135, 145
 366, 241, 234, 145
 367, 235, 199, 102
 368, 235, 135, 234
 369, 236, 140,  72
 370, 237,  73,  74
 371, 239,  74,  75
 372, 317, 170, 255
 373, 238, 139, 140
 374, 435, 347, 331
 375, 240, 239,  75
 376, 318, 202, 317
 377, 240,  75,  76
 378, 240, 147, 202
 379, 294, 203, 273
 380, 405, 261, 333
 381, 242, 203, 150
 382, 246, 156, 158
 383, 244, 243, 154
 384, 349, 274, 338
 385, 247, 158, 160
 386, 397, 350, 323
 387, 245, 157, 155
 388, 296, 206, 274
 389, 248, 160, 162
 390, 325, 248, 298
 391, 249, 162, 164
 392, 327, 249, 251
 393, 249, 248, 162
 394, 326, 209, 276
 395, 251, 164, 166
 396, 400, 327, 277
 397, 250, 165, 163
 398, 352, 328, 326
 399, 252, 166, 168
 400, 329, 250, 328
 401, 252, 168, 212
 402, 330, 212, 211
 403, 253, 168,  99
 404, 253,  99, 100
 405, 254, 171, 141
 406, 374, 293, 172
 407, 255, 170, 171
 408, 353, 171, 254
 409, 257, 179, 183
 410, 283, 182, 256
 411, 282, 176, 281
 412, 354, 257, 285
 413, 383, 355, 304
 414, 258, 175, 215
 415, 259, 219, 114
 416, 357, 265, 195
 417, 312, 194, 264
 418, 289, 230, 260
 419, 271, 233, 261
 420, 261, 233, 213
 421, 305, 173, 262
 422, 262, 173, 218
 423, 263, 213, 136
 424, 279, 263, 201
 425, 306, 221, 220
 426, 314, 220, 231
 427, 265, 259, 188
 428, 357, 195, 268
 429, 309, 290, 266
 430, 311, 226, 228
 431, 307, 225, 267
 432, 267, 225, 227
 433, 268, 195, 223
 434, 422, 370, 386
 435, 314, 306, 220
 436, 271, 269, 197
 437, 270, 238, 140
 438, 278, 236, 237
 439, 407, 359, 402
 440, 294, 196, 203
 441, 297, 242, 243
 442, 295, 204, 187
 443, 360, 319, 359
 444, 273, 203, 242
 445, 361, 322, 350
 446, 296, 207, 206
 447, 321, 297, 243
 448, 324, 246, 247
 449, 322, 208, 245
 450, 323, 322, 245
 451, 277, 251, 252
 452, 328, 210, 326
 453, 362, 318, 317
 454, 363, 270, 278
 455, 408, 353, 334
 456, 364, 263, 279
 457, 339, 262, 287
 458, 384, 262, 339
 459, 282, 177, 176
 460, 281, 176, 217
 461, 467, 340, 429
 462, 410, 287, 341
 463, 286, 184, 178
 464, 286, 178, 283
 465, 284, 107, 181
 466, 303, 180, 179
 467, 285, 257, 183
 468, 285, 177, 282
 469, 288, 184, 286
 470, 367, 283, 340
 471, 368, 366, 282
 472, 287, 262, 218
 473, 288, 185, 108
 474, 369, 286, 367
 475, 386, 370, 260
 476, 345, 308, 289
 477, 290, 191, 192
 478, 290, 193, 266
 479, 301, 214, 186
 480, 372, 219, 332
 481, 411, 364, 404
 482, 320, 234, 241
 483, 293, 254, 200
 484, 374, 172, 299
 485, 395, 297, 321
 486, 359, 319, 294
 487, 295, 272, 204
 488, 295, 187, 291
 489, 323, 207, 296
 490, 425, 335, 365
 491, 297, 273, 242
 492, 321, 243, 244
 493, 398, 298, 351
 494, 298, 248, 249
 495, 364, 279, 347
 496, 299, 172, 199
 497, 300, 199, 235
 498, 393, 320, 360
 499, 305, 186, 173
 500, 372, 280, 219
 501, 302, 258, 215
 502, 302, 174, 280
 503, 381, 257, 354
 504, 382, 284, 303
 505, 304, 217, 175
 506, 383, 258, 344
 507, 305, 301, 186
 508, 440, 413, 372
 509, 358, 306, 314
 510, 306, 264, 221
 511, 414, 336, 337
 512, 307, 260, 190
 513, 308, 268, 223
 514, 308, 189, 289
 515, 309, 267, 227
 516, 309, 227, 191
 517, 313, 222, 229
 518, 480, 452, 469
 519, 311, 266, 226
 520, 311, 228, 310
 521, 385, 264, 306
 522, 337, 336, 311
 523, 313, 310, 222
 524, 313, 229, 312
 525, 316, 231, 232
 526, 391, 269, 376
 527, 315, 232, 198
 528, 376, 269, 271
 529, 457, 392, 420
 530, 316, 314, 231
 531, 317, 202, 170
 532, 362, 317, 255
 533, 318, 237, 239
 534, 318, 239, 202
 535, 320, 241, 319
 536, 319, 241, 196
 537, 360, 320, 319
 538, 320, 292, 234
 539, 432, 416, 343
 540, 396, 275, 324
 541, 416, 338, 343
 542, 322, 276, 208
 543, 323, 245, 207
 544, 349, 296, 274
 545, 324, 275, 246
 546, 325, 324, 247
 547, 325, 247, 248
 548, 398, 378, 325
 549, 326, 210, 209
 550, 418, 352, 399
 551, 327, 298, 249
 552, 400, 351, 327
 553, 419, 400, 329
 554, 328, 250, 210
 555, 330, 329, 277
 556, 329, 211, 250
 557, 330, 252, 212
 558, 330, 211, 329
 559, 347, 238, 331
 560, 331, 238, 270
 561, 332, 219, 259
 562, 443, 390, 391
 563, 373, 299, 300
 564, 333, 261, 213
 565, 353, 254, 334
 566, 441, 374, 299
 567, 342, 291, 335
 568, 335, 291, 301
 569, 389, 310, 313
 570, 403, 309, 336
 571, 387, 313, 312
 572, 469, 414, 455
 573, 375, 359, 294
 574, 338, 274, 272
 575, 423, 354, 366
 576, 424, 423, 366
 577, 340, 283, 256
 578, 366, 285, 282
 579, 341, 287, 185
 580, 369, 288, 286
 581, 342, 295, 291
 582, 482, 431, 371
 583, 551, 492, 496
 584, 343, 338, 272
 585, 344, 258, 302
 586, 449, 384, 339
 587, 385, 306, 358
 588, 371, 308, 345
 589, 386, 307, 346
 590, 346, 307, 267
 591, 347, 279, 238
 592, 421, 408, 334
 593, 376, 271, 348
 594, 348, 271, 261
 595, 471, 466, 426
 596, 397, 296, 349
 597, 398, 325, 298
 598, 350, 322, 323
 599, 351, 298, 327
 600, 427, 419, 352
 601, 399, 326, 276
 602, 399, 276, 361
 603, 472, 255, 353
 604, 353, 255, 171
 605, 429, 354, 423
 606, 368, 281, 355
 607, 355, 281, 304
 608, 473, 448, 450
 609, 356, 332, 259
 610, 357, 356, 265
 611, 371, 268, 308
 612, 431, 268, 371
 613, 392, 314, 316
 614, 390, 316, 315
 615, 375, 294, 273
 616, 432, 342, 412
 617, 474, 407, 402
 618, 379, 373, 300
 619, 361, 276, 322
 620, 396, 394, 275
 621, 362, 278, 318
 622, 472, 428, 255
 623, 475, 428, 472
 624, 363, 331, 270
 625, 364, 333, 263
 626, 434, 331, 363
 627, 365, 335, 301
 628, 384, 305, 262
 629, 436, 369, 409
 630, 366, 354, 285
 631, 367, 286, 283
 632, 429, 381, 354
 633, 368, 282, 281
 634, 476, 423, 424
 635, 409, 369, 367
 636, 369, 341, 288
 637, 370, 345, 289
 638, 386, 260, 307
 639, 431, 357, 268
 640, 415, 385, 358
 641, 380, 280, 372
 642, 440, 332, 401
 643, 411, 405, 333
 644, 379, 300, 292
 645, 374, 334, 293
 646, 441, 299, 373
 647, 377, 375, 273
 648, 445, 377, 395
 649, 444, 376, 348
 650, 477, 405, 411
 651, 394, 321, 275
 652, 460, 395, 394
 653, 396, 324, 378
 654, 378, 324, 325
 655, 393, 379, 292
 656, 458, 360, 407
 657, 513, 413, 478
 658, 380, 344, 302
 659, 382, 340, 256
 660, 381, 303, 257
 661, 382, 256, 284
 662, 382, 303, 381
 663, 383, 304, 258
 664, 448, 430, 383
 665, 479, 448, 473
 666, 413, 380, 372
 667, 454, 451, 388
 668, 388, 264, 385
 669, 414, 346, 403
 670, 480, 422, 452
 671, 389, 313, 387
 672, 388, 387, 312
 673, 388, 312, 264
 674, 454, 385, 415
 675, 453, 387, 451
 676, 389, 337, 310
 677, 392, 316, 390
 678, 391, 390, 315
 679, 391, 315, 269
 680, 548, 490, 489
 681, 482, 371, 439
 682, 392, 358, 314
 683, 393, 292, 320
 684, 407, 360, 359
 685, 459, 426, 350
 686, 493, 442, 445
 687, 395, 377, 297
 688, 395, 321, 394
 689, 527, 465, 464
 690, 462, 396, 461
 691, 397, 323, 296
 692, 417, 349, 338
 693, 485, 427, 418
 694, 465, 398, 464
 695, 399, 352, 326
 696, 466, 399, 361
 697, 400, 277, 329
 698, 419, 329, 328
 699, 447, 365, 384
 700, 401, 332, 356
 701, 402, 359, 375
 702, 579, 495, 553
 703, 403, 267, 309
 704, 414, 403, 336
 705, 411, 333, 364
 706, 404, 364, 347
 707, 433, 373, 379
 708, 405, 348, 261
 709, 439, 345, 406
 710, 406, 345, 370
 711, 446, 433, 379
 712, 550, 470, 546
 713, 434, 363, 428
 714, 421, 334, 374
 715, 488, 423, 476
 716, 436, 341, 369
 717, 437, 355, 430
 718, 410, 339, 287
 719, 499, 435, 434
 720, 489, 404, 435
 721, 412, 342, 335
 722, 470, 438, 431
 723, 450, 380, 413
 724, 497, 412, 492
 725, 452, 422, 386
 726, 455, 337, 389
 727, 457, 358, 392
 728, 420, 392, 390
 729, 483, 474, 402
 730, 417, 338, 416
 731, 460, 445, 395
 732, 463, 349, 417
 733, 426, 361, 350
 734, 464, 351, 427
 735, 419, 328, 352
 736, 427, 352, 418
 737, 456, 391, 376
 738, 552, 493, 522
 739, 557, 533, 499
 740, 490, 404, 489
 741, 422, 406, 370
 742, 452, 386, 346
 743, 424, 366, 368
 744, 437, 368, 355
 745, 430, 355, 383
 746, 511, 410, 436
 747, 425, 412, 335
 748, 449, 447, 384
 749, 459, 350, 397
 750, 486, 416, 432
 751, 427, 351, 400
 752, 427, 400, 419
 753, 428, 363, 362
 754, 428, 362, 255
 755, 488, 467, 429
 756, 429, 340, 381
 757, 448, 383, 344
 758, 514, 447, 449
 759, 438, 357, 431
 760, 439, 371, 345
 761, 432, 343, 342
 762, 627, 594, 612
 763, 506, 405, 477
 764, 446, 379, 393
 765, 498, 421, 494
 766, 499, 434, 475
 767, 435, 404, 347
 768, 435, 331, 434
 769, 436, 410, 341
 770, 534, 488, 476
 771, 535, 512, 479
 772, 437, 424, 368
 773, 482, 470, 431
 774, 438, 401, 356
 775, 518, 454, 415
 776, 581, 480, 539
 777, 440, 372, 332
 778, 491, 401, 438
 779, 441, 421, 374
 780, 505, 373, 433
 781, 442, 402, 375
 782, 445, 442, 377
 783, 443, 420, 390
 784, 456, 376, 444
 785, 506, 348, 405
 786, 562, 540, 506
 787, 462, 394, 396
 788, 495, 417, 416
 789, 458, 446, 393
 790, 509, 507, 446
 791, 447, 425, 365
 792, 510, 339, 410
 793, 450, 344, 380
 794, 450, 448, 344
 795, 537, 476, 501
 796, 479, 430, 448
 797, 513, 450, 413
 798, 538, 537, 501
 799, 451, 387, 388
 800, 518, 415, 503
 801, 452, 346, 414
 802, 516, 469, 455
 803, 453, 389, 387
 804, 589, 581, 539
 805, 454, 388, 385
 806, 520, 503, 457
 807, 455, 414, 337
 808, 455, 389, 453
 809, 540, 444, 506
 810, 456, 443, 391
 811, 457, 415, 358
 812, 487, 420, 443
 813, 458, 393, 360
 814, 521, 407, 474
 815, 463, 459, 397
 816, 541, 471, 524
 817, 461, 396, 378
 818, 623, 576, 577
 819, 465, 378, 398
 820, 527, 464, 485
 821, 523, 462, 461
 822, 462, 460, 394
 823, 463, 397, 349
 824, 495, 416, 486
 825, 464, 398, 351
 826, 528, 527, 485
 827, 528, 466, 471
 828, 465, 461, 378
 829, 466, 418, 399
 830, 466, 361, 426
 831, 467, 409, 367
 832, 467, 367, 340
 833, 502, 406, 468
 834, 468, 406, 422
 835, 469, 452, 414
 836, 516, 455, 453
 837, 519, 487, 443
 838, 491, 438, 470
 839, 524, 426, 459
 840, 508, 459, 463
 841, 472, 353, 408
 842, 475, 472, 408
 843, 478, 413, 440
 844, 496, 492, 425
 845, 483, 402, 442
 846, 492, 412, 425
 847, 475, 434, 428
 848, 494, 421, 441
 849, 500, 409, 488
 850, 501, 424, 437
 851, 490, 411, 404
 852, 507, 433, 446
 853, 504, 440, 401
 854, 583, 481, 540
 855, 512, 430, 479
 856, 568, 513, 478
 857, 517, 451, 515
 858, 480, 468, 422
 859, 519, 456, 481
 860, 481, 456, 444
 861, 502, 439, 406
 862, 574, 559, 520
 863, 522, 445, 460
 864, 526, 417, 495
 865, 523, 461, 484
 866, 484, 461, 465
 867, 485, 464, 427
 868, 528, 418, 466
 869, 497, 432, 412
 870, 568, 555, 513
 871, 558, 482, 439
 872, 520, 420, 487
 873, 488, 409, 467
 874, 488, 429, 423
 875, 547, 498, 494
 876, 584, 533, 547
 877, 585, 548, 584
 878, 490, 477, 411
 879, 550, 549, 491
 880, 558, 546, 482
 881, 496, 425, 447
 882, 544, 496, 514
 883, 493, 483, 442
 884, 525, 460, 462
 885, 505, 441, 373
 886, 505, 494, 441
 887, 578, 522, 525
 888, 556, 553, 486
 889, 554, 449, 510
 890, 568, 478, 545
 891, 556, 497, 529
 892, 497, 486, 432
 893, 498, 475, 408
 894, 498, 408, 421
 895, 499, 489, 435
 896, 557, 475, 498
 897, 511, 436, 500
 898, 500, 436, 409
 899, 501, 476, 424
 900, 512, 501, 437
 901, 558, 439, 502
 902, 518, 515, 454
 903, 503, 415, 457
 904, 586, 559, 574
 905, 504, 478, 440
 906, 504, 401, 491
 907, 560, 433, 507
 908, 584, 547, 561
 909, 506, 444, 348
 910, 562, 477, 530
 911, 509, 446, 458
 912, 583, 573, 481
 913, 564, 524, 508
 914, 526, 508, 463
 915, 521, 509, 458
 916, 565, 563, 509
 917, 510, 449, 339
 918, 511, 510, 410
 919, 566, 500, 534
 920, 569, 511, 566
 921, 512, 437, 430
 922, 587, 566, 537
 923, 513, 473, 450
 924, 545, 478, 504
 925, 514, 496, 447
 926, 659, 634, 588
 927, 539, 480, 516
 928, 515, 451, 454
 929, 516, 480, 469
 930, 517, 516, 453
 931, 517, 453, 451
 932, 570, 539, 517
 933, 571, 503, 559
 934, 589, 539, 570
 935, 519, 443, 456
 936, 602, 531, 552
 937, 520, 457, 420
 938, 574, 487, 543
 939, 521, 458, 407
 940, 532, 474, 483
 941, 522, 493, 445
 942, 577, 523, 484
 943, 580, 577, 484
 944, 591, 564, 542
 945, 524, 471, 426
 946, 524, 459, 508
 947, 525, 522, 460
 948, 525, 462, 523
 949, 526, 463, 417
 950, 665, 602, 624
 951, 593, 528, 471
 952, 527, 484, 465
 953, 528, 485, 418
 954, 593, 580, 528
 955, 573, 519, 481
 956, 529, 497, 492
 957, 562, 506, 477
 958, 530, 477, 490
 959, 532, 483, 531
 960, 531, 483, 493
 961, 662, 590, 628
 962, 575, 474, 532
 963, 533, 489, 499
 964, 557, 547, 533
 965, 534, 500, 488
 966, 538, 512, 535
 967, 567, 479, 473
 968, 555, 473, 513
 969, 550, 491, 470
 970, 536, 502, 468
 971, 566, 534, 537
 972, 537, 534, 476
 973, 538, 501, 512
 974, 616, 554, 569
 975, 539, 516, 517
 976, 570, 517, 515
 977, 563, 507, 509
 978, 540, 481, 444
 979, 542, 508, 526
 980, 577, 576, 523
 981, 579, 526, 495
 982, 553, 495, 486
 983, 620, 574, 543
 984, 543, 487, 519
 985, 551, 496, 544
 986, 595, 514, 554
 987, 643, 596, 615
 988, 549, 504, 491
 989, 546, 470, 482
 990, 652, 598, 607
 991, 561, 494, 505
 992, 560, 505, 433
 993, 548, 530, 490
 994, 548, 489, 533
 995, 598, 550, 546
 996, 549, 545, 504
 997, 607, 558, 606
 998, 599, 550, 598
 999, 551, 529, 492
1000, 614, 600, 568
1001, 552, 531, 493
1002, 578, 525, 576
1003, 631, 591, 625
1004, 602, 590, 531
1005, 554, 514, 449
1006, 569, 510, 511
1007, 604, 555, 600
1008, 604, 603, 567
1009, 556, 486, 497
1010, 596, 529, 551
1011, 557, 499, 475
1012, 557, 498, 547
1013, 606, 502, 536
1014, 606, 558, 502
1015, 617, 608, 571
1016, 559, 503, 520
1017, 609, 561, 560
1018, 610, 507, 563
1019, 561, 547, 494
1020, 561, 505, 560
1021, 632, 609, 610
1022, 609, 584, 561
1023, 565, 509, 521
1024, 612, 594, 565
1025, 564, 541, 524
1026, 564, 508, 542
1027, 575, 565, 521
1028, 633, 583, 619
1029, 566, 511, 500
1030, 634, 537, 538
1031, 567, 535, 479
1032, 567, 473, 555
1033, 614, 545, 597
1034, 635, 544, 595
1035, 587, 569, 566
1036, 569, 554, 510
1037, 572, 515, 518
1038, 572, 518, 571
1039, 571, 518, 503
1040, 617, 559, 586
1041, 668, 645, 636
1042, 572, 570, 515
1043, 573, 543, 519
1044, 619, 540, 562
1045, 574, 520, 487
1046, 613, 543, 573
1047, 575, 521, 474
1048, 621, 532, 590
1049, 576, 525, 523
1050, 622, 564, 591
1051, 580, 484, 527
1052, 626, 623, 577
1053, 582, 578, 576
1054, 578, 552, 522
1055, 625, 579, 592
1056, 579, 542, 526
1057, 580, 527, 528
1058, 626, 541, 622
1059, 581, 536, 468
1060, 581, 468, 480
1061, 637, 591, 631
1062, 624, 578, 582
1063, 611, 530, 585
1064, 594, 563, 565
1065, 638, 609, 632
1066, 584, 548, 533
1067, 611, 562, 530
1068, 585, 530, 548
1069, 671, 644, 656
1070, 599, 549, 550
1071, 588, 538, 535
1072, 616, 569, 587
1073, 603, 535, 567
1074, 660, 616, 659
1075, 589, 536, 581
1076, 618, 570, 572
1077, 590, 532, 531
1078, 641, 592, 601
1079, 625, 542, 579
1080, 592, 579, 553
1081, 601, 553, 556
1082, 605, 529, 596
1083, 593, 471, 541
1084, 626, 593, 541
1085, 675, 642, 655
1086, 612, 565, 575
1087, 595, 544, 514
1088, 616, 595, 554
1089, 615, 551, 544
1090, 670, 648, 658
1091, 597, 545, 549
1092, 692, 630, 652
1093, 607, 546, 558
1094, 645, 618, 636
1095, 630, 599, 598
1096, 599, 597, 549
1097, 690, 647, 649
1098, 600, 555, 568
1099, 601, 592, 553
1100, 605, 601, 556
1101, 624, 552, 578
1102, 624, 602, 552
1103, 603, 588, 535
1104, 604, 567, 555
1105, 649, 600, 648
1106, 672, 640, 647
1107, 605, 556, 529
1108, 688, 650, 643
1109, 686, 589, 645
1110, 686, 668, 651
1111, 686, 651, 606
1112, 607, 598, 546
1113, 608, 572, 571
1114, 661, 586, 644
1115, 610, 560, 507
1116, 610, 609, 560
1117, 654, 563, 594
1118, 655, 585, 638
1119, 619, 562, 611
1120, 675, 669, 642
1121, 621, 612, 575
1122, 702, 681, 665
1123, 676, 629, 662
1124, 657, 573, 583
1125, 614, 568, 545
1126, 646, 597, 599
1127, 615, 596, 551
1128, 635, 615, 544
1129, 659, 588, 640
1130, 677, 659, 640
1131, 617, 571, 559
1132, 661, 653, 617
1133, 618, 589, 570
1134, 618, 572, 608
1135, 642, 633, 619
1136, 619, 583, 540
1137, 620, 586, 574
1138, 620, 543, 613
1139, 621, 575, 532
1140, 662, 621, 590
1141, 678, 631, 666
1142, 622, 541, 564
1143, 623, 582, 576
1144, 626, 577, 580
1145, 667, 663, 623
1146, 679, 637, 631
1147, 625, 591, 542
1148, 641, 601, 639
1149, 626, 580, 593
1150, 667, 626, 622
1151, 697, 657, 633
1152, 629, 612, 621
1153, 628, 590, 602
1154, 639, 601, 605
1155, 656, 644, 620
1156, 629, 627, 612
1157, 644, 586, 620
1158, 646, 599, 630
1159, 679, 631, 678
1160, 666, 625, 592
1161, 654, 610, 563
1162, 654, 632, 610
1163, 657, 583, 633
1164, 642, 619, 611
1165, 634, 587, 537
1166, 634, 538, 588
1167, 660, 595, 616
1168, 649, 604, 600
1169, 636, 618, 608
1170, 653, 608, 617
1171, 667, 637, 663
1172, 637, 622, 591
1173, 638, 585, 584
1174, 638, 584, 609
1175, 680, 639, 673
1176, 650, 605, 596
1177, 640, 588, 603
1178, 647, 603, 604
1179, 709, 681, 702
1180, 704, 682, 680
1181, 655, 642, 611
1182, 694, 683, 675
1183, 684, 615, 635
1184, 658, 648, 614
1185, 721, 668, 691
1186, 656, 620, 613
1187, 686, 536, 589
1188, 645, 589, 618
1189, 658, 597, 646
1190, 705, 685, 700
1191, 647, 640, 603
1192, 648, 600, 614
1193, 706, 672, 690
1194, 689, 648, 670
1195, 689, 649, 648
1196, 649, 647, 604
1197, 650, 639, 605
1198, 650, 596, 643
1199, 652, 607, 651
1200, 651, 607, 606
1201, 685, 653, 661
1202, 652, 630, 598
1203, 674, 652, 651
1204, 653, 636, 608
1205, 693, 594, 627
1206, 707, 675, 683
1207, 655, 611, 585
1208, 694, 638, 632
1209, 695, 613, 657
1210, 696, 629, 676
1211, 657, 613, 573
1212, 697, 633, 669
1213, 658, 614, 597
1214, 698, 658, 646
1215, 659, 616, 587
1216, 659, 587, 634
1217, 699, 684, 635
1218, 660, 635, 595
1219, 661, 617, 586
1220, 700, 644, 671
1221, 662, 629, 621
1222, 708, 682, 704
1223, 664, 582, 663
1224, 663, 582, 623
1225, 727, 708, 709
1226, 664, 624, 582
1227, 665, 628, 602
1228, 665, 624, 664
1229, 666, 631, 625
1230, 666, 592, 641
1231, 667, 623, 626
1232, 667, 622, 637
1233, 674, 651, 668
1234, 691, 636, 653
1235, 669, 633, 642
1236, 726, 676, 701
1237, 687, 646, 630
1238, 688, 673, 650
1239, 711, 656, 695
1240, 728, 704, 714
1241, 719, 690, 689
1242, 699, 660, 677
1243, 673, 639, 650
1244, 765, 760, 762
1245, 692, 652, 674
1246, 721, 713, 692
1247, 694, 675, 655
1248, 683, 632, 654
1249, 723, 711, 695
1250, 701, 662, 628
1251, 677, 660, 659
1252, 677, 640, 672
1253, 702, 664, 679
1254, 682, 666, 641
1255, 679, 664, 663
1256, 679, 663, 637
1257, 680, 641, 639
1258, 745, 744, 726
1259, 701, 628, 681
1260, 681, 628, 665
1261, 682, 678, 666
1262, 682, 641, 680
1263, 693, 654, 594
1264, 693, 683, 654
1265, 684, 643, 615
1266, 699, 635, 660
1267, 700, 661, 644
1268, 729, 728, 714
1269, 686, 606, 536
1270, 686, 645, 668
1271, 730, 721, 716
1272, 716, 685, 705
1273, 717, 643, 684
1274, 718, 698, 715
1275, 720, 719, 689
1276, 720, 670, 718
1277, 690, 672, 647
1278, 690, 649, 689
1279, 691, 668, 636
1280, 691, 653, 685
1281, 692, 687, 630
1282, 721, 691, 716
1283, 722, 627, 696
1284, 707, 683, 693
1285, 694, 655, 638
1286, 694, 632, 683
1287, 695, 656, 613
1288, 743, 731, 724
1289, 732, 723, 731
1290, 696, 627, 629
1291, 697, 695, 657
1292, 724, 723, 697
1293, 698, 670, 658
1294, 698, 646, 687
1295, 699, 677, 672
1296, 706, 699, 672
1297, 700, 685, 661
1298, 725, 700, 671
1299, 703, 701, 681
1300, 701, 676, 662
1301, 702, 665, 664
1302, 727, 679, 678
1303, 733, 704, 728
1304, 726, 701, 703
1305, 714, 680, 673
1306, 712, 673, 688
1307, 734, 671, 711
1308, 735, 730, 705
1309, 720, 689, 670
1310, 706, 684, 699
1311, 707, 669, 675
1312, 722, 707, 693
1313, 758, 745, 751
1314, 708, 678, 682
1315, 709, 703, 681
1316, 727, 709, 702
1317, 732, 711, 723
1318, 710, 696, 676
1319, 711, 671, 656
1320, 723, 695, 697
1321, 736, 688, 717
1322, 737, 713, 730
1323, 715, 687, 713
1324, 713, 687, 692
1325, 714, 704, 680
1326, 714, 673, 712
1327, 747, 730, 735
1328, 715, 698, 687
1329, 716, 691, 685
1330, 730, 716, 705
1331, 717, 688, 643
1332, 717, 684, 706
1333, 748, 741, 740
1334, 718, 670, 698
1335, 739, 706, 719
1336, 719, 706, 690
1337, 740, 718, 738
1338, 755, 738, 737
1339, 721, 692, 674
1340, 721, 674, 668
1341, 722, 693, 627
1342, 742, 696, 710
1343, 724, 697, 669
1344, 743, 742, 731
1345, 724, 669, 707
1346, 743, 707, 722
1347, 759, 752, 756
1348, 725, 705, 700
1349, 726, 710, 676
1350, 751, 709, 733
1351, 727, 702, 679
1352, 727, 678, 708
1353, 729, 714, 712
1354, 752, 750, 732
1355, 753, 736, 748
1356, 760, 757, 729
1357, 730, 713, 721
1358, 735, 705, 725
1359, 756, 742, 710
1360, 731, 723, 724
1361, 744, 710, 726
1362, 734, 711, 732
1363, 733, 709, 708
1364, 733, 708, 704
1365, 746, 744, 745
1366, 734, 725, 671
1367, 754, 725, 734
1368, 755, 737, 747
1369, 736, 712, 688
1370, 739, 736, 717
1371, 738, 715, 737
1372, 737, 715, 713
1373, 748, 736, 741
1374, 738, 718, 715
1375, 739, 717, 706
1376, 741, 719, 720
1377, 741, 736, 739
1378, 740, 720, 718
1379, 741, 739, 719
1380, 741, 720, 740
1381, 742, 722, 696
1382, 756, 752, 742
1383, 743, 724, 707
1384, 743, 722, 742
1385, 745, 726, 703
1386, 757, 728, 729
1387, 751, 745, 703
1388, 758, 728, 757
1389, 749, 748, 740
1390, 754, 734, 750
1391, 747, 737, 730
1392, 761, 750, 759
1393, 749, 740, 738
1394, 760, 753, 748
1395, 766, 759, 746
1396, 767, 766, 765
1397, 750, 734, 732
1398, 752, 732, 731
1399, 751, 703, 709
1400, 758, 733, 728
1401, 752, 731, 742
1402, 759, 756, 744
1403, 753, 729, 712
1404, 753, 712, 736
1405, 763, 754, 761
1406, 754, 735, 725
1407, 755, 749, 738
1408, 763, 735, 754
1409, 756, 710, 744
1410, 759, 744, 746
1411, 764, 758, 757
1412, 760, 729, 753
1413, 758, 746, 745
1414, 758, 751, 733
1415, 766, 764, 765
1416, 759, 750, 752
1417, 760, 748, 749
1418, 762, 749, 755
1419, 763, 762, 755
1420, 761, 754, 750
1421, 762, 760, 749
1422, 763, 755, 747
1423, 763, 747, 735
1424, 767, 763, 761
1425, 764, 746, 758
1426, 765, 764, 757
1427, 765, 757, 760
1428, 767, 765, 762
1429, 766, 761, 759
1430, 766, 746, 764
1431, 767, 762, 763
1432, 767, 761, 766
];

% number of nodes
Sz1=size(Node);             

% fixed nodes
Sz2=size(BNDY4);

% build ID array (degree of freedom)
ID_t1=ones(3,Sz1(1));
for i=1:Sz2(2)
    ID_t1(:,BNDY4(i))=0;                % the nodes on Boundary 4 are fixed (boundary condition)
end

Sz3=size(BNDY1);                        % X coordinates of Boundary 1 and Boundarys 2 are fixed
for i=1:Sz3(2)
    ID_t1(1,BNDY1(i))=0;
end

Sz4=size(BNDY2);                        % X coordinates of Boundary 1 and Boundarys 2 are fixed
for i=1:Sz4(2)
    ID_t1(1,BNDY2(i))=0;
end

Sz5=size(BNDY3);                        % Y coordinates of Boundary 3 are fixed
for i=1:Sz5(2)
    ID_t1(2,BNDY3(i))=0;
end

ID_t2=reshape(ID_t1,1,3*Sz1(1));        % reshape ID_t1 into a row vector

ID_t3=zeros(1,3*Sz1(1));
A=find(ID_t2);
B=sum(ID_t2);                           % B is the total number of degree of freedom  
for i=1:B
    ID_t3(1,A(i))=i;
end
ID=reshape(ID_t3,3,Sz1(1));             % construct ID array by reshaping ID_t3 back to 3xN matrix
% ID(:,1:400)

% LM array
Sz3=size(IEN);
LM=zeros(9,Sz3(1));
for i=1:Sz3(1)
    LM(1:3,i)=ID(:,IEN(i,2));
    LM(4:6,i)=ID(:,IEN(i,3));
    LM(7:9,i)=ID(:,IEN(i,4));
end
% LM

disp('Arrays OK!')