1 |
gforget |
1.1 |
function []=mitgcm_plot_plume(fil,cc); |
2 |
gforget |
1.2 |
%Example: |
3 |
gforget |
1.3 |
% dirExp='MITgcm/verification/tutorial_plume_on_slope/run/'; |
4 |
|
|
% grid_load(dirExp,1,'straight'); |
5 |
gforget |
1.1 |
% cc=[-0.1 0.02]; |
6 |
gforget |
1.3 |
% figure; |
7 |
gforget |
1.1 |
% for tt=300:300:8640; |
8 |
gforget |
1.3 |
% fil=sprintf('%sT.%010d',dirExp,tt); |
9 |
gforget |
1.1 |
% mitgcm_plot_plume(fil,cc); |
10 |
|
|
% pause(0.1); |
11 |
|
|
% end; |
12 |
gforget |
1.2 |
% |
13 |
|
|
%note: this requires that the gcmfaces software has been |
14 |
|
|
% installed and added to the matlab path as explained in |
15 |
gforget |
1.3 |
% MITgcm_contrib/gael/matlab_class/gcmfaces.pdf |
16 |
gforget |
1.1 |
|
17 |
|
|
gcmfaces_global; |
18 |
|
|
|
19 |
|
|
nx=320; nr=60; |
20 |
|
|
x=squeeze(mygrid.XC{1})*ones(1,nr); |
21 |
|
|
y=ones(nx,1)*mygrid.RC'; |
22 |
|
|
m=squeeze(mygrid.mskC{1}); |
23 |
|
|
|
24 |
|
|
c=rdmds2gcmfaces(fil); |
25 |
|
|
c=squeeze(c{1}); |
26 |
|
|
|
27 |
|
|
% figureL; |
28 |
|
|
depthStretchPlot('pcolor',{x,y,m.*c}); shading flat; |
29 |
|
|
if isempty(whos('cc')); caxis(cc); end; |
30 |
|
|
colorbar; title(fil); |
31 |
|
|
|