;*************************************************
; a default black and white plot with the default
; ncl labeling. All neg contour lines are dashed
; and the zero line is made to be double thickness.
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.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->TN                         ; read in data
  x     = halfzp("TN",x)                 ; convert to full zp

  t2plt = 0                              ; time index to plot
  l2plt = 5                              ; level value to plot    
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps" ,"polar")              ; open a ps file

  res                          = True         ; plot mods desired
  res@gsnPolar                 = "SH"         ; default is NH
  res@mpMinLatF                = -30          ; mpMaxLatF for NH
  res@cnLineLabelPlacementMode = "constant"   ; best line label mode
  title = "Perim lat = -30"+" zplev= "+in->lev({l2plt})+" UT = "+in->ut(t2plt)

  res@tiMainString             = title

  plot = gsn_csm_contour_map_polar(wks,x(t2plt,{l2plt},:,:), res)
end