Changeset 3560 for palm/trunk/SOURCE/lpm_init.f90

Timestamp:

Nov 23, 2018 9:20:21 AM (5 years ago)

Author:

raasch

Message:

bugfix to guarantee correct particle releases in case that the release interval is smaller than the model timestep

File:

• palm/trunk/SOURCE/lpm_init.f90 (modified) (3 diffs)

palm/trunk/SOURCE/lpm_init.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! set the first particle release time
+! 
+! 3524 2018-11-14 13:36:44Z raasch
 ! added missing working precision
 ! 
@@ -260 +263 @@
                 curvature_solution_effects, density_ratio, grid_particles,     &
                 isf,i_splitting_mode, initial_weighting_factor, ibc_par_b,     &
-                ibc_par_lr, ibc_par_ns, ibc_par_t, iran_part, log_z_z0,        &
+                ibc_par_lr, ibc_par_ns, ibc_par_t, iran_part,                  &
+                last_particle_release_time, log_z_z0,                          &
                 max_number_of_particle_groups, min_nr_particle,                &
                 number_concentration,                                          &
@@ -654 +658 @@
        iran_part = iran_part + myid
 
-       CALL lpm_create_particle (PHASE_INIT)
+!
+!--    Create the particle set, and set the initial particles
+       CALL lpm_create_particle( phase_init )
+       last_particle_release_time = particle_advection_start
+
 !
 !--    User modification of initial particles