s_examples/global_oce_biogeo/biogeochem.tex

\section{Biogeochemistry Tutorial}
%\label{www:tutorials}
\label{sec:eg-biogeochem_tutorial}
\begin{rawhtml}
<!-- CMIREDIR:eg-biogeochem_tutorial: -->
\end{rawhtml}
\begin{center}
(in directory: {\it verification/tutorial\_global\_oce\_biogeo/})
\end{center}

\subsection{Overview}
This model overlays the dissolved inorganic carbon biogeochemistry
model (the ``dic'' package) over a 2.8$^o$ global physical model. The
physical model has 15 levels, and is forced with a climatological
annual cycle of surface wind stresses \cite{Trenberth_etal_89},
surface heat and freshwater fluxes \cite{jiang99} with additional
relaxation toward climatological sea surface temperature and salinity
\cite{lev:94a,Levitus94}.  It uses the Gent and and McWilliams
\cite{gen-mcw:90} eddy parameterization scheme, has an implicit
free-surface, implicit vertical diffusion and uses the convective
adjustment scheme. The files for this experiment can be found in
the verification directory under tutorial\_global\_oce\_biogeo.

The biogeochemical model considers the coupled cycles of carbon,
oxygen, phosphorus and alkalinity.  A simplified parameterization of
biological production is used, limited by the availability of light
and phosphate.  A fraction of this productivity enters the dissolved
organic pool pool, which has an e-folding timescale for
remineralization of 6 months [following Yamanaka and Tajika, 1997].
The remaining fraction of this productivity is instantaneously
exported as particulate to depth [Yamanaka and Tajika, 1997] where it
is remineralized according to the empirical power law relationship
determined by Martin et al [1987].  The fate of carbon is linked to
that of phosphorus by the Redfield ratio. Carbonate chemistry is
explicitly solved \cite{Follows_etal_05} and the air-sea exchange of
CO$_2$ is parameterized with a uniform gas transfer coefficient
following \cite{Wanninkhof_92}. Oxygen is also linked to phosphorus by
the Redfield ratio, and oxygen air-sea exchange also follows
\cite{Wanninkhof_92}.  For more details see \cite{Dutkiewicz_etal_05}.

The example setup described here shows the physical model after 5900
years of spin-up and the biogeochemistry after 2900 years of spin-up.
The biogeochemistry is at a pre-industrial steady-state (atmospheric
ppmv is kept at 278). Five tracers are resolved: dissolved inorganic
carbon ($DIC$), alkalinity ($ALK$), phosphate ($PO4$), dissolved
organic phosphorus ($DOP$) and dissolved oxygen ($O2$).

\begin{figure} [tpb]
\begin{center}
\includegraphics[width=\textwidth,height=.3\textheight]{s_examples/global_oce_biogeo/co2flux.eps}
\caption{Modelled annual  mean air-sea CO$_2$ flux (mol C m$^{-2}$ y$^{-1}$)
for pre-industrial steady-state. Positive indicates  flux of CO$_2$
from ocean to the atmosphere (out-gassing),
contour interval is 1 mol C m$^{-2}$ y$^{-1}$.}
\label{lFcarflux}
\end{center}
\end{figure}


\subsection{Equations Solved}

The physical ocean model velocity and diffusivities are used to
redistribute the 5 tracers within the ocean. Additional redistribution
comes from chemical and biological sources and sinks.  For any tracer
$A$:
\begin{equation}
  \frac{\partial A}{\partial t}=-\nabla \cdot (\vec{u^{*}} A)+\nabla \cdot
  (\mathbf{K}\nabla A)+S_A \nonumber \label{lEtrac}
\end{equation}
where $\vec{u^{*}}$ is the transformed Eulerian mean circulation
(which includes Eulerian and eddy-induced advection), $\mathbf{K}$ is
the mixing tensor, and $S_A$ are the sources and sinks due to
biological and chemical processes.

