tenthalpy-basal-bc.tex - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
 (HTM) git clone git://src.adamsgaard.dk/pism
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
       tenthalpy-basal-bc.tex (2452B)
       ---
            1 \documentclass{standalone}
            2 \usepackage{amsmath}
            3 \usepackage{tikz}
            4 \usetikzlibrary{arrows,shapes}
            5 
            6 % Define block styles
            7 \tikzstyle{decision} = [ellipse, draw, text width=7em, text badly centered, inner sep=2pt]
            8 \tikzstyle{block} = [rectangle, draw, text width=5em, text badly centered, rounded corners, minimum height=4em]
            9 \tikzstyle{line} = [draw, -latex']
           10 
           11 \newcommand{\bn}{\mathbf{n}}
           12 \newcommand{\bq}{\mathbf{q}}
           13 
           14 % compile with pdflatex, then run
           15 % convert -density 300 enthalpy-basal-bc.pdf -quality 90 enthalpy-basal-bc.png
           16 
           17 \begin{document}
           18 \begin{tikzpicture}[node distance = 3cm, auto]
           19     % Place nodes
           20     \node (invisiblestart) {};
           21 
           22     \node [decision, below of=invisiblestart, text height=0.2cm] (coldvstemp) {$H<H_{\text s}(p)$ ?};
           23     \node [decision, left of=coldvstemp, xshift=-4em] (excludebad) {$\eta_{\text b}>0$ ?};
           24     \node [block, below of=excludebad, text width=6em] (fixbad) {$H := H_{\text s}(p)$};
           25 
           26     % edges
           27     \path [line] (invisiblestart) -- (coldvstemp);
           28     \path [line] (excludebad) -- node [text width=6em] {yes (consider base to be temperate)} (fixbad);
           29 
           30     % cold branch:
           31     \node [block, left of=fixbad, text width=7.5em] (coldmodeltype) {Eqn (49) is Neumann b.c.~for Eqn (67); $M_b=0$};
           32     % edges
           33     \path [line] (coldvstemp) -- node {yes} (excludebad);
           34     \path [line] (excludebad) -- node {no} (coldmodeltype);
           35 
           36     % temperate branch
           37     \node [block, below of=coldvstemp, text width=12em] (qtemperate) {$\nabla H \cdot \bn=0$ is Neumann b.c.~for Eqn (67)};
           38     \node [decision, below left of=qtemperate, text width=8em] (tempthick) {positive thickness of temperate ice at base?};
           39     \node [block, below right of=tempthick, text width=10em] (Mbforqtemperate) {$\bq = - k(H,p)\nabla T_{\text m}(p)$ \\ at ice base};
           40     \node [block, below left of=tempthick, text width=9em, xshift=-4em] (Mbforqcold) {$\bq = - K_{\text i}(H) \nabla H$ \\ at ice base};
           41     \node [block, below left of=Mbforqtemperate, text width=9em] (getMbtemp) {compute $M_b$ from Eqn (50) or (66)};
           42 
           43     % edges
           44     \path [line] (fixbad) -- (qtemperate);
           45     \path [line] (coldvstemp) -- node {no} (qtemperate);
           46     \path [line] (tempthick) -- node [above] {no} (Mbforqcold);
           47     \path [line] (tempthick) -- node {yes} (Mbforqtemperate);
           48     \path [line] (qtemperate) -- (tempthick);
           49     \path [line] (Mbforqcold) -- node {} (getMbtemp);
           50     \path [line] (Mbforqtemperate) -- node {} (getMbtemp);
           51 \end{tikzpicture}
           52