%
% Load data from files
%
filename=input('Enter the name of .sq4 file: ','s');
load(filename);
vsq4=eval(blank(filename(1:length(filename)-4)));
%
% Set data arrays
%
x=vsq4(:,1);
y=vsq4(:,2);
u=vsq4(:,3);
v=vsq4(:,4);
w=vsq4(:,5);
rho=vsq4(:,6);
T=vsq4(:,7);
S=vsq4(:,8);
q2=vsq4(:,9);
q2l=vsq4(:,10);
enzm=vsq4(:,11);
enzh=vsq4(:,12);
enzq=vsq4(:,13);
%
% Set fd coordinate arrays
%
nnv=21
size(x);
nn=ans(1);
ny=round(nnv);
nx=round(nn/nnv);
xx=flipud(reshape(x,ny,nx));
yy=flipud(reshape(y,ny,nx));
%
% Make plot
%
newplot='y';
while newplot == 'y',
normalfigure
scalarc=input('Enter variable you desire to contour (u,v,w,rho,T,S,q2,q2l,enzm,enzh,enzq): ','s');
scalar=flipud(reshape(eval(scalarc),ny,nx));;
[smin,smax,ibw]=colorband_fd_xy_cb(xx,yy,scalar,1);
axis('tight')
newplot=input('New variable plot? (y/n): ','s');
end