-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathOld_Parser_For_Aims_Clusters.nb
15820 lines (15715 loc) · 853 KB
/
Old_Parser_For_Aims_Clusters.nb
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 10.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 873767, 15812]
NotebookOptionsPosition[ 865963, 15681]
NotebookOutlinePosition[ 866316, 15697]
CellTagsIndexPosition[ 866273, 15694]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["Parser for FHI-AIMS LCAO coefficients file of cluster", "Title",
CellChangeTimes->{{3.756209670170629*^9,
3.7562097009027157`*^9}},ExpressionUUID->"98d6c56c-bf9b-44a9-a98a-\
14a73f125b99"],
Cell["\<\
If the AIMS calculation is not performed in periodic boundary conditions, \
then the whole LCAO coef. file is crazily parsed. This is an old way how to \
put it into \[OpenCurlyDoubleQuote]Fireball-like\[CloseCurlyDoubleQuote] \
files. For these calculations it was also necessary to print basis to get +/- \
coefficients of the basis far-away from nuclei.
Modern parser in PP-STM working for PBC-calculations do all the +/- coef. \
calculations automatically.\
\>", "Text",
CellChangeTimes->{{3.756209711804784*^9,
3.756209954254887*^9}},ExpressionUUID->"58188be7-044f-4982-8c14-\
fb2116141f5e"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"path", "=",
RowBox[{"(*",
RowBox[{"\"\<~/Desktop/Luna/AIMS_mol/TOAT_old/Cluster_lda/\>\"", ";"}],
"*)"}], "\"\<~/Desktop/WORK/probe_particle_code/TAP_cluster/\>\""}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"inp", "=",
RowBox[{"Import", "[",
RowBox[{
RowBox[{"path", "<>", "\"\</KS_eigenvectors_sh.out\>\""}], ",",
"\"\<Table\>\""}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"eigen", "=",
RowBox[{"Drop", "[",
RowBox[{
RowBox[{"Import", "[",
RowBox[{
RowBox[{"path", "<>", "\"\</all_states_sh.dat\>\""}], ",",
"\"\<Table\>\"", ",",
RowBox[{"HeaderLines", "\[Rule]", "1"}]}], "]"}], ",",
RowBox[{"-", "2"}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{"maxOrb", "=",
RowBox[{"Max", "[",
RowBox[{"eigen", "[",
RowBox[{"[",
RowBox[{"All", ",", "1"}], "]"}], "]"}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"geom", "=",
RowBox[{
RowBox[{"Import", "[",
RowBox[{
RowBox[{"path", "<>", "\"\</geometry.in\>\""}], ",", "\"\<Table\>\""}],
"]"}], "[",
RowBox[{"[",
RowBox[{"All", ",", "5"}], "]"}], "]"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.646656679692748*^9, 3.646656845000211*^9}, {
3.646656898548725*^9, 3.6466569216707153`*^9}, {3.6466569641974087`*^9,
3.646657013793148*^9}, {3.6466572475712967`*^9, 3.646657248016695*^9}, {
3.646659492904944*^9, 3.646659560185162*^9}, 3.648994406698945*^9,
3.649410462801566*^9, {3.6648684596985407`*^9, 3.6648684609473867`*^9}, {
3.670658962439204*^9, 3.670658978577338*^9}, {3.6706591496539307`*^9,
3.6706591551637793`*^9}, {3.671087081995823*^9, 3.67108708272377*^9}, {
3.67108765372051*^9, 3.6710876566929626`*^9}, {3.671088419788949*^9,
3.671088420246943*^9}, {3.671089692904976*^9, 3.671089693430003*^9}, {
3.671878060313167*^9, 3.671878084320608*^9}, 3.672041550249308*^9,
3.672061291419229*^9, {3.672062882429625*^9, 3.672062907931217*^9},
3.6720630443847923`*^9, {3.672137023635724*^9, 3.6721370288936567`*^9}, {
3.681708066170299*^9, 3.681708077601798*^9}, {3.68170890733045*^9,
3.6817089098263483`*^9}, {3.681709544507606*^9, 3.68170954661793*^9}, {
3.6817241504995327`*^9, 3.6817241629307613`*^9}, 3.6817243240774937`*^9, {
3.681726165324143*^9, 3.681726166835465*^9}, {3.756209621343478*^9,
3.756209628404077*^9}},ExpressionUUID->"f0c29b07-f192-41d0-8ed2-\
968acf8fd28f"],
Cell[BoxData["13400"], "Output",
CellChangeTimes->{
3.682755374898778*^9},ExpressionUUID->"178be7d2-8959-434a-ad92-\
584826c765f0"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"(*",
RowBox[{"maxOrb", "=", "1500"}], "*)"}]], "Input",
CellChangeTimes->{{3.6817245485198317`*^9, 3.6817245579221897`*^9},
3.681726157375598*^9, {3.6827552497271223`*^9, 3.682755251494546*^9}, {
3.682762389624846*^9,
3.682762391031801*^9}},ExpressionUUID->"fd210b32-9345-46ad-8b6c-\
bdf0856a2504"],
Cell[BoxData["1500"], "Output",
CellChangeTimes->{
3.6827553750407887`*^9},ExpressionUUID->"5eaa4a64-a671-477a-8bd9-\
06372368d059"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"numBasis", "=", "13552"}], ";"}],
RowBox[{"(*",
RowBox[{"|",
RowBox[{"Total", " ", "number", " ", "of", " ", "basis", " ",
RowBox[{"functions", ":"}]}]}], "*)"}]}]], "Input",
CellChangeTimes->{{3.646657038520666*^9, 3.646657075711282*^9},
3.649410488777108*^9, {3.670658518975626*^9, 3.670658532269808*^9}, {
3.672041600853631*^9, 3.6720416016292667`*^9}, {3.672061276229604*^9,
3.672061276930915*^9}, {3.672063171292104*^9, 3.672063171745957*^9}, {
3.67206324143327*^9, 3.672063242146289*^9}, {3.6817081111792727`*^9,
3.6817081116696033`*^9}, {3.6817089334407873`*^9, 3.681708934310667*^9}, {
3.6817095579476147`*^9, 3.681709558667387*^9}, {3.681724567497469*^9,
3.681724568585532*^9}, {3.682755302268689*^9,
3.682755304100643*^9}},ExpressionUUID->"4587a317-f23d-4eb1-b517-\
910dbf62acb6"],
Cell[BoxData[{
RowBox[{
RowBox[{"header", "=",
RowBox[{"inp", "[",
RowBox[{"[",
RowBox[{
RowBox[{"Range", "[",
RowBox[{"2", ",",
RowBox[{"numBasis", "+", "1"}]}], "]"}], ",",
RowBox[{"Range", "[",
RowBox[{"2", ",", "6"}], "]"}]}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"coefsi", "[", "i_", "]"}], ":=",
RowBox[{
RowBox[{
RowBox[{"Quotient", "[",
RowBox[{
RowBox[{"i", "-", "1"}], ",", "5"}], "]"}], "*",
RowBox[{"(",
RowBox[{"numBasis", "+", "2"}], ")"}]}], "+",
RowBox[{"Range", "[",
RowBox[{"2", ",",
RowBox[{"numBasis", "+", "1"}]}], "]"}]}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"coefsj", "[", "i_", "]"}], ":=",
RowBox[{"7", "+",
RowBox[{"Mod", "[",
RowBox[{
RowBox[{"i", "-", "1"}], ",", "5"}], "]"}]}]}]}], "Input",
CellChangeTimes->{{3.646657091987279*^9, 3.646657093098833*^9}, {
3.646657217818761*^9, 3.646657233201927*^9}, {3.646657284105816*^9,
3.6466573562596817`*^9}, {3.6466574128827953`*^9, 3.646657413289793*^9}, {
3.646657482506587*^9, 3.6466575582417393`*^9}, 3.646657698686228*^9, {
3.646657815729519*^9, 3.6466578697043743`*^9}, {3.646657916699012*^9,
3.646657922658945*^9}, {3.646657962398122*^9, 3.646658007161848*^9}, {
3.646658069051629*^9,
3.646658089059417*^9}},ExpressionUUID->"924ce0cd-4a3a-433f-b5dd-\
c1c0dd1a03b5"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"For", "[",
RowBox[{
RowBox[{"j", "=", "1"}], ",",
RowBox[{"j", "\[LessEqual]",
RowBox[{"Length", "[", "eigen", "]"}]}], ",",
RowBox[{"j", "++"}], ",",
RowBox[{"If", "[",
RowBox[{
RowBox[{
RowBox[{"eigen", "[",
RowBox[{"[",
RowBox[{"j", ",", "4"}], "]"}], "]"}], ">",
RowBox[{"-", "8.583"}]}],
RowBox[{"(*",
RowBox[{"-", "10"}], "*)"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"n1", "=", "j"}], ",",
RowBox[{"Print", "[", "j", "]"}], ",",
RowBox[{"Break", "[", "]"}]}], "}"}]}], "]"}]}],
"]"}], "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"j", "=", "n1"}], ",",
RowBox[{"j", "\[LessEqual]",
RowBox[{"Length", "[", "eigen", "]"}]}], ",",
RowBox[{"j", "++"}], ",",
RowBox[{"If", "[",
RowBox[{
RowBox[{
RowBox[{"eigen", "[",
RowBox[{"[",
RowBox[{"j", ",", "4"}], "]"}], "]"}], ">",
RowBox[{"-", "2.583"}]}],
RowBox[{"(*", "5", "*)"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"n2", "=", "j"}], ",",
RowBox[{"Print", "[", "j", "]"}], ",",
RowBox[{"Break", "[", "]"}]}], "}"}], ",",
RowBox[{"n2", "=", "numBasis"}]}], "]"}]}], "]"}], "\[IndentingNewLine]",
RowBox[{"n2", "-", "n1", "+", "1"}]}], "Input",
CellChangeTimes->{{3.646658161789647*^9, 3.6466582181114683`*^9}, {
3.6466582495041847`*^9, 3.646658343742758*^9}, {3.646658407879977*^9,
3.646658412368754*^9}, {3.646658459596148*^9, 3.646658493539097*^9}, {
3.646658980505591*^9, 3.6466589854564333`*^9}, {3.681724584339113*^9,
3.681724598299439*^9}, {3.6817246327158613`*^9, 3.681724639395372*^9}, {
3.68275550937494*^9,
3.682755514569496*^9}},ExpressionUUID->"37440310-43d9-4f0a-ad48-\
10f87d76c2af"],
Cell[CellGroupData[{
Cell[BoxData["459"], "Print",
CellChangeTimes->{
3.646658519433031*^9, 3.646658985900682*^9, 3.648380566277824*^9,
3.64899442309718*^9, 3.649059442011338*^9, 3.649410496211339*^9,
3.6494262470766277`*^9, 3.664868473924695*^9, 3.671087194773745*^9,
3.6710877609720297`*^9, 3.6710883537251577`*^9, 3.6710885165415173`*^9,
3.671089799605638*^9, 3.671878161201194*^9, 3.6718782425839787`*^9,
3.6720399862932453`*^9, 3.672040106391234*^9, 3.672041611715839*^9,
3.672061371644594*^9, 3.6720614565934553`*^9, {3.672063175713469*^9,
3.672063200117035*^9}, 3.6720632481345453`*^9, 3.672137296262094*^9,
3.6817081172658873`*^9, 3.681708953538452*^9, 3.681709641413719*^9,
3.6817245732933493`*^9, {3.681724641281295*^9, 3.681724649449624*^9},
3.682755519318467*^9},ExpressionUUID->"1ba0e049-b5da-4026-9e09-\
840c31e50c05"],
Cell[BoxData["1304"], "Print",
CellChangeTimes->{
3.646658519433031*^9, 3.646658985900682*^9, 3.648380566277824*^9,
3.64899442309718*^9, 3.649059442011338*^9, 3.649410496211339*^9,
3.6494262470766277`*^9, 3.664868473924695*^9, 3.671087194773745*^9,
3.6710877609720297`*^9, 3.6710883537251577`*^9, 3.6710885165415173`*^9,
3.671089799605638*^9, 3.671878161201194*^9, 3.6718782425839787`*^9,
3.6720399862932453`*^9, 3.672040106391234*^9, 3.672041611715839*^9,
3.672061371644594*^9, 3.6720614565934553`*^9, {3.672063175713469*^9,
3.672063200117035*^9}, 3.6720632481345453`*^9, 3.672137296262094*^9,
3.6817081172658873`*^9, 3.681708953538452*^9, 3.681709641413719*^9,
3.6817245732933493`*^9, {3.681724641281295*^9, 3.681724649449624*^9},
3.682755519319372*^9},ExpressionUUID->"31c80240-b856-49c2-8791-\
1a9cb5a049e8"]
}, Open ]],
Cell[BoxData["846"], "Output",
CellChangeTimes->{
3.646658985903287*^9, 3.648380566279875*^9, 3.6489944230987186`*^9,
3.649059442013249*^9, 3.649410496212945*^9, 3.6494262470786552`*^9,
3.664868473925993*^9, 3.671087194776494*^9, 3.671087760973447*^9,
3.6710883537286243`*^9, 3.6710885165438757`*^9, 3.671089799607712*^9,
3.671878161249506*^9, 3.671878242586255*^9, 3.6720399862961483`*^9,
3.672040106394456*^9, 3.672041611717671*^9, 3.672061371646139*^9,
3.672061456595167*^9, {3.672063175716152*^9, 3.6720632001189127`*^9},
3.6720632481367197`*^9, 3.67213729626373*^9, 3.6817081172674522`*^9,
3.681708953540711*^9, 3.681709641418215*^9, 3.681724573294794*^9, {
3.681724641284296*^9, 3.681724649451097*^9},
3.682755519320545*^9},ExpressionUUID->"0c2b6797-814b-49b5-8e69-\
7472c3437d5c"]
}, Open ]],
Cell[BoxData[{
RowBox[{
RowBox[{"lcaocoefs", "=",
RowBox[{"Array", "[",
RowBox[{"Null", ",",
RowBox[{"{",
RowBox[{"n2", "-", "n1", "+", "1"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"j", "=", "n1"}], ",",
RowBox[{"j", "\[LessEqual]", "n2"}], ",",
RowBox[{"j", "++"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"lcaocoefs", "[",
RowBox[{"[",
RowBox[{"j", "-", "n1", "+", "1"}], "]"}], "]"}], "=",
RowBox[{"inp", "[",
RowBox[{"[",
RowBox[{
RowBox[{"coefsi", "[", "j", "]"}], ",",
RowBox[{"coefsj", "[", "j", "]"}]}], "]"}], "]"}]}], "}"}]}],
"]"}]}], "Input",
CellChangeTimes->{{3.64665687562199*^9, 3.64665687644069*^9}, {
3.64665723817236*^9, 3.646657242085182*^9}, {3.646658531146142*^9,
3.646658545783949*^9}, {3.646658587198885*^9, 3.646658591905816*^9}, {
3.646658624400887*^9, 3.646658854490366*^9}, {3.6466589538299026`*^9,
3.6466589558570957`*^9}, {3.646659065876842*^9, 3.6466590732409353`*^9}, {
3.646659206733293*^9,
3.646659241508234*^9}},ExpressionUUID->"6bc03f13-a5f5-447f-93bc-\
21453af3400f"],
Cell[BoxData[
RowBox[{
RowBox[{"species", "=",
RowBox[{"{",
RowBox[{
"\"\<Au\>\"", ",", "\"\<N\>\"", ",", "\"\<C\>\"", ",", "\"\<H\>\"", ",",
"\"\<O\>\""}], "}"}]}], ";"}]], "Input",
CellChangeTimes->{{3.646659580673628*^9, 3.646659607257201*^9}, {
3.648382152394652*^9, 3.6483821656504793`*^9}, {3.649410508853694*^9,
3.649410512841061*^9}, {3.671087136724626*^9, 3.671087139524714*^9},
3.672063256807166*^9, {3.6817081290596323`*^9, 3.681708131195629*^9}, {
3.681724655777287*^9, 3.681724656120906*^9}, {3.681725944298409*^9,
3.681725997826*^9}},ExpressionUUID->"c2097d50-106e-4e2c-8966-7234aa79ddb1"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"lmax", "=",
RowBox[{"{",
RowBox[{"1", ",", "2", ",", "2", ",", "1", ",", "1", ",", "1", ",", "4"}],
"}"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"o", "[", "1", "]"}], ":=",
"\"\<_s.dat\>\""}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"o", "[", "2", "]"}], ":=",
"\"\<_p.dat\>\""}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"o", "[", "3", "]"}], ":=",
"\"\<_d.dat\>\""}], "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"i", "=", "1"}], ",",
RowBox[{"i", "\[LessEqual]",
RowBox[{"Length", "[", "species", "]"}]}], ",",
RowBox[{"i", "++"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"n", "[", "i_", "]"}], ":=",
RowBox[{"ElementData", "[",
RowBox[{
RowBox[{"species", "[",
RowBox[{"[", "i", "]"}], "]"}], ",", "\"\<Period\>\""}], "]"}]}],
"}"}]}], "]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"name", "[",
RowBox[{"i_", ",", "k_", ",", "l_"}], "]"}], ":=",
RowBox[{"path", "<>", "\"\</at_\>\"", "<>",
RowBox[{"ToString", "[", "i", "]"}], "<>", "\"\<_\>\"", "<>",
RowBox[{"ToString", "[", "k", "]"}], "<>", "\"\<_\>\"", "<>",
RowBox[{"ToString", "[",
RowBox[{"n", "[", "i", "]"}], "]"}], "<>",
RowBox[{"o", "[", "l", "]"}]}]}], ";"}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"i", "=", "1"}], ",",
RowBox[{"i", "\[LessEqual]",
RowBox[{"Length", "[", "species", "]"}]}], ",",
RowBox[{"i", "++"}], ",",
RowBox[{"{",
RowBox[{"For", "[",
RowBox[{
RowBox[{"l", "=", "1"}], ",",
RowBox[{"l", "\[LessEqual]",
RowBox[{"lmax", "[",
RowBox[{"[",
RowBox[{"n", "[", "i", "]"}], "]"}], "]"}]}], ",",
RowBox[{"l", "++"}], ",",
RowBox[{"{",
RowBox[{"(*",
RowBox[{
RowBox[{"Print", "[",
RowBox[{"name", "[",
RowBox[{"i", ",", "I", ",", "l"}], "]"}], "]"}], ","}], "*)"}],
RowBox[{"For", "[",
RowBox[{
RowBox[{"k", "=", "1"}], ",",
RowBox[{"k", "\[LessEqual]", "100"}], ",",
RowBox[{"k", "++"}], ",",
RowBox[{"If", "[",
RowBox[{
RowBox[{"FileExistsQ", "[",
RowBox[{"name", "[",
RowBox[{"i", ",", "k", ",", "l"}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"Print", "[",
RowBox[{"name", "[",
RowBox[{"i", ",", "k", ",", "l"}], "]"}], "]"}], ",",
RowBox[{
RowBox[{"R", "[",
RowBox[{"i", ",", "l"}], "]"}], "=",
RowBox[{"Import", "[",
RowBox[{
RowBox[{"name", "[",
RowBox[{"i", ",", "k", ",", "l"}], "]"}], ",",
"\"\<Table\>\"", ",",
RowBox[{"HeaderLines", "\[Rule]", "1"}]}], "]"}]}]}], "}"}]}],
"]"}]}], "]"}], "}"}]}], "]"}], "}"}]}], "]"}]}], "Input",
CellChangeTimes->{{3.648381114673232*^9, 3.648381122839981*^9}, {
3.648381254367094*^9, 3.6483813013916187`*^9}, {3.6483813562549267`*^9,
3.648381402792494*^9}, {3.648381445432322*^9, 3.648381448800462*^9}, {
3.648381521768354*^9, 3.6483815624325047`*^9}, {3.64838159559304*^9,
3.648381696985634*^9}, {3.6483820471227427`*^9, 3.6483820907433453`*^9}, {
3.6483822415233593`*^9, 3.648382272770693*^9}, {3.6483823097705193`*^9,
3.648382353995221*^9}, {3.648382400425809*^9, 3.648382643032461*^9}, {
3.648382687883821*^9, 3.648382819697509*^9}, {3.6483849191967154`*^9,
3.648385026650432*^9}, {3.6483850904040527`*^9, 3.648385112184224*^9}, {
3.648385198145749*^9, 3.648385209898244*^9}, {3.6483857344267683`*^9,
3.6483857605600452`*^9}, {3.648385909505146*^9, 3.6483859147318296`*^9},
3.671087143706547*^9, {3.6817246584078007`*^9, 3.68172465995295*^9}, {
3.681724967810999*^9, 3.681724976423284*^9}, 3.681725068491808*^9,
3.68172600927094*^9},ExpressionUUID->"60a3ddaa-7b6a-4de5-977b-\
ca1bdb6c11b0"],
Cell[CellGroupData[{
Cell[BoxData["\<\"~/Desktop/WORK/probe_particle_code/TAP_cluster//at_1_6_6_s.\
dat\"\>"], "Print",
CellChangeTimes->{
3.682755603295825*^9},ExpressionUUID->"6ffd8126-a43d-496d-9c93-\
5d4af457ff3a"],
Cell[BoxData["\<\"~/Desktop/WORK/probe_particle_code/TAP_cluster//at_2_16_2_s.\
dat\"\>"], "Print",
CellChangeTimes->{
3.682755603428369*^9},ExpressionUUID->"50070d30-d288-47b3-845f-\
93be7a092413"],
Cell[BoxData["\<\"~/Desktop/WORK/probe_particle_code/TAP_cluster//at_2_17_2_p.\
dat\"\>"], "Print",
CellChangeTimes->{
3.682755603567068*^9},ExpressionUUID->"0afaba39-b004-42f4-a12c-\
c13cb8426b9e"],
Cell[BoxData["\<\"~/Desktop/WORK/probe_particle_code/TAP_cluster//at_3_19_2_s.\
dat\"\>"], "Print",
CellChangeTimes->{
3.682755603705353*^9},ExpressionUUID->"75446690-f2d3-4442-8d7c-\
3037c2762dd0"],
Cell[BoxData["\<\"~/Desktop/WORK/probe_particle_code/TAP_cluster//at_3_20_2_p.\
dat\"\>"], "Print",
CellChangeTimes->{
3.682755603838647*^9},ExpressionUUID->"fec387a8-fab6-4f7a-bfab-\
ff9b0b407911"],
Cell[BoxData["\<\"~/Desktop/WORK/probe_particle_code/TAP_cluster//at_4_21_1_s.\
dat\"\>"], "Print",
CellChangeTimes->{
3.682755603977858*^9},ExpressionUUID->"80834163-2b65-4f56-89aa-\
3e0c020ef10a"],
Cell[BoxData["\<\"~/Desktop/WORK/probe_particle_code/TAP_cluster//at_5_23_2_s.\
dat\"\>"], "Print",
CellChangeTimes->{
3.682755604138974*^9},ExpressionUUID->"0378d341-5be3-4fc2-9b84-\
54523a18fb44"],
Cell[BoxData["\<\"~/Desktop/WORK/probe_particle_code/TAP_cluster//at_5_24_2_p.\
dat\"\>"], "Print",
CellChangeTimes->{
3.682755604306319*^9},ExpressionUUID->"af24904c-6424-42c4-9628-\
733a98ffc205"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"ListLinePlot", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"R", "[",
RowBox[{"i", ",", "l"}], "]"}], ",",
RowBox[{"{",
RowBox[{"l", ",", "1", ",",
RowBox[{"lmax", "[",
RowBox[{"[",
RowBox[{"n", "[", "i", "]"}], "]"}], "]"}]}], "}"}]}], "]"}],
",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "10"}], "}"}], ",", "All"}], "}"}]}], ",",
RowBox[{"PlotLabel", "\[Rule]",
RowBox[{"species", "[",
RowBox[{"[", "i", "]"}], "]"}]}], ",",
RowBox[{"LabelStyle", "\[Rule]",
RowBox[{"Directive", "[",
RowBox[{"Black", ",", "8"}], "]"}]}], ",",
RowBox[{"ImageSize", "\[Rule]", "Small"}], ",",
RowBox[{"PlotLegends", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\<s\>\"", ",", "\"\<p\>\"", ",", "\"\<d}\>\""}],
"}"}]}]}], "]"}], "}"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",",
RowBox[{"Length", "[", "species", "]"}]}], "}"}]}], "]"}], "//",
"TableForm"}]], "Input",
CellChangeTimes->{{3.648385042597542*^9, 3.648385046847447*^9}, {
3.6483851646559267`*^9, 3.648385175847562*^9}, {3.648385230030488*^9,
3.6483856063439913`*^9}, {3.648385639672269*^9, 3.648385642102378*^9}, {
3.648385927478169*^9, 3.64838594127734*^9}, {3.664868483309063*^9,
3.6648685534148207`*^9}, {3.6650572902898827`*^9,
3.665057300982132*^9}},ExpressionUUID->"e9f893e0-43ea-4242-9950-\
a337cd286aed"],
Cell[BoxData[
TagBox[GridBox[{
{
TemplateBox[{GraphicsBox[{{}, {{{}, {}, {
Hue[0.67, 0.6, 0.6],
Directive[
PointSize[0.008333333333333333],
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1.6]],
LineBox[CompressedData["
1:eJwUl3dYSH8bxlMSipLGj6JFQmnvnLu9tUtR2iLRliIaUmnQHkh7713KaU/a
g0qSkRSSiApv71/nuq/7XPd5zvUdz/PhsnE2sKekoKAYoKag+P9TO/H1aL5c
FRF4Wdql+sEBGKw5jYbeqiIScmZons0dgMfkaXvZ1ipCaZsbW5bYQUTdymZv
oakmGNOo0pQCDsLtw51kzlPVxK8Om+i8gYMo8L5wxyymmtjZu7z4goMDIrQs
950nqok9csXVw84cYOMdcLbnriGIJ29+P2jiwLcyhw8nL9UQcf4H3u7bwYkt
n0wb1ytriI+9pok6IpyoyvvQkfavhmjj6zn6+SwnondLTIto1RL5t7VyfgZx
4oGr41xZQi3xn5gLj1MpJy7XNvdzva8lni6GVehMcmLsX2lEkGgdET5Q5/x4
GxesUrdwzQRs6rtv3SxFuJAY+SVabLiOoL7ZNXHvHBc8WRV/Bh2qJ5poqLuE
w7iQuubuMH61nng9uO+cai0Xrr3tpOPvqSeMZdxb+95zYauiE0vwwSeEwDMX
rmd7ubF6O+npgvsTYrLXKROK3OCvLtOw6HlCNI+1Xj/myo1EH6ttr7kaiEPx
Uy8jUrnx/eqTUx4+DcSjIMEVxwFuXN1vG8g50kBcea5D1UrBg96BDOa5E41E
l9AXxQRhHiQ9TPDtvdtIHGJY+PbRhgdVbOoYnGskWPyMDOpieZD7c2B+XeUp
UeLhnrejkwd+DQGsullPicQtw5j8zQMD/m1X+qhJgrtE3VpQ4BDov9xK8XUg
ifvvNsxprA+h8kmlvVUPSdQ6uAfZxR0C918W1ZOsTYTmcoyQes8hsF4U/D6n
0UT0TZXXVPw7hGP9dsdEfJuI2tdJXrnihxG8wTh7oGxTf7o/cMTpMPyiWlmK
3jcRlxgdeEUzDuOcl0Lvm/3NRFfs9ZneicMIiy7c/lSvmXBd6bj4nZEXZeUR
00rBzcTu4++ECrR4Udz8ycbzaTPB+n1f7sptXuwVOVtu9LOZ2HPUTmXwKS8+
FbEvvT3RQnQp+wYo/uZF7TY+1cMXWgjtT3NTGmJHoOeg+n5/egtxIebX6Efn
I/Dp4qXunGohQvVfL3MWHgF9nvu74/+1Enq5ZQ3fPx5B2s75Li2jVkLE36PJ
hpcPghWO//iiWgltKd57znZ8YKjaP9Xe10pkt3Q1smXwofukS/2hXW2EzJdL
7y/M8iG7UOePmnYbMfREMsOU+yj2L7PTi4S3EYbDVxMXbY4iR+qV9ftnbcQp
xjRH3syjMHyh7Gy3u50wEz/fSvPhKALjRmrL9NoJ14e1etF8x3CY/1tnf0w7
EX/sRVvnpWMQ59sQaB5vJ/rlLcbySo7Bie7y7RD2DoJjVEleZuUYjnV9cuKz
6SBcs342+0gfR7SvfVRGbgfRcecMx+Vbx9GYKJS58bWDYBh8ybm34zi4pQ9Y
SUh1EvNzzjFXdvGD70tsiKF/J+E8U6nub8yPyc7KeZPeTuKy8jzLqRR+fG+4
YKrE0kV4fN/ycniOH3/PB5ez2XQRbVyVV1hFBLCT5cv0bHEXcc/Tv2GfrwBC
5R49f7DRRUwothVMdQngldiol4ZWN3Hkw5NjNswncM5adOZLcjdhSSTxF9qc
wOHcfsZ7C90Er2pYXkPpCVRqnTkscLKH2PK8NyXu3wn4c8lIPrvXQwToFe6W
1hVExZECV+e3PcTcj3SKvMeCkL7+ima/VC9hlXnaf25JEJolTJL9Eb0Ef6qp
/nUuIcQGd6pGv+slaG9/O7moKQQnlqaQ83LPiK07GmM7PYRAJS9nphv3jLh2
ter87sdCYMzczaH39RmBsnudzd1CYP5KajlpPidkTrQ8n10Rgk9s6dmc7OeE
yc+waDdOYWz7LjO2dWsfYdCwLuymLYyqd6OnQmz6CP6tqz2z14Rx5YmklXhL
H1FurnCrPUsYCwXe1fQ8/QRzQaA727Aw2HNlqdmC+olPIYb97ylEcDm/i+7s
x34iu1Km8bCgCCrrHj27yDtAPIl96jhlIQLvV/oKymYDROJzh+20ESIwXYgb
exo+QJCfurvLGzb9S49vTzYNED/2h70fWRTBSMS4xsMfA8THf5fi3A6I4mLZ
l2NUxwcJ5yhdlggdUSzzLYkwWg8SoUxKuYf9RPHypIl7b8IgocZ7J/JkuSiE
3A5uFe4fJGYmounH3oniCuU0lQrNEMH57bDnEqsYOnLGEnbIDxFm+9XoIrXE
MPaAdtHPe4igOa1/uPSWGG5Pm4gWVQwRthfW95pWimH06c2c+1+GiD2nrQSD
58UQl6MVKnBsmNhPv32I4BBHjlAK553zw8T5PUtmfsbiaPUManiYMUxQycby
nQoXR1qFcYX7m2HiaPpqenarOJT5SZtdnCNEeGiRcOS6OC4Zvhd0shwhZi7O
ylOLSeCz0tXg6McjxPoee/ldlyWwdrPlVcDMCJH3qiQjJ1sCnD9HisE9Styi
Wn4zOiMBKuuwm712o8Tg7/+uJrFJovJ18Qu+3FFi7/eX099NJJH0iU7/zOIo
8dL1ScqHaEm8un8/0EZ4jEheOSjn3i8J823mTMpeY0R8nabCAzopUATMlmx5
OkZ8T3aXs9KUQsZlStoU6nFCn6qTbA2Rgt5HjrdsOuOE1LqHdWenFBge4q9v
wjhhrFv/xIlGGvwfuY53vhkn/vV/cqpSk8aq1TmFNf4XxN7866wZIdI4mX2L
a7/3CyL25S4byR5pqEZxVB3teEFQVF9bdaeTQcYTqq/HmV4S/xJH7pzRlUGf
pUwvj+1LQuZxUOtitAwyYyV1mSpeEsuBqvb84zKgf78cSrF1gjiTIMrHxi6L
J8rUsQvGE8TBgdFn9dayELHtvPcid4IQ+SXw365cWTifc8nr2ZggiB/k8N6v
smAbE2Nt158kPvOnPB+QkEO61fl9vTmThAK15wuVW3IwtY1Ze/13kjD1j+/1
6JLDG50ceprTU0TbGnWgJeNJsCnQ1iqVThFNq4GvaSxOondXp0XSzlfE3YPJ
Hc65J6FVku5Je/4VIaiaw520chJ3ONVvPWp5RXzIXJgNUiAQlmZPbcg5TVxl
j1uXvEeA5sLAI8Fb00RLZq992SsCaw058xKvp4mekUaOFX7gX8i1OCf514T3
3SJ+Sl8gQM/abCT9NXGlbCP21XMgQ+/BhhTVzOb5OCDaRi8P/d+G6ZniMwSP
fYIZpbQ8Fm+U05ZcmCFuqP7QyLeRx1pdGv+ZhzPE7IePjI3h8jCrEPmXMzBD
0J0ZbZetkceHxZLAZOo3RJD56NUjs/Lgpi5rEpV9QyQdaZHz36WA0PATrZ6u
b4g5YxlhTWkF2LHFp9jkviHIF7kugfYKmMqedtp4/Ybw/Ku1/0S0AhhZHumo
/DdL9IbK6WmRCtj6375L0J8l4ic8JKcWFcAga7O0eHeWeHbkxeqb/Yp4rWq/
V699luDqlW6w0FBE5MstDJco3hK6KXbletcU8Uad86i83FtC0ZGbtiFHEdZd
nbmD194S8heU1lLHFSE+KjbDVb3pj2ZVU2xXQtTHr0dEV94ScQknfaaklJDt
+u41jeg7YusqW4C4oxLesF5zTnd7Rzx28Nq966ESnFXTVagr3hGaDhU2dn1K
mNmy9aXIyjtC78qFdrktyhjdTpErIPGeoP5+yS9BTBkV8sNSv669J9r8Xsx4
XlAGj3PRvriG94TeLLPky4fKSFTomdm25QPx2vLNWMugMjwy7vLrqX4g3jGY
sIrRqIDDLjLDPfwDMabww+DESRV0CPStuA9/IHZ/KZ4rd1fB1POqcUO2OWLp
1CeJ1nwVeA8QdCx2c8Q978gsq1kV2DjTnm0smiPec9A+DN+vii+j/nHav+aI
TwbWd9QNVVFfPxfXofSRiLrAPxIXrgpLvZs4cf8j8Xsp/5dnhyooI1Nu3nn1
kVho7nT8RKEGlsFR9aHj88QoK0/WiqwaWF9SPmK6Pk94PGPjjfdSw7OD7oF6
vfOEecYL54EKNUQGU+wOZf9E+Bbqb89aUkPYVUHT5iufiA8NZt7MJ9ShIHcq
eKP5E/H5T44ej5M6/vsQ2XOSZYGYin7DMZCvjr9XhM6GXFogXIwzbhz4pI4d
BrE3XjcvEGJ5YRS7jmmg0s/2huq+RWLXeWWNVEcNZDQWlTe5LhIUE+G/XhRo
4ILFzgD93kXi9qxSZe1nDUhdqTu9cfgz4XHAjVdBSBMK6u0Fbf6fibfCRn/c
3TVxK/jQluLpz8S6kP2CQY0mHrhMvX8i94Wgp7+Y83pdE6ItPN+/PfhC8Ip7
LbIpaCH6e9Bzw40vBMtxap8dwVqIn912+L3FV+LKmW97859roeGbalVe01fC
9lNEwD8mbdA9G043YFwiZioLomkttHFIpHbks/ISQWmcxTyQpQ1KmTVLvmtL
BJOu4xudr9oYEzKX2Vq4RJw1uDsQKn0K94yeWQbOLBG9sfItAbdPYVLx2lAO
8zfirfqtZNmBU5DOmkl21/pGTAusK1Sz62A72+fOj/7fCCcZyeyVizpQt6Fw
3lH3jeC2uV/3q0YH70m6wvGlb8QLTw33tm26OJabm2Z8dJk4wJHQb2aiC8Nc
Hu8gm2WCerWgvz5bF9+jHthcfLhMXK5t855f1cWfj3viqMaXiXXV5MXvh/Rw
Js5fx4TxO/G9UerfWR09VL5YGbbV+U7sdr//T+CaHjxDH7ieCPtOiOxjXHdM
10NfzBuf2q7vRGDr0Bam53oIbQ7S3EqzQsgaRvMe+6WH7O5SRSbVFSL3e/PN
Yh593Mqle/4paIXYyLzLna6rjxiFqTN3O1aIu6PD/Ntu6CNSa6/9D5ofxIsd
mk8mc/XRFXzF6YTmD+JWjP0XnjF9LKuE/5aM+EGsm7jOv6YywB8x6kdMgz+I
G6lRQwwiBtgXEPS1g/knMUzTO1RjZYAPqh+i9M/+JPJtd7CM3TNAgpqaQE3a
T2Kll7fXiTSAAuVR97WPP4mXBcs8vl8NcJQh+MBB4VXi4Fsp4+2chjg+sfzl
kM8qwbjaFU+vbwj+tOLyPW2rhI2hL09sgCEaJrvk5nb/IvySFaUTKg0xPC/p
mHvmF8G0jWPf/jlDjHD85DTN+UWcGpemObjfCENHD+r9XvlF2I/PamaeMsKp
LpqlSKXfxNGPSSeK/Y2gM2dHyR7zm5jPZKE4WW0EyT2W0RlvfxOFe/UYDBaM
oEBUVfKJrxEvz9+u+MJpjJEtE8HFwWuEW7+8AtNpY3xjYRSTmlwjvlheZGqL
MEaXdsxCl+A6kfdPPuxfmzFSRzI+2QatExbOrjRdG8ZwWa7zoptaJ4zP7f1y
UNwE+/c8fNsuukHE3q/Oob5ign5Z/9CY8A2Ci8Ux1TfHBHF/hPuvfdggvGLa
1ELemMAh9+hhb4U/BPeDRYYj7Kexg7p2W+KjP0Tpf2luZ0+fxs7LPGsv1v4Q
1M8V5nhjToNGr/c+zP4SIlPji3f7T6PusJZDf+1fosXpIHUonSm29D06HL7v
H9FkGlbFrWkKSyuFW97X/xEGS6eDzUJMweV0VvTB9D9CtS+URrDTFOdHr9zx
4aKA1NP1jsxtZvgebeyqYkSBH5XnzzSommHL/LRWaTAFAiYcz3gFm+Gzbvau
xicUcM1OtB3vMsMopXfl5SUKyOz2YZneeQak30XltsNbMND+gDXi1BkYa73t
bDm7BbbR8fQf7p9Bydunupeit+BKyX8Fn0bO4Ng2k+9k9xYoqX5NfrjvLAYP
Jw60bKHEh1X62D/nzsKBj4PuqgwlTnun69JnnUWGrfLgpDslAhV3hg0vnEXs
lLjBzyJK5Oou0WiLmuNVN21X50dKuP7QDve5bo4Gs2k3g0NUyO50nbVqM8f9
hsmkRCsq3I/e+5JitwX6Pc6FPXxEBWkHIyVzUwsofaZ6ZDVJBbX65788Mixg
ZtJ6fG7fVvTGJY9qf7WAZ+/WhyJmW9Hzyrr+vew5CDbXmqkkbcWTHKtwhdBz
UG90fsE1sRUvMz3EbcbPwcrMLaaHjRoXHvNmavBaQirq5S7Fc9QY3FnV+9vT
ErctcobuplFjLPtgpkeHJcQevvbIeU+NvmQW5mlaKwiqpd1LPrYNetxNgvNi
VnBaQ8d5l22w2qqdZX7OCi49Fm27a7YhgtflunSIFXicqIm4v9sQWG9R71du
BZqStq8bqjSw6epwEn1lhRbNzhT1+zQoann72GC7NfZtWf7v+gQNoni8bV+L
WsOS+a9YwuHtqNz7qnvC0ho6++OGHrtux+n9f+ZUwq3hy3p5IuHpdmwbtx3l
qrOGC409EUC3A4FtRKHzB2ucK0n7bX12B0Sp7UKFmGwwPJjPIFOwA4NSpuH2
ijYQawyJ2rWxA15njcbpXW1gITkb8/rUTqwUad0TSrWBZ+d+6YrUnbiY822y
q98GHJauhWErO8G8RM69/GcDGTsOHicNWpy233htJ2SLs+zXpk4/psUML7ni
aG0L4evbuQ1+0mK28Jfp52hbsBrXqJ3ToYP3nJfi5zZbbF19X3Arhw4L55fm
L63a4tjhqJq6LbvwM5wizPGYHc76u03TW+zClk5Ju08Wdqg7w/04qG4X9gyy
1yxG2eH5nHTgQdbdeLRFrMutww6XFj05pj13Y2sK9WvfdTtkPnu5t2t0NzIk
D6nSC9tDdvcw7YwEPaL3iGrxOdhjTMPi6dEkelRGLhI9j+xhF8H7JXeDHjq5
fJ4/R+wR8d9EZakQA146vVUupjuP4wxtQQW2DGAOVtr+Xfk8KLR5XqkmMIDv
eSp7h+95vBK4Unm7lwFT43dn+WrOo7DMgNaBYg+K0vob2b+dB+uWnrmv4nuQ
KxDAnHHcAYPK/CZcTntALbHFqPG8A1QrlS/+Sd8D9dk/H5zSHeDxzl7mzsQe
9N4OQMW0A56Y3JxvZGTEY8vfg1FsF2CnLvAoW4sRCS6iLDvNLqC8dTxAMYgR
RsrOwQcSLkD9e2NfDMkIiQHLzOejF1Bdy9DwcI0RVIaiX9mYL+LPHZoIS4m9
4JJL/rzD+CLm/X/enXbbC6F9fXxJcRex663dNvbSvXgXmijQO3YRGlk8p/d/
2dTSZ8LT/nNEiZfd1CQ/EyJn9rw8cNYReb8GNmycmGCgalugkOIILZl03vJC
JjzUTUnZPeuI4vu1/b2fmfCVrkP/Du8loKDCs1yQGXQDCQ75ly7hx++dKRfc
mPFTJfrOjbJLMLbkHVuqYsb1uA3Tv6uXIC3kVXFqjRnHAz6HC8MJOgFRw/5g
gX6MwQhrsBOUe1Sm7gWxYJRYfVfe74SsnvrrPs9YUMMafm3rvsswGwm+oMzE
irp/Rdq7bS/jmVm/9ZI5Kwjzz9xDRZt61EzQP5sVxadGm/R/XwYTTVna7yVW
XJtf/hqpcgVJ/rLR5nL/gfrqlG9Y9BXkxZu+LQz5D1G5UnIaM1fwLvaG3+fR
/+D+rp2i54Qz+K1PneM6tA+2162z6G86I+Mi6aPlvg/5X3ZvZ+9zxrPvIq+u
tO6DUuKtoGv0LljpiIqPZNqP8WUtTXUpF9CIXW4uOr8fxwwb1W5au+DfdfaA
/rr9EOVQ9WQLc8FU3OS21V1sUPzp+4yzygWvjT8FH7FlQ7HuYY17r11w/meB
ml0dG2z/OzrvsNMVd7d2RpcysIPQzS4sEHfF21KVJvqL7Gh5/DnurLUr9O5U
Cwe1sINt/U7etQhX5AafMmQ6cAC0guJfqOpdEVVF/bDp2gFE+fx2+vvBFUf/
tkREjB7ALw5W4YtMbrBJ8Ki6KXoQO1ZY1FUU3eDVLlmQEH0QPJ96mmNc3LCf
/jrj1PJBdD/ULzV+7Ab1Ga5kXSMOxA3t4wvvc4Pto9yNH9Uc+Pj3PaT+uuF3
zCG/i7s4cf5pP8PZE+7gXbgn+wOciNTRqv5xzh2TzrkcnO6c0KE4rrPjvjt4
OBK4FrM5wZmV/iu+2R17t5bInZ3kxNmfmi8ef3dHT6qfmzMDF1R+HmI8zOsB
/9jPPcdUuZCm+ej5ETMPhN6SM7p3gws89++J5YV74ONhAZ7HFVyIy9Zwz2/y
gPkgnaH5AheYfkTW8v/wgN/3HupuHm5QZj8QFz7miW07dtvNmXNDbeLIwTpL
T8Q9VX9UF8+NU8OOue1xnrA07Z9TGOSG2d6Du0yfeWLR/LNfAB0P0tz9gl2p
rsLg9nqhrzoP2gXNjPfKXsXlRc8a8SAe/Hq/Vk24X8X7aMbvhS08uBam+nWl
4Cr4WV7Uvac4hCyjYiXJ91fxqVHs8iwO4eR9qY1tHF6QesZ8NevWIez3DLC1
N/PCJ9Z/ysJNh/DiTUuHYawXgr9xG4RsOYzBXc89h/u9AEOfA0VKh1FzbG/r
e9prOGLx6kT6ncPAz1yqcPVruHmQncWx5zAYZVXvdQZdA0sR6yUael64yWs0
JbVeQ2HY5PMbRrwYPlL3nYrSG6JWNR+eJfMiXVI9lEbBG4LrpufX3vDiDqPJ
2yw/byyWfqDefuwInCa+nZ9p8kYv9f7gX25H4EBMXa/d4gPb4jOVfQ1H0OMW
Gi6o5IPn715aRdDwwZ+RZ5takA8+7i2+IGnIh48fQ00pujb9l19K+x/zQSnj
9E9z2uuwtxk8YvZ5078wY3RO9zrM2drbxmSPQsZMbJE69jpsW+6HaIUdRcQX
Tknjl9fxo8/Zt37yKH4oi1Sc4riBJ1JCNbwCx8CmXDmzZH8D9Q/0DGNuHcOP
/ApF+aIbMHpy0Jtq+BiCN/6eO/njBlztx8x8jhxHsPaRp/MnffGmuklo/cZx
jH3yW1AN9oV4i6FS6PBxtCzyeJoM+kI62eDNkeP8mI9eXNjHfhOOzZ/tJgL5
kXb3+ttoh5u4/O+tSOYrfhTvkFhuqbiJ7SkZ9SFSAuBtVSkroLiFwt1OV8Ni
BSBRmdigo3MLf347L5R8E0B650x50cNbSKlfzfutewIpOjHiPQu3EKh8+5RT
yQm84/VYGTrgBwmf5Ru0DIK42rcozaXqB9NvKZ8mXAXBsMjQNHDFDx1H7hpN
jwjC+/4Bi9kEP5jEES3h7EJgHi/+aNLsh3hTVj51XSGcZSJlBRf88MK4Ujsj
QAgsn7KEnJj9MXiI9UhqtRC65Vkzdsv7w0zWtVJ+QQi6h7ldWC/5o+2C5K87
nMKgvSt4NTjeHzcixTc8TITR9IkizrLFH03x6QN0EcLQHyouSv3iD4l8sXC9
NmGEKuSka7EF4HiWuLbihjB2/iw4ba0eAInnCWLvxEQQeX+gbNYzAMxdO12k
roggJOlawkBGAO6MMAgp5IpAxdZ46fhQALapMhVRvBVBYpdt1g+KQDAp+/J6
HxQFc+d86gmhQDD3dc4WnhHF8nf3vnHLQFzcmiyYkiCKgJ6h/cv3AtFTM2Gt
OyoKlWFWj5tNgYhH12jLXjFEpkd3+34LxL5E+q9rBmJIK3pO+437NpK63/63
Fi2Gkz/u8r80ug22JOaOlmExXL3ddFAy+DaKA4+7nmYWR0tL7Bht/W1oUfs8
rD8tjo8JZzRMP98GD7/c0Jdkcfy1irh8gCsIWt6FmT+nxfG28pmMuXEQWugP
Fo/wSEBwdTyP+W4Q6kRYosIuSgC1oZWaZBAkmc7/OlAqAZU3Q8YbK0F4v6Kc
ELkqgfL8qtDjx+9AMN7x6AwkscNc/tRL6zvI+nfUnTlUEsnDwSnbku6gIDFc
QXhYEjdH0m5UDtzB9YTLjhIHpTC47/Hc2+3BuNUR9+yIoxRyBIJm4hSCcTJt
XHd7jRQ6t6hf7PUJRkr38fcTW6XhpdF6O6gyGKtjz++kG0rjV1SFRPuXYBTd
MpOyyZDGm/9KAsOOhqDh4MdtB1akwXXuoOtLuxCIxZ/YMaYigwAXMerS1BA0
3ntgeT9RBn7LH1QZXoXAm/2qmO6iDK6MvT+5vi8USlUbT5gVZHGpP3rZ9XQo
gh7eFpuLl8Vbz1+XPeNC8TM96UPHZ1l0Hmiu2zYSCk29g4drVeRgai4ww7v3
LiRcnU49TZHDO9mdi2MGd/GB9mrn1C85BLPlfGSNuQvlhg/UrEYnMS7V++bT
8F2oeKvccSk9ieQ3+W/1mMOQ0MYyuLiLQITW1Jrm6TDMqPzzvudEYI1eW2wi
KQzbhj9wmD8jwK0h/OjfVBiaPrBdNhEAHHS3KDVwhmP4jhdrwD2gfq+t6B77
cJg+8J+ZXQa8mLlv/MsLh8xUj+bT4/Lg+sZ5PORrOL4IbA3JOyePxCAtnULx
CIzrZu8Xj5bfnOdDvjvdiMBP4Yh1mw559IUX8PW2RoD1Ws9R8XV5KDrYrXfS
RuIaTV1+kZACWumEve2MIkEzxHC377wCBLu1S9IeReLL67WXSY8U8Hb6cb7f
XCRkEw9X7Bnd1EHFXpQi92Cn8VJMapcidC8Hc/P53oPnMl0QvaoiBhvKqpe7
7qH81PWBpFuKcNeIeFy28z48DtRqvqxTRPpfCtuTwvfB3aByYmBFEUSl+7ED
pvdBs5Tf4CekhB0tMj76t+5jjiqU+OykhPFeibGJ7PvgPVxO9V++Ejx5upvr
+u7Daru8LdVHJdDPdJFfft6H+mj+kzJeZeTneRn4cESh4cyM2+HzyjgqeXXd
XD0Kr8PO/rTKVgaVb61sgmsURNTknpyfU4ZSmHfM0QdRCLlgSUgfVQHdPsEA
pvYotPF6bZ28pIIPT27cMf0aBcOnPMk6JSpg7+Y9+HNfNPQu0sXGf1eB/+ut
HZ+Vo6GS/SuySkoVhawrv+ESjZC+ecHim6rY+bJL7PODaPD+ETsU2K6KSn35
pR+d0bBSe8QluUsNjve+ZZ5ZiUaV2pUPz4zVkCNJ9LBxxSBoI15a7bEapqYU
x6R1YqC53vYpd14Nf750nKy6HoMbv28NrYipQ1OF3yMyLwY8hSpDAv7qOI1z
C+3jMRD/lzVi0qcOvetNe823xaIjy/K5G7sG7DkOmOuLx+L4Vt+625c0IFoZ
LZxtF4tjnq6P7z3RgECk/78zsbEwWpC4F0+niZzWa4IubbHI/3op/eE5TcQV
cR/9sBKLur+uG5llmohcEjJoPBwHbZXwnnJqLazRqu/+aRKHhca/St1ntLDs
9Onp/ZA4GJcK5C+UamGa0rY1oj4O3LGN6ge2a2OfYcm1z4txGDld4WltpQ3n
HRelSzniEbgrLupJnTaOUJVYjBnEw+3vXRpBplNYHNwla3UnfnN9/E+Rzqfg
VOlwSKc+Hrn7Nd45PzuFv1umLqZ+iYe0ScNp5WM6mN5WbHqaJwGOv6cJ+VAd
+Oscl3U3TYDRN0o5h3kdeJpFma5EJEDq/sJGg6Yu9jUsr0y0JiA+7puMSpEu
ajLY9I+sJeDVcMwr/+16kCM9818KJYIhFUf+yuhhV+iCzLJDIvZsW72w74oe
2L72SLk+TkT6G42E0TQ9PDbJGTMZT4SzzFKh2JgeWn/JKGXQJ2E9+1W+FK0+
zsWp5ZmoJ2F60Sjxrbw+2LyvCrr6JyF4ej1A0ksfAiOSv5brkzDg9fWmRLE+
6s6lcU2tJCFX+kX27Ht9SKtrt/MLJoOx4y8bcdAAE1dEfr27mIyRBqttWqcN
wDp1t2drVjIoJ9pjaaMM8CAwSy18JhmzJqMb/r0GuKXQ4Xmd/QGUI49fL91m
CBbmf8ZDpg/w7Z2zYZKiITjjdry7H/cAe1YdB4hbhmhP5mCpGXoAS9oDyhkN
hqDjHJ5SZngI0ecjNN1rhjiu+55bRuchFk+p3iiXNsJVkv15fPhDLOmk77Dx
NgL7q/31Z3ofYnKkiGOizgjt14RagnY+wgHqCfED60aYvC5cw6r5CAk65r0C
J40x9lTNjuHuI6hMqzrs9jOG9d+2COeeR2C2TktsbjXG6cXvS2K0KfDdG/BT
Y7sJls0+Kllpp4CH8whn/ikT9HjpiH+LSMFpHYbGdzEm4Dwq7TTXnwLbLft8
NyZMoOrwNViR8TEmQ6P/W+E+DYfyWg5K48cI7j6n13fpNJ6G57YcTnyMG8dp
v0RUbZ6uPxdWCycfY31poVWcwhSfaSlTeulTITwiWdSltakVxyiXhVMhRnko
UiPRFOvBilYaRqmYc6IwevLeFMlRFQ6DV1Nx4KrRD04xM3DzjbSHJqWC5se6
q2+gGaitnbndGlLxL155eHjYDNeMDDjvvk5F00Na/iOHz+CgtoDFIGUarvyY
TPTxOoP0ix6uarxp0Kilkx/sOYPFViwtaKThzN/X+ic4zuLN72MOTy+nQfFA
+vY4j7O4Li127mlUGryUu3O3956FdH3y6YWqNPzyZfYP5zZHR15Dv/JE2mYf
KvvEc90cs6Ercj1/0yB/Vl9jaMQcOX896a/zpEPt42naB4IW4KM7Vmaing52
vpU8/zALlLTcDTK/nI6m8yp5oR8t0Gx6cz4sOn2T5HwTq1XPQapcWPRdTTpa
vCgZ6bLPocNwuczuVTqkTsU+idhmCcUWmmI6qgxUHb+iiguW0D0VfmGSLwMC
Ly85H+q1xCdlFbYBnQwkbk+4lcpqBaskn1/zHhlg7OGw19a0woP3z4/xPchA
hvGjFwm+VnB/yzAa2pwB7SivgdtlVviswca4+2MGTlkoazB/sMJt5tSVit2Z
uLhjUkWLzRruMpTBPuKZiKO50M+nbw2TZr5eK/NMRE75/SgJtsYlHYN6+8BM
qOtK9c0+tca7vR+VQ/Iz8XFc5GrrT2tMCOoZdA5mgl73CYORoA1yrpJNXL8z
cXZ6vDP+gg0c83fax3Nm4amI5bOIdBuUcHrtPKSehdRvu0xkX9nAZ/yuW7dz
FpJcH2ak/WeLgbVJo6CELNysGVlrMbKF6dTrcyZkFp7U5uamRdki9lydqOxc
Fprfqm+R77fFvdFtTiL02TD/KmOVussOH8+VxMlIZuPrNxe+dm07uFaUiBpY
ZsOw2aqtNNwOhR4fGr2Cs2HxkC3b/rkd2n69Ts0rycb0wfNqH3bbozUh7eSH
8Www+w8cFde3x48bHvPHKXJwjad/zCjWHkP/Sj5f48tBq8Pwa7UX9qgODgzu
1cuB2Qu7aYYD56H5l5zj8s5BRYqMR7nNefiPpIr4pOXgmKaxlmDeeWinxN4a
6s5B09x/LHeXzqNDg/ol33IO5Jwc7rVLOeDUTw656/tzoe36xOedvwN+uTxM
7FTMxYqqUuZ8rwOCz94eoL2UC/4zGzPjLBdQSrVWoxKTi65vfbvLbS6gcukv
k/uTXNQp5jNdL7mAZgGZx1Fvc1E9lTor/mfTP/30XwptHq64ULu917oIT8/+
gUTRPMgtBbeHP7gIzfTcRzfO5qGZznD+xOJF7IlWmFULzMP+3ce+9510RNLq
Xtpf+XnQoQ2mdLnvCHJ0xuPuUB7q/Orkmd85wrXN02njdx5GVKnfNEtdAruo
2wktrnyEf2rm9I68hBvldHPu6vngveCjRry/BCW9F2Xuzvk4qbkrh+mkExLo
7yapxeeDFHj+4F+cE/JDnyYsNuTDpDAtnHLJCXEtZKzl2828Yr4FTq3LyCW0
HB7tKEAR7wvRczmX4W9vOpIuWABGCc4vjVuvYL/860gn4wIMNEYUKtlegXrq
Y8lVnwIomH+m+95yBQebPpojtQAdL67U9vM4I/eCtY9CewEY/nXZvgpyxnzJ
wK/f8wXAgkU/x7wz0s/NCdjvLoTwuxLb57wuOEi/0/yWSCH0t6zvyTNzwSoL
e7mySSHUvcreHYpwgZaOllaZdyGufmcaUWx2wVDn1Inmh4XQ+ZY8vfWnCx5G
rOm6kIWgUKTYcfm4Ky57BUTUvylE7qD5OV9rV/AKjZU9oCpCY77kkmSiKwz9
Ct3oDxfhjtSDsZR+V6j+eOuyR7UIZt1HhSpp3IANE6lH54twWV+E20feDal3
y1B8pwi7ndtfLHu74Z6T+kdkF2GjtOcOZ6Ubrhst+55qL0IXB8t56q9uGOqY
lRh6WwRK2ps1GcfcQff3jmj3lmLQsX1r+XPeHbLv7G4KcBTjSm74OHOmO3jl
N2jXZIvha1ah9eWNO4bE8uuPmhZDgeaQbyinB9JE7xnUuBfD6RPD4LylB6b8
vYIeRRZDfY9k5Z5UD6S7d3wZyy1GJMfPZKo3HpC7Zbf3XEsxNkwvbGnh9sQb
7VJrgcliZK/6XNK198T3OTdhfC/Gm75t9sW5npiWHG+N2FmCbtl5/zeLnuib
s1Pay1UCDbsD5xaEryLr/d/WbokSZLWPPurzuoqbae8VC7RKNvtn/OeIp1fh
KdWYXWZZAuULpyj5t3lBvjqBHHMrgUiXS1KBjhfykqnUmIJKULDK7E2f6IWs
gYMztvElEMze7XJm1gvBqsWsZHYJBr5QmdwTuAb30zJ2HNUlUBSMpi3zvoan
kkK8N9tKkCRh79PScQ2L8dtzRwZL8OSZUVwHkzd0C9wPsk+XoDbG37TJ1hsW
U3dzT82XYF3Jrbu8whuyitv1bL6XIFVZeDV1qw9mhG/vMNwowUKq41qEiQ/O
B5ufoaMsBbv53R838nxwoO/8bhPqUpCHJplc//jg+lbDDcrtpci7kBV52eA6
Dv9+5Em/sxQDPToZnrnXMfr8b7MXbSlEN7wfhf+7jhHtZA4pulJQ/+Ydqjp9
A5L/ybRqb+rLBcOpK6U3oDEc/KJs8/2vKw/dtWl9oe+gHWC7mfeNw4psOu8L
i6CKCYvN7wXXT1EZt/qCCI75/nCznvHJsV46rpsIHqobY9+s99KodNHCrZsI
bfaiq1gvwVN9WY0fr28igCv596NvJbB8c1NbQOEWdvWImC+/L8HWX9LOCRm3
8DxXc3fuixL4CBTWP6Dyg2OG/avq7hLs8ddL2C7hB2er02EsdSWoqb4sy3TR
D+3MPP1PN9enpmYtuemhH87kfTDOiS7BoxqmUvpBP5i5DrZ23iiBWRTjdapt
/qCZ5a1itS+BqXUiZYqsP6q8TpGRm/vFzKLA8I2rPxaOf3I9JFSCuk/t/v25
/vjTtt10gLEE1CsL2Y4z/virq/87dHN/Rtwk31b+FwCSSWtSfagY/X5xlqX6
AdhjVBSwtaQY3mt31S3DAiDN+dSjMrQY4zrJw63tAQh7TErp2RQjYLZUboYi
ENt3OJn1SxWD/U9bR7VcIPjy/c4f3VWM718Xnmh4B+KgiulP49dFcAhwufyg
OhAPmeNytEuK8C5zXrR4JRCT3RPc1L5FsLyk73xb9DYKGxLZvdSLcESPPZjT
/TZ4G96ZJO8pwgqp9+l2xW0E5pzPOP+iEMeXUugrV26DScD4/fDmfbWDwtuz
VCIIe5YC906cKwTl+ePpN7yD0Ho2UMb1YCEm8n7wsDUGYcNh6GbUZAEKnHii
Irfcge19bg7R+AI8Jw+HvFS9g/cni6/q6BTgUm9hyd/wO2ifTvs3SVWAz+3N
1RQjdyA5l2A8Up0PnpNvjGbZgrHj2xEpIYd8KJUza2bbBaNOIDTiI1M+/EIZ
LXWLg5FzRt3+J5mHDh1d+5lfwWCz3bNucCEP67o8kubKIcgbD69b3ZUHw8g3
Ne33Q8B3Tv/kdOlmf6T5McwxHQJjqfnTFPq52GX2x+8yfyjs1hYYTb/k4Fn9
QmPZ9VBsDeKveB2Sg3ZTzpivvaHo01t/EsWZA8dtIYzHDtxFKaPv8oWqbFzU
qz9t53wXVvXM1lZq2ZD9zBqS1nIXu9vnzriPZUE/Yn7oHUsYKMcdbj2yycIV
uluuwk5h4A/e+nRiIRN7qXVz7raE4QpmzY64ZSK8k6Xl+75wpNMWfL71IwNu
VWv7XNzC4cjn5vzGMwMEVc4eqmfh2LOvQEzjezqqZlxYSnkjULdF7Gfl5jwq
ufNaqk9ABCx+uY1wv0+DQGuXyPnXEWBaCTobZZYGGrmbZ6+ejITBz9vB/3pS
scfufUXxw0ic/TPa5SSdCr6jqTOMfyJxPujQnnKNx/Dm7Y7KPncPuR6/zEL3
p8CtlTXwcvM9jLvTylm/fwhVG8qnVYz30SjfGyOX/wBJOatHDFTuwzPimiHH
Jp8xCHN9jbh2H81zYeG7uZOwYX1qx6XC+9BaFfNmGE3Azcqrbgszm77ichB/
QDxW+cb597BEgTIEZy8djUP88rr8tFYUft8znRvricEn65eNRgFRyFXYfcHn
fDS4d7+uvVoXBTKqwMNm4z4aTtLpqH6LQkV4P73MvntIvLc1tf1oNN42dux1
PxGOH+8OvfhpEw2KWLPKKKlQWEck8k49jMatkZiKesk70K660uI1Ho1HEXc+
MfMHgvm/0vlexhhcFz9r0ce8ee/89HgzoRODJ5p1KYc+eeOZ/vqP7LAY9IS8
KuKNdQfHlUYfke4YyJRXte8zcsSO5JCeWzSx2HLVpIaz7AzYEasWpRqLQzGF
dx/orRP3ZJLOXrwTC/uKJHpHrjNk7K15r52dsfjtrPbZ7oYjuUYVwuy1PQ72
0irfc6k9yJO/5KvKNOMw65ewfmTah+zUanndEBEHSn+qqZC3fqTvzcrDDwbj
cL3y7dEZy9ukVtMwhy5LPORs3i+2vA0mf1nucJ48Gw+Duzo7xq6Eke9DgxoU
0+Nxq++hcD3FPZLC9GVd8Hw8ShnkFoK0o8ha0T1cRcIJaIp6GFbGFkOu3+2v
rPJJgFlV8Mu8z7HkuZcO8pltCYjXcPSVaIknz21bS7hOnwj+sWdv2JMSyef7
IhOJs4l4k/uXn9MtmbypJsS3lJMIqzZ+2S2nHpKqhrskYn4kwiNzMfHm0RRS
Q+de/THlJLgnciGfKpX8XexTUxeThBqd/+wPfk8lL7fRSim9S8LIUftbq7Np
ZOMdJv0e8WSon/h5kGE4nVxt289lEJKMlm2qGQqtGaRp+dWBV5PJMKA4K21X
kUkmD7tGuQg9gN4zzkmdjCwy5dT51F13HkCd4NYdj84mn8ZPK9VOPQCpncs3
5J9DzodVPXMXe4jpbKUuRpdcsvyTRLxSxEP8kewwv2SRR5YIy287PvcQUStl
7wo18kmrj9YRAoqPsKHFLZQjVkC+1d2So5PyCKoJ4m/5DhaShb4+L6PWH2F4
fP3aP+oisuaNYvWaWQpK9yypri4UkdbSdBPhdSnoD9ziP9lfTOouzS1q7d/k
8h/CxNWyErJlyTFd9sZjXOKaI9TVS8mSd4N1lq8fw4Mrg8W7upS8tsL47xVX
KthDjFiFecrI3RNCPPmbXF5u9F7QMLKMPJbUK8MSkoq2d7SDzatlZFnv3ZOH
Njm8wiry2znLctJ9Re3o2FIq0ihmGenay8lXFNyUJza5+60Hn2o+XwXJsT1q
SMA8DT+LI4T3360gb6j7FIxHpyFNejuL3HwFec9WOVOgJw11xM3RryqVpPrE
woQwZTr2qev/ZkqtJC27WT0+yKTjWLVB3d3VSvLw3IkHpzzSUTTITSWkXUWa
tzoFXCxOx7rqhv/bR1XkNs1VG7n5dBiQn7MDF6vI16ErN7oOZWBdMLDsi2Q1
6U+msjJYZ4A8PVNL6VdNjmTZ+TOnbHJv06mqR23V5Au6RfrJyQyUuUvffURd
QwbVyHHa7s/EK6dChhnFGrLK/zlLqVkmpFPEuZVu1JAzfY8s25I2OVlaKzOz
rIYsUTx1JX0iEzOl3LGzszWkepvoU232LHBU/WOcpK8l/3vuPdR+Lgtz4mtu
TlK1pN5WYQmG9Cz8FZZV9jSvJWuOXL5x4kMWDr4usRm5XkvSKVnI8xzPRssV
ka0X42vJagn7w99csnHI6UgnY2EtORspOZ5ck73ZzyKbyxpqyfWCrzt4/mUj
X/D2fd6uWtL/JL9PuFoOXm4Rv2DQV0uC4djnl/dzUNste4i1v5acF3ottHsy
B122Tqane2pJUTsWdn7eXATcpLb5SdaSxj/uRkq45eJ2n6r3i5JaUmRqq6sQ
mQuKJofW2aRaklXOs4p9s2/aX/xkt+5bS0rI7T39d5MLhZhVY3ds/s9XGe8z
LwrysFPtXeEP0Vpyn0XL86KNPIQVOCrmUW/m++R2+evko3heS4x+sIaMd0h3
Op2Wj4nBxOa98TUkaxrle+Ef+TD+uPNRplENOWbE4sCkWYATXfJbyuhqyP4+
ZRHKTS7LMD3UfvRpNXlMXfD+358FqNZwc9lysZoM3z84QqdbiMR2Mz3+XdXk
S9PUc0K5hVh9ciw8Nb+KvN9XHO9EWYSz6dtlTBSqyN0OvUttFkW4uBbvrjZY
SW7/m7EgV1+EPNtYR5czlaRr8hzDC9ZidLib+vZPVZCSzBd/JF0tRgU1n5a9
aQXJKKAofWesGM94DbvEnpeTvdtXg1MlS/CfAr2Bkmw5qXM0OXExqQSDD/41
J6aXkRvPlgSu/CnBj04Lf1CVkbWzWmFMwqVQ313xQNmilKywio1nsCvF6Bfe
dGfrEvI3tRB/UUIpLNjjdE+vF5F61mkqi72leLBneM0xpJBcak6f7acow4sT
2icmtxeQKb3SO85IlMHJ4pvN9I080k29tz3cqQzGoYMHyt7lkFTSgscdM8oA
mU+08fLZZCNPu8HKxKb/kZZ97H4meSIzxFxgbzleHyyKbxpNJxlkeZ1Ztcsx
ThcsWkCfRk5YMtRWBpXjaEKB2TvfFDJvJPYSTVM5Nsb4U5c6k0njk8tte9bL
Iekn/Z/O33jSe+3ytjGJChQvVP2a4oohTTQcvU3cK3D81Ugbs0MkeTbvxpW4
0goway6NZF24Q4qF3D4Z96UC1KfVvH/G3iBPireomApUgkc8lGa77AVS9sCX
zzNOlUhmm/i6qmAAaV29TOGiSsw0LFtpu3jCR+vFtPaXSry9aHLLiwiE+QFx
NnGhKrT73ox6+CIMPM2fpxfdqmDgQB25rBCNhDHvmqvVVTDfMp3WmhIPu7Qz
MqNrVcDpQhq9LQ8gbipGTStfjfClvB/7LjyGX+tc4/471djZ3dHf+zEN//gL
GWmeV8Pw+772fQczMZp6tH6IqQYvlGQ9eA2zYeYf4+5nUYP7lcUU7MG5SLzS
um1vTg22iC7octTlIzTASDniWw3ajOhV9T4WYuvbp9TLcrV4/Po2Tc/eEkTk
uwsoh9bC8vrIY/1HpVjn+519Z6wWKuJqEqtFZdjnS+PQcKgOgW6rrBZPymES
knXpo3sd/nCXDG1vr0DdSdWGXW11EHv08dzNnkp4a3+4KMhcD6vTinryvVVQ
O3rsnp5DPfZ75mwf6KiGi+cNPY/6erhZnPhS1LBZZzv/1KPdTyAclp9hVliL
7I5IxwHbJzjH8EnJIK4Oe/r2EfT1T8C2N/oxi1c9DjeHZFjsaQCDRIjvYf0n
2Hdx7dvTiw2QeTH46g93A6iOMYZJtDbgSdVzgmOxATrc+yY7DzTCiSXv3Y7C
Rmwv5TC76t2Ifrmv3LB5ChobjmCVsUYwGDBsjnQknhgeWJYSewqSvqEnroTE
7WVtEaOYp4ga9Fyc4W5CtrjoQvL3p5iplBE7cqEJ80RPL4MxidJZfyPbnCa8
ccmKrKshISmacdZ8pgn+v3f6PtvVhB9+Sb8mGJvRP3rExUy+CfnNMVl5RDOK
jDvvBbk3QdtseijWthkDIQcHDTZz9MN3Ml7yb0Zoebhg62QTOJylFP4kNONz
RXHhC4Zm7NE8X30wpxki5W56carN0OPLccwqbsbk2vGj6zeaoSbKddhxU0sc
ydGiq2zGu5iKLPXsZty8SQz3LDRDuIW7nyO+GVxJX95IHWqBqF/ewKBvMwKV
LqWaWrTgJcczHWWLzdxMMYI/oQU3ubTCrMSa0ZdR+K1qsAVe1XvC91I1gzqR
/9cqXSt4RTbkdbqb0CrKFr+i3goOCoNc6uAmODyhZC0OagWLZt83EbkmiPnl
POdtbcV93kPyOYMkZlMPMllvaQNvWKXe8D4STHvfa1nKt6FjMIdzVe8pnH78
fMXj14bfIoN7j3o1YvzonExZUxvCQ+v3hYY34IDs6sQOynYsxchMyUc9wc5I
QUJIuR0fHv83HxhUj8RXXou8we3I+eM2du9iHSSpSz2Xetphk51iXXiyFsb/
dbveo++AyoJkyHHKGvD0TPttM+4AVd+BBxdrqlDPo69l+qADW8b2aj23qIRO
+LR/8GwH/ox1FXf+KAc/z9NX8cc6cfPEVsWRW2WIdKrnDHXvBDWV75Nvv0oQ
ZKDNbN3YCcvVV/phm3w75VMdybW9C19/VFLlP85DwTLnxeeGXdg5fpmM182G
uUhVqENqF3IrV0yLN+eOZq/66eXPXago7j3TMp2C4uvllm5y3dBMPX7ZfnPu
zfO7ums+rBtj3aHsbsH3sdBgtGw21Q2neu9AV5MAXOI+wd0l0APSyfpnt8hZ
fAuIrZfw64EHZ/onTw8vcmPv/GzucA9C9L+YvaEKI9kNrz/j4uvFN4mMbG/J
OPKmZm55pm8v/lx+Y9v1+yEpLL02ITbSC3O6/7ZdtEknnR+px40ef4YHxjVe
0prZ5MHsF7J3bz+DqqCjnpRIPtlNc1PUZPoZtggSrMR/xaSHbMOCnPRzpO77
w77rTyl55NDHaYW451jQus7V86mcpOHUjL6w/BzjryUjJEYqSdZA+fhKvT7s
P7brqFBtNSnuZ3PnaGkfGquHG6Ljasm/NzN29DL0A78vN0pdqidTZq4OPXTr
x3+V3/+sSzaQyRIrKemj/aD/aWtbvdZIyv/wZV9iH8BWXdtkxXKS9OE/8aRQ
dwB9ieYOT4ubSMuy6+w7AgfAfuBkeZpvM+kn7SFJUzOAxuaQl7IqLWTqTeMd
eQsDSG7Ss/63tZXMqk+NW+IcxH+PJbdZN7SS7VYnPrw2GUReZDTHC8c20oBJ
kdY7YhDPLfqoXtO3k56vU9ietg2iKkLW+Xt+OymR+0CsbGMQA9toP8fLdpAq
c5I+RuJDeHV7SEG3uYOsTRFizb8yhPdM/mPZsp3kYgzzycrcITyLt+YWyO8k
H/SeO+T5dgj3btYzetF1kQEaAXRfDw6jR6H7CodNFznGtFPu0NlhdMQKS3bm
d5GiAfkUrInDmNxzjoLlQxfJ2Wpa0j06DMX0LLX6vd0kXaN1JsE0gpoVtF4U
7SYF+fyVrhmOIELqWNEb5W7y8Oqv1WsxI8grZ9sxp9ZNymzPE1MaGQFjedJ7
Vblu0oLriMdL5lEcTzmUOMbVTXpKjksRpqO4Rmez4v6ri7Sm3LrL48EoPISY
b1I3dZGS06sBvq9H4flfXpzPtS5ybo1Z3OLQGC6W0Tzu5eoi6Sa/tzM7jiE9
05N77Ukn6czjPpNXOgbNZIE5BrVOcqN+UpL91xgyE27QHGrtIKmJn8HO8uNQ
2Mb/V0+og6Ro443LDh1HsMzeoqeR7f8ruj7jqXD/N4CTnZk9s/kWMpOs+7Yz
UhJCpVAkq1JCCpVskSQiWUkqSkbiRigj68gIOY6DzIyQ/f89+L/O/Xn2fl1P
r+vBB+1ZfyBfT+gFFQXlp9N/1SPHRHqplr19wMRirURGrB657g52qPboA1aE
irfp1l9Q7Lf2ssyyPlAGim9bBtSh9QQmhet0/YCwdiEgLqEW2TncLNE92Q+C
q27PVz6rQQe8m62ocvoBbZ0I1y11hNTW5Dir/vaDz5XcIkxCVejUpCvLTeOf
QMdP0pf99yekr7R6Tu3pTxAdnTAan1GOklx0NJdnfoLMtJ/KaaAUHa+ObK3S
GwAs5k0eVq0fkPJdI9uEJwMg89W8w4hRMRKo4xC7NjcAGsZmO9kd36CpUNKt
i8aDYEPdsmsvUz4qs9zX7P18ELSHh7anp2ShOsUZn/i1QVASRl2Qy5+K+mlm
iC0nh0Ch7S/LS89i0e4jDamyxUPA5cVhy1PVHijdhFMmj+0XyJaNpxeJDAEN
5bG85l6/AAeHe2whVRIQaTjILPj9F8hjlnMz7PnfXxjq85b3wDBgXlIZ8x/O
BYTu/aWGD4fBKGOSCf/P16Clg578cmkYLI/oigxNF4F2ncpdEvJEkNVDC8Kn
PoAkhqUANici8Hiw3eU4UQoUbH3N0hKJ4Dj3caYiYgXo0xfyam4kgte+4yuE
7s+AXLNnJH2DCOSFrz7mr0NgQq7yNZ/yCAgK5ya5WdeCo1yDfYfdRoCy8PGc
vJo6EGFYEEifMQIaBO9PLMnWA2HukVd3foyAya3Ev7ciGsDb8R/3X7KSwIjb
o/tXRhrBq/UYtRBjEmBzn9OnVfkGbl5rIzGHkIC2+HvdywFNoDJPq97iEwmw
GHZ4EUubwWShK5fRMgk8ERvdmzXZAsoO/tqzqjQKgo4LeS9xfAfUL+wYL3iN
giqmlY/C8m3g3Zli44SCUaCg6rh5WbMdVMVxSd3+PQqWPhHfpIp0gPertDsK
smRgz7fxYTK5A2grh+lkXySDvD9cXwmMnWBShucoMY8M4iSmDrFe6wRlv2ay
JibIIENfTNjyRydoiF6urNw3Bk6d57lxQrEL+M/MCzt5jgHiXLn9z5AuUL4c
d6L73RhYhnpsBU1dILElfEHk7xgI3yyTj2IigFTGomTdw+Og9/SZShNIAIsn
hSc07oyDy2DRqMSLADqOKj9kbxwH+fkMS7kJBPC4us6kgXUC2F2+k7z6mgC8
5RVaT9tOAMvBq7URnwmg5aMvue/5BFA3Lm/X+UIA3MYW57SnJkD0sPbSDiIA
KGrPH3XwN3CYOatT9Z4AJD/JzNWH/gb+LC557mkE4L5t2/qn7TfQA7u51gMI
4ArhUx6ryCS4GqV+ytuSAEyZfQLFPSfB9Yzr2l/5CcA5f+3ogc+TwNhT3Z+5
vwuU3XFTPsQ6BbhHnN8Yx3cBNQ0PXeA0BXjMuS5HanWBoPwbKabvp0CdFY/g
5FAncEka9bSnnwYhbczLd250Ar4WxqUrjtPAVGZrxZOuEzx8edYzuXgaZGdr
fCVEdYDKpXSpZsYZgC7HvRik7wA3lIhXOM7PgBdn2q67L7UBaoOaEs9PMyCM
/NaQ7+t3YBAZdIrEPQuu5L/WUYpsBWG2XAQ/31lgR3zAJaPZAlT6vXNlv8+C
KPVti2s/m0DJ5WDH7f1zQFPKZe6pxzewz8o+ZjNyDgwCnb2S040gbNjjp/TU
HMgPMWvION0A/PljloLN/wBvP4/92egLcH8vn8D49g/YNTEW1cpdB86Yan0m
M86DO10x7z6frgFvzhzWcdaaBw4XTLbVlKoAjy3fp9ve8+DIbxURyfgKsCl8
RVY7ax5oK2yfze79CJrmdoc865kHg+zX1oWZ3wOFKnL9C+YFkJ6VlLthUwi6
249uHdVbAGMvw0SvXckBWVn/1LP9F0Anc1KvRG8KSN7FHp79duF/e3q4oyMS
Chw9jOmsxhdACEvvmWd1d9D9jL8/3uxdBHGqQnfyM1NQTZSU2Ge7RcDdTQ4t
rMpBdedec95OWARsVLbmNjpvEN+6zuhyyyIIj3r9n9T4e5R9K7xNmmEJLEWS
umpUy9AZA01+DoMlwPOGN5cqvBJxsCoxl91ZAhlO4w4XBxA643NjTbhq6X99
PsmlLVCHKhI8tA03lkD6MbIK0bAeTR9oV1PU/AuI409kJnwaET9x1IwY8Bdk
7/5inPb0G8orujVn/+kv6DS6sWpa24wCNwSeZGz8BfqH6BYNxluR92Oq+iKd
ZfAsaW/UEEM7eqV3YjMxZBmQz6v6lyV2IHM/mXaT+mUQWGuksv95J2p1ziZ0
Mq6ATXARruZ3IVUXvleqlitAe+eWfOc7AspL5Vz1SVoB0d/8+z0/dCOz7RCv
yIEV8O6r2ZPQDz9QRU59fZjkKti3//PBpnc96JhHR9VZz1WwrUjLQ5Xfi45O
6UuIla6CcbGjamNpfUipNbS+ieofuGMrpqsZ3Y+Ouv24c9biHzgu25745fpP
JFh1T2Mk5R9om5s4fsRxANm0LgycHP8HnqNvjolag+i4dum5z2pr4MUV53Pe
vENo07mmUfDeGmioA32vp4eQ7Utfvmvda0BLeGy3SOUvpKtbduGb9DoINM5g
fXJvGPUZq7SJ3FwHRssTxH+SRHTM0ML/Rss6qExe+H03jIherDkk/BDdAM47
vrYDQ0RUkF1tq3t9A2x5e8E/B0cQXQo1S0nLBtD23SNYHDWC9I3/7RyW3ATG
No1vuQZH0O/0kWhC0CYQ8PghwryfhM4PPOC692MT2M89/RHjR0KLn87THVPa
AkZvLePiK0mIi6fk56GYLaDcEmrIvENC+eQPfPqTW8BzXkN1TXcUyZXlm18x
2QanO9u/WweNoqBo9j1NedvAziW1grdkFIkLV7MeZdgB/F3WX3UnR9HgU/ol
qks7wK3YeaZJkIz8egeCf7fsgGCPda2yI2T0cTGwkoafCu51FCfSXyMj+4Sg
/iIzKugQm8X56SkZ9RCPRG4FU0GpMws8hM9kJDFiWztdTAXvfqc9YDP4v/y9
ye3gcSp4skLxqeE/MvJg2akrEaKG126sR+RxjKHi2Y8pj6yoYdsstVCo9Bhq
fMjGLhxBDQVShq0G1McQZ+p/2scRNUzWFWGrMhxDNx581FZZpYZ3utQZFI+N
oapzrGr1irsgUaKwXtN2DJW8NLBgvLQLipVIPpmyH0PD0dqlVNm7YNDN5wmm
DmMo3qnu1buhXVCm8EuYi90Y8rdKu8AnQAMvpinOGVqNoe1LzxU1bWggOXuC
bsNkDHm2tZ0UTKSB7Jbt2Y+0xtAtU/d9n9ppIGla2UBEfgwVClExCbPRwgPI
bj1PYAwNcT6x1T9KC7d0pf0g7RjKnSq+pRhLC9OKayfXpsno3orv7Nh3Wqj5
R0+Y3EFG53jzJdzZ6WDuN2tLjg9kRFxZSS+3ooM7h8d1HiaS0YH0kd/9SXQw
5UizyA0fMnLQm0xq66OD2fJH4rtNyajseq7WUxF6+PpQTlq1OBnpZAzd1nah
h5tNcW9MVkfRVFkWY+Urekh8Hhh1t3kUMV65FM23QA9J37vNw9NG0e4kqbqT
hxngzJzAKVePURTQ4OR7PYwBOpyod1Q5NIruPXZwvtXKABkusuvT7RpFsxOS
Xp58jNBRIVZhsZmENIs/eh9xYYQSJ6J92RJJ6CXzHkuuIka4BIzv+9iRUGqO
w0b7FiMsaJykkxMioei7qz6hFkzQ/I9Tju3QCLKVK3i7/xkTvEo1L86YMYJu
8RV9aZthgt28xcY2Z0aQUfJwva/ubnhd7UeftdAIyqVq7ORK2A3V5AiPBPqI
qNaqkauSvBt+OOIv+vkREa2KaFRfOswMK4TJ9FaWROQ7E84mGc8Mwxjj3s0y
ENHLRSuLqTFmmDoXmGkZOoxaVGOa6nRYoJyDZYf0wi/EQ59OepfMAv2tD3tG
Of1C77dCtkvmWWC1xB+n3NYhVAWm436Ys0L63/5/MzSGULuOGwdPPiucO8le
k5E9iE4PBu72p2ODJlalyS2sg0iOtol9x4UN/ohLCDL3H0CxmWrVRV/YYDmD
wiMD4k+0m4UWxUuxw0+tFbx9R36i4hnar8/D2eFIwKkj8sX96HdIUTJ5kh1q
qt5xdxboR+cU/9tSkuWAf/XCC9NC+1AaYc8YqwMHnHN2Vv832YtU7PSXo2I5
YFdOPOuTE71oTyyDTEEtB+z1kBNIruxBntqj171WOKDSwCEhTukeNJ2fMdkj
twc+r9oulIn/gexPnY35fX4PLHx1X3turRvd5y31yUnZAz+ep2YOu9CNLDSq
XrN07IExviry650EJLFX0kmSkRP6JSSEewECylg1eDcFOaGxtGL3wpsuxPVN
55trICck0cvkpIl0IZs3dEOPSjihz407p2/HdaLFS8GygX84ochzF+1Cqk4U
qhL5W0COC7Jf9T2r7teB9Guk3a+7ccG1kL037Lva0Xj+7eW4HC7oelzPRhu1
oVebtBPuJC7YpVrYpvTuO2KUD0+lE+eGpC0xkeNZrWhBzeC62zluWPBUw6A4
pQVFCTCSH2ZyQ+mLpsyXEpuRZ6aB470RbrjefpUuKr4J/XYStTOU5IG+M5F3
pR5+Q/3cKY6ECzyQRT8/0TjpK4LGyX9VX/HA1hXJxY20RmSmzlDqPssDWXur
GSxfNiABiRUtPxVeaDc03GFQWo+ed+Tut7/JC1FsUcfIty+o6GrBAUHEC820
Ntf3Ddeh9wSarQp6PmiV+j5NaK0WBQ4dstY6xgeTr81rfuCtRYo+XquZKXxw
0Mg3eOpQDZqNeFA+S+KDj+6H7vKtqEZzrdBb6gA/bChV9Nq7uwpZ2JmsmQTy
w5tKwpdqzlci89lTBo5f+eGZzi1OxZoK9Oi+EzjPIwAfinuzakuXI7g/bvK0
qwB0ladjK35Yir4cvG1iVSIAaUyFnI7t+og6pi3c9egEoXaBqk17wAfU4FPh
qmQnCH0CHtQMrxQjdgEDJ/ECQZhJr84nHVCEdGSP3OPbFoTPnAOEr1G9RS09
Urt5rYVgHLeztjPfa1TVe0JZ7JUQZPU8liuhno8OyEYaHqYShjTv8/ZNOuYi
v6ODKa6nhGG2rfDwvcgs5EI9HJtbLAz/A4q1HA+eowaVL7E7zCKQkylT0EY2
FWXdCmK94SYCze7ve/yW4zEqoxqKYakXgYb8PKPzrg9R2Zva5EbxvfCn8doe
UkQEer6wv/ZlyF54U7TvLHV8KFruc33zgbgX1mnqKX4Y90efn+1jnNMThRX2
GcG8yBsFZzZEOuSIwtwM/47+WGdEmBcJ2qYRg6i89XdB90mU4fkxzlNdDIbW
3u/jjTBCBydTvAIuiUFTyVGivJYamnZoYZZNF4O+lt9ylKnFkYnMDRf/TjH4
hWlTchcrE0oK3vTzpBeHfve+s50K/lN9e4PNhFFbHG5VPIvvt+yoPhWExi2u
ikOfRelVp6SC6j7Fx27ar8Shx42PUYdpQqrhWuZoH1EcRpmd1U8QM6v22EgI
khKQgA3d+m3yekzVz/yabcVPSEC5OHbOQqmyKrkev+zOaAm4c0CoINjdvMpG
9GnyoUYJSNzO9ShtnqiIClLxtKGWhFT/f9lSRm4KOtg0UX87awOw06uebnGX
YXvL/tTft4wdElJNTa0mRXEib/OLjGvY7XacibQl2KavbU0PLmNzJfAbqx6S
pniWo4SNOhD7yPWbyi+qsFWsaw15d8lQfJrEPuFsgs2QI4UiY7FtHQVuxXZj
I9ugH77CshSzs16MVr+AHceacXHkLXarYaGG3xp2aMrawILhfxT/65tSPZOA
Lcq8JV/xC5uFhfSZUWEfxa6x6STLW9g13ll5ca3Y767MC7aI7KfY8aiKI4Mv
tvTewRjjL9iPNcyGo/nkKP7Pnhj90xM7RnNfz8Ev2B4xK8wvBOUpbjwWel/0
Gvblp1G5Ja3YupW3+l1kFSi+2RPzVOkutnXPoPteIrZdb/oPJd0DFHutqwRc
TscWyTn5r3MLe05BjcbdSZHiPNpMG+U6bDWRpAF+LiWKxRpfBJcaYRt5fQpZ
uomtFVP+s7MQu7M/64H1CLZeBYoM5FWmuHUif/KIBXaYm1NpTRj2RL4THbkC
O+3vx8WSBWxNbq4XB/erUHxJZ+Swuwv2jo79rGU6ts70k8XpXmzr7g8P9LhU
KZY9TFiwOoY9YKwfLBqDbfM0MriwCduE9v7RDQY1ik/kBB5jM8HePKQyOxGO