begin

  a = addfile("wrfout_d01_2017-07-27_12:00:00.nc","r")

  hgt = a->HGT(0,:,:)                 ; terrain, 0 is the first time step

  lat = a->XLAT(0,:,0)                ; latitude
  lon = a->XLONG(0,0,:)               ; longitude 

  xi = lat
  yi = lon
  zi = -hgt/2000

;print(lat)
;print(lon)
;print(hgt)


  wks = gsn_open_wks("png","tdpack")   ; send graphics to PNG file
;
;  Draw plot.
;
  rho   = -3.
  theta = -0.
  phi   = -15.

  tdez2d(wks, xi,yi,zi, rho, theta, phi, -2)
  frame(wks)
end