;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"    
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$HAO_NCL_LIB/halfzp2full.ncl"
;************************************************
begin
;************************************************
; before the data can be plotted, there is a 
; sequence of things that have to be done to it:
;************************************************
  in       = addfile("/chewy/d/murphys/files/pyrj001.nc","r")
  x        = in->UN                         ; read in data
  x        = halfzp("UN",x)                 ; convert to full zp

  t2plt    = 1                              ; time index to plot
  lon2plt  = -180                           ; lon to plot
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps" ,"zplat")              ; open a ps file

  res                          = True         ; plot mods desired
  res@cnLineLabelPlacementMode = "constant"   ; best line label mode
  res@tiMainString             = "UT = "+in->ut(t2plt)+" Lon = "+lon2plt
  res@gsnDraw                  = False        ; don't draw yet
  res@gsnFrame                 = False        ; don't advance frame yet

  plot = gsn_csm_contour(wks,x(t2plt,:,:,{lon2plt}), res)
  plot = ZeroNegDashLineContour(plot)
  draw(plot)
  frame(wks)
end