begin
  ncout = addfile("/Users/isezen/Dropbox/wrf/wrfout_d02_2015-09-06_00_00_00","r")
  ; tgrid = ncout->W(0,:,:,:)
  ; ugrid = ncout->U(0,:,:,:)
  ; vgrid = ncout->V(0,:,:,:)
  u = wrf_user_getvar(ncout,"ua",-1) ; 3D U at mass points
  v = wrf_user_getvar(ncout,"va",-1) ; 3D V at mass points
  w = wrf_user_getvar(ncout,"wa",-1) ; 3D V at mass points
  tgrid = w(0,:,:,:)
  ugrid = u(0,:,:,:)
  vgrid = v(0,:,:,:)


wtype = "png"
;
; These lines were commented out due to the fact that they make the plot
; unmanageably large.
; wtype@wkWidth = 2500
; wtype@wkHeight = 2500
  wks = gsn_open_wks(wtype,"stream")
;
; The following line is for debugging purposes.  Uncomment it to debug.
;  wks = gsn_open_wks("x11","stream")

  mpres                       = True

  mpres@gsnMaximize           = True
  mpres = wrf_map_resources(ncout, mpres)
  mpres@tfDoNDCOverlay   = True

;---Zoom in on map
  ; mpres@mpMinLatF             = -12
  ; mpres@mpMaxLatF             =  68
  ; mpres@mpMinLonF             = -44
  ; mpres@mpMaxLonF             =  35

  mpres@mpOutlineOn           = True
  mpres@mpDataBaseVersion     = "MediumRes"   ; slightly better outlines
  mpres@mpLandFillColor = "gray91"

;---Change the streamline levels
  ;mpres@stLevelSpacingF       = 30.

  mpres@stLevelPalette        = "ncl_default" ; "MPL_Spectral"
  mpres@stLineThicknessF      = 5.0         ; default is 1.0
  ;mpres@stLevelSelectionMode  = "ExplicitLevels"
  ;mpres@stLevels = ispan(245,305,5)

  ;mpres@pmTickMarkDisplayMode = "Always"    ; more detailed tickmarks

  mpres@tiMainOn = False

  mpres@gsnDraw = False
  mpres@gsnFrame = False
  mpres@lbLabelBarOn = False

  plot = new(2,graphic)
  plot(0) = gsn_csm_streamline_scalar_map(wks, ugrid(0,:,:), vgrid(0,:,:), tgrid(0,:,:), mpres)
  plot(1) = gsn_csm_streamline_scalar_map(wks, ugrid(3,:,:), vgrid(3,:,:), tgrid(3,:,:), mpres)

  panres = True
  panres@gsnPanelLabelBar = True
  panres@lbOrientation = "vertical"

  gsn_panel(wks,plot,(/2,1/),panres)
end