대기과학지

이번에는 NCL을 이용해서 지도를 그리고 그 위에 기상전문을 입력시켜셔 일기 부호를 입력해보겠습니다 이번 포스팅은 기상전문에 대해 설명이 조금 필요합니다 NCL 공식 사이트에서 제공하는 전문 작성 방식입니다 영어로 되어 있어서 한글로 무슨 의미인지 뒤에 달아두겠습니다 character 0=iR-the precipitation data indicator #강수 자료 표시 character 1=iX-weather data and station type indicator #일기 자료 표시 character 2=h-height above ground of base of lowest cloud #운저고도 characters 3-4=VV-visibility in miles and fractions #시정 chara..

이번에는 NCL로 그린 WRF 결과를 파이썬으로 gif로 만들어보겠습니다 혼자 해본 그림중에 3일 분량으로 21km * 21km로 그린 그림입니다 시간 간격은 3시간 간격입니다 from PIL import Image as image import numpy as np cd /content/drive/MyDrive/사진/png from IPython.display import Image al = [] for i in np.arange(1, 10, 2): a1 = "Surface_multi.00000" + str(i) + ".png" al.append(a1) for j in np.arange(11, 35, 2): a2 = "Surface_multi.0000" + str(j) + ".png" al.append(..

이번에는 WRF 모델의 결과 기록을 1시간 단위로 나오게 해보겠습니다 저번 포스팅들은 전부 WRF의 결과를 3시간 단위로 기록을 했었습니다 그래서 이번에는 1시간 단위로 줄여서 해보겠습니다 자료의 기간을 기존에는 2016년 10월 6일 00:00:00부터 2016년 10월 8일 00:00:00까지였지만 https://rda.ucar.edu/datasets/ds083.2/ NCAR RDA Dataset ds083.2 Abstract: These NCEP FNL (Final) Operational Global Analysis data are on 1-degree by 1-degree grids prepared operationally every six hours. This product is from the..

이번에는 이전 포스팅과 같은 그림이지만 WRF 결과를 보면 3시간 단위로 2016년 10월 6일 00:00:00 부터 10월 8일 00:00:00까지 총 17개의 파일이 있습니다 아 파일들을 반복문을 이용해서 한 번에 그려지도록 해보겠습니다 begin DATADir = "/home/lsh/" FILES = systemfunc("ls -1 /home/lsh/wrfout*") ;ls -1 로 해당 디렉토리 안에 wrfout이 들어간 파일들을 전부 변수로 지정 numFILES = dimsizes(FILES) print("numFILES =" + numFILES) print(FILES) print(" ") type = "png" wks = gsn_open_wks(type, "Surface_multi") res ..

이번에는 저번 포스팅에서 이슬점 온도와 풍속 풍향 그림의 결과를 보여드리겠습니다 코드는 다른게 없으니 코드와 결과 사진만 올려드리겠습니다 begin a = addfile("./wrfout_d01_2016-10-06_00.nc", "r") type = "png" wks = gsn_open_wks(type, "Surface") res = True res@MainTitle = "REAL-TIME WRF" pltres = True mpres = True mpres@mpDataBaseVersion = "MediumRes" mpres@mpDataResolution = "FinestResolution" mpres@mpDataSetName = "Earth..4" mpres@mpGridAndLimbOn = True mp..

이번에는 조금더 복잡한 그림을 그려보겠습니다 이번 그림은 변수가 하나가 아니라 기온, 풍향, 풍속, 기압 이렇게 4가지 변수를 한 번에 그렸습니다 이슬점온도와 풍향, 풍속을 그린 코드도 함께 있지만 그 부분은 다음 포스팅에서 다루겠습니다 begin a = addfile("./wrfout_d01_2016-10-06_00.nc", "r") type = "png" wks = gsn_open_wks(type, "Surface") res = True res@MainTitle = "REAL-TIME WRF" pltres = True mpres = True mpres@mpDataBaseVersion = "MediumRes" mpres@mpDataResolution = "FinestResolution" mpres@mp..

이번에는 저번 포스팅에서 나왔던 결과를 시각화 해보겠습니다 시작하기 전에 저번과는 다른게 있습니다 저번에는 예제 그대로 미국 허리케인 matthew로 진행했지만 컴퓨터 사양 문제로 격자를 크게 잡을 수 없고 크게 잡는거랑은 별개로 육지의 비중이 너무 작아서 도메인을 조금 수정해서 한반도에서 진행했습니다 이전 결과는 이렇게 출력이 됩니다 수정한 도메인을 보여드리겠습니다 &geogrid e_we = 120, e_sn = 120, ref_lat = 37.541, ref_lon = 126.986, truelat1 = 35.0, truelat2 = 60.0, stand_lon = 127.0, 다른 부분은 바뀐게 없고 이렇게 7개만 변경했습니다 ref_lat, ref_lon의 위경도는 서울시청을 기준으로 했습니다 ..

이번에는 마지막으로 WRF를 실행해보겠습니다 이제 작업 디렉토리를 /home/lsh/WRF/WRF-4.1.2/test/em_real 로 옮기겠습니다 앞에서 WRF를 컴파일 할 때 em_real로 했기 때문에 test폴더 들 중 em_real로 가주셔야합니다 먼저 실행하기전에 링크라는 것을 해주어야합니다 ln -sf ../../../WPS/met_em.d01.2016-10* 온라인 튜토리얼에서는 이렇게 하라고 나와있었는데 저는 이렇게 하면 오류가 나서 각각 개별 파일마다 링크를 했습니다 ln -sf /home/lsh/WRF/WPS-4.1/met_em.d01.2016-10-06_00:00:00.nc ln -sf /home/lsh/WRF/WPS-4.1/met_em.d01.2016-10-06_06:00:00.n..

오늘은 WRF를 직접 실행해보겠습니다 https://www2.mmm.ucar.edu/wrf/OnLineTutorial/CASES/SingleDomain/ungrib.php ARW Online Tutorial You may have noticed that the SST data are only available once a day, but we set interval_seconds to 6 hourly. By doing this, ungrib will temporally interpolate the available SST data to a frequency of 6 hourly. Normally there is no advantage in temporally www2.mmm.ucar.edu 제가 사용한 ..

이번 포스팅은 WPS를 설치해보겠습니다 설치에 앞서 이번 포스팅은 꽤 의미 있는 포스팅입니다 제 블로그에 100번째 글로 올라가는 포스팅입니다 몇 개 안쓴거 같지만 벌써 100번째라는게 놀랍네요 ㅠㅠ 앞으로도 대학원을 가서도 열심히 남겨보겠습니다! cd $HOME/WRF/Downloads wget -c https://github.com/wrf-model/WPS/archive/v4.1.tar.gz tar -xvzf v4.1.tar.gz -C $HOME/WRF cd $HOME/WRF/WPS-4.1 ./configure #3 ./compile WPS도 코드는 간단합니다 다 넘어가고 역시 configure부터 보겠습니다 Will use NETCDF in dir: /home/lsh/WRF/Library Using..