%plot_roms_drift
% 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;
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);

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

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

% 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_a',i);
end

close(nc);
anim_make