1 |
% |
% |
2 |
% [] = compute_JFzx(SNAPSHOT) |
% [JFzx] = compute_JFzx(SNAPSHOT) |
3 |
% |
% |
4 |
% Here we compute the PV flux due to the zonal frictionnal force as |
% Here we compute the PV flux due to the zonal frictionnal force as |
5 |
% JFzx = ( TAUx * dSIGMATHETA/dy ) / RHO / EKL |
% JFzx = ( TAUx * dSIGMATHETA/dy ) / RHO / EKL |
24 |
% 06/04/12 |
% 06/04/12 |
25 |
% gmaze@mit.edu |
% gmaze@mit.edu |
26 |
|
|
27 |
function compute_JFzx(snapshot) |
function varargout = compute_JFzx(snapshot) |
28 |
|
|
29 |
global sla toshow |
global sla toshow |
30 |
global netcdf_suff netcdf_domain |
global netcdf_suff netcdf_domain |
171 |
close(ncRHO); |
close(ncRHO); |
172 |
close(ncH); |
close(ncH); |
173 |
close(nc); |
close(nc); |
174 |
|
|
175 |
|
% Output: |
176 |
|
output = struct('JFzx',JFz,'lat',STlat(2:ny-1),'lon',STlon(2:nx-1)); |
177 |
|
switch nargout |
178 |
|
case 1 |
179 |
|
varargout(1) = {output}; |
180 |
|
end |