Changeset 75 for palm/trunk/DOC


Ignore:
Timestamp:
Mar 22, 2007 9:54:05 AM (17 years ago)
Author:
raasch
Message:

preliminary update for changes concerning non-cyclic boundary conditions

Location:
palm/trunk/DOC/app
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/DOC/app/chapter_4.1.html

    r72 r75  
    8888<span style="font-weight: bold;">bc_lr</span> may
    8989also be
    90 assigned the values <span style="font-style: italic;">'dirichlet/neumann'</span>
    91 (inflow from left, outflow to the right) or <span style="font-style: italic;">'neumann/dirichlet'</span>
     90assigned the values <span style="font-style: italic;">'dirichlet/radiation'</span>
     91(inflow from left, outflow to the right) or <span style="font-style: italic;">'radiation/dirichlet'</span>
    9292(inflow from
    9393right, outflow to the left). This requires the multi-grid method to be
    9494used for solving the Poisson equation for perturbation pressure (see <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.2.html#psolver">psolver</a>)
    95 and it also requires cyclic boundary conditions along y (see<br> <a href="#bc_ns">bc_ns</a>).<br> <br>
     95and it also requires cyclic boundary conditions along y (see&nbsp;<a href="#bc_ns">bc_ns</a>).<br> <br>
    9696In case of these non-cyclic lateral boundaries, a Dirichlet condition
    9797is used at the inflow for all quantities (initial vertical profiles -
    9898see <a href="#initializing_actions">initializing_actions</a>
    9999- are fixed during the run) except u, to which a Neumann (zero
    100 gradient) condition is applied. At the outflow, a Neumann (zero
    101 gradient) condition is used for all quantities except v, which is set
    102 to its horizontal average along the outflow (e.g. v(k,:,nx+1) =
    103 average_along_y( v(k,:,nx)), and except w, which is set to zero
    104 (Dirichlet condition). These conditions ensure the velocity field to be
    105 free of divergence at the inflow and at the outflow. For perturbation
     100gradient) condition is applied. At the outflow, a radiation condition is used for all velocity components, while a Neumann (zero
     101gradient) condition is used for the scalars. For perturbation
    106102pressure Neumann (zero gradient) conditions are assumed both at the
    107103inflow and at the outflow.<br> <br>
     
    137133<span style="font-weight: bold;">bc_ns</span> may
    138134also be
    139 assigned the values <span style="font-style: italic;">'dirichlet/neumann'</span>
    140 (inflow from rear ("north"), outflow to the front ("south")) or <span style="font-style: italic;">'neumann/dirichlet'</span>
     135assigned the values <span style="font-style: italic;">'dirichlet/radiation'</span>
     136(inflow from rear ("north"), outflow to the front ("south")) or <span style="font-style: italic;">'radiation/dirichlet'</span>
    141137(inflow from front ("south"), outflow to the rear ("north")). This
    142138requires the multi-grid
     
    146142In case of these non-cyclic lateral boundaries, a Dirichlet condition
    147143is used at the inflow for all quantities (initial vertical profiles -
    148 see <a href="#initializing_actions">initializing_actions</a>
    149 - are fixed during the run) except v, to which a Neumann (zero
    150 gradient) condition is applied. At the outflow, a Neumann (zero
    151 gradient) condition is used for all quantities except u, which is set
    152 to its horizontal average along the outflow (e.g. u(k,ny+1,:) =
    153 average_along_x( u(k,ny,:)), and except w, which is set to zero
    154 (Dirichlet condition). These conditions ensure the velocity field to be
    155 free of divergence at the inflow and at the outflow. For perturbation
     144see <a href="chapter_4.1.html#initializing_actions">initializing_actions</a>
     145- are fixed during the run) except u, to which a Neumann (zero
     146gradient) condition is applied. At the outflow, a radiation condition is used for all velocity components, while a Neumann (zero
     147gradient) condition is used for the scalars. For perturbation
    156148pressure Neumann (zero gradient) conditions are assumed both at the
    157149inflow and at the outflow.<br> <br>
     
    428420Also cloud-top cooling by longwave radiation can be utilized (see <a href="#radiation">radiation</a>)<br> <b><br>
    429421cloud_physics =</b> <span style="font-style: italic;">.TRUE.
    430 </span>requires <a href="#moisture">moisture</a>
     422</span>requires&nbsp;<a href="#humidity">humidity</a>
    431423=<span style="font-style: italic;"> .TRUE.</span> .<br>
    432424Detailed information about the condensation scheme is given in the
     
    890882the so-called Blackadar mixing length is used (l = kappa * z / ( 1 +
    891883kappa * z / lambda ) with the limiting value lambda = 2.7E-4 * u_g / f).<br>
    892 </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="moisture"></a><b>moisture</b></p>
     884</td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="humidity"></a><b>humidity</b></p>
    893885</td> <td style="vertical-align: top;">L</td>
    894886<td style="vertical-align: top;"><i>.F.</i></td>
     
    953945<br> <span style="font-weight: bold;">Important: </span>The&nbsp;
    954946upstream-spline scheme is not implemented for humidity and passive
    955 scalars (see <a href="#moisture">moisture</a>
     947scalars (see&nbsp;<a href="#humidity">humidity</a>
    956948and <a href="#passive_scalar">passive_scalar</a>)
    957949and requires the use of a 2d-domain-decomposition. The last conditions
     
    12061198With <span style="font-weight: bold;">passive_scalar</span>
    12071199switched
    1208 on, the simultaneous use of humidity (see <a href="#moisture">moisture</a>)
     1200on, the simultaneous use of humidity (see&nbsp;<a href="#humidity">humidity</a>)
    12091201is impossible.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="phi"></a><b>phi</b></p>
    12101202</td> <td style="vertical-align: top;">R</td>
     
    13671359<br> <b>_level</b></p> </td> <td style="vertical-align: top;">R (10)</td> <td style="vertical-align: top;"> <p><i>10 *</i>&nbsp;
    13681360<i>0.0</i></p> </td> <td style="vertical-align: top;"> <p>Height level from
    1369 which on the moisture gradient defined by <a href="#q_vertical_gradient">q_vertical_gradient</a>
     1361which on the humidity gradient defined by <a href="#q_vertical_gradient">q_vertical_gradient</a>
    13701362is effective (in m).&nbsp; </p> <p>The height levels
    13711363are to be assigned in ascending order. The
     
    15331525</p> <p style="margin-left: 40px;"><span style="font-weight: bold;">Important: </span>The&nbsp;
    15341526upstream-spline scheme is not implemented for humidity and passive
    1535 scalars (see <a href="#moisture">moisture</a>
     1527scalars (see&nbsp;<a href="#humidity">humidity</a>
    15361528and <a href="#passive_scalar">passive_scalar</a>)
    15371529and requires the use of a 2d-domain-decomposition. The last conditions
     
    18231815<i>&nbsp;</i><a href="#alpha_surface">alpha_surface</a>
    18241816= 0.0, <a href="#bc_lr">bc_lr</a> = <a href="#bc_ns">bc_ns</a> = <span style="font-style: italic;">'cyclic'</span>,&nbsp;<a style="" href="#galilei_transformation">galilei_transformation</a>
    1825 = <span style="font-style: italic;">.F.</span>,&nbsp;<a href="#cloud_physics">cloud_physics&nbsp;</a> = <span style="font-style: italic;">.F.</span>,&nbsp; <a href="#cloud_droplets">cloud_droplets</a> = <span style="font-style: italic;">.F.</span>,&nbsp; <a href="#moisture">moisture</a> = <span style="font-style: italic;">.F.</span>, and <a href="#prandtl_layer">prandtl_layer</a> = .T..<br>
     1817= <span style="font-style: italic;">.F.</span>,&nbsp;<a href="#cloud_physics">cloud_physics&nbsp;</a> = <span style="font-style: italic;">.F.</span>,&nbsp; <a href="#cloud_droplets">cloud_droplets</a> = <span style="font-style: italic;">.F.</span>,&nbsp;&nbsp;<a href="#humidity">humidity</a> = <span style="font-style: italic;">.F.</span>, and <a href="#prandtl_layer">prandtl_layer</a> = .T..<br>
    18261818<font color="#000000"><br>
    18271819Note that an inclined model domain requires the use of <span style="font-weight: bold;">topography</span> = <span style="font-style: italic;">'flat'</span> and a
  • palm/trunk/DOC/app/chapter_4.2.html

    r72 r75  
    394394= <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">pt</span></td><td style="width: 196px; vertical-align: top;">potential
    395395temperature<br></td><td style="vertical-align: top;">K</td><td style="vertical-align: top;"></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">q</span></td><td style="width: 196px; vertical-align: top;">specific humidity
    396 (or total water content, if cloud physics is switched on)</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#moisture">moisture</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">ql</span></td><td style="width: 196px; vertical-align: top;">liquid water
     396(or total water content, if cloud physics is switched on)</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires&nbsp;<a href="chapter_4.1.html#humidity">humidity</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">ql</span></td><td style="width: 196px; vertical-align: top;">liquid water
    397397content</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#cloud_physics">cloud_physics</a>
    398398= <span style="font-style: italic;">.TRUE.</span>
     
    414414is allowed</td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">v</span></td><td style="width: 196px; vertical-align: top;">v-component of
    415415the velocity</td><td style="vertical-align: top;">m/s</td><td style="vertical-align: top;"></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">vpt</span></td><td style="width: 196px; vertical-align: top;">virtual potential
    416 temperature</td><td style="vertical-align: top;">K</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#moisture">moisture</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">w</span></td><td style="width: 196px; vertical-align: top;">w-component of
     416temperature</td><td style="vertical-align: top;">K</td><td style="vertical-align: top;">requires&nbsp;<a href="chapter_4.1.html#humidity">humidity</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">w</span></td><td style="width: 196px; vertical-align: top;">w-component of
    417417the velocity</td><td style="vertical-align: top;">m/s</td><td style="vertical-align: top;"></td></tr><tr><td style="vertical-align: top;"><span style="font-style: italic;">z0*</span></td><td style="vertical-align: top;">roughness length</td><td style="vertical-align: top;">m</td><td></td></tr></tbody></table><br>Multiple
    418418quantities can be assigned, e.g. <span style="font-weight: bold;">data_output</span>
  • palm/trunk/DOC/app/chapter_4.6.html

    r72 r75  
    652652<td style="vertical-align: middle;" width="57%"> <p>Variable
    653653to adjust the subdomain sizes in parallel runs.</p> </td> </tr>
    654 <tr> <td style="vertical-align: middle;"><b><a href="chapter_4.1.html#inflow_disturbance_begin"><b>inflow_disturbance_begin</b></a></b></td>
     654<tr><td><b><a href="chapter_4.1.html#humidity"><b>humidity</b></a></b></td><td>I</td><td>L</td><td><span style="font-style: italic;">.F.</span></td><td>Parameter
     655to switch on the prognostic equation for
     656specific
     657humidity q.</td></tr><tr> <td style="vertical-align: middle;"><b><a href="chapter_4.1.html#inflow_disturbance_begin"><b>inflow_disturbance_begin</b></a></b></td>
    655658<td style="vertical-align: middle;">I<br> </td>
    656659<td style="vertical-align: middle;">I<br> </td>
     
    749752which are
    750753to be created by the <span style="font-weight: bold;">dvrp</span>
    751 software. <br> </p> </td> </tr> <tr>
    752 <td style="vertical-align: middle;" width="15%"> <p><a href="chapter_4.1.html#moisture"><b>moisture</b></a></p>
    753 </td> <td style="vertical-align: middle;" width="5%">
    754 <p>I</p> </td> <td style="vertical-align: middle;" width="7%"> <p>L</p>
    755 </td> <td style="vertical-align: middle;" width="16%">
    756 <p><i>.F.</i></p> </td> <td style="vertical-align: middle;" width="57%"> <p>Parameter
    757 to switch on the prognostic equation for
    758 specific
    759 humidity q.</p> </td> </tr> <tr> <td style="vertical-align: middle;" width="15%"> <p><a href="chapter_4.1.html#momentum_advec"><b>momentum_advec</b></a></p>
     754software. <br> </p> </td> </tr>  <tr> <td style="vertical-align: middle;" width="15%"> <p><a href="chapter_4.1.html#momentum_advec"><b>momentum_advec</b></a></p>
    760755</td> <td style="vertical-align: middle;" width="5%">
    761756<p>I</p> </td> <td style="vertical-align: middle;" width="7%"> <p>C
     
    11051100(10)</p> </td> <td style="vertical-align: middle;" width="16%"> <p><i>10 *</i>&nbsp; <i>0.0</i></p>
    11061101</td> <td style="vertical-align: middle;" width="57%">
    1107 <p>Height level from which on the moisture gradient defined
     1102<p>Height level from which on the humidity gradient defined
    11081103by <a href="chapter_4.1.html#q_vertical_gradient">q_vertical_gradient</a>
    11091104is effective (in m).&nbsp;</p> </td> </tr> <tr>
Note: See TracChangeset for help on using the changeset viewer.