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