
function plot_stations(col,lonLat)

%% LD stations
lat_stat_LD   = -19;
lon_stat_LD(1) = 155;
lon_stat_LD(2) = 176.32;
lon_stat_LD(3) = 200;
label_stat_LD(1) = 'A';
label_stat_LD(2) = 'B';
label_stat_LD(3) = 'C';


%% SD stations
SD_stat = load('SD_stations_coord.txt');
lon_stat_SD = SD_stat(2:18,2);
lat_stat_SD = SD_stat(2:18,1);


%% Check geographical limits
indOK = lon_stat_LD >= lonLat(1) & lon_stat_LD <= lonLat(2);
lon_stat_LD = lon_stat_LD(indOK);

indOK = lon_stat_SD >= lonLat(1) & lon_stat_SD <= lonLat(2);
lon_stat_SD = lon_stat_SD(indOK);
lat_stat_SD = lat_stat_SD(indOK);
label_stat_SD = num2str(find(indOK));


%% Plot stations
deltaStat = (lonLat(4)-lonLat(3))/20;

for iLD = 1 : length(lon_stat_LD)
    plot(lon_stat_LD(iLD),lat_stat_LD,'+','Color',col,'Markersize',10);
    text(lon_sta(iLD),lat_stat_LD+deltaStat,label_stat_LD(iLD),'Color',col,'FontWeight','bold','FontSize',14);
end

plot(lon_stat_SD,lat_stat_SD,'+','Markersize',7,'Color',col);
text(lon_stat_SD,lat_stat_SD+deltaStat,label_stat_SD,'Color',col,'FontWeight','bold','FontSize',12);
% text(lon_stat_SD(9:end)-0.5,lat_stat_SD(9:end)+0.9,num2str([11:19]'),'Color',col);


%% Plot territorial waters
% New Caledonia
if lonLat(1) < 175
    load ../ZEE_NewCaledonia.txt
    plot([ZEE_NewCaledonia(:,2)],[ZEE_NewCaledonia(:,1)],'-','Color',col,'Linewidth',2);
    text(158.2,-22,'NEW','Color',col,'FontWeight','bold','FontSize',14);
    text(158.2,-22.5,'CALEDONIA','Color',col,'FontWeight','bold','FontSize',14);
end

% Vanuatu
if lonLat(1) < 175
    load ../ZEE_Vanuatu.txt
    plot([ZEE_Vanuatu(:,2)],[ZEE_Vanuatu(:,1)],'-','Color',col,'Linewidth',2);
    h = text(166.5,-17,'VANUATU','Color',col,'FontWeight','bold','FontSize',14);
    get(h)
end

% Tonga
if lonLat(1) > 188 && lonLat(2) < 195
    plot([189 189],[-19.7 -18.3],'-','Color',col,'Linewidth',5);
    plot([194 194],[-19.7 -18.3],'-','Color',col,'Linewidth',5);
    text(188.5,-21,'TONGA','Color',col,'FontWeight','bold','FontSize',14);
end

% Cook islands
if lonLat(1) > 198 && lonLat(2) < 210
    plot([199 199],[-19.7 -18.3],'-','Color',col,'Linewidth',5);
    plot([209 209],[-19.7 -18.3],'-','Color',col,'Linewidth',5);
    text(201,-21,'COOK','Color',col,'FontWeight','bold','FontSize',14);
    text(201,-22,'ISLANDS','Color',col,'FontWeight','bold','FontSize',14);
end

% Niue
if lonLat(1) > 193 && lonLat(2) < 200
    plot([194 194],[-19.7 -18.3],'-','Color',col,'Linewidth',5);
    plot([199 199],[-19.7 -18.3],'-','Color',col,'Linewidth',5);
    text(194.5,-17,'NIUE','Color',col,'FontWeight','bold','FontSize',14);
end

% Fidji
if lonLat(1) > 188
    plot([189 189],[-19.7 -18.3],'-','Color',col,'Linewidth',5);
    text(180.5,-21,'FIDJI','Color',col,'FontWeight','bold','FontSize',14);
end

