Reference
- 정지원(2007), 『인천 지역 해륙풍의 특성에 한 연구』, 경인교육대학교 석사 학위 논문
- 박재홍(2001), 『강릉 연안 지역 해륙풍의 기후학적 특성 연구』, 부산대학교 석사 학위 논문
- 정민희, 이희관(2006), "순수해륙풍의 순환구조에 대한 모델링 연구", 『대한환경공학회 학술발표 논문집』
- 지효은(2007), 『연안도시지역의 대기환경 예측을 위한 해양효과와 식생효과의 반영』, 부산대학교 대학원 석사학위 논문
- 문승의, 장광미, 이화운(1990), "해륙풍의 수치모델에 관한 연구", 『Asia-Pacific Journal of Atmospheric Sciences』
- 기상자료개방포털, "https://data.kma.go.kr/cmmn/main.do"
Code
Solar radiation
; Example script to produce standard plots for a WRF grav2x run
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "./WRFUserARW.ncl"
begin
;
; The WRF ARW input file.
; This needs to have a ".nc" appended, so just do it.
a = addfile("/home/nwp-1/1773/wrf/WRF-4.0.3/test/em_seabreeze2d_x/wrfout_d01_2018-06-01_00:00:00.nc","r")
type = "png"
wks = gsn_open_wks(type,"plt_Sea2x")
res = True
res@MainTitle = "WRF SEABREEZE"
res@InitTime = False
res@Footer = False
pltres = True
times = wrf_user_list_times(a) ; get times in the file
ntimes = dimsizes(times) ; number of times in the file
FirstTime=True
do it = 0,ntimes-1,1
time = it
res@TimeLabel = times(it)
SWDOWN=wrf_user_getvar(a,"SWDOWN",time)
gsn_define_colormap(wks,"WhiteYellowOrangeRed")
opts_SWDOWN = res
opts_SWDOWN@cnFillOn = True
opts_SWDOWN@ContourParameters = (/ 0., 1100., 10. /)
contour_SWDOWN = wrf_contour(a,wks,SWDOWN,opts_SWDOWN)
plot = wrf_overlays(a,wks,(/contour_SWDOWN/),pltres)
end do
end
Cloud
; Example script to produce standard plots for a WRF grav2x run
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "./WRFUserARW.ncl"
begin
;
; The WRF ARW input file.
; This needs to have a ".nc" appended, so just do it.
a = addfile("/home/nwp-1/1773/wrf/WRF-4.0.3/test/em_seabreeze2d_x/wrfout_d01_2018-06-01_00:00:00.nc","r")
; We generate plots, but what kind do we prefer?
; type = "x11"
type = "png"
; type = "ps"
; type = "ncgm"
wks = gsn_open_wks(type,"CLOUD")
; Set some basic resources
res = True
res@MainTitle = "WRF SEABREEZE"
res@InitTime = False
res@Footer = False
pltres = True
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; What times and how many time steps are in the data set?
times = wrf_user_list_times(a) ; get times in the file
ntimes = dimsizes(times) ; number of times in the file
; The specific plane we want to plot data on
plane = (/ 20., 0./) ; (x,y) point for vertical plane
angle = 90.0
pii = 3.14159
; This is the big loop over all of the time periods to process.
FirstTime=True
do it = 0,ntimes-1,1
time = it
res@TimeLabel = times(it)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; First get the variables we will need
CLDFRA = wrf_user_getvar(a,"CLDFRA",time)
z = wrf_user_getvar(a, "z",time) ; grid point height
CLDFRA_plane = wrf_user_intrp3d( CLDFRA,z,"v",plane,angle,False)
; get height info for labels ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
if (FirstTime) then
zmin = 0.
zmax = 10. ;(km)
nz = floattoint(zmax + 1)
zz = wrf_user_intrp3d(z,z,"v",plane,angle,False)
b = ind(zz(:,0) .gt. zmax*1000. )
zmax_pos = b(0) - 1
if ( abs(zz(zmax_pos,0)-zmax*1000.) .lt. abs(zz(zmax_pos+1,0)-zmax*1000.) ) then
zspan = b(0) - 1
else
zspan = b(0)
end if
zmax_pos = zspan
zmax = zz(zmax_pos,0)/1000.
delete(b)
FirstTime = False
end if
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
res@tmYLMode = "Explicit"
res@tmYLValues = fspan(0,zspan,nz)
res@tmYLLabels = sprintf("%.1f",fspan(0,zmax,nz))
gsn_define_colormap(wks,"WhiteBlue")
opts_CLDFRA = res
opts_CLDFRA@cnFillOn = True
opts_CLDFRA@ContourParameters = (/ 0., 1., .1 /)
contour_CLDFRA = wrf_contour(a,wks, CLDFRA_plane(0:zmax_pos,:),opts_CLDFRA)
plot = wrf_overlays(a,wks,(/contour_CLDFRA/),pltres)
; ************************************************************
end do ; end of the time loop
end