source: palm/trunk/TUTORIAL/SOURCE/mrun_data_analysis.tex @ 1033

% $Id: mrun_data_analysis.tex 973 2012-08-07 16:03:47Z letzel $
\input{header_tmp.tex}
%\input{../header_lectures.tex}

\usepackage[utf8]{inputenc}
\usepackage{ngerman}
\usepackage{pgf}
\usetheme{Dresden}
\usepackage{subfigure}
\usepackage{units}
\usepackage{multimedia}
\usepackage{hyperref}
\newcommand{\event}[1]{\newcommand{\eventname}{#1}}
\usepackage{xmpmulti}
\usepackage{tikz}
\usetikzlibrary{shapes,arrows,positioning}
\usetikzlibrary{decorations.markings}             %neues paket
\usetikzlibrary{decorations.pathreplacing}        %neues paket
\def\Tiny{\fontsize{4pt}{4pt}\selectfont}
\usepackage{amsmath}
\usepackage{amssymb}
\usepackage{multicol}
\usepackage{pdfcomment}
\usepackage{graphicx}
\usepackage{listings}
\lstset{showspaces=false,language=fortran,basicstyle=
        \ttfamily,showstringspaces=false,captionpos=b}

\institute{Institut fÌr Meteorologie und Klimatologie, Leibniz UniversitÀt Hannover}
\date{last update: \today}
\event{PALM Seminar}
\setbeamertemplate{navigation symbols}{}

\setbeamertemplate{footline}
  {
    \begin{beamercolorbox}[rightskip=-0.1cm]&
     {\includegraphics[height=0.65cm]{imuk_logo.pdf}\hfill \includegraphics[height=0.65cm]{luh_logo.pdf}}
    \end{beamercolorbox}
    \begin{beamercolorbox}[ht=2.5ex,dp=1.125ex,
      leftskip=.3cm,rightskip=0.3cm plus1fil]{title in head/foot}
      {\leavevmode{\usebeamerfont{author in head/foot}\insertshortauthor} \hfill \eventname \hfill \insertframenumber \; / \inserttotalframenumber}
    \end{beamercolorbox}
    \begin{beamercolorbox}[colsep=1.5pt]{lower separation line foot}
    \end{beamercolorbox}
  }
%\logo{\includegraphics[width=0.3\textwidth]{luhimuk_logo.pdf}}

\title[Data Analysis]{Data Analysis}
\author{Siegfried Raasch}

\begin{document}

% Folie 1
\begin{frame}
\titlepage
\end{frame}

\section{Data Analysis}
\subsection{Data Analysis}

% Folie 2
\begin{frame}
   \frametitle{PALM Data Analysis}
   
   \footnotesize
   PALM  offers different ways of standard output to analyze the simulation data: 
   \begin{enumerate}
      \item<2-> Run control informations (for certain timesteps)
      \item<3-> Mean vertical profiles
      \begin{itemize}
         \footnotesize
         \item[-]<3-> horizontally averaged, instantaneous or time-averaged
      \end{itemize}
      \item<4-> Time series
      \item<5-> Two-dimensional cross sections ($xy$, $xz$, $yz$)
      \begin{itemize}
         \footnotesize
         \item[-]<5-> instantaneous or time-averaged, + averaged along 3rd dimension
      \end{itemize}
      \item<6-> Complete three-dimensional data arrays
      \begin{itemize}
         \footnotesize
         \item[-]<6-> instantaneous or time-averaged
      \end{itemize}
   \end{enumerate}
   \quad \onslide<7-> Others (spectra, particle data)\\
   \vspace{3mm}
   \onslide<8-> First item is pure ASCII-data to give a short overview of the run. All other output is in netCDF format and can be graphically displayed.\\
   \vspace{2mm}
   Additionally, power spectra can be output using the spectra-package (use \texttt{mrun}-option \grqq \texttt{-p spectra}\grqq).
   
\end{frame}

% Folie 3
\begin{frame}
   \tikzstyle{yellow} = [rectangle, fill=yellow!20, font=\scriptsize, minimum height=1.2cm]
   \tikzstyle{purple} = [rectangle, fill=magenta!10, font=\scriptsize, minimum height=0.9cm]
   \tikzstyle{white} = [rectangle, font=\tiny, minimum height=0.9cm]
      
   \frametitle{Processing Standard Output Data with Graphics Software (Public Domain-Software)}
   
   \footnotesize
   The standard output files are in netCDF format which can be easily displayed with \texttt{ncview}, \texttt{ncl}, or \texttt{ferret}.
   
   \begin{tikzpicture}[auto, node distance=0]
   
      % Zeile 1
      \node [yellow, text width=0.2\textwidth] (1) {standard output};
      \node [yellow, text width=0.25\textwidth, right=0.0cm of 1] (2) {names of local files created by PALM};   
      \node [yellow, text width=0.25\textwidth, right=0.0cm of 2] (3) {suffix of permanent files as defined in \texttt{.mrun.config}};   
      \node [yellow, text width=0.25\textwidth, right=0.0cm of 3] (4) {call of graphics software};   
      
      % Zeile 2      
      \node [purple, text width=0.2\textwidth, below=0.0cm of 1] (12) {mean vertical profiles};
      \node [white, text width=0.25\textwidth, below=0.0cm of 2] (22) {\texttt{DATA\_1D\_PR\_NETCDF}};
      \node [white, text width=0.25\textwidth, below=0.0cm of 3] (32) {\texttt{\_pr.nc}};
      \node [white, text width=0.25\textwidth, below=0.0cm of 4] (42) {\texttt{ncview <filename>}\\
      \texttt{palmplot pr <options>}\\
      \texttt{ferret -gui}};
      
      % Zeile 3      
      \node [purple, text width=0.2\textwidth, below=0.0cm of 12] (13) {time series};
      \node [white, text width=0.25\textwidth, below=0.0cm of 22] (23) {\texttt{DATA\_1D\_TS\_NETCDF}};
      \node [white, text width=0.25\textwidth, below=0.0cm of 32] (33) {\texttt{\_ts.nc}};
      \node [white, text width=0.25\textwidth, below=0.0cm of 42] (43) {\texttt{ncview <filename>}\\
      \texttt{palmplot ts <options>}\\
      \texttt{ferret -gui}};
      
      % Zeile 4           
      \node [purple, text width=0.2\textwidth, minimum height=1.5cm, below=0.0cm of 13] (14) {2d-cross-sections};
      \node [white, text width=0.25\textwidth, minimum height=1.5cm, below=0.0cm of 23] (24) {\texttt{DATA\_2D\_XY\_NETCDF}\\ \texttt{DATA\_2D\_XZ\_NETCDF}\\ \texttt{DATA\_2D\_YZ\_NETCDF}};
      \node [white, text width=0.25\textwidth, minimum height=1.5cm, below=0.0cm of 33] (34) {\texttt{\_xy.nc}\\ \texttt{\_xz.nc}\\ \texttt{\_yz.nc}};
      \node [white, text width=0.25\textwidth, minimum height=1.5cm, below=0.0cm of 43] (44) {\texttt{ncview <filename>}\\
      \texttt{palmplot xy <options>}\\
      \texttt{palmplot xz <options>}\\
      \texttt{palmplot yz <options>}\\
      \texttt{ferret -gui}};
                  
      % Zeile 5      
      \node [purple, text width=0.2\textwidth, below=0.0cm of 14] (15) {3d-data}; 
      \node [white, text width=0.25\textwidth, below=0.0cm of 24] (25) {\texttt{DATA\_3D\_NETCDF}};
      \node [white, text width=0.25\textwidth, below=0.0cm of 34] (35) {\texttt{\_3d.nc}};  
      \node [white, text width=0.25\textwidth, below=0.0cm of 44] (45) {\texttt{ncview <filename>}\\
      \texttt{palmplot xy <options>}\\
      \texttt{ferret -gui}};
   
      \draw[-] (1.south west) -- (4.south east);   
      \draw[-] (13.south west) -- (43.south east); 
      \draw[-] (1.north east) -- (15.south east); 
   \end{tikzpicture}
\end{frame}

% Folie 4
\begin{frame}
   \frametitle{User Actions Necessary to Get Output}
   Example for the case of mean vertical profiles:
   \footnotesize
   \begin{enumerate}
   
      \item<2-> \textbf{Carefully think about, which variables you want to analyze}
      \begin{itemize}
         \scriptsize
              \item[]<2-> e.g. temperature, wind, etc., resolved- or subgrid-scale?, output interval?, should they be time averaged?
      \end{itemize}
      \vspace{1mm}
      \item<3-> \textbf{Find and set the appropriate output parameters in the parameter file}
      \begin{itemize}
         \scriptsize
              \item[]<3-> \texttt{data\_output\_pr = 'pt', 'u', 'w \dq pt\dq\, ',  dt\_dopr = 900.0,}
              \item[]<3-> \texttt{averaging\_interval\_pr = 600.0, dt\_averaging\_input\_pr = 10.0}
      \end{itemize}    
      \vspace{1mm}
      \item<4-> \textbf{Check, if the necessary file connection statements exist in the .mrun.config file}\\
      \begin{itemize}
         \scriptsize
              \item[]<4-> \texttt{DATA\_1D\_PR\_NETCDF   out:loc:tr  pr\#      \~{}/palm/current\_version/JOBS/\$fname/OUTPUT  \_pr  nc}
      \end{itemize}  
      \vspace{1mm}
      \item<5-> \textbf{Don't forget to activate the file connection statements by setting the correct activation string(s) in \texttt{mrun}-option \texttt{-r}}\\
      \begin{itemize}
         \scriptsize
              \item[]<5-> \texttt{mrun ... -r \dq ... pr\# ...\dq\, ...}
      \end{itemize}   
      \vspace{1mm}
      \item<6-> \textbf{After job completion, change to the OUTPUT-directory and call the appropriate graphics software}\\
      \begin{itemize}
         \scriptsize
              \item[]<6-> \texttt{ncview <filename>} \quad \textbf{or} \quad \texttt{palmplot pr <options>}
      \end{itemize}  
   \end{enumerate}   

\end{frame}

\end{document}