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	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