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;