1 |
adcroft |
1.2 |
load GRID |
2 |
|
|
|
3 |
|
|
nx=32;olx=5; |
4 |
|
|
j=0; |
5 |
|
|
|
6 |
|
|
cax=[1 1]; |
7 |
|
|
frame=0; |
8 |
|
|
%for t=48:48:192; |
9 |
|
|
for t=0:8:288; |
10 |
|
|
%for t=0:8:64; |
11 |
|
|
frame=frame+1; |
12 |
|
|
|
13 |
|
|
sf=sprintf('%10.10i',t); |
14 |
|
|
P=read_tiles(['Eta.' sf],nx*[1 1]); |
15 |
|
|
%T=read_tiles(['T.' sf],nx*[1 1]); |
16 |
|
|
%S=read_tiles(['S.' sf],nx*[1 1]); |
17 |
|
|
%W=read_tiles(['W.' sf],nx*[1 1]); |
18 |
|
|
%GT=read_tiles(['GT.' sf],nx*[1 1]); |
19 |
|
|
|
20 |
|
|
Q=P; |
21 |
|
|
|
22 |
|
|
%displaytiles(Q); |
23 |
|
|
|
24 |
|
|
%plotcube(XG,YG,ZG,Q) |
25 |
|
|
%axis square |
26 |
|
|
%axis off |
27 |
|
|
%view(130+0*t/256*360,20+0*sin(t/256*pi)); |
28 |
|
|
%axis([-1 1 -1 1 -1 1]*sqrt(2/3)) |
29 |
|
|
%cax=caxis; |
30 |
|
|
%caxis([-1 1]*max(abs(max(cax)),abs(min(cax)))) |
31 |
|
|
|
32 |
|
|
merccube(lonG,latG,Q);colorbar |
33 |
|
|
cax=caxis; |
34 |
|
|
caxis([-1 1]*max(abs(max(cax)),abs(min(cax))));colorbar |
35 |
|
|
|
36 |
|
|
title(sprintf('Time = %4.1f hours',t*450/3600)) |
37 |
|
|
drawnow |
38 |
|
|
|
39 |
|
|
%print('-djpeg60','-r60',sprintf('splash-%3.3i.jpg',frame)); |
40 |
|
|
end |