;================================================;
;  Popscal_1.ncl
;================================================;
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
; ================================================;
begin
  f         = addfile("TEMP.nc","r")

  t         = f->TEMP

  lat2d     = f->TLAT
  lon2d     = f->TLONG
  t@lon2d   = lon2d
  t@lat2d   = lat2d
;=================================================;
; Create plot
;=================================================;
  wks  = gsn_open_wks("ps","popscal")               ; open a ps file

  res                     = True
  res@cnLineLabelDensityF = 4.0                     ; increase dens of labels
  res@gsnAddCyclic        = True 
  plot = gsn_csm_contour_map_ce(wks,t(0,0,:,:),res)        
  end