/[MITgcm]/manual/s_getstarted/text/getting_started.tex
ViewVC logotype

Diff of /manual/s_getstarted/text/getting_started.tex

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.12 by adcroft, Wed Dec 5 15:49:39 2001 UTC revision 1.14 by edhill, Wed Jul 30 13:42:52 2003 UTC
# Line 28  and site, as well as links for downloadi Line 28  and site, as well as links for downloadi
28  documentation, to data-sources and other related sites.  documentation, to data-sources and other related sites.
29    
30  There is also a support news group for the model that you can email at  There is also a support news group for the model that you can email at
31  \texttt{support@mitgcm.org} or browse at:  \texttt{MITgcm-support@mitgcm.org} or browse at:
32  \begin{verbatim}  \begin{verbatim}
33  news://mitgcm.org/mitgcm.support  news://mitgcm.org/mitgcm.support
34  \end{verbatim}  \end{verbatim}
# Line 41  in the future. Line 41  in the future.
41    
42  MITgcm can be downloaded from our system by following  MITgcm can be downloaded from our system by following
43  the instructions below. As a courtesy we ask that you send e-mail to us at  the instructions below. As a courtesy we ask that you send e-mail to us at
44  \begin{rawhtml} <A href=mailto:support@mitgcm.org> \end{rawhtml}  \begin{rawhtml} <A href=mailto:MITgcm-support@mitgcm.org> \end{rawhtml}
45  support@mitgcm.org  MITgcm-support@mitgcm.org
46  \begin{rawhtml} </A> \end{rawhtml}  \begin{rawhtml} </A> \end{rawhtml}
47  to enable us to keep track of who's using the model and in what application.  to enable us to keep track of who's using the model and in what application.
48  You can download the model two ways:  You can download the model two ways:
# Line 330  calculation. Described in detail in sect Line 330  calculation. Described in detail in sect
330  \item \textit{global\_ocean.90x40x15} Global circulation with  \item \textit{global\_ocean.90x40x15} Global circulation with
331  GM, flux boundary conditions and poles.  GM, flux boundary conditions and poles.
332    
333    \item \textit{global\_ocean\_pressure} Global circulation in pressure
334      coordinate (non-Boussinesq ocean model). Described in detail in
335      section \ref{sect:eg-globalpressure}.
336    
337  \item \textit{solid-body.cs-32x32x1} Solid body rotation test for cube sphere  \item \textit{solid-body.cs-32x32x1} Solid body rotation test for cube sphere
338  grid.  grid.
339    
# Line 867  here I think. To come soon...) Line 871  here I think. To come soon...)
871  \item time-discretization  \item time-discretization
872  \end{itemize}  \end{itemize}
873    
874  The time steps are set through the real variables \textbf{deltaTMom }and  The time steps are set through the real variables \textbf{deltaTMom}
875  \textbf{deltaTtracer }(in s) which represent the time step for the momentum  and \textbf{deltaTtracer} (in s) which represent the time step for the
876  and tracer equations, respectively. For synchronous integrations, simply set  momentum and tracer equations, respectively. For synchronous
877  the two variables to the same value (or you can prescribe one time step only  integrations, simply set the two variables to the same value (or you
878  through the variable \textbf{deltaT}). The Adams-Bashforth stabilizing  can prescribe one time step only through the variable
879  parameter is set through the variable \textbf{abEps }(dimensionless). The  \textbf{deltaT}). The Adams-Bashforth stabilizing parameter is set
880  stagger baroclinic time stepping can be activated by setting the logical  through the variable \textbf{abEps} (dimensionless). The stagger
881  variable \textbf{staggerTimeStep }to '.\texttt{TRUE}.'.  baroclinic time stepping can be activated by setting the logical
882    variable \textbf{staggerTimeStep} to '.\texttt{TRUE}.'.
883    
884  \subsection{Equation of state}  \subsection{Equation of state}
885    
886  First, because the model equations are written in terms of perturbations, a  First, because the model equations are written in terms of
887  reference thermodynamic state needs to be specified. This is done through  perturbations, a reference thermodynamic state needs to be specified.
888  the 1D arrays \textbf{tRef}\textit{\ }and \textbf{sRef}. \textbf{tRef }%  This is done through the 1D arrays \textbf{tRef} and \textbf{sRef}.
889  specifies the reference potential temperature profile (in $^{o}$C for  \textbf{tRef} specifies the reference potential temperature profile
890  the ocean and $^{o}$K for the atmosphere) starting from the level  (in $^{o}$C for the ocean and $^{o}$K for the atmosphere) starting
891  k=1. Similarly, \textbf{sRef}\textit{\ }specifies the reference salinity  from the level k=1. Similarly, \textbf{sRef} specifies the reference
892  profile (in ppt) for the ocean or the reference specific humidity profile  salinity profile (in ppt) for the ocean or the reference specific
893  (in g/kg) for the atmosphere.  humidity profile (in g/kg) for the atmosphere.
894    
895  The form of the equation of state is controlled by the character variables  The form of the equation of state is controlled by the character
896  \textbf{buoyancyRelation}\textit{\ }and \textbf{eosType}\textit{. }\textbf{%  variables \textbf{buoyancyRelation} and \textbf{eosType}.
897  buoyancyRelation}\textit{\ }is set to '\texttt{OCEANIC}' by default and  \textbf{buoyancyRelation} is set to '\texttt{OCEANIC}' by default and
898  needs to be set to '\texttt{ATMOSPHERIC}' for atmosphere simulations. In  needs to be set to '\texttt{ATMOSPHERIC}' for atmosphere simulations.
899  this case, \textbf{eosType}\textit{\ }must be set to '\texttt{IDEALGAS}'.  In this case, \textbf{eosType} must be set to '\texttt{IDEALGAS}'.
900  For the ocean, two forms of the equation of state are available: linear (set  For the ocean, two forms of the equation of state are available:
901  \textbf{eosType}\textit{\ }to '\texttt{LINEAR}') and a polynomial  linear (set \textbf{eosType} to '\texttt{LINEAR}') and a polynomial
902  approximation to the full nonlinear equation ( set \textbf{eosType}\textit{\  approximation to the full nonlinear equation ( set
903  }to '\texttt{POLYNOMIAL}'). In the linear case, you need to specify the  \textbf{eosType}\textit{\ }to '\texttt{POLYNOMIAL}'). In the linear
904  thermal and haline expansion coefficients represented by the variables  case, you need to specify the thermal and haline expansion
905  \textbf{tAlpha}\textit{\ }(in K$^{-1}$) and \textbf{sBeta}\textit{\ }(in ppt$%  coefficients represented by the variables \textbf{tAlpha}\textit{\
906  ^{-1}$). For the nonlinear case, you need to generate a file of polynomial    }(in K$^{-1}$) and \textbf{sBeta} (in ppt$^{-1}$). For the nonlinear
907  coefficients called \textit{POLY3.COEFFS. }To do this, use the program  case, you need to generate a file of polynomial coefficients called
908  \textit{utils/knudsen2/knudsen2.f }under the model tree (a Makefile is  \textit{POLY3.COEFFS}. To do this, use the program
909  available in the same directory and you will need to edit the number and the  \textit{utils/knudsen2/knudsen2.f} under the model tree (a Makefile is
910  values of the vertical levels in \textit{knudsen2.f }so that they match  available in the same directory and you will need to edit the number
911  those of your configuration). \textit{\ }  and the values of the vertical levels in \textit{knudsen2.f} so that
912    they match those of your configuration).
913    
914    There there are also higher polynomials for the equation of state:
915    \begin{description}
916    \item['\texttt{UNESCO}':] The UNESCO equation of state formula of
917      Fofonoff and Millard \cite{fofonoff83}. This equation of state
918      assumes in-situ temperature, which is not a model variable; \emph{its use
919      is therefore discouraged, and it is only listed for completeness}.
920    \item['\texttt{JMD95Z}':] A modified UNESCO formula by Jackett and
921      McDougall \cite{jackett95}, which uses the model variable potential
922      temperature as input. The '\texttt{Z}' indicates that this equation
923      of state uses a horizontally and temporally constant pressure
924      $p_{0}=-g\rho_{0}z$.
925    \item['\texttt{JMD95P}':] A modified UNESCO formula by Jackett and
926      McDougall \cite{jackett95}, which uses the model variable potential
927      temperature as input. The '\texttt{P}' indicates that this equation
928      of state uses the actual hydrostatic pressure of the last time
929      step. Lagging the pressure in this way requires an additional pickup
930      file for restarts.
931    \item['\texttt{MDJWF}':] The new, more accurate and less expensive
932      equation of state by McDougall et~al. \cite{mcdougall03}. It also
933      requires lagging the pressure and therefore an additional pickup
934      file for restarts.
935    \end{description}
936    For none of these options an reference profile of temperature or
937    salinity is required.
938    
939  \subsection{Momentum equations}  \subsection{Momentum equations}
940    
# Line 1136  fields can be written out by setting the Line 1167  fields can be written out by setting the
1167  The precision with which to write the binary data is controlled by the  The precision with which to write the binary data is controlled by the
1168  integer variable w\textbf{riteBinaryPrec }(set it to \texttt{32} or \texttt{%  integer variable w\textbf{riteBinaryPrec }(set it to \texttt{32} or \texttt{%
1169  64}).  64}).
1170    
1171    %%% Local Variables:
1172    %%% mode: latex
1173    %%% TeX-master: t
1174    %%% End:
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22