The sources and sinks are:
\begin{eqnarray}
S_{DIC} & = &  F_{CO_2} + V_{CO_2} + r_{C:P} S_{PO_4}  + J_{Ca} \label{lEsdic} \\
S_{ALK} & = &  V_{ALK}-r_{N:P} S_{PO_4}  + 2 J_{Ca} \label{lEsalk} \\
S_{PO_4}& = &  -f_{DOP} J_{prod} - \frac{\partial F_P}{\partial z} + \kappa_{remin} [DOP]\\
S_{DOP} & = &  f_{DOP} J_{prod} -\kappa_{remin} [DOP] \\
S_{O_2} & = & \left\{ \begin{array}{ll}
               -r_{O:P} S_{PO_4} & \mbox{if $O_2>O_{2crit}$} \\
                0  & \mbox{if $O_2<O_{2crit}$}
                      \end{array}
              \right. 
\end{eqnarray}
where:
\begin{itemize}
\item $F_{CO_2}$ is the flux of CO$_2$ from the ocean to the
  atmosphere
\item $V_{CO_2}$ is ``virtual flux'' due to changes in $DIC$ due to
  the surface freshwater fluxes
\item $r_{C:P}$ is Redfield ratio of carbon to phosphorus
\item $J_{Ca}$ includes carbon removed from surface due to calcium
  carbonate formation and subsequent cumulation of the downward flux
  of CaCO$_3$
\item $V_{ALK}$ is ``virtual flux'' due to changes in alkalinity due
  to the surface freshwater fluxes
\item $r_{N:P}$ Redfield ratio is nitrogen to phosphorus
\item $f_{DOP}$ is fraction of productivity that remains suspended in
  the water column as dissolved organic phosphorus
\item $J_{prod}$ is the net community productivity
\item $\frac{\partial F_P}{\partial z}$ is the accumulation of
  remineralized phosphorus with depth
\item $\kappa_{remin}$ is rate with which $DOP$ remineralizes back to
  $PO_4$
\item $F_{O_2}$ is air-sea flux of oxygen
\item $r_{O:P}$ is Redfield ratio of oxygen to phosphorus
\item $O_{2crit}$ is a critical level below which oxygen consumption
  if halted
\end{itemize}

These terms (for the first four tracers) are described more in
\cite{Dutkiewicz_etal_05} and by \cite{McKinley_etal_04} for the terms
relating to oxygen.


\subsection{Code configuration}

The model configuration for this experiment resides in
verification/dic\_example. The modifications to the code (in {\it
  verification/dic\_example/code}) are:
\begin{itemize}
\item{{\bf SIZE.h}: which dictates the size of the model domain
    (128x64x15).}
\item{\bf PTRACERS\_SIZE.h}: which dictates how many tracers to assign
  how many tracers will be used (5).
\item{\bf GCHEM\_OPTIONS.h}: provides some compiler time options for
  the {\it /pkg/gchem}. In particular this example requires that {\it
    DIC\_BIOTIC} and {\it GCHEM\_SEPARATE\_FORCING} be defined.
\item{\bf GMREDI\_OPTIONS.h}: assigns the Gent-McWilliam eddy
  parameterization options.
\item{\bf DIAGNOSTICS\_SIZE.h}: assigns size information for the
  diagnostics package.
\item{\bf packages.conf}: which dictates which packages will be
  compiled in this version of the model - among the many that are used
  for the physical part of the model, this also includes {\it
    ptracers}, {\it gchem}, and {\it dic} which allow the
  biogeochemical part of this setup to function.
\end{itemize}

\vspace{1cm}
\noindent
The input fields needed for this run (in {\it
  verification/dic\_example/input}) are:
\begin{itemize}
\item {\bf data}: specifies the main parameters for the experiment,
  some parameters that may be useful to know: {\it nTimeSteps} number
  timesteps model will run, change to 720 to run for a year {\it
    taveFreq} frequency with which time averages are done, change to
  31104000 for annual averages.
\item {\bf data.diagnostics}: species details of diagnostic pkg output
\item {\bf data.gchem}: specifics files and other details needed in
  the biogeochemistry model run
\item {\bf data.gmredi}: species details for the GM parameterization
\item {\bf data.mnc}: specifies details for types of output, netcdf or
  binary
\item {\bf data.pkg}: set true or false for various packages to be
  used
\item {\bf data.ptracers}: details of the tracers to be used,
  including makes, diffusivity information and (if needed) initial
  files. Of particular importance is the {\it PTRACERS\_numInUse}
  which states how many tracers are used, and {\it PTRACERS\_Iter0}
  which states at which timestep the biogeochemistry model tracers
  were initialized.
\item {\bf depth\_g77.bin}: bathymetry data file
\item {\bf eedata}: This file uses standard default values and does
  not contain customizations for this experiment.
\item {\bf fice.bin}: ice data file, needed for the biogeochemistry
\item {\bf lev\_monthly\_salt.bin}: SSS values which model relaxes
  toward
\item {\bf lev\_monthly\_temp.bin}: SST values which model relaxes
  toward
\item {\bf pickup.0004248000.data}: variable and tendency values need
  to restart the physical part of the model
\item {\bf pickup\_cd.0004248000.data}: variable and tendency values
  need to restart the cd pkg
\item {\bf pickup\_ptracers.0004248000.data}: variable and tendency
  values need to restart the the biogeochemistry part of the model
\item {\bf POLY3.COEFFS}: coefficient for the non-linear equation of
  state
\item {\bf shi\_empmr\_year.bin}: freshwater forcing data file
\item {\bf shi\_qnet.bin}: heat flux forcing data file
\item {\bf sillev1.bin}: silica data file, need for the
  biogeochemistry
\item {\bf tren\_speed.bin}: wind speed data file, needed for the
  biogeochemistry
\item {\bf tren\_taux.bin}: meridional wind stress data file
\item {\bf tren\_tauy.bin}: zonal wind stress data file
\end{itemize}


\subsection{Running the example}

You will first need to download the MITgcm code. Instructions for
downloading the code can be found in section 3.2.

\begin{enumerate}
\item{go to the build directory in verification/dic\_example:\\
    \hspace{1cm}  {\it cd verification/dic\_example/build}}
\item{create the Makefile:\\
    \hspace{1cm}  {\it ../../../tools/genmake2 -mods=code}}
\item{create all the links:\\
    \hspace{1cm}  {\it make depend}}
\item{compile (the executable will be called mitgcmuv):\\
    \hspace{1cm}  {\it make}}
\item{move the executable to the directory with all the inputs:\\
    \hspace{1cm}  {\it mv mitgcmuv ../input/}}
\item{go to the input directory and run the model:\\
    \hspace{1cm}  {\it cd ../input}\\
    \hspace{1cm}  {\it ./mitgcmuv}}
\end{enumerate}
As the model is set up to run in the verification experiment, it only
runs for 4 timestep (2 days) and outputs data at the end of this short
run. For a more informative run, you will need to run longer. As set
up, this model starts from a pre-spun up state and initializes
physical fields and the biogeochemical tracers from the {\it pickup}
files.

Physical data (e.g. S,T, velocities etc) will be output as for any
regular ocean run. The biogeochemical output are:
\begin{itemize}
\item tracer snap shots: either netcdf, or older-style binary
  (depending on how {\it data.mnc} is set up). Look in {\it
    data.ptracers} to see which number matches which type of tracer
  (e.g. ptracer01 is DIC).
\item tracer time averages: either netcdf, or older-style binary
  (depending on how {\it data.mnc} is set up)
\item specific DIC diagnostics: these are averaged over {\it taveFreq}
  (set in {\it data}) and are specific to the dic package, and
  currently are only available in binary format:
  \begin{itemize}
  \item{\bf DIC\_Biotave}: 3-D biological community productivity (mol
    P m$^{-3}$ s$^{-1}$)
  \item{\bf DIC\_Cartave}: 3-D tendencies due to calcium carbonate
    cycle (mol C m$^{-3}$ s$^{-1}$)
  \item{\bf DIC\_fluxCO2ave}: 2-D air-sea flux of CO$_2$ (mol C
    m$^{-2}$ s$^{-1}$)
  \item{\bf DIC\_pCO2tave}: 2-D partial pressure of CO$_2$ in surface
    layer
  \item{\bf DIC\_pHtave}: 2-D pH in surface layer
  \item{\bf DIC\_SurOtave}: 2-D tendency due to air-sea flux of O$_2$
    (mol O m$^{-3}$ s$^{-1}$)
  \item{\bf DIC\_Surtave}: 2-D surface tendency of DIC due to air-sea
    flux and virtual flux (mol C m$^{-3}$ s$^{-1}$)
  \end{itemize}
\end{itemize}


%% \subsection{Reference Material}

%% \Hpar
%% Dutkiewicz. S., A. Sokolov, J.Scott and P. Stone, 2005:
%% A Three-Dimensional Ocean-Seaice-Carbon Cycle Model and its Coupling
%% to a Two-Dimensional Atmospheric Model: Uses in Climate Change Studies,
%% Report 122, Joint Program of the Science and Policy of Global Change,
%% M.I.T., Cambridge, MA.\\
%% (http://web.mit.edu/globalchange/www/MITJPSPGC\_Rpt122.pdf)
%% \Hpar
%% Follows, M., T. Ito and S. Dutkiewicz, 2005:
%% A Compact and Accurate Carbonate Chemistry Solver for Ocean
%% Biogeochemistry Models. {\it Ocean Modeling}, in press.
%% \Hpar
%% Gent, P. and J. McWilliams, 1990:
%% Isopycnal mixing in ocean circulation models.
%% {\it Journal of Physical Oceanography}, 20, 150 -- 155.
%% \Hpar
%% Jiang, S., P.H. Stone, and P. Malanotte-Rizzoli,
%% An assessment of the Geophysical Fluid Dynamics Laboratory
%% ocean model with coarse resolution: Annual-mean climatology,
%% {\it Journal of Geophysical Research}, 104, 25623 -- 25645, 1999.
%% \Hpar
%% Levitus, S. and T.P. Boyer, 1994:
%% {\it World Ocean Atlas 1994 Volume 4: Temperature},
%% NOAA Atlas NESDIS 4, U.S. Department of Commerce,
%% Washington, D.C., 117pp.
%% \Hpar
%% Levitus, S., R. Burgett, and T.P. Boyer, 1994:
%% {\it World Ocean Atlas 1994 Volume 3: Salinity},
%% NOAA Atlas NESDIS 3, U.S. Department of Commerce,
%% Washington, D.C., 99pp.
%% \Hpar
%% McKinley, G., M.J. Follows and J.C. Marshall, 2004:
%% Mechanisms of air-sea CO$_2$ flux variability in the Equatorial Pacific
%% and the North Atlantic.
%% {\it Global Biogeochemical Cycles}, 18, doi:10.1029/2003GB002179.
%% \Hpar
%% Trenberth, K., J. Olson, and W. Large, 1989:
%% {\it A global wind stress climatology based on ECMWF analyses,
%% Tech. Rep. NCAR/TN-338+STR},
%% National Center for Atmospheric Research, Boulder, Colorado.
%% \Hpar
%% Wanninkhof, R., 1992:
%% Relationship between wind speed and gas exchange over the ocean,
%% {\it Journal of Geophysical Research}, 97, 7373 -- 7382.

1	\section{Biogeochemistry Tutorial}
2	%\label{www:tutorials}
3	\label{sec:eg-biogeochem_tutorial}
4	\begin{rawhtml}
5	<!-- CMIREDIR:eg-biogeochem_tutorial: -->
6	\end{rawhtml}
7	\begin{center}
8	(in directory: {\it verification/tutorial\_global\_oce\_biogeo/})
9	\end{center}
10
11	\subsection{Overview}
12	This model overlays the dissolved inorganic carbon biogeochemistry
13	model (the ``dic'' package) over a 2.8$^o$ global physical model. The
14	physical model has 15 levels, and is forced with a climatological
15	annual cycle of surface wind stresses \cite{Trenberth_etal_89},
16	surface heat and freshwater fluxes \cite{jiang99} with additional
17	relaxation toward climatological sea surface temperature and salinity
18	\cite{lev:94a,Levitus94}. It uses the Gent and and McWilliams
19	\cite{gen-mcw:90} eddy parameterization scheme, has an implicit
20	free-surface, implicit vertical diffusion and uses the convective
21	adjustment scheme. The files for this experiment can be found in
22	the verification directory under tutorial\_global\_oce\_biogeo.
23
24	The biogeochemical model considers the coupled cycles of carbon,
25	oxygen, phosphorus and alkalinity. A simplified parameterization of
26	biological production is used, limited by the availability of light
27	and phosphate. A fraction of this productivity enters the dissolved
28	organic pool pool, which has an e-folding timescale for
29	remineralization of 6 months [following Yamanaka and Tajika, 1997].
30	The remaining fraction of this productivity is instantaneously
31	exported as particulate to depth [Yamanaka and Tajika, 1997] where it
32	is remineralized according to the empirical power law relationship
33	determined by Martin et al [1987]. The fate of carbon is linked to
34	that of phosphorus by the Redfield ratio. Carbonate chemistry is
35	explicitly solved \cite{Follows_etal_05} and the air-sea exchange of
36	CO$_2$ is parameterized with a uniform gas transfer coefficient
37	following \cite{Wanninkhof_92}. Oxygen is also linked to phosphorus by
38	the Redfield ratio, and oxygen air-sea exchange also follows
39	\cite{Wanninkhof_92}. For more details see \cite{Dutkiewicz_etal_05}.
40
41	The example setup described here shows the physical model after 5900
42	years of spin-up and the biogeochemistry after 2900 years of spin-up.
43	The biogeochemistry is at a pre-industrial steady-state (atmospheric
44	ppmv is kept at 278). Five tracers are resolved: dissolved inorganic
45	carbon ($DIC$), alkalinity ($ALK$), phosphate ($PO4$), dissolved
46	organic phosphorus ($DOP$) and dissolved oxygen ($O2$).
47
48	\begin{figure} [tpb]
49	\begin{center}
50	\includegraphics[width=\textwidth,height=.3\textheight]{s_examples/global_oce_biogeo/co2flux.eps}
51	\caption{Modelled annual mean air-sea CO$_2$ flux (mol C m$^{-2}$ y$^{-1}$)
52	for pre-industrial steady-state. Positive indicates flux of CO$_2$
53	from ocean to the atmosphere (out-gassing),
54	contour interval is 1 mol C m$^{-2}$ y$^{-1}$.}
55	\label{lFcarflux}
56	\end{center}
57	\end{figure}
58
59
60	\subsection{Equations Solved}
61
62	The physical ocean model velocity and diffusivities are used to
63	redistribute the 5 tracers within the ocean. Additional redistribution
64	comes from chemical and biological sources and sinks. For any tracer
65	$A$:
66	\begin{equation}
67	\frac{\partial A}{\partial t}=-\nabla \cdot (\vec{u^{*}} A)+\nabla \cdot
68	(\mathbf{K}\nabla A)+S_A \nonumber \label{lEtrac}
69	\end{equation}
70	where $\vec{u^{*}}$ is the transformed Eulerian mean circulation
71	(which includes Eulerian and eddy-induced advection), $\mathbf{K}$ is
72	the mixing tensor, and $S_A$ are the sources and sinks due to
73	biological and chemical processes.
74
75	The sources and sinks are:
76	\begin{eqnarray}
77	S_{DIC} & = & F_{CO_2} + V_{CO_2} + r_{C:P} S_{PO_4} + J_{Ca} \label{lEsdic} \\
78	S_{ALK} & = & V_{ALK}-r_{N:P} S_{PO_4} + 2 J_{Ca} \label{lEsalk} \\
79	S_{PO_4}& = & -f_{DOP} J_{prod} - \frac{\partial F_P}{\partial z} + \kappa_{remin} [DOP]\\
80	S_{DOP} & = & f_{DOP} J_{prod} -\kappa_{remin} [DOP] \\
81	S_{O_2} & = & \left\{ \begin{array}{ll}
82	-r_{O:P} S_{PO_4} & \mbox{if $O_2>O_{2crit}$} \\
83	0 & \mbox{if $O_2<O_{2crit}$}
84	\end{array}
85	\right.
86	\end{eqnarray}
87	where:
88	\begin{itemize}
89	\item $F_{CO_2}$ is the flux of CO$_2$ from the ocean to the
90	atmosphere
91	\item $V_{CO_2}$ is ``virtual flux'' due to changes in $DIC$ due to
92	the surface freshwater fluxes
93	\item $r_{C:P}$ is Redfield ratio of carbon to phosphorus
94	\item $J_{Ca}$ includes carbon removed from surface due to calcium
95	carbonate formation and subsequent cumulation of the downward flux
96	of CaCO$_3$
97	\item $V_{ALK}$ is ``virtual flux'' due to changes in alkalinity due
98	to the surface freshwater fluxes
99	\item $r_{N:P}$ Redfield ratio is nitrogen to phosphorus
100	\item $f_{DOP}$ is fraction of productivity that remains suspended in
101	the water column as dissolved organic phosphorus
102	\item $J_{prod}$ is the net community productivity
103	\item $\frac{\partial F_P}{\partial z}$ is the accumulation of
104	remineralized phosphorus with depth
105	\item $\kappa_{remin}$ is rate with which $DOP$ remineralizes back to
106	$PO_4$
107	\item $F_{O_2}$ is air-sea flux of oxygen
108	\item $r_{O:P}$ is Redfield ratio of oxygen to phosphorus
109	\item $O_{2crit}$ is a critical level below which oxygen consumption
110	if halted
111	\end{itemize}
112
113	These terms (for the first four tracers) are described more in
114	\cite{Dutkiewicz_etal_05} and by \cite{McKinley_etal_04} for the terms
115	relating to oxygen.
116
117
118	\subsection{Code configuration}
119
120	The model configuration for this experiment resides in
121	verification/dic\_example. The modifications to the code (in {\it
122	verification/dic\_example/code}) are:
123	\begin{itemize}
124	\item{{\bf SIZE.h}: which dictates the size of the model domain
125	(128x64x15).}
126	\item{\bf PTRACERS\_SIZE.h}: which dictates how many tracers to assign
127	how many tracers will be used (5).
128	\item{\bf GCHEM\_OPTIONS.h}: provides some compiler time options for
129	the {\it /pkg/gchem}. In particular this example requires that {\it
130	DIC\_BIOTIC} and {\it GCHEM\_SEPARATE\_FORCING} be defined.
131	\item{\bf GMREDI\_OPTIONS.h}: assigns the Gent-McWilliam eddy
132	parameterization options.
133	\item{\bf DIAGNOSTICS\_SIZE.h}: assigns size information for the
134	diagnostics package.
135	\item{\bf packages.conf}: which dictates which packages will be
136	compiled in this version of the model - among the many that are used
137	for the physical part of the model, this also includes {\it
138	ptracers}, {\it gchem}, and {\it dic} which allow the
139	biogeochemical part of this setup to function.
140	\end{itemize}
141
142	\vspace{1cm}
143	\noindent
144	The input fields needed for this run (in {\it
145	verification/dic\_example/input}) are:
146	\begin{itemize}
147	\item {\bf data}: specifies the main parameters for the experiment,
148	some parameters that may be useful to know: {\it nTimeSteps} number
149	timesteps model will run, change to 720 to run for a year {\it
150	taveFreq} frequency with which time averages are done, change to
151	31104000 for annual averages.
152	\item {\bf data.diagnostics}: species details of diagnostic pkg output
153	\item {\bf data.gchem}: specifics files and other details needed in
154	the biogeochemistry model run
155	\item {\bf data.gmredi}: species details for the GM parameterization
156	\item {\bf data.mnc}: specifies details for types of output, netcdf or
157	binary
158	\item {\bf data.pkg}: set true or false for various packages to be
159	used
160	\item {\bf data.ptracers}: details of the tracers to be used,
161	including makes, diffusivity information and (if needed) initial
162	files. Of particular importance is the {\it PTRACERS\_numInUse}
163	which states how many tracers are used, and {\it PTRACERS\_Iter0}
164	which states at which timestep the biogeochemistry model tracers
165	were initialized.
166	\item {\bf depth\_g77.bin}: bathymetry data file
167	\item {\bf eedata}: This file uses standard default values and does
168	not contain customizations for this experiment.
169	\item {\bf fice.bin}: ice data file, needed for the biogeochemistry
170	\item {\bf lev\_monthly\_salt.bin}: SSS values which model relaxes
171	toward
172	\item {\bf lev\_monthly\_temp.bin}: SST values which model relaxes
173	toward
174	\item {\bf pickup.0004248000.data}: variable and tendency values need
175	to restart the physical part of the model
176	\item {\bf pickup\_cd.0004248000.data}: variable and tendency values
177	need to restart the cd pkg
178	\item {\bf pickup\_ptracers.0004248000.data}: variable and tendency
179	values need to restart the the biogeochemistry part of the model
180	\item {\bf POLY3.COEFFS}: coefficient for the non-linear equation of
181	state
182	\item {\bf shi\_empmr\_year.bin}: freshwater forcing data file
183	\item {\bf shi\_qnet.bin}: heat flux forcing data file
184	\item {\bf sillev1.bin}: silica data file, need for the
185	biogeochemistry
186	\item {\bf tren\_speed.bin}: wind speed data file, needed for the
187	biogeochemistry
188	\item {\bf tren\_taux.bin}: meridional wind stress data file
189	\item {\bf tren\_tauy.bin}: zonal wind stress data file
190	\end{itemize}
191
192
193	\subsection{Running the example}
194
195	You will first need to download the MITgcm code. Instructions for
196	downloading the code can be found in section 3.2.
197
198	\begin{enumerate}
199	\item{go to the build directory in verification/dic\_example:\\
200	\hspace{1cm} {\it cd verification/dic\_example/build}}
201	\item{create the Makefile:\\
202	\hspace{1cm} {\it ../../../tools/genmake2 -mods=code}}
203	\item{create all the links:\\
204	\hspace{1cm} {\it make depend}}
205	\item{compile (the executable will be called mitgcmuv):\\
206	\hspace{1cm} {\it make}}
207	\item{move the executable to the directory with all the inputs:\\
208	\hspace{1cm} {\it mv mitgcmuv ../input/}}
209	\item{go to the input directory and run the model:\\
210	\hspace{1cm} {\it cd ../input}\\
211	\hspace{1cm} {\it ./mitgcmuv}}
212	\end{enumerate}
213	As the model is set up to run in the verification experiment, it only
214	runs for 4 timestep (2 days) and outputs data at the end of this short
215	run. For a more informative run, you will need to run longer. As set
216	up, this model starts from a pre-spun up state and initializes
217	physical fields and the biogeochemical tracers from the {\it pickup}
218	files.
219
220	Physical data (e.g. S,T, velocities etc) will be output as for any
221	regular ocean run. The biogeochemical output are:
222	\begin{itemize}
223	\item tracer snap shots: either netcdf, or older-style binary
224	(depending on how {\it data.mnc} is set up). Look in {\it
225	data.ptracers} to see which number matches which type of tracer
226	(e.g. ptracer01 is DIC).
227	\item tracer time averages: either netcdf, or older-style binary
228	(depending on how {\it data.mnc} is set up)
229	\item specific DIC diagnostics: these are averaged over {\it taveFreq}
230	(set in {\it data}) and are specific to the dic package, and
231	currently are only available in binary format:
232	\begin{itemize}
233	\item{\bf DIC\_Biotave}: 3-D biological community productivity (mol
234	P m$^{-3}$ s$^{-1}$)
235	\item{\bf DIC\_Cartave}: 3-D tendencies due to calcium carbonate
236	cycle (mol C m$^{-3}$ s$^{-1}$)
237	\item{\bf DIC\_fluxCO2ave}: 2-D air-sea flux of CO$_2$ (mol C
238	m$^{-2}$ s$^{-1}$)
239	\item{\bf DIC\_pCO2tave}: 2-D partial pressure of CO$_2$ in surface
240	layer
241	\item{\bf DIC\_pHtave}: 2-D pH in surface layer
242	\item{\bf DIC\_SurOtave}: 2-D tendency due to air-sea flux of O$_2$
243	(mol O m$^{-3}$ s$^{-1}$)
244	\item{\bf DIC\_Surtave}: 2-D surface tendency of DIC due to air-sea
245	flux and virtual flux (mol C m$^{-3}$ s$^{-1}$)
246	\end{itemize}
247	\end{itemize}
248
249
250	%% \subsection{Reference Material}
251
252	%% \Hpar
253	%% Dutkiewicz. S., A. Sokolov, J.Scott and P. Stone, 2005:
254	%% A Three-Dimensional Ocean-Seaice-Carbon Cycle Model and its Coupling
255	%% to a Two-Dimensional Atmospheric Model: Uses in Climate Change Studies,
256	%% Report 122, Joint Program of the Science and Policy of Global Change,
257	%% M.I.T., Cambridge, MA.\\
258	%% (http://web.mit.edu/globalchange/www/MITJPSPGC\_Rpt122.pdf)
259	%% \Hpar
260	%% Follows, M., T. Ito and S. Dutkiewicz, 2005:
261	%% A Compact and Accurate Carbonate Chemistry Solver for Ocean
262	%% Biogeochemistry Models. {\it Ocean Modeling}, in press.
263	%% \Hpar
264	%% Gent, P. and J. McWilliams, 1990:
265	%% Isopycnal mixing in ocean circulation models.
266	%% {\it Journal of Physical Oceanography}, 20, 150 -- 155.
267	%% \Hpar
268	%% Jiang, S., P.H. Stone, and P. Malanotte-Rizzoli,
269	%% An assessment of the Geophysical Fluid Dynamics Laboratory
270	%% ocean model with coarse resolution: Annual-mean climatology,
271	%% {\it Journal of Geophysical Research}, 104, 25623 -- 25645, 1999.
272	%% \Hpar
273	%% Levitus, S. and T.P. Boyer, 1994:
274	%% {\it World Ocean Atlas 1994 Volume 4: Temperature},
275	%% NOAA Atlas NESDIS 4, U.S. Department of Commerce,
276	%% Washington, D.C., 117pp.
277	%% \Hpar
278	%% Levitus, S., R. Burgett, and T.P. Boyer, 1994:
279	%% {\it World Ocean Atlas 1994 Volume 3: Salinity},
280	%% NOAA Atlas NESDIS 3, U.S. Department of Commerce,
281	%% Washington, D.C., 99pp.
282	%% \Hpar
283	%% McKinley, G., M.J. Follows and J.C. Marshall, 2004:
284	%% Mechanisms of air-sea CO$_2$ flux variability in the Equatorial Pacific
285	%% and the North Atlantic.
286	%% {\it Global Biogeochemical Cycles}, 18, doi:10.1029/2003GB002179.
287	%% \Hpar
288	%% Trenberth, K., J. Olson, and W. Large, 1989:
289	%% {\it A global wind stress climatology based on ECMWF analyses,
290	%% Tech. Rep. NCAR/TN-338+STR},
291	%% National Center for Atmospheric Research, Boulder, Colorado.
292	%% \Hpar
293	%% Wanninkhof, R., 1992:
294	%% Relationship between wind speed and gas exchange over the ocean,
295	%% {\it Journal of Geophysical Research}, 97, 7373 -- 7382.
296