작업10: WRF 모델 결과 NCL로 시각화하기3

이번에는 저번 포스팅에서 이슬점 온도와 풍속 풍향 그림의 결과를 보여드리겠습니다

 

코드는 다른게 없으니 코드와 결과 사진만 올려드리겠습니다

---

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
  mpres@mpPerimOn              = True

  mpres@mpGeophysicalLineColor     = "Black"
  mpres@mpGridLineColor            = "Black"
  mpres@mpLimbLineColor            = "Black"
  mpres@mpNationalLineColor        = "Black"
  mpres@mpPerimLineColor           = "Black"
  mpres@mpUSStateLineColor         = "Black"
  mpres@mpGeophysicalLineThicknessF   = 3.0

  times = wrf_user_getvar(a, "times", -1)
  ntimes = dimsizes(times)
  print(ntimes)

  do it = 0, ntimes-1, 2
    print("Working on time :" + times(it))
    res@TimeLabel = times(it)

;calculate Sea Level Temperature
    slp    = wrf_user_getvar(a, "slp", it)
      wrf_smooth_2d(slp, 3)
    tc     = wrf_user_getvar(a, "T2", it)
      tc = tc - 273.15
    td     = wrf_user_getvar(a, "td2", it)
    u      = wrf_user_getvar(a, "U", it)
    v      = wrf_user_getvar(a, "V", it)
    u10    = wrf_user_getvar(a, "U10", it)
    v10    = wrf_user_getvar(a, "V10", it)
    u10    = u10 * 1.94836
    v10    = v10 * 1.94836
      u10@units = "kts"
      v10@units = "kts"

;Plotting T
    opts = res
    opts@cnFillOn          = True
    opts@ContourParameters = (/-10., 30., 5./)
    opts@gsnSpreadColorEnd = -3
    contour_tc             = wrf_contour(a, wks, tc, opts)
    delete(opts)

;Plotting Td
    opts                   = res
    opts@cnFillOn          = True
    opts@cnLinesOn         = True
    opts@cnLineLabelsOn    = True
    opts@ContourParameters = (/-10., 30., 5./)
    opts@cnLineLabelBackgroundColor = -1
    opts@gsnSpreadColorEnd = -3
    contour_td             = wrf_contour(a, wks, td, opts)
    delete(opts)

;Plotting SLP
    opts = res
    opts@cnLineColor       = "Blue"
    opts@cnHighLabelsOn    = True
    opts@cnLowLabelsOn      = True
    opts@ContourParameters = (/1000., 1060., 2./)
    opts@cnLineLabelBackgroundColor = -1
    opts@gsnContourLineThicknessesScale = 2.0
    contour_psl = wrf_contour(a, wks, slp, opts)
    delete(opts)    

;Plotting Wind Vector
    opts = res
    opts@FieldTitle = "Wind"
    opts@NumVectors = 47
    vector = wrf_vector(a, wks, u10, v10, opts)
    delete(opts)
    

    plot = wrf_map_overlays(a, wks, (/contour_tc, contour_psl, vector/), pltres, mpres)
    plot = wrf_map_overlays(a, wks, (/contour_td, vector/), pltres, mpres)
  end do
end

 

결과 입니다