function hgrid = plot_eddiesgrid(aa)

axes(aa)

% load waypoint list
waypointsdir = '/Users/donglai/RO/EDDIES/scripts/';
c=load([waypointsdir 'coord_mtx_05new.txt']);
cw=fix(c(:,1)); % extract waypoint number only


%hold on

% if plot station positions (=1)
pgrd=1;
% plot station positions
if(pgrd==1)
    %load coord_mtx_05new.txt;
    %c=coord_mtx_05new;
    pp=axis
    % plot the only numbers inside axis
    ff=find(c(:,3)>pp(1) & c(:,3)<pp(2) & c(:,2)>pp(3) & c(:,2)<pp(4))
    for i=[1:length(ff)]         
        text(c(ff(i),3),c(ff(i),2),['   ',num2str(c(ff(i),1))],'fontsize',7);
    end 
    hold on
    hgrid = plot(c(ff,3),c(ff,2),'sk','markersize',7)  
end 
%hold off