;================================================;
;  Popscal_3.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
  gsn_define_colormap(wks,"BlAqGrYeOrRe")         ; choose colormap

  res                      = True                 ; plot mods desired
  res@cnFillOn             = True                 ; turn on color fill
  res@cnLinesOn            = False                ; turn off contour lines
  res@gsnSpreadColors      = True                 ; use full colormap
  res@gsnSpreadColorStart  = 10                   ; start at color 10
  res@gsnSpreadColorEnd    = 96                   ; end at color 96
  res@gsnAddCyclic         = True                 ; add cyclic point

  plot = gsn_csm_contour_map_ce(wks,t(0,0,:,:),res)        
  end