Skip to content

Instantly share code, notes, and snippets.

@pukpr
Created June 15, 2023 03:05
Show Gist options
  • Save pukpr/91241973cb825abd663484a884c8a2d8 to your computer and use it in GitHub Desktop.
Save pukpr/91241973cb825abd663484a884c8a2d8 to your computer and use it in GitHub Desktop.
function Tide_Sum (Template : in Data_Pairs;
Constituents : in Long_Periods_Amp_Phase;
Periods : in Long_Periods;
Ref_Time : in Long_Float := 0.0;
Scaling : in Long_Float := 1.0;
Cos_Phase : in Boolean := True;
Year_Len : in Long_Float := Year_Length;
Integ: in Long_Float := 0.0
) return Data_Pairs is
Pi : Long_Float := Ada.Numerics.Pi;
Time : Long_Float;
Res : Data_Pairs := Template;
Yr : Long_Float := Year_Len;
Drac : Long_Float := 2.0*pi*Yr/Draconic;
Trop : Long_Float := 2.0*pi*Yr/Tropical;
Anom : Long_Float := 2.0*pi*Yr/Anomalistic;
Nodal : Long_Float := Drac-Trop;
Peri : Long_Float := (Trop-Anom)/2.0;
K1 : Long_Float := Constituents(1).Value;
K2 : Long_Float := Constituents(2).Value;
TA : Long_Float := Constituents(3).Value;
TP : Long_Float := Constituents(4).Value;
AA : Long_Float := Constituents(5).Value;
AP : Long_Float := Constituents(6).Value;
T2A : Long_Float := Constituents(7).Value;
T2P : Long_Float := Constituents(8).Value;
SA : Long_Float := Constituents(9).Value;
SP : Long_Float := Constituents(10).Value;
A2A : Long_Float := Constituents(11).Value;
SYA : Long_Float := Constituents(12).Value;
SYP : Long_Float := Constituents(13).Value;
EA : Long_Float := Constituents(14).Value;
EP : Long_Float := Constituents(15).Value;
E2A : Long_Float := Constituents(16).Value;
YA : Long_Float := Constituents(17).Value;
YP : Long_Float := Constituents(18).Value;
Q1 : Long_Float := Constituents(19).Value;
Q2 : Long_Float := Constituents(20).Value;
DA : Long_Float := Constituents(21).Value;
DP : Long_Float := Constituents(22).Value;
K3 : Long_Float := Constituents(23).Value;
NA : Long_Float := Constituents(24).Value;
NP : Long_Float := Constituents(25).Value;
PA : Long_Float := Constituents(26).Value;
PP : Long_Float := Constituents(27).Value;
N2A : Long_Float := Constituents(28).Value;
N2P : Long_Float := Constituents(29).Value;
function Ds is new gsin (DA, 2.0*Drac, DP); -- = A*sin(2*Drac*(x+p));
function Ts is new gsin (TA, 2.0*Trop, TP); -- = A*sin(2*Trop*(x+p))
function T2s is new gsin (T2A, Trop+Drac, T2P); -- = A*sin((Trop+Drac)*(x+p))
function As is new gsin (AA, Anom, AP); -- = A*sin(Anom*(x+p))
function Semi is new gsin(SA, 4.0*Pi, SP); -- = A*sin(4*pi*(x+p))
function Evect is new gsin(EA, 2.0*Trop-Anom-4.0*Pi, EP); -- = A*sin((2*Trop-Anom-4*pi)*(x+p))
function Syn is new gsin(SYA, 2.0*Trop-4.0*Pi, SYP); -- = A*sin((2*Trop-4*pi)*(x+p))
function Annual is new gsin(YA, 2.0*Pi, YP); -- = A*cos(2*pi*(x+p))
function N is new gsin (NA, Nodal, NP); -- = A*sin(Nodal*(x+p))
function N2 is new gsin (N2A, Nodal/2.0, N2P); -- = A*sin(Nodal/2.0*(x+p))
function P is new gsin (PA, Peri, PP); -- = A*sin(Peri*(x+p))
-- # Define the model
function Model (Time : in Long_Float) return Long_Float is
MM : Long_Float := As(Time);
Msm : Long_Float := Evect(Time);
Mfp : Long_Float := T2s(Time + A2A*Mm + E2A*Msm);
Mf : Long_Float := Ts(Time + A2A*Mm + E2A*Msm);
Mfd : Long_Float := Ds(Time + A2A*Mm + E2A*Msm);
Msf : Long_Float := Syn(Time + Q1*Mf + Q2*Mfp);
begin
return Mf + Mfp + Mfd + Semi(Time) + Annual(Time) + Msm + Msf + N(Time) + N2(Time) + P(Time) +
As(Time + K1*(Mf+Mfp+Mfd) + K2*Mm + K3*Msf);
end Model;
begin
for I in Template'Range loop
Time := Template(I).Date - Ref_Time;
declare
TF : Long_Float := 0.0;
begin
Tf := Model (Time);
Res(I) := (Time, Scaling * TF);
end;
end loop;
return Res;
end Tide_Sum;
@pukpr
Copy link
Author

pukpr commented Jun 17, 2023

Tahiti -- same impulse as NINO4

image

   0.00116679885 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.17394464215:delA:
   0.27681578514:delB:
   9.17248760826:asym:
   0.15189805087:ann1:
  -0.97194502199:ann2:
  -0.03802357357:sem1:
   2.80640655167:sem2:
   0.00007646987:year:
  -0.11170347461:IR:
   0.00001248790 :mA:
   0.00075920093 :mP:
   0.00829059162 :shiftT:
  -0.81493725939 :init:
