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	dimitri	1.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)])