import math


def R_coat(delta_coat, k_coat, radius, Theta):
    """ thermal resistance due to conduction through coating layer
    
    Parameters
    ----------
    k_coat:     float
                thermal conductivity of the coating in W/(mK)
    delta_coat: float
                thickness of the coating in m
    radius:     float
                radius of drop in m
    Theta:      float    
                static contact angle in deg                   
    
    Returns
    ----------
    R_coat:     float
                thermal resistance due to conduction through coating layer in K/W
    """
    R_coat = delta_coat / (k_coat*math.pi*radius**2*(math.sin(math.radians(Theta)))**2)
    return R_coat