| 1 |
b150=loadgrid('/globec/flagg/realtime/bank150'); |
|---|
| 2 |
[a,b]=xy2ll(b150.x,b150.y); |
|---|
| 3 |
b150.x=b; |
|---|
| 4 |
b150.y=a; |
|---|
| 5 |
|
|---|
| 6 |
plotbnd(b150) |
|---|
| 7 |
drawelems(b150) |
|---|
| 8 |
lcontour(b150,'z',-[20 40 60 80 100]) |
|---|
| 9 |
|
|---|
| 10 |
b150.z=-b150.z; |
|---|
| 11 |
|
|---|
| 12 |
[m,n,p]=size(X); |
|---|
| 13 |
|
|---|
| 14 |
figure |
|---|
| 15 |
h1=slice(X,Y,Z,SAL,X(:,:,1),Y(:,:,1),Z(:,:,1),'linear');hold on |
|---|
| 16 |
h2=slice(X,Y,Z,SAL,X(:,:,1),Y(:,:,1),Z(:,:,4),'linear'); |
|---|
| 17 |
h3=slice(X,Y,Z,SAL,X(:,:,1),Y(:,:,1),Z(:,:,7),'linear'); |
|---|
| 18 |
h4=slice(X,Y,Z,SAL,X(:,:,1),Y(:,:,1),Z(:,:,10),'linear'); |
|---|
| 19 |
h5=slice(X,Y,Z,SAL,X(:,:,1),Y(:,:,1),Z(:,:,13),'linear'); |
|---|
| 20 |
h6=slice(X,Y,Z,SAL,X(:,:,1),Y(:,:,1),Z(:,:,16),'linear'); |
|---|
| 21 |
h7=slice(X,Y,Z,SAL,X(:,:,1),Y(:,:,1),Z(:,:,19),'linear'); |
|---|
| 22 |
h8=slice(X,Y,Z,SAL,X(:,:,1),Y(:,:,1),Z(:,:,22),'linear'); |
|---|
| 23 |
set(gca,'Clim',[32 33]) |
|---|
| 24 |
hp=colormesh3d(b150,b150.z) |
|---|
| 25 |
axis([min(X(:)) max(X(:)) min(Y(:)) max(Y(:)) -100 0]) |
|---|
| 26 |
colormap(jet(64)) |
|---|
| 27 |
colorbar |
|---|
| 28 |
rotate3d |
|---|
| 29 |
title('SALINITY','FontSize',20) |
|---|
| 30 |
xlabel('Long','FontSize',20);ylabel('Lat','FontSize',20) |
|---|
| 31 |
|
|---|
| 32 |
|
|---|
| 33 |
figure |
|---|
| 34 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,1),'linear');hold on |
|---|
| 35 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,3),'linear') |
|---|
| 36 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,5),'linear') |
|---|
| 37 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,7),'linear') |
|---|
| 38 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,9),'linear') |
|---|
| 39 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,11),'linear') |
|---|
| 40 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,13),'linear') |
|---|
| 41 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,15),'linear') |
|---|
| 42 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,17),'linear') |
|---|
| 43 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,19),'linear') |
|---|
| 44 |
slice(X,Y,Z,TEMP,X(:,:,1),Y(:,:,1),Z(:,:,21),'linear') |
|---|
| 45 |
hp=colormesh3d(b150,b150.z) |
|---|
| 46 |
axis([min(X(:)) max(X(:)) min(Y(:)) max(Y(:)) -100 0]) |
|---|
| 47 |
colormap(jet(64)) |
|---|
| 48 |
colorbar |
|---|
| 49 |
rotate3d |
|---|
| 50 |
title('TEMPERATURE','FontSize',20) |
|---|
| 51 |
xlabel('Long','FontSize',20);ylabel('Lat','FontSize',20) |
|---|
| 52 |
|
|---|
| 53 |
|
|---|
| 54 |
figure |
|---|
| 55 |
slice(X,Y,Z,TEMP,reshape(X(1,:,:),m,p),reshape(Y(1,:,:),m,p),reshape(Z(1,:,:),m,p));hold on |
|---|
| 56 |
slice(X,Y,Z,TEMP,reshape(X(1,:,:),m,p),reshape(Y(8,:,:),m,p),reshape(Z(1,:,:),m,p)) |
|---|
| 57 |
slice(X,Y,Z,TEMP,reshape(X(1,:,:),m,p),reshape(Y(15,:,:),m,p),reshape(Z(1,:,:),m,p)) |
|---|
| 58 |
slice(X,Y,Z,TEMP,reshape(X(1,:,:),m,p),reshape(Y(22,:,:),m,p),reshape(Z(1,:,:),m,p)) |
|---|
| 59 |
slice(X,Y,Z,TEMP,reshape(X(1,:,:),m,p),reshape(Y(29,:,:),m,p),reshape(Z(1,:,:),m,p)) |
|---|
| 60 |
slice(X,Y,Z,TEMP,reshape(X(1,:,:),m,p),reshape(Y(36,:,:),m,p),reshape(Z(1,:,:),m,p)) |
|---|
| 61 |
slice(X,Y,Z,TEMP,reshape(X(1,:,:),m,p),reshape(Y(43,:,:),m,p),reshape(Z(1,:,:),m,p)) |
|---|
| 62 |
title('TEMPERATURE','FontSize',20) |
|---|
| 63 |
xlabel('Long','FontSize',20);ylabel('Lat','FontSize',20) |
|---|
