function Initial_powerBased_powerHistory
input SI.Power[:, 2] history "Power history up to simulation time=0, [t,Q]";
input TRANSFORM.Units.InverseTime[:] lambdas "Decay constants for each precursor group";
input TRANSFORM.Units.NonDim[size(lambdas, 1)] alphas "Normalized precursor fractions [betas = alphas*Beta]";
input TRANSFORM.Units.NonDim Beta "Effective delayed neutron fraction";
input SI.Time Lambda "Prompt neutron generation time";
input TRANSFORM.Units.InverseTime[:] lambdas_dh "Decay constants for each group";
input Units.NonDim efs[size(lambdas_dh, 1)] "Decay-heat fraction of fission power";
input Boolean includeDH = false "=true if power history includes decay heat";
input SI.Power[size(lambdas, 1)] Cs_0 = fill(0, size(lambdas, 1)) "Precursor concentration at history time = 0";
input SI.Energy[size(lambdas_dh, 1)] Es_0 = fill(0, size(lambdas_dh, 1)) "Decay-heat concentration at history time = 0";
output SI.Power[size(lambdas, 1)] Cs;
output SI.Energy[size(lambdas_dh, 1)] Es;
end Initial_powerBased_powerHistory;