;-----------------------------------------------------------------------
pro plot_cont,fname,clevs,clipwin
;-----------------------------------------------------------------------

; called from main as:
; plot_cont,'/?dir?/lmich_cont.dat', $
; [10,20,30,40,50,75,100,150], [41.5,-88.5,43.5,-85.5]


spawn,'grep Cont '+fname+' | wc -l',a
reads,a,nc
tvlct,255,0,255,4  ; magenta
gname=' '
openr,1,fname
for i=0,nc-1 do begin
   readf,1,nopts,gname
   clev=strmid(gname,0,10)
   contour_line=fltarr(2,nopts)
   readf,1,contour_line,format='(2f9.5)'
   index=where(clevs eq clev,count)
   if(count ne 0) then begin
      oplot,contour_line(0,*),contour_line(1,*),linestyle=1,color=4
      minlat=min(contour_line(1,*),index)
      if contour_line(0,index) gt clipwin(1) and contour_line(0,index) lt clipwin(3) $
         and minlat gt clipwin(0) and minlat lt clipwin(2) then $
         xyouts,contour_line(0,index),minlat,strtrim(string(clev),2),charsize=.5,color=4
   end
end

; list contour intervals at bottom

; nlev=size(clevs)
; xyouts,0.3,0.04,'Depth Contours (m):',/norm,charsize=1
; for i=0,nlev(1)-1 do begin
;    xyouts,0.475+i*.05,0.04,strtrim(string(clevs(i))),/norm,charsize=1
; end

close,1
return
end