llc_1080/matlab/lookat.m

ts=354300;
pn='/nobackupp8/dmenemen/llc/llc_1080/MITgcm/run2/';
fld1={'Eta',
'KPPhbl',
'PhiBot',
'SIarea',
'SIheff',
'SIhsalt',
'SIhsnow',
'SIuice',
'SIvice',
'Salt',
'Theta',
'U',
'V',
'W',
'oceFWflx',
'oceQnet',
'oceQsw',
'oceSflux',
'oceTAUX',
'oceTAUY'};
fld2={'_ETAN',
'_KPPhbl',
'_PHIBOT',
'_SIarea',
'_SIheff',
'_SIhsalt',
'_SIhsnow',
'_SIuice',
'_SIvice',
'_SALTanom',
'_THETA',
'_UVELMASS',
'_VVELMASS',
'_WVELMASS',
'_oceFWflx',
'_oceQnet',
'_oceQsw',
'_oceSflux',
'_oceTAUX',
'_oceTAUY'};

for f=1:20
    f1=[pn fld1{f} '.' myint2str(ts-20,10) '.data'];
    f2=[pn fld2{f} '.' myint2str(ts,10) '.data'];
    fl1=quikread_llc(f1,1080);
    fl2=quikread_llc(f2,1080);
    if f==10, fl2=fl2+35; end
    clf
    subplot(311)
    quikplot_llc(fl1)
    thincolorbar
    title(fld1{f})
    subplot(312)
    quikplot_llc(fl2)
    thincolorbar
    title(fld2{f})
    subplot(313)
    quikplot_llc(fl2-fl1)
    thincolorbar
    pause
end

k=88;
for f=10:13
    f1=[pn fld1{f} '.' myint2str(ts-20,10) '.data'];
    f2=[pn fld2{f} '.' myint2str(ts,10) '.data'];
    fl1=quikread_llc(f1,1080,k);
    fl2=quikread_llc(f2,1080,k);
    if f==10, fl2=fl2+35; end
    clf
    subplot(311)
    quikplot_llc(fl1)
    thincolorbar
    title(fld1{f})
    subplot(312)
    quikplot_llc(fl2)
    thincolorbar
    title(fld2{f})
    subplot(313)
    quikplot_llc(fl2-fl1)
    thincolorbar
    pause
end


=============

ts1=430080;
ts2=431040;
u1=quikread_llc(['U.' myint2str(ts1,10) '.data'],1080);
v1=quikread_llc(['V.' myint2str(ts1,10) '.data'],1080);
e1=quikread_llc(['Eta.' myint2str(ts1,10) '.data'],1080);
u2=quikread_llc(['U.' myint2str(ts2,10) '.data'],1080);
v2=quikread_llc(['V.' myint2str(ts2,10) '.data'],1080);
e2=quikread_llc(['Eta.' myint2str(ts2,10) '.data'],1080);

clf
orient tall
subplot(311)
quikplot_llc(e1)
caxis([-1 1]*2.7)
thincolorbar
title(['Eta on ' ts2dte(ts1,90,2010)])
subplot(312)
quikplot_llc(e2)
caxis([-1 1]*2.7)
thincolorbar
title(['Eta on ' ts2dte(ts2,90,2010)])
subplot(313)
quikplot_llc(e2-e1)
caxis([-1 1]*.7)
thincolorbar
title('Difference')
eval(['print -dpsc Eta_day' int2str(ts2*90/60/60/24)])

