;---read in zonal winds
 ;a = addfile("MAM_GEOS_1anee.cam.h0.2000-12-07-00000.nc","r")
 ;u = a->U(0,35,:,:)    
  a = addfile("uwnd.2008.nc","r")
  u = a->uwnd(0,10,:,:)   
  printVarSummary (u)
 ;print(u)

  wks  = gsn_open_wks("png","adriana")  ; send graphics to PNG file

  res                 = True
  res@gsnMaximize     = True       ; maximize in frame
  res@cnLinesOn       = True
  res@cnFillOn        = True
  res@cnFillPalette   = "BlueYellowRed"   ; change the color palette
  res@tiMainString    = "Default"
  plot = gsn_csm_contour_map(wks,u,res)


  res@mpFillDrawOrder = "PostDraw" ; Make sure land fill is drawn
                                   ; on top of filled contours

  res@tiMainString    = "add gray land fill after contouring"

  plot = gsn_csm_contour_map(wks,u,res)