| 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 |