---- Tidal ----
   9.10846048884,    0.01604637681 1
  29.53063068368,   -0.01205477904 2
  27.09267692660,    0.35052067867 3
   7.08840411155,   -1.30045630029 4
   6.85248390316,    0.09965713020 5
  13.60611040750,    1.02820337889 6
  13.66083077700,    0.12976727508 7
3396.73824406533,    1.42389311091 8
   9.12068919638,    0.01007363343 9
  13.77727494300,   -0.02547242202 10
6793.47648813065,   -0.00382519574 11
  27.55454988600,    0.00153270698 12
  27.66676714572,    0.96599191529 13
  27.44323926226,    0.03752834353 14
1616.30271425126,   -0.82745095759 15
   6.85940288609,    0.00497079271 16
   7.09580797237,    0.00052601963 17
   9.55688260704,    0.48261834841 18
   9.13295078376,    0.12069498072 19
2190.35004466729,    0.11157935187 20
  27.21222081500,    0.01297859790 21
   9.18484996200,    0.93958775600 22
  14.76531534184,   -0.21577538545 23
  31.81203822360,    0.01339275114 24
   9.54345712458,  -27.02209230059 25
  13.63341568476,   -0.01448958000 26
  26.98505934729,    2.13981316731 27
 365.24996486987,    0.02730344078 28
 182.62498243493,   -8.37855224833 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
   0.06301900681 :K0:
  -0.22947329216 :level:
   0.08582897141,    0.23259069160,    0.34796345962 1
   0.17165794282,    0.07855607388,    1.92096384305 2
   0.25748691424,    0.07650430032,    2.17305847630 3
   0.34331588565,    0.11635141902,    1.22348903505 4
   0.42914485706,    0.07604287285,   -3.11723121523 5
   0.51497382847,    0.07909499727,   -2.52682920542 6
   0.60080279988,    0.11651441368,   -1.24066068129 7
   0.68663177130,    0.11354449515,   -2.08772971380 8
   0.77246074271,    0.15293990011,   -2.20358719524 9
   0.85828971412,    0.14987343157,   -2.00596920559 10
   0.94411868553,    0.08685228748,    0.12216479211 11
   1.02994765694,    0.09033262547,   -1.09401497683 12
   1.11577662836,    0.02617613654,   -1.04285486147 13
   1.20160559977,    0.05959286658,   -1.04430967033 14
   1.28743457118,    0.07677794637,    0.87984654910 15
   1.37326354259,    0.09663599903,    1.92072482584 16
   1.45909251400,    0.10929619745,    1.49427922013 17
   1.54492148542,    0.12231130737,    2.91468982240 18
   1.63075045683,    0.23602500874,    2.37135078512 19
   1.71657942824,    0.00908593044,   -1.74888194053 20
   1.80240839965,    0.15561901496,   -2.62134364650 21
   1.88823737106,    0.05140365621,    0.24320232669 22
   1.97406634248,    0.35878988835,   -2.00801475253 23
   2.05989531389,    0.04952875805,    0.16755698564 24
   2.14572428530,    0.12845411330,   -0.75195773003 25
   2.23155325671,    0.04857134050,   -1.60662021443 26
   2.31738222812,    0.10385822254,   -0.03203454268 27
   2.40321119954,    0.04072199222,    0.26184527743 28
CC   0.2742535755   0.6565710410   4 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\backup.bat  tahiti_constituents2crossvalid3posimpulseA

higher LTE factors
image

   0.00116562074 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.15329986950:delA:
   0.27681578514:delB:
   1.73377458700:asym:
   0.15636743052:ann1:
  -1.01598381811:ann2:
  -0.03940733478:sem1:
   2.72411687591:sem2:
   0.00007627082:year:
  -0.11392584361:IR:
   0.00001213764 :mA:
   0.00075380668 :mP:
   0.00829765333 :shiftT:
  -0.89486200845 :init:
---- Tidal ----
   9.10846048884,    0.01657544787 1
  29.53063068499,   -0.00886511689 2
  27.09267692660,    0.35135966126 3
   7.08840411170,   -1.30045548674 4
   6.85248390316,    0.10026505650 5
  13.60611040750,    1.02826861662 6
  13.66083077700,    0.12923748967 7
3396.73824406533,    1.42392607347 8
   9.12068919638,    0.00652150889 9
  13.77727494300,   -0.02780189713 10
6793.47648813065,   -0.00437229915 11
  27.55454988600,    0.00153240873 12
  27.66676714572,    0.96582795356 13
  27.44323926226,    0.03870375474 14
1616.30271425126,   -0.82782481062 15
   6.85940288609,    0.00408749592 16
   7.09580797252,    0.00052583688 17
   9.55688260731,    0.55002616720 18
   9.13295078376,    0.12055163718 19
