next up previous contents
Next: 39.7 local_potential_density_terms Up: 39. Diagnostics for physical Previous: 39.5 fct_netcdf

   
39.6 gyre_components

Option gyre_components computes instantaneous values of various components of northward tracer transport. The longitudinally averaged temperature and meridional velocity are constructed for the latitude of as a function of depth. Also, the vertically averaged temperature and meridional velocity are constructed for the latitude of as a function of longitude


$\displaystyle <\overline{T}^\phi_{k,j}>^\lambda$ = $\displaystyle \frac{1}{Volx^T_{k,j}}\sum_{i=2}^{imt-1} \overline{t_{i,k,j,n,\tau}}^{\phi}\Delta^T_i$ (39.33)
$\displaystyle <V_{k,j}>^\lambda$ = $\displaystyle \frac{1}{Volx^U_{k,j}}\sum_{i=2}^{imt-1} u_{i,k,j,2,\tau}\Delta^U_i$ (39.34)
$\displaystyle <\overline{T}^\phi_{i,j}>^z$ = $\displaystyle \frac{1}{Volz^T_{i,j}}\sum_{k=1}^{km} \overline{t_{i,k,j,n,\tau}}^{\phi}\Delta^T_k$ (39.35)
<Vi,j>z = $\displaystyle \frac{1}{Volz^U_{i,j}}\sum_{k=1}^{km} adv\_vnt_{i,k,j}\Delta^U_k$ (39.36)

Note that a factor of is built into $adv\_vnt_{i,k,j}$. The volume elements and volume of the latitude strip as a function of depth are


$\displaystyle \Delta^T_i$ = $\displaystyle tmask_{i,k,j}\cdot tmask_{i,k,j+1}\cdot\dxti$ (39.37)
$\displaystyle \Delta^U_i$ = $\displaystyle dxui\cdot \csuj$ (39.38)
$\displaystyle \Delta^T_k$ = $\displaystyle tmask_{i,k,j}\cdot tmask_{i,k,j+1}\cdot dzt_k$ (39.39)
$\displaystyle \Delta^U_k$ = dztk (39.40)
VolxTk,j = $\displaystyle \sum_{i=2}^{imt-1}\Delta^T_i$ (39.41)
VolxUk,j = $\displaystyle \sum_{i=2}^{imt-1}\Delta^U_i$ (39.42)
VolzTi,j = $\displaystyle \sum_{k=1}^{km}\Delta^T_k$ (39.43)
VolzUi,j = $\displaystyle \sum_{k=1}^{km}\Delta^U_k$ (39.44)

The canonical form of the northward components of tracer transport by various means and deviations is given as


ttn1,jrow,n = $\displaystyle \sum_{k=1}^{km} <\overline{T}^\phi_{k,j}>^\lambda \cdot <V_{k,j}>^\lambda\cdot dzt_k$ (39.45)
ttn2,jrow,n = ttn6,jrow,n - ttn1,jrow,n (39.46)
ttn3,jrow,n = $\displaystyle \sum_{i=2}^{imt-1} <\overline{T}^\phi_{i,j}>^z \cdot <V_{i,j}>^z\cdot \dxti$ (39.47)
ttn4,jrow,n = ttn6,jrow,n-ttn3,jrow,n-ttn5,jrow,n (39.48)
ttn5,jrow,n = $\displaystyle \sum_{i=2}^{imt-1} -(\overline{smf_{i-1,j,1}\cdot\dxuim}^\lambda)...
..._{i,1,j,n,\tau}}^\phi - <\overline{T}^\phi_{i,j}>^z)\cdot\frac{\csuj}{f_{jrow}}$ (39.49)
ttn6,jrow,n = $\displaystyle \sum_{k=1}^{km}\sum_{i=2}^{imt-1} 0.5\cdot adv\_fn_{i,k,j}\cdot \Delta^T_i\cdot dzt_k$ (39.50)
ttn7,jrow,n = $\displaystyle \sum_{k=1}^{km}\sum_{i=2}^{imt-1} diff\_fn_{i,k,j}\cdot \Delta^T_i\cdot dzt_k$ (39.51)
ttn8,jrow,n = ttn6,jrow,n + ttn7,jrow,n (39.52)

Note the factor of 0.5 which is needed to correct the advective flux of tracer as described in Section 22.8.2. These terms may also be broken down as a function of latitude within .

The output from this diagnostic may be written as ascii to the model printout or as 32 bit IEEE unformatted data to file gyre_comp.yyyyyy.mm.dd.dta. If option netcdf or gyre_components_netcdf is enabled, data is written in NetCDF format to file gyre_comp.yyyyyy.mm.dd.dta.nc rather than in unformatted IEEE. The ``yyyyyy.mm.dd'' is a place holder for year, month, and day and this naming convention is explained further in Section 38.2. The interval between output is specified by variable gyreint and the control is specified by variable iogyre.


next up previous contents
Next: 39.7 local_potential_density_terms Up: 39. Diagnostics for physical Previous: 39.5 fct_netcdf
RC Pacanowski and SM Griffies, GFDL, Jan 2000