/[MITgcm]/MITgcm_contrib/articles/ceaice/ceaice.tex
ViewVC logotype

Diff of /MITgcm_contrib/articles/ceaice/ceaice.tex

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

revision 1.7 by dimitri, Fri Jan 18 02:37:22 2008 UTC revision 1.10 by dimitri, Mon Feb 25 16:50:56 2008 UTC
# Line 1  Line 1 
1    % $Header$
2    % $Name$
3  \documentclass[12pt]{article}  \documentclass[12pt]{article}
4    
5  \usepackage[]{graphicx}  \usepackage[]{graphicx}
# Line 50  Line 52 
52  \maketitle  \maketitle
53    
54  \begin{abstract}  \begin{abstract}
55    Some blabla  
56    As part of ongoing efforts to obtain a best possible synthesis of most
57    available, global-scale, ocean and sea ice data, dynamic and thermodynamic
58    sea-ice model components have been incorporated in the Massachusetts Institute
59    of Technology general circulation model (MITgcm).  Sea-ice dynamics use either
60    a visco-plastic rheology solved with a line successive relaxation (LSR)
61    technique, reformulated on an Arakawa C-grid in order to match the oceanic and
62    atmospheric grids of the MITgcm, and modified to permit efficient and accurate
63    automatic differentiation of the coupled ocean and sea-ice model
64    configurations.
65    
66  \end{abstract}  \end{abstract}
67    
68  \section{Introduction}  \section{Introduction}
# Line 134  The maximum ice pressure $P_{\max}$, a m Line 146  The maximum ice pressure $P_{\max}$, a m
146  both thickness $h$ and compactness (concentration) $c$:  both thickness $h$ and compactness (concentration) $c$:
147  \begin{equation}  \begin{equation}
148    P_{\max} = P^{*}c\,h\,e^{[C^{*}\cdot(1-c)]},    P_{\max} = P^{*}c\,h\,e^{[C^{*}\cdot(1-c)]},
149  \label{icestrength}  \label{eq:icestrength}
150  \end{equation}  \end{equation}
151  with the constants $P^{*}$ and $C^{*}$. The nonlinear bulk and shear  with the constants $P^{*}$ and $C^{*}$. The nonlinear bulk and shear
152  viscosities $\eta$ and $\zeta$ are functions of ice strain rate  viscosities $\eta$ and $\zeta$ are functions of ice strain rate
# Line 342  five different cases at steady state (af Line 354  five different cases at steady state (af
354  \end{description}  \end{description}
355  \ml{[We have not implemented the EVP solver on a B-grid.]}  \ml{[We have not implemented the EVP solver on a B-grid.]}
356  \begin{figure*}[htbp]  \begin{figure*}[htbp]
357  %GET  \includegraphics[width=\widefigwidth]{\fpath/all_086280}    \includegraphics[width=\widefigwidth]{\fpath/all_086280}
358    \caption{Ice concentration, effective thickness [m], and ice    \caption{Ice concentration, effective thickness [m], and ice
359      velocities [m/s]      velocities [m/s]
360      for 5 different numerical solutions.}      for 5 different numerical solutions.}
# Line 412  contraint limits ice strength and viscos Line 424  contraint limits ice strength and viscos
424  time scale, resolution and EVP-time step, effectively allowing the  time scale, resolution and EVP-time step, effectively allowing the
425  elastic waves to damp out more quickly \citep{hunke01}.  elastic waves to damp out more quickly \citep{hunke01}.
426  \begin{figure*}[htbp]  \begin{figure*}[htbp]
427  %GET  \includegraphics[width=\widefigwidth]{\fpath/hun12days}    \includegraphics[width=\widefigwidth]{\fpath/hun12days}
428    \caption{Ice flow, divergence and bulk viscosities of three    \caption{Ice flow, divergence and bulk viscosities of three
429      experiments with \citet{hunke01}'s test case: C-LSRns (top),      experiments with \citet{hunke01}'s test case: C-LSRns (top),
430      C-EVPns (middle), and C-EVPns with damping described in      C-EVPns (middle), and C-EVPns with damping described in
# Line 435  These reduced viscosities lead to small Line 447  These reduced viscosities lead to small
447  which in turn can have a strong effect on solutions in the limit of  which in turn can have a strong effect on solutions in the limit of
448  nearly rigid regimes (arching and blocking, not shown).  nearly rigid regimes (arching and blocking, not shown).
449    
450    \ml{[Say something about performance? This is tricky, as the
451      perfomance depends strongly on the configuration. A run with slowly
452      changing forcing is favorable for LSR, because then only very few
453      iterations are required for convergences while EVP uses its fixed
454      number of internal timesteps. If the forcing in changing fast, LSR
455      needs far more iterations while EVP still uses the fixed number of
456      internal timesteps. I have produces runs where for slow forcing LSR
457      is much faster than EVP and for fast forcing, LSR is much slower
458      than EVP. EVP is certainly more efficient in terms of vectorization
459      and MFLOPS on our SX8, but is that a criterion?]}
460    
461  \subsection{C-grid}  \subsection{C-grid}
462  \begin{itemize}  \begin{itemize}
463  \item no-slip vs. free-slip for both lsr and evp;  \item no-slip vs. free-slip for both lsr and evp;
# Line 795  We thank Jinlun Zhang for providing the Line 818  We thank Jinlun Zhang for providing the
818  helpful discussions. ML thanks Elizabeth Hunke for multiple explanations.  helpful discussions. ML thanks Elizabeth Hunke for multiple explanations.
819    
820  \bibliography{bib/journal_abrvs,bib/seaice,bib/genocean,bib/maths,bib/mitgcmuv,bib/fram}  \bibliography{bib/journal_abrvs,bib/seaice,bib/genocean,bib/maths,bib/mitgcmuv,bib/fram}
 %\bibliography{journal_abrvs,seaice,genocean,maths,mixing,mitgcmuv,bib/fram}  
821    
822  \end{document}  \end{document}
823    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22