| 1 | cmin=20; | 
| 2 | cmax=30, | 
| 3 | V=[cmin cmax]; | 
| 4 | figure | 
| 5 | caxis('manual') | 
| 6 |  | 
| 7 | cmin=20; | 
| 8 | cmax=30, | 
| 9 | V=[cmin cmax]; | 
| 10 | figure | 
| 11 | caxis('manual') | 
| 12 |  | 
| 13 | subplot(1,2,1) | 
| 14 | title='in-situ time mean temperature at z=2'; | 
| 15 | imagesc(lat,long,squeeze(meantheta(:,:,2))');shading flat;axis image;caxis(V);colorbar('horizontal'); | 
| 16 | text(0,-30,descriptor); | 
| 17 | text(0,140,title); | 
| 18 | text(0,-20,'timestep');text(50,-20,num2str(it)) | 
| 19 | subplot(1,2,2) | 
| 20 | contour(flipud(squeeze(meantheta(:,Ny/2,:))'),50);caxis(V);colorbar;grid | 
| 21 | text(0,-2,'in-situ time mean central temperature east-west') |