Diff of /MITgcm_contrib/gmaze_pv/compute_JFz.m

Parent Directory | Revision Log | Revision Graph | Patch

--- MITgcm_contrib/gmaze_pv/compute_JFz.m	2006/10/06 21:44:52	1.2
+++ MITgcm_contrib/gmaze_pv/compute_JFz.m	2007/02/01 17:02:02	1.3
@@ -1,5 +1,5 @@
 %
-% [] = compute_JFz(SNAPSHOT)
+% [JFz] = compute_JFz(SNAPSHOT)
 %
 % Here we compute the PV flux due to frictionnal forces as
 % JFz = ( TAUx * dSIGMATHETA/dy - TAUy * dSIGMATHETA/dx ) / RHO / EKL
@@ -25,7 +25,7 @@
 % 06/27/06
 % gmaze@mit.edu
 
-function compute_JFz(snapshot)
+function varargout = compute_JFz(snapshot)
 
 global sla toshow
 global netcdf_suff netcdf_domain
@@ -216,3 +216,9 @@
 nc=close(nc);
 
 
+% Output:
+output = struct('JFz',JFz,'lat',STlat(2:ny-1),'lon',STlon(2:nx-1));
+switch nargout
+ case 1
+  varargout(1) = {output};
+end