def calc_burnup_rho(values):
    """Calculate burnup and reactivity.

    This function calculates burnup, reactivity and standard deviation
    for reactivity.

    Arguments:
        values (dict)[-]: Main dictionary storing processed values.

    Returns:
        out_tup (tuple)[-]: (Burnup[MWd], Reactivity[%dk/k], rho_stddev[%dk/k]).

    """
    round_digit = 5
    mU_UWNR = 62.25  # Estimated total mass of U in UWNR.
    keff_list = values['burnup_summary']['keff']
    burnup_list = values['burnup_summary']['burnup']
    stddev_list = values['burnup_summary']['std_dev']

    bu = [round(bp*mU_UWNR, round_digit) for bp in burnup_list]
    rho = [round((k-1)/k*100.0, round_digit) for k in keff_list]
    rho_stddev = [round((std_dev/keff**2)*100.0, round_digit)
                  for keff, std_dev in zip(keff_list, stddev_list)]
    return (bu, rho, rho_stddev)