begin

  ymStrt = 197901                              ; climatology start year/month
  ymLast = 201512                              ;      "      last  year/month


  a = addfile("adaptor.mars.internal-1572848815.075499-2543-15-38ef6f8b-52bf-43ff-9911-3827e1b62589.nc","r")


  time   = a->time         
  yyyymm = cd_calendar(time, -1)
  printVarSummary(yyyymm)
  print(yyyymm)
  ntStrt = ind(yyyymm.eq.ymStrt)                ; start time index
  ntLast = ind(yyyymm.eq.ymLast)                ; last  time index
  u        = a->u10            ; (time,lat,lon)
  v        = a->v10            ; (time,lat,lon)
  ntim = dimsizes(time)

   printVarSummary(u)

  months = (/"January", "February", "March", "April" \
            ,"May", "June", "July", "August"         \
            ,"September", "October", "November"      \
            ,"December" /)
;************************************************
; Compute the climatology using a function in contributed.ncl
;************************************************
  uClm = clmMonTLL (u)
  vClm = clmMonTLL (v)                      ; monthly climatology
  printVarSummary(uClm)
  printVarSummary(vClm)

 wspd = wind_speed(uClm,vClm)
;************************************************
; create colors
;*************************************************
  wks = gsn_open_wks("png","climo_ghcn")             ; send graphics to PNG file

;************************************************
; create panel plots
;*************************************************
  plot   = new (12, graphic)                    ; create graphical array

  res                      = True               ; plot options desired
  res@gsnDraw              = False
  res@gsnFrame             = False
  res@cnFillOn             = True               ; turn on color fill
  res@cnFillPalette        = "amwg256"          ; set color map
  res@cnInfoLabelOn        = False              ; turn off contour info label
  res@cnLinesOn            = False              ; turn off contour lines
  res@cnLineLabelsOn       = False              ; turn off line labels

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

  res@mpFillOn             = False              ; turn off gray continents
  res@mpCenterLonF         = 0                  ; Centers the plot at 0 
  res@lbLabelBarOn         = False              ; No single label bar

  resP                     = True               ; panel options
  
  resP@gsnMaximize         = True               ; maximize image
  resP@gsnPanelLabelBar    = True               ; Add common label bar 



  



  do nmo=0,11                                   ; loop over the months
     res@gsnLeftString     = months(nmo)
     plot(nmo) = gsn_csm_contour_map(wks,wspd(nmo,:,:), res)  ; create plot

      end do

  gsn_panel(wks,plot,(/4,3/),resP)
 

end