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 16, 2023

PDO odd harmonics

This is a set of odd harmonics that has a bug which multiplies the harmonics when saved and reloaded

image

0.00114616290 :offset:
0.00000000000 :bg:
0.00000000000 :impA:
0.00000000000 :impB:
0.00000000000 :impC:
-9.17224981172:delA:
0.22493884587:delB:
9.17248760826:asym:
0.09973887381:ann1:
-0.09646558041:sem1:
1.36226903406:sem2:
0.00007647650:year:
0.00538081717:IR:
0.00000890175 :mA:
0.00077698555 :mP:
0.00852871790 :shiftT:
-0.81788767220 :init:
---- Tidal ----
9.10846048884, 0.01752333734 1
29.53063068364, -0.02002212147 2
27.09267692660, 0.34985859282 3
6.85248390316, 0.09756523556 5
13.60611040750, 1.02746072649 6
13.66083077700, 0.13221326457 7
3396.73824406533, 1.42391150261 8
9.12068919638, 0.00542846536 9
13.77727494300, -0.01467729158 10
6793.47648813065, -0.00290652324 11
27.55454988600, 0.00151101429 12
27.66676714572, 0.97038310177 13
27.44323926226, 0.03796950012 14
1616.30271425126, 0.82446003780 15
6.85940288609, 0.00468342807 16
7.09580797237, 0.00052852283 17
9.55688260703, 0.56404751589 18
9.13295078376, 0.12099760718 19
2190.35004466729, 0.10577973110 20
27.21222081500, 0.01072397926 21
9.18484996200, 0.93951687242 22
14.76531534182, -0.21590169054 23
31.81203822350, 0.00363248539 24
9.54345712458, -11.23658033725 25
13.63341568476, -0.01408683466 26
26.98505934729, 2.51127206939 27
365.24996487650, 0.01144691421 28
182.62498243825, -4.97404975554 29
---- LTE ----
0.00000000000 :trend:
0.00000000000 :accel:
-0.01046520108 :K0:
0.06301108998 :level:
-0.08152632502, 0.07299352634, -0.71548622404 1
-3.66868462591, 0.35146328825, 2.76947311292 45
-0.24457897506, 0.24187970722, -1.06435274701 3
-2.20121077555, 0.48258023103, 1.05679118245 27
-2.69036872567, 0.42561621781, -1.55095510429 33
-1.38594752534, 0.22586338102, -0.41216858503 17
-3.99478992599, 0.28123552370, -2.58792170326 49
-4.64700052615, 0.09712888747, 2.41436351587 57
-5.29921112631, 0.15466050829, -2.88877956546 65
-4.32089522607, 0.13084753256, -2.78074667879 53
-6.60363232664, 0.22505873688, -2.59933176718 81
-7.25584292680, 0.05374704802, 0.41306254568 89
-2.52731607563, 0.13733955568, -2.20417350582 31
-8.56026412712, 0.17582788528, -1.45464260261 105
-2.03815812551, 0.08876777616, 0.92946854038 25
-9.86468532745, 0.16740216348, -3.02924210343 121
-10.51689592761, 0.03723743249, -1.63686424431 129
-6.60470046958 -3.34257932583, 0.22582343838, 2.10205688508 41
-11.82131712793, 0.02407856892, 0.27908796988, 0.24918698882, 1.05520504094 81
145
-7.25701656534, -5.95142172648, 0.07572121195, 0.05047891244, -0.18216166912 89
1.63959209990 73
-7.90933266110, 0.10050573759, -1.85090136778 97
-8.56164875686, 0.14684389680, 2.20693866305 105
-1.06001365561, 0.17600482605, 0.60998810796 13
-2.36464584713, 0.14875931698, -0.23108334356 29
-10.51859704414, 0.13573195681, -1.23313586723 129
-3.34311999077, 0.22756406739, 0.94721958610 41
-11.82322923566, 0.10867949332, 0.21487995729 145
-1.22309267955, 0.02654897021, -2.13647511113 15
-0.24461853591, 0.18172146780, -1.29298901032 3
-1.38617170349, 0.23304807719, -0.84135108879 17
-8.88780680474, 0.23159997155, -0.76254821707 109
-0.40769755985, 0.13351676655, -2.64420287900 5
-1.87540877531, 0.22673736821, -3.13429787637 23
-16.38944190599, 0.09940980237, -1.03327575324 201
-4.32159413441, 0.17245392866, 2.47061137051 53
-2.20156682319, 0.42907696545, 0.35666586304 27
-13.12573832825, 0.18893497736, 1.48911290437 161
-13.77794892842, 0.15633771139, -1.35848101623 169
-8.88636942720, 0.18908809438, 1.44686234963 109
-0.40763162510, 0.12227231760, -2.39601729576 5
-0.73373692518, 0.10680768601, 2.60956197785 9
-16.38679132906, 0.22663094873, -2.54234548196 201
-2.36426342559, 0.09671592509, 1.16278441680 29
-1.87510547546, 0.18495463229, -2.32352376132 23
CC 0.5041092212 0.8464354657 3 1
0.00000000000:dLOD:CC 0.0550207915 0.8464432213 1 1
0.00000000000:dLOD:

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

