def print_stats(label, means, rms, ndigits=3):
    """Print statistics about the force coefficients."""
    def r(v):
        return round(v, ndigits=ndigits)
    print(f'\n{label}:')
    print(f'<C_T> = {r(means.ct)}, (C_T)_rms = {r(rms.ct)}')
    print(f'<C_L> = {r(means.cl)}, (C_L)_rms = {r(rms.cl)}')
    print(f'<C_Z> = {r(means.cz)}, (C_Z)_rms = {r(rms.cz)}')