def calculate_viscosity_Sutherland(T):
    """Calculates the viscosity of the air

    Parameters
    -----------
    T : float
        Temperature (K)

    Returns
    -----------
    visc : float
        viscosity of the air (kg/(m s))

    Notes
    -----------
    Acoording to [1] the limis for this function are:

    p < p_c =36 Atm (3.65 MPa)
    T < 2000 K

    According to [2] the limits for this function are:

    T < 550 K


    """
    visc_0 = 1.176 * 1e-05
    T_0 = 273.1
    b = 0.4042
    return visc_0 * (T / T_0) ** (3 / 2) * ((1 + b) / (T / T_0 + b))
