cmin=20;
cmax=30,
V=[cmin cmax];
figure
caxis('manual')

cmin=20;
cmax=30,
V=[cmin cmax];
figure
caxis('manual')

subplot(1,2,1)
        title='in-situ time mean temperature at z=2';
        imagesc(lat,long,squeeze(meantheta(:,:,2))');shading flat;axis image;caxis(V);colorbar('horizontal');
        text(0,-30,descriptor);
        text(0,140,title);
        text(0,-20,'timestep');text(50,-20,num2str(it))
subplot(1,2,2)
        contour(flipud(squeeze(meantheta(:,Ny/2,:))'),50);caxis(V);colorbar;grid
        text(0,-2,'in-situ time mean central temperature east-west')