--- MITgcm_contrib/gmaze_pv/compute_QEk.m	2006/07/10 15:09:01	1.1
+++ MITgcm_contrib/gmaze_pv/compute_QEk.m	2007/02/01 17:02:02	1.2
@@ -1,5 +1,5 @@
 %
-% [] = compute_QEk(SNAPSHOT)
+% [QEk] = compute_QEk(SNAPSHOT)
 %
 % Here we compute the lateral heat flux induced by Ekman currents
 % from JFz, the PV flux induced by frictional forces:
@@ -29,7 +29,7 @@
 % 06/27/06
 % gmaze@mit.edu
 
-function QEk = compute_QEk(snapshot)
+function varargout = compute_QEk(snapshot)
 
 global sla toshow
 global netcdf_suff netcdf_domain
@@ -151,3 +151,10 @@
 nc=close(nc);
 
 
+
+% Output:
+output = struct('QEk',QEk,'lat',JFzlat,'lon',JFzlon);
+switch nargout
+ case 1
+  varargout(1) = {output};
+end