--- MITgcm_contrib/gmaze_pv/compute_EKLx.m	2007/01/30 22:10:10	1.1
+++ MITgcm_contrib/gmaze_pv/compute_EKLx.m	2007/02/01 17:02:02	1.2
@@ -1,5 +1,5 @@
 %
-% [] = compute_EKLx(SNAPSHOT)
+% [EKL] = compute_EKLx(SNAPSHOT)
 %
 % Here we compute the Ekman Layer Depth as:
 % EKL = 0.7 sqrt( TAUx/RHO )/f 
@@ -23,7 +23,7 @@
 % 12/04/06
 % gmaze@mit.edu
 
-function compute_EKLx(snapshot)
+function varargout = compute_EKLx(snapshot)
 
 global sla toshow
 global netcdf_suff netcdf_domain
@@ -144,3 +144,10 @@
 close(nc);
 
 
+
+% Output:
+output = struct('EKL',EKL,'lat',RHOlat,'lon',RHOlon);
+switch nargout
+ case 1
+  varargout(1) = {output};
+end