2190.35004466729,    0.11480054214 20
  27.21222081500,    0.01392263213 21
   9.18484996200,    0.93973979420 22
  14.76531534249,   -0.24438358806 23
  31.81203822662,    0.01417441386 24
   9.54345712486,  -26.90196032780 25
  13.63341568476,   -0.01419544058 26
  26.98505934729,    1.91736609980 27
 365.24996467082,    0.02796334217 28
 182.62498233541,   -8.24330058115 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
   0.06224618152 :K0:
  -0.23551962389 :level:
   0.08584797576,    0.14065436213,    0.28267354958 1
   0.17169595152,    0.08436494113,    2.01775546914 2
   3.86315890916,    0.10445889712,   -0.42156884493 45
   0.34339190304,    0.12252367217,    1.40919397720 4
   4.37824676372,    0.11386330226,   -0.20968803337 51
   3.17637510309,    0.16290125735,    0.38124875346 37
   0.60093583031,    0.09013005094,   -0.86746114230 7
   0.68678380607,    0.11907672582,   -2.15111112782 8
   0.77263178183,    0.17505228714,   -2.05784372659 9
   0.85847975759,    0.10638356134,   -1.71978364412 10
   0.94432773335,    0.14321863637,    0.84951474059 11
   3.26222307885,    0.13354196823,   -1.72105671865 38
   2.91883117581,    0.14027756213,    0.75309132961 34
   3.51976700613,    0.04695908511,    1.37959400508 41
   1.28771963639,    0.07830212367,    1.05200244382 15
   1.37356761215,    0.06131667364,    2.41258397030 16
   4.63579069100,    0.11556153608,   -2.17622408113 54
   1.54526356367,    0.09336454004,    3.12233045098 18
   1.63111153942,    0.25365758641,    2.85128612423 19
   1.20187166063,    0.10016618792,   -0.67792987796 14
   1.80280749094,    0.19726052358,   -1.75948032718 21
   1.88865546670,    0.07752315795,    0.48558989375 22
   1.97450344246,    0.35010197363,   -1.30486671670 23
   3.60561498189,    0.03973666187,    2.87704008747 42
   2.14619939398,    0.05694087302,   -0.29028344619 25
   2.31789534550,    0.06938569694,    1.18363830809 27
   1.11602368487,    0.09273698808,    0.69971102785 13
   2.06035141822,    0.02991741679,    2.07643967339 24
CC   0.1736748707   0.7043405732   4 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\backup.bat  tahiti_constituents2crossvalid3negimpulsehighC

@pukpr
Copy link
Author

pukpr commented Jun 17, 2023

image

   0.00112126556 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.18890383581:delA:
   0.27681578514:delB:
   9.17248760826:asym:
   0.05360211484:ann1:
  -0.96085610488:ann2:
  -0.01044401578:sem1:
   1.05828816888:sem2:
   0.00007989622:year:
   0.08580660342:IR:
   0.00001453882 :mA:
   0.00071155368 :mP:
   0.00854255279 :shiftT:
  -0.81998523023 :init:
---- Tidal ----
   9.10846048884,    0.01748549331 1
  29.53063066129,   -0.02018159033 2
  27.09267692660,    0.34868276650 3
   7.08840410897,    1.01123207657 4
   6.85248390316,    0.09936532330 5
  13.60611040750,    1.02748375223 6
  13.66083077700,    0.13106870771 7
3396.73824406533,    1.42390192500 8
   9.12068919638,    0.00584141877 9
  13.77727494300,   -0.01718472889 10
6793.47648813065,   -0.00267408983 11
  27.55454988600,    0.00149075807 12
  27.66676714572,    0.97096979711 13
  27.44323926226,    0.04213981653 14
1616.30271425126,    0.82489556989 15
   6.85940288609,    0.00490422906 16
   7.09580796979,    0.00052717599 17
   9.55688260235,    0.54879352021 18
   9.13295078376,    0.11783445094 19
2190.35004466729,    0.06643929930 20
  27.21222081500,    0.01112262594 21
   9.18484996200,    0.93953611481 22
  14.76531533064,   -0.22450532033 23
  31.81203817161,    0.00254211429 24
   9.54345711991,  -27.25267202085 25
  13.63341568476,   -0.01315177760 26
  26.98505934729,    2.34353228557 27
 365.24996829622,    0.01087019440 28
 182.62498414811,   -5.96600309916 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
  -0.02867808431 :K0:
   0.02884550175 :level:
   0.08162359490,    0.37178720159,   -1.64675809211 1
   0.16324718980,    0.31751291509,    1.16293105284 2
   0.24487078470,    0.26831141095,   -2.27119500859 3
   0.32649437960,    0.23552090241,    0.27555674656 4
   0.40811797450,    0.23424731913,    2.57001624976 5
   0.48974156940,    0.24603156191,   -0.54289781733 6
   0.57136516430,    0.32061375094,    2.49057309253 7
   0.65298875920,    0.21371699358,   -0.83120699858 8
   0.73461235410,    0.17765644429,    1.38594016125 9
   0.81623594900,    0.10202361508,   -1.72584862962 10
   0.89785954390,    0.12088322125,    0.33506106939 11
   0.97948313880,    0.18173245181,   -2.86760205859 12
   1.06110673370,    0.35051091499,    0.02455914592 13
   1.14273032860,    0.18992293431,   -3.13073225142 14
   1.22435392350,    0.16863413075,   -0.70734677928 15
   1.30597751840,    0.18692734714,    2.26075388250 16
   1.38760111330,    0.14710344088,   -1.27715403526 17
   1.46922470820,    0.30653264494,    1.44186513813 18
   1.55084830310,    0.23042889124,   -2.19089344575 19
   1.63247189800,    0.27113452760,    0.63177108504 20
   1.71409549290,    0.24772587024,   -2.66221069322 21
   1.79571908780,    0.10125637934,   -0.19542700303 22
   1.87734268270,    0.21350429932,    3.12565587225 23
   1.95896627760,    0.24007185356,   -0.58564999738 24
   2.04058987250,    0.11062383769,    1.84876277611 25
   2.12221346740,    0.22334019353,   -0.93810870826 26
   2.20383706230,    0.14732107267,    2.14760581184 27
   2.28546065720,    0.23228863121,   -1.07179539354 28
CC   0.8285277243   0.8243587173   8 1
   0.00000000000:dLOD:Main exiting, flushing other tasks
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\restore.bat  nino4_constituents2crossvalid3

image

   0.00113059441 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.17398407088:delA:
   0.27682190605:delB:
   9.17248760826:asym:
   0.05957717659:ann1:
  -0.83662065238:ann2:
  -0.01779593080:sem1:
   1.47276156802:sem2:
   0.00008009524:year:
   0.06559494673:IR:
   0.00001318147 :mA:
   0.00076326358 :mP:
   0.00852912452 :shiftT:
  -0.82761238960 :init:
