202304_10days location

pathi='/data3/yiwen/matlab_coawst/particle_tracking/0926_tracking/'
load([pathi,'spun_up_10days.mat'])  %north-east(-,-)

lon=lon(:,1);lat=lat(1,:);
%======================== case 1 :: harbor ===============================%
x1=linspace(121.9148,121.9163,30);y1=linspace(25.1221,25.1252,30);
[x11,y11]=meshgrid(x1,y1);

s2=linstace(121.9149,121.9158,30);y2=linspace(25.1260,25.1269,30);
[x22,y22]=meshgrid(x2,y2);
%======================== case 2 :: boundary =============================%
% x1=linspace(121.9125,121.9130,30);y1=linspace(25.1242,25.1285,30);
% [x11,y11]=meshgrid(x1,y1);
% % 
% x2=linspace(121.9130,121.9175,30);y2=linspace(25.1280,25.1285,30);
% [x22,y22]=meshgrid(x2,y2);


    %(0hrs:1152 / 6hrs=1212 / 12hrs=1272 / 18hrs=1332) %
    dstart=1152; 
    m=2881;mon1=1; 

for k=1:5:size(xx,2); %timestep
   % sperating area
   plot(xx(1:9:270,k), yy(1:9:270,k), 'o','Color', m_map(k,:),'Markersize',Msize(k));hold on
   plot(xx(270:12:end,k), yy(270:12:end,k), 'o','Color', m_map(k,:),'Markersize',Msize(k));hold on 
   
   plot(xx2(1:8*2:300,k), yy2(1:8*2:300,k), 'o','Color', m_map(k,:),'Markersize',Msize(k));hold on
   plot(xx2(300:8*3:end,k), yy2(300:8*3:end,k), 'o','Color', m_map(k,:),'Markersize',Msize(k));hold on
end