Home

Context Navigation

user_defined_code.tex @ 1734

Last change on this file since 1734 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

Rev	Line
[973]	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}
[1226]	16	%\usetikzlibrary{decorations.markings} %neues paket
	17	%\usetikzlibrary{decorations.pathreplacing} %neues paket
[973]	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
[1515]	28	\institute{Institute of Meteorology and Climatology, Leibniz UniversitÃ€t Hannover}
	29	\selectlanguage{english}
[973]	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}
[1515]	49	\author{PALM group}
[973]	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}
[1515]	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.
[973]	73	\end{itemize}
	74	\vspace{3mm}
	75	\onslide<4->
	76	\textbf{Solution:}
	77	\begin{itemize}
[1515]	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.
[973]	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}
[1226]	87	\tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\Tiny]
[973]	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}.\\
	92	Example: \texttt{user\underline{ }last\underline{ }actions.f90}
	93	\end{itemize}
	94	\onslide<3->
[1226]	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}\\
[1515]	107	~~~~\texttt{USE kinds}\\
[1226]	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}\\
[973]	132	\end{frame}
	133
	134	% Folie 4
	135	\begin{frame}
	136	\frametitle{Embedding of User-Interface Routines}
[1226]	137	\tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\Tiny]
[973]	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->
[1226]	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}\\
[973]	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]
[1226]	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]
[973]	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
[1226]	208	\node[yellow, text width=5em] (PALM) at (0,0.1) {PALM};
[973]	209
	210	\coordinate (A) at (0,-0.7);
[1226]	211	\node[yellow, right=0.6cm of A] (parin) {parin};
[973]	212
[1226]	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);
[973]	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]
[1226]	289	\tikzstyle{green1} = [draw, fill=green!30, text width=10em, text centered, minimum height=1.6em]
[973]	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
[1226]	307	\node[green, right=0.6cm of B, text width=14em] (useractions1) {\textcolor{blue}{user\_actions\_(\textit{'before\_timestep'})}};
[973]	308
[1234]	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
[1226]	312	\node[green, right=0.3cm of progframe, text width=14em] (useractions2) {\textcolor{blue}{user\_actions\_(\textit{'*\_tendency'})}};
[973]	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
[1226]	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};
[973]	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
[1226]	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};
[973]	342
[1226]	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};
[973]	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
[1226]	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};
[973]	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);
[1234]	367	\draw[-, thick] (D) -- (E);
	368	\draw[->, thick] (E) -- (progeqfast);
[973]	369	\draw[-, thick] (progeq) -- (D);
	370	\draw[-, thick] (progframe) -- (useractions2);
[1226]	371	\draw[-, thick] (F) -- (lpm);
	372	\draw[-, thick] (userlpmadvec) -- (lpm);
[973]	373	\draw[-, thick] (G) -- (exchangehoriz);
	374	\draw[-, thick] (H) -- (swaptimelevel);
	375	\draw[-, thick] (I) -- (pres);
	376	\draw[dotted, thick] (pres) -- (J);
[1226]	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);
[973]	384	\draw[-, thick] (O) -- (useractions4);
[1226]	385	% \draw[-, thick] (O) -- (timestep);
[973]	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
[1226]	417	\node[yellow, right=0.6cm of D] (write3dbinary) {write\_3d\_binary};
	418	\node[yellow, right=0.9cm of write3dbinary] (writevarlist) {write\_var\_list};
[973]	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);
[1226]	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};
[973]	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
[1226]	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\\
[973]	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
[1226]	455	\texttt{user\_check\_data\_output} & \texttt{variable}, \texttt{unit} & \texttt{check\_parameters} + \texttt{init\_masks}& check the user-defined output quantities\\
[973]	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
[1226]	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\\
[973]	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\\
[1226]	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
[973]	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
[1226]	492	\texttt{user\_init\_grid} & \texttt{gls} & \texttt{init\_grid} & defining a special topography\\
[973]	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\\
[1226]	497	\hline
	498	\texttt{user\_lpm\_advec} & - - - & \texttt{lpm} & modification of particles after advection \\
	499	\hline
