| 1 |
% |
% |
| 2 |
% [] = compute_EKL(SNAPSHOT) |
% [EKL] = compute_EKL(SNAPSHOT) |
| 3 |
% |
% |
| 4 |
% Here we compute the Ekmal Layer Depth as: |
% Here we compute the Ekmal Layer Depth as: |
| 5 |
% EKL = 0.7 sqrt( |TAU|/RHO )/f |
% EKL = 0.7 sqrt( |TAU|/RHO )/f |
| 24 |
% 08/16/06 |
% 08/16/06 |
| 25 |
% gmaze@mit.edu |
% gmaze@mit.edu |
| 26 |
|
|
| 27 |
function compute_EKL(snapshot) |
function varargout = compute_EKL(snapshot) |
| 28 |
|
|
| 29 |
global sla toshow |
global sla toshow |
| 30 |
global netcdf_suff netcdf_domain |
global netcdf_suff netcdf_domain |
| 145 |
nc=close(nc); |
nc=close(nc); |
| 146 |
|
|
| 147 |
|
|
| 148 |
|
|
| 149 |
|
% Output: |
| 150 |
|
output = struct('EKL',EKL,'lat',RHOlat,'lon',RHOlon); |
| 151 |
|
switch nargout |
| 152 |
|
case 1 |
| 153 |
|
varargout(1) = {output}; |
| 154 |
|
end |