-
Notifications
You must be signed in to change notification settings - Fork 0
/
FAST3D-parabolaTest7.nb
8560 lines (8416 loc) · 447 KB
/
FAST3D-parabolaTest7.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 11.3' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 457089, 8552]
NotebookOptionsPosition[ 445110, 8399]
NotebookOutlinePosition[ 445474, 8415]
CellTagsIndexPosition[ 445431, 8412]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"Needs", "[", "\"\<NDSolve`FEM`\>\"", "]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SetSystemOptions", "[",
RowBox[{"\"\<FiniteElementOptions\>\"", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\<DefaultNumberOfElements\>\"", "\[Rule]", "10"}], "}"}]}],
"]"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SetSystemOptions", "[",
RowBox[{"\"\<FiniteElementOptions\>\"", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\<InterpolationToleranceFactor\>\"", "\[Rule]", "1."}],
"}"}]}], "]"}], ";"}], "\[IndentingNewLine]",
RowBox[{"SetSystemOptions", "[",
RowBox[{"\"\<FiniteElementOptions\>\"", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\<SymbolicProcessing\>\"", "\[Rule]", "100."}], "}"}]}],
"]"}]}], "Input",
CellChangeTimes->{{3.7592312870016127`*^9, 3.759231317855836*^9}, {
3.7592313546623707`*^9, 3.759231380174076*^9}, 3.759232259776616*^9, {
3.7592323029217033`*^9, 3.759232325208775*^9}, 3.759232385911024*^9,
3.759232459586309*^9, {3.759232503878454*^9, 3.759232504156756*^9}, {
3.75923255782311*^9, 3.759232558139102*^9}},
CellLabel->"In[2]:=",ExpressionUUID->"509f3a06-88e9-4e4a-971f-f484c01d33d3"],
Cell[BoxData[
RowBox[{"\<\"FiniteElementOptions\"\>", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"\<\"CacheInterpolationElements\"\>", "\[Rule]", "True"}], ",",
RowBox[{"\<\"DefaultBounds\"\>", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"-", "1"}], ",", "1"}], "}"}]}], ",",
RowBox[{"\<\"DefaultElementMeshOrder\"\>", "\[Rule]", "2"}], ",",
RowBox[{"\<\"DefaultExtrapolationHandler\"\>", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"Indeterminate", "&"}], ",",
RowBox[{"\<\"WarningMessage\"\>", "\[Rule]", "True"}]}], "}"}]}], ",",
RowBox[{"\<\"DefaultNumberOfElements\"\>", "\[Rule]", "10"}], ",",
RowBox[{"\<\"InterpolationToleranceFactor\"\>", "\[Rule]", "1.`"}], ",",
RowBox[{"\<\"MinimumElementMeshQuality\"\>", "\[Rule]", "0.`"}], ",",
RowBox[{"\<\"SymbolicProcessing\"\>", "\[Rule]", "100.`"}]}],
"}"}]}]], "Output",
CellChangeTimes->{
3.759231238809537*^9, {3.759231289671607*^9, 3.759231318392024*^9}, {
3.759231368057869*^9, 3.759231381198244*^9}, 3.75923226061362*^9, {
3.7592323036147833`*^9, 3.759232325867476*^9}, 3.759232386482686*^9,
3.759232461482924*^9, 3.759232505633773*^9, 3.759232558738502*^9,
3.7592330475105677`*^9, 3.75924659575353*^9, 3.759247011760198*^9,
3.759247241427747*^9, 3.759255627094892*^9, 3.7596186586256037`*^9},
CellLabel->"Out[5]=",ExpressionUUID->"81fd9c28-f82d-420d-bcf6-89e08f4c05fa"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"hull1", "=",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[LessEqual]", "10.`"}], "&&",
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"2.215384615384615`", "\[VeryThinSpace]", "+",
RowBox[{"1.159647461884313`", " ",
SuperscriptBox["x", "2."]}]}]}], "&&",
RowBox[{"0.", "<", "y", "\[LessEqual]", " ", "10."}]}], ")"}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"hull2", "=",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[LessEqual]", "10.`"}], "&&",
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.7538461538461538`", "\[VeryThinSpace]", "+",
RowBox[{"0.3102990715481003`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{
"10.", "\[LessEqual]", " ", "y", "\[LessEqual]", " ", "20."}]}], ")"}],
",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"hull3", "=",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[LessEqual]", "10.`"}], "&&",
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.061538461538461535`", "\[VeryThinSpace]", "+",
RowBox[{"0.16809777948184113`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{
"20.", "\[LessEqual]", " ", "y", "\[LessEqual]", " ", "30."}]}], ")"}],
",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"hull4", "=",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[LessEqual]", "10.`"}], "&&",
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.13846153846153844`", "\[VeryThinSpace]", "+",
RowBox[{"0.11225568226713982`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{
"30.", "\[LessEqual]", " ", "y", "\[LessEqual]", " ", "40."}]}], ")"}],
",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"hull5", "=",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[LessEqual]", "10.`"}], "&&",
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.9846153846153846`", "\[VeryThinSpace]", "+",
RowBox[{"0.08302998710242387`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{
"40.", "\[LessEqual]", " ", "y", "\[LessEqual]", " ", "50."}]}], ")"}],
",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"hull6", "=",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[LessEqual]", "10.`"}], "&&",
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"2.6`", "\[VeryThinSpace]", "+",
RowBox[{"0.06526115558326687`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{"50.", "\[LessEqual]", " ", "y", "<", "58."}]}], ")"}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"hullreg", " ", "=",
RowBox[{"RegionUnion", "[",
RowBox[{
"hull1", ",", "hull2", ",", "hull3", ",", " ", "hull4", ",", " ", "hull5",
",", " ", "hull6"}], "]"}]}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{"waterdensity", " ", "=", " ", "1."}], "\[IndentingNewLine]",
RowBox[{"Volume", "[", "hullreg", "]"}], "\[IndentingNewLine]"}], "Input",
CellChangeTimes->{{3.7590709082786827`*^9, 3.759070976680941*^9}, {
3.7590710164343243`*^9, 3.759071086100382*^9}, 3.759188116258915*^9, {
3.75918816398985*^9, 3.7591882026296964`*^9}, {3.7591882759664655`*^9,
3.7591882853554068`*^9}, {3.759188332776496*^9, 3.7591883396282296`*^9}, {
3.759188374776904*^9, 3.759188376003834*^9}, {3.759188516054896*^9,
3.759188527537196*^9}, {3.7591887519323597`*^9, 3.759188765469929*^9}, {
3.759188800325141*^9, 3.759188816436474*^9}, 3.7591888606400876`*^9, {
3.759188918694499*^9, 3.759188938406581*^9}, {3.7591889842913303`*^9,
3.7591891330553503`*^9}, {3.7591891851592536`*^9, 3.759189211156124*^9}, {
3.7591899694242325`*^9, 3.759189973666709*^9}, {3.759190973952363*^9,
3.759191033811428*^9}, {3.7591914538635054`*^9, 3.759191456949845*^9}, {
3.7591928346226873`*^9, 3.7591928374629154`*^9}, {3.759192905840386*^9,
3.7591929087606525`*^9}, 3.7591929773381405`*^9, {3.7591930843332806`*^9,
3.759193084785929*^9}, {3.75919324958808*^9, 3.7591932831106415`*^9}, {
3.7591933485838013`*^9, 3.759193355292145*^9}, {3.7591934156245823`*^9,
3.7591934493357754`*^9}, {3.7591958377415714`*^9,
3.7591958855074897`*^9}, {3.7591959261356907`*^9, 3.759195936437623*^9}, {
3.759195970779643*^9, 3.7591960015091343`*^9}, {3.7591960625351934`*^9,
3.759196143677266*^9}, {3.759231202778841*^9, 3.759231219233143*^9}, {
3.759231257964735*^9, 3.759231272368237*^9}, {3.759231406906158*^9,
3.759231487761743*^9}, {3.75923153295437*^9, 3.7592315661025667`*^9}, {
3.759231600891844*^9, 3.759231601020112*^9}, {3.759231657919354*^9,
3.759231665954249*^9}, {3.7592319961869783`*^9, 3.7592320084834127`*^9}, {
3.7592430216919603`*^9, 3.7592432108054333`*^9}, {3.759246608030121*^9,
3.759246685871416*^9}},
CellLabel->"In[6]:=",ExpressionUUID->"2f6ff26d-771a-4b03-a8f4-d7ffe4270b88"],
Cell[BoxData[
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{
RowBox[{"z", "\[LessEqual]", "10.`"}], "&&",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"2.215384615384615`", "\[VeryThinSpace]", "+",
RowBox[{"1.159647461884313`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{"y", ">", "0"}], "&&",
RowBox[{"y", "\[LessEqual]", "10.`"}]}], ")"}], "||",
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"2.6`", "\[VeryThinSpace]", "+",
RowBox[{"0.06526115558326687`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{"y", "<", "58.`"}], "&&",
RowBox[{"y", "\[GreaterEqual]", "50.`"}]}], ")"}], "||",
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.9846153846153846`", "\[VeryThinSpace]", "+",
RowBox[{"0.08302998710242387`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{"y", "\[GreaterEqual]", "40.`"}], "&&",
RowBox[{"y", "\[LessEqual]", "50.`"}]}], ")"}], "||",
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.13846153846153844`", "\[VeryThinSpace]", "+",
RowBox[{"0.11225568226713982`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{"y", "\[GreaterEqual]", "30.`"}], "&&",
RowBox[{"y", "\[LessEqual]", "40.`"}]}], ")"}], "||",
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.061538461538461535`", "\[VeryThinSpace]", "+",
RowBox[{"0.16809777948184113`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{"y", "\[GreaterEqual]", "20.`"}], "&&",
RowBox[{"y", "\[LessEqual]", "30.`"}]}], ")"}], "||",
RowBox[{"(",
RowBox[{
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.7538461538461538`", "\[VeryThinSpace]", "+",
RowBox[{"0.3102990715481003`", " ",
SuperscriptBox["x", "2.`"]}]}]}], "&&",
RowBox[{"y", "\[GreaterEqual]", "10.`"}], "&&",
RowBox[{"y", "\[LessEqual]", "20.`"}]}], ")"}]}], ")"}]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.759188120365203*^9,
3.759188207894758*^9, 3.759188289719022*^9, 3.7591883455011263`*^9,
3.7591883839906178`*^9, 3.759188530724363*^9, {3.759188754722724*^9,
3.7591888171367264`*^9}, {3.7591890517720194`*^9,
3.7591890797836294`*^9}, {3.759189133713914*^9, 3.759189159968935*^9}, {
3.7591892115466604`*^9, 3.7591892223785925`*^9}, 3.7591899768159018`*^9,
3.7591901754738593`*^9, 3.7591902152142916`*^9, 3.7591903647086954`*^9,
3.759190466463532*^9, {3.759191036086466*^9, 3.7591910410587077`*^9},
3.759191100511945*^9, 3.7591914608430557`*^9, 3.759191613308679*^9,
3.7591928416955366`*^9, 3.7591929119777985`*^9, 3.759193090535079*^9, {
3.759193267218831*^9, 3.7591932917473707`*^9}, 3.7591933593193345`*^9, {
3.7591934166276407`*^9, 3.759193456845788*^9}, 3.7591945536547327`*^9,
3.7591948491173277`*^9, 3.7591950772002106`*^9, 3.7591952715910254`*^9,
3.759195479261425*^9, {3.7591958472870445`*^9, 3.7591960020553913`*^9},
3.7591960854514437`*^9, 3.759196144367815*^9, 3.759231243787692*^9, {
3.7592312769900103`*^9, 3.7592313199027443`*^9}, 3.7592313697604237`*^9,
3.7592314147513237`*^9, 3.7592314882450323`*^9, 3.759231606674178*^9,
3.759231667079802*^9, 3.759231710336154*^9, 3.759232012178594*^9,
3.759232512186304*^9, 3.759232560956871*^9, 3.7592330487104273`*^9,
3.759243044881804*^9, {3.759243131207652*^9, 3.759243211425343*^9},
3.759246597189707*^9, 3.759246687913784*^9, 3.75924701347337*^9,
3.759247244078658*^9, 3.7592556297720203`*^9, 3.759618661353381*^9},
CellLabel->"Out[12]=",ExpressionUUID->"aa8f3bce-1b06-437f-b193-88ab8a24b6cb"],
Cell[BoxData["1.`"], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.759188120365203*^9,
3.759188207894758*^9, 3.759188289719022*^9, 3.7591883455011263`*^9,
3.7591883839906178`*^9, 3.759188530724363*^9, {3.759188754722724*^9,
3.7591888171367264`*^9}, {3.7591890517720194`*^9,
3.7591890797836294`*^9}, {3.759189133713914*^9, 3.759189159968935*^9}, {
3.7591892115466604`*^9, 3.7591892223785925`*^9}, 3.7591899768159018`*^9,
3.7591901754738593`*^9, 3.7591902152142916`*^9, 3.7591903647086954`*^9,
3.759190466463532*^9, {3.759191036086466*^9, 3.7591910410587077`*^9},
3.759191100511945*^9, 3.7591914608430557`*^9, 3.759191613308679*^9,
3.7591928416955366`*^9, 3.7591929119777985`*^9, 3.759193090535079*^9, {
3.759193267218831*^9, 3.7591932917473707`*^9}, 3.7591933593193345`*^9, {
3.7591934166276407`*^9, 3.759193456845788*^9}, 3.7591945536547327`*^9,
3.7591948491173277`*^9, 3.7591950772002106`*^9, 3.7591952715910254`*^9,
3.759195479261425*^9, {3.7591958472870445`*^9, 3.7591960020553913`*^9},
3.7591960854514437`*^9, 3.759196144367815*^9, 3.759231243787692*^9, {
3.7592312769900103`*^9, 3.7592313199027443`*^9}, 3.7592313697604237`*^9,
3.7592314147513237`*^9, 3.7592314882450323`*^9, 3.759231606674178*^9,
3.759231667079802*^9, 3.759231710336154*^9, 3.759232012178594*^9,
3.759232512186304*^9, 3.759232560956871*^9, 3.7592330487104273`*^9,
3.759243044881804*^9, {3.759243131207652*^9, 3.759243211425343*^9},
3.759246597189707*^9, 3.759246687913784*^9, 3.75924701347337*^9,
3.759247244078658*^9, 3.7592556297720203`*^9, 3.75961866135559*^9},
CellLabel->"Out[13]=",ExpressionUUID->"7d8bc7ac-4222-4910-bbc2-377bebf305bf"],
Cell[BoxData["5286.276962145505`"], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.759188120365203*^9,
3.759188207894758*^9, 3.759188289719022*^9, 3.7591883455011263`*^9,
3.7591883839906178`*^9, 3.759188530724363*^9, {3.759188754722724*^9,
3.7591888171367264`*^9}, {3.7591890517720194`*^9,
3.7591890797836294`*^9}, {3.759189133713914*^9, 3.759189159968935*^9}, {
3.7591892115466604`*^9, 3.7591892223785925`*^9}, 3.7591899768159018`*^9,
3.7591901754738593`*^9, 3.7591902152142916`*^9, 3.7591903647086954`*^9,
3.759190466463532*^9, {3.759191036086466*^9, 3.7591910410587077`*^9},
3.759191100511945*^9, 3.7591914608430557`*^9, 3.759191613308679*^9,
3.7591928416955366`*^9, 3.7591929119777985`*^9, 3.759193090535079*^9, {
3.759193267218831*^9, 3.7591932917473707`*^9}, 3.7591933593193345`*^9, {
3.7591934166276407`*^9, 3.759193456845788*^9}, 3.7591945536547327`*^9,
3.7591948491173277`*^9, 3.7591950772002106`*^9, 3.7591952715910254`*^9,
3.759195479261425*^9, {3.7591958472870445`*^9, 3.7591960020553913`*^9},
3.7591960854514437`*^9, 3.759196144367815*^9, 3.759231243787692*^9, {
3.7592312769900103`*^9, 3.7592313199027443`*^9}, 3.7592313697604237`*^9,
3.7592314147513237`*^9, 3.7592314882450323`*^9, 3.759231606674178*^9,
3.759231667079802*^9, 3.759231710336154*^9, 3.759232012178594*^9,
3.759232512186304*^9, 3.759232560956871*^9, 3.7592330487104273`*^9,
3.759243044881804*^9, {3.759243131207652*^9, 3.759243211425343*^9},
3.759246597189707*^9, 3.759246687913784*^9, 3.75924701347337*^9,
3.759247244078658*^9, 3.7592556297720203`*^9, 3.759618661659857*^9},
CellLabel->"Out[14]=",ExpressionUUID->"ad0a7f51-d224-49bf-a82f-7abba4b44e08"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"Mass", "[",
RowBox[{"boatregion_", ",", " ", "boatdensity_"}], "]"}], ":=",
"1400."}], "\[IndentingNewLine]",
RowBox[{"boatmass", "=", "1000."}], "\[IndentingNewLine]",
RowBox[{"boatcom", " ", "=", " ",
RowBox[{"{",
RowBox[{"0", ",", "34.16", ",", " ", "2.5"}], "}"}]}]}], "Input",
CellChangeTimes->{{3.75907084715088*^9, 3.7590709012747087`*^9},
3.759071138603026*^9, {3.759071259499156*^9, 3.759071266928605*^9}, {
3.759071300947072*^9, 3.759071307409132*^9}, 3.759071376767318*^9, {
3.7591899840798655`*^9, 3.7591899848487515`*^9}, {3.759190121162506*^9,
3.759190126490387*^9}, {3.759231845123139*^9, 3.7592318602254477`*^9}, {
3.759251439590783*^9, 3.759251476239232*^9}, {3.759253306941372*^9,
3.759253307045931*^9}, {3.759253448364478*^9, 3.7592534484828157`*^9}, {
3.759253537476088*^9, 3.7592535375769997`*^9}, {3.759255457258191*^9,
3.759255457331435*^9}, {3.7592592631705847`*^9, 3.759259263352057*^9}, {
3.759262744551395*^9, 3.759262744589581*^9}, {3.7593434501253347`*^9,
3.759343503661457*^9}, 3.7593457199445267`*^9, {3.75934639664704*^9,
3.759346399249103*^9}, {3.7593472843614893`*^9, 3.759347287275352*^9}, {
3.759619478431678*^9, 3.759619494941098*^9}},
CellLabel->
"In[102]:=",ExpressionUUID->"148f209f-1c4c-4a67-bd1b-e73e94ce68a7"],
Cell[BoxData["1000.`"], "Output",
CellChangeTimes->{3.759259263929592*^9, 3.7593435060707407`*^9,
3.759345720642428*^9, 3.75934640035091*^9, 3.759347288269042*^9,
3.759618662637723*^9, 3.759619495736557*^9},
CellLabel->
"Out[103]=",ExpressionUUID->"6f662560-c5dd-4e09-b144-1a34a799698b"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",", "34.16`", ",", "2.5`"}], "}"}]], "Output",
CellChangeTimes->{3.759259263929592*^9, 3.7593435060707407`*^9,
3.759345720642428*^9, 3.75934640035091*^9, 3.759347288269042*^9,
3.759618662637723*^9, 3.759619495742299*^9},
CellLabel->
"Out[104]=",ExpressionUUID->"778973bd-fab2-41f4-93f8-582b01e60ee0"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData["1400.`"], "Input",
CellChangeTimes->{{3.759259253589293*^9, 3.759259253654072*^9}},
CellLabel->"In[46]:=",ExpressionUUID->"45362d03-ccfa-46ee-833d-2447bb7bda85"],
Cell[BoxData["1400.`"], "Output",
CellChangeTimes->{
3.7591901270370865`*^9, 3.7591901842152624`*^9, 3.7591902170444508`*^9,
3.7591903679788785`*^9, 3.7591904682806416`*^9, 3.7591910413266835`*^9,
3.75919110077754*^9, 3.759191461129629*^9, 3.7591916135899086`*^9,
3.7591928435437527`*^9, 3.759192912277458*^9, 3.7591930948705883`*^9,
3.75919329235915*^9, 3.759193457457187*^9, 3.7591945543443193`*^9,
3.7591948497920103`*^9, 3.7591950778719225`*^9, 3.7591952722982407`*^9,
3.759195479933139*^9, 3.759196015161483*^9, 3.7591961520562525`*^9, {
3.759231845446526*^9, 3.759231861540543*^9}, 3.759233055095726*^9,
3.759246692000321*^9, 3.759247016850753*^9, 3.759247245832082*^9, {
3.759251453654491*^9, 3.759251476764043*^9}, 3.7592533084729443`*^9,
3.7592533755081673`*^9, 3.759253449016714*^9, 3.7592535381777782`*^9,
3.759255458483305*^9, 3.759255634325502*^9, 3.759259255471898*^9},
CellLabel->"Out[46]=",ExpressionUUID->"91f981d3-7d8a-4226-94ce-dc3c03bc8bda"]
}, Open ]],
Cell[BoxData[{
RowBox[{
RowBox[{"Waterline3D", "[",
RowBox[{"\[Theta]_", ",", " ", "\[Phi]_"}], "]"}], ":=", " ",
RowBox[{"Quiet", "[",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"y", ",", " ",
RowBox[{"\[Phi]", "\[Equal]", "90"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"x", ",", " ",
RowBox[{"\[Theta]", "==", "90"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-", "z"}], ",", " ",
RowBox[{
RowBox[{"\[Theta]", ">", "90"}], "||",
RowBox[{"\[Phi]", ">", "90"}]}]}], "}"}]}], "}"}], ",", "z"}],
"]"}], "\[LessEqual]",
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"draft", ",", " ",
RowBox[{
RowBox[{"\[Theta]", "\[Equal]", "90"}], "||",
RowBox[{"\[Phi]", "\[Equal]", "90"}]}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{"-", "1"}], "*",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "*", "x"}], "+",
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "*", "y"}], "+",
"draft"}], ")"}]}], ",", " ",
RowBox[{
RowBox[{"\[Theta]", ">", "90"}], "||",
RowBox[{"\[Phi]", ">", "90"}]}]}], "}"}]}], "}"}], ",",
RowBox[{
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "*", "x"}], "+",
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "*", "y"}], "+",
"draft"}]}], "]"}]}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",",
RowBox[{"-", "15"}], ",", "15"}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"y", ",", "0", ",", "60"}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"z", ",", "0", ",", "30"}], "}"}]}], "}"}]}], "]"}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"CalculateDraft3D", "[",
RowBox[{"boatregion_", ",", " ", "\[Theta]_", ",", " ", "\[Phi]_"}], "]"}],
":=",
RowBox[{"Quiet", "[",
RowBox[{
RowBox[{"FindRoot", "[",
RowBox[{
RowBox[{
RowBox[{"Integrate", "[",
RowBox[{"1.", ",", " ",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}], "\[Element]",
RowBox[{"RegionIntersection", "[",
RowBox[{"boatregion", ",", " ",
RowBox[{"Waterline3D", "[",
RowBox[{"\[Theta]", ",", " ", "\[Phi]"}], "]"}]}], "]"}]}]}],
"]"}], "\[Equal]",
RowBox[{"boatmass", "/", "1."}]}], ",",
RowBox[{"{",
RowBox[{"draft", ",", " ",
RowBox[{"1", "-",
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "/", "3"}], "-",
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "/", "3"}]}]}], "}"}]}],
"]"}], "[",
RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"DisplacedWater3D", "[",
RowBox[{"boatregion_", ",", " ", "\[Theta]_", ",", " ", "\[Phi]_"}], "]"}],
":=", " ",
RowBox[{"Quiet", "[",
RowBox[{
RowBox[{"RegionIntersection", "[",
RowBox[{"boatregion", ",", " ",
RowBox[{"Waterline3D", "[",
RowBox[{"\[Theta]", ",", " ", "\[Phi]"}], "]"}]}], "]"}], "/.",
RowBox[{"CalculateDraft3D", "[",
RowBox[{"boatregion", ",", " ", "\[Theta]", ",", " ", "\[Phi]"}],
"]"}]}], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"CenterOfBuoyancy3D", "[",
RowBox[{"boatregion_", ",", "\[Theta]_", ",", " ", "\[Phi]_"}], "]"}], ":=",
RowBox[{"Quiet", "[",
RowBox[{"RegionCentroid", "[",
RowBox[{"DisplacedWater3D", "[",
RowBox[{"boatregion", ",", " ", "\[Theta]", ",", " ", "\[Phi]"}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"RightingMoment3D", "[",
RowBox[{"boatregion_", ",", " ", "\[Theta]_", ",", " ", "\[Phi]_"}], "]"}],
":=",
RowBox[{"Quiet", "[",
RowBox[{"Cross", "[",
RowBox[{
RowBox[{
RowBox[{"CenterOfBuoyancy3D", "[",
RowBox[{"boatregion", ",", " ", "\[Theta]", ",", " ", "\[Phi]"}],
"]"}], "-", "boatcom"}], ",", " ",
RowBox[{"boatmass", "*", "9.8", "*", "0.001", "*",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"-",
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}], ",", " ",
RowBox[{"-",
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}], ",", "1"}], "}"}], "/",
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{
RowBox[{"-", "1"}], ",", " ",
RowBox[{
RowBox[{"\[Theta]", ">", "90"}], "||",
RowBox[{"\[Phi]", ">", "90"}]}]}], "}"}], "}"}], ",", "1"}],
"]"}]}], "/",
RowBox[{"Norm", "[",
RowBox[{"{",
RowBox[{
RowBox[{"-",
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}], ",", " ",
RowBox[{"-",
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}], ",", "1"}], "}"}],
"]"}]}]}]}], "]"}], "]"}]}]}], "Input",
CellChangeTimes->{{3.75907150036087*^9, 3.75907151420553*^9}, {
3.759071689998314*^9, 3.759071692052773*^9}, {3.759071729200437*^9,
3.7590717906234617`*^9}, {3.7590718315538177`*^9, 3.75907184764012*^9}, {
3.759190139449133*^9, 3.7591901402817287`*^9}, {3.7591945137812686`*^9,
3.7591945507377234`*^9}, {3.7591962012005296`*^9,
3.7591962031752744`*^9}, {3.7592318052310867`*^9, 3.759231837116363*^9}, {
3.7592320308419333`*^9, 3.759232031476864*^9}, {3.759232595207486*^9,
3.7592326043528748`*^9}, {3.75923265937434*^9, 3.759232662472035*^9}, {
3.75924696984558*^9, 3.759246985264812*^9}, 3.759253436990367*^9, {
3.759255578825295*^9, 3.7592555854337683`*^9}},
CellLabel->
"In[105]:=",ExpressionUUID->"0440646f-917b-4090-ad63-896091dbb2a8"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"f", "[",
RowBox[{"angle1_", ",", " ", "angle2_"}], "]"}], ":=",
RowBox[{"N", "[",
RowBox[{"RightingMoment3D", "[",
RowBox[{"hullreg", ",", " ", "angle1", ",", " ", "angle2"}], "]"}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"angleTable1", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "0.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "0."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}], "//",
"AbsoluteTiming"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"angleTable2", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "2.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "2."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}], "//",
"AbsoluteTiming"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"angleTable3", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "4.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "4."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}], "//",
"AbsoluteTiming"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"angleTable4", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "6.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "6."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}], "//",
"AbsoluteTiming"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"angleTable5", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "8.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "8."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}], "//",
"AbsoluteTiming"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"angleTable6", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "10.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "10."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}], "//",
"AbsoluteTiming"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"angleTable7", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "12.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "12."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}], "//",
"AbsoluteTiming"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"angleTable8", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "14.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "14."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}], "//", "AbsoluteTiming"}],
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"angleTable9", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "16.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "16."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}], "//",
RowBox[{"AbsoluteTiming", "\[IndentingNewLine]",
RowBox[{
RowBox[{"angleTable10", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "18.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "18."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}]}]}], "//",
RowBox[{"AbsoluteTiming", "\[IndentingNewLine]",
RowBox[{
RowBox[{"angleTable11", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "20.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "20."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}],
"]"}], "]"}], "]"}]}], ";"}]}]}], "//", "AbsoluteTiming"}],
"\[IndentingNewLine]", "*)"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"data", " ", "=", " ",
RowBox[{"Union", "[",
RowBox[{
"angleTable1", ",", " ", "angleTable2", ",", " ", "angleTable3", ",", " ",
"angleTable4", ",", " ", "angleTable5", ",", " ", "angleTable6", ",",
"angleTable7", ",", " ", "angleTable8"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"xdata", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "1"}], "]"}], "]"}], ",", " ",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "2"}], "]"}], "]"}], ",", " ",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "3", ",", " ", "1"}], "]"}], "]"}]}], "}"}],
",", " ",
RowBox[{"{",
RowBox[{"i", ",", " ", "1", ",", " ",
RowBox[{"Length", "[", "data", "]"}], ",", "1"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"ydata", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "1"}], "]"}], "]"}], ",", " ",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "2"}], "]"}], "]"}], ",", " ",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "3", ",", " ", "2"}], "]"}], "]"}]}], "}"}],
",", " ",
RowBox[{"{",
RowBox[{"i", ",", " ", "1", ",", " ",
RowBox[{"Length", "[", "data", "]"}], ",", "1"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"zdata", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "1"}], "]"}], "]"}], ",", " ",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "2"}], "]"}], "]"}], ",", " ",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "3", ",", " ", "3"}], "]"}], "]"}]}], "}"}],
",", " ",
RowBox[{"{",
RowBox[{"i", ",", " ", "1", ",", " ",
RowBox[{"Length", "[", "data", "]"}], ",", "1"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"normdata", " ", "=", " ",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "1"}], "]"}], "]"}], ",", " ",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "2"}], "]"}], "]"}], ",", " ",
RowBox[{"Quiet", "[",
RowBox[{"Norm", "[",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "3"}], "]"}], "]"}], "]"}], "]"}]}],
"}"}], ",", " ",
RowBox[{"{",
RowBox[{"i", ",", " ", "1", ",", " ",
RowBox[{"Length", "[", "data", "]"}], ",", "1"}], "}"}]}], "]"}],
"]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{"Quiet", "[",
RowBox[{"ListPlot3D", "[",
RowBox[{"normdata", ",",
RowBox[{"Mesh", "\[Rule]", "None"}], ",",
RowBox[{"InterpolationOrder", "\[Rule]", "3"}], ",",
RowBox[{"ColorFunction", "\[Rule]", "\"\<SouthwestColors\>\""}]}], "]"}],
"]"}], "\[IndentingNewLine]",
RowBox[{"Show", "[",
RowBox[{
RowBox[{"ListPlot3D", "[",
RowBox[{"ydata", ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"120", ",", "140"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "15"}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}], ",",
RowBox[{"Mesh", "\[Rule]", "None"}], ",",
RowBox[{"InterpolationOrder", "\[Rule]", "3"}], ",",
RowBox[{"ColorFunction", "\[Rule]", "\"\<SouthwestColors\>\""}]}], "]"}],
",",
RowBox[{"ListPlot3D", "[",
RowBox[{"xdata", ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"120", ",", "140"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "15"}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}], ",",
RowBox[{"Mesh", "\[Rule]", "None"}], ",",
RowBox[{"InterpolationOrder", "\[Rule]", "3"}], ",",
RowBox[{"ColorFunction", "\[Rule]", "\"\<SouthwestColors\>\""}]}],
"]"}]}], "]"}], "\[IndentingNewLine]",
RowBox[{"ListPlot3D", "[",
RowBox[{"xdata", ",",
RowBox[{"Mesh", "\[Rule]", "None"}], ",",
RowBox[{"InterpolationOrder", "\[Rule]", "3"}], ",",
RowBox[{"ColorFunction", "\[Rule]", "\"\<SouthwestColors\>\""}]}],
"]"}], "\[IndentingNewLine]",
RowBox[{"ListPlot3D", "[",
RowBox[{"ydata", ",",
RowBox[{"Mesh", "\[Rule]", "None"}], ",",
RowBox[{"InterpolationOrder", "\[Rule]", "3"}], ",",
RowBox[{"ColorFunction", "\[Rule]", "\"\<SouthwestColors\>\""}]}],
"]"}], "\[IndentingNewLine]"}], "Input",
CellChangeTimes->{{3.759071972719798*^9, 3.759072000466646*^9}, {
3.7590728096669416`*^9, 3.7590728370525837`*^9}, {3.759072925626684*^9,
3.759073043280994*^9}, {3.7590733389563932`*^9, 3.759073395969692*^9}, {
3.759073737971199*^9, 3.759073745068993*^9}, {3.759196286946718*^9,
3.759196350019063*^9}, {3.7591968986470656`*^9, 3.7591969019949427`*^9}, {
3.759196953944845*^9, 3.7591969568600082`*^9}, 3.7591970221001887`*^9, {
3.7591973124142914`*^9, 3.7591973418859725`*^9}, {3.759197599361195*^9,
3.759197599533065*^9}, {3.759197840707432*^9, 3.759197840879014*^9}, {
3.7592327247668533`*^9, 3.75923278672547*^9}, {3.759233074147759*^9,
3.759233089106287*^9}, {3.759233164075914*^9, 3.7592331641651173`*^9}, {
3.7592332034436283`*^9, 3.7592332264305553`*^9}, {3.759233329409387*^9,
3.759233330562262*^9}, {3.759233374509461*^9, 3.759233375405281*^9}, {
3.759235205597747*^9, 3.7592352096624928`*^9}, {3.759246554546398*^9,
3.7592465783867483`*^9}, {3.75924720921434*^9, 3.759247230869266*^9}, {
3.7592486515739813`*^9, 3.759248754141222*^9}, {3.759248788809889*^9,
3.759248800420772*^9}, {3.759248841941307*^9, 3.759248904059903*^9},
3.759248974456441*^9, {3.759249700079706*^9, 3.759249784582097*^9}, {
3.759255657163827*^9, 3.759255662485325*^9}, {3.759344774221476*^9,
3.759344784093149*^9}},
CellLabel->
"In[110]:=",ExpressionUUID->"b677da27-0834-4a54-b1eb-70a1b1a91310"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"26.017744`", ",", "Null"}], "}"}]], "Output",
CellChangeTimes->{3.7592334028077707`*^9, 3.759235241768835*^9,
3.759246468304467*^9, 3.759246590194846*^9, 3.759246733316518*^9,
3.759247099012374*^9, 3.7592472377982893`*^9, 3.759247281489523*^9,
3.759247979162485*^9, 3.7592489209219418`*^9, 3.759249033223909*^9,
3.759249688861747*^9, 3.759249815976777*^9, 3.759250612075322*^9,
3.759251662967119*^9, 3.759255707007372*^9, 3.759343547648552*^9,
3.759345768277549*^9, 3.759346447240469*^9, 3.759347329004874*^9,
3.759618708362731*^9, 3.759619545764861*^9},
CellLabel->
"Out[111]=",ExpressionUUID->"afa328e4-b6e9-4695-88be-794aa61132cb"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"79.865405`", ",", "Null"}], "}"}]], "Output",
CellChangeTimes->{3.7592334028077707`*^9, 3.759235241768835*^9,
3.759246468304467*^9, 3.759246590194846*^9, 3.759246733316518*^9,
3.759247099012374*^9, 3.7592472377982893`*^9, 3.759247281489523*^9,
3.759247979162485*^9, 3.7592489209219418`*^9, 3.759249033223909*^9,
3.759249688861747*^9, 3.759249815976777*^9, 3.759250612075322*^9,
3.759251662967119*^9, 3.759255707007372*^9, 3.759343547648552*^9,
3.759345768277549*^9, 3.759346447240469*^9, 3.759347329004874*^9,
3.759618708362731*^9, 3.759619625629786*^9},
CellLabel->
"Out[112]=",ExpressionUUID->"6c789a9a-7725-4409-807f-cd9168c1ba7e"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"80.896742`", ",", "Null"}], "}"}]], "Output",
CellChangeTimes->{3.7592334028077707`*^9, 3.759235241768835*^9,
3.759246468304467*^9, 3.759246590194846*^9, 3.759246733316518*^9,
3.759247099012374*^9, 3.7592472377982893`*^9, 3.759247281489523*^9,
3.759247979162485*^9, 3.7592489209219418`*^9, 3.759249033223909*^9,
3.759249688861747*^9, 3.759249815976777*^9, 3.759250612075322*^9,
3.759251662967119*^9, 3.759255707007372*^9, 3.759343547648552*^9,
3.759345768277549*^9, 3.759346447240469*^9, 3.759347329004874*^9,
3.759618708362731*^9, 3.759619706527018*^9},
CellLabel->
"Out[113]=",ExpressionUUID->"591eae45-bde0-4acc-857f-d20ed63b6505"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"77.283469`", ",", "Null"}], "}"}]], "Output",
CellChangeTimes->{3.7592334028077707`*^9, 3.759235241768835*^9,
3.759246468304467*^9, 3.759246590194846*^9, 3.759246733316518*^9,
3.759247099012374*^9, 3.7592472377982893`*^9, 3.759247281489523*^9,
3.759247979162485*^9, 3.7592489209219418`*^9, 3.759249033223909*^9,
3.759249688861747*^9, 3.759249815976777*^9, 3.759250612075322*^9,