clear all
file_name='sylt.2d.nc';
bathy=getnc(file_name,'bathymetry');
elev=getnc(file_name,'elev');
u1=getnc(file_name,'u');
v1=getnc(file_name,'v');
tmax=2;
%tmax=20;
scale=3;
m=2;
r=datenum(2007,01,01);


for j=2:tmax
    h=figure(1);
    clf
    elev_t=squeeze(elev(j,:,:));
    elev_tn=elev_t(1:m:160,1:m:135);
    u=squeeze(u1(j,:,:));
    v=squeeze(v1(j,:,:));
    pcolor(elev_tn)
    un=u(1:m:160,1:m:135);
    vn=v(1:m:160,1:m:135);
    shading interp;
    hold on
    quiver(un,vn,scale,'color','k');
    get(gca,'XTick')
    get(gca,'YTick')
    set(gca,'XTickLabel',[4;8;12;16;20;24]);
    set(gca,'YTickLabel',[4;8;12;16;20;24;30;32]);
    hugo=strcat('','elevation in sylt test case [',datestr(r+(j*0.0052)),']','')
    title((hugo),'fontsize',16)
    xlabel('x / km','fontsize',14);
    ylabel('y / km','fontsize',14);
    caxis([-1 1]);
    cb=colorbar;
    set(get(cb,'ylabel'),'string','Surface Elevation / m','fontsize',14)
    tmp=strcat('print -dpng',' /yourfolder/getm_data/',num2str(j+1000),'.png');
    eval(tmp);
end