pukpr/LTE_Partial.adb

## LTE_Partial.adb
         if MLR_On or not (Forcing_Only or Is_Minimum_Entropy) then -- MLR_On
              if Climate_Trend then
                 Secular_Trend := 1.0;
              else
                 Secular_Trend := 0.0;
              end if;
              if Derivative then
                  DR(Data_Records'First).Value := Data_Records(Data_Records'First).Value;
                  for I in Data_Records'First+1 .. Data_Records'Last loop
                     DR(I).Value := DR(I-1).Value + Data_Records(I).Value;
                  end loop;
              else
                  if Partial then
                     DR(Data_Records'First) := Data_Records(Data_Records'First);
                     for I in Data_Records'First+1 .. Data_Records'Last loop
                        DR(I).Value := Data_Records(I-1).Value + D.B.ShiftT*(Data_Records(I).Value - Data_Records(I-1).Value);
                     end loop;
                  else
                     DR := Data_Records;
                  end if;
              end if;
              DR := Annual_Add(DR, -1.0);

              Regression_Factors (Data_Records => Excluded(DR), -- Time series
                                  --First => First,
                                  --Last => Last,  -- Training Interval
                                  Forcing => Excluded(Forcing),  -- Value @ Time
                                  NM => NM + NH, -- # modulations
                                  DBLT => M, --D.B.LT,
                                  DALTAP => MAP, --D.A.LTAP,
                                  DALEVEL => D.A.LEVEL,
                                  DAK0 => D.A.K0,
                                  Secular_Trend => Secular_Trend,
                                  Accel => Accel,
                                  Singular => Singular
                                 );
         else
            Singular := False;
         end if;

      CorrCoeff := 0.0;
      Pareto_Scale := 1.0;
      if Pareto then
            Pareto_Start := 1;
      else
            Pareto_Start := NH+1;
      end if;
      for Pareto_Index in Pareto_Start .. NH+1 loop
         if Forcing_Only or (Is_Minimum_Entropy and not MLR_On ) then
            Model := Forcing;
         else

            Model := LTE(Forcing => Forcing,
                         Wave_Numbers => M(1 .. NM - 1 + Pareto_Index), --D.B.Lt(1..NM),
                         Amp_Phase => MAP(1 .. NM - 1 + Pareto_Index), --D.A.LTAP,
                         Offset => D.A.level,
                         K0 => D.A.K0,
                         Trend => Secular_Trend,
                         Accel => Accel,
                         NonLin => NonLin+D.B.IR);  -- D.LT GEM.LTE.LT0

            if Derivative then
                for I in reverse Model'First+1 .. Model'Last loop
                   Model(I).Value := Model(I).Value - Model(I-1).Value;
                end loop;
                Model(Model'First).Value := Model(Model'First+1).Value;
            end if;

            if Partial then
                for I in Model'First+1 .. Model'Last loop
                   Impulses(I).Value := Model(I).Value - Model(I-1).Value;
                end loop;
                Impulses(Model'First).Value := 0.0;
                for I in Model'Range loop
                   Model(I).Value := Model(I).Value + D.B.ShiftT* Impulses(I).Value;
                end loop;
            end if;

            Model := Annual_Add(Model);

         end if;
	if MLR_On or not (Forcing_Only or Is_Minimum_Entropy) then -- MLR_On
	if Climate_Trend then
	Secular_Trend := 1.0;
	else
	Secular_Trend := 0.0;
	end if;
	if Derivative then
	DR(Data_Records'First).Value := Data_Records(Data_Records'First).Value;
	for I in Data_Records'First+1 .. Data_Records'Last loop
	DR(I).Value := DR(I-1).Value + Data_Records(I).Value;
	end loop;
	else
	if Partial then
	DR(Data_Records'First) := Data_Records(Data_Records'First);
	for I in Data_Records'First+1 .. Data_Records'Last loop
	DR(I).Value := Data_Records(I-1).Value + D.B.ShiftT*(Data_Records(I).Value - Data_Records(I-1).Value);
	end loop;
	else
	DR := Data_Records;
	end if;
	end if;
	DR := Annual_Add(DR, -1.0);

	Regression_Factors (Data_Records => Excluded(DR), -- Time series
	--First => First,
	--Last => Last, -- Training Interval
	Forcing => Excluded(Forcing), -- Value @ Time
	NM => NM + NH, -- # modulations
	DBLT => M, --D.B.LT,
	DALTAP => MAP, --D.A.LTAP,
	DALEVEL => D.A.LEVEL,
	DAK0 => D.A.K0,
	Secular_Trend => Secular_Trend,
	Accel => Accel,
	Singular => Singular
	);
	else
	Singular := False;
	end if;

	CorrCoeff := 0.0;
	Pareto_Scale := 1.0;
	if Pareto then
	Pareto_Start := 1;
	else
	Pareto_Start := NH+1;
	end if;
	for Pareto_Index in Pareto_Start .. NH+1 loop
	if Forcing_Only or (Is_Minimum_Entropy and not MLR_On ) then
	Model := Forcing;
	else

	Model := LTE(Forcing => Forcing,
	Wave_Numbers => M(1 .. NM - 1 + Pareto_Index), --D.B.Lt(1..NM),
	Amp_Phase => MAP(1 .. NM - 1 + Pareto_Index), --D.A.LTAP,
	Offset => D.A.level,
	K0 => D.A.K0,
	Trend => Secular_Trend,
	Accel => Accel,
	NonLin => NonLin+D.B.IR); -- D.LT GEM.LTE.LT0

	if Derivative then
	for I in reverse Model'First+1 .. Model'Last loop
	Model(I).Value := Model(I).Value - Model(I-1).Value;
	end loop;
	Model(Model'First).Value := Model(Model'First+1).Value;
	end if;

	if Partial then
	for I in Model'First+1 .. Model'Last loop
	Impulses(I).Value := Model(I).Value - Model(I-1).Value;
	end loop;
	Impulses(Model'First).Value := 0.0;
	for I in Model'Range loop
	Model(I).Value := Model(I).Value + D.B.ShiftT* Impulses(I).Value;
	end loop;
	end if;

	Model := Annual_Add(Model);

	end if;