timour_matlab/mscripts/plotmeanqgpv.m

clear path 
whos
global Nx Ny Nz
global lat long dz dm mdep
global delt_su su_its t_su delt
global descriptor this_path
global f deltaf Q beta r_expt r_heat H
global time rots it
global g Cp rho_bar alpha
global u v t w 
global iterations

cmin=0.75;
cmax=1.75,
V=[cmin cmax];
figure
caxis('manual')

subplot(1,2,1)
        title='in-situ time mean qgpv at z=2';
        imagesc(lat,long,squeeze(meanq(:,:,2))');shading flat;axis image;caxis(V);colorbar('horizontal');
        text(0,-30,descriptor);
        text(0,140,title);
        text(0,-20,'timestep');text(40,-20,int2str(eval(itstart)));
        text(40,-12,'to');text(40,-4,int2str(eval(itend)));
subplot(1,2,2)
        contour(flipud(squeeze(meanq(:,Ny/2,:))'),100);caxis(V);colorbar;grid
        text(0,-2,'in-situ cross-tank time mean qgpv')


1	edhill	1.1	clear path
2			whos
3			global Nx Ny Nz
4			global lat long dz dm mdep
5			global delt_su su_its t_su delt
6			global descriptor this_path
7			global f deltaf Q beta r_expt r_heat H
8			global time rots it
9			global g Cp rho_bar alpha
10			global u v t w
11			global iterations
12
13			cmin=0.75;
14			cmax=1.75,
15			V=[cmin cmax];
16			figure
17			caxis('manual')
18
19			subplot(1,2,1)
20			title='in-situ time mean qgpv at z=2';
21			imagesc(lat,long,squeeze(meanq(:,:,2))');shading flat;axis image;caxis(V);colorbar('horizontal');
22			text(0,-30,descriptor);
23			text(0,140,title);
24			text(0,-20,'timestep');text(40,-20,int2str(eval(itstart)));
25			text(40,-12,'to');text(40,-4,int2str(eval(itend)));
26			subplot(1,2,2)
27			contour(flipud(squeeze(meanq(:,Ny/2,:))'),100);caxis(V);colorbar;grid
28			text(0,-2,'in-situ cross-tank time mean qgpv')
29
30