1 |
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 |