---- Tidal ----
   9.10846048884,    0.01661376084 1
  29.53063065999,   -0.02120668030 2
  27.09267692660,    0.34911168549 3
   7.08840410882,    1.01123212064 4
   6.85248390316,    0.09768884388 5
  13.60611040750,    1.02748037345 6
  13.66083077700,    0.13158577093 7
3396.73824406533,    1.42388941594 8
   9.12068919638,    0.00657408355 9
  13.77727494300,   -0.01485265883 10
6793.47648813065,   -0.00259306420 11
  27.55454988600,    0.00152528442 12
  27.66676714572,    0.97045860642 13
  27.44323926226,    0.04134510011 14
1616.30271425126,    0.82528346299 15
   6.85940288609,    0.00479502383 16
   7.09580796964,    0.00052755533 17
   9.55688260208,    0.54643398258 18
   9.13295078376,    0.11772032980 19
2190.35004466729,    0.06660572719 20
  27.21222081500,    0.01093988874 21
   9.18484996200,    0.93940726777 22
  14.76531532999,   -0.23016042560 23
  31.81203816860,    0.00067343714 24
   9.54345711963,  -23.57552455229 25
  13.63341568476,   -0.01296690989 26
  26.98505934729,    2.61543611520 27
 365.24996849524,    0.01220426926 28
 182.62498424762,   -7.46616159518 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
  -0.01173551851 :K0:
  -0.14407904240 :level:
   0.08163028353,    0.02737346280,   -2.11689265924 1
   0.16326056706,    0.06741123792,    0.88661788227 2
   3.26521134126,    0.08432484695,    2.49093456930 40
   2.85705992360,    0.05106004422,   -0.48534733872 35
   2.61216907301,    0.05210070290,    1.21578076644 32
   1.71423595416,    0.06862975819,   -2.03281118782 21
   0.57141198472,    0.11990328177,    2.23315514056 7
   0.24489085059,    0.07293539327,   -1.67452252915 3
   0.73467255178,    0.10972803346,    0.52351373190 9
   3.59173247539,    0.08724243593,   -1.33976191401 44
   0.89793311885,    0.02386146153,   -1.37026317514 11
   3.18358105773,    0.05056809726,   -1.80588281917 39
   1.06119368591,    0.18476456196,    0.16760347637 13
   1.14282396944,    0.08193358361,   -2.30420213805 14
   2.53053878948,    0.11009671091,   -1.15668811957 31
   1.30608453650,    0.11056319148,    2.97715930569 16
   1.38771482004,    0.08420062403,   -0.18639276308 17
   1.46934510357,    0.17497232217,    2.07971445203 18
   2.93869020713,    0.05772752013,    2.56972541642 36
   4.48966559423,    0.08665281835,    2.07047974701 55
   3.10195077420,    0.07099231659,    1.31427458009 38
   1.79586623769,    0.05333833910,   -2.66489281094 22
   1.87749652122,    0.13018494112,   -2.38083858571 23
   1.95912680476,    0.10558078629,   -0.17342278719 24
   2.04075708829,    0.05279364145,   -0.72688481972 25
   2.12238737182,    0.11739493843,   -0.55457835671 26
   2.20401765535,    0.09192614015,    2.73717450124 27
   2.28564793888,    0.17260239695,   -0.89478916826 28
CC   0.8006006534   0.8973098111   4 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\backup.bat  nino4_constituents2crossvalid3highB

@pukpr
Copy link
Author

pukpr commented Jun 17, 2023

PDO -- even+odd harmonics

image

   0.00111229585 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.66774568652:delA:
   0.27681578514:delB:
   9.17248760826:asym:
   0.09933802625:ann1:
  -0.75828040131:ann2:
  -0.06476110840:sem1:
   1.76185813085:sem2:
   0.00007836156:year:
  -0.12238173591:IR:
   0.00001232060 :mA:
   0.00056019599 :mP:
   0.00856490117 :shiftT:
  -0.58007765362 :init:
---- Tidal ----
   9.10846048884,    0.01696567027 1
  29.53063067132,   -0.02087015593 2
  27.09267692660,    0.34112943919 3
   7.08840411013,    1.01125148121 4
   6.85248390316,    0.10492569380 5
  13.60611040750,    1.02721399346 6
  13.66083077700,    0.13198784495 7
3396.73824406533,    1.43047225374 8
   9.12068919638,    0.01618355821 9
  13.77727494300,   -0.02442838169 10
6793.47648813065,   -0.00291652307 11
  27.55454988600,    0.00144214042 12
  27.66676714572,    0.97086004517 13
  27.44323926226,    0.04214960282 14
1616.30271425126,    0.82580224942 15
   6.85940288609,    0.00534498140 16
   7.09580797095,    0.00052146730 17
   9.55688260445,    0.46024100877 18
   9.13295078376,    0.11749213750 19