@pukpr
Copy link
Author

pukpr commented Jun 17, 2023

AMO annual

image

  0.00112098511 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.35083453519:delA:
   0.27681578514:delB:
   0.00000000000:asym:
   0.03403812760:ann1:
  -3.01729778071:ann2:
  -0.00994197148:sem1:
   1.14606692255:sem2:
   0.00007553123:year:
   0.08168885333:IR:
   0.00001256002 :mA:
   0.00070828290 :mP:
   0.00855636784 :shiftT:
  -0.43369020857 :init:
---- Tidal ----
   9.10846048884,    0.01747994186 1
  29.53063068982,   -0.02316971812 2
  27.09267692660,    0.35410170786 3
   7.08840411226,    1.01122625285 4
   6.85248390316,    0.09935469410 5
  13.60611040750,    1.02815474586 6
  13.66083077700,    0.13276177492 7
3396.73824406533,    1.49613244952 8
   9.12068919638,    0.00589034986 9
  13.77727494300,   -0.01746610710 10
6793.47648813065,   -0.00266811744 11
  27.55454988600,    0.00142840434 12
  27.66676714572,    0.97261064484 13
  27.44323926226,    0.03369991542 14
1616.30271425126,    0.82907339430 15
   6.85940288609,    0.00562422194 16
   7.09580797308,    0.00050343830 17
   9.55688260833,    0.56581892268 18
   9.13295078376,    0.12139797627 19
2190.35004466729,    0.06217722682 20
  27.21222081500,    0.01024142276 21
   9.18484996200,    0.93905066261 22
  14.76531534491,   -0.22730385748 23
  31.81203823784,    0.00204794799 24
   9.54345712587,  -27.77146657826 25
  13.63341568476,   -0.01372217132 26
  26.98505934729,    2.35257905076 27
 365.24996393123,    0.01138571821 28
 182.62498196561,   -5.99502745787 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
-424.55945986242 :K0:
-2238.97991847792 :level:
   0.06374847351, 2110.70013385188,    2.11236197280 1
   0.12749694702, 1041.54259522066,    1.08303136059 2
   0.19124542054,  679.17079923676,    0.05363982615 3
   0.25499389405,  493.67726700503,   -0.97572883459 4
   0.31874236756,  379.12843205949,   -2.00500219494 5
   0.38249084107,  300.48547403737,   -3.03440719269 6
   0.44623931458,  242.66021511678,    2.21936881828 7
   0.50998778810,  198.01337895480,    1.18955896613 8
   0.57373626161,  162.45417347215,    0.15849622500 9
   0.63748473512,  133.51563609626,   -0.87225830614 10
   0.70123320863,  109.60765904114,   -1.90330447767 11
   0.76498168214,   89.78181859692,   -2.93537727002 12
   0.82873015566,   73.27065794103,    2.31579981411 13
   0.89247862917,   59.46465456039,    1.28233127793 14
   0.95622710268,   47.93624232939,    0.24837550015 15
   1.01997557619,   38.41093958673,   -0.78893699402 16
   1.08372404970,   30.38682958023,   -1.82906701959 17
   1.14747252322,   23.78693407049,   -2.87085363863 18
   1.21122099673,   18.40319113865,    2.37204208875 19
   1.27496947024,   13.86517840250,    1.33442512271 20
   1.33871794375,   10.35876683819,    0.29331868416 21
   1.40246641726,    7.57618662425,   -0.74678566603 22
   1.46621489078,    5.22497776841,   -1.79395745011 23
   1.52996336429,    3.45942460194,   -2.84058923452 24
   1.59371183780,    2.17754616777,    2.40922255339 25
   1.65746031131,    1.25365724743,    1.35198597544 26
   1.72120878482,    0.63654046014,    0.24293185203 27
   1.78495725834,    0.22018402397,   -0.92481860598 28
