source: palm/trunk/TUTORIAL/SOURCE/user_defined_code.tex @ 1703

Last change on this file since 1703 was 1649, checked in by knoop, 9 years ago

Latex lecture update for PALM seminar in Hong Kong

  • Property svn:keywords set to Id
File size: 46.4 KB
Line 
1% $Id: user_defined_code.tex 1649 2015-09-15 16:34:42Z raasch $
2\input{header_tmp.tex}
3%\input{../header_lectures.tex}
4
5\usepackage[utf8]{inputenc}
6\usepackage{ngerman}
7\usepackage{pgf}
8\usepackage{subfigure}
9\usepackage{units}
10\usepackage{multimedia}
11\usepackage{hyperref}
12\newcommand{\event}[1]{\newcommand{\eventname}{#1}}
13\usepackage{xmpmulti}
14\usepackage{tikz}
15\usetikzlibrary{shapes,arrows,positioning}
16%\usetikzlibrary{decorations.markings}             %neues paket
17%\usetikzlibrary{decorations.pathreplacing}        %neues paket
18\def\Tiny{\fontsize{4pt}{4pt}\selectfont}
19\usepackage{amsmath}
20\usepackage{amssymb}
21\usepackage{multicol}
22\usepackage{pdfcomment}
23\usepackage{graphicx}
24\usepackage{listings}
25\lstset{showspaces=false,language=fortran,basicstyle=
26        \ttfamily,showstringspaces=false,captionpos=b}
27
28\institute{Institute of Meteorology and Climatology, Leibniz UniversitÀt Hannover}
29\selectlanguage{english}
30\date{last update: \today}
31\event{PALM Seminar}
32\setbeamertemplate{navigation symbols}{}
33
34\setbeamertemplate{footline}
35  {
36    \begin{beamercolorbox}[rightskip=-0.1cm]&
37     {\includegraphics[height=0.65cm]{imuk_logo.pdf}\hfill \includegraphics[height=0.65cm]{luh_logo.pdf}}
38    \end{beamercolorbox}
39    \begin{beamercolorbox}[ht=2.5ex,dp=1.125ex,
40      leftskip=.3cm,rightskip=0.3cm plus1fil]{title in head/foot}
41      {\leavevmode{\usebeamerfont{author in head/foot}\insertshortauthor} \hfill \eventname \hfill \insertframenumber \; / \inserttotalframenumber}
42    \end{beamercolorbox}
43    \begin{beamercolorbox}[colsep=1.5pt]{lower separation line foot}
44    \end{beamercolorbox}
45  }
46%\logo{\includegraphics[width=0.3\textwidth]{luhimuk_logo.pdf}}
47
48\title[The PALM User-Interface]{The PALM User-Interface}
49\author{PALM group}
50
51\begin{document}
52
53% Folie 1
54\begin{frame}
55\titlepage
56\end{frame}
57
58\section{The PALM User-Interface}
59\subsection{The PALM User-Interface}
60
61% Folie 2
62\begin{frame}
63   \frametitle{Purpose of the User-Interface}
64   \scriptsize
65   \begin{itemize}
66      \item<2-> The standard (default) PALM code cannot account for every specific demand of a user. In order to include these specific demands, the user would have to modify the standard code.
67   \end{itemize}
68   \vspace{2mm}
69   \onslide<3->
70   \textbf{Problem:}
71   \begin{itemize}
72      \item<3-> New releases of PALM (current release is 3.10) would require the user to add his modifications to the new release again.
73   \end{itemize}   
74   \vspace{3mm}
75   \onslide<4->
76   \textbf{Solution:}
77   \begin{itemize}
78      \item<4-> PALM offers a \grqq user-interface\grqq, i.e. a set of subroutines, where the user can add his modifications, and which can be re-used for future releases of the standard PALM code.
79      \item<5-> By using the user-interface, the standard code does not have to be modified by the user in most of the cases.
80      \item<6-> The user-interface subroutines are almost \grqq empty\grqq\, by default. They are called from the standard PALM code but (with some very minor exceptions) do not contain any executable code.
81   \end{itemize}   
82\end{frame}
83
84% Folie 3
85\begin{frame}
86   \frametitle{General Structure of the User-Interface}
87   \tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\Tiny]
88   \scriptsize
89   \begin{itemize}
90      \item<1-> All routines can be found under \texttt{.../trunk/SOURCE}.
91      \item<2-> There is one file for each routine. Filenames are \texttt{user\underline{ }*.f90}.\\
92Example: \texttt{user\underline{ }last\underline{ }actions.f90}
93   \end{itemize}
94   \onslide<3->
95   \begin{center}
96   \begin{tikzpicture}[auto]
97      \node[yellow]{~\texttt{SUBROUTINE user\_last\_actions} \\
98                    \par\medskip
99                    \texttt{!------------------------------------------------------------------------------!}\\
100                    \texttt{!}\\
101                    \texttt{! Description:}\\
102                    \texttt{! ------------}\\
103                    \texttt{! Execution of user-defined actions at the end of a job.}\\
104                    \texttt{!------------------------------------------------------------------------------!}\\
105                    \par\medskip
106                    ~~~~\texttt{USE control\_parameters}\\
107                    ~~~~\texttt{USE kinds}\\
108                    ~~~~\texttt{USE user}\\
109                    \par\medskip
110                    ~~~~\texttt{IMPLICIT NONE}\\
111                    \par\medskip
112                    \texttt{!}\\
113                    \texttt{!-- Here the user-defined actions at the end of a job follow.}\\
114                    \texttt{!-- Sample for user-defined output:}\\
115                    \vspace{0.5mm}
116                    ~~~~\texttt{IF ( write\_binary(1:4) == 'true' )  THEN}\\
117                    \vspace{0.5mm}
118                    \texttt{!~~~~~~~IF ( ALLOCATED( u2\_av ) )  THEN}\\
119                    \vspace{0.5mm}
120                    \texttt{!~~~~~~~~~~WRITE ( 14 )  'u2\_av~~~~~~~~~~~~~~~';  WRITE ( 14 )  u2\_av}\\
121                    \vspace{0.5mm}
122                    \texttt{!~~~~~~~ENDIF}\\
123                    \par\medskip
124                    ~~~~~~~\texttt{WRITE ( 14 )  '*** end user ***~~~~'}\\
125                    \par\medskip
126                    ~~~~\texttt{ENDIF}\\
127                    \par\medskip
128                    ~\texttt{END SUBROUTINE user\_last\_actions}};
129   \end{tikzpicture}
130   \end{center}
131%   \includegraphics[scale=0.3]{user_defined_code_figures/user_last_actions.png}\\   
132\end{frame}
133
134% Folie 4
135\begin{frame}
136   \frametitle{Embedding of User-Interface Routines}
137   \tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\Tiny]
138   \scriptsize
139   \begin{itemize}
140      \item<1-> The user-interface routines are called from specific, well-defined locations in the standard PALM code.\\
141   \ \\
142   Example from \texttt{palm.f90}:
143   \end{itemize}
144   \onslide<2->
145   \begin{center}
146      \begin{tikzpicture}[auto]
147         \node[yellow]{~\texttt{...}\\
148                       \vspace{0.5mm}
149                       \texttt{!}\\
150                       \texttt{!-- If required, final user-defined actions, and}\\
151                       \texttt{!-- last actions on the open files and close files. Unit 14 was opened}\\
152                       \texttt{!-- in write\_3d\_binary but it is closed here, to allow writing on this}\\
153                       \texttt{!-- unit in routine user\_last\_actions.}\\
154                       \vspace{0.5mm}
155                       ~~~~\texttt{CALL cpu\_log( log\_point(4), 'last actions', 'start' )}\\
156                       ~~~~\texttt{DO  i = 0, io\_blocks-1}\\
157                       ~~~~~~~\texttt{IF ( i == io\_group )  THEN}\\
158                       ~~~~~~~~~~\texttt{\textcolor{blue}{CALL user\_last\_actions}}\\
159                       \vspace{0.5mm}
160                       ~~~~~~~~~~\texttt{IF ( write\_binary(1:4) == 'true' )  CALL close\_file( 14 )}\\
161                       ~~~~~~~\texttt{ENDIF}\\
162                       \vspace{0.5mm}
163                       \texttt{\#if defined( \_\_parallel )}\\
164                       ~~~~~~~\texttt{CALL MPI\_BARRIER( comm2d, ierr )}\\
165                       \vspace{0.5mm}
166                       \texttt{\#endif}\\
167                       \vspace{0.5mm}
168                       ~~~~\texttt{ENDDO}\\
169                       \vspace{0.5mm}
170                       ~~~~\texttt{CALL close\_file( 0 )}\\
171                       \vspace{0.5mm}
172                       ~~~~\texttt{CALL close\_dvrp}\\
173                       ~~~~\texttt{CALL cpu\_log( log\_point(4), 'last actions', 'stop' )}\\
174                       \vspace{1mm}
175                       ~\texttt{...}\\
176                       \vspace{1mm}
177                       \texttt{!}\\
178                       \texttt{!-- Take final CPU-time for CPU-time analysis}\\
179                       ~~~~\texttt{CALL cpu\_log( log\_point(1), 'total', 'stop' )}\\
180                       ~~~~\texttt{CALL cpu\_statistics}\\
181                       \texttt{\#if defined( \_\_parallel )}\\
182                       ~~~~\texttt{CALL MPI\_FINALIZE( ierr )}\\
183                       \vspace{0.5mm}
184                       \texttt{\#endif}\\
185                       \par\medskip
186                       ~\texttt{END PROGRAM palm}};
187      \end{tikzpicture} 
188   \end{center}
189
190%   \includegraphics[scale=0.33]{user_defined_code_figures/embedding.png}\\
191\end{frame}
192
193
194% Folie 5
195\begin{frame}
196   \frametitle{Embedding of User-Interface Routines\\
197   Flow Chart Overview (I): Initial Steps}
198   \scriptsize
199   \begin{tikzpicture}[scale=0.8, transform shape]
200      % Define colours
201      \tikzstyle{green} = [draw, fill=green!30, text width=7.1em, text centered, minimum height=1.6em]
202      \tikzstyle{green1} = [draw, fill=green!30, text width=12em, text centered, minimum height=1.6em]
203      \tikzstyle{green2} = [draw, fill=green!30, text width=10em, text centered, minimum height=1.6em]
204      \tikzstyle{yellow} = [draw, fill=yellow!30, text width=7.1em, text centered, minimum height=1.6em]
205      \tikzstyle{orange} = [draw, fill=orange!60, text width=7.1em, text centered]
206     
207      % Define coordinates and draw nodes     
208      \node[yellow, text width=5em] (PALM) at (0,0.1) {PALM};
209
210      \coordinate (A) at (0,-0.7);
211      \node[yellow, right=0.6cm of A] (parin) {parin};
212
213      \node[green, right=1mm of parin] (userparin) {user\_parin};
214      \draw[-] (parin) -- (userparin);
215
216      \coordinate[below=0.6cm of A] (B);
217      \node[yellow, right=0.6cm of B] (initgrid) {init\_grid};
218
219      \node[green, right=1mm of initgrid] (userinigrid) {user\_init\_grid};
220      \draw[-] (initgrid) -- (userinigrid);
221
222      \coordinate[below=0.6cm of B] (C);
223      \node[yellow, right=0.6cm of C] (checkparam) {check\_parameters};
224
225      \coordinate[below=0.3cm of checkparam] (D);
226      \node[green1, right=0.6cm of D] (usercheckdataout) {user\_check\_data\_output/\_pr};
227
228      \coordinate[below=0.55cm of D] (E);
229      \node[green1, right=0.6cm of E] (usercheckparam) {user\_check\_parameters};
230
231      \coordinate[below=1.7cm of C] (F);
232      \node[yellow, right=0.6cm of F] (init3dmodel) {init\_3d\_model};
233
234      \coordinate[below=0.3cm of init3dmodel] (G);
235      \node[green2, right=0.6cm of G] (userinit3d) {user\_init\_3d\_model};
236
237      \coordinate[below=0.55cm of G] (H);
238      \node[yellow, right=0.6cm of H] (read3dbin) {read\_3d\_binary};
239
240      \coordinate[below=0.3cm of read3dbin] (I);
241      \node[green1, right=0.6cm of I] (userreadrestart) {user\_read\_restart\_data};
242
243      \coordinate[below=1.15cm of H] (J);
244      \node[green1, right=0.6cm of J] (userinitplan) {user\_init\_plant\_canopy};
245
246      \coordinate[below=0.6cm of J] (K);
247      \node[yellow, right=0.6cm of K] (lpminit) {lpm\_init};
248
249      \node[green, right=1mm of lpminit] (userlpminit) {user\_lpm\_init};
250      \draw[-] (lpminit) -- (userlpminit);     
251
252      \coordinate[below=0.6cm of K] (L);
253      \node[green, right=0.6cm of L] (userinit) {user\_init};
254
255      \coordinate[below=4.05cm of F] (M);
256      \node[yellow, right=0.6cm of M] (header) {header};
257
258      \node[green, right=1mm of header] (userheader) {user\_header};
259      \draw[-] (header) -- (userheader);     
260%     
261       % Draw lines
262       \draw[-, thick] (PALM) -- (0,-7.8);
263       \draw[dotted, thick] (0,-8.2) -- (0,-7.8);
264       \draw[-] (A) -- (parin);
265       \draw[-] (B) -- (initgrid);
266       \draw[-] (C) -- (checkparam);
267       \draw[-] (checkparam.south) |- (usercheckdataout.west);
268       \draw[-] (checkparam.south) |- (usercheckparam.west);
269       \draw[-] (F) -- (init3dmodel);
270       \draw[-,dashed] (init3dmodel.south) |- (userinit3d.west);
271       \draw[-,dashed] (init3dmodel.south) |- (read3dbin.west);
272       \draw[-] (read3dbin.south) |- (userreadrestart.west);
273       \draw[-,dashed] (init3dmodel.south) |- (userinitplan);
274       \draw[-,dashed] (init3dmodel.south) |- (lpminit);
275       \draw[-] (init3dmodel.south) |- (userinit);
276       \draw[-] (M) -- (header);   
277     
278   \end{tikzpicture}
279\end{frame}
280
281% Folie 6
282\begin{frame}
283   \frametitle{Embedding of User-Interface Routines\\
284   Flow Chart Overview (II): Time Integration Loop}
285   \scriptsize
286   \begin{tikzpicture}[scale=0.7, transform shape]
287      % Define colours
288      \tikzstyle{green} = [draw, fill=green!30, text width=7.1em, text centered, minimum height=1.6em]
289      \tikzstyle{green1} = [draw, fill=green!30, text width=10em, text centered, minimum height=1.6em]
290      \tikzstyle{yellow} = [draw, fill=yellow!30, text width=7.1em, text centered, minimum height=1.6em]
291      \tikzstyle{orange} = [draw, fill=orange!60, text width=7.1em, text centered]
292      \tikzstyle{orange1} = [draw, fill=orange!80, text width=11em, text centered]     
293      \definecolor{grey}{rgb}{0.25,0.25,0.25}
294
295      % Define coordinates and draw nodes     
296      \node[draw, text width=25em,minimum height=4.0em,fill=black!15] (progframe) at (5.75,-1.5) {};
297
298      \coordinate (A) at (0,0);
299
300      \node[orange, right=0.6cm of A] (timeint) {time\_integration};       
301     
302      \coordinate[below=0.3cm of timeint] (B);
303      \coordinate[below=0.8cm of timeint] (C);   
304      \coordinate[right=2cm of C] (D);   
305      \coordinate[right=3.8cm of D] (E);
306                 
307      \node[green, right=0.6cm of B, text width=14em] (useractions1) {\textcolor{blue}{user\_actions\_(\textit{'before\_timestep'})}};   
308         
309      \node[orange1, below=0.2cm of D] (progeq) {prognostic\_equations\_vector};   
310      \node[orange1, below=0.2cm of E] (progeqfast) {prognostic\_equations\_cache};
311           
312      \node[green, right=0.3cm of progframe, text width=14em] (useractions2) {\textcolor{blue}{user\_actions\_(\textit{'*\_tendency'})}};
313     
314     
315      \coordinate[below=2.3cm of timeint] (F);
316      \coordinate[below=2.9cm of timeint] (G);
317      \coordinate[below=3.5cm of timeint] (H);
318      \coordinate[below=4.1cm of timeint] (I);
319
320                             
321      \node[yellow, right=0.6cm of F] (lpm) {lpm};
322      \node[green, right=0.9cm of lpm, text width=10em] (userlpmadvec) {user\_lpm\_advec};
323      \node[yellow, right=0.6cm of G] (exchangehoriz) {exchange\_horiz}
324%      \node[yellow, right=0.0cm of exchangehoriz] (asselin) {asselin\_filter};     
325      \node[yellow, right=0.6cm of H] (swaptimelevel) {swap\_timelevel};
326      \node[yellow, right=0.6cm of I] (pres) {pres};
327     
328      \coordinate[below=0.3cm of pres] (J);   
329      \coordinate[below=5.1cm of timeint] (K);     
330      \coordinate[below=5.7cm of timeint] (L);
331      \coordinate[below=6.3cm of timeint] (M);
332      \coordinate[below=6.9cm of timeint] (N);     
333      \coordinate[below=8.2cm of timeint] (O);
334      \coordinate[below=8.8cm of timeint] (P);     
335                 
336%      \node[yellow, right=0.6cm of K] (prandtlfluxes) {prandtl\_fluxes};   
337      \node[green, right=0.6cm of K, text width=14em] (useractions3) {\textcolor{blue}{user\_actions\_(\textit{'after\_integration'})}};             
338%      \node[yellow, right=0.6cm of L] (diffusivities) {diffusivities};             
339%      \node[yellow, right=0.6cm of L] (checkrestart) {check\_for\_restart};     
340      \node[yellow, right=0.6cm of L] (flowstat) {flow\_statistics};
341      \node[green, right=0.6cm of flowstat] (userstat2) {user\_statistics}
342
343      \node[yellow, right=0.6cm of M] (sum3ddata) {sum\_up\_3d\_data}
344      \node[green1, right=0.6cm of sum3ddata] (user3ddataav) {user\_3d\_data\_averaging};
345
346      \node[green, right=0.6cm of O, text width=14em] (useractions4) {\textcolor{blue}{user\_actions\_(\textit{'after\_timestep'})}};
347%      \node[yellow, right=0.6cm of P] (timestep) {timestep};             
348
349      \node[yellow, right=0.2cm of user3ddataav] (dataoutput) {data\_output\_***};
350      \node[green, right=0.2cm of dataoutput] (userdvrpcoltab) {user\_dvrp\_coltab};
351      \coordinate[below=0.2cm of dataoutput] (Q);   
352      \coordinate[right=0.4cm of Q] (R);   
353      \coordinate[right=3.0cm of R] (S);
354         
355      \node[green, below=0.2cm of R, text width=9em] (output2d) {user\_data\_output\_2d};
356      \node[green, below=0.2cm of S, text width=9em] (output3d) {user\_data\_output\_3d};
357     
358
359
360      \draw[dotted, thick] (0,0.05) -- (0,0.3);
361      \draw[-, thick] (0,0.05) -- (0,-8.8);
362      \draw[dotted, thick] (0,-9.2) -- (0,-8.8);
363      \draw[-, thick] (A) -- (timeint);
364      \draw[-, thick] (timeint) -- (P);
365      \draw[-, thick] (B) -- (useractions1);
366      \draw[-, thick] (C) -- (D);
367      \draw[-, thick] (D) -- (E);
368      \draw[->, thick] (E) -- (progeqfast);
369      \draw[-, thick] (progeq) -- (D);     
370      \draw[-, thick] (progframe) -- (useractions2);
371      \draw[-, thick] (F) -- (lpm);   
372      \draw[-, thick] (userlpmadvec) -- (lpm);   
373      \draw[-, thick] (G) -- (exchangehoriz);
374      \draw[-, thick] (H) -- (swaptimelevel);
375      \draw[-, thick] (I) -- (pres);     
376      \draw[dotted, thick] (pres) -- (J);
377%      \draw[-, thick] (K) -- (prandtlfluxes);   
378      \draw[-, thick] (K) -- (useractions3);     
379      \draw[-, thick] (L) -- (flowstat); 
380      \draw[-, thick] (M) -- (sum3ddata); 
381      \draw[-, thick] (sum3ddata) -- (user3ddataav);
382
383      \draw[-, thick] (userstat2) -- (flowstat); 
384      \draw[-, thick] (O) -- (useractions4);
385%      \draw[-, thick] (O) -- (timestep);   
386      \draw[-, thick] (dataoutput) -- (Q);
387      \draw[-, thick] (S) -- (Q);
388      \draw[-, thick] (R) -- (output2d);
389      \draw[-, thick] (S) -- (output3d);
390   \end{tikzpicture}
391\end{frame}
392
393
394% Folie 7
395\begin{frame}
396   \frametitle{Embedding of User-Interface Routines\\
397   Flow Chart Overview (III): Final Steps}
398   \scriptsize
399   \begin{tikzpicture}[scale=0.8, transform shape]
400      % Define colours
401      \tikzstyle{green} = [draw, fill=green!30, text width=7.1em, text centered, minimum height=1.6em]
402      \tikzstyle{yellow} = [draw, fill=yellow!30, text width=7.1em, text centered, minimum height=1.6em]
403      \tikzstyle{orange} = [draw, fill=orange!60, text width=7.1em, text centered]
404     
405      % Define coordinates and draw nodes     
406      \node[green, text width=5em] (PALM) at (0,0.1) {PALM};
407
408      \coordinate (A) at (0,-0.7);
409
410
411      \node[orange, right=0.6cm of A] (timeint) {time\_integration};       
412      \coordinate[below=0.5cm of timeint] (B);
413      \coordinate[below=1.0cm of timeint] (C);   
414     
415      \coordinate[below=2.5cm of PALM] (D); 
416     
417      \node[yellow, right=0.6cm of D] (write3dbinary) {write\_3d\_binary};
418      \node[yellow, right=0.9cm of write3dbinary] (writevarlist) {write\_var\_list};   
419       
420      \coordinate[below=3.5cm of PALM] (E);
421      \coordinate[below=4.1cm of PALM] (F);
422      \coordinate[below=4.7cm of PALM] (G);
423      \node[yellow, right=0.6cm of E] (header) {header};     
424      \node[green, right=0.6cm of F] (userlastactions) {user\_last\_actions};   
425      \node[yellow, right=0.6cm of G] (cpustat) {cpu\_statistics};               
426      % Draw lines
427      \draw[-, thick] (PALM) -- (G);
428      \draw[-, thick] (A) -- (timeint);
429      \draw[-, thick] (B) -- (timeint);
430      \draw[dotted, thick] (B) -- (C);
431      \draw[dashed, thick] (D) -- (write3dbinary);
432      \draw[-, thick] (write3dbinary) -- (writevarlist);
433      \draw[-, thick] (E) -- (header);
434      \draw[-, thick] (F) -- (userlastactions);
435      \draw[-, thick] (G) -- (cpustat);
436
437   \end{tikzpicture}
438\end{frame}
439
440% Folie 8
441\begin{frame}
442   \frametitle{Complete List of User-Interface Routines (I)}
443   \tiny
444   \begin{tabular}{|p{2.5cm}|p{2.4cm}|p{2.1cm}|p{3cm}|}
445      \hline
446           \textbf{Name} & \textbf{Arguments} & \textbf{Called from} & \textbf{Task}\\
447           \hline
448           \texttt{user\_3d\_data\_averaging} & \texttt{mode}, \texttt{variable} & \texttt{average\_3d\_data} + \texttt{sum\_up\_3d\_data} & temporal averaging for user-defined quantities\\
449           \hline
450           \texttt{user\_actions} \texttt{user\_actions} & \texttt{location} \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\, \texttt{i}, \texttt{j}, \texttt{location} & \texttt{time\_integration} +
451\texttt{prognostic\_equations} & e.g. additional forces to be included in the prognostic equations\\
452           \hline
453           \texttt{user\_dummy (user\_additional \_routines.f90)} & - - - & - - - & for additional subroutines defined by the user\\
454           \hline
455           \texttt{user\_check\_data\_output} & \texttt{variable}, \texttt{unit} & \texttt{check\_parameters} + \texttt{init\_masks}& check the user-defined output quantities\\
456           \hline
457           \texttt{user\_check\_data\_output\_pr} & \texttt{variable}, \texttt{var\_count}, \texttt{unit} & \texttt{check\_parameters} & check the user-defined profile output quantities\\
458           \hline
459           \texttt{user\_check\_parameters} & - - - & \texttt{check\_parameters} & check user-defined variables\\
460           \hline         
461           \texttt{user\_data\_output\_2d} & \texttt{av}, \texttt{variable}, \texttt{found}, \texttt{grid}, \texttt{local\_pf}, \texttt{two\_d} & \texttt{data\_output\_2d} & output/calculation of additional user-defined quantities\\       
462           \hline         
463           \texttt{user\_data\_output\_3d} & \texttt{av}, \texttt{variable}, \texttt{found}, \texttt{local\_pf}, \texttt{nz\_do} & \texttt{data\_output\_3d} & output/calculation of additional user-defined quantities\\
464           \hline         
465           \texttt{user\_data\_output\_dvrp} & \texttt{output\_variable}, \texttt{local\_pf} & \texttt{data\_output\_dvrp} & output of additional user-defined quantities\\
466           \hline
467           \texttt{user\_data\_output\_mask} & \texttt{av}, \texttt{variable}, \texttt{found}, \texttt{local\_pf} & \texttt{data\_output\_mask} & output of additional masked user-defined quantities\\
468           \hline
469           \texttt{user\_define\_netcdf\_grid} & \texttt{variable}, \texttt{found}, \texttt{grid\_x}, \texttt{grid\_y}, \texttt{grid\_z} & \texttt{netcdf} & defining the grid for additional output quantities\\
470           \hline
471           \texttt{user\_dvrp\_coltab} & \texttt{mode}, \texttt{variable} & \texttt{data\_output\_dvrp} & defining color tables for particles\\   
472           \hline
473           \texttt{user\_header} & \texttt{io} & \texttt{header} & output user variables to header\\
474           \hline
475           \texttt{user\_init} & - - - & \texttt{init\_3d\_model} & e.g. reading from restart file\\
476           \hline
477   \end{tabular}
478
479\end{frame}
480
481% Folie 9
482\begin{frame}
483   \frametitle{Complete List of User-Interface Routines (II)}
484   \tiny
485   \ \\
486   \begin{tabular}{|p{2.5cm}|p{2.4cm}|p{2.1cm}|p{3cm}|}
487      \hline
488           \textbf{Name} & \textbf{Arguments} & \textbf{Called from} & \textbf{Task}\\
489           \hline
490           \texttt{user\_init\_3d\_model} & - - - & \texttt{init\_3d\_model} & special initializations\\
491           \hline
492           \texttt{user\_init\_grid} & \texttt{gls} & \texttt{init\_grid} & defining a special topography\\
493           \hline         
494           \texttt{user\_init\_plant\_canopy} & - - - & \texttt{init\_3d\_model} & setting of leaf area density and canopy drag coefficient\\
495           \hline
496           \texttt{user\_last\_actions} & - - - & \texttt{palm} & e.g. output for restart runs\\
497           \hline
498           \texttt{user\_lpm\_advec} & - - - & \texttt{lpm} & modification of particles after advection \\
499           \hline
500           \texttt{user\_lpm\_init} & - - - & \texttt{lpm\_init} & defining initial particle sources \\
501           \hline
502           \texttt{user\_lpm\_set\_attributes} & - - - & \texttt{lpm} & defining particles attributes \\
503           \hline
504           \texttt{MODULE user (user\_module.f90)} & - - - & - - - & contains user defined variables\\
505           \hline                 
506           \texttt{user\_parin} &  & \texttt{parin} & reading user variables\\     
507           \hline 
508           \texttt{user\_read\_restart\_data} & \texttt{i}, \texttt{nxlfa}, \texttt{nxl\_on\_file}, \texttt{nxrfa}, \texttt{nxr\_on\_file}, \texttt{nynfa}, \texttt{nyn\_on\_file}, \texttt{nysfa}, \texttt{nys\_on\_file}, \texttt{offset\_xa}, \texttt{offset\_ya}, \texttt{overlap\_count}, \texttt{tmp\_2d}, \texttt{tmp\_3d} & \texttt{read\_3d\_binary} & reading user-defined 2d/3d-arrays from the restart file\\         
509           \hline                 
510           \texttt{user\_spectra} & \texttt{mode}, \texttt{m}, \texttt{pr} & \texttt{calc\_spectra} + \texttt{data\_output\_spectra} & output/calculation of additional user-defined quantities\\         
511           \hline         
512           \texttt{user\_statistics} & \texttt{mode}, \texttt{sr}, \texttt{tn} & \texttt{flow\_statistics} & calculating additional horizontal averages + time series quantities\\             
513           \hline             
514   \end{tabular}
515   \ \\
516   \ \\
517   See PALM online documentation under \\
518   \textbf{http://palm.muk.uni-hannover.de/trac/wiki/doc/app/userint}\\
519   for detailed explanations.
520
521\end{frame}
522
523% Folie 10
524\begin{frame}
525   \frametitle{Data Access / Exchange}
526   \tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.6\textwidth, font=\tiny]
527   \footnotesize
528   \begin{itemize}
529      \item<2-> \textbf{Between the standard PALM code and the user-interface:}\\
530      \begin{itemize}
531         \scriptsize
532         \item<3-> by including the respective \textcolor{blue}{PALM modules} in the user-interface subroutines.\\
533         \begin{center}
534%           \includegraphics[scale=0.3]{user_defined_code_figures/data_access_exchange.png}
535            \begin{tikzpicture}[auto]
536               \node[yellow]{\texttt{SUBROUTINE user\_actions( location )} \\
537                             \par\medskip
538                             ~~~~~\texttt{\textcolor{blue}{USE arrays\_3d}} \\
539                             ~~~~~\texttt{\textcolor{blue}{USE control\_parameters}} \\
540                             ~~~~~\texttt{\textcolor{blue}{USE cpulog}} \\
541                             ~~~~~\texttt{\textcolor{blue}{USE indices}} \\
542                             ~~~~~\texttt{\textcolor{blue}{USE interfaces}} \\
543                             ~~~~~\texttt{\textcolor{blue}{USE pegrid}} \\
544                             ~~~~~\texttt{\textcolor{red}{USE user}} \\
545                             \par\medskip
546                             ~~~~~\texttt{IMPLICIT NONE} \\
547                             \par\medskip
548                             ~~~~~\texttt{CHARACTER (LEN=$^\ast$) ::  location} \\
549                             \par\medskip
550                             ~~~~~\texttt{INTEGER ::  i, j, k}};
551            \end{tikzpicture}
552         \end{center} 
553      \end{itemize}
554      \item<4-> \textbf{Within the user-interface:}
555      \begin{itemize}
556         \scriptsize
557         \item by the module \texttt{\textcolor{red}{user}} (file \texttt{user\underline{ }module.f90}),  which is used in every subroutine included in the interface.\\ \textbf{This module is (and should be) never used in the standard PALM code (otherwise, the default code would depend on the user\\ interface).}
558      \end{itemize}
559   \end{itemize}
560
561\end{frame}
562
563% Folie 11
564\begin{frame}
565   \frametitle{Usage of \texttt{user\underline{ }actions} (I)}
566   \tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\tiny]
567   \scriptsize
568   \begin{itemize}
569      \item<1-> \texttt{user\underline{ }actions} is designed to add additional terms to the prognostic equations or to carry out special actions at the beginning or the end of each timestep.
570      \vspace{1.8mm}
571      \item<2-> Therefore, several calls of \texttt{user\underline{ }actions} can be found in the default PALM routines \texttt{time\underline{ }integration} and \texttt{prognostic\underline{ }equations}. The place from which it is called is communicated to the routine by a string-argument, e.g.\\
572   \ \\
573   \quad \texttt{CALL user\underline{ }actions( 'u-tendency' )}\\
574   \ \\
575        \onslide<3-> It means that this call is from a line within \texttt{prognostic\underline{ }equations}, where the tendencies for the u-component are calculated and integrated:
576   \end{itemize}
577   \vspace{-3mm}
578   \begin{center}
579   \begin{tikzpicture}[auto]
580      \node[yellow] {\texttt{DO  i = nxl, nxr} \\
581                     ~~~\texttt{DO  j = nys, nyn} \\
582                     \texttt{...} \\
583                     ~~~~~~\texttt{CALL diffusion\_u( i, j )}\\
584                     ~~~~~~\texttt{CALL coriolis( i, j, 1 )} \\
585                     \texttt{...}\\
586                     ~~~~~~\texttt{\textcolor{blue}{CALL user\_actions( i, j, 'u-tendency' )}}\\
587                     \texttt{!} \\
588                     \texttt{!-- Prognostic equation for u-velocity component} \\
589                     ~~~~~~\texttt{DO  k = nzb\_u\_inner(j,i)+1, nzt} \\
590                     ~~~~~~~~~\texttt{u\_p(k,j,i) = u(k,j,i) + dt\_3d * ( tsc(2) * tend(k,j,i) + \&} \\
591                     \hspace{4.08cm} \texttt{tsc(3) * tu\_m(k,j,i) )       \&} \\
592                     \hspace{2.93cm} \texttt{- tsc(5) * rdf(k) * ( u(k,j,i) - ug(k) )}\\
593                     ~~~~~~\texttt{ENDDO} \\
594                     \texttt{...}};
595   \end{tikzpicture}
596   \end{center}
597%    \begin{center}
598%       \includegraphics[scale=0.33]{user_defined_code_figures/user_actions_1.png}
599%    \end{center}
600\end{frame}
601
602% Folie 12
603\begin{frame}
604   \frametitle{Usage of \texttt{user\underline{ }actions} (II)}
605   \tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\tiny]
606   \scriptsize
607   \begin{itemize}
608      \item Additional tendencies have to be included by the user at the respective code line in \texttt{user\underline{ }actions}:
609   \end{itemize}
610   \begin{center}
611      \begin{tikzpicture}[auto]
612         \node[yellow]{~\texttt{SUBROUTINE user\_actions( location )}\\
613                       \texttt{...}\\
614                       \texttt{!}\\
615                       \texttt{!--    Here the user-defined actions follow}\\
616                       \texttt{!--    No calls for single grid points are allowed at locations before and}\\
617                       \texttt{!--    after the timestep, since these calls are not within an i,j-loop}\\
618                       ~~~~~~\texttt{SELECT CASE ( location )}\\
619                       \texttt{...}\\
620                       ~~~~~~~~~~\texttt{CASE ( 'after\_timestep' )}\\
621                       \texttt{!}\\
622                       \texttt{!-- \hspace{6mm}Enter actions to be done after every timestep here}\\
623                       \par\medskip
624                       ~~~~~~~~~~\texttt{CASE ( 'u-tendency' )}\\
625                       \texttt{!}\\
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                       \par\bigskip
635                       ~~~~~~~~~~\texttt{CASE ( 'v-tendency' )}\\
636                       \texttt{...}};
637      \end{tikzpicture}
638%      \includegraphics[scale=0.35]{user_defined_code_figures/user_actions_2.png}
639   \end{center}
640\end{frame}
641
642% Folie 13
643\begin{frame}[t]
644   \frametitle{Usage of \texttt{user\underline{ }actions} (III)}
645   \tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.5\textwidth, font=\Tiny]
646   \scriptsize
647   \begin{columns}[T]
648      \begin{column}{1.15\textwidth}
649         \begin{itemize}
650            \item<1-> The different versions of \texttt{prognostic\underline{ }equations} (\texttt{prognostic\underline{ }equations\underline{ }cache}, \texttt{prognostic\underline{ }equations\underline{ }vector}) contain different calls of \texttt{user\underline{ }actions}:
651            \item<2-> From \texttt{prognostic\underline{ }equations\underline{ }vector}: \hspace{0.6em}\texttt{CALL user\underline{ }actions( 'u-tendency' )}
652            \item<3-> From \texttt{prognostic\underline{ }equations},\\
653            \texttt{prognostic\underline{ }equations\underline{ }cache}:\hspace{4em}\texttt{CALL user\underline{ }actions( i, j, 'u-tendency' )}
654            \vspace{3mm}
655            \item<4-> In case that \texttt{prognostic\underline{ }equations}\\
656            \texttt{prognostic\underline{ }equations\underline{ }cache} are\\
657            used, the user has to add his code in\\
658            the interface routine\\
659            \texttt{user\underline{ }actions\underline{ }ij}:
660            \vspace{3mm}
661            \item<5-> Here, only the \texttt{k}-loop (vertical direction)\\
662            has to be used, because loops over\\
663            i and j are executed in\\
664            \texttt{prognostic\underline{ }equations\underline{ }cache}.
665         \end{itemize}
666      \end{column}
667   \end{columns}
668 
669   \uncover<4->{%
670   \begin{tikzpicture}[remember picture, overlay]
671      \node [yellow,shift={(9.2 cm,3.0cm)}]  at (current page.south west){~\texttt{SUBROUTINE user\_actions\_ij( i, j, location )}\\
672                                                                          \par\medskip
673                                                                          ~~~~~~\texttt{USE control\_parameters}\\
674                                                                          ~~~~~~\texttt{USE pegrid}\\
675                                                                          ~~~~~~\texttt{USE user}\\
676                                                                          \par\medskip
677                                                                          ~~~~~~\texttt{IMPLICIT NONE}\
678                                                                          \par\medskip
679                                                                          ~~~~~~\texttt{CHARACTER (LEN=*) ::  location}\\
680                                                                          \par\medskip
681                                                                          ~~~~~~\texttt{INTEGER(iwp) ::  i, idum, j}\\
682                                                                          \par\medskip
683                                                                          \texttt{!}\\
684                                                                          \texttt{!--    Here the user-defined actions follow}\\
685                                                                          \vspace{0.5mm}
686                                                                          ~~~~~~\texttt{SELECT CASE ( location )}\\
687                                                                          \par\medskip
688                                                                          ~~~~~~~~~\texttt{CASE ( 'u-tendency' )}\\
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}}\\
694                                                                          \par\medskip
695                                                                          ~~~~~~~~~\texttt{CASE ( 'v-tendency' )}};
696%          {%
697%          \begin{tikzpicture}[remember picture, overlay]
698% %            \node[text width=14em] at (0,0) {\includegraphics[scale=0.3]{user_defined_code_figures/user_actions_3.png}};
699%          \end{tikzpicture}
700%          };
701    \end{tikzpicture}}
702   
703\end{frame}
704
705% Folie 14
706\begin{frame}
707   \frametitle{Steering the User-Interface}
708   \tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\Tiny]
709   \scriptsize
710   \begin{columns}[T]
711      \begin{column}{1.1\textwidth}
712         For steering the user-interface code, the user may want to add some additional variables and set their respective values within the parameter-file (e.g. \texttt{example\underline{ }cbl\underline{ }p3d}). This requires the following actions (example for a variable named \texttt{foo}):
713      \end{column}
714   \end{columns}
715   \vspace{1mm}
716      \begin{enumerate}
717         \item<1-> Add the variable name to module \texttt{user} in order to define it and to make it available in all user-interface subroutines. Set a default value for this variable.
718         \begin{center}
719%            \includegraphics[scale=0.3]{user_defined_code_figures/steering_1.png}
720            \begin{tikzpicture}
721               \node[yellow]{~\texttt{MODULE user}\\
722                             \vspace{0.4mm}
723                             \texttt{...}\\
724                             \vspace{0.4mm}
725                             ~\texttt{REAL(wp) \hspace{0.4cm}:: ~~\textcolor{blue}{foo = 0.0}}\\
726                             \vspace{0.4mm}
727                             \texttt{...}\\
728                             \vspace{0.4mm}
729                             ~\texttt{END MODULE user}};
730            \end{tikzpicture}
731         \end{center}
732         \vspace{1mm}
733         \item<2-> Add the variable to the NAMELIST \texttt{/userpar/}. This NAMELIST already contains four predefined variables.
734         \begin{center}
735%            \includegraphics[scale=0.3]{user_defined_code_figures/steering_2.png}
736            \begin{tikzpicture}
737               \node[yellow]{~\texttt{SUBROUTINE user\_parin}\\
738                             \vspace{0.4mm}
739                             \texttt{...}\\
740                             \vspace{0.4mm}
741                             ~~~~~\texttt{NAMELIST ~~/userpar/ ~~ data\_output\_pr\_user, data\_output\_user,}\\
742                             \vspace{0.4mm}
743                             \hspace{3cm}\texttt{\textcolor{blue}{foo}, region}\\
744                             \vspace{0.4mm}
745                             \texttt{...}\\
746                             \vspace{0.4mm}
747                             ~\texttt{END SUBROUTINE user\_parin}};
748            \end{tikzpicture}
749         \end{center}
750         \vspace{1mm}
751         \item<3-> Add the NAMELIST \texttt{\&userpar} to the parameter file (e.g. \texttt{example\underline{ }cbl\underline{ }p3d}) and assign a value to this variable.         
752         \begin{center}
753%            \includegraphics[scale=0.3]{user_defined_code_figures/steering_3.png}
754            \begin{tikzpicture}
755               \node[yellow]{~\texttt{\&inipar ~~~ nx = ... ~~/}\\
756                             \vspace{0.4mm}
757                             ~\texttt{\&d3par ~~~ end\_time = 3600.0, ... ~~/}\\
758                             \par\smallskip
759                             ~\texttt{\textcolor{blue}{\&userpar ~~foo = 12345.6 ~~/}}};
760            \end{tikzpicture}
761         \end{center}
762         \vspace{1mm}
763         \item<4-> Output the variable's value using \texttt{user\underline{ }header}.
764      \end{enumerate}
765     
766
767\end{frame}
768
769% Folie 15
770\begin{frame}
771   \frametitle{User-Defined Output}
772   \scriptsize
773   \begin{itemize}
774      \item<2-> A very typical request of users is the calculation and output of quantities which are not part of PALM's standard output (e.g. a 3D-array of the resolved-scale vertical heat (temperature) flux).
775      \vspace{2mm}
776      \item<3-> The default user interface includes a number of subroutines which allow the calculation of user-defined quantities and output of these quantities as profiles, timeseries, 2d cross section or 3d volume data. These are e.g.\\
777      \texttt{user\underline{ }check\underline{ }data\underline{ }output}, \texttt{user\underline{ }check\underline{ }data\underline{ }output\underline{ }pr},\\
778      \texttt{user\underline{ }define\underline{ }netcdf\underline{ }grid}, \texttt{user\underline{ }statistics},\\
779      \texttt{user\underline{ }3d\underline{ }data\underline{ }averaging}, \texttt{user\underline{ }data\underline{ }output\underline{ }2d},\\
780      \texttt{user\underline{ }data\underline{ }output\underline{ }3d}.
781      \vspace{2mm}
782      \item<4-> The respective subroutines contain, as an example, code lines (written as comment lines) for calculating and output the square of the u-component velocity.
783      \vspace{2mm}
784      \item<5-> These quantities are output to PALM‘s standard netCDF files, i.e.\\ \texttt{DATA\underline{ }1D\underline{ }PR\underline{ }NETCDF}, \texttt{DATA\underline{ }1D\underline{ }TS\underline{ }NETCDF}, \texttt{DATA\underline{ }2D\underline{ }XY\underline{ }NETCDF} or \texttt{DATA\underline{ }3D\underline{ }NETCDF}.
785      \vspace{2mm}
786      \item<6-> The online documentation gives very detailed instructions about how to modify the interface in order to output user-defined quantities under\\
787      \vspace{1mm}
788      \textbf{http://palm.muk.uni-hannover.de/trac/wiki/doc/app/userint/output}
789   \end{itemize}
790\end{frame}
791
792% Folie 16
793\begin{frame}
794   \frametitle{User-Defined Data for Restart Runs (I)}
795   \tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\tiny]
796   \scriptsize
797   \begin{itemize}
798      \item It might be neccessary to save the values of user-defined variables at the end of a model run in order to use them for a restart run.\\
799      \ \\
800      This can be done using the routine \texttt{user\underline{ }last\underline{ }actions}.\\
801      \grqq \texttt{14}\grqq\, is the file-id for the restart file (local filename \texttt{BINOUT}):
802   \end{itemize}
803   \begin{center}
804%      \includegraphics[scale=0.3]{user_defined_code_figures/restart_1.png}
805      \begin{tikzpicture}[auto]
806         \node[yellow]{~\texttt{SUBROUTINE user\_last\_actions}\\
807                       \texttt{...}\\
808                       ~~~~~\texttt{\textcolor{blue}{WRITE ( 14 )  'foo \hspace{2cm}';  WRITE ( 14 )  foo}}\\
809                       ~~~~~\texttt{\textcolor{blue}{WRITE ( 14 )  'bar \hspace{2cm}';  WRITE ( 14 )  bar}}
810                       \par\smallskip
811                       ~~~~~\texttt{WRITE ( 14 ) '*** end user *** \hspace{0.66cm}'}
812                       \par\smallskip
813                       ~\texttt{END SUBROUTINE user\_last\_actions}};
814      \end{tikzpicture}
815   \end{center}
816\end{frame}
817
818% Folie 17
819\begin{frame}
820   \frametitle{User-Defined Data for Restart Runs (II)}
821   \tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.95\textwidth, font=\tiny]
822   \scriptsize
823   \begin{itemize}
824      \item Additionally, these variables have to be read from the restart file (file-id \grqq \texttt{13}\grqq\,, local filename \texttt{BININ}) by adding code to the routine \texttt{user\underline{ }read\underline{ }restart\underline{ }data}:
825   \end{itemize}
826   \vspace{-2.7mm}
827   \begin{center}
828%      \includegraphics[scale=0.3]{user_defined_code_figures/restart_2.png}
829      \begin{tikzpicture}
830         \node[yellow]{\texttt{~SUBROUTINE user\_read\_restart\_data( i, nxlfa, nxl\_on\_file, nxrfa, nxr\_on\_file, \&} \\
831                       \hspace{3.84cm}\texttt{nynfa, nyn\_on\_file, nysfa, nys\_on\_file,    \&} \\
832                       \hspace{3.84cm}\texttt{offset\_xa, offset\_ya, overlap\_count,       \&} \\
833                       \hspace{3.84cm}\texttt{tmp\_2d, tmp\_3d )}\\
834                       \texttt{...} \\
835                       ~~~~~\texttt{IF ( initializing\_actions == 'read\_restart\_data' )  THEN}\\
836                       ~~~~~~~~~\texttt{READ ( 13 )  field\_char} \\
837                       ~~~~~~~~~\texttt{DO  WHILE ( TRIM( field\_char )  '*** end user ***' )}\\
838                       \par\smallskip
839                       ~~~~~~~~~~~~~\texttt{nxlf = nxlfa(i,k)}
840                       \texttt{...} \\
841                       ~~~~~~~~~~~~~\texttt{SELECT CASE ( TRIM( field\_char ) )}\\
842                       \par\smallskip
843                       ~~~~~~~~~~~~~~~~~\texttt{CASE ( '\textcolor{blue}{foo}' )} \\
844                       ~~~~~~~~~~~~~~~~~~~~~\texttt{IF ( .NOT. ALLOCATED( \textcolor{blue}{foo} ) ) THEN}\\
845                       ~~~~~~~~~~~~~~~~~~~~~~~~~\texttt{ALLOCATE( \textcolor{blue}{foo}(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )}\\
846                       ~~~~~~~~~~~~~~~~~~~~~\texttt{ENDIF}\\
847                       ~~~~~~~~~~~~~~~~~~~~~\texttt{IF ( k == 1 )  READ ( 13 )  tmp\_3d}\\
848                       ~~~~~~~~~~~~~~~~~~~~~\texttt{\textcolor{blue}{foo}(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = \&} \\
849                       \hspace{4.5cm}\texttt{tmp\_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)}\\
850                       \texttt{...}\\
851                       ~~~~~~~~~~~~~\texttt{END SELECT}\\
852                       \par\smallskip
853                       ~~~~~~~~~\texttt{ENDDO}
854                       \par\smallskip
855                       ~~~~~~~~~\texttt{READ ( 13 )  field\_char}\\
856                       \par\smallskip
857                       ~~~~~~~~~\texttt{ENDDO}\\
858                       ~~~~~\texttt{ENDIF}\\
859                       \par\smallskip
860                       ~\texttt{END SUBROUTINE user\_read\_restart\_data}};
861      \end{tikzpicture}
862   \end{center}
863\end{frame}
864
865% Folie 18
866\begin{frame}
867   \frametitle{Using the User-Interface with \texttt{mrun}}
868   \tiny
869   \begin{columns}[T]
870      \begin{column}{1.1\textwidth}
871         \vspace{-1mm}
872         \textbf{Users can add their own (modified) user-interface to a PALM-run by carrying out the following steps:}
873         \begin{enumerate}
874            \item<2-> Copy the default (empty) user-interface files that you need (e.g. user\underline{ }module.f90, user\underline{ }parin.f90, user\underline{ }actions.f90) to a directory of your choice, e.g.:\\
875            \scriptsize
876            \par\smallskip
877                 \quad \texttt{cd \~{ }/palm/current\underline{ }version}\\
878                 \quad \texttt{mkdir -p USER\underline{ }CODE/example\underline{ }cbl}\\
879                 \quad \texttt{cp trunk/SOURCE/user\underline{ }module.f90  USER\underline{ }CODE/example\underline{ }cbl/user\underline{ }module.f90}\\
880                 \quad \texttt{cp trunk/SOURCE/user\underline{ }parin.f90\hspace{1em}USER\underline{ }CODE/example\underline{ }cbl/user\underline{ }parin.f90}\\
881            \quad \texttt{...}
882           
883            \item<3-> Set an additional path in the configuration file \texttt{.mrun.config} to allow \texttt{mrun} to find and include this file:\\
884            \tiny
885            \par\smallskip
886                 \scriptsize \quad \texttt{\%add\underline{ }source\underline{ }path \hspace{1em} \$base\underline{ }directory/USER\underline{ }CODE/\$fname}
887            \par\smallskip
888                 \item<4-> Modify the interface routines according to your needs.
889                 \item<5-> Start a PALM run by executing\\
890            \par\smallskip
891                 \scriptsize \quad      \texttt{mrun -d example\underline{ }cbl ...}\\
892            \par\smallskip
893                 \tiny
894                 The files user\underline{ }*.f90 will be automatically compiled within the job / interactive run and will replace the respective PALM default user-interface files.
895              \end{enumerate}
896         \begin{itemize}
897                 \item<6-> \textbf{The modified user-interface file cannot be pre-compiled by using} \texttt{mbuild}!
898                 \item<7-> The above method allows to use different user-interfaces for different runs. Just store the respective interface-files in subdirectories \texttt{USER\underline{ }CODE/abcd}, \texttt{USER\underline{ }CODE/cdef}, etc. and start \texttt{mrun} with option \grqq \texttt{-d abcd}\grqq, \grqq \texttt{-d cdef}\grqq, etc.
899         \end{itemize}
900
901      \end{column}
902   \end{columns} 
903   
904\end{frame}
905
906\end{document}
Note: See TracBrowser for help on using the repository browser.