1 |
nx=4320; |
2 |
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 |
18 |
ts=279360; |
19 |
fn=[pn 'Theta.' myint2str(ts,10) '.data']; |
20 |
f=11; |
21 |
ix=3800:4100; |
22 |
iy=3100:3300; |
23 |
T=quikread_llc(fn,nx,1,'real*4',f); |
24 |
clf reset |
25 |
mypcolor(rot90(T(ix,iy),2)') |
26 |
colormap(cmap) |
27 |
caxis([13.7 16.8]) |
28 |
thincolorbar |
29 |
set(gca,'xtick',[],'ytick',[]) |
30 |
print -dpsc LA |
31 |
|
32 |
%%%%%%%%%%%% |
33 |
% compare SST off Peru on November 1, 2011 |
34 |
cx=[15 20]; |
35 |
n1=1080; |
36 |
n2=2160; |
37 |
n3=4320; |
38 |
ix1=380:590; |
39 |
iy1=300:450; |
40 |
ix2=760:1180; |
41 |
iy2=600:900; |
42 |
ix3=1520:2360; |
43 |
iy3=1200:1800; |
44 |
dte='1-Nov-2011'; |
45 |
ts1=dte2ts(dte,90,2010,1,1); |
46 |
ts2=dte2ts(dte,45,2011,1,17); |
47 |
ts3=dte2ts(dte,25,2011,9,10); |
48 |
p1='/nobackupp8/dmenemen/llc/llc_1080/MITgcm/run_2011/'; |
49 |
p2='/nobackupp5/dmenemen/llc_2160/MITgcm/run_day49_on/'; |
50 |
p3='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run/'; |
51 |
f1=[p1 'T.' myint2str(ts1,10) '.data']; |
52 |
f2=[p2 'Theta.' myint2str(ts2,10) '.data']; |
53 |
f3=[p3 'Theta.' myint2str(ts3,10) '.data']; |
54 |
f=12; |
55 |
T1=quikread_llc(f1,n1,1,'real*4',f); |
56 |
T2=quikread_llc(f2,n2,1,'real*4',f); |
57 |
T3=quikread_llc(f3,n3,1,'real*4',f); |
58 |
clf reset |
59 |
subplot(131) |
60 |
mypcolor(rot90(T1(ix1,iy1),2)'); |
61 |
caxis(cx) |
62 |
thincolorbar('horiz') |
63 |
title('1/12 deg') |
64 |
subplot(132) |
65 |
mypcolor(rot90(T2(ix2,iy2),2)'); |
66 |
caxis(cx) |
67 |
thincolorbar('horiz') |
68 |
title('1/24 deg') |
69 |
subplot(133) |
70 |
mypcolor(rot90(T3(ix3,iy3),2)'); |
71 |
caxis(cx) |
72 |
thincolorbar('horiz') |
73 |
title('1/48 deg') |
74 |
|
75 |
print -dpsc -r600 SST_Peru |
76 |
|
77 |
%%%%%%%%%%%% |
78 |
nm='SIarea'; |
79 |
n=1; |
80 |
clf reset |
81 |
colormap(cmap) |
82 |
for ts=10368:144:138384; |
83 |
mydisp(ts) |
84 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
85 |
if exist(fn) |
86 |
% fld=[rot90(quikread_llc(fn,nx,1,'real*4',7),1) quikread_llc(fn,nx,1,'real*4',8)]; |
87 |
fld=[ rot90(read_llc_fkij(fn,nx,3),1) read_llc_fkij(fn,nx,4,1,1:nx,(2*nx+1):(3*nx))]; |
88 |
clf |
89 |
quikpcolor(rot90(fld(:,1:5000),1)'); |
90 |
caxis([0.9 1]) |
91 |
title('llc4320 fractional sea ice concentration') |
92 |
thincolorbar |
93 |
text(3400,760,ts2dte(ts,25,2011,9,10),'color','w') |
94 |
eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]); |
95 |
n=n+1; |
96 |
end |
97 |
end |
98 |
|
99 |
%%%%%%%%%%%% |
100 |
nm='SIarea'; |
101 |
n=1; |
102 |
clf reset |
103 |
colormap(cmap) |
104 |
ts=3600; |
105 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
106 |
fld=[rot90(quikread_llc(fn,nx,1,'real*4',7),1) quikread_llc(fn,nx,1,'real*4',8)]; |
107 |
clf |
108 |
quikpcolor(rot90(fld(:,1:5000),1)') |
109 |
caxis([0.9 1]) |
110 |
title('llc4320, sea ice concentration (m)') |
111 |
thincolorbar |
112 |
|
113 |
%%%%%%%%%%%% |
114 |
p1='/nobackupp8/dmenemen/llc/llc_1080/MITgcm/run_2011/'; |
115 |
n1=1080; |
116 |
t1=599040; |
117 |
f1=[p1 'T.' myint2str(t1,10) '.data']; |
118 |
fld1=zeros(n1*4,n1); |
119 |
fld1(1:n1,:)=quikread_llc(f1,n1,1,'real*4',2); |
120 |
fld1((n1+1):(2*n1),:)=quikread_llc(f1,n1,1,'real*4',5); |
121 |
fld1((2*n1+1):(3*n1),:)=rot90(quikread_llc(f1,n1,1,'real*4',9),2); |
122 |
fld1((3*n1+1):(4*n1),:)=rot90(quikread_llc(f1,n1,1,'real*4',12),2); |
123 |
|
124 |
p2='/nobackupp5/dmenemen/llc_2160/MITgcm/run_day49_73/'; |
125 |
n2=2160; |
126 |
t2=466080; |
127 |
f2=[p2 'Theta.' myint2str(t2,10) '.data']; |
128 |
fld2=zeros(n2*4,n2); |
129 |
fld2(1:n2,:)=quikread_llc(f2,n2,1,'real*4',2); |
130 |
fld2((n2+1):(2*n2),:)=quikread_llc(f2,n2,1,'real*4',5); |
131 |
fld2((2*n2+1):(3*n2),:)=rot90(quikread_llc(f2,n2,1,'real*4',9),2); |
132 |
fld2((3*n2+1):(4*n2),:)=rot90(quikread_llc(f2,n2,1,'real*4',12),2); |
133 |
|
134 |
p3='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run/'; |
135 |
n3=4320; |
136 |
t3=26928; |
137 |
f3=[p3 'Theta.' myint2str(t3,10) '.data']; |
138 |
fld3=zeros(n3*4,n3); |
139 |
fld3(1:n3,:)=quikread_llc(f3,n3,1,'real*4',2); |
140 |
fld3((n3+1):(2*n3),:)=quikread_llc(f3,n3,1,'real*4',5); |
141 |
fld3((2*n3+1):(3*n3),:)=rot90(quikread_llc(f3,n3,1,'real*4',9),2); |
142 |
fld3((3*n3+1):(4*n3),:)=rot90(quikread_llc(f3,n3,1,'real*4',12),2); |
143 |
|
144 |
cx=[0 30]; |
145 |
clf reset |
146 |
colormap(cmap) |
147 |
subplot(311) |
148 |
mypcolor(fld1'); |
149 |
caxis(cx) |
150 |
thincolorbar |
151 |
title('llc_1080 SST, 17 September 2011') |
152 |
subplot(312) |
153 |
mypcolor(fld2'); |
154 |
caxis(cx) |
155 |
thincolorbar |
156 |
title('llc_2160 SST, 16 September 2011') |
157 |
subplot(313) |
158 |
mypcolor(fld3'); |
159 |
caxis(cx) |
160 |
thincolorbar |
161 |
title('llc_4320 SST, 17 September 2011') |
162 |
|
163 |
print -dpsc -r600 Theta_17sep2011 |
164 |
|
165 |
%%%%%%%%%%%% |
166 |
nx=4320; |
167 |
p1='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run_day1/'; |
168 |
pn='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run/'; |
169 |
nm='Theta'; |
170 |
f=2; iy=1300:1900; ix=2600:3400; |
171 |
clf reset |
172 |
colormap(cmap) |
173 |
ts=180; |
174 |
fn=[p1 nm '.' myint2str(ts,10) '.data']; |
175 |
fld=quikread_llc(fn,nx,1,'real*4',f); |
176 |
subplot(211) |
177 |
mypcolor(fld(ix,iy)'); |
178 |
caxis([10 20]) |
179 |
thincolorbar |
180 |
ts=14976; |
181 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
182 |
fld=quikread_llc(fn,nx,1,'real*4',f); |
183 |
subplot(212) |
184 |
mypcolor(fld(ix,iy)'); |
185 |
caxis([10 20]) |
186 |
thincolorbar |
187 |
|
188 |
print -dpsc -r600 Theta |
189 |
|
190 |
%%%%%%%%%%%% |
191 |
pn='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run/'; |
192 |
nm='Theta'; |
193 |
nx=4320; |
194 |
clf reset |
195 |
colormap(cmap) |
196 |
ts=136368; |
197 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
198 |
fld=quikread_llc(fn,nx); |
199 |
fld(find(~fld))=-3; |
200 |
quikplot_llc(fld) |
201 |
caxis([-3 31]) |
202 |
axis([0 17281 1600 15121]) |
203 |
title(['llc4320 ' ts2dte(ts,45,2011,1,17) ', Sea Surface Temperature (deg C)']) |
204 |
thincolorbar |
205 |
eval(['print -dpsc -r600 Theta4320_' int2str(ts)]) |
206 |
|
207 |
%%%%%%%%%%%% |
208 |
nm='KPPhbl'; |
209 |
n=1; |
210 |
clf reset |
211 |
colormap(cmap) |
212 |
for ts=92240:80:140720; |
213 |
mydisp(ts) |
214 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
215 |
fld=quikread_llc(fn,nx); |
216 |
clf |
217 |
quikplot_llc(log10(fld)) |
218 |
caxis([0 3]) |
219 |
axis([0 8641 800 7561]) |
220 |
title('llc2160, log10(PBL)') |
221 |
thincolorbar |
222 |
text(1000,1000,ts2dte(ts,45,2011,1,17),'color','w') |
223 |
eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]); |
224 |
n=n+1; |
225 |
end |
226 |
|
227 |
%%%%%%%%%%%% |
228 |
nm='Eta'; |
229 |
n=1; |
230 |
clf reset |
231 |
for ts=92240:80:140720; |
232 |
mydisp(ts) |
233 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
234 |
fld=quikread_llc(fn,nx); |
235 |
clf |
236 |
quikplot_llc(fld) |
237 |
caxis([-1 1]*3) |
238 |
axis([0 8641 800 7561]) |
239 |
title('llc2160, sea surface height (m)') |
240 |
thincolorbar |
241 |
text(1000,1000,ts2dte(ts,45,2011,1,17)) |
242 |
eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]); |
243 |
n=n+1; |
244 |
end |
245 |
|
246 |
%%%%%%%%%%%% |
247 |
nm='Theta'; |
248 |
n=1; |
249 |
clf reset |
250 |
colormap(cmap) |
251 |
ts=4464; |
252 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
253 |
fld=quikread_llc(fn,nx); |
254 |
clf |
255 |
quikplot_llc(fld) |
256 |
caxis([-2.5 30.5]) |
257 |
axis([0 17281 1600 15121]) |
258 |
title('llc4320, sea surface temperature (deg C)') |
259 |
thincolorbar |
260 |
print -dpsc -r600 Theta4320_4464 |
261 |
|
262 |
%%%%%%%%%%%% |
263 |
nm='W'; |
264 |
kx=68; |
265 |
n=1; |
266 |
clf reset |
267 |
colormap(cmap) |
268 |
for ts=92240:80:140720; |
269 |
mydisp(ts) |
270 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
271 |
fld=quikread_llc(fn,nx,kx); |
272 |
clf |
273 |
quikplot_llc(fld*1000) |
274 |
caxis([-1 1]*8) |
275 |
axis([0 8641 800 7561]) |
276 |
title('llc2160, Wvel at 2 km depth (mm/s)') |
277 |
thincolorbar |
278 |
text(1000,1000,ts2dte(ts,45,2011,1,17),'color','k') |
279 |
eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]); |
280 |
n=n+1; |
281 |
end |
282 |
|
283 |
%%%%%%%%%%%% |
284 |
% W2000 in EqWPac |
285 |
nm='W'; |
286 |
kx=68; |
287 |
n=1; |
288 |
clf reset |
289 |
colormap(cmap) |
290 |
for ts=92240:80:140720; |
291 |
mydisp(ts) |
292 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
293 |
fld=quikread_llc(fn,nx,kx,'real*4',9); |
294 |
clf |
295 |
quikpcolor(1000*rot90(fld,2)') |
296 |
caxis([-1 1]*5) |
297 |
title('llc2160, Wvel at 2 km depth (mm/s)') |
298 |
thincolorbar |
299 |
text(620,540,ts2dte(ts,45,2011,1,17),'color','k') |
300 |
eval(['print -dtiff -r135 ' pn 'figs/' nm 'WPac' myint2str(n,4)]); |
301 |
n=n+1; |
302 |
end |
303 |
|
304 |
%%%%%%%%%%%% |
305 |
nm='Salt'; |
306 |
n=1; |
307 |
clf reset |
308 |
colormap(cmap) |
309 |
for ts=92240:80:140720; |
310 |
mydisp(ts) |
311 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
312 |
fld=quikread_llc(fn,nx); |
313 |
clf |
314 |
quikplot_llc(fld) |
315 |
caxis([28 38]) |
316 |
axis([0 8641 800 7561]) |
317 |
title('llc2160, sea surface salinity (g/kg)') |
318 |
thincolorbar |
319 |
text(1000,1000,ts2dte(ts,45,2011,1,17),'color','w') |
320 |
eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]); |
321 |
n=n+1; |
322 |
end |
323 |
|
324 |
%%%%%%%%%%%% |
325 |
nm='PhiBot'; |
326 |
fld1=0; |
327 |
n=1; |
328 |
for ts=92240:80:140720; |
329 |
mydisp(ts) |
330 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
331 |
fld1=fld1+quikread_llc(fn,nx); |
332 |
n=n+1; |
333 |
end |
334 |
fld1=fld1/n; |
335 |
n=1; |
336 |
clf reset |
337 |
colormap(cmap) |
338 |
for ts=92240:80:140720; |
339 |
mydisp(ts) |
340 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
341 |
fld=quikread_llc(fn,nx)-fld1; |
342 |
clf |
343 |
quikplot_llc(fld/9.81) |
344 |
caxis([-2 2]) |
345 |
axis([0 8641 800 7561]) |
346 |
title('llc2160, bottom pressure (m)') |
347 |
thincolorbar |
348 |
text(1000,1000,ts2dte(ts,45,2011,1,17)) |
349 |
eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]); |
350 |
n=n+1; |
351 |
end |
352 |
|
353 |
%%%%%%%%%%%% |
354 |
% Eta perturbation |
355 |
nm='Eta'; |
356 |
fld1=0; |
357 |
fld2=0; |
358 |
n=1; |
359 |
for ts=92240:80:140720; |
360 |
mydisp(ts) |
361 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
362 |
fld=quikread_llc(fn,nx); |
363 |
fld1=fld1+fld; |
364 |
fld2=fld2+fld.^2; |
365 |
n=n+1; |
366 |
end |
367 |
fld1=fld1/n; |
368 |
fld2=fld2/n; |
369 |
writebin([pn 'EtaMean.data'],fld1); |
370 |
writebin([pn 'EtaMeanSquare.data'],fld2); |
371 |
clf reset |
372 |
colormap(cmap) |
373 |
quikplot_llc(sqrt(fld2-fld1.^2)*2*sqrt(2)) |
374 |
caxis([0 3]) |
375 |
thincolorbar |
376 |
axis([0 8641 800 7561]) |
377 |
title('Global Eta range (m): 2*sqrt(2)*std(Eta)') |
378 |
eval(['print -dtiff -r135 ' pn 'figs/' nm 'Range']); |
379 |
n=1; |
380 |
for ts=92240:80:140720; |
381 |
mydisp(ts) |
382 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
383 |
fld=quikread_llc(fn,nx)-fld1; |
384 |
clf |
385 |
quikplot_llc(fld) |
386 |
caxis([-1 1]*2) |
387 |
axis([0 8641 800 7561]) |
388 |
title('llc2160, Eta perturbation (m)') |
389 |
thincolorbar |
390 |
text(1000,1000,ts2dte(ts,45,2011,1,17)) |
391 |
eval(['print -dtiff -r135 ' pn 'figs/' nm 'Pert' myint2str(n,4)]); |
392 |
n=n+1; |
393 |
end |
394 |
|
395 |
%%%%%%%%%%%% |
396 |
% Theta perturbation |
397 |
nm='Theta'; |
398 |
fld1=0; |
399 |
n=1; |
400 |
for ts=92240:80:140720; |
401 |
mydisp(ts) |
402 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
403 |
fld1=fld1+quikread_llc(fn,nx); |
404 |
n=n+1; |
405 |
end |
406 |
fld1=fld1/n; |
407 |
n=1; |
408 |
clf reset |
409 |
colormap(cmap) |
410 |
for ts=92240:80:140720; |
411 |
mydisp(ts) |
412 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
413 |
fld=quikread_llc(fn,nx)-fld1; |
414 |
clf |
415 |
quikplot_llc(fld) |
416 |
caxis([-1 1]*1.5) |
417 |
axis([0 8641 800 7561]) |
418 |
title('llc2160, Theta perturbation (deg C)') |
419 |
thincolorbar |
420 |
text(1000,1000,ts2dte(ts,45,2011,1,17)) |
421 |
eval(['print -dtiff -r135 ' pn 'figs/' nm 'Pert' myint2str(n,4)]); |
422 |
n=n+1; |
423 |
end |
424 |
|
425 |
%%%%%%%%%%%% |
426 |
% Theta 12-hour difference |
427 |
nm='Theta'; |
428 |
n=1; |
429 |
clf reset |
430 |
colormap(cmap) |
431 |
for ts=93200:80:140720; |
432 |
mydisp(ts) |
433 |
fn1=[pn nm '.' myint2str(ts-960,10) '.data']; |
434 |
fn2=[pn nm '.' myint2str(ts,10) '.data']; |
435 |
fld=quikread_llc(fn2,nx)-quikread_llc(fn1,nx); |
436 |
clf |
437 |
quikplot_llc(fld) |
438 |
caxis([-1 1]) |
439 |
axis([0 8641 800 7561]) |
440 |
title('llc2160, Theta 12-hour difference (deg C)') |
441 |
thincolorbar |
442 |
text(1000,1000,ts2dte(ts,45,2011,1,17)) |
443 |
eval(['print -dtiff -r135 ' pn 'figs/' nm 'Diff' myint2str(n,4)]); |
444 |
n=n+1; |
445 |
end |
446 |
|
447 |
%%%%%%%%%%%% |
448 |
nm1='U'; nm2='V'; |
449 |
clf reset |
450 |
colormap(cmap) |
451 |
ts=140720; |
452 |
fn1=[pn nm1 '.' myint2str(ts,10) '.data']; |
453 |
fn2=[pn nm2 '.' myint2str(ts,10) '.data']; |
454 |
fl1=quikread_llc(fn1,nx); |
455 |
fl2=quikread_llc(fn2,nx); |
456 |
fld=sqrt(fl1.^2+fl2.^2); |
457 |
clf |
458 |
quikplot_llc(fld) |
459 |
caxis([0 1]) |
460 |
axis([0 8641 800 7561]) |
461 |
title('llc2160, ocean surface speed (m/s)') |
462 |
thincolorbar |
463 |
text(1000,1000,ts2dte(ts,45,2011,1,17),'color','w') |
464 |
eval(['print -dpsc ' pn 'figs/HiResSpeed']); |
465 |
|
466 |
%%%%%%%%%%%% |
467 |
nm1='U'; nm2='V'; |
468 |
n=1; |
469 |
clf reset |
470 |
colormap(cmap) |
471 |
for ts=92240:80:140720; |
472 |
mydisp(ts) |
473 |
fn1=[pn nm1 '.' myint2str(ts,10) '.data']; |
474 |
fn2=[pn nm2 '.' myint2str(ts,10) '.data']; |
475 |
fl1=quikread_llc(fn1,nx); |
476 |
fl2=quikread_llc(fn2,nx); |
477 |
fld=sqrt(fl1.^2+fl2.^2); |
478 |
clf |
479 |
quikplot_llc(fld) |
480 |
caxis([0 1]) |
481 |
axis([0 8641 800 7561]) |
482 |
title('llc2160, ocean surface speed (m/s)') |
483 |
thincolorbar |
484 |
text(1000,1000,ts2dte(ts,45,2011,1,17),'color','w') |
485 |
eval(['print -dtiff -r135 ' pn 'figs/Speed' myint2str(n,4)]); |
486 |
n=n+1; |
487 |
end |
488 |
|
489 |
%%%%%%%%%%%% |
490 |
% circum-Greenland speed |
491 |
nm1='U'; nm2='V'; |
492 |
n=1; |
493 |
clf reset |
494 |
orient landscape |
495 |
wysiwyg |
496 |
colormap(cmap) |
497 |
for ts=92240:80:140720; |
498 |
mydisp(ts) |
499 |
fn1=[pn nm1 '.' myint2str(ts,10) '.data']; |
500 |
fn2=[pn nm2 '.' myint2str(ts,10) '.data']; |
501 |
fl1=sqrt(quikread_llc(fn1,nx,1,'real*4',7).^2+quikread_llc(fn2,nx,1,'real*4',7).^2); |
502 |
fl2=sqrt(quikread_llc(fn1,nx,1,'real*4',11).^2+quikread_llc(fn2,nx,1,'real*4',11).^2); |
503 |
fl3=sqrt(quikread_llc(fn1,nx,1,'real*4',3).^2+quikread_llc(fn2,nx,1,'real*4',3).^2); |
504 |
fld=zeros(nx*2,'single'); |
505 |
fld(1:nx,(nx+1):end)=rot90(fl1,2); |
506 |
fld(1:nx,1:nx)=rot90(fl2,2); |
507 |
fld((nx+1):end,(nx+1):end)=rot90(fl1,1); |
508 |
fld((nx+1):end,1:nx)=fl3; |
509 |
clf |
510 |
quikpcolor(fld') |
511 |
caxis([0 1]*.8) |
512 |
axis([940 3080 1490 2620]) |
513 |
title('llc2160, ocean surface speed (m/s)') |
514 |
thincolorbar |
515 |
text(1900,2200,ts2dte(ts,45,2011,1,17),'color','w') |
516 |
eval(['print -dtiff -r135 ' pn 'figs/GreenSpeed' myint2str(n,4)]); |
517 |
n=n+1; |
518 |
end |
519 |
|
520 |
%%%%%%%%%%%% |
521 |
% circum-Greenland speed, 15-m |
522 |
nm1='U'; nm2='V'; |
523 |
k=9; % 14.68 m depth |
524 |
n=1; |
525 |
clf reset |
526 |
orient landscape |
527 |
wysiwyg |
528 |
colormap(cmap) |
529 |
for ts=92240:80:140720; |
530 |
mydisp(ts) |
531 |
fn1=[pn nm1 '.' myint2str(ts,10) '.data']; |
532 |
fn2=[pn nm2 '.' myint2str(ts,10) '.data']; |
533 |
fl1=sqrt(quikread_llc(fn1,nx,k,'real*4',7).^2+quikread_llc(fn2,nx,k,'real*4',7).^2); |
534 |
fl2=sqrt(quikread_llc(fn1,nx,k,'real*4',11).^2+quikread_llc(fn2,nx,k,'real*4',11).^2); |
535 |
fl3=sqrt(quikread_llc(fn1,nx,k,'real*4',3).^2+quikread_llc(fn2,nx,k,'real*4',3).^2); |
536 |
fld=zeros(nx*2,'single'); |
537 |
fld(1:nx,(nx+1):end)=rot90(fl1,2); |
538 |
fld(1:nx,1:nx)=rot90(fl2,2); |
539 |
fld((nx+1):end,(nx+1):end)=rot90(fl1,1); |
540 |
fld((nx+1):end,1:nx)=fl3; |
541 |
clf |
542 |
quikpcolor(fld') |
543 |
caxis([0 1]*.8) |
544 |
axis([940 3080 1490 2620]) |
545 |
title('llc2160, current speed at 15 m depth (m/s)') |
546 |
thincolorbar |
547 |
text(1900,2200,ts2dte(ts,45,2011,1,17),'color','w') |
548 |
eval(['print -dtiff -r135 ' pn 'figs/GreenSpeed15m' myint2str(n,4)]); |
549 |
n=n+1; |
550 |
end |
551 |
|
552 |
%%%%%%%%%%%% |
553 |
nm1='oceTAUX'; nm2='oceTAUY'; |
554 |
n=1; |
555 |
clf reset |
556 |
colormap(cmap) |
557 |
for ts=92240:80:140720; |
558 |
mydisp(ts) |
559 |
fn1=[pn nm1 '.' myint2str(ts,10) '.data']; |
560 |
fn2=[pn nm2 '.' myint2str(ts,10) '.data']; |
561 |
fl1=quikread_llc(fn1,nx); |
562 |
fl2=quikread_llc(fn2,nx); |
563 |
fld=sqrt(fl1.^2+fl2.^2); |
564 |
clf |
565 |
quikplot_llc(fld) |
566 |
caxis([0 1]/2) |
567 |
axis([0 8641 800 7561]) |
568 |
title('llc2160, ocean surface stress (Pa)') |
569 |
thincolorbar |
570 |
text(1000,1000,ts2dte(ts,45,2011,1,17),'color','w') |
571 |
eval(['print -dtiff -r135 ' pn 'figs/Stress' myint2str(n,4)]); |
572 |
n=n+1; |
573 |
end |
574 |
|
575 |
%%%%%%%%%%%% |
576 |
nm='SIheff'; |
577 |
n=1; |
578 |
clf reset |
579 |
colormap(cmap) |
580 |
for ts=92240:80:140720; |
581 |
mydisp(ts) |
582 |
fn=[pn nm '.' myint2str(ts,10) '.data']; |
583 |
fld=quikread_llc(fn,nx,1,'real*4',7); |
584 |
clf |
585 |
quikpcolor(rot90(fld,2)') |
586 |
caxis([0 10]) |
587 |
title('llc2160, effective sea ice thickness (m)') |
588 |
thincolorbar |
589 |
text(1500,400,ts2dte(ts,45,2011,1,17),'color','w') |
590 |
eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]); |
591 |
n=n+1; |
592 |
end |
593 |
|
594 |
%%%%%%%%%%%% |
595 |
ts=140704; |
596 |
pn='/nobackupp5/dmenemen/llc_2160/MITgcm/run_crash/'; |
597 |
fn=[pn 'Theta.' myint2str(ts,10) '.data']; |
598 |
for k=1:90 |
599 |
fld=quikread_llc(fn,nx); |
600 |
disp([k minmax(fld)]) |
601 |
end |