plot_sq4_11.m

Revision 495 (checked in by cbc, 11 years ago)	Initial import of Stark code.

Line
1	%
2	% Load data from files
3	%
4	!ls .sq4 .SQ4
5	filename=input('Enter the name of .sq4 file: ','s');
6	load(filename);
7	vsq4=eval(blank(filename(1:length(filename)-4)));
8	%
9	% Set data arrays
10	%
11	variable=['U ';'V ';'W ';'RHO ';'T ';'S ';'Q2 ';'Q2L ';'ENZM';'ENZH';'ENZQ']
12	ip=[1 2 3 4 5 6 7 8 10 11 12]
13	x=vsq4(:,1);
14	y=vsq4(:,2);
15	u=vsq4(:,3);
16	v=vsq4(:,4);
17	w=vsq4(:,5);
18	rho=vsq4(:,6);
19	T=vsq4(:,7);
20	S=vsq4(:,8);
21	q2=vsq4(:,9);
22	q2l=vsq4(:,10);
23	enzm=vsq4(:,11);
24	enzh=vsq4(:,12);
25	enzq=vsq4(:,13);
26	%
27	% Create mesh and create boundary
28	%
29	% nnv=input('Enter nnv: ');
30	nnv=21
31	size(x);
32	nn=ans(1);
33	nseg=nn/nnv-1;
34	femgen
35	bnd=detbndy(in);
36	xmin=min(x(:,1));
37	xmax=max(x(:,1));
38	ymin=min(y(:,1));
39	ymax=max(y(:,1));
40	%
41	% Create Plot
42	%
43	tallfigure
44	for i=1:11
45	subplot(4,3,ip(i))
46	bndo=plotbnd(x,y,bnd);
47	set(bndo,'Color','k')
48	axis([xmin xmax ymin ymax])
49	axis('off')
50	title(variable(i,:))
51	hold on
52	colormeshm(in,x,y,vsq4(:,i+2))
53	colorbar
54	end
55	%
56	%
57	%
58	break
59	!ls iter_0_.sq4
60	filename=input('Enter the name of reference .sq4 file: ','s');
61	load(filename);
62	vsq4r=eval(blank(filename(1:length(filename)-4)));
63	vsq4(:,5)=vsq4r(:,7);
64	vsq4(:,6)=vsq4r(:,8);
65	%
66	% Create Plot
67	%
68	massvar=['initial T';'initial S';'detided T';'detided S'];
69	tallfigure
70	for i=1:4
71	subplot(2,2,i)
72	bndo=plotbnd(x,y,bnd);
73	set(bndo,'Color','k')
74	axis([xmin xmax ymin ymax])
75	axis('off')
76	title(massvar(i,:))
77	hold on
78	colormeshm(in,x,y,vsq4(:,i+4))
79	colorbar
80	end

Note: See TracBrowser for help on using the browser.