%plot_roms_drift2
% plotting ROMS currents, simulated drifters and tracked ones

% loading currents file
ncfile='/sux1/models/toms/new/my25/my_adria03_avg.nc';
nc=netcdf(ncfile);
jd=(nc{'ocean_time'}(:)/86400)+2440000;

% loading coastline
load /d4/rps/p/adriatic/data/coast/coast
lim=[12.1798   15.7260   42.5   45.8602];
set(gcf,'pos',[200 50 600 700]);
set(gcf,'color','white');

% currents at 3 m (?)
load /d4/rps/p/adriatic/data/moorings/w3m.mat loni lati w1hi jdi
[wlp,jdlp]=plfilt(w1hi,jdi);
wlp=interp1(jdlp,wlp,jd);

% loading argos drifters
d=load('/d4/rps/p/adriatic/data/drifters/argos_10min.mat');

% loading floaters file
nc2file='/sux1/models/toms/new/my25/my_adria03_flt.nc';
timesub=15; %subsampling time of acquisition
driftsub=5; %subsampling # od drifters
[lonm,latm,depthm,jdm]=rfloat(nc2file,timesub,driftsub);


% plotting last 30 days
for i=(length(jd)-20):length(jd);
% plotting all values
%for i=1:length(jd);
 clf
% draw velocities 
 [w,lon,lat]=zsliceuv(ncfile,i,-3);
 pslice(lon,lat,abs(w),[0 0.5],'ROMS current (m/s)');dasp(43);
 psliceuv(lon,lat,w,3,.5,'black');
% draw coastline 
 line(coast(:,1),coast(:,2),'color','white');

% draw modeled drifters
 idriftm=find(jdm>(jd(i)-2) & jdm<(jd(i)));
 line(lonm(idriftm,:),latm(idriftm,:),'linewidth',2,'color','magenta');

% draw drifter tracks over past two days
 idrift=find(d.jdi>(jd(i)-2) & d.jdi<(jd(i)));
 line(d.loni(idrift,:),d.lati(idrift,:),'linewidth',2,'color','black');

 title(greg2str(gregorian(jd(i))));
 axis(lim);
 anim_frame('drift_b',i);
end

close(nc);
anim_make