Changeset 89


Ignore:
Timestamp:
May 25, 2007 12:08:31 PM (18 years ago)
Author:
raasch
Message:

further changes concerning user-defined profiles

Location:
palm/trunk
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/DOC/app/chapter_3.5.4.html

    r88 r89  
    1616section or 3d volume data. The respective subroutines
    1717contain sample code lines (written as comment lines) for defining, calculating and
    18 output of such quantities.<br><br>The
     18output of such quantities.<br><br>Output times, averaging intervals, etc. are steered by the same variables as used for the standard PALM output quantities, e.g. <a href="chapter_4.2.html#dt_data_output">dt_data_output</a>.<br><br>The
    1919rest of this chapter explains step-by-step how to modify/extend the
    20 default file user_interface.f90 in order to generate the respective output.<br><br><br><h4>Output of vertical profiles</h4><br><br><h4><span style="font-weight: bold;">Output of 2d cross sections or 3d volume data</span></h4><br>This example shows the output of the
     20default file user_interface.f90 in order to generate the respective output.<br><br><br><h4>Output of vertical profiles</h4>This example shows the output of the
     21quantity "turbulent resolved-scale horizontal momentum flux" (u*v*). If more than one user-defined
     22quantity shall be output, the following steps have to be carried out in the
     23same way for each of the quantities.<br><br><ol><li>The
     24quantity has to be given a unique string identifier, e.g. <span style="font-style: italic;">'u*v*'</span>.
     25This identifier must be different from the identifiers used for the
     26PALM standard output (see list in description of parameter <a href="chapter_4.2.html#data_output_pr">data_output_pr</a>).
     27To switch on output of this quantity, the user has to assign the string
     28identifier to the parameter <a href="chapter_4.3.html#data_output_pr_user">data_output_pr_user</a>,
     29eg.:<br><br><span style="font-family: monospace;">&nbsp;
     30&nbsp; data_output_pr_user</span> = <span style="font-style: italic;">'u*v*'</span>,<span style="font-style: italic;"></span><span style="font-style: italic;"></span><span style="font-style: italic;"></span><span style="font-style: italic;"></span><span style="font-style: italic;"></span><span style="font-style: italic;"></span><br><br></li><li>For
     31the
     32quantity, an identification number, a physical unit, and the vertical
     33grid on which it is defined (u- or w-grid), has to be assigned (subroutine <a href="chapter_3.5.1.html#user_check_data_output_pr"><span style="font-family: monospace;">user_check_data_output_pr</span></a>):<span style="font-family: monospace;"><br><br>&nbsp; &nbsp; CASE (
     34'u*v*' )<br></span><span style="font-family: monospace;">&nbsp; &nbsp; &nbsp; &nbsp;index = pr_palm + <span style="color: rgb(255, 0, 0);">1<span style="color: rgb(0, 0, 0);">
     35&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
     36&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
     37&nbsp; &nbsp; ! identification number</span></span><br>&nbsp; &nbsp; &nbsp; &nbsp;dopr_index(var_count) = index<br>&nbsp; &nbsp; &nbsp;&nbsp; dopr_unit(var_count) &nbsp;= '<span style="color: rgb(255, 0, 0);">m2/s2</span>'
     38&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
     39&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ! physical unit<br>&nbsp; &nbsp; &nbsp; &nbsp;hom(:,2,index,:) &nbsp; &nbsp; &nbsp;= SPREAD( <span style="color: rgb(255, 0, 0);">zu</span>, 2, statistic_regions+1 ) &nbsp; &nbsp;! vertical grid<br></span><br>Here only the those parts in <span style="color: rgb(255, 0, 0);">red</span> color have to be given by the user appropriately.<br>The identification number (<span style="font-family: Courier New,Courier,monospace;">index</span>) must be within the range [<span style="font-family: Courier New,Courier,monospace;"> pr_palm+1 , pr_palm+max_pr_user</span> ], where <span style="font-family: Courier New,Courier,monospace;">max_pr_user</span> is the number of&nbsp;user-defined profiles as given by parameter <a href="chapter_4.3.html#data_output_pr_user">data_output_pr_user</a> in the respective PALM run.&nbsp;The physical unit has to be given with respect to the NetCDF conventions. If no unit is given,
     40PALM will abort. The vertical grid has to be either <span style="font-family: Courier New,Courier,monospace;">zu</span> (u-grid) or <span style="font-family: Courier New,Courier,monospace;">zw</span> (w-grid).<br><br></li><li>The quantity has to be calculated for all gridpoints (subroutine <a href="chapter_3.5.1.html#user_statistics">user_statistics</a>):<br><br><span style="font-family: monospace;">&nbsp; &nbsp; !$OMP DO</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp; DO&nbsp; i = nxl, nxr</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DO&nbsp; j = nys, nyn</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DO&nbsp; k = nzb_s_outer(j,i)+1, nzt</span><br style="font-family: monospace;"><span style="font-family: monospace;"></span><span style="font-family: monospace;"></span><span style="font-family: monospace;"></span><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sums_l(k,pr_palm+<span style="color: rgb(255, 0, 0);">1</span>,tn) = sums_l(k,pr_palm+<span style="color: rgb(255, 0, 0);">1</span>,tn) +&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &amp;</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <span style="color: rgb(255, 0, 0);">( 0.5 * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) * &amp;</span></span><br style="font-family: monospace; color: rgb(255, 0, 0);"><span style="font-family: monospace; color: rgb(255, 0, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
     41( 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) ) * &amp;</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
     42* rmask(j,i,sr)<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ENDDO<br>&nbsp; &nbsp; &nbsp; &nbsp;ENDDO<br>&nbsp; &nbsp; ENDDO<br></span><br style="font-family: monospace;">Once again, only those parts in <span style="color: rgb(255, 0, 0);">red</span> have to be adjusted by the user.<br>The
     43turbulent resolved-scale momentum flux u*v* is defined as the product
     44of the deviations of the horizontal velocities from their respective
     45horizontally averaged mean values. These mean values are stored in
     46array <span style="font-family: Courier New,Courier,monospace;">hom(..,1,1,sr)</span> and <span style="font-family: Courier New,Courier,monospace;">hom(..,1,2,sr)</span> for the u- and v-component, respectively. Since due to the staggered grid, <span style="font-family: Courier New,Courier,monospace;">u</span> and <span style="font-family: Courier New,Courier,monospace;">v</span>
     47are not defined at the same gridpoints, they have to be interpolated
     48appropriately (here to the center of the gridbox). The result of the
     49calculation is stored in array <span style="font-family: Courier New,Courier,monospace;">sums_l</span>.
     50The second index of this array is the identification number of the
     51profile which must match the one given in the previous step 2.</li></ol><br><h4>Output of timeseries</h4>still to be added<br>&nbsp;<br><h4><span style="font-weight: bold;">Output of 2d cross sections or 3d volume data</span></h4><br>This example shows the output of the
    2152quantity "square of the u-component" (<span style="font-weight: bold;">Note:</span> this quantity
    2253could of course easily be calculated from the u-component by
    2354postprocessing the PALM output so that calculation within PALM is not
    2455necessarily required). If more than one user-defined
    25 quantity shall be output, these steps have to be carried out in the
     56quantity shall be output, the following steps have to be carried out in the
    2657same way for each of the quantities.<br><br><ol><li>The
    2758quantity has to be given a unique string identifier, e.g. <span style="font-style: italic;">'u2'</span>.
  • palm/trunk/DOC/app/chapter_4.2.html

    r75 r89  
    788788</ul> <p>These initial profiles have been either set by
    789789the user or
    790 have been calculated by a 1d-model prerun.<br> </p>In case
     790have been calculated by a 1d-model prerun.</p>The
     791user is allowed to extend the above list of quantities by defining his
     792own output quantities (see the user-parameter <a href="chapter_4.3.html#data_output_pr_user">data_output_pr_user</a>).<br><br>In case
    791793of ASCII data output to local file PLOT1D_DATA,
    792794PALM additionally creates a NAMELIST parameter file (local name <a href="chapter_3.4.html#PLOT1D_PAR">PLOT1D_PAR</a>)
  • palm/trunk/DOC/app/chapter_4.3.html

    r62 r89  
    11<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
    22<html><head>
    3 <meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM
    4 chapter 4.3</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20041013;13364831"> <meta name="CHANGED" content="20041015;15380006"> <meta name="KEYWORDS" content="parallel LES model"> <style>
     3<meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM chapter 4.3</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20041013;13364831"> <meta name="CHANGED" content="20041015;15380006"> <meta name="KEYWORDS" content="parallel LES model"> <style>
    54<!--
    65@page { size: 21cm 29.7cm }
     
    1413<br> <b><font size="4">value</font></b></p>
    1514</td> <td> <p><font size="4"><b>Explanation</b></font></p>
    16 </td> </tr> <tr><td style="vertical-align: top;"><a name="data_output_user"></a><span style="font-weight: bold;">data_output_user</span></td><td style="vertical-align: top;">C * 10 (100)</td><td style="vertical-align: top;"><span style="font-style: italic;">100 * ' '</span></td><td style="vertical-align: top;">User defined quantities for
     15</td> </tr> <tr><td style="vertical-align: top;"><span style="font-weight: bold;"><a name="data_output_pr_user"></a>data_output_pr_user</span></td><td style="vertical-align: top;">C * 10<br>(200)</td><td style="vertical-align: top;"><span style="font-style: italic;">200 * ' '</span></td><td style="vertical-align: top;">User defined quantities for
     16which horizontally averaged profile data is to be output.<br><br>Beside
     17the PALM standard output quantities (which can be selected via
     18parameter <a href="chapter_4.2.html#data_output_pr">data_output_pr</a>),
     19the user can define (without any restriction) additional output
     20quantities (e.g. turbulent resolved-scale horizontal momentum fluxes, etc.). Each of these
     21quantities has to be given a unique identifier (string) which must be
     22different from the strings defining the standard output quantities (see
     23list from the description of <a href="chapter_4.2.html#data_output_pr">data_output_pr</a>).
     24Data output can be switched on by assigning the respective strings to <span style="font-weight: bold;">data_output_user</span>.<br><br>The
     25user has to calculate/provide the respective data array(s) by
     26appropriately extending the user interface (see <a href="chapter_3.5.4.html">3.5.4</a>).<br><br>For
     27further steering of the user defined output (output times, cross
     28section levels, time averaging, etc.), the steering parameters for <a href="chapter_4.2.html#data_output_pr">data_output_pr</a>
     29apply accordingly.</td></tr><tr><td style="vertical-align: top;"><a name="data_output_user"></a><span style="font-weight: bold;">data_output_user</span></td><td style="vertical-align: top;">C * 10 (100)</td><td style="vertical-align: top;"><span style="font-style: italic;">100 * ' '</span></td><td style="vertical-align: top;">User defined quantities for
    1730which 2d cross section and/or 3d volume data are to be output.<br><br>Beside
    1831the PALM standard output quantities (which can be selected via
  • palm/trunk/DOC/app/chapter_4.6.html

    r83 r89  
    283283<td>Quantities for which vertical profiles (horizontally
    284284averaged)
    285 are to be output.</td> </tr> <tr> <td><a href="chapter_4.2.html#data_output_sp"><span style="font-weight: bold;">data_output_sp</span></a></td>
     285are to be output.</td> </tr> <tr><td><a href="chapter_4.3.html#data_output_pr_user"><span style="font-weight: bold;">data_output_pr_user</span></a></td><td>U</td><td>C * 10<br>(200)</td><td><span style="font-style: italic;">200 * ' '</span></td><td>User defined quantities for which horizontally averaged profile data is to be output.</td></tr><tr> <td><a href="chapter_4.2.html#data_output_sp"><span style="font-weight: bold;">data_output_sp</span></a></td>
    286286<td>P</td> <td>C * 10 (10)</td> <td><span style="font-style: italic;">10 * ' '</span></td>
    287287<td>Quantities for which horizontal spectra are to be calculated
  • palm/trunk/DOC/app/index.html

    r83 r89  
    154154of contents</span><span style="font-family: Thorndale;" lang="EN-GB"> <br>
    155155&nbsp; <o:p></o:p></span></p>
    156 <p style="margin-left: 36pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_1.0.html"><span style="" lang="EN-GB">1.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Introduction <o:p></o:p></span></p>
    157 <p style="margin-left: 36pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_2.0.html"><span style="" lang="EN-GB">2.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Basic
     156<p style="margin-left: 36pt;"><span style="font-family: Thorndale;"><a href="chapter_1.0.html"><span style="" lang="EN-GB">1.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Introduction <o:p></o:p></span></p>
     157<p style="margin-left: 36pt;"><span style="font-family: Thorndale;"><a href="chapter_2.0.html"><span style="" lang="EN-GB">2.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Basic
    158158techniques of the LES model and its parallelization <o:p></o:p></span></p>
    159 <p style="margin: 5pt 0cm 0.0001pt 36pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.0.html"><span style="" lang="EN-GB">3.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Execution of
     159<p style="margin: 5pt 0cm 0.0001pt 36pt;"><span style="font-family: Thorndale;"><a href="chapter_3.0.html"><span style="" lang="EN-GB">3.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Execution of
    160160model runs<o:p></o:p></span></p>
    161 <p style="margin: 0cm 0cm 0.0001pt 72pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.1.html"><span style="" lang="EN-GB">3.1</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Job control
     161<p style="margin: 0cm 0cm 0.0001pt 72pt;"><span style="font-family: Thorndale;"><a href="chapter_3.1.html"><span style="" lang="EN-GB">3.1</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Job control
    162162with <span class="SpellE"><b>mrun</b></span> <br>
    163 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.2.html"><span style="" lang="EN-GB">3.2</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Example of a
     163</span><span style="font-family: Thorndale;"><a href="chapter_3.2.html"><span style="" lang="EN-GB">3.2</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Example of a
    164164minimum configuration file <br>
    165 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.3.html"><span style="" lang="EN-GB">3.3</span></a></span><span style="font-family: Thorndale;" lang="EN-GB">
     165</span><span style="font-family: Thorndale;"><a href="chapter_3.3.html"><span style="" lang="EN-GB">3.3</span></a></span><span style="font-family: Thorndale;" lang="EN-GB">
    166166Initialization and restart runs <br>
    167 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.4.html"><span style="" lang="EN-GB">3.4</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Input and
     167</span><span style="font-family: Thorndale;"><a href="chapter_3.4.html"><span style="" lang="EN-GB">3.4</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Input and
    168168output files <br>
    169 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.5.html"><span style="" lang="EN-GB">3.5</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Programming
     169</span><span style="font-family: Thorndale;"><a href="chapter_3.5.html"><span style="" lang="EN-GB">3.5</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Programming
    170170of user-defined code extensions <o:p></o:p></span></p>
    171 <p style="margin: 0cm 0cm 0.0001pt 108pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.5.1.html"><span style="" lang="EN-GB">3.5.1</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Interfaces for
     171<p style="margin: 0cm 0cm 0.0001pt 108pt;"><span style="font-family: Thorndale;"><a href="chapter_3.5.1.html"><span style="" lang="EN-GB">3.5.1</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Interfaces for
    172172user-defined&nbsp; code <br>
    173 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.5.2.html"><span style="" lang="EN-GB">3.5.2</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Definition of
     173</span><span style="font-family: Thorndale;"><a href="chapter_3.5.2.html"><span style="" lang="EN-GB">3.5.2</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Definition of
    174174user-defined parameters <br>
    175 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.5.3.html"><span style="" lang="EN-GB">3.5.3</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Definition of
     175</span><span style="font-family: Thorndale;"><a href="chapter_3.5.3.html"><span style="" lang="EN-GB">3.5.3</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Definition of
    176176<span class="SpellE">user-defined subdomains</span></span></p><p style="margin: 0cm 0cm 0.0001pt 108pt;"><span style="font-family: Thorndale;" lang="EN-GB"><span class="SpellE"></span> <a href="chapter_3.5.4.html">3.5.4</a> User-defined output quantities<br>
    177 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.5.5.html"><span style="" lang="EN-GB">3.5.5</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Compiling and linking
     177</span><span style="font-family: Thorndale;"><a href="chapter_3.5.5.html"><span style="" lang="EN-GB">3.5.5</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Compiling and linking
    178178user-defined code <o:p></o:p></span></p>
    179 <p style="margin: 0cm 0cm 5pt 72pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.6.html"><span style="" lang="EN-GB">3.6</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Interactive mode and
     179<p style="margin: 0cm 0cm 5pt 72pt;"><span style="font-family: Thorndale;"><a href="chapter_3.6.html"><span style="" lang="EN-GB">3.6</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Interactive mode and
    180180debugging <br>
    181 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_3.7.html"><span style="" lang="EN-GB">3.7</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Optional
     181</span><span style="font-family: Thorndale;"></span><span style="font-family: Thorndale;" lang="EN-GB"><a href="chapter_3.7.html">3.7</a> Optional
    182182software packages<o:p></o:p></span></p>
    183 <p style="margin: 12pt 0cm 0.0001pt 36pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.0.html"><span style="" lang="EN-GB">4.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Steering
     183<p style="margin: 12pt 0cm 0.0001pt 36pt;"><span style="font-family: Thorndale;"><a href="chapter_4.0.html"><span style="" lang="EN-GB">4.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Steering
    184184parameters <o:p></o:p></span></p>
    185 <p style="margin: 0cm 0cm 0.0001pt 72pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html"><span style="" lang="EN-GB">4.1</span></a></span><span style="font-family: Thorndale;" lang="EN-GB">
     185<p style="margin: 0cm 0cm 0.0001pt 72pt;"><span style="font-family: Thorndale;"><a href="chapter_4.1.html"><span style="" lang="EN-GB">4.1</span></a></span><span style="font-family: Thorndale;" lang="EN-GB">
    186186Initialization parameters <br>
    187 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.2.html"><span style="" lang="EN-GB">4.2</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Runtime
     187</span><span style="font-family: Thorndale;"><a href="chapter_4.2.html"><span style="" lang="EN-GB">4.2</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Runtime
    188188parameters and package parameters <br>
    189 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.3.html"><span style="" lang="EN-GB">4.3</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> User-defined
     189</span><span style="font-family: Thorndale;"><a href="chapter_4.3.html"><span style="" lang="EN-GB">4.3</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> User-defined
    190190parameters <br>
    191 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.4.html"><span style="" lang="EN-GB">4.4</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Example of a
     191</span><span style="font-family: Thorndale;"><a href="chapter_4.4.html"><span style="" lang="EN-GB">4.4</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Example of a
    192192minimum parameter set <br>
    193 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.5.html"><span style="" lang="EN-GB">4.5</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Data analysis and
     193</span><span style="font-family: Thorndale;"><a href="chapter_4.5.html"><span style="" lang="EN-GB">4.5</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Data analysis and
    194194visualization <o:p></o:p></span></p>
    195195
    196 <p style="margin: 0cm 0cm 0.0001pt 108pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.5.1.html"><span style="" lang="EN-GB">4.5.1</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> NetCDF data output<b></b><br>
     196<p style="margin: 0cm 0cm 0.0001pt 108pt;"><span style="font-family: Thorndale;"><a href="chapter_4.5.1.html"><span style="" lang="EN-GB">4.5.1</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> NetCDF data output<b></b><br>
    197197<a href="chapter_4.5.2.html">4.5.2</a> </span><span style="font-family: Thorndale;" lang="EN-GB">Plots of
    198198vertical profiles with <b>profil</b></span><br>
    199 <span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.5.3.html"><span style="" lang="EN-GB">4.5.3</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of time
     199<span style="font-family: Thorndale;"><a href="chapter_4.5.3.html"><span style="" lang="EN-GB">4.5.3</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of time
    200200series with <b>profil</b> <br>
    201 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.5.4.html"><span style="" lang="EN-GB">4.5.4</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of 2D
     201</span><span style="font-family: Thorndale;"><a href="chapter_4.5.4.html"><span style="" lang="EN-GB">4.5.4</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of 2D
    202202cross sections with <b>iso2d</b> <br>
    203 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.5.5.html"><span style="" lang="EN-GB">4.5.5</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of 3D
     203</span><span style="font-family: Thorndale;"><a href="chapter_4.5.5.html"><span style="" lang="EN-GB">4.5.5</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of 3D
    204204volume data with <b>AVS</b> <br>
    205 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.5.6.html"><span style="" lang="EN-GB">4.5.6</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of
     205</span><span style="font-family: Thorndale;"><a href="chapter_4.5.6.html"><span style="" lang="EN-GB">4.5.6</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of
    206206compressed 3D volume data with <b>AVS</b><br>
    207 </span><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.5.7.html"><span style="" lang="EN-GB">4.5.7</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of <span class="SpellE">isosurfaces</span>, 2D cross sections and particles
     207</span><span style="font-family: Thorndale;"><a href="chapter_4.5.7.html"><span style="" lang="EN-GB">4.5.7</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Plots of <span class="SpellE">isosurfaces</span>, 2D cross sections and particles
    208208with <span class="SpellE"><b>dvrp</b></span><o:p></o:p></span></p>
    209 <p style="margin: 0cm 0cm 5pt 72pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.6.html"><span style="" lang="EN-GB">4.6</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Listing of the steering
     209<p style="margin: 0cm 0cm 5pt 72pt;"><span style="font-family: Thorndale;"><a href="chapter_4.6.html"><span style="" lang="EN-GB">4.6</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Listing of the steering
    210210parameters in alphabetical&nbsp; order<o:p></o:p></span></p>
    211 <p style="margin: 12pt 0cm 0.0001pt 36pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_5.0.html"><span style="" lang="EN-GB">5.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Installation
     211<p style="margin: 12pt 0cm 0.0001pt 36pt;"><span style="font-family: Thorndale;"><a href="chapter_5.0.html"><span style="" lang="EN-GB">5.0</span></a></span><span style="font-family: Thorndale;" lang="EN-GB"> Installation
    212212of the model <o:p></o:p></span></p>
    213 <p style="margin: 0cm 0cm 5.1pt 72pt;"><span style="font-family: Thorndale;"><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_5.1.html">5.1</a>
     213<p style="margin: 0cm 0cm 5.1pt 72pt;"><span style="font-family: Thorndale;"><a href="chapter_5.1.html">5.1</a>
    214214Installation of <span class="SpellE">new</span><span class="SpellE"></span> / other
    215215<span class="SpellE">versions, version update</span><o:p></o:p></span></p><font color="#000000">&nbsp;
  • palm/trunk/SOURCE/check_parameters.f90

    r87 r89  
    14051405          CASE ( 'Sw' )
    14061406             dopr_index(i) = 39
     1407             dopr_unit(i)  = 'none'
    14071408             hom(:,2,39,:) = SPREAD( zw, 2, statistic_regions+1 )
    14081409
  • palm/trunk/SOURCE/modules.f90

    r87 r89  
    55! Actual revisions:
    66! -----------------
    7 ! +data_output_pr_user, max_pr_user, num_user_profiles, size of data_output_pr,
    8 ! dopr_index, dopr_initial_index and dopr_unit enlarged,
     7! +data_output_pr_user, max_pr_user, size of data_output_pr, dopr_index,
     8! dopr_initial_index and dopr_unit enlarged,
    99! var_hom and var_sum renamed pr_palm
    1010!
     
    275275                mg_switch_to_pe0_level = 0, ngsrb = 2, nsor = 20, &
    276276                nsor_ini = 100, n_sor, normalizing_region = 0, &
    277                 num_user_profiles = 0, nz_do1d, nz_do3d = -9999, &
    278                 outflow_damping_width = -1, prt_time_count = 0, runnr = 0, &
    279                 skip_do_avs = 0, timestep_count = 0
     277                nz_do1d, nz_do3d = -9999, outflow_damping_width = -1, &
     278                prt_time_count = 0, runnr = 0, skip_do_avs = 0, &
     279                timestep_count = 0
    280280
    281281    INTEGER ::  dist_nxl(0:1), dist_nxr(0:1), dist_nyn(0:1), dist_nys(0:1), &
  • palm/trunk/SOURCE/user_interface.f90

    r87 r89  
    7070
    7171
    72     NAMELIST /userpar/  data_output_pr_user, data_output_user, max_pr_user, &
    73                         region
     72    NAMELIST /userpar/  data_output_pr_user, data_output_user, region
    7473
    7574!
     
    8988    READ ( 11, userpar )
    9089    user_defined_namelist_found = .TRUE.
    91 
    92 !
    93 !-- Determine the maximum number of user-profiles allowed to be calculated.
    94 !-- This number may be (much) larger than num_user_profiles which only gives
    95 !-- the number of user-defined profiles to be output.
    96     IF ( data_output_pr_user(1) /= ' '  .AND.  max_pr_user == 0 )  THEN
    97        max_pr_user = 100
    98     ENDIF
    9990
    10091!
     
    109100       DO  WHILE ( data_output_pr_user(j) /= ' '  .AND.  j <= 100 )
    110101          data_output_pr(i) = data_output_pr_user(j)
    111           num_user_profiles = num_user_profiles + 1
     102          max_pr_user       = max_pr_user + 1
    112103          i = i + 1
    113104          j = j + 1
    114105       ENDDO
    115     ENDIF
    116 
    117     IF ( num_user_profiles > max_pr_user )  THEN
    118        IF ( myid == 0 )  THEN
    119           PRINT*, '+++ user_parin: number of user-defined profiles given by '
    120           PRINT*, '                data_output_pr_user (', num_user_profiles, &
    121                                    ') is larger than the allowed maximum'
    122           PRINT*, '                number of profiles: max_pr_user(', &
    123                                    max_pr_user, ')'
    124        ENDIF
    125        CALL local_stop
    126106    ENDIF
    127107
     
    529509!             DO  k = nzb_s_outer(j,i)+1, nzt
    530510!!
    531 !!--             Sample on how to calculate the profile for vertical velocity
    532 !                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +       &
    533 !                                                 w(k,j,i)               &
     511!!--             Sample on how to calculate the profile of the resolved-scale
     512!!--             horizontal momentum flux u*v*
     513!                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +           &
     514!                      ( 0.5 * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) * &
     515!                      ( 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) ) * &
    534516!                                                 * rmask(j,i,sr)
    535517!!
     
    770752!--    the number "pr_palm+1", the second one "pr_palm+2", etc. The respective
    771753!--    user-profile-numbers have also to be used in routine user_statistics!
    772 !       CASE ( 'w_user' )                  ! quantity string as given in
    773 !                                          ! data_output_pr_user
     754!       CASE ( 'u*v*' )                      ! quantity string as given in
     755!                                            ! data_output_pr_user
    774756!          index = pr_palm + 1
    775 !          dopr_index(var_count)  = index  ! the quantities' user-profile-number
    776 !          dopr_unit(var_count)   = 'm/s'  ! quantity unit
    777 !          hom(:,2,index,:)       = SPREAD( zw, 2, statistic_regions+1 )
    778 !                                          ! grid on which the quantity is
    779 !                                          ! defined (use zu or zw)
     757!          dopr_index(var_count)  = index    ! quantities' user-profile-number
     758!          dopr_unit(var_count)   = 'm2/s2'  ! quantity unit
     759!          hom(:,2,index,:)       = SPREAD( zu, 2, statistic_regions+1 )
     760!                                            ! grid on which the quantity is
     761!                                            ! defined (use zu or zw)
    780762
    781763       CASE DEFAULT
Note: See TracChangeset for help on using the changeset viewer.