;************************************************
;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"
;load "$HAO_NCL_LIB/hao_contrib.ncl"
;load "$HAO_NCL_LIB/axis.ncl"
;************************************************
begin
;************************************************
; before the data can be plotted, there is a 
; sequence of things that have to be done to it:
;************************************************
  a  = addfile("/home/ew/NCL/bin/Task/temp.nc","r")
  x        = a->air                         ; read in data
  x        = halfzp("air",x)                 ; convert to full zp

  t2plt    = 0                              ; time index to plot
  lon2plt  = 96                           ; lon to plot
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps" ,"zplat")           ; open a ps file
  gsn_define_colormap(wks,"BlAqGrYeOrRevi200"); choose colormap


  res                          = True         ; plot mods desired
  res@cnLineLabelPlacementMode = "constant"   ; best line label mode
  ;res@tiMainString             = "UT = "+in->ut(t2plt)+" Lon = "+lon2plt
  res@cnFillOn                 = True         ; add color
  res@gsnSpreadColors          = True         ; use full range of color map
  res@lbLabelAutoStride        = True         ; nice lb labels
  res@cnLineLabelsOn           = True         ; turn on line labels

; interpolate to height
  arg      = True
  arg@bott = 1000          ; bottom height to start at
  arg@top  = 1600          ; top height to end at
  arg@dz   = 1000           ; height interval
  z        = a->level        ; get Z
  x_int    = zplev_to_height(z,x,arg)

  plot = gsn_csm_contour(wks,x_int(t2plt,:,:,{lon2plt}), res)

end