function []=mitgcm_plot_held_suarez(fil,cc); %Example: % dirExp='MITgcm/verification/tutorial_held_suarez_cs/run/'; % grid_load(dirExp,6,'cube'); % cc=[260:2:310]; suf=[dirExp 'T']; rec=1; % figure; set(gcf,'Renderer','zbuffer'); % for tt=282240:192*30:345600; % fil={sprintf('%s.%010d',suf,tt),rec}; % mitgcm_plot_held_suarez(fil,cc); % shading interp; % pause(0.1); % end; % %note: this requires that the gcmfaces software has been % installed and added to the matlab path as explained in % MITgcm_contrib/gael/matlab_class/gcmfaces.pdf typePlot=0; if ~isempty(which('gcmfaces_sphere')); typePlot=1; end; gcmfaces_global; m=squeeze(mygrid.mskC(:,:,1)); if ~iscell(fil); c=rdmds2gcmfaces(fil); else; c=rdmds2gcmfaces(fil{1},'rec',fil{2}); end; c=c(:,:,1); clf; if isempty(whos('cc')); cc=[]; end; if typePlot==1; gcmfaces_sphere(c.*m,cc,[],[],1); else; qwckplot(c.*m); gcmfaces_cmap_cbar(cc); end; title(fil);