--- MITgcm_contrib/gmaze_pv/compute_JFz.m 2006/10/06 21:44:52 1.2 +++ MITgcm_contrib/gmaze_pv/compute_JFz.m 2007/09/19 15:37:38 1.5 @@ -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