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
load vhov.out
load uhov.out
load hhov.out
load shov.out
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   9) Surface theta_e    \n \n  0) Exit \n ');
if it == 1
pcolor(x,y,s)
shading interp
colorbar
hold
quiver(x,y,u,v)
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
set(gca,'fontweight','bold')
b=max(max(s));
a=min(min(s));
title(['\theta_e, from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
hold off
elseif it == 2
clf
pcolor(x,y,m)
shading interp
colorbar
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
set(gca,'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
colorbar
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
set(gca,'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
colorbar
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
set(gca,'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
colorbar
xlabel('x (km)','fontweight','bold')
ylabel('y (km)','fontweight','bold')
set(gca,'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,t,uhov)
colorbar
xlabel('x (km)','fontweight','bold')
ylabel('Time (days)','fontweight','bold')
set(gca,'fontweight','bold')
b=max(max(uhov));
a=min(min(uhov));
title(['u (m/s), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it ==  7
clf
contourf(x,t,vhov)
colorbar
xlabel('x (km)','fontweight','bold')
ylabel('Time (days)','fontweight','bold')
set(gca,'fontweight','bold')
b=max(max(vhov));
a=min(min(vhov));
title(['v (m/s), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 8
clf
contourf(x,t,whov)
colorbar
xlabel('x (km)','fontweight','bold')
ylabel('Time (days)','fontweight','bold')
set(gca,'fontweight','bold')
b=max(max(whov));
a=min(min(whov));
title(['w (cm/s), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 9
clf
contourf(x,t,shov)
colorbar
xlabel('x (km)','fontweight','bold')
ylabel('Time (days)','fontweight','bold')
set(gca,'fontweight','bold')
b=max(max(shov));
a=min(min(shov));
title(['Surface \theta_e (K), from ',num2str(a),' to ',num2str(b)],'fontweight','bold')
elseif it == 0
break
end
end
