;****************************************************
; native_2.ncl
;*****************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;*****************************************************
begin
 f1 = addfile("1994_256_FSD.nc","r")
 fsd=f1->FSD

 lat  = fsd&lat
 lon  = fsd&lon
 nlat = dimsizes(lat)
 nlon = dimsizes(lon)
;**************************************
;   Create plot
;***************************************
 wks  = gsn_open_wks("ps", "native")
 gsn_define_colormap(wks,"BlGrYeOrReVi200")         ; choose color map
 i = NhlNewColor(wks,0.7,0.7,0.7)                   ; add gray to colormap

 res                             = True             ; plot mods desired

 res@cnFillOn                    = True             ; turn on color 
 res@cnLinesOn                   = False            ; no contour lines

 res@cnLevelSelectionMode        = "ManualLevels"   ; set manual contour levels
 res@cnMinLevelValF              = 0                ; set min contour level
 res@cnMaxLevelValF              = 70               ; set max contour level
 res@cnLevelSpacingF             = 5                ; contour spacing


 res@gsnSpreadColors             = True             ; use full color map
 res@gsnSpreadColorEnd           = -3               ; don't use added gray
 res@lbOrientation               ="vertical"        ; vertical label bar

 res@mpDataBaseVersion           = "Ncarg4_1"       ; use finer database
 res@pmTickMarkDisplayMode       = "Always"         ; turn on tickmarks

 res@gsnAddCyclic                = False            ; regional data, don't add pt

 res@mpProjection                = "mercator"       ; projection
 res@mpLimitMode                 = "Corners"        ; method to zoom
 res@mpLeftCornerLatF            = min(lat)
 res@mpLeftCornerLonF            = min(lon)
 res@mpRightCornerLatF           = max(lat)       
 res@mpRightCornerLonF           = max(lon)
 res@tfDoNDCOverlay              = True             ; do not transform data

 res@tiMainString                ="Native Mercator Projection"
 plot = gsn_csm_contour_map (wks,fsd(0,:,:),res)    ; create plot

end