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