def func_SH_Erbs_4d(TH:float,
                    KT:float)->float:
    """天空日射量の推計式 Erbsモデル 4次式（KTが0.22を超えて0.80以下）
    Args:
      TH(float): 水平面全天日射量(MJ/m2)
      KT(float): 晴天指数(-)

    Returns:
      SH(float)): 水平面天空日射量(MJ/m2)
    """
    return TH * (0.9511 - 0.1604*KT + 4.388*KT**2 - 16.638*KT**3 + 12.336*KT**4)