pukpr/Tide_sum_minus_offset.adb

## Tide_sum_minus_offset.adb

   -- Conventional tidal series summation or superposition of cycles
   function Tide_Sum_Diff (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;
      Cumulative : Long_Float := 0.0;
   begin
      for I in Template'Range loop
         Time := Template(I).Date + Ref_Time;
         declare
            TF1, TF2 : Long_Float := 0.0;
            use Ada.Numerics.Long_Elementary_Functions;
         begin
            for J in Constituents'First .. Constituents'First + 3 loop
               declare
                  L : Amp_Phase renames Constituents(J);
                  Period : Long_Float := Periods(J);
               begin
                  if Cos_Phase then
                     TF1 := TF1 + L.Amplitude*(Cos(2.0*Pi*Year_Len/Period*Time + L.Phase));
                  else
                     TF1 := TF1 + L.Amplitude*(Sin(2.0*Pi*Year_Len/Period*Time + L.Phase));
                  end if;
               end;
            end loop;
            for J in Constituents'First + 4 .. Constituents'Last loop
               declare
                  L : Amp_Phase renames Constituents(J);
                  Period : Long_Float := Periods(J);
               begin
                  if Cos_Phase then
                     TF2 := TF2 + L.Amplitude*(Cos(2.0*Pi*Year_Len/Period*Time + L.Phase));
                  else
                     TF2 := TF2 + L.Amplitude*(Sin(2.0*Pi*Year_Len/Period*Time + L.Phase));
                  end if;
               end;
            end loop;
            Cumulative := Cumulative + cos(Scaling + TF1);
            Res(I) := (Time, cos(Scaling + TF1)*(Integ + TF2) );
         end;
      end loop;
      Cumulative := Cumulative/Long_Float(Template'Length);
      for I in Template'Range loop
         Res(I).Value := Res(I).Value - Cumulative*Integ;
      end loop;
     return Res;
   end Tide_Sum_Diff;

   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
      --Res2,
      Res1 : Data_Pairs := Template;
   begin
      Res1 := Tide_Sum_Diff (Template, Constituents, Periods, 0.0, Scaling, Cos_Phase, Year_Len, Integ);
      --Res2 := Tide_Sum_Diff (Template, Constituents, Periods, Ref_Time, Scaling, Cos_Phase, Year_Len, Integ);
      --for I in Res1'Range loop
      --   Res1(I).Value := Res1(I).Value - Res2(I).Value;
      --end loop;
      return Res1;
   end Tide_Sum;

	-- Conventional tidal series summation or superposition of cycles
	function Tide_Sum_Diff (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;
	Cumulative : Long_Float := 0.0;
	begin
	for I in Template'Range loop
	Time := Template(I).Date + Ref_Time;
	declare
	TF1, TF2 : Long_Float := 0.0;
	use Ada.Numerics.Long_Elementary_Functions;
	begin
	for J in Constituents'First .. Constituents'First + 3 loop
	declare
	L : Amp_Phase renames Constituents(J);
	Period : Long_Float := Periods(J);
	begin
	if Cos_Phase then
	TF1 := TF1 + L.Amplitude(Cos(2.0PiYear_Len/PeriodTime + L.Phase));
	else
	TF1 := TF1 + L.Amplitude(Sin(2.0PiYear_Len/PeriodTime + L.Phase));
	end if;
	end;
	end loop;
	for J in Constituents'First + 4 .. Constituents'Last loop
	declare
	L : Amp_Phase renames Constituents(J);
	Period : Long_Float := Periods(J);
	begin
	if Cos_Phase then
	TF2 := TF2 + L.Amplitude(Cos(2.0PiYear_Len/PeriodTime + L.Phase));
	else
	TF2 := TF2 + L.Amplitude(Sin(2.0PiYear_Len/PeriodTime + L.Phase));
	end if;
	end;
	end loop;
	Cumulative := Cumulative + cos(Scaling + TF1);
	Res(I) := (Time, cos(Scaling + TF1)*(Integ + TF2) );
	end;
	end loop;
	Cumulative := Cumulative/Long_Float(Template'Length);
	for I in Template'Range loop
	Res(I).Value := Res(I).Value - Cumulative*Integ;
	end loop;
	return Res;
	end Tide_Sum_Diff;

	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
	--Res2,
	Res1 : Data_Pairs := Template;
	begin
	Res1 := Tide_Sum_Diff (Template, Constituents, Periods, 0.0, Scaling, Cos_Phase, Year_Len, Integ);
	--Res2 := Tide_Sum_Diff (Template, Constituents, Periods, Ref_Time, Scaling, Cos_Phase, Year_Len, Integ);
	--for I in Res1'Range loop
	-- Res1(I).Value := Res1(I).Value - Res2(I).Value;
	--end loop;
	return Res1;
	end Tide_Sum;