--- MITgcm_contrib/gmaze_pv/compute_JFzx.m	2007/01/30 22:10:10	1.1
+++ MITgcm_contrib/gmaze_pv/compute_JFzx.m	2007/02/01 17:02:02	1.2
@@ -1,5 +1,5 @@
 %
-% [] = compute_JFzx(SNAPSHOT)
+% [JFzx] = compute_JFzx(SNAPSHOT)
 %
 % Here we compute the PV flux due to the zonal frictionnal force as
 % JFzx = ( TAUx * dSIGMATHETA/dy ) / RHO / EKL
@@ -24,7 +24,7 @@
 % 06/04/12
 % gmaze@mit.edu
 
-function compute_JFzx(snapshot)
+function varargout = compute_JFzx(snapshot)
 
 global sla toshow
 global netcdf_suff netcdf_domain
@@ -171,3 +171,10 @@
 close(ncRHO);
 close(ncH);
 close(nc);
+
+% Output:
+output = struct('JFzx',JFz,'lat',STlat(2:ny-1),'lon',STlon(2:nx-1));
+switch nargout
+ case 1
+  varargout(1) = {output};
+end