|
Revision 495
(checked in by cbc, 12 years ago)
|
Initial import of Stark code.
|
| Line | |
|---|
| 1 |
% |
|---|
| 2 |
% Load data from files |
|---|
| 3 |
% |
|---|
| 4 |
ls *.f5d |
|---|
| 5 |
filename=input('Enter the name of .f5d file: ','s'); |
|---|
| 6 |
load(filename); |
|---|
| 7 |
vf5d=eval(blank(filename(1:length(filename)-4))); |
|---|
| 8 |
% |
|---|
| 9 |
% Set data arrays |
|---|
| 10 |
% |
|---|
| 11 |
x=vf5d(:,1)/1000.0; |
|---|
| 12 |
y=vf5d(:,2); |
|---|
| 13 |
u=vf5d(:,3); |
|---|
| 14 |
v=vf5d(:,4); |
|---|
| 15 |
w=vf5d(:,5); |
|---|
| 16 |
sigt=vf5d(:,6); |
|---|
| 17 |
T=vf5d(:,7); |
|---|
| 18 |
S=vf5d(:,8); |
|---|
| 19 |
% |
|---|
| 20 |
% Create mesh and create boundary |
|---|
| 21 |
% |
|---|
| 22 |
% nnv=input('Enter nnv: '); |
|---|
| 23 |
nnv=21 |
|---|
| 24 |
size(x); |
|---|
| 25 |
nn=ans(1); |
|---|
| 26 |
nseg=nn/nnv-1; |
|---|
| 27 |
femgen |
|---|
| 28 |
bnd=detbndy(in); |
|---|
| 29 |
xmin=min(x(:,1)); |
|---|
| 30 |
xmax=max(x(:,1)); |
|---|
| 31 |
ymin=min(y(:,1)); |
|---|
| 32 |
ymax=max(y(:,1)); |
|---|
| 33 |
% |
|---|
| 34 |
% Plot boundary |
|---|
| 35 |
% |
|---|
| 36 |
newplot='y'; |
|---|
| 37 |
while newplot == 'y', |
|---|
| 38 |
figure |
|---|
| 39 |
whitebg('w') |
|---|
| 40 |
hold on |
|---|
| 41 |
bndo=plotbnd(x,y,bnd); |
|---|
| 42 |
set(bndo,'Color','k') |
|---|
| 43 |
axis([xmin xmax ymin ymax]) |
|---|
| 44 |
% axis('off') |
|---|
| 45 |
scalarc=input('Enter variable you desire to contour (u,v,w,sigt,T,S): ','s'); |
|---|
| 46 |
scalar=eval(scalarc); |
|---|
| 47 |
scrange(scalar) |
|---|
| 48 |
[smin,smax,ibw]=colorband_fe(in,x,y,bnd,scalar) |
|---|
| 49 |
newplot=input('New variable plot? (y/n): ','s'); |
|---|
| 50 |
end |
|---|
