/[MITgcm]/MITgcm/verification/tutorial_global_oce_latlon/diags_matlab/mit_loadglobal.m
ViewVC logotype

Diff of /MITgcm/verification/tutorial_global_oce_latlon/diags_matlab/mit_loadglobal.m

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

revision 1.3 by jmc, Sat Aug 12 20:25:13 2006 UTC revision 1.4 by heimbach, Thu Sep 15 22:12:09 2016 UTC
# Line 105  mit_plotstreamfunctions Line 105  mit_plotstreamfunctions
105  mit_plotzonalvelocity  mit_plotzonalvelocity
106    
107  if meanfields  if meanfields
108      % plot some mean fields
109      figure('PaperPosition',[0.31 0.25 10.5 7.88],'PaperOrientation','landscape')
110      k=kmax;
111      mySST =   t(:,:,1,k);
112      mySSS =   s(:,:,1,k);
113      myETA = eta(:,:,k);
114      myVEL = sqrt(u(:,:,1,k).*u(:,:,1,k)+v(:,:,1,k).*v(:,:,1,k));
115    %
116      clear sh;
117      sh(1) = subplot(2,2,1);imagesc(grd.lonc,grd.latc,mySST'.*grd.hfacc(:,:,1)');
118      title(['Sea Surface Temperature [degC]'])
119      sh(2) = subplot(2,2,2);imagesc(grd.lonc,grd.latc,mySSS'.*grd.hfacc(:,:,1)');
120      title(['Sea Surface Salinity [psu]'])
121      sh(3) = subplot(2,2,3);imagesc(grd.lonc,grd.latc,myETA');
122      title(['Sea Surface Height [m]'])
123      sh(4) = subplot(2,2,4);imagesc(grd.lonc,grd.latc,myVEL'.*grd.hfacc(:,:,1)');
124      title(['Sea Surface Speed [m/s]'])
125    %
126      axis(sh,'image'); axis(sh,'xy')
127      set(gcf,'currentAxes',sh(1));colorbar('h')
128      set(gcf,'currentAxes',sh(2));colorbar('h')
129      set(gcf,'currentAxes',sh(3));colorbar('h')
130      set(gcf,'currentAxes',sh(4));colorbar('h')
131    end
132    
133    if meanfields
134    % compute actual heat and E-P fluxes from time averages    % compute actual heat and E-P fluxes from time averages
135    % this only makes sense when you load the time averaged fields    % this only makes sense when you load the time averaged fields
136    tauTheta = mit_getparm('data','tauThetaClimRelax');    tauTheta = mit_getparm('data','tauThetaClimRelax');
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22