[1515]	500	\texttt{user\_lpm\_init} & - - - & \texttt{lpm\_init} & defining initial particle sources \\
[1226]	501	\hline
	502	\texttt{user\_lpm\_set\_attributes} & - - - & \texttt{lpm} & defining particles attributes \\
	503	\hline
[973]	504	\texttt{MODULE user (user\_module.f90)} & - - - & - - - & contains user defined variables\\
	505	\hline
	506	\texttt{user\_parin} & & \texttt{parin} & reading user variables\\
[1226]	507	\hline
[973]	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
[1226]	510	\texttt{user\_spectra} & \texttt{mode}, \texttt{m}, \texttt{pr} & \texttt{calc\_spectra} + \texttt{data\_output\_spectra} & output/calculation of additional user-defined quantities\\
[973]	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 \\
[1649]	518	\textbf{http://palm.muk.uni-hannover.de/trac/wiki/doc/app/userint}\\
[973]	519	for detailed explanations.
	520
	521	\end{frame}
	522
	523	% Folie 10
	524	\begin{frame}
	525	\frametitle{Data Access / Exchange}
[1226]	526	\tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.6\textwidth, font=\tiny]
[973]	527	\footnotesize
	528	\begin{itemize}
	529	\item<2-> \textbf{Between the standard PALM code and the user-interface:}\\
	530	\begin{itemize}
	531	\scriptsize
[1515]	532	\item<3-> by including the respective \textcolor{blue}{PALM modules} in the user-interface subroutines.\\
[973]	533	\begin{center}
[1226]	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}
[973]	552	\end{center}
	553	\end{itemize}
	554	\item<4-> \textbf{Within the user-interface:}
	555	\begin{itemize}
	556	\scriptsize
[1515]	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).}
[973]	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)}
[1226]	566	\tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\tiny]
[973]	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.
[1226]	570	\vspace{1.8mm}
[973]	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	\ \\
[1515]	573	\quad \texttt{CALL user\underline{ }actions( 'u-tendency' )}\\
[973]	574	\ \\
[1515]	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:
[973]	576	\end{itemize}
[1226]	577	\vspace{-3mm}
[973]	578	\begin{center}
[1226]	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}
[973]	596	\end{center}
[1226]	597	% \begin{center}
	598	% \includegraphics[scale=0.33]{user_defined_code_figures/user_actions_1.png}
	599	% \end{center}
[973]	600	\end{frame}
	601
	602	% Folie 12
	603	\begin{frame}
	604	\frametitle{Usage of \texttt{user\underline{ }actions} (II)}
[1226]	605	\tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\tiny]
[973]	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}
[1226]	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}\\
[1531]	618	~~~~~~\texttt{SELECT CASE ( location )}\\
[1226]	619	\texttt{...}\\
	620	~~~~~~~~~~\texttt{CASE ( 'after\_timestep' )}\\
	621	\texttt{!}\\
[1531]	622	\texttt{!-- \hspace{6mm}Enter actions to be done after every timestep here}\\
[1226]	623	\par\medskip
	624	~~~~~~~~~~\texttt{CASE ( 'u-tendency' )}\\
	625	\texttt{!}\\
[1531]	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}}\\
[1226]	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}
[973]	639	\end{center}
	640	\end{frame}
	641
	642	% Folie 13
	643	\begin{frame}[t]
	644	\frametitle{Usage of \texttt{user\underline{ }actions} (III)}
[1226]	645	\tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.5\textwidth, font=\Tiny]
[973]	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}:
[1515]	651	\item<2-> From \texttt{prognostic\underline{ }equations\underline{ }vector}: \hspace{0.6em}\texttt{CALL user\underline{ }actions( 'u-tendency' )}
[973]	652	\item<3-> From \texttt{prognostic\underline{ }equations},\\
[1515]	653	\texttt{prognostic\underline{ }equations\underline{ }cache}:\hspace{4em}\texttt{CALL user\underline{ }actions( i, j, 'u-tendency' )}
[973]	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\\
[1515]	663	i and j are executed in\\
[973]	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]
[1226]	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
[1531]	673	~~~~~~\texttt{USE control\_parameters}\\
	674	~~~~~~\texttt{USE pegrid}\\
	675	~~~~~~\texttt{USE user}\\
