Changes between Version 457 and Version 458 of doc/app/initialization_parameters

Timestamp:

Feb 12, 2020 12:59:19 PM (5 years ago)

Author:

scharf

Comment:

--

doc/app/initialization_parameters

@@ -10 +10 @@
 ==== [#others Others] ====
 \\\\
-'''NAMELIST group name:''' [=#inipar '''{{{initialization_parameters}}}''']
+'''NAMELIST group name:''' [=#initialization_parameters '''{{{initialization_parameters}}}''']
 ----
 
@@ -534 +534 @@
 If it is set to .F., the large-scale subsidence velocity w_subs is used to advect the prognostic variables (subroutine  {{{subsidence}}}).
 If it is set to .T., the subroutine {{{subsidence}}} is not called and the subsidence tendencies from [../../app/iofiles#LSF_DATA LSF_DATA] are applied to the prognostic variables in the new subroutine {{{ls_advec}}}.
-The usage of use_subsidence_tendencies requires [wiki:doc/app/inipar#large_scale_forcing large_scale_forcing] = .T. as well as [wiki:doc/app/inipar#large_scale_subsidence large_scale_subsidence] = .T..
+The usage of use_subsidence_tendencies requires [#large_scale_forcing large_scale_forcing] = .T. as well as [#large_scale_subsidence large_scale_subsidence] = .T..
 }}}
 |----------------
@@ -875 +875 @@
 '''Remark 1''' \\
 
-In case of very complex topography, including narrow street canyons or complex shapes resolved by only a few grid points, the use of [../../app/inipar#psolver psolver] = 'multigrid_noopt' in combination with [../../app/inipar#masking_method masking_method] = .F. can lead to a blow up of wind velocity. So far, this never happened if masking_method = .T. or [../../app/inipar#psolver psolver] = 'multigrid' (where topography is always masked). 
+In case of very complex topography, including narrow street canyons or complex shapes resolved by only a few grid points, the use of [#psolver psolver] = 'multigrid_noopt' in combination with [#masking_method masking_method] = .F. can lead to a blow up of wind velocity. So far, this never happened if masking_method = .T. or [#psolver psolver] = 'multigrid' (where topography is always masked). 
 A possible workaround to prevent such velocity blow-up is to preprocess the topography by filling holes and removing complex shapes on the small scale before running the LES. (This approach should have no significant effect on the flow field, as in LES the flow in such regions is poorly resolved and do not yield to reliable physical information.)
 The reason for this velocity blow-up is unclear so far. 
@@ -882 +882 @@
 '''Remark 2''' \\
 
-If topography is prescribed and [../../app/inipar#psolver psolver] = 'multigrid_noopt', the mean velocity-divergence reduction (see RUN_CONTROL) is smaller in case of [../../app/inipar#masking_method masking_method] = .T. compared to [../../app/inipar#masking_method masking_method] = .F. (about one-order of magnitude).
+If topography is prescribed and [#psolver psolver] = 'multigrid_noopt', the mean velocity-divergence reduction (see RUN_CONTROL) is smaller in case of [#masking_method masking_method] = .T. compared to [#masking_method masking_method] = .F. (about one-order of magnitude).
 This smaller reduction can be attributed to wall-bounded grid points, where the divergence after pressure correction remains significantly larger in case of masking, while the divergence reduction for non-wall-bounded grid points is similar in both cases. 
 However, validation tests revealed that this smaller divergence reduction has no effect on the resulting flow.  
@@ -2064 +2064 @@
 '''In case of cyclic boundary conditions:'''\\
 Shifts the left/right boundary by multiples of a subdomain size along y.\\
-Requires  [../inipar#bc_lr bc_lr] = 'cyclic',  [../inipar#bc_ns bc_ns] = 'cyclic' and [../inipar#psolver psolver] = 'multigrid'.\\\\
+Requires  [#bc_lr bc_lr] = 'cyclic',  [#bc_ns bc_ns] = 'cyclic' and [#psolver psolver] = 'multigrid'.\\\\
 This parameter introduces a shift by multiples of a subdomain size in y-direction at the left/right domain boundary. The shift is such that a point (coordinate y=yr) on the right domain boundary no longer corresponds with the point with the same y-coordinate on the left boundary, but instead with the point with the coordinate {{{y=yr+y_shift*(ny+1)/npey}}}. Conversely, a point on the left boundary with coordinate {{{y=yl}}} corresponds with the point on the right boundary with the coordinate {{{y=yl-y_shift*(ny+1)/npey}}}.\\
 Negative values and also those that exceed npey are allowed. They are transformed to the range {{{[0,npey-1]}}} internally.\\\\
 This feature has been implemented to alleviate the occurrence of streak-like structures that appear during simulations with neutral stratification, cyclic boundary conditions, and driving velocity with no v-component. For a description of these structures, see Munters (2016; ​http://www.dx.doi.org/10.1063/1.4941912).\\\\
 '''Schematic of the y_shift method:'''\\
-(Model domain with [../inipar#npex npex] = [../inipar#npey npey] = 5 and y_shift = 1)\\
+(Model domain with [#npex npex] = [#npey npey] = 5 and y_shift = 1)\\
 [[Image(y_shift_pic.png,260px)]]\\\\
-'''In case of [../initialization_parameters#bc_lr bc_lr] = 'dirichlet/radiation' and [../initialization_parameters#turbulent_inflow turbulent_inflow] = .TRUE.:'''\\
+'''In case of [#bc_lr bc_lr] = 'dirichlet/radiation' and [#turbulent_inflow turbulent_inflow] = .TRUE.:'''\\
 Parameter to define a y-shift for the recycled inflow turbulence, given in multiples of PE.
 E.g. with y_shift = 3 the turbulent fluctuations, obtained at the recycling plane, will be shifted by 3 processors in positive y-direction before being imposed on the inflow. The parameter y_shift replaces the deprecated parameter recycling_yshift since revision r4301.
@@ -3106 +3106 @@
 By setting '''alpha_surface''' unequal to zero the model domain can be inclined in x-direction with respect to the horizontal. In this way flows over inclined surfaces (e.g. drainage flows, gravity flows) can be simulated. In case of '''alpha_surface''' /= ''0'' the buoyancy term appears both in the equation of motion of the u-component and of the w-component.\\\\
 An inclination is only possible in case of cyclic horizontal boundary conditions along x '''AND''' y (see [#bc_lr bc_lr] and [#bc_ns bc_ns]) and [#topography topography] = '' 'flat'.''\\\\
-In the case of [../inipar#alpha_surface alpha_surface] /= 0.0 the simultaneous use of [../inipar#humidity humidity] = ''.TRUE.'' is not allowed.\\\\
+In the case of [#alpha_surface alpha_surface] /= 0.0 the simultaneous use of [#humidity humidity] = ''.TRUE.'' is not allowed.\\\\
 Runs with inclined surface still require additional [../userint user-defined code] as well as modifications to the default code. Please ask the PALM developer group.
 }}}