clear all
file_name='ns_6nm.2d.2007-05.nc';
bathy=getnc(file_name,'bathymetry');
elev=getnc(file_name,'elev');
lat=getnc(file_name,'latc');
lon=getnc(file_name,'lonc');
u=getnc(file_name,'u');
v=getnc(file_name,'v');
tmax=401;
%tmax=2;
scale=2;
m=2;
r=datenum(1997,05,01);


for j=1:tmax;
    h=figure(1);
    clf
    elev_t=squeeze(elev(j,:,:));
    elev_tn=elev_t(1:87,1:111);
    u1=squeeze(u(j,:,:));
    v1=squeeze(v(j,:,:));
    pcolor(lon,lat,elev_tn)
    un=u1(1:m:87,1:m:111);
    vn=v1(1:m:87,1:m:111);
    lon1=lon(1:m:87,1:m:111);
    lat1=lat(1:m:87,1:m:111);
    shading flat;
    hold on
    quiver(lon1,lat1,un,vn,scale,'color','k');       
    xlabel('lat /°E','fontsize',14);
    ylabel('lon /°N','fontsize',14);
    axis([-4 16 50.5 60])
    caxis([-2.39369 2.01685]);
    hugo=strcat('','surface elevation in North Sea test case [',datestr(r+(j*1/24)),']','')
    title((hugo),'fontsize',16)  
    cb=colorbar('location','southoutside','fontsize',14);    
    set(get(cb,'xlabel'),'string','Elevation/ m','fontsize',14)
    tmp=strcat('print -dpng',' /yourfolder/',num2str(j+1000),'.png');
    eval(tmp);
end