def AirDensity(RH, Tc, P=101.2):
    Rd =  286.9  # Specific gas constant of dry air in J/(kg*K)    
    q  =  0.622*P/(RH*SatVapor(Tc)) # Mixing ratio in kg/kg
    Tv = (Tc + 273.15)*(1.0 + 0.61*q) # Compute virtual temperature in K
    P  *= 1000.0 # Convert pressure from kPa to Pa (= Nm/m^3 = J/m^3)
    rho_a = P/(Rd*Tv) # Compute air density in kg/
    return rho_a
dataset/ETHPy150Open LEAF-BoiseState/SPEED/Module05/PenmanMonteithEx.py/AirDensity/swapped operands of `/`