/[MITgcm]/MITgcm/utils/matlab/dens_poly3.m
ViewVC logotype

Diff of /MITgcm/utils/matlab/dens_poly3.m

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

revision 1.1.6.1 by heimbach, Mon Apr 8 20:10:39 2002 UTC revision 1.3 by jmc, Sat Feb 17 23:49:43 2007 UTC
# Line 16  function [dens] = dens_poly3(poly3,t,s) Line 16  function [dens] = dens_poly3(poly3,t,s)
16  % or to work within a single model level  % or to work within a single model level
17  % >> D=dens_poly3(P(3),T(:,:,3),S(:,:,3));  % >> D=dens_poly3(P(3),T(:,:,3),S(:,:,3));
18    
19    % $Header$
20    % $Name$
21    
22  if size(t) ~= size(s)  if size(t) ~= size(s)
23   error('T and S must be the same shape and size')   error('T and S must be the same shape and size')
24  end  end
# Line 51  for k=1:nz, Line 54  for k=1:nz,
54  end  end
55    
56  dens=reshape(dens,n);  dens=reshape(dens,n);
57    dens( find(t==0 & s==0) )=0;

Legend:
Removed from v.1.1.6.1  
changed lines
  Added in v.1.3

  ViewVC Help
Powered by ViewVC 1.1.22