--- MITgcm_contrib/gmaze_pv/compute_EKL.m	2006/10/06 19:23:50	1.1
+++ MITgcm_contrib/gmaze_pv/compute_EKL.m	2007/02/01 17:02:02	1.2
@@ -1,5 +1,5 @@
 %
-% [] = compute_EKL(SNAPSHOT)
+% [EKL] = compute_EKL(SNAPSHOT)
 %
 % Here we compute the Ekmal Layer Depth as:
 % EKL = 0.7 sqrt( |TAU|/RHO )/f 
@@ -24,7 +24,7 @@
 % 08/16/06
 % gmaze@mit.edu
 
-function compute_EKL(snapshot)
+function varargout = compute_EKL(snapshot)
 
 global sla toshow
 global netcdf_suff netcdf_domain
@@ -145,3 +145,10 @@
 nc=close(nc);
 
 
+
+% Output:
+output = struct('EKL',EKL,'lat',RHOlat,'lon',RHOlon);
+switch nargout
+ case 1
+  varargout(1) = {output};
+end