== Particle Parameters ==
[[TracNav(doc/app/partoc|nocollapse)]]
\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
'''NAMELIST group name: particles_par''' \\

||='''Parameter Name'''  =||='''[../fortrantypes FORTRAN Type]'''  =||='''Default Value'''  =||='''Explanation'''  =||
|----------------
{{{#!td style="vertical-align:top; text-align:left;width: 150px"
[=#alloc_factor '''alloc_factor''']
}}}
{{{#!td style="vertical-align:top; text-align:left;style="width: 50px"
R
}}}
{{{#!td style="vertical-align:top; text-align:left;style="width: 75px"
20.0
}}}
{{{#!td
Factor (in percent) describing the memory allocated additionally to the memory needed for initial particles at a given grid cell. For example, 100 particles are initialized at a grid cell. Hence, an array for 120 particles (default value of '''alloc_factor''') is allocated providing sufficient memory for particles being transported to this grid cell during the simulation. Note that this array is automatically reallocated if more than 120 particles (following the example above) are transported to this array in order to provide sufficient memory. Thus, it is not necessary to choose alloc_factor too high.
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#bc_par_b '''bc_par_b''']
}}}
{{{#!td style="vertical-align:top"
C*15
}}}
{{{#!td style="vertical-align:top"
'reflect'
}}}
{{{#!td
Bottom boundary condition for particle transport. \\
By default, particles are reflected at the bottom boundary. Alternatively, a particle absorption can set by '''bc_par_b''' = '' 'absorb'. ''
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#bc_par_lr '''bc_par_lr''']
}}}
{{{#!td style="vertical-align:top"
C*15
}}}
{{{#!td style="vertical-align:top"
'cyclic'
}}}
{{{#!td
Lateral boundary condition (x-direction) for particle transport. \\
By default, cyclic boundary conditions are used along x. Alternatively, reflection ('''bc_par_lr''' = '' 'reflect' '') or absorption ('''bc_par_lr''' = '' 'absorb' '') can be set. \\
This lateral boundary conditions should correspond to the lateral boundary condition used for the flow (see [../inipar#bc_lr bc_lr]).
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#bc_par_ns '''bc_par_ns''']
}}}
{{{#!td style="vertical-align:top"
C*15
}}}
{{{#!td style="vertical-align:top"
'cyclic'
}}}
{{{#!td
Lateral boundary condition (y-direction) for particle transport.\\
By default, cyclic boundary conditions are used along y. Alternatively, reflection ('''bc_par_ns''' = '' 'reflect' '') or absorption (bc_par_ns = '' 'absorb' '') can be set. \\
This lateral boundary conditions should correspond to the lateral boundary condition used for the flow (see [../inipar#bc_ns bc_ns]).
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#bc_par_t '''bc_par_t''']
}}}
{{{#!td style="vertical-align:top"
C*15
}}}
{{{#!td style="vertical-align:top"
'absorb'
}}}
{{{#!td
Top boundary condition for particle transport.\\
By default, particles are absorbed at the top boundary. Alternatively, a reflection condition can be set by '''bc_par_t''' = '' 'reflect'. ''
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#collision_kernel '''collision_kernel''']
}}}
{{{#!td style="vertical-align:top"
C*15
}}}
{{{#!td style="vertical-align:top"
'none'
}}}
{{{#!td
Parameter to steer cloud droplet growth by collision processes.\\\\

The growth of cloud droplets due to collision is parameterized using the so-called collision kernel. By default, collision is switched off. The user can choose between the following kernels:\\\\

'' 'hall' ''
      Collision kernel from Hall (1980, J. Atmos. Sci., 2486-2507), which considers collision due to pure gravitational effects. Larger droplets have a higher terminal fall velocity and are collecting smaller ones. Only terminal droplet velocities are considered in this kernel (not their effective velocities).

'' 'hall_fast' ''
      Same as '' 'hall' '', but a collision efficiency table is calculated only once (at the beginning of the simulation) for fixed radius classes in the range [1.0E-6,1000.0E-6] m. The number of classes to be used (i.e. the resolution of the kernel) can be set by parameter [#radius_classes radius_classes]. This method significantly reduces the total CPU time of a job.

'' 'none' ''
      Droplet collision is switched off.

'' 'wang' ''
      Beside gravitational effects (treated with the Hall-kernel) also the effects of turbulence on the collision are considered using parameterizations of Ayala et al. (2008, New J. Phys., 10, 075015) and Wang and Grabowski (2009, Atmos. Sci. Lett., 10, 1-8). This kernel includes three possible effects of turbulence: the modification of the relative velocity between the droplets, the effect of preferential concentration, and the enhancement of collision efficiencies.

'' 'wang_fast' ''
      Same as '' 'wang' '', but a collision efficiency table is calculated only once (at the beginning of the simulation) for fixed radius- and dissipation classes in the ranges [1.0E-6,1000.0E-6] m, and [0.0,600.0] cm**2/s**3 respectively. The number of classes to be used (i.e. the resolution of the kernel) can be set by parameters [#radius_classes radius_classes], and [#dissipation_classes dissipation_classes]. This method significantly reduces the total CPU time of a job.

'''Attention:''' Switching on the collision process drastically increases the CPU time of jobs.
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#collision_algorithm '''collision_algorithm''']
}}}
{{{#!td style="vertical-align:top"
C*15
}}}
{{{#!td style="vertical-align:top"
'all_or_nothing'
}}}
{{{#!td
Parameter to steer the algorithm for cloud droplet growth by collision.\\\\

By default, the collision algorithm is set to 'all_or_nothing'. The user can choose between the following algorithms:\\\\

'' 'all_or_nothing' ''
      Probabilistic collision algorithm based on the ideas of Shima et al. (2009) and Sölch and Kärcher (2010).  Each particles represented by one superdroplet grows by the collection of one particle of another superdroplet if the probability for this event if larger than a random number.

'' 'average_impact' ''
      Original PALM collision algorithm (Riechelmann et al, 2012), in which the average grow of every superdroplet is calculated. In contrast to the  'all_or_nothing' algorithm, the number of collected particles is equally distributed over the collecting particles, i.e., a particle might grow by collecting a certain fraction of particles.  

}}}
|----------------
{{{#!td style="vertical-align:top;width: 150px"
[=#curvature_solution_effects '''curvature_solution_effects''']
}}}
{{{#!td style="vertical-align:top;width: 50px"
L
}}}
{{{#!td style="vertical-align:top;width: 75px"
.F.
}}}
{{{#!td
Parameter to consider solution and curvature effects on the equilibrium vapor pressure of cloud droplets. For the initialization of the corresponding dry aerosol spectrum, see [#init_aerosol_probabilistic init_aerosol_probabilistic]. In this case, the radius growth equation is a stiff o.d.e, which is integrated in time using a Rosenbrock method (see Numerical Recipes in FORTRAN, 2nd Edition, p.731). '''curvature_solution_effects''' = ''.T.'' may significantly increase CPU time of jobs.
}}}
|----------------
{{{#!td style="vertical-align:top;width: 150px"
[=#deallocate_memory '''deallocate_memory''']
}}}
{{{#!td style="vertical-align:top;width: 50px"
L
}}}
{{{#!td style="vertical-align:top;width: 75px"
.T.
}}}
{{{#!td
Parameter to enable deallocation of unused memory. If the number of particles in a grid box exceeds the allocated memory, new memory is allocated. However, in case the number of particles per grid box in only temporarily high, most of the memory remains unused. If [#deallocate_memory deallocate_memory] = .True., 
the allocated memory used for particles will be dynamically adjusted with respect to the current number of particles every [#step_dealloc step_dealloc] 's timestep.  
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#density_ratio '''density_ratio''']
}}}
{{{#!td style="vertical-align:top"
R(10)
}}}
{{{#!td style="vertical-align:top"
0.0, 9 *
9999999.9
}}}
{{{#!td
Ratio of the density of the fluid and the density of the particles. \\
With the default value the particles are weightless and transported passively with the resolved scale flow. In case of '''density_ratio''' ''/= 0.0'' particles have a mass and hence inertia so that their velocity deviates more or less from the velocity of the surrounding flow. Particle velocity is calculated analytically and depends on (besides the density ratio and the current velocity difference between particles and surrounding fluid) the particle radius which is determined via [#radius radius] as well as on the molecular viscosity (assumed as 1.461E-5 m^2^/s).

If '''density_ratio''' ''= 1.0'', the particle density corresponds to the density of the surrounding fluid and the particles do not feel any buoyancy. Otherwise, particles will be accelerated upwards ('''density_ratio''' > ''1.0'') or downwards ('''density_ratio''' < ''1.0'').

With several groups of particles (see [#number_of_particle_groups number_of_particle_groups]), each group can be assigned a different value. If the number of values given for '''density_ratio''' is less than the number of groups defined by [#number_of_particle_groups number_of_particle_groups]), then the last assigned value is used for all remaining groups. This means that by default the particle density ratio for all groups will be ''0.0''.
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#dissipation_classes '''dissipation_classes''']
}}}
{{{#!td style="vertical-align:top"
I
}}}
{{{#!td style="vertical-align:top"
10
}}}
{{{#!td
Number of dissipation classes to be used in the collision efficiency table.\\
This parameter comes into effect if the parameter [#collision_kernel collision_kernel] is set to '' 'wang_fast' ''. It defines the number of dissipation classes which spawn the collision efficiency table. The interval [1.0,1000.0] cm**2/s**3 is divided into n (= '''dissipation_classes''') equidistant parts. 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#dt_dopts '''dt_dopts''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
value of [../d3par#dt_data_output dt_data] \\ [../d3par#dt_data_output _output] 
}}}
{{{#!td
Temporal interval at which time series data of particle quantities shall be output (in s). 
If particle advection is switched on (see [#dt_prel dt_prel]) this parameter can be used to assign the temporal interval at which time series of particle quantities shall be output. Output is written in netCDF format on local file [../io_file#DATA_1D_PTS_NETCDF DATA_1D_PTS_NETCDF].

The following list gives a short description of the quantities available. Most quantities are averaged over all particles. The quantity name given in the first column is identical to the respective name of the variable on the netCDF file.

In case of using more than one particle group (see [#number_of_particle_groups number_of_particle_groups]), separate time series are output for each of the groups. The long names of the variables in the netCDF file containing the respective time series all end with the string ''PG'' ##, where ## is the number of the particle group ''(01, 02, etc.)''. \\\\

||='''netCDF Variable Name''' =||='''Explanation''' =||
||tnpt   || total number of particles  ||
||x_    || particle x-coordinate with respect to the particle origin (in m)  ||
||y_    || particle y-coordinate with respect to the particle origin (in m)  ||
||z_    || particle z-coordinate with respect to the particle origin (in m)  ||
||z_abs    || absolute particle z-coordinate (in m)  ||
||u     || u particle velocity component (in m/s)  ||
||v     || v particle velocity component (in m/s)  ||
||w     || w particle velocity component (in m/s)  ||
||u''    || subgrid-scale u particle velocity component (in m/s)  ||
||v''    || subgrid-scale v particle velocity component (in m/s)  || 
||w''    || subgrid-scale w particle velocity component (in m/s)  ||
||npt_up    || total number of upward moving particles  ||
||w_up    || vertical velocity of the upward moving particles (in m/s)  ||
||w_down    || vertical velocity of the downward moving particles (in m/s)  ||
||npt_max    || maximum number of particles in a subdomain (=tnpt for non-parallel runs)  ||
||npt_min    || minimum number of particles in a subdomain (=tnpt for non-parallel runs)  ||
||x*2    || variance of the particle x-coordinate with respect to x_ (in m^2^)  ||
||y*2    || variance of the particle y-coordinate with respect to y_ (in m^2^)  ||
||z*2    || variance of the particle z-coordinate with respect to z_ (in m^2^)  ||
||u*2    || variance of the u particle velocity component with respect to u (in m^2^/s^2^)  ||
||v*2    || variance of the v particle velocity component with respect to v (in m^2^/s^2^)  ||
||w*2    || variance of the w particle velocity component with respect to w (in m^2^/s^2^)  ||
||u"2    || variance of the subgrid-scale u particle velocity component with respect to u" (in m^2^/s^2^)  ||
||v"2    || variance of the subgrid-scale v particle velocity component with respect to v" (in m^2^/s^2^)  ||
||w"2    || variance of the subgrid-scale w particle velocity component with respect to w" (in m^2^/s^2^)  ||
||npt*2    || variance of the number of particles with respect to the average number of particles per subdomain  ||                
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#dt_min_part '''dt_min_part''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
0.0002
}}}
{{{#!td
Minimum value for the particle time step when SGS velocities are used (in s).\\
For a further explanation see package parameter [#use_sgs_for_particles use_sgs_for_particles].
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#dt_prel '''dt_prel''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
9999999.9
}}}
{{{#!td
Temporal interval at which particles are to be released from a particle source (in s).\\  
By default, particles are released only at the beginning of a simulation (t_init=0). The time of the first release (t_init) can be changed with package parameter [#particle_advection_start particle_advection_start]. The time of the last release can be set with the package parameter [#end_time_prel end_time_prel]. If '''dt_prel''' has been set, additional releases will be at ''t = t_init+'''dt_prel''', t_init+2*'''dt_prel''', t_init+3*'''dt_prel''', etc..'' Actual release times may slightly deviate from these values (see e.g. [../d3par#dt_dopr dt_dopr]).\\\\
The domain of the particle source, as well as the distance of  released particles within this source, are determined via package parameters [#pst pst], [#psl psl], [#psr psr], [#pss pss], [#psn psn], [#psb psb], [#pdx pdx], [#pdy pdy] and [#pdz pdz]. By default, one particle is released at all points defined by these parameters. The package parameter [#particles_per_point particles_per_point] can be used to start more than one particle per point.\\\\
Up to 10 different groups of particles can be released at the same time (see [#number_of_particle_groups number_of_particle_groups]) where each group may have a different source. All particles belonging to one group have the same density ratio and the same radius. All other particle features (e.g. location of the source) are identical for all groups of particles.\\\\
Subgrid scale velocities can (optionally) be included for calculating the particle advection, using the method of Weil et al. (2004, JAS, 61, 2877-2887). This method is switched on by the package parameter [#use_sgs_for_particles use_sgs_for_particles]. This also forces the Euler/upstream method to be used for time advancement of the TKE (see initialization parameter [../inipar#use_upstream_for_tke use_upstream_for_tke]). The minimum time step during the sub-time steps is controlled by package parameter [#dt_min_part dt_min_part]. \\\\
By default, particles are weightless and transported passively with the resolved scale flow. Particles can be given a mass and thus an inertia by assigning the package parameter density_ratio a non-zero value (it defines the ratio of the density of the fluid and the density of the particles). In this case, their radius must also be defined, which affects their flow resistance. \\\\
Boundary conditions for the particle transport can be defined with package parameters [#bc_par_t bc_par_t], [#bc_par_lr bc_par_lr], [#bc_par_ns bc_par_ns] and [#bc_par_b bc_par_b].\\\\
Time series of particle quantities in netCDF format can be output to local file [../iofiles#DATA_1D_PTS_NETCDF DATA_1D_PTS_NETCDF] by using package parameter [#dt_dopts dt_dopts].\\\\
For analysis, additional output of particle information in equidistant temporal intervals can be carried out using [#dt_write_particle_data dt_write_particle_data] (file [../iofiles#PARTICLE_DATA PARTICLE_DATA]).\\\\
Statistical information (e.g. the total number of particles used, the number of particles exchanged between the PEs, etc.) are output to the local file [../iofiles#PARTICLE_INFOS PARTICLE_INFOS], if switched on by the parameter [#write_particle_statistics write_particle_statistics]. \\\\
If a job chain is to be carried out, particle information for the restart run (e.g. current location of all particles at the end of the run) is output to the local file [../iofiles#PARTICLE_RESTART_DATA_OUT PARTICLE_RESTART_DATA_OUT], which must be saved at the end of the run and given as an input file to the restart run under local file name [../iofiles#PARTICLE_RESTART_DATA_IN PARTICLE_RESTART_DATA_IN] using respective file connection statements in the '''mrun''' configuration file. \\\\
'''So far, the particle transport realized in PALM does only work duly in case of a constant vertical grid spacing! '''
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#dt_write_particle_data '''dt_write_particle_data''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
9999999.9
}}}
{{{#!td
Temporal interval for output of particle data (in s).\\
This parameter can be used to assign the temporal interval at which particle data shall be output. Data are output to the local file [../iofiles#PARTICLE_DATA PARTICLE_DATA]. See the file description for more details about its format.\\
By default, no particle data are output.
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#end_time_prel '''end_time_prel''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
9999999.9
}}}
{{{#!td
Time of the last release of particles (in s).\\
See also [#particle_advection_start particle_advection_start].
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#initial_weighting_factor '''initial_weighting\\_factor''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
1.0
}}}
{{{#!td
Factor to define the real number of initial droplets in a grid box.\\
In case of explicitly simulating cloud droplets (see [../inipar#cloud_droplets cloud_droplets]), the real number of initial droplets in a grid box is equal to the initial number of droplets in this box (defined by the particle source parameters [#pst pst], [#psl psl], [#psr psr], [#pss pss], [#psn psn], [#psb psb], [#pdx pdx], [#pdy pdy] and [#pdz pdz]) times the '''initial_weighting_factor'''.
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#init_aerosol_probabilistic '''init_aerosol_probabilistic''']
}}}
{{{#!td style="vertical-align:top"
L
}}}
{{{#!td style="vertical-align:top"
.FALSE.
}}}
{{{#!td
A logical which steers the initialization of the aerosol spectrum (only necessary if [#curvature_solution_effects curvature_solution_effects] are activated). Up to 3 log-normal distributions can be predefined to initialize the aerosol spectrum via [#n1 n1], [#n2 n2], [#n3 n3], [#rm1 rm1], [#rm2 rm2], [#rm3 rm3], [#s1 s1], [#s2 s2], [#s3 s3] in the dry aerosol radius range from 0.01 to 1.0 microns. In a subsaturated environment, the initial radius of the haze particle is computed by a parametrization, i.e., the parameter [#radius radius] does not affect the initial radius of the particles. If no aerosol spectrum is desired, see [#monodisperse_aerosols monodisperse_aerosols].

Options:
* .TRUE.: The aerosol dry radius is initialized by a random number generator. The weighting factor is not changed by this initialization.
* .FALSE.: The aerosol spectrum is divided in logarithmically-spaced bins (the number of bins equals the number of super-droplets per grid box). The dry aerosol radius of the super-droplet is set to the mean dry aerosol radius of the super-droplet's bin. The weighting factor is adjusted to be proportional to the number of aerosols in the bin, but the mean weighting factors still matches the [#initial_weighting_factor initial_weighting_factor].
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#max_number_particles_per_gridbox '''max_number_particles_per_gridbox''']
}}}
{{{#!td style="vertical-align:top"
I
}}}
{{{#!td style="vertical-align:top"
100
}}}
{{{#!td
Threshold for splitting of super droplets: If the number of particles in a grid box reached this value no more splitting operations are conducted. 

'''Remark:'''\\
It is important to set this threshold to an appropriate value (if [#splitting splitting]=.T.). Otherwise the unlimited creation of new particles in single grid boxes may significantly increase CPU time of jobs. 
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#min_nr_particle '''min_nr_particle''']
}}}
{{{#!td style="vertical-align:top"
I
}}}
{{{#!td style="vertical-align:top"
50
}}}
{{{#!td
Minimum number of particles for which memory is allocated at every grid cell.\\

'''Remark:'''\\
In case you are using larger number of gridpoints and if you release particles in a limited area of the total domain only, you may think about reducing the default value to a small number, e.g. 2. For example, if you release particles only near the surface, you may never need to allocate memory for particles well above the boundary layer. Furthermore, please be aware that you may run out of memory if you are using the default value of '''min_nr_particle''' for a larger number of grid points.
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#monodisperse_aerosols '''monodisperse_aerosols''']
}}}
{{{#!td style="vertical-align:top"
L
}}}
{{{#!td style="vertical-align:top"
.FALSE.
}}}
{{{#!td
* .TRUE.: Initializes a monodisperse aerosol spectrum, i.e., the dry aerosol radius is 0.1 micron for each super-droplet. 
* .FALSE.: A log-nomal distributed aerosol spectrum is initialized (see [#init_aerosol_probabilistic init_aerosol_probabilistic]).
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#number_concentration '''number_concentration'']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
-1.0
}}}
{{{#!td
Initial particle number concentration (in units of 1/cm³).\\
If this value is set, the initial weighting-factor is calculated that the prescribed initial number concentration is ensured in every grid box that contains particles.
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#number_of_particle_groups '''number_of_particle\\_groups''']
}}}
{{{#!td style="vertical-align:top"
I
}}}
{{{#!td style="vertical-align:top"
1
}}}
{{{#!td
Number of particle groups to be used.\\
Each particle group can be assigned its own source region (see [#pdx pdx], [#psl psl], [#psr psr], etc.), particle diameter ([#radius radius]) and particle density ratio ([#density_ratio density_ratio]).\\
If fewer values are given for [#pdx pdx], [#psl psl], etc. than the number of particle groups, then the last value is used for the remaining values (or the default value, if the user did not set the parameter).\\
The maximum allowed number of particle groups is limited to ''10.''
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#number_particles_per_gridbox '''number_particles_per\\_gridbox''']
}}}
{{{#!td style="vertical-align:top"
I
}}}
{{{#!td style="vertical-align:top"
-1
}}}
{{{#!td
Number of particles which are created in every grid box.\\ 
If a value >=1 is chosen, [#pdx pdx], [#pdy pdy], and [#pdz pdz] are calculated assuming an isotropic particle distribution. The default value -1 requires initialization by the user with parameters [#pdx pdx], [#pdy pdy], and [#pdz pdz].
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#n1 '''n1''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
100.0
}}}
{{{#!td
Number concentration of the first log-normal distribution steering the initial dry aerosol spectrum. See [#init_aerosol_probabilistic init_aerosol_probabilistic] for more details. 

'''n1''' can be given in arbitrary units, since the final number concentration is still steered via [#initial_weighting_factor initial_weighting_factor].
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#n2 '''n2''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
0.0
}}}
{{{#!td
Number concentration of the second log-normal distribution steering the initial dry aerosol spectrum. See [#n1 n1].
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#n3 '''n3''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
0.0
}}}
{{{#!td
Number concentration of the third log-normal distribution steering the initial dry aerosol spectrum. See [#n1 n1].
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#particles_per_point '''particles_per_point''']
}}}
{{{#!td style="vertical-align:top"
I
}}}
{{{#!td style="vertical-align:top"
1
}}}
{{{#!td
Number of particles to be started per point.\\
By default, one particle is started at all points of the particle source, defined by the package parameters [#pst pst], [#psl psl], [#psr psr], [#pss pss], [#psn psn], [#psb psb], [#pdx pdx], [#pdy pdy] and [#pdz pdz].
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#particle_advection_start '''particle_advection\\_start''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
0.0
}}}
{{{#!td
Time of the first release of particles (in s).\\
If particles are not to be released at the beginning of the run, the release time can be set via '''particle_advection_start'''.
If particle transport is switched on in a restart run, then [#read_particles_from_restartfile read_particles_from_restartfile] = ''.F.'' is also required. \\
See also [#end_time_prel end_time_prel]. 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#particle_maximum_age '''particle_maximum_age''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
9999999.9
}}}
{{{#!td
Maximum allowed age of particles (in s). \\
If the age of a particle exceeds the time set by '''particle_maximum_age''', the particle is deleted.
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#pdx '''pdx''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
10*[../inipar#dx dx]
}}}
{{{#!td
Distance along x between particles within a particle source (in m). \\
If the particle source shall be confined to one grid point, the distances given by '''pdx''', [#pdy pdy] and [#pdz pdz] must be set larger than the respective domain size or [#psl psl] = [#psr psr] has to be set alternatively.\\
'''pdx''' can be assigned a different value for each particle group (see [#number_of_particle_groups number_of_particle_groups]).
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#pdy '''pdy''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
10*[../inipar#dy dy]
}}}
{{{#!td
Distance along y between particles within a particle source (in m). 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#pdz '''pdz''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
10*\\(zu(2)-zu(1))
}}}
{{{#!td
Distance along z between particles within a particle source (in m). 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#psb '''psb''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
10*zu([../inipar#nz nz]/2)
}}}
{{{#!td
Bottom edge of a particle source (in m).\\
In case of [#seed_follows_topography seed_follows_topography] switched on, the bottom edge height is interpreted as relative to the given topography. 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#psl '''psl''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
10*0.0
}}}
{{{#!td
Left edge of a particle source (in m). 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#psn '''psn''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
10*([#ny ny]*[../inipar#dy dy])
}}}
{{{#!td
Rear ("north") edge of a particle source (in m). 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#psr '''psr''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
10*([#nx nx]*[../inipar#dx dx])
}}}
{{{#!td
Right edge of a particle source (in m). 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#pss '''pss''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
10*0.0
}}}
{{{#!td
Front ("south") edge of a particle source (in m). 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#pst '''pst''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
10*zu([../inipar#nz nz]/2)
}}}
{{{#!td
Top edge of a particle source (in m).\\
In case of [#seed_follows_topography seed_follows_topography] switched on, the top edge height is interpreted as relative to the given topography.  
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#radius '''radius''']
}}}
{{{#!td style="vertical-align:top"
R (10)
}}}
{{{#!td style="vertical-align:top"
0.0, 9*\\9999999.9
}}}
{{{#!td
Particle radius (in m).\\
The viscous friction (in case of a velocity difference between particles and surrounding fluid) depends on the particle radius which must be assigned as soon as [#density_ratio density_ratio] /= ''0.0.'' \\
With several groups of particles (see [#number_of_particle_groups number_of_particle_groups]), each group can be assigned a different value. If the number of values given for '''radius''' is less than the number of groups defined by [#number_of_particle_groups number_of_particle_groups]), then the last assigned value is used for all remaining groups. This means that by default the particle radius for all groups will be ''0.0.''
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#radius_classes '''radius_classes''']
}}}
{{{#!td style="vertical-align:top"
I
}}}
{{{#!td style="vertical-align:top"
20
}}}
{{{#!td
Number of radius classes to be used in the collision efficiency table.\\
This parameter comes into effect if parameter [#collision_kernel collision_kernel] is set to '' 'hall_fast' '' or '' 'wang_fast' ''. It defines the number of radius classes which spawn the collision efficiency table. The interval [1.0E-6,2.0E-4] m is divided into n (= '''radius_classes''') logarithmic equidistant parts. 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#radius_split '''radius_split''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
40.0E-6
}}}
{{{#!td
Radius threshold for splitting of super droplets. Only super droplets with a larger radius than [#radius_split radius_split] are considered for cloning (If  [#splitting splitting]=.T.) . Sensitivity studies showed, that a radius of 40µm seems to be sufficient to achieve better statistics of the rain drop distribution. However, those results are based on simulations of shallow cumulus clouds. Therefore, this parameter may be adjusted for different applications. 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#random_start_position '''random_start_position''']
}}}
{{{#!td style="vertical-align:top"
L
}}}
{{{#!td style="vertical-align:top"
.F.
}}}
{{{#!td
Initial position of the particles is varied randomly within certain limits. \\
By default, the initial positions of particles within the source exactly correspond with the positions given by [#psl psl], [#psr psr], [#psn psn], [#pss pss], [#psb psb], [#pst pst], [#pdx pdx], [#pdy pdy], and [#pdz pdz]. With '''random_start_position''' = ''.T.'' the initial positions of the particles are allowed to randomly vary from these positions within certain limits. \\\\
'''Very important:'''
In case of '''random_start_position''' = ''.T.'', the random-number generators on the individual PEs no longer  run synchronously. If random disturbances are applied to the velocity field (see [../d3par#create_disturbances create_disturbances]), then as a consequence for parallel runs the realizations of the turbulent flow fields will deviate between runs which used different numbers of PEs!
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#read_particles_from_restartfile '''read_particles\\_from_restartfile''']
}}}
{{{#!td style="vertical-align:top"
L
}}}
{{{#!td style="vertical-align:top"
.T.
}}}
{{{#!td
Read particle data from the previous run. \\
By default, with restart runs particle data is read from file [../iofiles#PARTICLE_RESTART_DATA_IN PARTICLE_RESTART_DATA_IN], which is created by the preceding run. If this is not requested or if in a restart run particle transport is switched on for the first time (see [#particle_advection_start particle_advection_start]), then '''read_particles_from_restartfile''' = ''.F.'' is required.
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#rm1 '''rm1''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
0.05E-6
}}}
{{{#!td
Mode radius of the first log-normal distribution steering the initial dry aerosol spectrum. See [#init_aerosol_probabilistic init_aerosol_probabilistic] for more details. 

'''rm1''' should be entered in meters.
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#rm2 '''rm2''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
0.05E-6
}}}
{{{#!td
Mode radius of the second log-normal distribution steering the initial dry aerosol spectrum. See [#rm1 rm1].
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#rm3 '''rm3''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
0.05E-6
}}}
{{{#!td
Mode radius of the third log-normal distribution steering the initial dry aerosol spectrum. See [#rm1 rm1].
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#seed_follows_topography '''seed_follows_\\topography''']
}}}
{{{#!td style="vertical-align:top"
L
}}}
{{{#!td style="vertical-align:top"
.F.
}}}
{{{#!td
Heights of initial particles are interpreted relative to the given topography.\\
In case of topography, the heights of the initial particle sources set by [#psb psb] and [#pst pst] are interpreted as relative to the given topography if [#seed_follows_topography seed_follows_topography]=.TRUE. . Otherwise, if [#seed_follows_topography seed_follows_topography]=.FALSE. and the local particle source is below the surface, no particle are released at this location.
}}}
|----------------
{{{#!td style="vertical-align:top;width: 150px"
[=#splitting '''splitting''']
}}}
{{{#!td style="vertical-align:top;width: 50px"
L
}}}
{{{#!td style="vertical-align:top;width: 75px"
F
}}}
{{{#!td
Switch on/off splitting algorithm for cloud droplets (if [../inipar#cloud_droplets cloud_droplets]=.T.). If [#splitting splitting] is set .TRUE. every time step the splitting algorithm is executed. The algorithm split particles which fulfill certain criterion's into several super droplets with a reduced number of represented particles of every super droplet. The splitting algorithm can be steered by several parameters (a critical radius, a critical weighting factor,..). This mechanism allows an improved representation of the right tail of the drop size distribution with a feasible amount of computational costs. The limits of particle creation should be chosen carefully! The idea of this algorithm is based on Unterstrasser and Soelch, 2014.  
}}}
|----------------
{{{#!td style="vertical-align:top;width: 150px"
[=#splitting_factor '''splitting_factor''']
}}}
{{{#!td style="vertical-align:top;width: 50px"
I
}}}
{{{#!td style="vertical-align:top;width: 75px"
2
}}}
{{{#!td
This parameter is used for [#splitting_mode splitting_mode] = 'const'. It defines the number into how many partiles one particle is divided. This means that, if all other splitting criterion's for one super droplet are fulfilled, splitting_factor-1 new particles (with the same properties as the origin super droplet) are created. The weighting factor of each of the super droplets is divided by the splitting_factor. For [#splitting_mode splitting_mode] = 'cl_av' and  'gb_av' this parameter is calculated automatically and given by the ratio of the actual weighting factor of a super droplet and the approximated number concentration for cloud droplets with a similar radius in a certain volume (grid volume). 
}}}
|----------------
{{{#!td style="vertical-align:top;width: 150px"
[=#splitting_factor_max '''splitting_factor_max''']
}}}
{{{#!td style="vertical-align:top;width: 50px"
I
}}}
{{{#!td style="vertical-align:top;width: 75px"
5
}}}
{{{#!td
This parameter restricts the [#splitting_factor splitting_factor] to an upper bound. It is used for the  [#splitting_mode splitting_mode] 'cl_av_ and 'gb_av' where the [#splitting_factor splitting_factor] is calculated automatically and can get very large values which are not feasible.
}}}
|----------------
{{{#!td style="vertical-align:top;width: 150px"
[=#splitting_function '''splitting_function''']
}}}
{{{#!td style="vertical-align:top;width: 50px"
C
}}}
{{{#!td style="vertical-align:top;width: 75px"
'gamma'
}}}
{{{#!td
Parameter to steer splitting algorithm.\\\\

For the  [#splitting_mode splitting_mode] ='cl_av' and 'gb_av' different functions can be used to approximate the drop size distribution. The user can choose between the following functions:\\\\

'' 'exp' ''
      An exponential distribution.

'' 'gamma' ''
      A gamma size distribution.
'' 'log' ''
      A lognormal distribution.
}}}
|----------------
{{{#!td style="vertical-align:top;width: 150px"
[=#splitting_mode '''splitting_mode''']
}}}
{{{#!td style="vertical-align:top;width: 50px"
C
}}}
{{{#!td style="vertical-align:top;width: 75px"
'const'
}}}
{{{#!td
Parameter to steer splitting algorithm.\\\\

The splitting of cloud droplets can be done in using different modes. The user can choose between the following modes:\\\\

'' 'const' ''
      In splitting mode 'const' a critical radius  [#radius_split radius_split] a critical weighting factor [#weight_factor_split weight_factor_split] and a [#splitting_factor splitting_factor] must  be prescribed. Super droplets which have a larger radius and larger weighting factor are split into [#splitting_factor splitting_factor] super droplets. Therefore, the weighting factor of the super droplet and all created clones is reduced by the factor of [#splitting_factor splitting_factor].

'' 'cl_av' ''
      In splitting mode 'cl_av'  a critical radius [#radius_split radius_split] and a splitting function must be prescribed. The critical weighting factor is calculated while approximating a 'gamma', 'log' or 'exp'- drop size distribution. In this mode a drop size distribution (discretized in 100 bins in a range from 1.0 µm to 1 mm) is calculated as an average over all cloudy grid boxes. Super droplets which have a larger radius and larger weighting factor are split into 'splitting_factor' super droplets. In this case the [#splitting_factor splitting_factor] is calculated of weighting factor of the super droplet and the approximated number concentration for droplet of such a size. Due to the splitting, the weighting factor of the super droplet and all created clones is reduced by the factor of 'splitting_facor'.

'' 'gb_av' ''
      Same as for 'cl_av' but a drop size distribution is calculated for every grid box.

}}}
|----------------
{{{#!td style="vertical-align:top;width: 150px"
[=#step_dealloc '''step_dealloc''']
}}}
{{{#!td style="vertical-align:top;width: 50px"
I
}}}
{{{#!td style="vertical-align:top;width: 75px"
100
}}}
{{{#!td
Number of timesteps after which particle memory is deallocated. The parameter has only an effect if [#dealloc_memory dealloc_memory] = .True. .  
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#s1 '''s1''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
2.0
}}}
{{{#!td
Geometric standard deviation of the first log-normal distribution steering the initial dry aerosol spectrum. See [#init_aerosol_probabilistic init_aerosol_probabilistic] for more details. 

'''s1''' has no units.
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#s2 '''s2''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
2.0
}}}
{{{#!td
Geometric standard deviation of the second log-normal distribution steering the initial dry aerosol spectrum. See [#s1 s1].
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#s3 '''s3''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
2.0
}}}
{{{#!td
Geometric standard deviation of the third log-normal distribution steering the initial dry aerosol spectrum. See [#s1 s1].
}}} 
|----------------
{{{#!td style="vertical-align:top"
[=#use_sgs_for_particles '''use_sgs_for_particles''']
}}}
{{{#!td style="vertical-align:top"
L
}}}
{{{#!td style="vertical-align:top"
.F.
}}}
{{{#!td
Use subgrid-scale velocities for particle advection.\\
These velocities are calculated from the resolved and subgrid-scale TKE using the Monte-Carlo random-walk method described by Weil et al. (2004, JAS, 61, 2877-2887). When using this method, the time step for the advancement of the particles is limited by the so-called Lagrangian time scale. This may be smaller than the current LES time step so that several particle (sub-) time steps have to be carried out within one LES time step. In order to limit the number of sub-time steps (and to limit the CPU-time), the minimum value for the particle time step is defined by the package parameter [#dt_min_part dt_min_part].
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#vertical_particle_advection '''vertical_particle\\_advection''']
}}}
{{{#!td style="vertical-align:top"
L
}}}
{{{#!td style="vertical-align:top"
.T.
}}}
{{{#!td
Switch on/off vertical particle transport.\\
By default, particles are transported along all three directions in space. With '''vertical_particle_advection''' = ''.F.'', the particles will only be transported horizontally.
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#weight_factor_split '''weight_factor_split''']
}}}
{{{#!td style="vertical-align:top"
R
}}}
{{{#!td style="vertical-align:top"
0.1 * [#initial_weighting_factor initial_weighting_factor]
}}}
{{{#!td
For [#splitting_mode splitting_mode]='const' a critical weighting factor must be prescribed. Super droplets with a larger weighting factor are considered for cloning. 
}}}
|----------------
{{{#!td style="vertical-align:top"
[=#write_particle_statistics '''write_particle_statistics''']
}}}
{{{#!td style="vertical-align:top"
L
}}}
{{{#!td style="vertical-align:top"
.F.
}}}
{{{#!td
Switch on/off output of particle information.\\
For '''write_particle_statistics''' = ''.T.'' statistical information (e.g. the total number of particles used, the number of particles exchanged between the PEs, etc.) which may be used for debugging are output to the local file [../iofiles#PARTICLE_INFOS PARTICLE_INFOS]. \\\\
'''Note:''' For parallel runs files may become very large and performance of PALM may decrease.
}}}