CC   0.1026013846   0.8008152775   6 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> $env:CC_START
1880
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\backup.bat amo_annual_1880_1932

image

   0.00111285555 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.56714148611:delA:
   0.27681578514:delB:
   0.00000000000:asym:
   0.02627188835:ann1:
  -3.01933241452:ann2:
  -0.01541132495:sem1:
   2.15945145374:sem2:
   0.00007080127:year:
   0.05481033024:IR:
   0.00000985933 :mA:
  -0.00071327720 :mP:
   0.00856661400 :shiftT:
   0.64473092418 :init:
---- Tidal ----
   9.10846048884,    0.01751524502 1
  29.53063072074,   -0.02884814349 2
  27.09267692660,    0.35462631489 3
   7.08840411582,    1.01131217622 4
   6.85248390316,    0.09941748308 5
  13.60611040750,    1.02809933479 6
  13.66083077700,    0.13696309407 7
3396.73824406533,    1.49609420887 8
   9.12068919638,    0.00578893533 9
  13.77727494300,   -0.01778767674 10
6793.47648813065,   -0.00264868799 11
  27.55454988600,    0.00149194998 12
  27.66676714572,    0.93495064500 13
  27.44323926226,    0.03093884606 14
1616.30271425126,    0.99969287678 15
   6.85940288609,    0.00821180617 16
   7.09580797665,    0.00047907795 17
   9.55688261480,    0.57527880621 18
   9.13295078376,    0.13001592702 19
2190.35004466729,    0.06476519571 20
  27.21222081500,    0.01167179304 21
   9.18484996200,    0.71568344120 22
  14.76531536037,   -0.23678089561 23
  31.81203830960,    0.00231517806 24
   9.54345713232,  -27.56095478406 25
  13.63341568476,   -0.01390588229 26
  26.98505934729,    2.49103430935 27
 365.24995920127,    0.01134600205 28
 182.62497960063,   -5.99111354263 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
   0.09139539678 :K0:
   0.38022412009 :level:
 505.64442648742,    0.04963310857,    2.60722431591 7489
   0.06751445367,    0.45640869402,   -1.55758990631 1
   0.13502890734,    0.18540665283,   -2.86189114977 2
   0.20254336100,    0.08303687090,    2.04582650985 3
   0.94520235136,    0.06258102157,    2.77818554941 14
   0.33757226834,    0.09977350148,    0.35996310938 5
   0.40508672201,    0.09462530431,   -0.90032828925 6
   0.47260117568,    0.09269064491,   -2.55665589595 7
   0.54011562935,    0.07479091331,    1.38433108973 8
   2.16046251739,    0.02733741595,   -0.22628751306 32
   0.67514453668,    0.05584658032,    1.77755229776 10
   3.78080940542,    0.04528615027,    0.66171774142 56
   2.22797697105,    0.03330520985,   -2.14985783482 33
   0.27005781467,    0.05641360080,    1.70666586485 4
   3.71329495176,    0.01709463889,    2.38456460861 55
   1.01271680502,    0.05782849162,    0.64382257749 15
   2.97063596140,    0.01444049767,    1.87129204279 44
   1.14774571236,    0.04529977445,   -3.05389462875 17
   1.21526016603,    0.09192528606,    1.77071607287 18
   1.28277461970,    0.05655968770,   -0.98833788800 19
   2.43052033206,    0.01612491733,   -0.08480607274 36
   1.41780352703,    0.08871973253,   -1.55538214546 21
   1.48531798070,    0.05953374447,    3.02711342791 22
   1.55283243437,    0.03136344878,    1.08678624307 23
   1.62034688804,    0.03108223818,   -1.57046051670 24
   3.57826604442,    0.01881657722,   -2.66942878570 53
   3.10566486874,    0.02050041211,    2.02894762309 46
   1.82289024904,    0.07458943198,   -3.02898667178 27
   1.89040470271,    0.07068548578,    0.99699652738 28
