clear all
file_name='box_curvilinear.3d.nc';
lon=getnc(file_name,'xc')/1000;
lat=getnc(file_name,'yc')/1000;
bathy=getnc(file_name,'bathymetry');
temp=getnc(file_name,'temp');
u1=getnc(file_name,'uu');
v1=getnc(file_name,'vv');
%tmax=1;
tmax=121;
scale=3;


for j=1:tmax
    h=figure(1);
    clf
    temp_t=temp(j,2,:,:); %2 for bottom temp
    pcolor(lon,lat,squeeze(temp_t))
    uu=squeeze(u1(j,2,:,:));
    vv=squeeze(v1(j,2,:,:));
    uun=uu(1:2:30,1:2:100);
    vvn=vv(1:2:30,1:2:100);
    lon1=lon(1:2:30,1:2:100);
    lat1=lat(1:2:30,1:2:100);
    shading interp;
    hold on   
    quiver(lon1,lat1,uun,vvn,scale,'color','k');
    title('bottom temperature in box curvilinear test case','fontsize',16)
    xlabel('x /km','fontsize',14);
    ylabel('y /km','fontsize',14);
    axis equal
    axis([-100 100 70 130]);
    caxis([4. 20.]);
    cb=colorbar('location','southoutside');
    set(get(cb,'xlabel'),'string','Temperature / °C','fontsize',14)    
    tmp=strcat('print -dpng',' /yourfolder/getm_data/',num2str(j+1000),'.png');
    eval(tmp);
end