== Initialization parameters == ==== [#mode Mode] ==== ==== [#grid Grid] ==== ==== [#num Numerics] ==== ==== [#phys Physics] ==== ==== [#bc Boundary Conditions] ==== ==== [#ini Initialization] ==== ==== [#topo Topography] ==== ==== [#canopy Canopy] ==== ==== [#others Others] ==== \\\\ NAMELIST group name: '''inipar'''\\ ---- [=#mode '''Mode:]\\ ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN]\\[../fortrantypes Type]''' =||='''Default\\Value''' =||='''Explanation''' =|| |---------------- {{{#!td style="vertical-align:top;width: 150px" [=#cloud_droplets '''cloud_droplets'''] }}} {{{#!td style="vertical-align:top;width: 50px" L }}} {{{#!td style="vertical-align:top;width: 75px" .F. }}} {{{#!td Parameter to switch on usage of cloud droplets.\\\\ Cloud droplets require to use particles (i.e. the NAMELIST group [../parpar particles_par] has to be included in the parameter file). Then each particle is a representative for a certain number of droplets. The droplet features (number of droplets, initial radius, etc.) can be steered with the respective particle parameters (see e.g. [../parpar#radius radius]). The real number of initial droplets in a grid cell is equal to the initial number of droplets (defined by the particle source parameters [../parpar#pst pst], [../parpar#psl psl], [../parpar#psr psr], [../parpar#pss pss], [../parpar#psn psn], [../parpar#psb psb], [../parpar#pdx pdx], [../parpar#pdy pdy] and [../parpar#pdz pdz]) times the [../parpar#initial_weighting_factor initial_weighting_factor].\\\\ In case of using cloud droplets, the default condensation scheme in PALM cannot be used, i.e. [#cloud_physics cloud_physics] must be set ''.F.''. }}} |---------------- {{{#!td style="vertical-align:top" [=#cloud_physics '''cloud_physics'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .F. }}} {{{#!td Parameter to switch on the condensation scheme.\\\\ For '''cloud_physics''' = ''.T.'', equations for the liquid water content and the liquid water potential temperature are solved instead of those for specific humidity and potential temperature. Note that a grid volume is assumed to be either completely saturated or completely unsaturated (0%-or-100%-scheme). A simple precipitation scheme can additionally be switched on with parameter [#precipitation precipitation]. Also cloud-top cooling by longwave radiation can be utilized (see [#radiation radiation]).\\\\ '''cloud_physics''' = ''.T.'' requires [#humidity humidity] = ''.T.''.\\\\ Detailed information about the condensation scheme is given in the description of the [[cloud physics module]] (pdf-file, only in German).\\\\ This condensation scheme is not allowed if cloud droplets are simulated explicitly (see [#cloud_droplets cloud_droplets]). }}} |---------------- {{{#!td style="vertical-align:top" [=#conserve_volume_flow '''conserve_volume_flow'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .F. }}} {{{#!td Conservation of volume flow in x- and y-direction.\\\\ '''conserve_volume_flow''' = ''.T.'' guarantees that the volume flow through the xz- and yz-cross-sections of the total model domain remains constant throughout the run depending on the chosen [#conserve_volume_flow_mode conserve_volume_flow_mode].\\\\ Note that '''conserve_volume_flow''' = ''.T.'' requires [#dp_external dp_external] = ''.F.''. }}} |---------------- {{{#!td style="vertical-align:top" [=#conserve_volume_flow_mode '''conserve_volume_flow_mode'''] }}} {{{#!td style="vertical-align:top" C*16 }}} {{{#!td style="vertical-align:top" 'default' }}} {{{#!td Modus of volume flow conservation.\\\\ The following values are allowed:\\\\ '' 'default' ''\\\\ Per default, PALM uses '' 'initial_profiles' '' for cyclic lateral boundary conditions ([#bc_lr bc_lr] = '' 'cyclic' '' and [#bc_ns bc_ns] = '' 'cyclic{{{'}}}'') and '' 'inflow_profile' '' for non-cyclic lateral boundary conditions (bc_lr /= '' 'cyclic' '' or bc_ns /= '' 'cyclic{{{'}}}'').\\\\ '' 'initial_profiles' ''\\\\ The target volume flow is calculated at t=0 from the initial profiles of u and v. This setting is only allowed for cyclic lateral boundary conditions (bc_lr = '' 'cyclic' '' and bc_ns = '' 'cyclic{{{'}}}'').\\\\ '' 'inflow_profile' ''\\\\ The target volume flow is calculated at every timestep from the inflow profile of u or v, respectively. This setting is only allowed for non-cyclic lateral boundary conditions (bc_lr /= '' 'cyclic' '' or bc_ns /= '' 'cyclic{{{'}}}'').\\\\ '' 'bulk_velocity' ''\\\\ The target volume flow is calculated from a predefined bulk velocity (see [#u_bulk u_bulk] and [#v_bulk v_bulk]). This setting is only allowed for cyclic lateral boundary conditions (bc_lr = '' 'cyclic' '' and bc_ns = '' 'cyclic{{{'}}}'').\\\\ Note that '''conserve_volume_flow_mode''' only comes into effect if [#conserve_volume_flow conserve_volume_flow] = ''.T.''. }}} |---------------- {{{#!td style="vertical-align:top" [=#coupling_start_time '''coupling_start_time'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Simulation time of precursor run.\\\\ Sets the time period a precursor run shall run uncoupled. This parameter is used to set up the precursor run control for atmosphere-ocean-coupled runs. It has to be set individually to the atmospheric / oceanic precursor run. The time in the data output will show negative values during the precursor run. See documentation for further information. }}} |---------------- {{{#!td style="vertical-align:top" [=#dp_external '''dp_external'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .F. }}} {{{#!td External pressure gradient switch.\\\\ This parameter is used to switch on/off an external pressure gradient as driving force. The external pressure gradient is controlled by the parameters [#dp_smooth dp_smooth], [#dp_level_b dp_level_b] and [#dpdxy dpdxy].\\\\ Note that '''dp_external''' = ''.T.'' requires [#conserve_volume_flow conserve_volume_flow] = ''.F.''. It is normally recommended to disable the Coriolis force by setting [#omega omega] = ''0.0''. }}} |---------------- {{{#!td style="vertical-align:top" [=#dp_smooth '''dp_smooth'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .F. }}} {{{#!td Vertically smooth the external pressure gradient using a sinusoidal smoothing function.\\\\ This parameter only applies if [#dp_external dp_external] = ''.T.''. It is useful in combination with [#dp_level_b dp_level_b] >> 0 to generate a non-accelerated boundary layer well below dp_level_b. }}} |---------------- {{{#!td style="vertical-align:top" [=#dp_level_b '''dp_level_b'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Lower limit of the vertical range for which the external pressure gradient is applied (in m).\\\\ This parameter only applies if [#dp_external dp_external] = ''.T.''. It must hold the condition zu(0) <= '''dp_level_b''' <= zu([#nz nz]). It can be used in combination with [#dp_smooth dp_smooth] = ''.T.'' to generate a non-accelerated boundary layer well below '''dp_level_b''' if '''dp_level_b''' >> 0.\\\\ Note that there is no upper limit of the vertical range because the external pressure gradient is always applied up to the top of the model domain. }}} |---------------- {{{#!td style="vertical-align:top" [=#dpdxy '''dpdxy'''] }}} {{{#!td style="vertical-align:top" R(2) }}} {{{#!td style="vertical-align:top" 2 * 0.0 }}} {{{#!td Values of the external pressure gradient applied in x- and y-direction, respectively (in Pa/m).\\\\ This parameter only applies if [#dp_external dp_external] = ''.T.''. It sets the pressure gradient values. Negative values mean an acceleration, positive values mean deceleration. For example, '''dpdxy''' = ''-0.0002, 0.0,'' drives the flow in positive x-direction. }}} |---------------- {{{#!td style="vertical-align:top" [=#dt_coupling '''dt_coupling'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 9999999.9 }}} {{{#!td Temporal interval for the data exchange in case of runs with [[coupled models]] (e.g. atmosphere - ocean) (in s).\\\\ This parameter has an effect only in case of a run with coupled models. It is available starting from version 3.3a.\\\\ This parameter specifies the temporal interval at which data are exchanged at the interface between coupled models (currently: interface between atmosphere and ocean). If this parameter is not explicitly specified in the parameter files for both coupled models, or if there is an inconsistency between its values for both coupled models, the execution will terminate and an informative [[error message]] will be given. In order to ensure synchronous coupling throughout the simulation, '''dt_coupling''' should be chosen larger than [../d3par#dt_max dt_max]. }}} |---------------- {{{#!td style="vertical-align:top" [=#e_init '''e_init'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Initial subgrid-scale TKE in m^2^s^-2^.\\\\ This option prescribes an initial subgrid-scale TKE from which the initial diffusion coefficients K,,m,, and K,,h,, will be calculated if '''e_init''' is positive. This option only has an effect if [#km_constant km_constant] is not set. }}} |---------------- {{{#!td style="vertical-align:top" [=#e_min '''e_min'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Minimum subgrid-scale TKE in m^2^s^-2^.\\\\ This option adds artificial viscosity to the flow by ensuring that the subgrid-scale TKE does not fall below the minimum threshold '''e_min'''. }}} |---------------- {{{#!td style="vertical-align:top" [=#galilei_transformation '''galilei_transformation'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .F. }}} {{{#!td Application of a Galilei-transformation to the coordinate system of the model.\\\\ With '''galilei_transformation''' = ''.T.,'' a so-called Galilei-transformation is switched on which ensures that the coordinate system of the model is moved along with the geostrophical wind. Alternatively, the model domain can be moved along with the averaged horizontal wind (see [#use_ug_for_galilei_tr use_ug_for_galilei_tr], this can and will naturally change in time). With this method, numerical inaccuracies of the Piascek - Williams - scheme (concerns in particular the momentum advection) are minimized. Beyond that, in the majority of cases the lower relative velocities in the moved system permit a larger time step ([#dt dt]). Switching the transformation on is only worthwhile if the geostrophical wind (ug, vg) and the averaged horizontal wind clearly deviate from the value 0. In each case, the distance the coordinate system has been moved is written to the file [../iofiles#RUN_CONTROL RUN_CONTROL].\\\\ Non-cyclic lateral boundary conditions (see [#bc_lr bc_lr] and [#bc_ns bc_ns]), the specification of a gestrophic wind that is not constant with height as well as e.g. stationary inhomogeneities at the bottom boundary do not allow the use of this transformation. }}} |---------------- {{{#!td style="vertical-align:top" [=#humidity '''humidity'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .F. }}} {{{#!td Parameter to switch on the prognostic equation for specific humidity q.\\\\ The initial vertical profile of q can be set via parameters [#q_surface q_surface], [#q_vertical_gradient q_vertical_gradient] and [#q_vertical_gradient_level q_vertical_gradient_level]. Boundary conditions can be set via [#q_surface_initial_change q_surface_initial_change] and [#surface_waterflux surface_waterflux].\\\\ If the condensation scheme is switched on ([#cloud_physics cloud_physics] = ''.T.''), q becomes the total liquid water content (sum of specific humidity and liquid water content). }}} |---------------- {{{#!td style="vertical-align:top" [=#km_constant '''km_constant'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" variable (computed from TKE) }}} {{{#!td Constant eddy diffusivities are used (laminar simulations).\\\\ If this parameter is specified, both in the 1d and in the 3d-model constant values for the eddy diffusivities are used in space and time with K,,m,, = '''km_constant''' and K,,h,, = K,,m,, / [[prandtl_number]]. The prognostic equation for the subgrid-scale TKE is switched off. Constant eddy diffusivities are only allowed with the Prandtl layer ([#prandtl_layer prandtl_layer]) switched off. }}} |---------------- {{{#!td style="vertical-align:top" [=#km_damp_max '''km_damp_max'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.5*([#dx dx] or [#dy dy]) }}} {{{#!td Maximum diffusivity used for filtering the velocity field in the vicinity of the outflow (in m^2^/s).\\\\ When using non-cyclic lateral boundaries (see [#bc_lr bc_lr] or [#bc_ns bc_ns]), a smoothing has to be applied to the velocity field in the vicinity of the outflow in order to suppress any reflections of outgoing disturbances. Smoothing is done by increasing the eddy diffusivity along the horizontal direction which is perpendicular to the outflow boundary. Only velocity components parallel to the outflow boundary are filtered (e.g. v and w, if the outflow is along x). Damping is applied from the bottom to the top of the domain.\\\\ The horizontal range of the smoothing is controlled by [#outflow_damping_width outflow_damping_width] which defines the number of gridpoints (counted from the outflow boundary) from where on the smoothing is applied. Starting from that point, the eddy diffusivity is linearly increased (from zero to its maximum value given by '''km_damp_max''') until half of the damping range width, from where it remains constant up to the outflow boundary. If at a certain grid point the eddy diffusivity calculated from the flow field is larger than as described above, it is used instead.\\\\ The default value of '''km_damp_max''' has been empirically proven to be sufficient. }}} |---------------- {{{#!td style="vertical-align:top" [=#ocean '''ocean'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .F. }}} {{{#!td Parameter to switch on ocean runs.\\\\ By default PALM is configured to simulate atmospheric flows. However, starting from version 3.3, '''ocean''' = ''.T.'' allows simulation of ocean turbulent flows. Setting this switch has several effects:\\\\ * An additional prognostic equation for salinity is solved.[[BR]] * Potential temperature in buoyancy and stability-related terms is replaced by potential density.[[BR]] * Potential density is calculated from the equation of state for seawater after each timestep, using the algorithm proposed by Jackett et al. (2006, J. Atmos. Oceanic Technol., '''23''', 1709-1728).[[BR]] So far, only the initial hydrostatic pressure is entered into this equation.[[BR]] * z=0 (sea surface) is assumed at the model top (vertical grid index k=[#nzt nzt] on the w-grid), with negative values of z indicating the depth.[[BR]] * Initial profiles are constructed (e.g. from [#pt_vertical_gradient pt_vertical_gradient] / [#pt_vertical_gradient_level pt_vertical_gradient_level]) starting from the sea surface, using surface values given by [#pt_surface pt_surface], [#sa_surface sa_surface], [#ug_surface ug_surface], and [#vg_surface vg_surface].[[BR]] * Zero salinity flux is used as default boundary condition at the bottom of the sea.[[BR]] * If switched on, random perturbations are by default imposed to the upper model domain from zu(nzt*2/3) to zu(nzt-3).\\\\ Relevant parameters to be exclusively used for steering ocean runs are [#bc_sa_t bc_sa_t], [#bottom_salinityflux bottom_salinityflux], [#sa_surface sa_surface], [#sa_vertical_gradient sa_vertical_gradient], [#sa_vertical_gradient_level sa_vertical_gradient_level], and [#top_salinityflux top_salinityflux].\\\\ Section [[4.4.2]] gives an example for appropriate settings of these and other parameters neccessary for ocean runs.\\\\ '''ocean''' = ''.T.'' does not allow settings of [#timestep_scheme timestep_scheme] = '' 'leapfrog' '' or '' 'leapfrog+euler' '' as well as [#scalar_advec scalar_advec] = '' 'ups-scheme'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#passive_scalar '''passive_scalar'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .F. }}} {{{#!td Parameter to switch on the prognostic equation for a passive scalar.\\\\ The initial vertical profile of ''s'' can be set via parameters [#s_surface s_surface], [#s_vertical_gradient s_vertical_gradient] and [#s_vertical_gradient_level s_vertical_gradient_level]. Boundary conditions can be set via [#s_surface_initial_change s_surface_initial_change] and [#surface_scalarflux surface_scalarflux].\\\\ '''Note:'''\\ With '''passive_scalar''' switched on, the simultaneous use of humidity (see [#humidity humidity]) is impossible. }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} [[BR]] [=#grid '''Grid:]\\ ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN]\\[../fortrantypes Type]''' =||='''Default\\Value''' =||='''Explanation''' =|| |---------------- {{{#!td style="vertical-align:top;width: 150px" [=#dx '''dx'''] }}} {{{#!td style="vertical-align:top;width: 50px" R }}} {{{#!td style="vertical-align:top;width: 75px" 1.0 }}} {{{#!td Horizontal grid spacing along the x-direction (in m).\\\\ Along x-direction only a constant grid spacing is allowed.\\\\ For coupled runs (see [[3.8]]) this parameter must be equal in both parameter files [../iofiles#PARIN PARIN] and [../iofiles#PARIN_O PARIN_O]. }}} |---------------- {{{#!td style="vertical-align:top" [=#dy '''dy'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 1.0 }}} {{{#!td Horizontal grid spacing along the y-direction (in m).\\\\ Along y-direction only a constant grid spacing is allowed.\\\\ For coupled runs (see [[3.8]]) this parameter must be equal in both parameter files [../iofiles#PARIN PARIN] and [../iofiles#PARIN_O PARIN_O]. }}} |---------------- {{{#!td style="vertical-align:top" [=#dz '''dz'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" }}} {{{#!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 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\\\\ The w-levels lie half between them:\\\\ zw(k) = ( zu(k) + zu(k+1) ) * 0.5 }}} |---------------- {{{#!td style="vertical-align:top" [=#dz_max '''dz_max'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 9999999.9 }}} {{{#!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. }}} |---------------- {{{#!td style="vertical-align:top" [=#dz_stretch_factor '''dz_stretch_factor'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 1.08 }}} {{{#!td Stretch factor for a vertically stretched grid (see [#dz_stretch_level dz_stretch_level]).\\\\ The stretch factor should not exceed a value of approx. ''1.10 - 1.12'', otherwise the discretization errors due to the stretched grid are not negligible any more. (refer Kalnay de Rivas) }}} |---------------- {{{#!td style="vertical-align:top" [=#dz_stretch_level '''dz_stretch_level'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 100000.0 }}} {{{#!td Height level above/below which the grid is to be stretched vertically (in m).\\\\ For [#ocean ocean] = ''.F.,'' '''dz_stretch_level''' is the height level (in m) above which the grid is to be stretched vertically. The vertical grid spacings [#dz dz] above this level are calculated as\\\\ dz(k+1) = dz(k) * [#dz_stretch_factor dz_stretch_factor]\\\\ 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''' 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]. }}} |---------------- {{{#!td style="vertical-align:top" [=#grid_matching '''grid_matching'''] }}} {{{#!td style="vertical-align:top" C*6 }}} {{{#!td style="vertical-align:top" 'strict' }}} {{{#!td Variable to adjust the subdomain sizes in parallel runs.\\\\ For '''grid_matching''' = '' 'strict' '', the subdomains are forced to have an identical size on all processors. In this case the processor numbers in the respective directions of the virtual processor net must fulfill certain divisor conditions concerning the grid point numbers in the three directions (see [#nx nx], [#ny ny] and [#nz nz]). Advantage of this method is that all PEs bear the same computational load.\\\\ There is no such restriction by default, because then smaller subdomains are allowed on those processors which form the right and/or north boundary of the virtual processor grid. On all other processors the subdomains are of same size. Whether smaller subdomains are actually used, depends on the number of processors and the grid point numbers used. Information about the respective settings are given in file [../iofiles#RUN_CONTROL RUN_CONTROL].\\\\ When using a multi-grid method for solving the Poisson equation (see [#psolver psolver]) only '''grid_matching''' = '' 'strict' '' is allowed.\\\\ '''Note:'''\\ In some cases for small processor numbers there may be a very bad load balancing among the processors which may reduce the performance of the code. }}} |---------------- {{{#!td style="vertical-align:top" [=#npex '''npex'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" }}} {{{#!td Number of processors along x-direction of the virtual processor net.\\\\ For parallel runs, the total number of processors to be used is given by the '''mrun'''-option -X. By default, depending on the type of the parallel computer, PALM generates a 1d processor net (domain decomposition along x, [#npey npey] = ''1'') or a 2d-net (this is favored on machines with fast communication network and/or large number of processors (>256)). In case of a 2d-net, it is tried to make it more or less square-shaped. If, for example, 16 processors are assigned (-X 16), a 4 * 4 processor net is generated ('''npex''' = 4, '''npey''' = 4). This choice is optimal for square total domains ([#nx nx] = [#ny ny]), since then the number of ghost points at the lateral boundarys of the subdomains reaches a minimum. If nx and ny differ extremely, the processor net should be manually adjusted using adequate values for '''npex''' and '''npey'''.\\\\ '''Important:'''\\\\ The value of '''npex''' * '''npey''' must exactly match the value assigned by the '''mrun'''-option -X. Otherwise the model run will abort with a corresponding error message.\\\\ Additionally, the specification of '''npex''' and '''npey''' may of course override the default setting for the domain decomposition (1d or 2d) which may have a significant (negative) effect on the code performance. }}} |---------------- {{{#!td style="vertical-align:top" [=#npey '''npey'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" }}} {{{#!td Number of processors along y-direction of the virtual processor net.\\\\ For further information see [#npex npex]. }}} |---------------- {{{#!td style="vertical-align:top" [=#nx '''nx'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" }}} {{{#!td Number of grid points in x-direction.\\\\ A value for this parameter must be assigned. Since the lower array bound in PALM starts with i = 0, the actual number of grid points is equal to '''nx+1'''. In case of cyclic boundary conditions along x, the domain size is ('''nx+1''')* [#dx dx].\\\\ For parallel runs, in case of [#grid_matching grid_matching] = '' 'strict','' '''nx+1''' must be an integral multiple of the processor numbers (see [#npex npex] and [#npey npey]) along x- as well as along y-direction (due to data transposition restrictions).\\\\ For coupled runs this parameter must be equal in both parameter files [#../iofiles#PARIN PARIN] and [../iofiles#PARIN_O PARIN_O]. }}} |---------------- {{{#!td style="vertical-align:top" [=#ny '''ny'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" }}} {{{#!td Number of grid points in y-direction.\\\\ A value for this parameter must be assigned. Since the lower array bound in PALM starts with j = 0, the actual number of grid points is equal to '''ny+1'''. In case of cyclic boundary conditions along y, the domain size is ('''ny+1''') * [#dy dy].\\\\ For parallel runs, in case of [#grid_matching grid_matching] = '' 'strict','' '''ny+1''' must be an integral multiple of the processor numbers (see [#npex npex] and [#npey npey]) along y- as well as along x-direction (due to data transposition restrictions).\\\\ For coupled runs this parameter must be equal in both parameter files [#../iofiles#PARIN PARIN] and [../iofiles#PARIN_O PARIN_O]. }}} |---------------- {{{#!td style="vertical-align:top" [=#nz '''nz'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" }}} {{{#!td Number of grid points in z-direction.\\\\ A value for this parameter must be assigned. Since the lower array bound in PALM starts with k = 0 and since one additional grid point is added at the top boundary (k = '''nz+1'''), the actual number of grid points is '''nz+2'''. However, the prognostic equations are only solved up to '''nz''' (u, v) or up to '''nz-1''' (w, scalar quantities). The top boundary for u and v is at k = '''nz+1''' (u, v) while at k = '''nz''' for all other quantities.\\\\ For parallel runs, in case of [#grid_matching grid_matching] = '' 'strict','' '''nz''' must be an integral multiple of the number of processors in x-direction (due to data transposition restrictions). }}} |---------------- {{{#!td style="vertical-align:top" [=#nz_do3d '''nz_do3d'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" nz+1 }}} {{{#!td Limits the output of 3d volume data along the vertical direction (grid point index k).\\\\ By default, data for all grid points along z are output. The parameter '''nz_do3d''' can be used to limit the output up to a certain vertical grid point (e.g. in order to reduce the amount of output data). It affects all output of volume data ("normal" output to file, see [../d3par#data_output data_output], as well as output for '''dvrp'''-software, see [../d3par#mode_dvrp mode_dvrp]). }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} [[BR]] [=#num '''Numerics:]\\ ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN]\\[../fortrantypes Type]''' =||='''Default\\Value''' =||='''Explanation''' =|| |---------------- {{{#!td style="vertical-align:top;width: 150px" [=#call_psolver_at_all_substeps '''call_psolver_at_all_substeps'''] }}} {{{#!td style="vertical-align:top;width: 50px" L }}} {{{#!td style="vertical-align:top;width: 75px" ''.T.'' }}} {{{#!td Switch to steer the call of the pressure solver.\\\\ In order to speed-up performance, the Poisson equation for perturbation pressure (see [#psolver psolver]) can be called only at the last substep of multistep Runge-Kutta timestep schemes (see [#timestep_scheme timestep_scheme]) by setting '''call_psolver_at_all_substeps''' = ''.F.''. In many cases, this sufficiently reduces the divergence of the velocity field. Nevertheless, small-scale ripples (2-delta-x) may occur. In this case and in case of non-cyclic lateral boundary conditions, '''call_psolver_at_all_timesteps''' = ''.T.'' should be used. }}} |---------------- {{{#!td style="vertical-align:top" [=#cfl_factor '''cfl_factor'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.1, 0.8 or 0.9 (see right) }}} {{{#!td Time step limiting factor.\\\\ In the model, the maximum allowed time step according to CFL and diffusion-criterion dt_max is reduced by [#dt dt] = dt_max * '''cfl_factor''' in order to avoid stability problems which may arise in the vicinity of the maximum allowed timestep. The condition ''0.0'' < '''cfl_factor''' < ''1.0'' applies.\\\\ The default value of '''cfl_factor''' depends on the [#timestep_scheme timestep_scheme] used:\\\\ For the third order Runge-Kutta scheme it is '''cfl_factor''' = ''0.9.''\\\\ In case of the leapfrog scheme a quite restrictive value of '''cfl_factor''' = ''0.1'' is used because for larger values the velocity divergence significantly effects the accuracy of the model results. Possibly larger values may be used with the leapfrog scheme but these are to be determined by appropriate test runs.\\\\ The default value for the Euler scheme is '''cfl_factor''' = ''0.8.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#cycle_mg '''cycle_mg'''] }}} {{{#!td style="vertical-align:top" C*1 }}} {{{#!td style="vertical-align:top" 'w' }}} {{{#!td Type of cycle to be used with the multi-grid method.\\\\ This parameter determines which type of cycle is applied in the multi-grid method used for solving the Poisson equation for perturbation pressure (see [#psolver psolver]). It defines in which way it is switched between the fine and coarse grids. So-called v- and w-cycles are realized (i.e. '''cycle_m'''g may be assigned the values '' 'v' '' or '' 'w' ''). The computational cost of w-cycles is much higher than that of v-cycles, however, w-cycles give a much better convergence. }}} |---------------- {{{#!td style="vertical-align:top" [=#fft_method '''fft_method'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'system-specific' }}} {{{#!td FFT-method to be used.\\\\ The fast fourier transformation (FFT) is used for solving the perturbation pressure equation with a direct method (see [#psolver psolver]) and for calculating power spectra (see optional software packages, [[4.2]]).\\\\ By default, system-specific, optimized routines from external vendor libraries are used. However, these are available only on certain computers and there are more or less severe restrictions concerning the number of gridpoints to be used with them.\\\\ There are two other PALM internal methods available on every machine (their respective source code is part of the PALM source code):\\\\ 1.: The '''Temperton'''-method from Clive Temperton (ECWMF) which is computationally very fast and switched on with '''fft_method''' = '' 'temperton-algorithm'.'' The number of horizontal gridpoints ([#nx nx]+1, [#ny ny]+1) to be used with this method must be composed of prime factors 2, 3 and 5.\\\\ 2.: The '''Singleton'''-method which is very slow but has no restrictions concerning the number of gridpoints to be used with, switched on with '''fft_method''' = '' 'singleton-algorithm'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#long_filter_factor '''long_filter_factor'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Filter factor for the so-called Long-filter.\\\\ This filter very efficiently eliminates 2-delta-waves sometimes caused by the upstream-spline scheme (see Mahrer and Pielke, 1978: Mon. Wea. Rev., '''106''', 818-830). It works in all three directions in space. A value of '''long_filter_factor''' = ''0.01'' sufficiently removes the small-scale waves without affecting the longer waves.\\\\ By default, the filter is switched off (= ''0.0''). It is exclusively applied to the tendencies calculated by the upstream-spline scheme (see [#momentum_advec momentum_advec] and [#scalar_advec scalar_advec]), not to the prognostic variables themselves. At the bottom and top boundary of the model domain the filter effect for vertical 2-delta-waves is reduced. There, the amplitude of these waves is only reduced by approx. 50%, otherwise by nearly 100%.\\\\ Filter factors with values > ''0.01'' also reduce the amplitudes of waves with wavelengths longer than 2-delta (see the paper by Mahrer and Pielke, quoted above). }}} |---------------- {{{#!td style="vertical-align:top" [=#loop_optimization '''loop_optimization'''] }}} {{{#!td style="vertical-align:top" C*16 }}} {{{#!td style="vertical-align:top" see right }}} {{{#!td Method used to optimize loops for solving the prognostic equations.\\\\ By default, the optimization method depends on the host on which PALM is running. On machines with vector-type CPUs, single 3d-loops are used to calculate each tendency term of each prognostic equation, while on all other machines, all prognostic equations are solved within one big loop over the two horizontal indices {{{i}}} and {{{j}}} (giving a good cache uitilization).\\\\ The default behaviour can be changed by setting either '''loop_optimization''' = '' 'vector' '' or '''loop_optimization''' = '' 'cache'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#mg_cycles '''mg_cycles'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" -1 }}} {{{#!td Number of cycles to be used with the multi-grid scheme.\\\\ This parameter determines the number of cycles to be carried out in the multi-grid method used for solving the Poisson equation for perturbation pressure (see [#psolver psolver]). The type of the cycles can be set with [#cycle_mg cycle_mg].\\\\ By default ('''mg_cyles''' = -''1''), the number of cycles depends on the requested accuracy of the scheme (see [#residual_limit residual_limit]) and may vary from time step to time step. In this case, the CPU time for a run will be difficult to estimate, since it heavily depends on the total number of the cycles to be carried out.\\\\ By assigning '''mg_cycles''' a value (>= ''1''), the number of cycles can be fixed so that the CPU time can be clearly estimated.\\\\ '''Note:'''\\ When using a fixed number of cycles, the user must examine the local file [../iofiles#RUN_CONTROL RUN_CONTROL] regularly to check whether the divergence of the velocity field is sufficiently reduced by the pressure solver. It should be reduced at least by two orders of magnitude. For cyclic boundary conditions along both horizontal directions (see [#bc_lr bc_lr] and [#bc_ns bc_ns]) '''mg_cycles''' = ''2'' is typically a good choice, for non-cyclic lateral boundary conditions '''mg_cycles''' = ''4'' may be sufficient. }}} |---------------- {{{#!td style="vertical-align:top" [=#mg_switch_to_pe0_level '''mg_switch_to_pe0_level'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" }}} {{{#!td Grid level at which data shall be gathered on PE0.\\\\ In case of a run using several PEs and the multigrid method for solving the Poisson equation for perturbation pressure (see [#psolver psolver]), the value of this parameter defines on which grid level the data are gathered on PE0 in order to allow for a further coarsening of the grid. The finest grid defines the largest grid level. By default, the gathering level is determined automatically and displayed in file [../iofiles#RUN_CONTROL RUN_CONTROL]. It is only possible to gather data from a level larger than the one determined automatically. A test run may be neccessary to determine this level.\\\\ Setting of '''mg_switch_to_pe0_level''' = -''1'' prevents that data are collected on PE0 at all, i.e. coarsening of grids is limited by the subdomains. }}} |---------------- {{{#!td style="vertical-align:top" [=#momentum_advec '''momentum_advec'''] }}} {{{#!td style="vertical-align:top" C*10 }}} {{{#!td style="vertical-align:top" 'pw-scheme' }}} {{{#!td Advection scheme to be used for the momentum equations.\\\\ The user can choose between the following schemes:\\\\ '' 'pw-scheme' ''\\\\ The scheme of Piascek and Williams (1970, J. Comp. Phys., 6, 392-405) with central differences in the form C3 is used. If intermediate Euler-timesteps are carried out in case of [#timestep_scheme timestep_scheme] = '' 'leapfrog+euler' '' the advection scheme is - for the Euler-timestep - automatically switched to an upstream-scheme.\\\\ '' 'ups-scheme' ''\\\\ The upstream-spline scheme is used (see Mahrer and Pielke, 1978: Mon. Wea. Rev., 106, 818-830). In opposite to the Piascek-Williams scheme, this is characterized by much better numerical features (less numerical diffusion, better preservation of flow structures, e.g. vortices), but computationally it is much more expensive. In addition, the use of the Euler-timestep scheme is mandatory ([#timestep_scheme timestep_scheme] = '' 'euler' ''), i.e. the timestep accuracy is only of first order. For this reason the advection of scalar variables (see [#scalar_advec scalar_advec]) should then also be carried out with the upstream-spline scheme, because otherwise the scalar variables would be subject to large numerical diffusion due to the upstream scheme.\\\\ Since the cubic splines used tend to overshoot under certain circumstances, this effect must be adjusted by suitable filtering and smoothing (see [#cut_spline_overshoot cut_spline_overshoot], [#long_filter_factor long_filter_factor], [#ups_limit_pt ups_limit_pt], [#ups_limit_u ups_limit_u], [#ups_limit_v ups_limit_v], [#ups_limit_w ups_limit_w]). This is always neccessary for runs with stable stratification, even if this stratification appears only in parts of the model domain.\\\\ With stable stratification the upstream-spline scheme also produces gravity waves with large amplitude, which must be suitably damped (see [../d3par#rayleigh_damping_factor rayleigh_damping_factor]).\\\\ '''Important:'''\\ The upstream-spline scheme is not implemented for humidity and passive scalars (see [#humidity humidity] and [#passive_scalar passive_scalar]) and requires the use of a 2d-domain-decomposition. The last conditions severely restricts code optimization on several machines leading to very long execution times! The scheme is also not allowed for non-cyclic lateral boundary conditions (see [#bc_lr bc_lr] and [#bc_ns bc_ns]). }}} |---------------- {{{#!td style="vertical-align:top" [=#ngsrb '''ngsrb'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" 2 }}} {{{#!td Grid level at which data shall be gathered on PE0.\\\\ In case of a run using several PEs and the multigrid method for solving the Poisson equation for perturbation pressure (see [#psolver psolver]), the value of this parameter defines on which grid level the data are gathered on PE0 in order to allow for a further coarsening of the grid. The finest grid defines the largest grid level. By default, the gathering level is determined automatically and displayed in file [../iofiles#RUN_CONTROL RUN_CONTROL]. It is only possible to gather data from a level larger than the one determined automatically. A test run may be neccessary to determine this level. }}} |---------------- {{{#!td style="vertical-align:top" [=#nsor '''nsor'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" 20 }}} {{{#!td Number of iterations to be used with the SOR-scheme.\\\\ This parameter is only effective if the SOR-scheme is selected as pressure solver ([#psolver psolver] = '' 'sor' ''). The number of iterations necessary for a sufficient convergence of the scheme depends on the grid point numbers and is to be determined by appropriate test runs (the default value will not at all be sufficient for larger grid point numbers). The number of iterations used for the first call of the SOR-scheme (t = 0) is determined via the parameter [#nsor_ini nsor_ini]. }}} |---------------- {{{#!td style="vertical-align:top" [=#nsor_ini '''nsor_ini'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" 100 }}} {{{#!td Initial number of iterations with the SOR algorithm.\\\\ This parameter is only effective if the SOR algorithm was selected as the pressure solver scheme ([#psolver psolver] = '' 'sor' '') and specifies the number of initial iterations of the SOR scheme (at t = 0). The number of subsequent iterations at the following timesteps is determined with the parameter [#nsor nsor]. Usually nsor < '''nsor_ini''', since in each case subsequent calls to [#psolver psolver] use the solution of the previous call as initial value. Suitable test runs should determine whether sufficient convergence of the solution is obtained with the default value and if necessary the value of '''nsor_ini''' should be changed. }}} |---------------- {{{#!td style="vertical-align:top" [=#omega_sor '''omega_sor'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 1.8 }}} {{{#!td Convergence factor to be used with the the SOR-scheme.\\\\ If the SOR-scheme is selected ([#psolver psolver] = '' 'sor' ''), this parameter determines the value of the convergence factor, where ''1.0'' <= '''omega_sor''' < ''2.0.'' The optimum value of '''omega_sor''' depends on the number of grid points along the different directions in space. For non-equidistant grids it can only be determined by appropriate test runs. }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} [[BR]] [=#phys '''Physics:]\\ ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN]\\[../fortrantypes Type]''' =||='''Default\\Value''' =||='''Explanation''' =|| |---------------- {{{#!td style="vertical-align:top;width: 150px" [=#omega '''omega'''] }}} {{{#!td style="vertical-align:top;width: 50px" R }}} {{{#!td style="vertical-align:top;width: 75px" 7.29212E-5 }}} {{{#!td Angular velocity of the rotating system (in rad s^-1^).\\\\ The angular velocity of the earth is set by default. The values of the Coriolis parameters are calculated as:\\\\ f = 2.0 * '''omega''' * sin([#phi phi])\\ f* = 2.0 * '''omega''' * cos(phi) }}} |---------------- {{{#!td style="vertical-align:top" [=#phi '''phi'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 55.0 }}} {{{#!td Geographical latitude (in degrees).\\\\ The value of this parameter determines the value of the Coriolis parameters f and f*, provided that the angular velocity (see [#omega omega]) is non-zero. }}} |---------------- {{{#!td style="vertical-align:top" [=#prandtl_number '''prandtl_number'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 1.0 }}} {{{#!td Ratio of the eddy diffusivities for momentum and heat (K,,m,,/K,,h,,).\\\\ For runs with constant eddy diffusivity (see [#km_constant km_constant]), this parameter can be used to assign the Prandtl number (ratio K,,m,, / K,,h,,). }}} [[BR]] [=#bc '''Boundary Conditions:]\\ ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN]\\[../fortrantypes Type]''' =||='''Default\\Value''' =||='''Explanation''' =|| |---------------- {{{#!td style="vertical-align:top;width: 150px" [=#adjust_mixing_length '''adjust_mixing_length'''] }}} {{{#!td style="vertical-align:top;width: 50px" L }}} {{{#!td style="vertical-align:top;width: 75px" .F. }}} {{{#!td Near-surface adjustment of the mixing length to the Prandtl-layer law.\\\\ Usually the mixing length in LES models l,,LES,, depends (as in PALM) on the grid size and is possibly restricted further in case of stable stratification and near the lower wall (see parameter [#wall_adjustment wall_adjustment]). With '''adjust_mixing_length''' = ''.T.'' the Prandtl' mixing length l,,PR,, = kappa * z/phi is calculated and the mixing length actually used in the model is set l = MIN (l,,LES,, , l,,PR,,). This usually gives a decrease of the mixing length at the bottom boundary and considers the fact that eddy sizes decrease in the vicinity of the wall.\\\\ '''Warning:''' So far, there is no good experience with '''adjust_mixing_length''' = ''.T.''{{{!}}}\\\\ With '''adjust_mixing_length''' = ''.T.'' and the Prandtl-layer being switched on (see [#prandtl_layer prandtl_layer]) '' '(u*){{{**}}} 2+neumann' '' should always be set as the lower boundary condition for the TKE (see [#bc_e_b bc_e_b]), otherwise the near-surface value of the TKE is not in agreement with the Prandtl-layer law (Prandtl-layer law and Prandtl-Kolmogorov-Ansatz should provide the same value for K,,m,,). A warning is given, if this is not the case. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_e_b '''bc_e_b'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'neumann' }}} {{{#!td Bottom boundary condition of the TKE.\\\\ '''bc_e_b''' may be set to '' 'neumann' '' or '' '(u*){{{**}}}2+neumann' ''. '''bc_e_b''' = '' 'neumann' '' yields to e(k=0)=e(k=1) (Neumann boundary condition), where e(k=1) is calculated via the prognostic TKE equation. Choice of '' '(u*){{{**}}}2+neumann' '' also yields to e(k=0)=e(k=1), but the TKE at the Prandtl-layer top (k=1) is calculated diagnostically by e(k=1)=(us/0.1){{{**}}}2. However, this is only allowed if a Prandtl-layer is used ([#prandtl_layer prandtl_layer]). If this is not the case, a warning is given and '''bc_e_b''' is reset to '' 'neumann' ''.\\\\ At the top boundary a Neumann boundary condition is generally used: (e(nz+1) = e(nz)). }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_lr '''bc_lr'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'cyclic' }}} {{{#!td Boundary condition along x (for all quantities).\\\\ By default, a cyclic boundary condition is used along x.\\\\ '''bc_lr''' may also be assigned the values '' 'dirichlet/radiation' '' (inflow from left, outflow to the right) or '' 'radiation/dirichlet' '' (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 [../d3par#psolver psolver]) and it also requires cyclic boundary conditions along y (see [#bc_ns bc_ns]).\\\\ In case of these non-cyclic lateral boundaries, a Dirichlet condition is used at the inflow for all quantities (initial vertical profiles - see [#initializing_actions initializing_actions] - 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.\\\\ When using non-cyclic lateral boundaries, a filter is applied to the velocity field in the vicinity of the outflow in order to suppress any reflections of outgoing disturbances (see [#km_damp_max km_damp_max] and [#outflow_damping_width outflow_damping_width]).\\\\ In order to maintain a turbulent state of the flow, it may be neccessary to continuously impose perturbations on the horizontal velocity field in the vicinity of the inflow throughout the whole run. This can be switched on using [../d3par#create_disturbances create_disturbances]. The horizontal range to which these perturbations are applied is controlled by the parameters [#inflow_disturbance_begin inflow_disturbance_begin] and [#inflow_disturbance_end inflow_disturbance_end]. The vertical range and the perturbation amplitude are given by [../d3par#disturbance_level_b disturbance_level_b], [../d3par#disturbance_level_t disturbance_level_t], and [../d3par#disturbance_amplitude disturbance_amplitude]. The time interval at which perturbations are to be imposed is set by [../d3par#dt_disturb dt_disturb].\\\\ In case of non-cyclic horizontal boundaries [../d3par#call_psolver_at_all_substeps call_psolver_at_all_substeps] = ''.T.'' should be used. '''Note:''' Using non-cyclic lateral boundaries requires very sensitive adjustments of the inflow (vertical profiles) and the bottom boundary conditions, e.g. a surface heating should not be applied near the inflow boundary because this may significantly disturb the inflow. Please check the model results very carefully. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_ns '''bc_ns'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'cyclic' }}} {{{#!td Boundary condition along y (for all quantities).\\\\ By default, a cyclic boundary condition is used along y.\\\\ '''bc_ns''' may also be assigned the values '' 'dirichlet/radiation' '' (inflow from rear ("north"), outflow to the front ("south")) or '' 'radiation/dirichlet' '' (inflow from front ("south"), outflow to the rear ("north")). This requires the multi-grid method to be used for solving the Poisson equation for perturbation pressure (see [../d3par#psolver psolver]) and it also requires cyclic boundary conditions along x (see [#bc_lr bc_lr]).\\\\ In case of these non-cyclic lateral boundaries, a Dirichlet condition is used at the inflow for all quantities (initial vertical profiles - see [#initializing_actions initializing_actions] - 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.\\\\ For further details regarding non-cyclic lateral boundary conditions see [#bc_lr bc_lr]. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_p_b '''bc_p_b'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'neumann' }}} {{{#!td Bottom boundary condition of the perturbation pressure.\\\\ Allowed values are '' 'dirichlet' '', '' 'neumann' '' and '' 'neumann+inhomo' ''. '' 'dirichlet' '' sets p(k=0)=0.0, '' 'neumann' '' sets p(k=0)=p(k=1). '' 'neumann+inhomo' '' corresponds to an extended Neumann boundary condition where heat flux or temperature inhomogeneities near the surface (pt(k=1)) are additionally regarded (see Shen and LeClerc (1995, Q.J.R. Meteorol. Soc., 1209)). This condition is only permitted with the Prandtl-layer switched on ([#prandtl_layer prandtl_layer]), otherwise the run is terminated.\\\\ Since at the bottom boundary of the model the vertical velocity disappears (w(k=0) = 0.0), the consistent Neumann condition ('' 'neumann' '' or '' 'neumann+inhomo' '') dp/dz = 0 should be used, which leaves the vertical component w unchanged when the pressure solver is applied. Simultaneous use of the Neumann boundary conditions both at the bottom and at the top boundary ([#bc_p_t bc_p_t]) usually yields no consistent solution for the perturbation pressure and should be avoided. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_p_t '''bc_p_t'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'dirichlet' }}} {{{#!td Top boundary condition of the perturbation pressure.\\\\ Allowed values are '' 'dirichlet' '' (p(k=nz+1)= 0.0) or '' 'neumann' '' (p(k=nz+1)=p(k=nz)).\\\\ Simultaneous use of Neumann boundary conditions both at the top and bottom boundary ([#bc_p_b bc_p_b]) usually yields no consistent solution for the perturbation pressure and should be avoided. Since at the bottom boundary the Neumann condition is a good choice (see [#bc_p_b bc_p_b]), a Dirichlet condition should be set at the top boundary. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_pt_b '''bc_pt_b'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'dirichlet' }}} {{{#!td Bottom boundary condition of the potential temperature.\\\\ Allowed values are '' 'dirichlet' '' (pt(k=0) = const. = [#pt_surface pt_surface] + [#pt_surface_initial_change pt_surface_initial_change]; the user may change this value during the run using [#user-defined] code) and '' 'neumann' '' (pt(k=0)=pt(k=1)). When a constant surface sensible heat flux is used ([#surface_heatflux surface_heatflux]), '''bc_pt_b''' = '' 'neumann' '' must be used, because otherwise the resolved scale may contribute to the surface flux so that a constant value cannot be guaranteed.\\\\ In the coupled atmosphere executable, '''bc_pt_b''' is internally set and does not need to be prescribed. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_pt_t '''bc_pt_t'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'initial_gradient' }}} {{{#!td Top boundary condition of the potential temperature.\\\\ Allowed are the values '' 'dirichlet' '' (pt(k=nz+1) does not change during the run), '' 'neumann' '' (pt(k=nz+1)=pt(k=nz)), and '' 'initial_gradient' ''. With the '' 'initial_gradient' ''-condition the value of the temperature gradient at the top is calculated from the initial temperature profile (see [#pt_surface pt_surface], [#pt_vertical_gradient pt_vertical_gradient]) by bc_pt_t_val = (pt_init(k=nz+1) - pt_init(k=nz)) / dzu(nz+1). Using this value (assumed constant during the run) the temperature boundary values are calculated as pt(k=nz+1) = pt(k=nz) + bc_pt_t_val * dzu(nz+1) (up to k=nz the prognostic equation for the temperature is solved). When a constant sensible heat flux is used at the top boundary ([#top_heatflux top_heatflux]), '''bc_pt_t''' = '' 'neumann' '' must be used, because otherwise the resolved scale may contribute to the top flux so that a constant value cannot be guaranteed. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_q_b '''bc_q_b'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'dirichlet' }}} {{{#!td Bottom boundary condition of the specific humidity / total water content.\\\\ Allowed values are '' 'dirichlet' '' (q(k=0) = const. = [#q_surface q_surface] + [#q_surface_initial_change q_surface_initial_change]; the user may change this value during the run using [#user-defined] code) and '' 'neumann' '' (q(k=0)=q(k=1)). When a constant surface latent heat flux is used ([#surface_waterflux surface_waterflux]), '''bc_q_b''' = '' 'neumann' '' must be used, because otherwise the resolved scale may contribute to the surface flux so that a constant value cannot be guaranteed. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_q_t '''bc_q_t'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'neumann' }}} {{{#!td Top boundary condition of the specific humidity / total water content.\\\\ Allowed are the values '' 'dirichlet' '' (q(k=nz) and q(k=nz+1) do not change during the run) and '' 'neumann' ''. With the Neumann boundary condition the value of the humidity gradient at the top is calculated from the initial humidity profile (see [#q_surface q_surface], [#q_vertical_gradient q_vertical_gradient]) by: bc_q_t_val = ( q_init(k=nz) - q_init(k=nz-1)) / dzu(nz). Using this value (assumed constant during the run) the humidity boundary values are calculated as q(k=nz+1) =q(k=nz) + bc_q_t_val * dzu(nz+1) (up tp k=nz the prognostic equation for q is solved). }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_s_b '''bc_s_b'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'dirichlet' }}} {{{#!td Bottom boundary condition of the scalar concentration.\\\\ Allowed values are '' 'dirichlet' '' (s(k=0) = const. = [#s_surface s_surface] + [#s_surface_initial_change s_surface_initial_change]; the user may change this value during the run using [#user-defined] code) and '' 'neumann' '' (s(k=0) = s(k=1)). When a constant surface concentration flux is used ([#surface_scalarflux surface_scalarflux]), '''bc_s_b''' = '' 'neumann' '' must be used, because otherwise the resolved scale may contribute to the surface flux so that a constant value cannot be guaranteed. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_s_t '''bc_s_t'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'neumann' }}} {{{#!td Top boundary condition of the scalar concentration.\\\\ Allowed are the values '' 'dirichlet' '' (s(k=nz) and s(k=nz+1) do not change during the run) and '' 'neumann' ''. With the Neumann boundary condition the value of the scalar concentration gradient at the top is calculated from the initial scalar concentration profile (see [#s_surface s_surface], [#s_vertical_gradient s_vertical_gradient]) by: bc_s_t_val = (s_init(k=nz) - s_init(k=nz-1)) / dzu(nz). Using this value (assumed constant during the run) the concentration boundary values are calculated as s(k=nz+1) = s(k=nz) + bc_s_t_val * dzu(nz+1) (up to k=nz the prognostic equation for the scalar concentration is solved). }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_sa_t '''bc_sa_t'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'neumann' }}} {{{#!td Top boundary condition of the salinity.\\\\ This parameter only comes into effect for ocean runs (see parameter [#ocean ocean]).\\\\ Allowed are the values '' 'dirichlet' '' (sa(k=nz+1) does not change during the run) and '' 'neumann' '' (sa(k=nz+1)=sa(k=nz)).\\\\ When a constant salinity flux is used at the top boundary ([#top_salinityflux top_salinityflux]), '''bc_sa_t''' = '' 'neumann' '' must be used, because otherwise the resolved scale may contribute to the top flux so that a constant value cannot be guaranteed. }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_uv_b '''bc_uv_b'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'dirichlet' }}} {{{#!td Bottom boundary condition of the horizontal velocity components u and v.\\\\ Allowed values are '' 'dirichlet' '' and '' 'neumann' ''. '''bc_uv_b''' = '' 'dirichlet' '' yields the no-slip condition with u=v=0 at the bottom. Due to the staggered grid u(k=0) and v(k=0) are located at z = - 0,5 * [#dz dz] (below the bottom), while u(k=1) and v(k=1) are located at z = +0,5 * dz. u=v=0 at the bottom is guaranteed using mirror boundary condition: u(k=0) = - u(k=1) and v(k=0) = - v(k=1) The Neumann boundary condition yields the free-slip condition with u(k=0) = u(k=1) and v(k=0) = v(k=1). With Prandtl - layer switched on (see [#prandtl_layer prandtl_layer]), the free-slip condition is not allowed (otherwise the run will be terminated). }}} |---------------- {{{#!td style="vertical-align:top" [=#bc_uv_t '''bc_uv_t'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'dirichlet' }}} {{{#!td Top boundary condition of the horizontal velocity components u and v.\\\\ Allowed values are '' 'dirichlet' '', '' 'dirichlet_0' '' and '' 'neumann' ''. The Dirichlet condition yields u(k=nz+1) = ug(nz+1) and v(k=nz+1) = vg(nz+1), Neumann condition yields the free-slip condition with u(k=nz+1) = u(k=nz) and v(k=nz+1) = v(k=nz) (up to k=nz the prognostic equations for the velocities are solved). The special condition '' 'dirichlet_0' '' can be used for channel flow, it yields the no-slip condition u(k=nz+1) = ug(nz+1) = 0 and v(k=nz+1) = vg(nz+1) = 0. In the coupled ocean executable, bc_uv_t is internally set ('neumann') and does not need to be prescribed. }}} |---------------- {{{#!td style="vertical-align:top" [=#bottom_salinityflux '''bottom_salinityflux'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Kinematic salinity flux near the surface (in psu m/s).\\\\ This parameter only comes into effect for ocean runs (see parameter [#ocean ocean]).\\\\ The respective salinity flux value is used as bottom (horizontally homogeneous) boundary condition for the salinity equation. This additionally requires that a Neumann condition must be used for the salinity, which is currently the only available condition. }}} |---------------- {{{#!td style="vertical-align:top" [=#inflow_damping_height '''inflow_damping_height'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" from precursor run }}} {{{#!td Height below which the turbulence signal is used for turbulence recycling (in m).\\\\ In case of a turbulent inflow (see [#turbulent_inflow turbulent_inflow]), this parameter defines the vertical thickness of the turbulent layer up to which the turbulence extracted at the recycling plane (see [#recycling_width recycling_width]) shall be imposed to the inflow. Above this level the turbulence signal is linearly damped to zero. The transition range within which the signal falls to zero is given by the parameter [#inflow_damping_width inflow_damping_width].\\\\ By default, this height is set as the height of the convective boundary layer as calculated from a precursor run. See chapter [[3.9]] about proper settings for getting this CBL height from a precursor run. }}} |---------------- {{{#!td style="vertical-align:top" [=#inflow_damping_width '''inflow_damping_width'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.1 * [#inflow_damping_height inflow_damping]\\ [#inflow_damping_height _height] }}} {{{#!td Transition range within which the turbulance signal is damped to zero (in m).\\\\ See [#inflow_damping_height inflow_damping_height] for explanation. }}} |---------------- {{{#!td style="vertical-align:top" [=#inflow_disturbance_begin '''inflow_disturbance_begin'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" MIN(10, nx/2 or ny/2) }}} {{{#!td Lower limit of the horizontal range for which random perturbations are to be imposed on the horizontal velocity field (gridpoints).\\\\ If non-cyclic lateral boundary conditions are used (see [#bc_lr bc_lr] or [#bc_ns bc_ns]), this parameter gives the gridpoint number (counted horizontally from the inflow) from which on perturbations are imposed on the horizontal velocity field. Perturbations must be switched on with parameter [../d3par#create_disturbances create_disturbances]. }}} |---------------- {{{#!td style="vertical-align:top" [=#inflow_disturbance_end '''inflow_disturbance_end'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" MIN(100, 3/4*nx or 3/4*ny) }}} {{{#!td Upper limit of the horizontal range for which random perturbations are to be imposed on the horizontal velocity field (gridpoints).\\\\ If non-cyclic lateral boundary conditions are used (see [#bc_lr bc_lr] or [#bc_ns bc_ns]), this parameter gives the gridpoint number (counted horizontally from the inflow) up to which perturbations are imposed on the horizontal velocity field. Perturbations must be switched on with parameter [../d3par#create_disturbances create_disturbances]. }}} |---------------- {{{#!td style="vertical-align:top" [=#outflow_damping_width '''outflow_damping_width'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" MIN(20, [#nx nx]/2 or [#ny ny]/2) }}} {{{#!td Width of the damping range in the vicinity of the outflow (gridpoints).\\\\ When using non-cyclic lateral boundaries (see [#bc_lr bc_lr] or [#bc_ns bc_ns]), a smoothing has to be applied to the velocity field in the vicinity of the outflow in order to suppress any reflections of outgoing disturbances. This parameter controlls the horizontal range to which the smoothing is applied. The range is given in gridpoints counted from the respective outflow boundary. For further details about the smoothing see parameter [#km_damp_max km_damp_max], which defines the magnitude of the damping. }}} |---------------- {{{#!td style="vertical-align:top" [=#prandtl_layer '''prandtl_layer'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .T. }}} {{{#!td Parameter to switch on a Prandtl layer.\\\\ By default, a Prandtl layer is switched on at the bottom boundary between z = 0 and z = 0.5 * [#dz dz] (the first computational grid point above ground for u, v and the scalar quantities). In this case, at the bottom boundary, free-slip conditions for u and v (see [#bc_uv_b bc_uv_b]) are not allowed. Likewise, laminar simulations with constant eddy diffusivities ([#km_constant km_constant]) are forbidden.\\\\ With Prandtl-layer switched off, the TKE boundary condition [#bc_e_b bc_e_b] = '' '(u*){{{**}}}2+neumann' '' must not be used and is automatically changed to '' 'neumann' '' if necessary. Also, the pressure boundary condition [#bc_p_b bc_p_b] = '' 'neumann+inhomo' '' is not allowed.\\\\ If the Prandtl-layer is switched off and fluxes shall be prescribed at the surface (by setting [#surface_heatflux surface_heatflux]), it is required to set the parameter [#use_surface_fluxes use_surface_fluxes] = ''.T.''.\\\\ The roughness length is declared via the parameter [#roughness_length roughness_length]. }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} [[BR]] [=#ini '''Initialization:]\\ ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN]\\[../fortrantypes Type]''' =||='''Default\\Value''' =||='''Explanation''' =|| |---------------- {{{#!td style="vertical-align:top;width: 150px" [=#damp_level_1d '''damp_level_1d'''] }}} {{{#!td style="vertical-align:top;width: 50px" R }}} {{{#!td style="vertical-align:top;width: 75px" zu(nz+1) }}} {{{#!td Height where the damping layer begins in the [[1d-model]] (in m).\\\\ This parameter is used to switch on a damping layer for the 1d-model, which is generally needed for the damping of inertia oscillations. Damping is done by gradually increasing the value of the eddy diffusivities about 10% per vertical grid level (starting with the value at the height given by '''damp_level_1d''', or possibly from the next grid pint above), i.e. K,,m,,(k+1) = 1.1 * K,,m,,(k). The values of K,,m,, are limited to 10 m^2^/s at maximum.\\\\ This parameter only comes into effect if the 1d-model is switched on for the initialization of the 3d-model using [#initializing_actions initializing_actions] = '' 'set_1d-model_profiles'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#dissipation_1d '''dissipation_1d'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'as_in_3d_model' }}} {{{#!td Calculation method for the energy dissipation term in the TKE equation of the [[1d-model]].\\\\ By default the dissipation is calculated as in the 3d-model using diss = (0.19 + 0.74 * l / l_grid) * e{{{**}}}1.5 / l.\\\\ Setting '''dissipation_1d''' = '' 'detering' '' forces the dissipation to be calculated as diss = 0.064 * e{{{**}}}1.5 / l. }}} |---------------- {{{#!td style="vertical-align:top" [=#dt_pr_1d '''dt_pr_1d'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 9999999.9 }}} {{{#!td Temporal interval of vertical profile output of the [[1d-model]] (in s).\\\\ Data are written in ASCII format to file [../iofiles#LIST_PROFIL_1D LIST_PROFIL_1D]. This parameter is only in effect if the 1d-model has been switched on for the initialization of the 3d-model with [#initializing_actions initializing_actions] = '' 'set_1d-model_profiles'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#dt_run_control_1d '''dt_run_control_1d'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 60.0 }}} {{{#!td Temporal interval of runtime control output of the [[1d-model]] (in s).\\\\ Data are written in ASCII format to file [../iofiles#RUN_CONTROL RUN_CONTROL]. This parameter is only in effect if the 1d-model is switched on for the initialization of the 3d-model with [#initializing_actions initializing_actions] = '' 'set_1d-model_profiles'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#end_time_1d '''end_time_1d'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 864000.0 }}} {{{#!td Time to be simulated for the [[1d-model]] (in s).\\\\ The default value corresponds to a simulated time of 10 days. Usually, after such a period the inertia oscillations have completely decayed and the solution of the 1d-model can be regarded as stationary (see [#damp_level_1d damp_level_1d]). This parameter is only in effect if the 1d-model is switched on for the initialization of the 3d-model with [#initializing_actions initializing_actions] = '' 'set_1d-model_profiles'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#initializing_actions '''initializing_actions'''] }}} {{{#!td style="vertical-align:top" C*100 }}} {{{#!td style="vertical-align:top" }}} {{{#!td Initialization actions to be carried out.\\\\ This parameter does not have a default value and therefore must be assigned with each model run. For restart runs '''initializing_actions''' = '' 'read_restart_data' '' must be set. For the initial run of a job chain the following values are allowed:\\\\ '' 'set_constant_profiles' ''\\\\ A horizontal wind profile consisting of linear sections (see [#ug_surface ug_surface], [#ug_vertical_gradient ug_vertical_gradient], [#ug_vertical_gradient_level ug_vertical_gradient_level] and [#vg_surface vg_surface], [#vg_vertical_gradient vg_vertical_gradient], [#vg_vertical_gradient_level vg_vertical_gradient_level], respectively) as well as a vertical temperature (humidity) profile consisting of linear sections (see [#pt_surface pt_surface], [#pt_vertical_gradient pt_vertical_gradient], [#q_surface q_surface] and [#q_vertical_gradient q_vertical_gradient]) are assumed as initial profiles. The subgrid-scale TKE is set to 0 but K,,m,, and K,,h,, are set to very small values because otherwise no TKE would be generated.\\\\ '' 'set_1d-model_profiles' ''\\\\ The arrays of the 3d-model are initialized with the (stationary) solution of the [[1d-model]]. These are the variables e, K,,h,,{{{,}}} K,,m,,{{{,}}} u, v and with Prandtl layer switched on rif, us, usws, vsws. The temperature (humidity) profile consisting of linear sections is set as for '' 'set_constant_profiles' '' and assumed as constant in time within the 1d-model. For steering of the 1d-model a set of parameters with suffix "_1d" (e.g. [#end_time_1d end_time_1d], [#damp_level_1d damp_level_1d]) is available.\\\\ '' 'by_user' ''\\\\ The initialization of the arrays of the 3d-model is under complete control of the user and has to be done in routine {{{user_init_3d_model}}} of the user-interface (see [[3.5.1]]).\\\\ '' 'initialize_vortex' ''\\\\ The initial velocity field of the 3d-model corresponds to a Rankine-vortex with vertical axis. This setting may be used to test advection schemes. Free-slip boundary conditions for u and v (see [#bc_uv_b bc_uv_b], [#bc_uv_t bc_uv_t]) are necessary. In order not to distort the vortex, an initial horizontal wind profile constant with height is necessary (to be set by [#initializing_actions initializing_actions] = '' 'set_constant_profiles{{{'}}}'') and some other conditions have to be met (neutral stratification, diffusion must be switched off, see [#km_constant km_constant]). The center of the vortex is located at jc = ([#nx nx]+1)/2. It extends from k = 0 to k = [#nz nz]+1. Its radius is 8 * [#dx dx] and the exponentially decaying part ranges to 32 * dx (see {{{init_rankine.f90}}}).\\\\ '' 'initialize_ptanom' ''\\\\ A 2d-Gauss-like shape disturbance (x,y) is added to the initial temperature field with radius 10.0 * [#dx dx] and center at jc = ([#nx nx]+1)/2. This may be used for tests of scalar advection schemes (see [#scalar_advec scalar_advec]). Such tests require a horizontal wind profile constant with hight and diffusion switched off (see '' 'initialize_vortex' ''). Additionally, the buoyancy term must be switched of in the equation of motion for w (this requires the user to comment out the call of buoyancy in the source code of {{{prognostic_equations.f90}}}).\\\\ '' 'cyclic_fill' ''\\\\ Here, 3d-data from a precursor run are read by the initial (main) run. The precursor run is allowed to have a smaller domain along x and y compared with the main run. Also, different numbers of processors can be used for these two runs. Limitations are that the precursor run must use cyclic horizontal boundary conditions and that the number of vertical grid points, [#nz nz], must be same for the precursor run and the main run. If the total domain of the main run is larger than that of the precursor run, the domain is filled by cyclic repetition of the (cyclic) precursor data. This initialization method is recommended if a turbulent inflow is used (see [#turbulent_inflow turbulent_inflow]). 3d-data must be made available to the run by activating an appropriate file connection statement for local file [../iofiles#BININ BININ]. See chapter [[3.9]] for more details, where usage of a turbulent inflow is explained.\\\\ Values may be combined, e.g. '''initializing_actions''' = '' 'set_constant_profiles initialize_vortex' '', but the values of '' 'set_constant_profiles' '', '' 'set_1d-model_profiles' '' , and '' 'by_user' '' must not be given at the same time. }}} |---------------- {{{#!td style="vertical-align:top" [=#mixing_length_1d '''mixing_length_1d'''] }}} {{{#!td style="vertical-align:top" C*20 }}} {{{#!td style="vertical-align:top" 'as_in_3d_model' }}} {{{#!td Mixing length used in the [[1d-model]].\\\\ By default the mixing length is calculated as in the 3d-model (i.e. it depends on the grid spacing).\\\\ By setting '''mixing_length_1d''' = '' 'blackadar','' 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). }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} [[BR]] [=#topo '''Topography:]\\ ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN]\\[../fortrantypes Type]''' =||='''Default\\Value''' =||='''Explanation''' =|| |---------------- {{{#!td style="vertical-align:top; text-align:left;width: 150px" [=#building_height '''building_height'''] }}} {{{#!td style="vertical-align:top; text-align:left;style="width: 50px" R }}} {{{#!td style="vertical-align:top; text-align:left;style="width: 75px" 50.0 }}} {{{#!td Height of a single building in m.\\\\ '''building_height''' must be less than the height of the model domain. This parameter requires the use of [#topography topography] = '' 'single_building'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#building_length_x '''building_length_x'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 50.0 }}} {{{#!td Width of a single building in m.\\\\ Currently, '''building_length_x''' must be at least 3 * [#dx dx] and no more than ( [#nx nx] - 1 ) * dx - [#building_wall_left building_wall_left]. This parameter requires the use of [#topography topography] = '' 'single_building'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#building_length_y '''building_length_y'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 50.0 }}} {{{#!td Depth of a single building in m.\\\\ Currently, '''building_length_y''' must be at least 3 * [#dy dy] and no more than ( [#ny ny] - 1 ) * dy - [#building_wall_south building_wall_south]. This parameter requires the use of [#topography topography] = '' 'single_building'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#building_wall_left '''building_wall_left'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" building centered in x-direction }}} {{{#!td x-coordinate of the left building wall (distance between the left building wall and the left border of the model domain) in m.\\\\ Currently, '''building_wall_left''' must be at least 1 * [#dx dx] and less than ( [#nx nx] - 1 ) * dx - [#building_length_x building_length_x]. This parameter requires the use of [#topography topography] = '' 'single_building'.''\\\\ The default value '''building_wall_left''' = ( ( nx + 1 ) * dx - building_length_x ) / 2 centers the building in x-direction. Due to the staggered grid the building will be displaced by -0.5 dx in x-direction and -0.5 [#dy dy] in y-direction. }}} |---------------- {{{#!td style="vertical-align:top" [=#building_wall_south '''building_wall_south'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" building centered in y-direction }}} {{{#!td y-coordinate of the South building wall (distance between the South building wall and the South border of the model domain) in m.\\\\ Currently, '''building_wall_south''' must be at least 1 * [#dy dy] and less than ( [#ny ny] - 1 ) * dy - [#building_length_y building_length_y]. This parameter requires the use of [#topography topography] = '' 'single_building'.''\\\\ The default value '''building_wall_south''' = ( ( ny + 1 ) * dy - building_length_y ) / 2 centers the building in y-direction. Due to the staggered grid the building will be displaced by -0.5 [#dx dx] in x-direction and -0.5 dy in y-direction. }}} |---------------- {{{#!td style="vertical-align:top" [=#canyon_height '''canyon_height'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 50.0 }}} {{{#!td Street canyon height in m.\\\\ '''canyon_height''' must be less than the height of the model domain. This parameter requires [#topography topography] = '' 'single_street_canyon'.'' }}} |---------------- {{{#!td style="vertical-align:top" [=#canyon_width_x '''canyon_width_x'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 9999999.9 }}} {{{#!td Street canyon width in x-direction in m.\\\\ Currently, '''canyon_width_x''' must be at least 3 * [#dx dx] and no more than ( [#nx nx] - 1 ) * dx - [#canyon_wall_left canyon_wall_left]. This parameter requires [#topography topography] = '' 'single_street_canyon'.'' A non-default value implies a canyon orientation in y-direction. }}} |---------------- {{{#!td style="vertical-align:top" [=#canyon_width_y '''canyon_width_y'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 9999999.9 }}} {{{#!td Street canyon width in y-direction in m.\\\\ Currently, '''canyon_width_y''' must be at least 3 * [#dy dy] and no more than ( [#ny ny] - 1 ) * dy - [#canyon_wall_south canyon_wall_south]. This parameter requires [#topography topography] = '' 'single_street_canyon'.'' A non-default value implies a canyon orientation in x-direction. }}} |---------------- {{{#!td style="vertical-align:top" [=#canyon_wall_left '''canyon_wall_left'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" canyon centered in x-direction }}} {{{#!td x-coordinate of the left canyon wall (distance between the left canyon wall and the left border of the model domain) in m.\\\\ Currently, '''canyon_wall_left''' must be at least 1 * [#dx dx] and less than ( [#nx nx] - 1 ) * dx - [#canyon_width_x canyon_width_x]. This parameter requires [#topography topography] = '' 'single_street_canyon'.''\\\\ The default value '''canyon_wall_left''' = ( ( nx + 1 ) * dx - canyon_width_x ) / 2 centers the canyon in x-direction. }}} |---------------- {{{#!td style="vertical-align:top" [=#canyon_wall_south '''canyon_wall_south'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" canyon centered in y-direction }}} {{{#!td y-coordinate of the South canyon wall (distance between the South canyon wall and the South border of the model domain) in m.\\\\ Currently, '''canyon_wall_south''' must be at least 1 * [#dy dy] and less than ( [#ny ny] - 1 ) * dy - [#canyon_width_y canyon_width_y]. This parameter requires [#topography topography] = '' 'single_street_canyon'.''\\\\ The default value '''canyon_wall_south''' = ( ( ny + 1 ) * dy - canyon_width_y ) / 2 centers the canyon in y-direction. }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} |---------------- {{{#!td style="vertical-align:top" [=#topography '''topography'''] }}} {{{#!td style="vertical-align:top" C*40 }}} {{{#!td style="vertical-align:top" '' 'flat' '' }}} {{{#!td Topography mode.\\\\ The user can choose between the following modes:\\\\ '' 'flat' '' Flat surface. '' 'single_building' '' Flow around a single rectangular building mounted on a flat surface. The building size and location can be specified by the parameters [#building_height building_height], [#building_length_x building_length_x], [#building_length_y building_length_y], [#building_wall_left building_wall_left] and [#building_wall_south building_wall_south]. '' 'single_street_canyon' '' Flow over a single, quasi-2D street canyon of infinite length oriented either in x- or in y-direction. The canyon size, orientation and location can be specified by the parameters [#canyon_height canyon_height] plus '''either''' [#canyon_width_x canyon_width_x] and [#canyon_wall_left canyon_wall_left] '''or''' [#canyon_width_y canyon_width_y] and [#canyon_wall_south canyon_wall_south]. '' 'read_from_file' '' Flow around arbitrary topography. This mode requires the input file [#TOPOGRAPHY_DATA TOPOGRAPHY_DATA]. This file contains the arbitrary topography height information in m. These data must exactly match the horizontal grid.\\\\ Alternatively, the user may add code to the user interface subroutine [#user_init_grid user_init_grid] to allow further topography modes. These require to explicitly set the [#topography_grid_convention topography_grid_convention] to either '' 'cell_edge' '' or '' 'cell_center' ''.\\\\ Non-flat '''topography''' modes may assign a kinematic sensible [#wall_heatflux wall_heatflux] and a kinematic [#wall_humidityflux wall_humidityflux] (requires [#humidity humidity] = .T.) or a [#wall_scalarflux wall_scalarflux] (requires [#passive_scalar passive_scalar] = .T.) at the five topography faces.\\\\ All non-flat '''topography''' modes require the use of [#momentum_advec momentum_advec] = [#scalar_advec scalar_advec] = '' 'pw-scheme' '', [#psolver psolver] /= '' 'sor' '', [#alpha_surface alpha_surface] = 0.0, [#galilei_transformation galilei_transformation] = ''.F.'', [#cloud_physics cloud_physics] = ''.F.'', [#cloud_droplets cloud_droplets] = ''.F.'', and [#prandtl_layer prandtl_layer] = ''.T.''.\\\\ Note that an inclined model domain requires the use of [#topography topography] = '' 'flat' '' and a nonzero alpha_surface. }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} [[BR]] [=#canopy '''Canopy:]\\ ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN]\\[../fortrantypes Type]''' =||='''Default\\Value''' =||='''Explanation''' =|| |---------------- {{{#!td style="vertical-align:top;width: 150px" [=#canopy_mode '''canopy_mode'''] }}} {{{#!td style="vertical-align:top;width: 50px" C*20 }}} {{{#!td style="vertical-align:top;width: 75px" 'block' }}} {{{#!td Canopy mode.\\\\ Besides using the default value, that will create a horizontally homogeneous plant canopy that extends over the total horizontal extension of the model domain, the user may add code to the user interface (see [#3.5.1 3.5.1]) subroutine {{{user_init_plant_canopy}}} to allow further canopy modes.\\\\ The setting of '''canopy_mode''' becomes only active, if [#plant_canopy plant_canopy] has been set ''.T.'' and a non-zero [#drag_coefficient drag_coefficient] has been defined. }}} |---------------- {{{#!td style="vertical-align:top" [=#cthf '''cthf'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Average heat flux that is prescribed at the top of the plant canopy.\\\\ If [#plant_canopy plant_canopy] is set ''.T.'', the user can prescribe a heat flux at the top of the plant canopy. It is assumed that solar radiation penetrates the canopy and warms the foliage which, in turn, warms the air in contact with it.\\\\ '''Note:''' Instead of using the value prescribed by [#surface_heatflux surface_heatflux], the near surface heat flux is determined from an exponential function that is dependent on the cumulative leaf_area_index (Shaw and Schumann (1992, Boundary Layer Meteorol., '''61''', 47-64)). }}} |---------------- {{{#!td style="vertical-align:top" [=#drag_coefficient '''drag_coefficient'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Drag coefficient used in the plant canopy model.\\\\ This parameter has to be non-zero, if the parameter [#plant_canopy plant_canopy] is set ''.T.''. }}} |---------------- {{{#!td style="vertical-align:top" [=#lad_surface '''lad_surface'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Surface value of the leaf area density (in m^2^/m^3^).\\\\ This parameter assigns the value of the leaf area density '''lad''' at the surface (k=0). Starting from this value, the leaf area density profile is constructed with [#lad_vertical_gradient lad_vertical_gradient] and [#lad_vertical_gradient_level lad_vertical_gradient_level]. }}} |---------------- {{{#!td style="vertical-align:top" [=#lad_vertical_gradient '''lad_vertical_gradient'''] }}} {{{#!td style="vertical-align:top" R(10) }}} {{{#!td style="vertical-align:top" 10 * 0.0 }}} {{{#!td Gradient(s) of the leaf area density (in m^2^/m^4^).\\\\ This leaf area density gradient holds starting from the height level defined by [#lad_vertical_gradient_level lad_vertical_gradient_level] (precisely: for all uv levels k where zu(k) > lad_vertical_gradient_level, lad(k) is set: lad(k) = lad(k-1) + dzu(k) * '''lad_vertical_gradient''') up to the level defined by [#pch_index pch_index]. Above that level lad(k) will automatically be set to 0.0. A total of 10 different gradients for 11 height intervals (10 intervals if lad_vertical_gradient_level(1) = 0.0) can be assigned. The leaf area density at the surface is assigned via [#lad_surface lad_surface]. }}} |---------------- {{{#!td style="vertical-align:top" [=#lad_vertical_gradient_level '''lad_vertical_gradient_level'''] }}} {{{#!td style="vertical-align:top" R(10) }}} {{{#!td style="vertical-align:top" 10 * 0.0 }}} {{{#!td Height level from which on the gradient of the leaf area density defined by [#lad_vertical_gradient lad_vertical_gradient] is effective (in m).\\\\ The height levels have to be assigned in ascending order. The default values result in a leaf area density that is constant with height up to the top of the plant canopy layer defined by [#pch_index pch_index]. For the piecewise construction of temperature profiles see [#lad_vertical_gradient lad_vertical_gradient]. }}} |---------------- {{{#!td style="vertical-align:top" [=#leaf_surface_concentration '''leaf_surface_concentration'''] }}} {{{#!td style="vertical-align:top" R }}} {{{#!td style="vertical-align:top" 0.0 }}} {{{#!td Concentration of a passive scalar at the surface of a leaf (in K m/s).\\\\ This parameter is only of importance in cases in that both, [#plant_canopy plant_canopy] and [#passive_scalar passive_scalar], are set ''.T..'' The value of the concentration of a passive scalar at the surface of a leaf is required for the parametrisation of the sources and sinks of scalar concentration due to the canopy. }}} |---------------- {{{#!td style="vertical-align:top" [=#pch_index '''pch_index'''] }}} {{{#!td style="vertical-align:top" I }}} {{{#!td style="vertical-align:top" 0 }}} {{{#!td Grid point index (scalar) of the upper boundary of the plant canopy layer.\\\\ Above '''pch_index''' the arrays of leaf area density and [#drag_coeffient drag_coeffient] are automatically set to zero in case of [#plant_canopy plant_canopy] = ''.T.''. Up to '''pch_index''' a leaf area density profile can be prescribed by using the parameters [#lad_surface lad_surface], [#lad_vertical_gradient lad_vertical_gradient] and [#lad_vertical_gradient_level lad_vertical_gradient_level]. }}} |---------------- {{{#!td style="vertical-align:top" [=#plant_canopy '''plant_canopy'''] }}} {{{#!td style="vertical-align:top" L }}} {{{#!td style="vertical-align:top" .F. }}} {{{#!td Switch for the plant_canopy_model.\\\\ If '''plant_canopy''' is set ''.T.'', the plant canopy model of Watanabe (2004, BLM 112, 307-341) is used.\\ The impact of a plant canopy on a turbulent flow is considered by an additional drag term in the momentum equations and an additional sink term in the prognostic equation for the subgrid-scale TKE. These additional terms are dependent on the leaf drag coefficient (see [#drag_coefficient drag_coefficient]) and the leaf area density (see [#lad_surface lad_surface], [#lad_vertical_gradient lad_vertical_gradient], [#lad_vertical_gradient_level lad_vertical_gradient_level). The top boundary of the plant canopy is determined by the parameter [#pch_index pch_index]. For all heights equal to or larger than zw(k=pch_index) the leaf area density is 0 (i.e. there is no canopy at these heights!).\\ By default, a horizontally homogeneous plant canopy is prescribed, if '''plant_canopy''' is set ''.T.''. However, the user can define other types of plant canopies (see [#canopy_mode canopy_mode]).\\\\ If '''plant_canopy''' and [#passive_scalar passive_scalar] are set ''.T.'', the canopy acts as an additional source or sink, respectively, of scalar concentration. The source/sink strength is dependent on the scalar concentration at the leaf surface, which is generally constant with time in PALM and which can be specified by specifying the parameter [#leaf_surface_concentration leaf_surface_concentration].\\\\ Additional heating of the air by the plant canopy is taken into account, when the default value of the parameter [#cthf cthf] is altered in the parameter file. In that case the value of [#surface_heatflux surface_heatflux] specified in the parameter file is not used in the model. Instead the near-surface heat flux is derived from an exponential function that is dependent on the cumulative leaf area index.\\\\ '''plant_canopy''' = ''.T.'' is only allowed together with a non-zero [#drag_coefficient drag_coefficient]. }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}} [[BR]] [=#others '''Others:]\\ ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN]\\[../fortrantypes Type]''' =||='''Default\\Value''' =||='''Explanation''' =|| |---------------- {{{#!td style="vertical-align:top; text-align:left;width: 150px" [=# ''''''] }}} {{{#!td style="vertical-align:top; text-align:left;style="width: 50px" }}} {{{#!td style="vertical-align:top; text-align:left;style="width: 75px" }}} {{{#!td }}} |---------------- {{{#!td style="vertical-align:top" [=# ''''''] }}} {{{#!td style="vertical-align:top" }}} {{{#!td style="vertical-align:top" }}} {{{#!td }}}