/[MITgcm]/MITgcm_contrib/gmaze_pv/compute_JFz.m
ViewVC logotype

Diff of /MITgcm_contrib/gmaze_pv/compute_JFz.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by gmaze, Fri Oct 6 21:44:52 2006 UTC revision 1.5 by gmaze, Wed Sep 19 15:37:38 2007 UTC
# Line 1  Line 1 
1  %  %
2  % [] = compute_JFz(SNAPSHOT)  % [JFz] = compute_JFz(SNAPSHOT)
3  %  %
4  % Here we compute the PV flux due to frictionnal forces as  % Here we compute the PV flux due to frictionnal forces as
5  % JFz = ( TAUx * dSIGMATHETA/dy - TAUy * dSIGMATHETA/dx ) / RHO / EKL  % JFz = ( TAUx * dSIGMATHETA/dy - TAUy * dSIGMATHETA/dx ) / RHO / EKL
# Line 25  Line 25 
25  % 06/27/06  % 06/27/06
26  % gmaze@mit.edu  % gmaze@mit.edu
27    
28  function compute_JFz(snapshot)  function varargout = compute_JFz(snapshot)
29    
30  global sla toshow  global sla toshow
31  global netcdf_suff netcdf_domain  global netcdf_suff netcdf_domain
# Line 216  nc{ncid}(:,:,:)        = JFz; Line 216  nc{ncid}(:,:,:)        = JFz;
216  nc=close(nc);  nc=close(nc);
217    
218    
219    % Output:
220    output = struct('JFz',JFz,'lat',STlat(2:ny-1),'lon',STlon(2:nx-1));
221    switch nargout
222     case 1
223      varargout(1) = {output};
224    end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22