대기과학지

저번 포스팅에 이어서 대기 상단 일사량을 계산하고 시각화 해보겠습니다 저희가 해야할 건 이 공식을 계산하는 코드를 구현하면 됩니다 대기 상단 일사량은 시간각하고 똑같게 위도별로 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..

정말 오랜만에 파이썬 글을 쓰네요.. 시험 기간이기도 하고 연구실 구성하는거랑 이것 저것 일이 많아서 좀 걸렸네요 ㅠㅠ 그동안 했던 것들을 다시 하나씩 올려보려고 합니다 오늘은 대기 상단 연평균 일사량을 계산해보고 시각화해보려고 합니다 기후학이라는 과목에서 과제로 나왔습니다 다음에 기후학을 듣는 학생들도 도움이 됐으면 좋겠습니다 먼저 코드 공유전에 몇가지 공식들 설명이 필요합니다 포트란으로 적위와 이심률을 계산했던 포스팅이 있어서 보셨을만한 공식입니다 순서대로 태양의 적위, 감마, 이심률 공식입니다 포트란으로 계산해 둔 테이블이 있지만 이번에는 파이썬으로도 새로 계산해서 사용해보려고 합니다 다음은 시간각입니다 정확히는 일출 시간각을 구해서 사용합니다 이렇게 구해진 3가지 값을 이용해서 아래 공식을 계산합..