load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
;
begin
;
  wks = gsn_open_wks("x11","weather_sym")
  cmap = (/(/1.,1.,1./),(/0.,0.,0./),(/1.,0.,0./) /)  ; color map
  gsn_define_colormap(wks,cmap)
;
;  Plot some station model data for some cities in Rwanda 
;
; location of map in frame. Size is the same as before.
;
  mpres              = True
;
  mpres@vpXF         =   0.05
  mpres@vpYF         =   1.0
  mpres@vpWidthF     =   0.90
  mpres@vpHeightF    =   0.80
  mpres@mpLimitMode  = "LatLon"
  mpres@mpMinLonF    =   28.85
  mpres@mpMaxLonF    =   30.9
  mpres@mpMinLatF    =   -3.20
  mpres@mpMaxLatF    =   -1.0
  ;mpres@mpCenterLatF =   -2.0
  ;mpres@mpCenterLonF =   30.02
  mpres@mpPerimOn    = True
  mpres@mpOutlineBoundarySets   = "AllBoundaries"
  mpres@mpUSStateLineThicknessF = 1.25

  mpres@tiMainString = "Rwandan stations"
  mpres@tiMainFont   = "Helvetica-Bold"
  mpres@mpFillOn = False
  map = gsn_map(wks,"Satellite",mpres)
;
;  Flag the station model procedure that the wind barbs are
;  being drawn over a map.
;
  wmsetp("ezf",1)
;
;  Draw in the foreground color.
;
  wmsetp("col", 1)
;
;  Plot station model data at selected cities.
;
  lons = (/  30.13, 28.92, 29.25, 29.57, 30.5, 30.05, 29.55/)
  lats = (/  -1.97, -2.47, -1.67, -2.48, -2.15, -1.6, -1.58/)
  imdat= (/ "11721700181008020050300004983052026604007289086925",  \
            "11060032571033020380300004033056030610507808089258",  \
            "11854813511029020330300004025054016609507726087036",  \
            "11000022751126021360300004955054054600007757087712",  \
            "11515500121004020000300004975050034603017207084703",  \
            "11751718481027020310300004021053012609017685086925",  \
            "11206227031102021040300004963056046601517084081470" /)

  wmstnm(wks,lats,lons,imdat)
  wks = gsn_open_wks("x11","weather_sym")

end