CC  -0.4155325601   0.8656196543   6 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\backup.bat amo_annual_1880_1932c

@pukpr
Copy link
Author

pukpr commented Jun 17, 2023

Darwin from NINO4

image

   0.00109679653 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
  -9.07759355841:delA:
   0.28482498794:delB:
   9.17098647185:asym:
   0.09218781368:ann1:
  -2.80295183685:ann2:
  -0.16390860780:sem1:
   0.93380332335:sem2:
   0.00008001455:year:
   0.09120358694:IR:
   0.00001423192 :mA:
  -0.00084014605 :mP:
   0.00852138525 :shiftT:
  -0.69506096917 :init:
---- Tidal ----
   9.10846048884,    0.01780193986 1
  29.53063066051,   -0.02618876943 2
  27.09267692660,    0.34403047539 3
   7.08840410888,    1.04863253062 4
   6.85248390316,    0.09729712617 5
  13.60611040750,    1.02538228569 6
  13.66083077700,    0.12733805414 7
3396.73824406533,    1.42389083326 8
   9.12068919638,    0.00732280496 9
  13.77727494300,   -0.06385257650 10
6793.47648813065,   -0.00192941472 11
  27.55454988600,    0.00156193507 12
  27.66676714572,    1.32989278112 13
  27.44323926226,    0.04885687199 14
1616.30271425126,    1.52041696912 15
   6.85940288609,    0.00450355767 16
   7.09580796970,    0.00054845036 17
   9.55688260219,    0.57647055022 18
   9.13295078376,    0.11908277253 19
2190.35004466729,    0.07470718227 20
  27.21222081500,    0.01154428152 21
   9.18484996200,    1.20039667414 22
  14.76531533026,   -0.28412560903 23
  31.81203816982,    0.00293250235 24
   9.54345711974,  -33.70230950695 25
  13.63341568476,   -0.02491031318 26
  26.98505934729,    2.74038381957 27
 365.24996841455,    0.00687047953 28
 182.62498420727,   -1.07028516539 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
  -0.04221279836 :K0:
   0.33645601493 :level:
   0.09951886453,    0.23217948415,    2.72420566404 1
   1.89085842611,    0.19636738212,    0.50427761983 19
   0.29855659360,    0.09479244603,   -0.97468061231 3
   0.59711318719,    0.04847954656,   -1.12114414247 6
   0.49759432266,    0.10415069810,    1.07704123512 5
   4.47834890394,    0.06329167970,    2.20201900253 45
   5.17498095567,    0.14017565191,    0.59761840829 52
   1.19422637439,    0.07654797875,   -2.53395583517 12
   0.79615091626,    0.07025452438,    2.19760817153 8
   0.99518864532,    0.20691321757,    1.43247373996 10
   0.39807545813,    0.05091239301,    0.67840525244 4
   3.28412252956,    0.17681475346,   -0.28230042718 33
   2.68700934237,    0.20528399214,   -2.51079463751 27
   2.58749047783,    0.13853096643,   -2.27462053169 26
   2.88604707143,    0.13992373151,    1.16178736094 29
   1.59230183251,    0.16511269771,   -2.02171682590 16
   2.48797161330,    0.17790728294,   -2.69178985400 25
   1.39326410345,    0.10465806136,   -0.87726726573 14
   3.38364139409,    0.12373966407,   -1.96081745121 34
   5.37401868473,    0.14467388742,   -2.94952323220 54
   1.69182069705,    0.10928181704,   -1.20321919156 17
   3.98075458128,    0.15081949216,    1.76945843383 40
   4.17979231035,    0.09365039494,    0.07539475915 42
   2.38845274877,    0.16769231179,    0.95164309536 24
   2.18941501971,    0.09671047843,   -1.20065157167 22
   1.99037729064,    0.23228402962,   -3.01957974476 20
   4.37883003941,    0.10898029597,    2.19626022871 44
   1.49278296798,    0.09604388330,   -0.49569266883 15
CC   0.0905655313   0.8910796732   3 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\backup.bat  darwin_constituents2crossvalid3bhigh

Tahiti -- reversed the impulse