2190.35004466729,    0.06350821389 20
  27.21222081500,    0.01310201239 21
   9.18484996200,    0.93925625527 22
  14.76531533566,   -0.31091937687 23
  31.81203819490,    0.00192327513 24
   9.54345712200,  -27.05779153651 25
  13.63341568476,   -0.01106672959 26
  26.98505934729,    2.45868003102 27
 365.24996676156,    0.00955950651 28
 182.62498338078,   -5.34039381766 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
  -0.00897849710 :K0:
  -0.01389384533 :level:
   0.08167378705,    0.23723136943,   -1.74446041842 1
   0.16334757410,    0.13706613834,    0.82133942872 2
   0.24502136115,    0.03150504098,    1.42999204147 3
   0.32669514820,    0.12799837791,    2.05353893047 4
   0.40836893525,    0.13908810640,    0.02089201168 5
   0.49004272230,    0.12918667595,    0.71096059518 6
   0.57171650935,    0.12549509901,   -0.68176236046 7
   0.65339029640,    0.22805940938,    1.97510861151 8
   0.73506408345,    0.08940470167,   -1.79455158788 9
   0.81673787049,    0.02230549667,   -2.47121164999 10
   0.89841165754,    0.39711667655,    0.69384407892 11
   0.98008544459,    0.03133888324,   -1.61512676674 12
   1.06175923164,    0.25282191433,    2.83701005717 13
   1.14343301869,    0.32273378097,    0.25533692720 14
   1.22510680574,    0.22538985610,   -1.78636149959 15
   1.30678059279,    0.08487381346,    1.40676494663 16
   1.38845437984,    0.15776289258,   -2.57453495079 17
   1.47012816689,    0.17223977599,   -0.50741584747 18
   1.55180195394,    0.02762044639,   -2.40589067727 19
   1.63347574099,    0.39012119575,   -2.38858284797 20
   1.71514952804,    0.18938810551,    1.54111201497 21
   1.79682331509,    0.27888529516,   -0.25712844192 22
   1.87849710214,    0.27240262909,   -1.02708766903 23
   1.96017088919,    0.15233449696,    3.07440186859 24
   2.04184467624,    0.41928278459,    0.85973870374 25
   2.12351846329,    0.28406652464,    0.06245131008 26
   2.20519225034,    0.32892575529,   -2.58335200455 27
   2.28686603739,    0.16866099845,    0.67373474692 28
CC  -0.1747694735   0.7922659281   8 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\backup.bat  pdo_constituents2crossvalid3

harmonics
image

   0.00111180439 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.64994091477:delA:
   0.27630863247:delB:
   5.95458639215:asym:
   0.22609252674:ann1:
  -0.45346158850:ann2:
  -0.08550753219:sem1:
   1.72689847276:sem2:
   0.00008228029:year:
  -0.14722222217:IR:
   0.00000013322 :mA:
   0.00043700738 :mP:
   0.00856067229 :shiftT:
  -0.54036412374 :init:
---- Tidal ----
   9.10846048884,    0.01666869981 1
  29.53063064570,   -0.01916857338 2
  27.09267692660,    0.34143087163 3
   7.08840410718,    1.01125238068 4
   6.85248390316,    0.10515235192 5
  13.60611040750,    1.02745565845 6
  13.66083077700,    0.13181920882 7
3396.73824406533,    1.43043496265 8
   9.12068919638,    0.01655445879 9
  13.77727494300,   -0.03386159850 10
6793.47648813065,   -0.00289520863 11
  27.55454988600,    0.00144478807 12
  27.66676714572,    0.97104704322 13
  27.44323926226,    0.04214680061 14
1616.30271425126,    0.82614491137 15
   6.85940288609,    0.00480132633 16
   7.09580796799,    0.00052359452 17
   9.55688259908,    0.46050148501 18
   9.13295078376,    0.11795164516 19
2190.35004466729,    0.06165977212 20
  27.21222081500,    0.01328755889 21
   9.18484996200,    0.75279367102 22
  14.76531532285,   -0.34435482862 23
  31.81203813544,    0.00424949678 24
   9.54345711665,  -10.39511941352 25
  13.63341568476,   -0.01340906778 26
  26.98505934729,    1.28831743026 27
 365.24997068029,    0.01089873641 28
 182.62498534015,   -6.46485169453 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
  -0.03424358065 :K0:
  -0.05102121934 :level:
   0.08153101077,    0.24729278211,   -1.50329223998 1
   3.42430245242,    0.17005038646,    0.96259203078 42
   2.36439931239,    0.09133692393,    0.14119297100 29
   0.32612404309,    0.11129055077,    2.62682004567 4
   0.40765505386,    0.16934227834,   -0.39009445176 5
   3.58736447397,    0.20670689628,    0.48154052332 44
   4.56573660323,    0.19917531954,   -2.84174777553 56
   0.65224808618,    0.18769505029,    1.90872902450 8
   1.54908920467,    0.19535095972,    0.26959577190 19
   1.95674425853,    0.06556443742,   -3.10451840805 24
   0.89684111849,    0.36845368784,    0.82452689458 11
   4.40267458169,    0.22734266160,   -0.33743294882 54
   1.05990314004,    0.27516087470,    2.61029540592 13
   1.14143415081,    0.30447537291,    0.21302334066 14
   1.22296516158,    0.18106689147,   -1.98947847733 15
   4.48420559246,    0.06284027792,   -2.69382427367 55
   1.38602718312,    0.15308769754,    2.87455638175 17
   1.46755819390,    0.12363234891,   -0.83172507264 18
   4.23961256014,    0.18471127931,    0.08060991925 52
   1.63062021544,    0.38325125364,   -2.42161155439 20
   1.71215122621,    0.15340376977,    1.35110569900 21
   1.79368223698,    0.23906272719,   -0.20898536394 22
   1.87521324776,    0.29453514486,   -1.24239103837 23
   3.66889548474,    0.18346549962,   -2.42479170868 45
   2.03827526930,    0.34094942254,    0.54849950737 25
   2.11980628007,    0.38219199743,   -0.11713529470 26
   2.20133729084,    0.34650142338,   -2.83065831943 27
   2.28286830162,    0.15269295978,    0.73700920294 28
CC  -0.4066941619   0.8643025720   8 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\backup.bat  pdo_constituents2crossvalid3high

@pukpr
Copy link
Author

pukpr commented Jun 20, 2023

Using a sin^N delta of power N= 15

PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ......\GeoEnergyMath\io\backup.bat darwin_constituents2crossvalid3bSINcomplement

Delta variant
2390 ......\GeoEnergyMath\io\restore.bat darwin_constituents2crossvalid3b
2392 ......\GeoEnergyMath\io\restore.bat darwin_constituents2crossvalid3b
2393 $env:CC_START="1880"
2394 $env:CC_END="1932"
2400 ......\GeoEnergyMath\io\backup.bat darwin_constituents2crossvalid3bSIN
2404 ......\GeoEnergyMath\io\backup.bat darwin_constituents2crossvalid3bSINhigh
2406 ......\GeoEnergyMath\io\backup.bat darwin_constituents2crossvalid3bSINhigh-pt8
2408 ......\GeoEnergyMath\io\backup.bat darwin_constituents2crossvalid3bSINhigh-pt85worse
2409 ......\GeoEnergyMath\io\restore.bat darwin_constituents2crossvalid3bSIN

2410 $env:CC_END="2015"
2411 $env:CC_START="1935"
2414 ......\GeoEnergyMath\io\backup.bat darwin_constituents2crossvalid3bSINcomplement
2415 .\enso_opt.exe

image

  0.00111916648 :offset:
   0.00000000000 :bg:
  -9.10072106846 :impA:
   0.28438035101 :impB:
  14.99048476157 :impC:
  -9.09960468396:delA:
   0.27676443679:delB:
   9.17134999257:asym:
   0.07094617114:ann1:
  -2.81079818916:ann2:
  -0.15325271031:sem1:
   0.82592576109:sem2:
   0.00009296189:year:
   0.16091717257:IR:
   0.00001536702 :mA:
  -0.00084083884 :mP:
   0.00860963428 :shiftT:
  -0.59680237727 :init:
---- Tidal ----
   9.10846048884,    0.01735147656 1
  29.53063057588,   -0.02841717728 2
  27.09267692660,    0.34212651950 3
   7.08840409913,    0.98902914949 4
   6.85248390316,    0.10124629564 5
  13.60611040750,    1.02538819660 6
  13.66083077700,    0.12810273401 7
3396.73824406533,    1.42401647831 8
   9.12068919638,    0.01132689320 9
  13.77727494300,   -0.00129361749 10
6793.47648813065,   -0.00141824121 11
  27.55454988600,    0.00154127154 12
  27.66676714572,    1.04013588153 13
  27.44323926226,    0.04801111162 14
1616.30271425126,    0.82350083257 15
   6.85940288609,    0.00477189186 16
   7.09580795993,    0.00053293671 17
   9.55688258446,    0.47846435153 18
   9.13295078376,    0.11822503299 19
2190.35004466729,    0.08081744484 20
  27.21222081500,    0.01239904500 21
   9.18484996200,    0.94334717135 22
  14.76531528794,   -0.36577320399 23
  31.81203797339,    0.01780510042 24
   9.54345710207,  -35.06822898616 25
  13.63341568476,   -0.02476631047 26
  26.98505934729,    3.07001049660 27
 365.24998136189,    0.00486549849 28
 182.62499068094,   -4.22040260892 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
   0.02960125178 :K0:
   0.26312236581 :level:
   0.09829301324,    0.12809910700,   -0.43097963143 1
   0.19658602648,    0.11265463505,   -1.90715906435 2
   0.29487903972,    0.12366899916,    2.53306413829 3
   0.39317205296,    0.13115016410,   -1.00123312594 4
   0.49146506620,    0.09832753942,   -2.41960389136 5
   0.58975807944,    0.05829593852,    0.87636664162 6
   0.68805109268,    0.06254261335,    0.76857327396 7
   0.78634410592,    0.13035683320,   -0.53500952585 8
   0.88463711916,    0.07085795920,    2.69514441824 9
   0.98293013240,    0.10572920781,    1.50357447397 10
   1.08122314564,    0.09659980480,    0.27879539477 11
   1.17951615888,    0.07741724101,   -1.55555411855 12
   1.27780917212,    0.10110073886,   -0.92953866256 13
   1.37610218536,    0.09717414843,    2.82266415713 14
   1.47439519860,    0.20512998393,    0.54275260650 15
   1.57268821184,    0.10603572259,   -0.97981187486 16
   1.67098122508,    0.05799658566,   -0.89035815650 17
   1.76927423832,    0.07720518497,    0.46976104366 18
   1.86756725156,    0.12519995028,   -2.14973644482 19
   1.96586026480,    0.09429808065,   -1.16526740575 20
   2.06415327804,    0.16079515683,    1.56047312721 21
   2.16244629128,    0.14044981063,   -1.51094107709 22
   2.26073930452,    0.10259274726,   -2.43552323023 23
   2.35903231776,    0.05793799112,    1.90795379664 24
   2.45732533100,    0.14164047023,   -1.80620850855 25
   2.55561834424,    0.00420607273,    2.89569948071 26
   2.65391135748,    0.07049823326,   -1.37420827228 27
   2.75220437072,    0.13149912518,    0.08227794618 28
CC   0.4425475022   0.6975592971   6 1

image

  0.00110985420 :offset:
   0.00000000000 :bg:
  -9.10596735501 :impA:
   0.28457288762 :impB:
  14.90969085924 :impC:
  -9.09960468396:delA:
   0.27676443679:delB:
   9.17134999257:asym:
   0.09724201034:ann1:
  -3.04928964789:ann2:
  -0.14742732792:sem1:
   0.86686044813:sem2:
   0.00009287011:year:
   0.17721500732:IR:
   0.00001526087 :mA:
  -0.00083037102 :mP:
   0.00861491148 :shiftT:
  -0.60991844542 :init:
