저번 포스팅에 이어서 대기 상단 일사량을 계산하고 시각화 해보겠습니다 저희가 해야할 건 이 공식을 계산하는 코드를 구현하면 됩니다 대기 상단 일사량은 시간각하고 똑같게 위도별로 365개를 계산해줘야합니다 과제에는 위도별 평균 일사량과 적도, 북위 30, 60, 90도에서의 그래프 이렇게 2가지를 요구하고 있기에 일반적으로 계산된 값이 나열된 자료와 평균값 이렇게 2개가 필요합니다 바로 가보겠습니다 def cal_e(phi): ii = [] s = 1361 d1 = de d2 = be for i in range(len(d1)): E = (s/np.pi) * d2[i] * (globals()['a{}'.format(phi)][i] * np.sin(np.radians(phi)) * np.sin(np.radia..
