대기과학지

오늘은 grib2 파일을 ungrib할 때 발생하는 Data not found가 발생할 때 해결 방법을 알아보겠습니다 먼저 에러 구문 입니다 길어서 중간 부분은 생략했습니다 중간 부분은 상층 하층을 확인하는 부분으로 o x 이렇게 표시 되는 부분입니다 *** Starting program ungrib.exe *** Start date = 2020-11-28_00:00:00 , End date = 2020-11-30_00:00:00 output format is WPS ######################################################################### Inventory for date = 2020-11-28_00:00:00 PRES TT UU VV RH HG..

이번 포스팅은 파이썬 라이브러리인 WRF-PYTHON을 사용해서 제가 만든 WRF 결과물을 시각화해보겠습니다 먼저 코드 입니다 1 import pygrib 2 import numpy as np 3 import matplotlib.pyplot as plt 4 from netCDF4 import Dataset 5 from wrf import to_np, getvar, smooth2d, get_basemap, latlon_coords 6 7 from matplotlib.cm import get_cmap 8 9 from mpl_toolkits.basemap import Basemap 10 11 f = Dataset('./wrfout_d01_2016-10-06_03.nc') 12 13 t = getvar(f, "..

이번에는 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 제가 사용한 ..