clf
orient tall
subplot(311)
quikplot_llc(sqrt(u1.^2+v1.^2))
caxis([0 1]*.7)
thincolorbar
title(['Surface speed on ' ts2dte(ts1,90,2010)])
subplot(312)
quikplot_llc(sqrt(u2.^2+v2.^2))
caxis([0 1]*.7)
thincolorbar
title(['Surface speed on ' ts2dte(ts2,90,2010)])
subplot(313)
quikplot_llc(sqrt(u2.^2+v2.^2)-sqrt(u1.^2+v1.^2))
caxis([-1 1]*.25)
thincolorbar
title('Difference')
eval(['print -dpsc Speed_day' int2str(ts2*90/60/60/24)])
1	ts=354300;
2	pn='/nobackupp8/dmenemen/llc/llc_1080/MITgcm/run2/';
3	fld1={'Eta',
4	'KPPhbl',
5	'PhiBot',
6	'SIarea',
7	'SIheff',
8	'SIhsalt',
9	'SIhsnow',
10	'SIuice',
11	'SIvice',
12	'Salt',
13	'Theta',
14	'U',
15	'V',
16	'W',
17	'oceFWflx',
18	'oceQnet',
19	'oceQsw',
20	'oceSflux',
21	'oceTAUX',
22	'oceTAUY'};
23	fld2={'_ETAN',
24	'_KPPhbl',
25	'_PHIBOT',
26	'_SIarea',
27	'_SIheff',
28	'_SIhsalt',
29	'_SIhsnow',
30	'_SIuice',
31	'_SIvice',
32	'_SALTanom',
33	'_THETA',
34	'_UVELMASS',
35	'_VVELMASS',
36	'_WVELMASS',
37	'_oceFWflx',
38	'_oceQnet',
39	'_oceQsw',
40	'_oceSflux',
41	'_oceTAUX',
42	'_oceTAUY'};
43
44	for f=1:20
45	f1=[pn fld1{f} '.' myint2str(ts-20,10) '.data'];
46	f2=[pn fld2{f} '.' myint2str(ts,10) '.data'];
47	fl1=quikread_llc(f1,1080);
48	fl2=quikread_llc(f2,1080);
49	if f==10, fl2=fl2+35; end
50	clf
51	subplot(311)
52	quikplot_llc(fl1)
53	thincolorbar
54	title(fld1{f})
55	subplot(312)
56	quikplot_llc(fl2)
57	thincolorbar
58	title(fld2{f})
59	subplot(313)
60	quikplot_llc(fl2-fl1)
61	thincolorbar
62	pause
63	end
64
65	k=88;
66	for f=10:13
67	f1=[pn fld1{f} '.' myint2str(ts-20,10) '.data'];
68	f2=[pn fld2{f} '.' myint2str(ts,10) '.data'];
69	fl1=quikread_llc(f1,1080,k);
70	fl2=quikread_llc(f2,1080,k);
71	if f==10, fl2=fl2+35; end
72	clf
73	subplot(311)
74	quikplot_llc(fl1)
75	thincolorbar
76	title(fld1{f})
77	subplot(312)
78	quikplot_llc(fl2)
79	thincolorbar
80	title(fld2{f})
81	subplot(313)
82	quikplot_llc(fl2-fl1)
83	thincolorbar
84	pause
85	end
86
87
88	=============
89
90	ts1=430080;
91	ts2=431040;
92	u1=quikread_llc(['U.' myint2str(ts1,10) '.data'],1080);
93	v1=quikread_llc(['V.' myint2str(ts1,10) '.data'],1080);
94	e1=quikread_llc(['Eta.' myint2str(ts1,10) '.data'],1080);
95	u2=quikread_llc(['U.' myint2str(ts2,10) '.data'],1080);
96	v2=quikread_llc(['V.' myint2str(ts2,10) '.data'],1080);
97	e2=quikread_llc(['Eta.' myint2str(ts2,10) '.data'],1080);
98
99	clf
100	orient tall
101	subplot(311)
102	quikplot_llc(e1)
103	caxis([-1 1]*2.7)
104	thincolorbar
105	title(['Eta on ' ts2dte(ts1,90,2010)])
106	subplot(312)
107	quikplot_llc(e2)
108	caxis([-1 1]*2.7)
109	thincolorbar
110	title(['Eta on ' ts2dte(ts2,90,2010)])
111	subplot(313)
112	quikplot_llc(e2-e1)
113	caxis([-1 1]*.7)
114	thincolorbar
115	title('Difference')
116	eval(['print -dpsc Eta_day' int2str(ts2*90/60/60/24)])
117
118	clf
119	orient tall
120	subplot(311)
121	quikplot_llc(sqrt(u1.^2+v1.^2))
122	caxis([0 1]*.7)
123	thincolorbar
124	title(['Surface speed on ' ts2dte(ts1,90,2010)])
125	subplot(312)
126	quikplot_llc(sqrt(u2.^2+v2.^2))
127	caxis([0 1]*.7)
128	thincolorbar
129	title(['Surface speed on ' ts2dte(ts2,90,2010)])
130	subplot(313)
131	quikplot_llc(sqrt(u2.^2+v2.^2)-sqrt(u1.^2+v1.^2))
132	caxis([-1 1]*.25)
133	thincolorbar
134	title('Difference')
135	eval(['print -dpsc Speed_day' int2str(ts2*90/60/60/24)])