def func_TH(P:float,
            IN0:float,
            Sinh:float,
            SH:float)->float:
    """水平面全天日射量の式
    Args:
      P(float): 大気透過率(-)      
      IN0(float): 大気外法線面日射量(MJ/m2)
      Sinh(float): 太陽高度角のサイン(-)
      SH(float): 水平面天空日射量(MJ/m2)

    Returns:
      TH(float)): 水平面全天日射量(MJ/m2)
    """
    return IN0 * P**(1/Sinh) * Sinh + SH