Revision 495
(checked in by cbc, 11 years ago)
|
Initial import of Stark code.
|
Line | |
---|
1 |
% |
---|
2 |
% Load data from files |
---|
3 |
% |
---|
4 |
filename=input('Enter the name of .sq4 file: ','s'); |
---|
5 |
load(filename); |
---|
6 |
vsq4=eval(blank(filename(1:length(filename)-4))); |
---|
7 |
% |
---|
8 |
% Set data arrays |
---|
9 |
% |
---|
10 |
x=vsq4(:,1); |
---|
11 |
y=vsq4(:,2); |
---|
12 |
u=vsq4(:,3); |
---|
13 |
v=vsq4(:,4); |
---|
14 |
w=vsq4(:,5); |
---|
15 |
rho=vsq4(:,6); |
---|
16 |
T=vsq4(:,7); |
---|
17 |
S=vsq4(:,8); |
---|
18 |
q2=vsq4(:,9); |
---|
19 |
q2l=vsq4(:,10); |
---|
20 |
enzm=vsq4(:,11); |
---|
21 |
enzh=vsq4(:,12); |
---|
22 |
enzq=vsq4(:,13); |
---|
23 |
% |
---|
24 |
% Set fd coordinate arrays |
---|
25 |
% |
---|
26 |
nnv=21 |
---|
27 |
size(x); |
---|
28 |
nn=ans(1); |
---|
29 |
ny=round(nnv); |
---|
30 |
nx=round(nn/nnv); |
---|
31 |
xx=flipud(reshape(x,ny,nx)); |
---|
32 |
yy=flipud(reshape(y,ny,nx)); |
---|
33 |
% |
---|
34 |
% Make plot |
---|
35 |
% |
---|
36 |
newplot='y'; |
---|
37 |
while newplot == 'y', |
---|
38 |
normalfigure |
---|
39 |
scalarc=input('Enter variable you desire to contour (u,v,w,rho,T,S,q2,q2l,enzm,enzh,enzq): ','s'); |
---|
40 |
scalar=flipud(reshape(eval(scalarc),ny,nx));; |
---|
41 |
[smin,smax,ibw]=colorband_fd_xy_cb(xx,yy,scalar,1); |
---|
42 |
axis('tight') |
---|
43 |
newplot=input('New variable plot? (y/n): ','s'); |
---|
44 |
end |
---|