1 |
nx=4320; |
nx=4320; |
2 |
pn='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run/'; |
pn='~dmenemen/llc_4320/MITgcm/run/'; |
3 |
|
|
4 |
|
%%%%%%%%%%%% |
5 |
|
nx=4320; |
6 |
|
pn='~dmenemen/llc_4320/MITgcm/run/'; |
7 |
|
nm='U'; |
8 |
|
ts=468432; |
9 |
|
fn=[pn nm '.' myint2str(ts,10) '.data']; |
10 |
|
fld=read_llc_fkij(fn,nx,3); |
11 |
|
clf |
12 |
|
quikpcolor(rot90(fld)'); |
13 |
|
caxis([-1 1]/4) |
14 |
|
thincolorbar |
15 |
|
|
16 |
%%%%%%%%%%%%%% |
%%%%%%%%%%%%%% |
17 |
% SST of San Pedro basin |
% SST of San Pedro basin |
80 |
clf reset |
clf reset |
81 |
colormap(cmap) |
colormap(cmap) |
82 |
for ts=10368:144:138384; |
for ts=10368:144:138384; |
83 |
mydisp(ts) |
mydisp(ts) |
84 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
85 |
if exist(fn) |
if exist(fn) |
86 |
fld=[rot90(quikread_llc(fn,nx,1,'real*4',7),1) quikread_llc(fn,nx,1,'real*4',8)]; |
% fld=[rot90(quikread_llc(fn,nx,1,'real*4',7),1) quikread_llc(fn,nx,1,'real*4',8)]; |
87 |
clf |
fld=[ rot90(read_llc_fkij(fn,nx,3),1) read_llc_fkij(fn,nx,4,1,1:nx,(2*nx+1):(3*nx))]; |
88 |
quikpcolor(rot90(fld(:,1:5000),1)'); |
clf |
89 |
caxis([0.9 1]) |
quikpcolor(rot90(fld(:,1:5000),1)'); |
90 |
title('llc4320 fractional sea ice concentration') |
caxis([0.9 1]) |
91 |
thincolorbar |
title('llc4320 fractional sea ice concentration') |
92 |
text(3400,760,ts2dte(ts,25,2011,9,10),'color','w') |
thincolorbar |
93 |
eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]); |
text(3400,760,ts2dte(ts,25,2011,9,10),'color','w') |
94 |
n=n+1; |
eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]); |
95 |
end |
n=n+1; |
96 |
|
end |
97 |
end |
end |
98 |
|
|
99 |
%%%%%%%%%%%% |
%%%%%%%%%%%% |