| 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 |