---- Tidal ----
   9.10846048884,    0.01846597017 1
  29.53063057648,   -0.03394950807 2
  27.09267692660,    0.34203486752 3
   7.08840409920,    0.98902914949 4
   6.85248390316,    0.09944457395 5
  13.60611040750,    1.02523107011 6
  13.66083077700,    0.12861229547 7
3396.73824406533,    1.42401647831 8
   9.12068919638,    0.01270231406 9
  13.77727494300,   -0.00126833170 10
6793.47648813065,   -0.00134537377 11
  27.55454988600,    0.00155005248 12
  27.66676714572,    1.04000442849 13
  27.44323926226,    0.04962952737 14
1616.30271425126,    0.82351149714 15
   6.85940288609,    0.00454258727 16
   7.09580795999,    0.00053422352 17
   9.55688258458,    0.48212442132 18
   9.13295078376,    0.11812462122 19
2190.35004466729,    0.07932590507 20
  27.21222081500,    0.01287833701 21
   9.18484996200,    0.94434959805 22
  14.76531528824,   -0.36575811194 23
  31.81203797478,    0.02360084843 24
   9.54345710219,  -35.06931620105 25
  13.63341568476,   -0.02518887575 26
  26.98505934729,    2.96969415112 27
 365.24998127011,    0.00427406948 28
 182.62499063506,   -4.35248443878 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
   0.02810361321 :K0:
   0.26042392472 :level:
   0.09835892336,    0.10525471578,   -0.80460442819 1
   0.19671784672,    0.10439747689,   -1.37300148009 2
   0.29507677007,    0.12683353431,    2.67708204217 3
   0.39343569343,    0.11129386927,   -1.12267685907 4
   4.81958724454,    0.11890742254,   -3.05129070775 49
   0.59015354015,    0.03436432985,    1.20037230725 6
   0.68851246351,    0.08199904050,    0.26497494707 7
   5.21302293797,    0.00479807367,   -0.39225118556 53
   2.95076770074,    0.08866196614,    0.91687894118 30
   0.98358923358,    0.15147976695,    1.49927932639 10
   1.08194815694,    0.06397850461,    0.53860107841 11
   1.18030708029,    0.08787268363,   -1.21055053995 12
   3.73763908760,    0.05950573892,    2.93847148446 38
   1.37702492701,    0.10244065999,   -2.64154640632 14
   1.47538385037,    0.13721345635,    1.04845284335 15
   1.57374277373,    0.07196677228,   -0.78516733690 16
   1.67210169708,    0.09032295539,    0.27646572433 17
   3.54092124088,    0.02439358742,    2.88964533141 36
   3.93435693431,    0.03646694797,    0.82434107520 40
   1.96717846716,    0.09618113945,   -0.08917390613 20
   2.06553739052,    0.19205911552,    2.44930592731 21
   2.16389631387,    0.11084363883,   -1.09774459741 22
   2.26225523723,    0.07339497541,   -1.13121489440 23
   2.36061416059,    0.10706568280,    2.19062642602 24
   2.45897308395,    0.17359820240,   -1.17766558835 25
   4.42615155110,    0.16680404172,   -1.96922413621 45
   2.65569093066,    0.03948719888,   -0.98296303915 27
   2.75404985402,    0.11549894513,    0.81862339353 28
CC   0.5313022260   0.8067687281   4 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\restore.bat  darwin_constituents2crossvalid3bSINhigh-pt8

complement

image

PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> $env:CC_START="1935"
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> $env:CC_END="2015"
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\restore.bat  darwin_constituents2crossvalid3bSINcomplement

   0.00111410351 :offset:
   0.00000000000 :bg:
  -9.10159665685 :impA:
   0.29100278069 :impB:
  15.00099325490 :impC:
  -9.09960468396:delA:
   0.27676443679:delB:
   9.17239635150:asym:
   0.05359581574:ann1:
  -2.40364058872:ann2:
  -0.05127054039:sem1:
   1.86746940586:sem2:
   0.00010058079:year:
   0.08634772665:IR:
   0.00001720056 :mA:
   0.00084345706 :mP:
   0.00859911890 :shiftT:
  -0.43543304884 :init:
---- Tidal ----
   9.10846048884,    0.01803272019 1
  29.53063052608,   -0.03658112291 2
  27.09267692660,    0.34112909106 3
   7.08840409339,    0.98910415544 4
   6.85248390316,    0.10175600791 5
  13.60611040750,    1.02592115418 6
  13.66083077700,    0.12791759738 7
3396.73824406533,    1.42197611951 8
   9.12068919638,    0.01788987735 9
  13.77727494300,   -0.00605340828 10
6793.47648813065,   -0.00145833085 11
  27.55454988600,    0.00151198479 12
  27.66676714572,    1.03895065084 13
  27.44323926226,    0.05075911753 14
1616.30271425126,    0.82355871052 15
   6.85940288609,    0.00320896705 16
   7.09580795417,    0.00053410739 17
   9.55688257403,    0.57306744100 18
   9.13295078376,    0.12435688738 19