image

   0.00111706818 :offset:
   0.00000000000 :bg:
   0.00000000000 :impA:
   0.00000000000 :impB:
   0.00000000000 :impC:
   9.13892040543:delA:
   0.27665632225:delB:
   9.17519751825:asym:
   0.10027624074:ann1:
   4.55747889060:ann2:
   0.03528371114:sem1:
  -0.49024695247:sem2:
   0.00007277788:year:
   0.12330662698:IR:
   0.00001516005 :mA:
   0.00066746186 :mP:
   0.00851180548 :shiftT:
  -1.17659597372 :init:
---- Tidal ----
   9.10846048884,    0.01970808414 1
  29.53063070782,   -0.02203620164 2
  27.09267692660,    0.33447176700 3
   7.08840411433,    1.01125247893 4
   6.85248390316,    0.09469597033 5
  13.60611040750,    1.02766984654 6
  13.66083077700,    0.13859839206 7
3396.73824406533,    1.42279510825 8
   9.12068919638,    0.00201351092 9
  13.77727494300,   -0.17054252454 10
6793.47648813065,   -0.00203828469 11
  27.55454988600,    0.00104589368 12
  27.66676714572,    0.73075200142 13
  27.44323926226,    0.04068132090 14
1616.30271425126,    0.82393126909 15
   6.85940288609,    0.00615891450 16
   7.09580797516,    0.00051625333 17
   9.55688261210,    0.53766230007 18
   9.13295078376,   -0.11484822577 19
2190.35004466729,    0.03536548818 20
  27.21222081500,    0.01262896773 21
   9.18484996200,    0.94010408306 22
  14.76531535391,   -0.19583917000 23
  31.81203827961,    0.01607464046 24
   9.54345712963,  -30.24704333995 25
  13.63341568476,    0.00007822884 26
  26.98505934729,   17.83364711651 27
 365.24996117788,    0.00723345340 28
 182.62498058894,  -11.39884305816 29
---- LTE ----
   0.00000000000 :trend:
   0.00000000000 :accel:
  -0.05533755098 :K0:
  -0.08202414309 :level:
   0.08039511787,    0.29430572700,    1.43023638980 1
   1.12553165016,    0.19739572744,    2.91882259713 14
   0.16079023574,    0.19535145077,    0.65493248880 2
   3.45699006835,    0.08932445955,    1.15411821837 43
   1.76869259311,    0.19692380843,   -1.65586049907 22
   0.40197558934,    0.06139977409,    2.64796083410 5
   0.56276582508,    0.13399723447,   -1.39344723857 7
   0.64316094295,    0.06570586215,   -3.11760248161 8
   1.92948282884,    0.06389084305,    2.96609606505 24
   3.69817542195,    0.18196535512,   -1.89291916400 46
   0.88434629655,    0.08627067225,    1.99116700569 11
   4.01975589343,    0.14479478474,    1.52498565030 50
   1.04513653229,    0.11922067069,   -2.67847283604 13
   2.49224865392,    0.04589770421,    0.41813639708 31
   3.13540959687,    0.04862964783,   -1.73348793951 39
   4.10015101129,    0.09953147249,   -0.84148993099 51
   2.97461936113,    0.08632914077,   -0.22381971517 37
   1.44711212163,    0.14875394609,    2.26935194494 18
   0.72355606082,    0.16180866948,    2.13508187894 9
   1.20592676803,    0.27017702753,    2.50231354967 15
   3.53738518621,    0.08081196350,   -1.44081105469 44
   3.21580471474,    0.08846133830,   -2.30359655524 40
   1.52750723950,    0.04725254576,   -0.48604894115 19
   1.36671700376,    0.10433480865,   -1.17267587894 17
   2.00987794671,    0.21430822796,    1.80459887492 25
   2.09027306458,    0.21769542085,    0.51808455670 26
   1.84908771098,    0.04510899274,    2.33653157074 23
   2.25106330032,    0.14386429998,   -1.74257275909 28
CC  -0.1156376618   0.7069699849   7 1
   0.00000000000:dLOD:
PS C:\Users\pp\github\pukpr\GEM2\GeoEnergyMath\eqlong> ..\..\..\GeoEnergyMath\io\backup.bat  tahiti_constituents2crossvalid3negimpulsehigh

image
image

@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