;; -*- mode: emacs-lisp; coding: emacs-mule; -*-;  -*- mode: ncl;-*-
;***************************************************************************
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/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"


ncol = 5
ana  = readAsciiTable("Book4.csv", ncol, "float", 1)
lon1=ana(:,0)
lat1=ana(:,1)
no=ana(:,2)
no2=ana(:,3)
bc=ana(:,4)

total= new((/5,21488/),typeof(bc),bc@_FillValue)
total(0,:)=lat1
total(1,:)=lon1
total(2,:)=no
total(3,:)=no2
total(4,:)=bc

lat=lat1
lon=lon1
lat@longname = "latitude"
lat@units = "degrees_north"
lon@longname ="longitude"
lon@units = "degrees_east"

total!0       = "lat"
total!1       = "lon"

;total&lon     =  lon
;total&lat     =  lat
  


wks = gsn_open_wks("eps","NO")
;gsn_define_colormap(wks,"BlWhRe")
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  res                     = True               ; plot mods desired
  res@cnFillOn            = True              ; turn on color   
  res@gsnSpreadColors     = True               ; use full color map
  res@cnLinesOn           = False              ; no contour lines
  res@cnLineLabelsOn      = False              ; no line labels

  ;res@mpMaxLatF =  50.                         ; specify the plot domain
  ;res@mpMinLatF = 39.                         ;                         
  ;res@mpMinLonF = 267                          ;
  ;res@mpMaxLonF = 285    
  ;res@mpCenterLonF= 276 
  ;res@gsnAddCyclic        = False 


;res@lbTitleOn        =  True                ; turn on title
;res@lbTitleString    = "Pa/s"                ; title string
;res@lbTitlePosition  = "Right"              ; title position
;res@lbTitleFontHeightF= .02                ; make title smaller
;res@lbTitleDirection = "Across"             ; title direction


 res@cnLevelSelectionMode =  "AutomaticLevels"   ; AutomaticLevels

; res@cnLevelSelectionMode =  "ManualLevels"   ; AutomaticLevels
;  res@cnMinLevelValF       = 0.              ; min level
;  res@cnMaxLevelValF       =  1.              ; max level
;  res@cnLevelSpacingF      =  0.05               ; interval

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
res@tiMainString        = " NO"

res@gsnDraw              = False           ; Do not draw plot
res@gsnFrame             = False           ; Do not advance frome


plot= gsn_csm_contour_map(wks,total(2,:),res)

  draw (plot)
  frame(wks)
delete(wks)