Changeset 75 for palm/trunk/DOC/app/chapter_4.1.html

Timestamp:

Mar 22, 2007 9:54:05 AM (18 years ago)

Author:

raasch

Message:

preliminary update for changes concerning non-cyclic boundary conditions

File:

• palm/trunk/DOC/app/chapter_4.1.html (modified) (10 diffs)

palm/trunk/DOC/app/chapter_4.1.html

@@ -88 +88 @@
 <span style="font-weight: bold;">bc_lr</span> may
 also be
-assigned the values <span style="font-style: italic;">'dirichlet/neumann'</span>
-(inflow from left, outflow to the right) or <span style="font-style: italic;">'neumann/dirichlet'</span>
+assigned the values <span style="font-style: italic;">'dirichlet/radiation'</span>
+(inflow from left, outflow to the right) or <span style="font-style: italic;">'radiation/dirichlet'</span>
 (inflow from
 right, outflow to the left). This requires the multi-grid method to be
 used 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>)
-and it also requires cyclic boundary conditions along y (see<br> <a href="#bc_ns">bc_ns</a>).<br> <br>
+and it also requires cyclic boundary conditions along y (see&nbsp;<a href="#bc_ns">bc_ns</a>).<br> <br>
 In case of these non-cyclic lateral boundaries, a Dirichlet condition
 is used at the inflow for all quantities (initial vertical profiles -
 see <a href="#initializing_actions">initializing_actions</a>
 - are fixed during the run) except u, to which a Neumann (zero
-gradient) condition is applied. At the outflow, a Neumann (zero
-gradient) condition is used for all quantities except v, which is set
-to its horizontal average along the outflow (e.g. v(k,:,nx+1) =
-average_along_y( v(k,:,nx)), and except w, which is set to zero
-(Dirichlet condition). These conditions ensure the velocity field to be
-free of divergence at the inflow and at the outflow. For perturbation
+gradient) condition is applied. At the outflow, a radiation condition is used for all velocity components, while a Neumann (zero
+gradient) condition is used for the scalars. For perturbation
 pressure Neumann (zero gradient) conditions are assumed both at the
 inflow and at the outflow.<br> <br>
@@ -137 +133 @@
 <span style="font-weight: bold;">bc_ns</span> may
 also be
-assigned the values <span style="font-style: italic;">'dirichlet/neumann'</span>
-(inflow from rear ("north"), outflow to the front ("south")) or <span style="font-style: italic;">'neumann/dirichlet'</span>
+assigned the values <span style="font-style: italic;">'dirichlet/radiation'</span>
+(inflow from rear ("north"), outflow to the front ("south")) or <span style="font-style: italic;">'radiation/dirichlet'</span>
 (inflow from front ("south"), outflow to the rear ("north")). This
 requires the multi-grid
@@ -146 +142 @@
 In case of these non-cyclic lateral boundaries, a Dirichlet condition
 is used at the inflow for all quantities (initial vertical profiles -
-see <a href="#initializing_actions">initializing_actions</a>
-- are fixed during the run) except v, to which a Neumann (zero
-gradient) condition is applied. At the outflow, a Neumann (zero
-gradient) condition is used for all quantities except u, which is set
-to its horizontal average along the outflow (e.g. u(k,ny+1,:) =
-average_along_x( u(k,ny,:)), and except w, which is set to zero
-(Dirichlet condition). These conditions ensure the velocity field to be
-free of divergence at the inflow and at the outflow. For perturbation
+see <a href="chapter_4.1.html#initializing_actions">initializing_actions</a>
+- are fixed during the run) except u, to which a Neumann (zero
+gradient) condition is applied. At the outflow, a radiation condition is used for all velocity components, while a Neumann (zero
+gradient) condition is used for the scalars. For perturbation
 pressure Neumann (zero gradient) conditions are assumed both at the
 inflow and at the outflow.<br> <br>
@@ -428 +420 @@
 Also cloud-top cooling by longwave radiation can be utilized (see <a href="#radiation">radiation</a>)<br> <b><br>
 cloud_physics =</b> <span style="font-style: italic;">.TRUE.
-</span>requires <a href="#moisture">moisture</a>
+</span>requires&nbsp;<a href="#humidity">humidity</a>
 =<span style="font-style: italic;"> .TRUE.</span> .<br>
 Detailed information about the condensation scheme is given in the
@@ -890 +882 @@
 the so-called Blackadar mixing length is used (l = kappa * z / ( 1 +
 kappa * z / lambda ) with the limiting value lambda = 2.7E-4 * u_g / f).<br>
-</td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="moisture"></a><b>moisture</b></p>
+</td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="humidity"></a><b>humidity</b></p>
 </td> <td style="vertical-align: top;">L</td>
 <td style="vertical-align: top;"><i>.F.</i></td>
@@ -953 +945 @@
 <br> <span style="font-weight: bold;">Important: </span>The&nbsp;
 upstream-spline scheme is not implemented for humidity and passive
-scalars (see <a href="#moisture">moisture</a>
+scalars (see&nbsp;<a href="#humidity">humidity</a>
 and <a href="#passive_scalar">passive_scalar</a>)
 and requires the use of a 2d-domain-decomposition. The last conditions
@@ -1206 +1198 @@
 With <span style="font-weight: bold;">passive_scalar</span>
 switched
-on, the simultaneous use of humidity (see <a href="#moisture">moisture</a>)
+on, the simultaneous use of humidity (see&nbsp;<a href="#humidity">humidity</a>)
 is impossible.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="phi"></a><b>phi</b></p>
 </td> <td style="vertical-align: top;">R</td>
@@ -1367 +1359 @@
 <br> <b>_level</b></p> </td> <td style="vertical-align: top;">R (10)</td> <td style="vertical-align: top;"> <p><i>10 *</i>&nbsp;
 <i>0.0</i></p> </td> <td style="vertical-align: top;"> <p>Height level from
-which on the moisture gradient defined by <a href="#q_vertical_gradient">q_vertical_gradient</a>
+which on the humidity gradient defined by <a href="#q_vertical_gradient">q_vertical_gradient</a>
 is effective (in m).&nbsp; </p> <p>The height levels
 are to be assigned in ascending order. The
@@ -1533 +1525 @@
 </p> <p style="margin-left: 40px;"><span style="font-weight: bold;">Important: </span>The&nbsp;
 upstream-spline scheme is not implemented for humidity and passive
-scalars (see <a href="#moisture">moisture</a>
+scalars (see&nbsp;<a href="#humidity">humidity</a>
 and <a href="#passive_scalar">passive_scalar</a>)
 and requires the use of a 2d-domain-decomposition. The last conditions
@@ -1823 +1815 @@
 <i>&nbsp;</i><a href="#alpha_surface">alpha_surface</a>
 = 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>
-= <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>
+= <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>
 <font color="#000000"><br>
 Note 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