*
Initializes particles equally distributed over the first release location * specified in file RELEASES. This box is assumed to be the domain for doing * domain-filling trajectory calculations. * All particles carry the same amount of mass which alltogether comprises the* mass of air within the box. * * Author: A. Stohl * * 15 October 2002 * * Changes * 2022, L. Bakels: OpenMP parallelisation * 2023, L. Bakels: smooth vertical particle distribution instead of * distributing particles on fixed vertical layers * *
*
Variables: * * numparticlecount consecutively counts the number of particles released * nx_we(2) grid indices for western and eastern boundary of domain- * filling trajectory calculations * ny_sn(2) grid indices for southern and northern boundary of domain- * filling trajectory calculations * *
First for the south pole
$OMP DO
Do the following steps only if particles are not read in from previous model run
!$OMP DO