Changes between Version 389 and Version 390 of doc/app/initialization_parameters

Timestamp:

Jun 12, 2018 11:18:03 AM (7 years ago)

Author:

Giersch

Comment:

--

doc/app/initialization_parameters

@@ -659 +659 @@
 }}}
 {{{#!td style="vertical-align:top"
-R
-}}}
-{{{#!td style="vertical-align:top"
-
+R(10)
+}}}
+{{{#!td style="vertical-align:top"
+10 * -1.0
 }}}
 {{{#!td
 Vertical grid spacing (in m).\\\\
-This parameter must be assigned by the user because no default value is given.\\\\
-By default, the model uses constant grid spacing along the z-direction, but it can be stretched using the parameters [#dz_stretch_level dz_stretch_level] and [#dz_stretch_factor dz_stretch_factor]. In case of stretching, a maximum allowed grid spacing can be given by [#dz_max dz_max].\\\\
-Assuming a constant '''dz''', the scalar levels (zu) are calculated directly by:\\\\
-      zu(0) = - '''dz''' * 0.5\\\\ 
-      zu(1) = '''dz''' * 0.5\\\\
+At least one value must be assigned by the user because no realistic value is given.\\\\
+By default, the model uses constant grid spacing along the z-direction, but it can be stretched using either the old parameters [#dz_stretch_level dz_stretch_level], [#dz_stretch_factor dz_stretch_factor], [#dz_max dz_max] or the new ones [#dz_stretch_level_start dz_stretch_level_start], [#dz_stretch_level_end dz_stretch_level_end] which allow grid stretching to a finer or coarser grid and the definition of several stretching regions. If the new stretching method is used different dz values have to be specified by the user (equal to the number of defined [#dz_stretch_level_end dz_stretch_level_end] + 1). In that way PALM knows which grid spacing shall be used before and after the grid stretching region limited through [#dz_stretch_level_start dz_stretch_level_start] and [#dz_stretch_level_end dz_stretch_level_end].\\\\
+Assuming a constant dz, the scalar levels (zu) are calculated directly by:\\\\
+      zu(0) = - dz(1) * 0.5    or    zu(0) = 0.0 depending on the boundary condition you choose,\\\\ 
+      zu(1) = dz(1) * 0.5,\\\\
+      zu(k) = zu(k-1) + dz(1).\\\\
 The w-levels lie half between them:\\\\
-      zw(k) = ( zu(k) + zu(k+1) ) * 0.5
+      zw(k) = ( zu(k) + zu(k+1) ) * 0.5.
 }}}
 |----------------
@@ -682 +683 @@
 }}}
 {{{#!td style="vertical-align:top"
-9999999.9
+1000.0
 }}}
 {{{#!td
 Allowed maximum vertical grid spacing (in m).\\\\
-If the vertical grid is stretched (see [#dz_stretch_factor dz_stretch_factor] and [#dz_stretch_level dz_stretch_level]), '''dz_max''' can be used to limit the vertical grid spacing.\\\\
-'''Attention:'''\\
-So far, vertical grid stretching cannot be applied when simultaneously using PALM's particle model (see list of [../parpar Particle Parameters]). 
+If the vertical grid is stretched (see [#dz_stretch_factor dz_stretch_factor] and [#dz_stretch_level dz_stretch_level]), '''dz_max''' can be used to limit the vertical grid spacing. The default value is 1000 to prevent unrealistic high vertical grid spacings.
 }}}
 |----------------
@@ -712 +711 @@
 }}}
 {{{#!td style="vertical-align:top"
-100000.0
+-9999999.9
 }}}
 {{{#!td
@@ -721 +720 @@
       zw(k) = ( zu(k) + zu(k+1) ) * 0.5.\\\\
 For [#ocean ocean] = ''.T.,'' '''dz_stretch_level''' is the height level (in m, negative) below which the grid is to be stretched vertically. The vertical grid spacings [#dz dz] below this level are calculated correspondingly as\\\\
-      dz(k-1) = dz(k) * [#dz_stretch_factor dz_stretch_factor].\\\\
-'''Attention:'''\\
-So far, vertical grid stretching cannot be applied when simultaneously using PALM's particle model (see list of [../parpar Particle Parameters]). 
+      dz(k-1) = dz(k) * [#dz_stretch_factor dz_stretch_factor].
+}}}
+|----------------
+{{{#!td style="vertical-align:top"
+[=#dz_stretch_level_end '''dz_stretch_level_end''']
+}}}
+{{{#!td style="vertical-align:top"
+R(9)
+}}}
+{{{#!td style="vertical-align:top"
+9 * 9999999.9
+}}}
+{{{#!td
+Height level until which the grid is to be stretched vertically (in m). Up to 9 heights/separate stretching regions are possible.\\\\
+For [#ocean ocean] = ''.F.,'' '''dz_stretch_level_end''' is the height level (in m) until which the grid is to be stretched vertically. The vertical grid spacings [#dz dz] between this level and the corresponding [#dz_stretch_level_start dz_stretch_level_start] are calculated as\\\\
+      dz(k+1) = dz(k) * dz_stretch_factor_array\\\\
+and used as spacings for the scalar levels (zu). The w-levels are then defined as:\\\\
+      zw(k) = ( zu(k) + zu(k+1) ) * 0.5.\\\\
+For [#ocean ocean] = ''.T.,'' '''dz_stretch_level_end''' is also the height level (in m) until which the grid is to be stretched vertically but it is defined negative. The vertical grid spacings [#dz dz] between this level and the corresponding [#dz_stretch_level_start dz_stretch_level_start] are calculated as\\\\
+      dz(k-1) = dz(k) * dz_stretch_factor_array.\\\\
+For each '''dz_stretch_level_end''' a corresponding [=#dz_stretch_level_start dz_stretch_level_start] must be defined. \\\\
+'''Hint:'''\\
+dz_stretch_factor_array is an array in the code which accounts for different stretching factors due to different stretching regions that are possible. The factors are calculated internally except for the case when the old stretching for the free atmosphere with [#dz_stretch_factor dz_stretch_factor], [#dz_stretch_level dz_stretch_level] and [#dz_max dz_max] is used. 
+}}}
+|----------------
+{{{#!td style="vertical-align:top"
+[=#dz_stretch_level_start '''dz_stretch_level_start''']
+}}}
+{{{#!td style="vertical-align:top"
+R(9)
+}}}
+{{{#!td style="vertical-align:top"
+9 * -9999999.9
+}}}
+{{{#!td
+Height level above/below which the grid is to be stretched vertically (in m). Up to 9 heights/separate stretching regions are possible.\\\\
+For [#ocean ocean] = ''.F.,'' '''dz_stretch_level_start''' is the height level (in m) above which the grid is to be stretched vertically. The vertical grid spacings [#dz dz] between this level and the corresponding [#dz_stretch_level_end dz_stretch_level_end] are calculated as\\\\
+      dz(k+1) = dz(k) * dz_stretch_factor_array\\\\
+and used as spacings for the scalar levels (zu). The w-levels are then defined as:\\\\
+      zw(k) = ( zu(k) + zu(k+1) ) * 0.5.\\\\
+For [#ocean ocean] = ''.T.,'' '''dz_stretch_level_start''' is the height level (in m, negative) below which the grid is to be stretched vertically. The vertical grid spacings [#dz dz] between this level and the corresponding [#dz_stretch_level_end dz_stretch_level_end] are calculated as\\\\
+      dz(k-1) = dz(k) * dz_stretch_factor_array.\\\\
+For each '''dz_stretch_level_start''' a corresponding [#dz_stretch_level_end dz_stretch_level_end] must be defined except for the last level. Here, it is possible to omit the value for [#dz_stretch_level_end dz_stretch_level_end] to consider 'endless' stretching until the value of [=dz_max dz_max]) is reached. In that case the stretching factor can not be calculated and is set to the value of [=dz_stretch_factor '''dz_stretch_factor''']. \\\\
+'''Hint:'''\\
+dz_stretch_factor_array is an array in the code which accounts for different stretching factors due to different stretching regions that can be defined. The factors are calculated internally except for the case when the old stretching for the free atmosphere with [#dz_stretch_factor dz_stretch_factor], [#dz_stretch_level dz_stretch_level] and [#dz_max dz_max] is used. 
 }}}
 |----------------