/[MITgcm]/MITgcm_contrib/gmaze_pv/compute_MLD.m
ViewVC logotype

Diff of /MITgcm_contrib/gmaze_pv/compute_MLD.m

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

revision 1.2 by gmaze, Fri Oct 6 21:44:53 2006 UTC revision 1.3 by gmaze, Thu Feb 1 17:02:02 2007 UTC
# Line 1  Line 1 
1  %  %
2  % [] = compute_MLD(SNAPSHOT)  % [MLD] = compute_MLD(SNAPSHOT)
3  %  %
4  % Here we compute the Mixed Layer Depth as:  % Here we compute the Mixed Layer Depth as:
5  % MLD = min depth for which : ST > ST(SSS,SST-0.8,p0)    % MLD = min depth for which : ST > ST(SSS,SST-0.8,p0)  
# Line 28  Line 28 
28  % 09/20/06  % 09/20/06
29  % gmaze@mit.edu  % gmaze@mit.edu
30    
31  function compute_MLD(snapshot)  function varargout = compute_MLD(snapshot)
32    
33  global sla toshow  global sla toshow
34  global netcdf_suff netcdf_domain  global netcdf_suff netcdf_domain
# Line 155  nc{ncid}.uniquename    = ncchar(uniquena Line 155  nc{ncid}.uniquename    = ncchar(uniquena
155  nc{ncid}(:,:,:)        = MLD;  nc{ncid}(:,:,:)        = MLD;
156    
157  nc=close(nc);  nc=close(nc);
158    close(ncST);
159    close(ncS);
160    close(ncT);
161    
162    
163    % Output:
164    output = struct('MLD',MLD,'lat',DTlat,'lon',STlon);
165    switch nargout
166     case 1
167      varargout(1) = {output};
168    end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22