[1226]	676	\par\medskip
[1531]	677	~~~~~~\texttt{IMPLICIT NONE}\
[1226]	678	\par\medskip
[1531]	679	~~~~~~\texttt{CHARACTER (LEN=*) :: location}\\
[1226]	680	\par\medskip
[1531]	681	~~~~~~\texttt{INTEGER(iwp) :: i, idum, j}\\
[1226]	682	\par\medskip
	683	\texttt{!}\\
	684	\texttt{!-- Here the user-defined actions follow}\\
	685	\vspace{0.5mm}
[1531]	686	~~~~~~\texttt{SELECT CASE ( location )}\\
[1226]	687	\par\medskip
[1531]	688	~~~~~~~~~\texttt{CASE ( 'u-tendency' )}\\
[1226]	689	\texttt{!}\\
[1531]	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}}\\
[1226]	694	\par\medskip
[1531]	695	~~~~~~~~~\texttt{CASE ( 'v-tendency' )}};
[1226]	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	% };
[973]	701	\end{tikzpicture}}
	702
	703	\end{frame}
	704
	705	% Folie 14
	706	\begin{frame}
	707	\frametitle{Steering the User-Interface}
[1226]	708	\tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\Tiny]
[973]	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}
[1226]	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}
[1515]	725	~\texttt{REAL(wp) \hspace{0.4cm}:: ~~\textcolor{blue}{foo = 0.0}}\\
[1226]	726	\vspace{0.4mm}
	727	\texttt{...}\\
	728	\vspace{0.4mm}
	729	~\texttt{END MODULE user}};
	730	\end{tikzpicture}
[973]	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}
[1226]	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}
[973]	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}
[1226]	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}
[973]	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}
[1649]	788	\textbf{http://palm.muk.uni-hannover.de/trac/wiki/doc/app/userint/output}
[973]	789	\end{itemize}
	790	\end{frame}
	791
	792	% Folie 16
	793	\begin{frame}
	794	\frametitle{User-Defined Data for Restart Runs (I)}
[1226]	795	\tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.7\textwidth, font=\tiny]
[973]	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}
[1226]	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}
[973]	815	\end{center}
	816	\end{frame}
	817
	818	% Folie 17
	819	\begin{frame}
	820	\frametitle{User-Defined Data for Restart Runs (II)}
[1226]	821	\tikzstyle{yellow} = [rectangle, draw, fill=yellow!30, text width=0.95\textwidth, font=\tiny]
[973]	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}
[1226]	826	\vspace{-2.7mm}
[973]	827	\begin{center}
[1226]	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}
[973]	862	\end{center}
	863	\end{frame}
	864
	865	% Folie 18
	866	\begin{frame}
	867	\frametitle{Using the User-Interface with \texttt{mrun}}
[1515]	868	\tiny
[973]	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.:\\
[1515]	875	\scriptsize
[1226]	876	\par\smallskip
[973]	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{...}
[1515]	882
[973]	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
[1515]	885	\par\smallskip
	886	\scriptsize \quad \texttt{\%add\underline{ }source\underline{ }path \hspace{1em} \$base\underline{ }directory/USER\underline{ }CODE/\$fname}
	887	\par\smallskip
[973]	888	\item<4-> Modify the interface routines according to your needs.
	889	\item<5-> Start a PALM run by executing\\
[1515]	890	\par\smallskip
	891	\scriptsize \quad \texttt{mrun -d example\underline{ }cbl ...}\\
	892	\par\smallskip
[973]	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
[1234]	906	\end{document}

Note: See TracBrowser for help on using the repository browser.

Context Navigation

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

Download in other formats: