Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/lpm_read_restart_file.f90

    r3049 r3560  
    2525! -----------------
    2626! $Id$
     27! time_prel replaced by last_particle_release_time
     28!
     29! 3049 2018-05-29 13:52:36Z Giersch
    2730! Error messages revised
    2831!
     
    9396    USE kinds
    9497
    95     USE lpm_pack_and_sort_mod,                                                   &
     98    USE lpm_pack_and_sort_mod,                                                 &
    9699        ONLY:  lpm_sort_in_subboxes
    97100
    98101    USE particle_attributes,                                                   &
    99102        ONLY:  alloc_factor, bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,         &
    100                grid_particles, min_nr_particle, number_of_particles,           &
    101                number_of_particle_groups, particle_groups, particle_type,      &
    102                prt_count, time_prel, time_write_particle_data, zero_particle
     103               grid_particles, last_particle_release_time, min_nr_particle,    &
     104               number_of_particles, number_of_particle_groups, particle_groups,&
     105               particle_type, prt_count, time_write_particle_data, zero_particle
    103106
    104107    USE pegrid
     
    153156!-- allocate them and read their contents.
    154157    READ ( 90 )  bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,                     &
    155                  number_of_particle_groups, particle_groups, time_prel,        &
    156                  time_write_particle_data
     158                 last_particle_release_time, number_of_particle_groups,        &
     159                 particle_groups, time_write_particle_data
    157160
    158161    ALLOCATE( prt_count(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                        &
Note: See TracChangeset for help on using the changeset viewer.