Changeset 1531
- Timestamp:
- Jan 26, 2015 1:58:29 PM (10 years ago)
- Location:
- palm/trunk/TUTORIAL/SOURCE
- Files:
-
- 7 edited
-
basic_equations.tex (modified) (7 diffs)
-
final_remarks.tex (modified) (1 diff)
-
non_cyclic_boundary_conditions.tex (modified) (1 diff)
-
parallelization.tex (modified) (3 diffs)
-
parallelization_figures/perf_4.png (modified) (previous)
-
program_structure.tex (modified) (1 diff)
-
user_defined_code.tex (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/TUTORIAL/SOURCE/basic_equations.tex
r1413 r1531 1 1 % $Id$ 2 2 \input{header_tmp.tex} 3 %\input{../header_ lectures.tex}3 %\input{../header_LECTURE.tex} 4 4 5 5 \usepackage[utf8]{inputenc} 6 6 \usepackage{ngerman} 7 7 \usepackage{pgf} 8 \usetheme{Dresden}9 8 \usepackage{subfigure} 10 9 \usepackage{units} … … 14 13 \usepackage{xmpmulti} 15 14 \usepackage{tikz} 16 \usetikzlibrary{shapes,arrows,positioning }15 \usetikzlibrary{shapes,arrows,positioning,decorations.pathreplacing} 17 16 \def\Tiny{\fontsize{4pt}{4pt}\selectfont} 17 18 %---------- neue Pakete 18 19 \usepackage{amsmath} 19 20 \usepackage{amssymb} 20 21 \usepackage{multicol} 21 22 23 \institute{Institut fÃŒr Meteorologie und Klimatologie, Leibniz UniversitÀt Hannover} 22 \usepackage{pdfcomment} 23 \usepackage{xcolor} 24 25 \institute{Institute of Meteorology and Climatology, Leibniz UniversitÀt Hannover} 26 \selectlanguage{english} 24 27 \date{last update: \today} 25 28 \event{PALM Seminar} 26 29 \setbeamertemplate{navigation symbols}{} 27 30 \setbeamersize{text margin left=.5cm,text margin right=.2cm} 28 31 \setbeamertemplate{footline} 29 { 32 {% 30 33 \begin{beamercolorbox}[rightskip=-0.1cm]& 31 34 {\includegraphics[height=0.65cm]{imuk_logo.pdf}\hfill \includegraphics[height=0.65cm]{luh_logo.pdf}} 32 35 \end{beamercolorbox} 33 \begin{beamercolorbox}[ht=2.5ex,dp=1.125ex, 34 leftskip=.3cm,rightskip=0.3cm plus1fil]{title in head/foot} 35 {\leavevmode{\usebeamerfont{author in head/foot}\insertshortauthor} \hfill \eventname \hfill \insertframenumber \; / \inserttotalframenumber} 36 \end{beamercolorbox} 37 \begin{beamercolorbox}[colsep=1.5pt]{lower separation line foot} 38 \end{beamercolorbox} 39 } 40 %\logo{\includegraphics[width=0.3\textwidth]{luhimuk_logo.pdf}} 36 \begin{beamercolorbox}[ht=2.5ex,dp=1.125ex,% 37 leftskip=.3cm,rightskip=0.3cm plus1fil]{title in head/foot}% 38 {\leavevmode{\usebeamerfont{author in head/foot}\insertshortauthor} \hfill \eventname \hfill \insertframenumber \; / \inserttotalframenumber}% 39 \end{beamercolorbox}% 40 % \begin{beamercolorbox}[colsep=1.5pt]{lower separation line foot}% 41 % \end{beamercolorbox} 42 }%\logo{\includegraphics[width=0.3\textwidth]{luhimuk_logo.eps}} 41 43 42 44 \title[Basic Equations]{Basic Equations} 43 \author{ Siegfried Raasch}45 \author{PALM group} 44 46 45 47 \begin{document} … … 159 161 \end{equation*} 160 162 \end{itemize} 161 \onslide<6-> \tikzstyle{plain} = [rectangle, draw, text width=0.2 7\textwidth, font=\small]163 \onslide<6-> \tikzstyle{plain} = [rectangle, draw, text width=0.255\textwidth, font=\small] 162 164 163 165 \begin{tikzpicture}[remember picture, overlay] … … 218 220 \end{column} 219 221 \begin{column}{0.32\textwidth} 220 \vspace{45mm} 222 \vspace{45mm} 223 \hspace{-1.75cm} 221 224 \begin{footnotesize} 222 225 \onslide<3->$R_{ki}$: \textbf{Reynolds-stress} \\ 223 $C_{ki}$: \textbf{cross-stress} \\ 224 $L_{ki}$: \textbf{Leonard-stress} \\ 225 $\tau_{ki}$: \textbf{total stress-tensor} 226 \hspace*{-1.5cm}$C_{ki}$: \textbf{cross-stress} \\ 227 \hspace*{-1.5cm}$L_{ki}$: \textbf{Leonard-stress} \\ 228 \hspace*{-1.5cm}$\tau_{ki}$: \textbf{total stress-tensor}\\ 229 \hspace*{-1.05cm} \textbf{generalized Reynolds stress}\\ 226 230 \end{footnotesize} 227 231 \end{column} 228 232 \end{columns} 229 \onslide<4->\tikzstyle{plain} = [rectangle, draw, text width=0.2 7\textwidth, font=\small]233 \onslide<4->\tikzstyle{plain} = [rectangle, draw, text width=0.25\textwidth, font=\small] 230 234 \begin{tikzpicture}[remember picture, overlay] 231 235 \node at (current page.north west){ … … 261 265 \end{equation*} 262 266 \end{itemize} 263 \onslide<1->\tikzstyle{plain} = [rectangle, draw, text width=0.2 7\textwidth, font=\small]267 \onslide<1->\tikzstyle{plain} = [rectangle, draw, text width=0.25\textwidth, font=\small] 264 268 \begin{tikzpicture}[remember picture, overlay] 265 269 \node at (current page.north west){ … … 311 315 312 316 \begin{picture}(0.0,0.0) 313 \put(140, 13){\uncover<6->{\includegraphics[width=0.6\textwidth]{basic_equations_figures/explicit_implicit.png}}}317 \put(140,0){\uncover<6->{\includegraphics[width=0.6\textwidth]{basic_equations_figures/explicit_implicit.png}}} 314 318 \end{picture} 315 319 \end{frame} … … 351 355 \onslide<2-> 352 356 \begin{flalign*} 353 &\frac{\partial \overline{u_i}}{\partial t} = - \frac{\partial \overline{u_k}\,\overline{u_i}}{\partial x_k} - \frac{1}{\rho_0} \frac{\partial \overline{\pi}^*}{\partial x_i} - \varepsilon_{ijk}f_j \overline{u_k} + \varepsilon_{i3k} f_3 \overline{u}_{k_\mathrm{g}} + g \frac{\overline{ T}-T_0}{T_0} \delta_{i3} + \nu \frac{\partial^2 \overline{u_i}}{\partial x_k^2} - \frac{\partial \tau_{ki}^r}{\partial x_k}&357 &\frac{\partial \overline{u_i}}{\partial t} = - \frac{\partial \overline{u_k}\,\overline{u_i}}{\partial x_k} - \frac{1}{\rho_0} \frac{\partial \overline{\pi}^*}{\partial x_i} - \varepsilon_{ijk}f_j \overline{u_k} + \varepsilon_{i3k} f_3 \overline{u}_{k_\mathrm{g}} + g \frac{\overline{\theta}-\theta_0}{\theta_0} \delta_{i3} + \nu \frac{\partial^2 \overline{u_i}}{\partial x_k^2} - \frac{\partial \tau_{ki}^r}{\partial x_k}& 354 358 \end{flalign*} 355 359 \item<4-> First principle (using potential\\ temperature): -
palm/trunk/TUTORIAL/SOURCE/final_remarks.tex
r1515 r1531 151 151 \vspace{5mm} 152 152 \item<3->{Ready for deep convection using anelastic system of equations} 153 \vspace{5mm} 154 \item<4->{Nesting methods are under development} 153 155 \end{itemize} 154 156 \end{frame} -
palm/trunk/TUTORIAL/SOURCE/non_cyclic_boundary_conditions.tex
r1515 r1531 385 385 \item<2->{So far, these methods have been applied only to a few special cases (cold air outbreaks, urban canopy layer for neutral stratification). Other setups may require modifications.} 386 386 \item<3->{Biggest problems are caused by gravity waves in capping inversions. Simulations with pure neutral stratification cause less problems.} 387 \item<4->{A synthetic turbulence generator will be available around mid 2015 as an additional option for creating inflow turbulence} 387 388 \end{itemize} 388 389 \end{frame} -
palm/trunk/TUTORIAL/SOURCE/parallelization.tex
r1515 r1531 249 249 \ \\ 250 250 \ \\ 251 \vspace*{-0.05cm} 251 252 \textcolor{blue}{columns of i}\\ 253 \vspace*{-0.04cm} 252 254 \textcolor{red}{no contiguous data at all}\\ 255 \vspace*{0.17cm} 253 256 \onslide<9-> \textcolor{blue}{columns of k}\\ 257 \vspace*{-0.04cm} 254 258 \textcolor{red}{planes of k,j (all data contiguous)} 255 259 \end{column} … … 262 266 \vspace{2mm} 263 267 \item<12-> OpenMP parallelization as well as mixed usage of OpenMP and 264 MPI is alsorealized.268 MPI is realized. 265 269 \end{itemize} 266 270 \end{frame} … … 570 574 \frametitle{Performance Examples (III)} 571 575 \begin{itemize} 572 \item Simulation with $ 2160^3$ grid points ($\sim$ 2TByte memory)576 \item Simulation with $4320^3$ grid points ($\sim$ 13 TByte memory) 573 577 \end{itemize} 574 578 \begin{columns}[T] 575 579 \begin{column}{0.5\textwidth} 576 \includegraphics[scale=0. 3]{parallelization_figures/perf_4.png} \\580 \includegraphics[scale=0.5]{parallelization_figures/perf_4.png} \\ 577 581 \scriptsize 578 \quad Cray-XC 30, HLRN-III, Hannover\\582 \quad Cray-XC40, HLRN-III, Hannover\\ 579 583 \quad (2D-domain decomposition) 580 584 \end{column} -
palm/trunk/TUTORIAL/SOURCE/program_structure.tex
r1515 r1531 916 916 \frametitle{Important Variables and Their Declaration} 917 917 \begin{itemize} 918 \item<2->{3D-arrays of prognostic variables are named $\Psi$, and $\Psi_p$ for time level $t$, and $t+\Delta t$, respectively, with $\Psi = u$, $v$, $w$, $pt$, $q$, $e$, $sa$ }918 \item<2->{3D-arrays of prognostic variables are named $\Psi$, and $\Psi_p$ for time level $t$, and $t+\Delta t$, respectively, with $\Psi = u$, $v$, $w$, $pt$, $q$, $e$, $sa$, $u\_p$, $v\_p$, ...} 919 919 \vspace{2.5mm} 920 920 \item<3->{They are by default declared as $\Psi$($z$,$y$,$x$) or $\Psi$($k$,$j$,$i$), e.g.\\ -
palm/trunk/TUTORIAL/SOURCE/user_defined_code.tex
r1515 r1531 616 616 \texttt{!-- No calls for single grid points are allowed at locations before and}\\ 617 617 \texttt{!-- after the timestep, since these calls are not within an i,j-loop}\\ 618 ~~~~~~ ~\texttt{SELECT CASE ( location )}\\618 ~~~~~~\texttt{SELECT CASE ( location )}\\ 619 619 \texttt{...}\\ 620 620 ~~~~~~~~~~\texttt{CASE ( 'after\_timestep' )}\\ 621 621 \texttt{!}\\ 622 \texttt{!-- Enter actions to be done after every timestep here}\\622 \texttt{!-- \hspace{6mm}Enter actions to be done after every timestep here}\\ 623 623 \par\medskip 624 624 ~~~~~~~~~~\texttt{CASE ( 'u-tendency' )}\\ 625 625 \texttt{!}\\ 626 \texttt{!-- Enter actions to be done in the u-tendency term here}\\627 ~~~~~~~~~~~~~ \texttt{\textcolor{blue}{DO i = nxl, nxr}}\\628 ~~~~~~~~~~~~~~~~ \texttt{\textcolor{blue}{DO j = nys, nyn}}\\629 ~~~~~~~~~~~~~~~~~~~ \texttt{\textcolor{blue}{DO k = nxb+1, nzt}}\\630 ~~~~~~~~~~~~~~~~~~~~~~~~~ \texttt{\textcolor{blue}{tend(k,j,i) = tend(k,j,i) - const * u(k,j,i) ...}}\\631 ~~~~~~~~~~~~~~~~~~~ \texttt{\textcolor{blue}{ENDDO}}\\632 ~~~~~~~~~~~~~~~~ \texttt{\textcolor{blue}{ENDDO}}\\633 ~~~~~~~~~~~~~ \texttt{\textcolor{blue}{ENDDO}}\\626 \texttt{!-- \hspace{6mm}Enter actions to be done in the u-tendency term here}\\ 627 ~~~~~~~~~~~~~~\texttt{\textcolor{blue}{DO i = nxl, nxr}}\\ 628 ~~~~~~~~~~~~~~~~~\texttt{\textcolor{blue}{DO j = nys, nyn}}\\ 629 ~~~~~~~~~~~~~~~~~~~~\texttt{\textcolor{blue}{DO k = nxb+1, nzt}}\\ 630 ~~~~~~~~~~~~~~~~~~~~~~~~~~\texttt{\textcolor{blue}{tend(k,j,i) = tend(k,j,i) - const * u(k,j,i) ...}}\\ 631 ~~~~~~~~~~~~~~~~~~~~\texttt{\textcolor{blue}{ENDDO}}\\ 632 ~~~~~~~~~~~~~~~~~\texttt{\textcolor{blue}{ENDDO}}\\ 633 ~~~~~~~~~~~~~~\texttt{\textcolor{blue}{ENDDO}}\\ 634 634 \par\bigskip 635 635 ~~~~~~~~~~\texttt{CASE ( 'v-tendency' )}\\ … … 671 671 \node [yellow,shift={(9.2 cm,3.0cm)}] at (current page.south west){~\texttt{SUBROUTINE user\_actions\_ij( i, j, location )}\\ 672 672 \par\medskip 673 ~~~ \texttt{USE control\_parameters}\\674 ~~~ \texttt{USE pegrid}\\675 ~~~ \texttt{USE user}\\673 ~~~~~~\texttt{USE control\_parameters}\\ 674 ~~~~~~\texttt{USE pegrid}\\ 675 ~~~~~~\texttt{USE user}\\ 676 676 \par\medskip 677 ~~~ \texttt{IMPLICIT NONE}\677 ~~~~~~\texttt{IMPLICIT NONE}\ 678 678 \par\medskip 679 ~~~ \texttt{CHARACTER (LEN=*) :: location}\\679 ~~~~~~\texttt{CHARACTER (LEN=*) :: location}\\ 680 680 \par\medskip 681 ~~~ \texttt{INTEGER(iwp) :: i, idum, j}\\681 ~~~~~~\texttt{INTEGER(iwp) :: i, idum, j}\\ 682 682 \par\medskip 683 683 \texttt{!}\\ 684 684 \texttt{!-- Here the user-defined actions follow}\\ 685 685 \vspace{0.5mm} 686 ~~~ \texttt{SELECT CASE ( location )}\\686 ~~~~~~\texttt{SELECT CASE ( location )}\\ 687 687 \par\medskip 688 ~~~~~~ \texttt{CASE ( 'u-tendency' )}\\688 ~~~~~~~~~\texttt{CASE ( 'u-tendency' )}\\ 689 689 \texttt{!}\\ 690 \texttt{!--Enter actions to be done in the u-tendency term here}\\691 ~~~~~~~~~ \texttt{\textcolor{blue}{DO k = nzb+1, nzt-1}}\\692 ~~~~~~~~~~~~~~~ \texttt{\textcolor{blue}{tend(k,j,i) = tend(k,j,i) + ...}}\\693 ~~~~~~~~~ \texttt{\textcolor{blue}{ENDDO}}\\690 ~~~~~~\texttt{!-- Enter actions to be done in the u-tendency term here}\\ 691 ~~~~~~~~~~~~\texttt{\textcolor{blue}{DO k = nzb+1, nzt-1}}\\ 692 ~~~~~~~~~~~~~~~~~~\texttt{\textcolor{blue}{tend(k,j,i) = tend(k,j,i) + ...}}\\ 693 ~~~~~~~~~~~~\texttt{\textcolor{blue}{ENDDO}}\\ 694 694 \par\medskip 695 ~~~~~~ \texttt{CASE ( 'v-tendency' )}};695 ~~~~~~~~~\texttt{CASE ( 'v-tendency' )}}; 696 696 % {% 697 697 % \begin{tikzpicture}[remember picture, overlay]
Note: See TracChangeset
for help on using the changeset viewer.