2190.35004466729,    0.10023747605 20
  27.21222081500,    0.01365874938 21
   9.18484996200,    0.94249414112 22
  14.76531526304,   -0.33085511500 23
  31.81203785780,    0.01437373677 24
   9.54345709166,  -35.15834522678 25
  13.63341568476,   -0.02561841048 26
  26.98505934729,    2.66228966886 27
 365.24998898079,    0.01042655187 28
 182.62499449039,   -2.06877620682 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
   0.01627185181 :K0:
   0.24607100842 :level:
   0.09821737513,    0.12982921648,    0.37864452113 1
   0.19643475026,    0.09495160464,    0.53876964029 2
   0.29465212539,    0.05200158278,   -1.73246801232 3
   0.39286950052,    0.15327436844,    2.93128993582 4
   0.49108687565,    0.12116171640,   -0.61026199805 5
   0.58930425078,    0.16826850214,    2.72831836262 6
   0.68752162591,    0.16201386758,    0.50713272508 7
   0.78573900104,    0.09778749475,   -2.54260767908 8
   0.88395637617,    0.10042271599,   -2.82285366034 9
   0.98217375130,    0.06308774705,   -2.65791921994 10
   1.08039112643,    0.05307354995,    3.01091203331 11
   1.17860850156,    0.10959776635,    0.89530367263 12
   1.27682587669,    0.12269611081,    2.05908982701 13
   1.37504325182,    0.14001996377,   -0.75310837014 14
   1.47326062695,    0.05450698278,    0.00570651443 15
   1.57147800208,    0.05791597434,    0.85654285798 16
   1.66969537721,    0.15447576611,    0.12500089638 17
   1.76791275234,    0.15176765519,   -2.67725363810 18
   1.86613012747,    0.06050471512,    2.24851720732 19
   1.96434750260,    0.01117398781,   -0.33024610614 20
   2.06256487773,    0.06161422446,   -2.38631185710 21
   2.16078225286,    0.22410084932,   -1.63859807521 22
   2.25899962799,    0.06850570208,    0.75543081880 23
   2.35721700312,    0.13480200795,    1.64377951810 24
   2.45543437825,    0.08348829256,   -0.55605911354 25
   2.55365175338,    0.15922125348,   -1.34910826202 26
   2.65186912851,    0.20931998383,    2.10276326277 27
   2.75008650364,    0.28512382344,    0.52740751024 28
CC   0.3048770993   0.7029153436   2 1

compare to delta version
image

PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ......\GeoEnergyMath\io\restore.bat darwin_constituents2crossvalid3b
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> $env:CC_START="1880"
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> $env:CC_END="1932"

   0.00109088247 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.09960468396:delA:
   0.27676443679:delB:
   9.17134999257:asym:
   0.06889137968:ann1:
  -2.27376860093:ann2:
  -0.15715015236:sem1:
   0.92513196376:sem2:
   0.00008029612:year:
   0.06650862460:IR:
   0.00001532610 :mA:
  -0.00082255133 :mP:
   0.00852819665 :shiftT:
  -0.72843188072 :init:
---- Tidal ----
   9.10846048884,    0.01755781878 1
  29.53063065867,   -0.02837080741 2
  27.09267692660,    0.34519971062 3
   7.08840410867,    1.01123278425 4
   6.85248390316,    0.09753091315 5
  13.60611040750,    1.02538182651 6
  13.66083077700,    0.12773101598 7
3396.73824406533,    1.42389799800 8
   9.12068919638,    0.00964048632 9
  13.77727494300,   -0.03355442307 10
6793.47648813065,   -0.00183211491 11
  27.55454988600,    0.00154686936 12
  27.66676714572,    1.04662152018 13
  27.44323926226,    0.04871832111 14
1616.30271425126,    0.82350083257 15
   6.85940288609,    0.00506536922 16
   7.09580796949,    0.00053185233 17
   9.55688260180,    0.56622942023 18
   9.13295078376,    0.11822295260 19
2190.35004466729,    0.07507394684 20
  27.21222081500,    0.01181143551 21
   9.18484996200,    0.93943288590 22
  14.76531532934,   -0.28076102127 23
  31.81203816555,    0.00468638492 24
   9.54345711936,  -34.69940596576 25
  13.63341568476,   -0.02367995865 26
  26.98505934729,    2.61832577375 27
 365.24996869612,    0.00632525027 28
 182.62498434806,   -4.35727303373 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
  -0.04154440309 :K0:
   0.32272740167 :level:
   0.09945049539,    0.14233495792,    2.82293339007 1
   0.19890099078,    0.04074860781,   -1.18095774486 2
   0.29835148617,    0.07736963942,    0.22509833347 3
   0.39780198156,    0.01801059313,    0.95928086552 4
   0.49725247695,    0.09577604899,    1.01439357459 5
   0.59670297234,    0.06282921208,   -0.62233294256 6
   0.69615346773,    0.12472337575,   -1.60605415489 7
   0.79560396312,    0.13116440049,    2.61309540850 8
   0.89505445851,    0.07699326678,    1.18684004096 9
   0.99450495390,    0.17352383322,    1.51341038983 10
   1.09395544929,    0.09750132014,    0.59655667129 11
   1.19340594468,    0.06280150681,   -1.79792613936 12
   1.29285644007,    0.13311197174,   -2.57370048704 13
   1.39230693546,    0.07227354075,   -0.75952872923 14
   1.49175743085,    0.15196195110,   -0.77599806381 15
   1.59120792624,    0.18134615860,   -2.24338876599 16
   1.69065842163,    0.16815721263,   -0.72839865759 17
   1.79010891702,    0.05043092479,    2.62748007942 18
   1.88955941241,    0.15429751798,    0.67546858462 19
   1.98900990780,    0.18528345314,    3.03613883875 20
   2.08846040319,    0.07829199998,   -1.87386005882 21
   2.18791089858,    0.07744109973,   -2.46782574971 22
   2.28736139397,    0.00467212630,    1.31626589276 23
   2.38681188936,    0.15861118711,    0.72053297557 24
   2.48626238475,    0.15835312136,   -2.98697579417 25
   2.58571288014,    0.17772235712,   -1.98248117218 26
   2.68516337553,    0.26258090863,   -2.78194228824 27
   2.78461387092,    0.11031703379,    1.66223720371 28
CC   0.0675539728   0.8460361054   2 1


Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment