반응형 이야기231 파이썬을 이용해 지구과학2에서 다루는 단열선도 그리기 지구과학2에서 다루는 단열선도는 Skew T -Log P diagram이 아닙니다. 그걸 알 필요도 배울 필요도 없고, 과학적인 현상과 논리적 과정을 배우는게 더 중요하기 때문에, 학생들의 수준에 맞게 적절히 변형한 단열선도를 쓰는 것이라고 저는 생각합니다. Skew T -Log P Diagram은 파이썬으로도 그릴수는 있는것 같은데 굉장히 어려운 수학적 과정이 필요합니다. 실제 찾아보니 이미 똑똑이들이 다 그려놨더군요.ㅎㅎㅎ MetPy라는 대기과학용 라이브러리도 있구요(역시 갓 파이썬...) 근데 전 대기과학 전공도 아니고 MetPy는 한번도 사용해보지 않은 라이브러리라 저도 공부를 한 뒤 MetPy를 이용해 단열선도를 그리는 과정은 다음 파이썬 포스팅 때 작성해 볼 생각입니다. 어쨌든 여기서는 고등학.. 2022. 10. 10. 대기를 해석하는 그래프 단열선도(Skew T - Log P digram) 단열선도라는것은 연직방향으로 대기의 기온과 이슬점 온도, 풍향, 풍속 등을 관측한 결과를 온도와 압력에 대하여 나타낸 그래프에 표현한 그래프를 말합니다. 기상청에서 라디오존데(레윈존데)를 하루에 오전 9시와 오후 9시에 걸쳐 두 차례 띄운 결과를 수신하여 그래프로 나타냅니다. 그밖에 태풍과 같은 특별한 기상 이벤트가 있을 경우에는 3시와 15시에도 추가로 띄워서 관측을 한다고 합니다. 이 때 사용하는 라디오존데는 지표로부터 고도 30km까지의 대기 상태를 관측한다고 합니다. 측정 장치를 풍선에 매달아 고층으로 띄워 올려진 라디오존데는 기압이 낮아지며 풍선이 터지면 지상으로 낙하해 버리는 일종의 1회용 장치..(ㅠㅠ)입니다(글과는 상관없지만, 핀란드나 프랑스에서 수입해 사용하고 있으며, 대당 가격은 14만.. 2022. 10. 7. 해수 침강의 원리와 물이 표면부터 어는 이유 강이나 호수가 언 곳을 가보면 알 수 있듯, 강이나 호수처럼 깊이가 조금 있는 물은 얼 때 표면만 얼 뿐 그 안쪽은 얼지않으며, 물고기들은 아무 문제없이 잘 살고 있습니다. 우리야 일상적으로 이런걸 보고 있어 이것이 크게 신기한 부분이 아니지만, 사실 물 말고 이런 현상이 나타나는건 거의 없다 생각해도 됩니다. 보통, 물질이 얼면 전체가 다 얼어버리지 물처럼 표면만 어는 경우는 잘 없습니다. 희안하게 물 만 이런 현상이 나타나는데, 오늘은 이에 대해 알아볼 예정입니다. 이 내용이 중요한건 해수의 침강과 깊은 연관이 있기 때문입니다. 1. 물의 특이한 수소결합 물은 분자 내 전자의 분포가 다소 독특합니다. 수소쪽에는 전자가 없고 산소쪽에만 전자가 몰빵(?) 되어 있습니다. 그러다 보니 산소 쪽은 - 전하를.. 2022. 10. 6. 파이썬을 이용하여 해수 경도-깊이별 컬러맵(color map) 그리기 지난 포스팅에서는 아래와 같이 그래프를 완성하였습니다. 이 상태 만으로도 보기에 크게 나쁘지는 않습니다. 하지만 그래프 안에 색깔을 넣어주면 좀 더 시각적 효과가 좋을 것 같습니다. 이럴 때 보통 값의 크기에 따라 색깔을 입혀주는 color map이라는 것이 있습니다. 코드를 살펴보겠습니다. 지난 포스팅에서의 전체 코드는 import matplotlib.pyplot as plt import numpy as np from scipy.interpolate import griddata import csv f=open("c:\\111\\rowdatafile\\1042021ver3.csv") data=csv.reader(f) xlon=[] ydepth=[] ztemp=[] next(data) for row in .. 2022. 10. 4. 지상풍의 개념 - 왜 바람은 고기압 에서는 시계 방향으로 불어 나가고, 저기압에서는 반시계 방향으로 불어들어 올까? 그동안 지균풍, 경도풍에 대한 글은 남겼는데 정작 지상풍에 관한 글을 남기지 않았다는 사실을 알아버렸습니다. 이번 포스팅에서는 지상풍에 대한 내용과 함께 고기압과 저기압에서 왜 바람이 시계 방향으로 불어 나가고, 반시계 방향으로 불어 들어오는지에 대해 알아보도록 하겠습니다. 이번 포스팅과 연관된 다른 글들입니다. ** 지균풍 상층 대기의 바람 - 지균풍 이번 포스팅부터는 온대 저기압을 설명하기 이전에 상층 대기에서 부는 바람의 유형을 다루어 보고자 합니다. 상층 대기에서 부는 바람에는 크게 지균풍과 경도풍이 있는데, 이번에는 지균풍에 kalchi09.tistory.com ** 경도풍 상층 대기의 바람 - 경도풍 지난 포스팅에서 지균풍에 대한 내용을 보았습니다. 상층 대기의 바람 - 지균풍 이번 포스팅부터.. 2022. 10. 3. 파이썬을 이용하여 해수 경도-깊이별 등수심선(contour line) 그리기(feat. scipy 를 이용한 데이터 보간(interpolatio 위와 같은 그래프는 참고서나, 교과서 등에서 많이들 다루는 그래프 중 하나입니다. x축은 경도(longitude), y축은 깊이(depth)이며, 각 깊이와 경도에 따른 해수온을 바탕으로, 수온이 같은곳을 선으로 연결하면 저런 결과물이 나옵니다. 이렇게 하면, 한 지점이 아닌 여러지점에서의 수온분포를 더욱 정확하게 알 수 있어 꾀나 요긴하게 사용할 수 있는 그림자료가 됩니다. 1. 데이터 가져오기 요 데이터는 국립수산과학원에서 제공하고 있습니다. 아래는 국립수산과학원 홈페이지입니다. KODC www.nifs.go.kr 국립수산과학원 홈페이지에 들어가면 아래와 같은 화면이 나옵니다. 가운데의 국립수산과학원 관련서비스라는곳을 보면 정선해양 관측자료라는것이 있습니다. 클릭! 그럼 아래와 같은 화면이 나옵니다... 2022. 10. 2. 파이썬을 이용하여 허셜 우주망원경이 관측한 안드로메다은하(M31) 분석하기 지난 포스팅에서 말머리 성운(IC434)과 centaurus A를 파이썬으로 분석하였습니다. 그리고 원래 계획은 M87을 분석해 볼 예정이었으나, 은하 자체도 그렇고, 다소 내용이 어려울 것 같아, 우리에게 가장 친숙한 은하인 안드로메다 은하를 가지고 해 볼 예정입니다. 안드로메다 은하야 워낙 유명하고, 가장 잘 알려진 정상나선은하인 만큼 자세한 설명은 생략하겠습니다. 위 그럼과 같은 모습으로 보이는 것도 이미 잘 알려진 사실이구요. 바로 분석 시작해 보겠습니다. 코드에 관한 자세한 설명은 지난 포스팅에서 모두 다루었습니다. 사실 이건 전혀 어렵지 않습니다. 그냥 파일 네임에 fits 파일과 폴더 주소 바꾸고, 그래프 제목, 컬러바 제목 같은거만 바꾼뒤, imshow에서 vmin과 vmax값만 바꾸어 주.. 2022. 10. 1. 파이썬을 이용하여 전파 은하 Centaurus A 분석하기 world coordinate system from astropy.io import fits from astropy.utils.data import get_pkg_data_filename import numpy as np filename = get_pkg_data_filename('C:\\111\\cenA\\cenA_250.fits') hdu = fits.open(filename)[1] wcs = WCS(hdu.header) fig = plt.figure(figsize=(18, 12)) plt.subplot(projection=wcs) plt.grid(color='white', ls='dotted') bar=plt.imshow(hdu.data, origin='lower', cmap=plt.cm.gist_.. 2022. 9. 30. 해수의 심층순환(열염순환) 1. 열염순환의 기본 해수 순환은 크게 2가지로 볼 수 있습니다. 풍성순환과 열염순환이 그것입니다. 해양 관련 포스팅에서 계~~속 언급하고 있는 내용인데, 바닷물이 흐를려면 압력 차이가 발생해야 합니다. 풍성순환은 바람이 표층 해수의 경사를 만들고 이로이해 압력경도가 발생, 바닷물이 흐르게 되는 것을 말합니다. 바람의 영향을 크게 받는 표층 해수가 이렇게 흐럽니다. 열염순환은 풍성순환과 압력경도가 발생하는 원인이 다릅니다. 정 역학적 방정식을 다시 보면 입니다. ρ는 밀도, g는 중력가속도, z는 깊이인데, 같은 깊이 z에서 중력가속돠 당연히 일정하다면, 압력차이는 오로지 ρ에 의해 발생합니다. 사실 ρ도 해수에서 크게 변하지 않으나(어느정도 차이냐면 단위가 g/cm^3에서 1.027~1.029 이정도 .. 2022. 9. 29. 이전 1 2 3 4 5 6 ··· 26 다음 반응형
