/[MITgcm]/MITgcm_contrib/gael/comm/course-idma2016/matlab/mitgcm_plot_plume.m
ViewVC logotype

Annotation of /MITgcm_contrib/gael/comm/course-idma2016/matlab/mitgcm_plot_plume.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.3 - (hide annotations) (download)
Fri Jan 22 20:30:55 2016 UTC (9 years, 6 months ago) by gforget
Branch: MAIN
CVS Tags: HEAD
Changes since 1.2: +5 -5 lines
- revise example
- switch to gcmfaces_sphere if available

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    

  ViewVC Help
Powered by ViewVC 1.1.22