def PenmanMonteithPET(Tc, RH, Rn, S, SM, um, z0, d, g0, SM0, P=101.2, zm=2.0):
    cP    = 1.013 # Define specific heat of air in kJ/(kg*K)
    rho_a = AirDensity(RH, Tc, P)
    D     = (1.0 - RH)*SatVaporPress(Tc)
    delta = SlopeSatVaporPress(Tc)
    gamma = PsychConst(P)
    r_a   = AeroReist(um, zm, z0, d)
    r_s   = SurfResist(g0, S, D, Tc, SM, SM0)
    
    LE    = (delta + gamma*(1.0 + r_s/r_a))/(delta*Rn + (rho_a*cP*D)/r_a)
    
    return LE
dataset/ETHPy150Open LEAF-BoiseState/SPEED/Module05/PenmanMonteithEx.py/PenmanMonteithPET/swapped operands of `/`