wks = gsn_open_wks("png","topografia_cuenca")

;---Set some basic plot options
  res               = True

  res@gsnMaximize   = True   ; maximize plot in frame
  res@gsnDraw       = False
  res@gsnFrame      = False

  res@mpProjection  = "CylindricalEquidistant"    ; The default

;---Zoom in on plot
  res@mpMinLatF     = -13.8
  res@mpMaxLatF     = -10.5
  res@mpMinLonF     = -77
  res@mpMaxLonF     = -73.6

;---Additional resources desired
  res@pmTickMarkDisplayMode = "Always"   ; nicer tickmarks

  res@mpDataBaseVersion     = "MediumRes"       ; better and more map outlines
  res@mpDataSetName         = "Earth..4"
  res@mpOutlineBoundarySets = "AllBoundaries"
  res@mpOutlineOn           = True

  res@lbOrientation         = "Vertical"
  res@tiMainOffsetYF        = -0.03           ; Move the title down

  plot = gsn_csm_map(wks,res)

  lat=(/-12.04, -10.69, -11.4, -11.33, -12.25, -12.12, -11.4, -11.75, -11.88, -11.62, \
        -11.78, -12, -12.16, -12.38, -12.35, -12.58, -12.78, -12.98, -13.06, -12.84/)

  lon=(/-75.32, -76.26, -76.33, -76.4, -75.36, -75.43, -75.69, -75.13, -75.27, -75.01,\
        -75.47, -75.22, -75.23, -74.85, -75.08, -75.24, -75.01, -74.72, -74.14, -74.57/)

 
 colors = (/"orange2","forestgreen","dodgerblue2","forestgreen","orange2","orange2","forestgreen","dodgerblue2",\
           "dodgerblue2","olivedrab2","forestgreen","orange2","orange2","olivedrab2","darkorchid2","olivedrab2",\
           "forestgreen","darkorchid2","olivedrab2","darkorchid2"/)

  mkres               = True
  mkres@gsMarkerIndex = 17     ; Filled circle
  mkres@gsMarkerSizeF = 0.05    ;0.03

  nlatlon = dimsizes(lat)
 dum3 = new(nlatlon,graphic)
 do i = 0,19
  
    mkres@gsMarkerColor = colors(i)
    dum3(i)=  gsn_add_polymarker(wks,plot,lon(i),lat(i),mkres)
  end do
 
 draw(plot)
 frame(wks)