-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathRayTracingPostProcessingRayPosition.nb
More file actions
16554 lines (15901 loc) · 728 KB
/
RayTracingPostProcessingRayPosition.nb
File metadata and controls
16554 lines (15901 loc) · 728 KB
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
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 0, 0]
NotebookDataLength[ 745556, 16553]
NotebookOptionsPosition[ 684050, 15893]
NotebookOutlinePosition[ 684449, 15909]
CellTagsIndexPosition[ 684406, 15906]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"Ray", " ", "Trace", " ", "Through", " ", "the", " ", "IONOSPHERE", " ",
"Mathematica", " ", "Notebook", " ", "Post", " ", "Processsing"}], " ",
"*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"Copyright", " ",
RowBox[{"(", "C", ")"}], "2018", " ", "Sasan", " ", "Ardalan"}],
" ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"This", " ", "program", " ", "is", " ", "free", " ",
RowBox[{"software", ":",
RowBox[{"you", " ", "can", " ", "redistribute", " ", "it", " ",
RowBox[{"and", "/", "or"}], " ", "modify"}]}]}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"it", " ", "under", " ", "the", " ", "terms", " ", "of", " ", "the", " ",
"GNU", " ", "General", " ", "Public", " ", "License", " ", "as", " ",
"published", " ", "by"}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"the", " ", "Free", " ", "Software", " ", "Foundation"}], ",",
RowBox[{
"either", " ", "version", " ", "3", " ", "of", " ", "the", " ",
"License"}], ",", "or"}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{" ",
RowBox[{"at", " ", "your", " ", "option"}], ")"}], " ", "any", " ",
"later", " ",
RowBox[{"version", ".", "This"}], " ", "program", " ", "is", " ",
"distributed"}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{
"in", " ", "the", " ", "hope", " ", "that", " ", "it", " ", "will", " ",
"be", " ", "useful"}], ",",
RowBox[{"but", " ", "WITHOUT", " ", "ANY", " ", "WARRANTY"}]}],
" ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"without", " ", "even", " ", "the", " ", "implied", " ", "warranty", " ",
"of"}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"MERCHANTABILITY", " ", "or", " ", "FITNESS", " ", "FOR", " ", "A", " ",
"PARTICULAR", " ",
RowBox[{"PURPOSE", ".", "See"}], " ", "the"}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"GNU", " ", "General", " ", "Public", " ", "License", " ", "for", " ",
"more", " ",
RowBox[{"details", ".", "You"}], " ", "should", " ", "have"}],
" ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"received", " ", "a", " ", "copy", " ", "of", " ", "the", " ", "GNU", " ",
"General", " ", "Public", " ", "License"}], " ",
"*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"along", " ", "with", " ", "this", " ",
RowBox[{"program", ".", "If"}], " ", "not"}], ",",
RowBox[{
RowBox[{"see", "<",
RowBox[{"http", ":"}]}], "//",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"www", ".", "gnu", ".", "org"}], "/", "licenses"}], "/"}], ">",
"."}]}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"Author", ":",
RowBox[{"Sasan", " ", "Ardalan"}]}],
" ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"Date", ":",
RowBox[{"December", " ", "28"}]}], ",", "2018"}],
" ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"http", ":"}], "//",
RowBox[{"www", ".", "radiocalc", ".", "com"}]}],
" ", "*)"}],
"\[IndentingNewLine]",
RowBox[{
"(*", " \
", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"The", " ", "Mathematica", " ", "Code", " ", "processes", " ", "the", " ",
"Ray", " ", "Matrix", " ", "computed", " ", "using"}], " ",
"*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"The", " ", "Mathematica", " ", "Notebook", " ", "for", " ", "Ray", " ",
"Tracing", " ", "which", " ", "is", " ", "based", " ", "on"}],
" ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"the", " ", "FORTRAN", " ", "Code", " ",
RowBox[{"by", " ", ":"}]}],
" ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{
RowBox[{"R", ".", " ", "Michael"}], " ", "Jones"}], ",", " ",
RowBox[{"Judith", " ",
RowBox[{"J", ".", " ", "Stephenson"}]}], ",", " ",
RowBox[{
RowBox[{"A", " ", "Versatile", " ", "Three"}], "-", "Dimensional"}]}],
" ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{
"Ray", " ", "Tracing", " ", "Compute", " ", "Program", " ", "for", " ",
"Radio", " ", "Waves", " ", "in", " ", "the", " ", "Ionosphere"}], ","}],
" ", "*)"}], " ", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{
RowBox[{"OT", " ", "Report", " ", "75"}], "-", "75"}], ",", " ",
RowBox[{"US", " ", "Department", " ", "of", " ", "Commerce"}], ","}],
" ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"October", " ", "1975"}],
" ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"Visit", " ",
RowBox[{"http", ":"}]}], "//",
RowBox[{
RowBox[{"www", ".", "radiocalc", ".", "com"}], " ", "for", " ",
"details"}]}], " ", "*)"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"The", " ", "program", " ", "reads", " ", "in", " ", "the", " ", "Ray", " ",
"Matrix", " ",
RowBox[{"rayMatrix", ".", "wdx"}]}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"It", " ", "also", " ", "reads", " ", "in", " ", "and", " ", "uses", " ",
"the", " ", "parameters", " ", "stored", " ", "in", " ",
RowBox[{"W", ".", "wdx"}]}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"Plots", " ", "Ray", " ", "Points", " ", "in", " ", "3", "D", " ", "with",
" ", "Origin", " ", "at", " ", "Earth", " ", "Center"}],
" ", "*)"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"Uses", " ", "Plot", " ", "Range", " ", "to", " ", "Show", " ", "Ray", " ",
"Paths"}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**",
"**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**",
"**", "**", "**", "**", "**", "**", "**", "**"}], "******)"}],
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{
"**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", " ",
"RAY_PATH"}], " ", "Ray", " ", "Vector", " ", "3", "D", " ",
RowBox[{"Plot", " ", "**", "**", "**"}], "*",
RowBox[{"**", "**"}]}], "******)"}], "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**",
"**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**",
"**", "**", "**", "**", "**", "**", "**", "**"}], "******)"}],
"\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"RayMatrix", " ", "=",
RowBox[{"Import", "[", "\"\<rayMatrix.wdx\>\"", "]"}]}], ";"}],
"\[IndentingNewLine]", " ", "\[IndentingNewLine]",
RowBox[{"lengthMatrix", "=",
RowBox[{"Length", "[", "RayMatrix", "]"}]}], "\n",
RowBox[{
RowBox[{"Print", " ", "[",
RowBox[{"\"\<Length Matrix=\>\"", ",", "lengthMatrix"}], "]"}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{"WPost", "=",
RowBox[{"Import", "[", "\"\<W.wdx\>\"", "]"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{"CoordinatesRay", "=",
RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ", "\[IndentingNewLine]", " ",
RowBox[{"*", " ", "Initialize"}], "\[IndentingNewLine]", " ", "*)"}], "\[IndentingNewLine]",
RowBox[{"PIT2", "=",
RowBox[{"2.", "*", "Pi"}]}], "\n",
RowBox[{"PID2", "=",
RowBox[{"Pi", "/", "2."}]}], "\n",
RowBox[{"DEGS", "=",
RowBox[{"180.", "/", "Pi"}]}], "\n",
RowBox[{"RAD", "=",
RowBox[{"Pi", "/", "180."}]}], "\n",
RowBox[{"CC", "=",
RowBox[{"2.997925", "*",
RowBox[{"10", "^", "5"}]}]}], "\n",
RowBox[{"K", "=",
RowBox[{
RowBox[{"2.81785", "E"}], "-",
RowBox[{"15", "*",
RowBox[{
RowBox[{"CC", "^", "2"}], "/", "Pi"}]}]}]}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"TLON", "=",
RowBox[{"WPost", "[",
RowBox[{"[", "5", "]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"TLAT", "=",
RowBox[{"WPost", "[",
RowBox[{"[", "4", "]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"Print", " ", "[",
RowBox[{"\"\<XMTRH=\>\"", ",", "XMTRH", ",", "\"\< PLAT=\>\"", ",",
"PLAT", ",", "\"\< PLON=\>\"", ",", "PLON"}], "]"}], "\[IndentingNewLine]",
RowBox[{"Print", " ", "[",
RowBox[{"\"\< TLAT=\>\"", ",", "TLAT", ",", "\"\< TLON=\>\"", ",",
"TLON"}], "]"}], "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ",
RowBox[{
RowBox[{"EARTHR", "=",
RowBox[{"WPost", "[",
RowBox[{"[", "2", "]"}], "]"}]}], ";"}], " ",
RowBox[{"(*", " ",
RowBox[{"Earth", " ", "Radious", " ", "Km"}], " ", "*)"}],
"\[IndentingNewLine]", " ",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"lengthMatrix", "=", "362"}], ";"}], " ", "*)"}], " ", "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"ii", "=", "1"}], ",",
RowBox[{"ii", "<", " ", "lengthMatrix"}], ",", " ",
RowBox[{"ii", "++"}], ",", "\[IndentingNewLine]", " ",
RowBox[{
RowBox[{"Ray", "=",
RowBox[{"Re", "[",
RowBox[{"RayMatrix", "[",
RowBox[{"[", "ii", "]"}], "]"}], " ", "]"}]}], ";"}]}],
"\[IndentingNewLine]", " ",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"RayPath", "[", "]"}], ";"}], "*)"}], "\[IndentingNewLine]",
"\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]",
"\[IndentingNewLine]", "]"}], " ",
RowBox[{"(*", " ",
RowBox[{"End", " ", "For", " ", "Loop"}], " ", "*)"}],
"\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"Store", " ", "in", " ", "column", " ", "format", " ", "in", " ", "text",
" ", "file"}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"sfile", "=",
RowBox[{"OpenWrite", "[",
RowBox[{"\"\<earth_center_spherical.txt\>\"", ",",
RowBox[{"FormatType", "\[Rule]", "StandardForm"}]}], "]"}]}],
"\[IndentingNewLine]", "\n", "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"ii", "=", "1"}], ",",
RowBox[{"ii", "<", "lengthMatrix"}], ",", " ",
RowBox[{"ii", "++"}], ",", "\[IndentingNewLine]", " ",
"\[IndentingNewLine]", " ",
RowBox[{
RowBox[{"Ray", "=",
RowBox[{"RayMatrix", "[",
RowBox[{"[", "ii", "]"}], "]"}]}], " ", ";", "\[IndentingNewLine]",
RowBox[{"Write", "[",
RowBox[{"sfile", ",",
RowBox[{"Re", "[",
RowBox[{"Ray", "[",
RowBox[{"[", "1", "]"}], "]"}], "]"}], ",", "\"\< \>\"", ",",
RowBox[{"Re", "[",
RowBox[{"Ray", "[",
RowBox[{"[", "2", "]"}], "]"}], "]"}], ",", "\"\< \>\"", ",",
RowBox[{"Re", "[",
RowBox[{"Ray", "[",
RowBox[{"[", "3", "]"}], "]"}], "]"}]}], "]"}], ";"}]}],
"\[IndentingNewLine]", "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Close", "[", "sfile", "]"}], ";"}]}]}]], "Input",
CellChangeTimes->{{3.755460160174726*^9, 3.755460200851099*^9}, {
3.755460867164393*^9, 3.7554608696101847`*^9}, {3.755461109551591*^9,
3.755461168134016*^9}, {3.755461222559847*^9, 3.755461225073848*^9}, {
3.755461266072176*^9, 3.755461289662408*^9}, {3.755461340848538*^9,
3.755461352157701*^9}, {3.7554613900973577`*^9, 3.755461449990507*^9}, {
3.755461589566959*^9, 3.755461619871396*^9}, {3.755463118441763*^9,
3.755463126184739*^9}, 3.75546315757937*^9, 3.7554636478915443`*^9, {
3.755464749687929*^9, 3.7554647897473803`*^9}, {3.755466771693224*^9,
3.755466774660864*^9}, {3.755466831213698*^9, 3.7554668354604*^9}, {
3.755467133843219*^9, 3.755467145920827*^9}, {3.755467324005324*^9,
3.755467334735742*^9}, {3.7554673830966597`*^9, 3.755467418515942*^9}, {
3.7554674620970993`*^9, 3.755467466636093*^9}, {3.755467638108069*^9,
3.755467638523624*^9}, {3.7554677940417347`*^9, 3.755467794707505*^9}, {
3.7554679523070908`*^9, 3.755467955262204*^9}, {3.755468165531446*^9,
3.755468166386797*^9}, {3.755468265416511*^9, 3.7554683358189363`*^9}, {
3.755468422710059*^9, 3.755468463299807*^9}, {3.755468666343663*^9,
3.755468739513462*^9}, 3.755468770988674*^9, {3.755468815015856*^9,
3.755468821191153*^9}, {3.755469030255931*^9, 3.755469033941498*^9},
3.75546930322577*^9, {3.7555121188523073`*^9, 3.7555121247926283`*^9}, {
3.755513114731526*^9, 3.755513121472557*^9}, {3.755555558730844*^9,
3.755555590016902*^9}, {3.7555556365887957`*^9, 3.755555727513669*^9},
3.75556328754053*^9, {3.7555637437099257`*^9, 3.755563992239903*^9}, {
3.7558540251641273`*^9, 3.755854042142034*^9}, {3.7558541294245243`*^9,
3.755854180837818*^9}, {3.755855284898017*^9, 3.755855469860422*^9}, {
3.755855502029087*^9, 3.755855536897591*^9}, {3.75586584972805*^9,
3.7558658655082073`*^9}, {3.7558660745526543`*^9, 3.755866177455607*^9}, {
3.755866292608629*^9, 3.75586630756979*^9}, 3.7558678586088123`*^9},
CellLabel->"In[1]:=",ExpressionUUID->"f3bf48f0-c859-41ce-8d10-9071e7993b36"],
Cell[BoxData["604"], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.967434461133341*^9},
CellLabel->"Out[2]=",ExpressionUUID->"48abf5a1-cfd4-7e4b-9143-17b3f9f4967d"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"Length Matrix=\"\>", "\[InvisibleSpace]", "604"}],
SequenceForm["Length Matrix=", 604],
Editable->False]], "Print",
CellChangeTimes->{3.9674343335148125`*^9, 3.967434374408081*^9,
3.967434461136343*^9},
CellLabel->
"During evaluation of \
In[1]:=",ExpressionUUID->"1389afee-fe4f-8941-95c9-690148c392dd"],
Cell[BoxData["6.283185307179586`"], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.967434461138342*^9},
CellLabel->"Out[6]=",ExpressionUUID->"ee394879-3655-e04b-932e-d13479ca6947"],
Cell[BoxData["1.5707963267948966`"], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.967434461140339*^9},
CellLabel->"Out[7]=",ExpressionUUID->"3f832850-27f3-f34b-b44b-96833184d961"],
Cell[BoxData["57.29577951308232`"], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.9674344611413403`*^9},
CellLabel->"Out[8]=",ExpressionUUID->"80b66d63-64ee-0146-9e24-2a04aee2b5da"],
Cell[BoxData["0.017453292519943295`"], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.967434461143339*^9},
CellLabel->"Out[9]=",ExpressionUUID->"2ff0274d-6a3c-5043-addb-efe60e9c1841"],
Cell[BoxData["299792.5`"], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.9674344611443405`*^9},
CellLabel->"Out[10]=",ExpressionUUID->"ea7c05e5-f779-964e-84c3-a7693a8e26dc"],
Cell[BoxData[
RowBox[{"-", "4.2912410820646`*^11"}]], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.9674344611463394`*^9},
CellLabel->"Out[11]=",ExpressionUUID->"e232362c-1848-ea45-8b9e-9d0f7db54069"],
Cell[BoxData[
RowBox[{"-", "105"}]], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.96743446114834*^9},
CellLabel->"Out[12]=",ExpressionUUID->"d2b3fece-4b8b-2f42-acfc-a6e7a81d5e4c"],
Cell[BoxData["40.`"], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.967434461150343*^9},
CellLabel->"Out[13]=",ExpressionUUID->"d871e84b-78f0-ec41-ad9e-75fd4fc37722"],
Cell[CellGroupData[{
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"XMTRH=\"\>", "\[InvisibleSpace]", "XMTRH",
"\[InvisibleSpace]", "\<\" PLAT=\"\>", "\[InvisibleSpace]", "PLAT",
"\[InvisibleSpace]", "\<\" PLON=\"\>", "\[InvisibleSpace]", "PLON"}],
SequenceForm[
"XMTRH=", $CellContext`XMTRH, " PLAT=", $CellContext`PLAT,
" PLON=", $CellContext`PLON],
Editable->False]], "Print",
CellChangeTimes->{3.9674343335148125`*^9, 3.967434374408081*^9,
3.96743446115234*^9},
CellLabel->
"During evaluation of \
In[1]:=",ExpressionUUID->"513e08ee-b395-e342-b490-3940b4ace7ca"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\" TLAT=\"\>", "\[InvisibleSpace]", "40.`",
"\[InvisibleSpace]", "\<\" TLON=\"\>", "\[InvisibleSpace]",
RowBox[{"-", "105"}]}],
SequenceForm[" TLAT=", 40., " TLON=", -105],
Editable->False]], "Print",
CellChangeTimes->{3.9674343335148125`*^9, 3.967434374408081*^9,
3.967434461154339*^9},
CellLabel->
"During evaluation of \
In[1]:=",ExpressionUUID->"810eaeea-ac47-364b-8162-dfb79af262b4"]
}, Open ]],
Cell[BoxData[
InterpretationBox[
RowBox[{
TagBox["OutputStream",
"SummaryHead"], "[",
DynamicModuleBox[{Typeset`open$$ = False, Typeset`embedState$$ = "Ready"},
TemplateBox[{
PaneSelectorBox[{False -> GridBox[{{
PaneBox[
ButtonBox[
DynamicBox[
FEPrivate`FrontEndResource["FEBitmaps", "SummaryBoxOpener"]],
ButtonFunction :> (Typeset`open$$ = True), Appearance -> None,
BaseStyle -> {}, Evaluator -> Automatic, Method ->
"Preemptive"], Alignment -> {Center, Center}, ImageSize ->
Dynamic[{
Automatic,
3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification])}]],
GraphicsBox[{
Thickness[0.0016806722689075631`], {
FaceForm[{
GrayLevel[0.93],
Opacity[1.]}],
FilledCurveBox[{{{1, 4, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {
1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}}}, {{{
25.499999999999996`, 2.5}, {25.499999999999996`,
1.3953100000000003`}, {24.604699999999998`,
0.49999999999999994`}, {23.5, 0.49999999999999994`}, {2.5,
0.49999999999999994`}, {1.3953100000000003`,
0.49999999999999994`}, {0.49999999999999994`,
1.3953100000000003`}, {0.49999999999999994`, 2.5}, {
0.49999999999999994`, 23.5}, {0.49999999999999994`,
24.604699999999998`}, {1.3953100000000003`,
25.499999999999996`}, {2.5, 25.499999999999996`}, {23.5,
25.499999999999996`}, {24.604699999999998`,
25.499999999999996`}, {25.499999999999996`,
24.604699999999998`}, {25.499999999999996`, 23.5}, {
25.499999999999996`, 2.5}}}]}, {
RGBColor[0.699951, 0.699951, 0.699951],
Opacity[1.],
JoinForm[{ElisionsDump`Miter, 10.}],
JoinedCurveBox[{{{1, 4, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {
1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}}}, {{{
25.499999999999996`, 2.5}, {25.499999999999996`,
1.3953100000000003`}, {24.604699999999998`,
0.49999999999999994`}, {23.5, 0.49999999999999994`}, {2.5,
0.49999999999999994`}, {1.3953100000000003`,
0.49999999999999994`}, {0.49999999999999994`,
1.3953100000000003`}, {0.49999999999999994`, 2.5}, {
0.49999999999999994`, 23.5}, {0.49999999999999994`,
24.604699999999998`}, {1.3953100000000003`,
25.499999999999996`}, {2.5, 25.499999999999996`}, {23.5,
25.499999999999996`}, {24.604699999999998`,
25.499999999999996`}, {25.499999999999996`,
24.604699999999998`}, {25.499999999999996`, 23.5}, {
25.499999999999996`, 2.5}}}]}, {
FaceForm[{
RGBColor[0.5, 0.5, 0.5],
Opacity[1.]}],
FilledCurveBox[{{{0, 2, 0}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {
0, 1, 0}, {1, 3, 3}, {1, 3, 3}}}, {{{3.52539,
1.0035199999999997`}, {2.5250000000000004`,
1.0035199999999997`}, {1.6828099999999997`,
1.0035199999999997`}, {0.9999999999999999,
1.6863299999999997`}, {0.9999999999999999,
2.5285199999999994`}, {0.9999999999999999,
23.474999999999998`}, {0.9999999999999999,
24.317199999999996`}, {1.6828099999999997`,
24.999999999999996`}, {2.5250000000000004`,
24.999999999999996`}, {3.52539, 24.999999999999996`}, {
3.52539, 24.999999999999993`}, {5.523440000000001,
22.421099999999996`}, {5.523440000000001, 13.0039}, {
5.523440000000001, 3.5867199999999992`}, {3.52539,
1.0035199999999997`}, {3.52539, 1.0035199999999997`}}}],
FilledCurveBox[{{{0, 2, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {
0, 1, 0}, {0, 1, 0}}}, {{{22.4977, 12.9016}, {
17.740599999999997`, 16.8398}, {17.740599999999997`,
14.856599999999998`}, {19.980900000000002`,
12.947299999999997`}, {17.740599999999997`,
10.946100000000001`}, {17.740599999999997`, 8.96289}, {
22.4977, 12.9016}}}]}, {
FaceForm[{
RGBColor[
0.46093800000000007`, 0.46093800000000007`,
0.46093800000000007`],
Opacity[1.]}],
FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3,
3}}}, {{{9.202339999999998, 12.8969}, {9.202339999999998,
12.335499999999998`}, {8.747660000000002,
11.880899999999999`}, {8.186329999999998,
11.880899999999999`}, {7.624999999999999,
11.880899999999999`}, {7.169920000000001,
12.335499999999998`}, {7.169920000000001, 12.8969}, {
7.169920000000001, 13.458200000000003`}, {7.624999999999999,
13.9133}, {8.186329999999998, 13.9133}, {8.747660000000002,
13.9133}, {9.202339999999998, 13.458200000000003`}, {
9.202339999999998, 12.8969}}}],
FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3,
3}}}, {{{12.6227, 12.8969}, {12.6227, 12.335499999999998`}, {
12.167599999999998`, 11.880899999999999`}, {
11.606200000000001`, 11.880899999999999`}, {11.0449,
11.880899999999999`}, {10.589799999999999`,
12.335499999999998`}, {10.589799999999999`, 12.8969}, {
10.589799999999999`, 13.458200000000003`}, {11.0449,
13.9133}, {11.606200000000001`, 13.9133}, {
12.167599999999998`, 13.9133}, {12.6227,
13.458200000000003`}, {12.6227, 12.8969}}}],
FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3,
3}}}, {{{16.042600000000004`, 12.8969}, {16.042600000000004`,
12.335499999999998`}, {15.587499999999999`,
11.880899999999999`}, {15.026199999999998`,
11.880899999999999`}, {14.464799999999997`,
11.880899999999999`}, {14.010199999999998`,
12.335499999999998`}, {14.010199999999998`, 12.8969}, {
14.010199999999998`, 13.458200000000003`}, {
14.464799999999997`, 13.9133}, {15.026199999999998`,
13.9133}, {15.587499999999999`, 13.9133}, {
16.042600000000004`, 13.458200000000003`}, {
16.042600000000004`, 12.8969}}}]}}, AspectRatio -> 1, Axes ->
False, Background -> GrayLevel[0.5], Frame -> True, FrameStyle ->
Directive[
Thickness[Tiny],
GrayLevel[0.7]], FrameTicks -> None, ImageSize ->
Dynamic[{
Automatic,
3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification])}]],
GridBox[{{
RowBox[{
TagBox["\"Name: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox[
StyleBox[
PaneBox[
"\"earth_center_spherical.txt\"", ContentPadding -> False,
FrameMargins -> 0, StripOnInput -> True, BaselinePosition ->
Baseline, ImageSize -> {{1, 300}, Automatic}],
LineBreakWithin -> False], "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Unique ID: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["3", "SummaryItem"]}]}},
GridBoxAlignment -> {
"Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete ->
False, GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
GridBoxSpacings -> {"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
BaseStyle -> {
ShowStringCharacters -> False, NumberMarks -> False,
PrintPrecision -> 3, ShowSyntaxStyles -> False}]}},
GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}},
AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
BaselinePosition -> {1, 1}], True -> GridBox[{{
PaneBox[
ButtonBox[
DynamicBox[
FEPrivate`FrontEndResource["FEBitmaps", "SummaryBoxCloser"]],
ButtonFunction :> (Typeset`open$$ = False), Appearance -> None,
BaseStyle -> {}, Evaluator -> Automatic, Method ->
"Preemptive"], Alignment -> {Center, Center}, ImageSize ->
Dynamic[{
Automatic,
3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification])}]],
GraphicsBox[{
Thickness[0.0016806722689075631`], {
FaceForm[{
GrayLevel[0.93],
Opacity[1.]}],
FilledCurveBox[{{{1, 4, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {
1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}}}, {{{
25.499999999999996`, 2.5}, {25.499999999999996`,
1.3953100000000003`}, {24.604699999999998`,
0.49999999999999994`}, {23.5, 0.49999999999999994`}, {2.5,
0.49999999999999994`}, {1.3953100000000003`,
0.49999999999999994`}, {0.49999999999999994`,
1.3953100000000003`}, {0.49999999999999994`, 2.5}, {
0.49999999999999994`, 23.5}, {0.49999999999999994`,
24.604699999999998`}, {1.3953100000000003`,
25.499999999999996`}, {2.5, 25.499999999999996`}, {23.5,
25.499999999999996`}, {24.604699999999998`,
25.499999999999996`}, {25.499999999999996`,
24.604699999999998`}, {25.499999999999996`, 23.5}, {
25.499999999999996`, 2.5}}}]}, {
RGBColor[0.699951, 0.699951, 0.699951],
Opacity[1.],
JoinForm[{ElisionsDump`Miter, 10.}],
JoinedCurveBox[{{{1, 4, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {
1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}}}, {{{
25.499999999999996`, 2.5}, {25.499999999999996`,
1.3953100000000003`}, {24.604699999999998`,
0.49999999999999994`}, {23.5, 0.49999999999999994`}, {2.5,
0.49999999999999994`}, {1.3953100000000003`,
0.49999999999999994`}, {0.49999999999999994`,
1.3953100000000003`}, {0.49999999999999994`, 2.5}, {
0.49999999999999994`, 23.5}, {0.49999999999999994`,
24.604699999999998`}, {1.3953100000000003`,
25.499999999999996`}, {2.5, 25.499999999999996`}, {23.5,
25.499999999999996`}, {24.604699999999998`,
25.499999999999996`}, {25.499999999999996`,
24.604699999999998`}, {25.499999999999996`, 23.5}, {
25.499999999999996`, 2.5}}}]}, {
FaceForm[{
RGBColor[0.5, 0.5, 0.5],
Opacity[1.]}],
FilledCurveBox[{{{0, 2, 0}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {
0, 1, 0}, {1, 3, 3}, {1, 3, 3}}}, {{{3.52539,
1.0035199999999997`}, {2.5250000000000004`,
1.0035199999999997`}, {1.6828099999999997`,
1.0035199999999997`}, {0.9999999999999999,
1.6863299999999997`}, {0.9999999999999999,
2.5285199999999994`}, {0.9999999999999999,
23.474999999999998`}, {0.9999999999999999,
24.317199999999996`}, {1.6828099999999997`,
24.999999999999996`}, {2.5250000000000004`,
24.999999999999996`}, {3.52539, 24.999999999999996`}, {
3.52539, 24.999999999999993`}, {5.523440000000001,
22.421099999999996`}, {5.523440000000001, 13.0039}, {
5.523440000000001, 3.5867199999999992`}, {3.52539,
1.0035199999999997`}, {3.52539, 1.0035199999999997`}}}],
FilledCurveBox[{{{0, 2, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {
0, 1, 0}, {0, 1, 0}}}, {{{22.4977, 12.9016}, {
17.740599999999997`, 16.8398}, {17.740599999999997`,
14.856599999999998`}, {19.980900000000002`,
12.947299999999997`}, {17.740599999999997`,
10.946100000000001`}, {17.740599999999997`, 8.96289}, {
22.4977, 12.9016}}}]}, {
FaceForm[{
RGBColor[
0.46093800000000007`, 0.46093800000000007`,
0.46093800000000007`],
Opacity[1.]}],
FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3,
3}}}, {{{9.202339999999998, 12.8969}, {9.202339999999998,
12.335499999999998`}, {8.747660000000002,
11.880899999999999`}, {8.186329999999998,
11.880899999999999`}, {7.624999999999999,
11.880899999999999`}, {7.169920000000001,
12.335499999999998`}, {7.169920000000001, 12.8969}, {
7.169920000000001, 13.458200000000003`}, {7.624999999999999,
13.9133}, {8.186329999999998, 13.9133}, {8.747660000000002,
13.9133}, {9.202339999999998, 13.458200000000003`}, {
9.202339999999998, 12.8969}}}],
FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3,
3}}}, {{{12.6227, 12.8969}, {12.6227, 12.335499999999998`}, {
12.167599999999998`, 11.880899999999999`}, {
11.606200000000001`, 11.880899999999999`}, {11.0449,
11.880899999999999`}, {10.589799999999999`,
12.335499999999998`}, {10.589799999999999`, 12.8969}, {
10.589799999999999`, 13.458200000000003`}, {11.0449,
13.9133}, {11.606200000000001`, 13.9133}, {
12.167599999999998`, 13.9133}, {12.6227,
13.458200000000003`}, {12.6227, 12.8969}}}],
FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3,
3}}}, {{{16.042600000000004`, 12.8969}, {16.042600000000004`,
12.335499999999998`}, {15.587499999999999`,
11.880899999999999`}, {15.026199999999998`,
11.880899999999999`}, {14.464799999999997`,
11.880899999999999`}, {14.010199999999998`,
12.335499999999998`}, {14.010199999999998`, 12.8969}, {
14.010199999999998`, 13.458200000000003`}, {
14.464799999999997`, 13.9133}, {15.026199999999998`,
13.9133}, {15.587499999999999`, 13.9133}, {
16.042600000000004`, 13.458200000000003`}, {
16.042600000000004`, 12.8969}}}]}}, AspectRatio -> 1, Axes ->
False, Background -> GrayLevel[0.5], Frame -> True, FrameStyle ->
Directive[
Thickness[Tiny],
GrayLevel[0.7]], FrameTicks -> None, ImageSize ->
Dynamic[{
Automatic,
3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification])}]],
GridBox[{{
RowBox[{
TagBox["\"Name: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox[
PaneBox[
"\"earth_center_spherical.txt\"",
ImageSize -> {{1, 500}, Automatic}, BaselinePosition ->
Baseline, ContentPadding -> False, FrameMargins -> 0,
StripOnInput -> True], "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Unique ID: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["3", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Binary: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["False", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Open: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox[
DynamicBox[
ToBoxes[Options[
OutputStream["earth_center_spherical.txt", 3]] =!= {},
StandardForm], UpdateInterval -> 1], "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Encoding: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["Automatic", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Format: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["StandardForm", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Page width: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["78", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Number marks: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["Automatic", "SummaryItem"]}]}},
GridBoxAlignment -> {
"Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete ->
False, GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
GridBoxSpacings -> {"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
BaseStyle -> {
ShowStringCharacters -> False, NumberMarks -> False,
PrintPrecision -> 3, ShowSyntaxStyles -> False}]}},
GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}},
AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
BaselinePosition -> {1, 1}]},
Dynamic[Typeset`open$$], ImageSize -> Automatic]},
"SummaryPanel"],
DynamicModuleValues:>{}], "]"}],
OutputStream["earth_center_spherical.txt", 3],
Editable->False,
SelectWithContents->True,
Selectable->False]], "Output",
CellChangeTimes->{3.967434333513813*^9, 3.9674343744040813`*^9,
3.9674344611943398`*^9},
CellLabel->"Out[18]=",ExpressionUUID->"868cd156-008f-d64d-9e1c-8e498693327d"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"a", "=",
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{"arrowArray", "=",
RowBox[{"{", "}"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"AppendTo", "[",
RowBox[{"arrowArray", ",",
RowBox[{"{",
RowBox[{"Green", ",",
RowBox[{"Ball", "[",
RowBox[{"a", ",", "6370"}], "]"}]}], "}"}]}], "]"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Graphics3D", "[",
RowBox[{"arrowArray", ",",
RowBox[{"Axes", "\[Rule]", "True"}], ",",
RowBox[{"Boxed", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}]}], " ", "]"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Print", "[",
RowBox[{"\"\<lengthMatrix=\>\"", ",", "lengthMatrix"}], "]"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"For", "[",
RowBox[{
RowBox[{"ii", "=", "3"}], ",",
RowBox[{"ii", "<", "lengthMatrix"}], ",", " ",
RowBox[{"ii", "++"}], ",", "\[IndentingNewLine]", "\n", " ",
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Ray", "=",
RowBox[{"RayMatrix", "[",
RowBox[{"[", "ii", "]"}], "]"}]}], " ", ";", "\[IndentingNewLine]",
RowBox[{"spherical", "=",
RowBox[{"{",
RowBox[{
RowBox[{"Re", "[",
RowBox[{"Ray", "[",
RowBox[{"[", "1", "]"}], "]"}], "]"}], ",",
RowBox[{"Re", "[",
RowBox[{"Ray", "[",
RowBox[{"[", "2", "]"}], "]"}], "]"}], ",",
RowBox[{"Re", "[",
RowBox[{"Ray", "[",
RowBox[{"[", "3", "]"}], "]"}], "]"}]}], "}"}]}], ";",
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"Print", "[",
RowBox[{"\"\<Spherical+\>\"", ",", "spherical"}], "]"}], ";"}], " ",
"*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ", "FromSphericalCoordinates", " ", "*)"}],
"\[IndentingNewLine]", " ",
RowBox[{"b", "=",
RowBox[{"FromSphericalCoordinates", "[", "spherical", "]"}]}], ";",
" ", "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"bb", "=",
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{"spherical", "[",
RowBox[{"[", "1", "]"}], "]"}], "*",
RowBox[{"Cos", "[",
RowBox[{"spherical", "[",
RowBox[{"[", "3", "]"}], "]"}], "]"}], "*",
RowBox[{"Sin", "[",
RowBox[{"spherical", "[",
RowBox[{"[", "2", "]"}], "]"}], "]"}]}], ",",
RowBox[{
RowBox[{"spherical", "[",
RowBox[{"[", "1", "]"}], "]"}], "*",
RowBox[{"Sin", "[",
RowBox[{"spherical", "[",
RowBox[{"[", "3", "]"}], "]"}], "]"}], "*",
RowBox[{"Sin", "[",
RowBox[{"spherical", "[",
RowBox[{"[", "2", "]"}], "]"}], "]"}]}], ",",
RowBox[{
RowBox[{"spherical", "[",
RowBox[{"[", "1", "]"}], "]"}], "*",
RowBox[{"Cos", "[",
RowBox[{"spherical", "[",
RowBox[{"[", "2", "]"}], "]"}], "]"}]}]}], "}"}]}], ";"}], " ",
"*)"}], "\[IndentingNewLine]",
RowBox[{"R", "=",
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{
RowBox[{"b", "[",
RowBox[{"[", "1", "]"}], "]"}], "^", "2"}], "+",
RowBox[{
RowBox[{"b", "[",
RowBox[{"[", "2", "]"}], "]"}], "^", "2"}], "+",
RowBox[{
RowBox[{"b", "[",
RowBox[{"[", "3", "]"}], "]"}], "^", "2"}]}], "]"}]}], ";",
"\[IndentingNewLine]", " ",
RowBox[{"Print", "[",
RowBox[{"\"\<ii=\>\"", ",", "ii", ",", "\"\< b=\>\"", ",", "b", ",",
" ", "\"\< magnitude b=\>\"", ",", "R", " ", ",",
" ", "\"\< (R - EARTHR)= \>\"", ",",
RowBox[{"R", "-", "EARTHR"}]}], " ", "]"}], ";", " ",
"\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ",
RowBox[{"AppendTo", "[",
RowBox[{"arrowArray", ",",
RowBox[{"{",
RowBox[{"Red", ",",
RowBox[{"Point", "[",
RowBox[{"{",
RowBox[{"a", ",", "b"}], "}"}], " ", "]"}]}], "}"}]}], "]"}],
";"}]}], " ", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"AppendTo", "[",
RowBox[{"arrowArray", ",",
RowBox[{"{",
RowBox[{"Red", ",",
RowBox[{"Arrow", "[",
RowBox[{"{",
RowBox[{"a", ",", "b"}], "}"}], " ", "]"}]}], "}"}]}], "]"}],
";"}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
"\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{"Print", " ", "[",
RowBox[{"\"\<arrowArray=\>\"", ",", "arrowArray"}], "]"}], ";"}],
"*)"}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"Graphics3D", "[",
RowBox[{"arrowArray", ",",
RowBox[{"Axes", "\[Rule]", "True"}], ",",
RowBox[{"Boxed", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}]}], "]"}], ";"}], " ",
"*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]"}]}]], "Input",
CellChangeTimes->{{3.755460160174726*^9, 3.755460200851099*^9}, {
3.755460867164393*^9, 3.7554608696101847`*^9}, {3.755461109551591*^9,
3.755461168134016*^9}, {3.755461222559847*^9, 3.755461225073848*^9}, {
3.755461266072176*^9, 3.755461289662408*^9}, {3.755461340848538*^9,
3.755461352157701*^9}, {3.7554613900973577`*^9, 3.755461449990507*^9}, {
3.755461589566959*^9, 3.755461619871396*^9}, {3.755463118441763*^9,
3.755463126184739*^9}, 3.75546315757937*^9, 3.7554636478915443`*^9, {
3.755464749687929*^9, 3.7554647897473803`*^9}, {3.755466771693224*^9,
3.755466774660864*^9}, {3.755466831213698*^9, 3.7554668354604*^9}, {
3.755467133843219*^9, 3.755467145920827*^9}, {3.755467324005324*^9,
3.755467334735742*^9}, {3.7554673830966597`*^9, 3.755467418515942*^9}, {
3.7554674620970993`*^9, 3.755467466636093*^9}, {3.755467638108069*^9,
3.755467638523624*^9}, {3.7554677940417347`*^9, 3.755467794707505*^9}, {
3.7554679523070908`*^9, 3.755467955262204*^9}, {3.755468165531446*^9,
3.755468166386797*^9}, {3.755468265416511*^9, 3.7554683358189363`*^9}, {
3.755468422710059*^9, 3.755468463299807*^9}, {3.755468666343663*^9,
3.755468739513462*^9}, 3.755468770988674*^9, {3.755468815015856*^9,
3.755468821191153*^9}, {3.755469030255931*^9, 3.755469033941498*^9},
3.75546930322577*^9, {3.7555121188523073`*^9, 3.7555121247926283`*^9}, {
3.755513114731526*^9, 3.755513121472557*^9}, {3.755555558730844*^9,
3.755555590016902*^9}, {3.7555556365887957`*^9, 3.755555727513669*^9},
3.75556328754053*^9, {3.7555637437099257`*^9, 3.7555641418577538`*^9}, {
3.755564184927292*^9, 3.755564199731595*^9}, {3.755564330562236*^9,
3.755564560287678*^9}, {3.7555645911477337`*^9, 3.755564625043189*^9}, {
3.755564676613998*^9, 3.755564724401286*^9}, {3.755564814474552*^9,
3.7555648439775953`*^9}, {3.7555651937514067`*^9, 3.755565194710245*^9}, {
3.755565249275519*^9, 3.755565262466814*^9}, {3.755565313838909*^9,
3.755565349472348*^9}, {3.755565721227339*^9, 3.7555657253073473`*^9}, {
3.755566054914033*^9, 3.755566057758024*^9}, {3.755566104710037*^9,
3.755566131974083*^9}, {3.755566208619033*^9, 3.755566209475679*^9},
3.755580940857008*^9, {3.7555810838111897`*^9, 3.7555810989217997`*^9}, {
3.755581233092634*^9, 3.755581233926701*^9}, {3.755581456726119*^9,
3.755581461196878*^9}, {3.755583512081643*^9, 3.7555835143683853`*^9}, {
3.7556346840056*^9, 3.755634685434472*^9}, {3.755634815077203*^9,
3.755634832453389*^9}, {3.755854621460703*^9, 3.755854740375824*^9}, {
3.755854855788068*^9, 3.7558549353133717`*^9}, {3.755855148810149*^9,
3.755855152583824*^9}, {3.755858534656389*^9, 3.7558585586944113`*^9}, {
3.7558644296283083`*^9, 3.755864444714732*^9}, {3.755864924971488*^9,
3.755864944233095*^9}, {3.7558662057476254`*^9, 3.755866209383171*^9}},
CellLabel->"In[21]:=",ExpressionUUID->"3446ef60-298e-42b6-8d5d-36899cdd63e6"],
Cell[BoxData[
RowBox[{"{", "}"}]], "Output",
CellChangeTimes->{3.967434333619812*^9, 3.9674343747132816`*^9,
3.9674344619963017`*^9},
CellLabel->"Out[22]=",ExpressionUUID->"70d6c3e4-6771-a54c-9431-08c0913bd1e4"],
Cell[CellGroupData[{
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"lengthMatrix=\"\>", "\[InvisibleSpace]", "604"}],
SequenceForm["lengthMatrix=", 604],
Editable->False]], "Print",
CellChangeTimes->{3.967434333631321*^9, 3.967434374714283*^9,
3.967434461997303*^9},
CellLabel->
"During evaluation of \
In[21]:=",ExpressionUUID->"1df84846-a2bc-c34e-9da3-104b3ce49513"],
Cell[BoxData[