;*************************************************
; proj_3.ncl
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" 
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
  a = addfile("atmos.nc","r")
;************************************************
; read in data
;************************************************
  t = a->TS(0,:,:)                                
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps","proj")          ; open a ncgm file
  gsn_define_colormap(wks,"gui_default")


  res                   = True
  res@mpProjection      = "Orthographic"       ; choose projection
  res@mpPerimOn         = False             ; turn off box around plot
  res@mpFillOn          = False
  res@mpCenterLonF      = -120.             ; choose center lon
  res@mpCenterLatF      = 50.               ; choose center lat

  res@cnFillOn          = True              ; color plot desired
  res@cnLineLabelsOn    = False             ; turn off contour line labels
  res@cnLinesOn         = False             ; turn off contour lines
;  res@txFontHeightF     = 0.015 

  res@gsnSpreadColors     = True            ; use full range of color map
;  res@gsnSpreadColorStart = 2 
;  res@gsnSpreadColorEnd   = -3 
;  res@lbLabelFontHeightF  = 0.015           ; label bar font height
 
  res@tiMainString       = "Example of Orthographic Projection"   ; add a title
;  res@tiMainFontHeightF  = .018                               ; font height

  plot = gsn_csm_contour_map(wks,t,res)  ; create the plot
end