clear
clf
load x.out
load y.out
load t.out
load v.out
load u.out
load w.out
load s.out
load s0.out
load m.out
load whov.out
a=size(whov);
n=round(a(1)/2);
w2hov=whov(n:a(1),:);
t2=t(n:a(1));
load vhov.out
v2hov=vhov(n:a(1),:);
load uhov.out
u2hov=uhov(n:a(1),:);
i=1;
while i == 1
it=input('Plot choices: \n \n Final Time: \n   1) s, u & v \n   2) Mass flux \n   3) w \n   4) u   \n   5) s0 \n \n Hovmueller:  \n   6) Zonal velocity \n   7) Meridional velocity \n   8) Vertical velocity \n \n  0) Exit \n ');
if it == 1
pcolor(x,y,s)
shading interp
hold
quiver(x,y,u,v)
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
b=max(max(s));
a=min(min(s));
title(['Theta_e, from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 2
clf
pcolor(x,y,m)
shading interp
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
b=max(max(m));
a=min(min(m));
title(['Mass flux (cm/s), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 3
clf   
pcolor(x,y,w)
shading interp
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
b=max(max(w));
a=min(min(w));
title(['w (cm/s), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 4
clf
pcolor(x,y,u)
shading interp
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
b=max(max(u));
a=min(min(u));
title(['u (m/s), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 5
clf
pcolor(x,y,s0)
shading interp
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
b=max(max(s0));
a=min(min(s0));
title(['Theta_{e0} (K), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 6
clf
contourf(x,t2,u2hov)
xlabel('x (km)','fontweight','bold')
ylabel('Time (days)','fontweight','bold')
b=max(max(u2hov));
a=min(min(u2hov));
title(['u (m/s), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 7
clf
contourf(x,t2,v2hov)
xlabel('x (km)','fontweight','bold')
ylabel('Time (days)','fontweight','bold')
b=max(max(v2hov));
a=min(min(v2hov));
title(['v (m/s), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 8
clf
contourf(x,t2,w2hov)
xlabel('x (km)','fontweight','bold')
ylabel('Time (days)','fontweight','bold')
b=max(max(w2hov));
a=min(min(w2hov));
title(['w (cm/s), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 0
break
end
end
