대기과학지

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

보호되어 있는 글입니다.

이번에 한 내용은 이심률 구하는 작업입니다 ​ 이심률이란 원운동하는 물체가 궤도에서 벗어난 정도를 뜻합니다 ​ 이걸 지구에 적용하면 태양을 도는 공전 궤도의 변화 정도를 말합니다 ​ 먼저 이심률 공식을 설명하기 전에 Day Number 라는 숫자를 알아보겠습니다 1년은 총 365일로 구성이 되어있습니다 1월 1일을 0으로 두고 12월 31일을 365로 계산합니다 ​ 2월은 4년마다 29일이 되지만 이는 포함하지 않습니다 이심률 공식입니다 ​ 2가지가 있는데 위에 공식이 조금 더 정밀하게 값이 나옵니다 ​ 위에 식을 가지고 만든 코드입니다 fortran 90까지는 한 줄에 입력 가능한 글자수가 제한이 있기 때문에 자동으로 줄바꿈이 됩니다 그럴때는 넘어간 줄 맨 앞으로 커서를 옮겨서 백스페이스바를 한 번 눌..