1 |
nx=4320; |
nx=4320; |
2 |
pn='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run/'; |
pn='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run/'; |
3 |
|
|
4 |
|
%%%%%%%%%%%%%% |
5 |
|
% SST of San Pedro basin |
6 |
|
ts=279360; |
7 |
|
fn=[pn 'Theta.' myint2str(ts,10) '.data']; |
8 |
|
f=11; |
9 |
|
ix=3800:4100; |
10 |
|
iy=3100:3300; |
11 |
|
T=quikread_llc(fn,nx,1,'real*4',f); |
12 |
|
clf reset |
13 |
|
mypcolor(rot90(T(ix,iy),2)') |
14 |
|
colormap(cmap) |
15 |
|
caxis([13.7 16.8]) |
16 |
|
thincolorbar |
17 |
|
set(gca,'xtick',[],'ytick',[]) |
18 |
|
print -dpsc LA |
19 |
|
|
20 |
%%%%%%%%%%%% |
%%%%%%%%%%%% |
21 |
% compare SST off Peru on November 1, 2011 |
% compare SST off Peru on November 1, 2011 |
22 |
cx=[15 20]; |
cx=[15 20]; |