Changeset 2265 for palm/trunk/SOURCE/lpm_init.f90
- Timestamp:
- Jun 8, 2017 4:58:28 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/lpm_init.f90
r2263 r2265 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Unused variables removed. 28 ! 29 ! 2263 2017-06-08 14:59:01Z schwenkel 27 30 ! Implemented splitting and merging algorithm 28 31 ! … … 177 180 ONLY: alloc_factor, bc_par_b, bc_par_lr, bc_par_ns, bc_par_t, & 178 181 block_offset, block_offset_def, collision_kernel, & 179 curvature_solution_effects, & 180 density_ratio, grid_particles, isf,i_splitting_mode, & 181 initial_weighting_factor, ibc_par_b, ibc_par_lr, ibc_par_ns, & 182 ibc_par_t, iran_part, log_z_z0, & 183 max_number_of_particle_groups, maximum_number_of_particles, & 184 min_nr_particle, mpi_particle_type, & 185 number_concentration, number_particles_per_gridbox, & 186 number_of_particles, & 182 curvature_solution_effects, density_ratio, grid_particles, & 183 isf,i_splitting_mode, initial_weighting_factor, ibc_par_b, & 184 ibc_par_lr, ibc_par_ns, ibc_par_t, iran_part, log_z_z0, & 185 max_number_of_particle_groups, min_nr_particle, & 186 mpi_particle_type, number_concentration, & 187 number_particles_per_gridbox, number_of_particles, & 187 188 number_of_particle_groups, number_of_sublayers, & 188 189 offset_ocean_nzt, offset_ocean_nzt_m1, & 189 190 particles, particle_advection_start, particle_groups, & 190 191 particle_groups_type, particles_per_point, & 191 particle_type, pdx, pdy, pdz, & 192 prt_count, psb, psl, psn, psr, pss, pst, & 193 radius, random_start_position, read_particles_from_restartfile,& 194 seed_follows_topography, sgs_wf_part, sort_count, & 195 splitting_function, splitting_mode, total_number_of_particles, & 196 use_sgs_for_particles, & 197 write_particle_statistics, uniform_particles, zero_particle, & 198 z0_av_global 192 particle_type, pdx, pdy, pdz, prt_count, psb, psl, psn, psr, & 193 pss, pst, radius, random_start_position, & 194 read_particles_from_restartfile, seed_follows_topography, & 195 sgs_wf_part, sort_count, splitting_function, splitting_mode, & 196 total_number_of_particles, use_sgs_for_particles, & 197 write_particle_statistics, zero_particle, z0_av_global 199 198 200 199 USE pegrid … … 532 531 grid_particles(nzb+1:nzt,nys:nyn,nxl:nxr) ) 533 532 534 maximum_number_of_particles = 0535 533 number_of_particles = 0 536 534 … … 592 590 CALL check_open( 80 ) 593 591 WRITE ( 80, 8000 ) current_timestep_number, simulated_time, & 594 number_of_particles, & 595 maximum_number_of_particles 592 number_of_particles 596 593 CALL close_file( 80 ) 597 594 ENDIF … … 956 953 957 954 ! 958 !-- Determine maximum number of particles (i.e., all possible particles that 959 !-- have been allocated) and the current number of particles 955 !-- Determine the current number of particles 960 956 DO ip = nxl, nxr 961 957 DO jp = nys, nyn 962 958 DO kp = nzb+1, nzt 963 maximum_number_of_particles = maximum_number_of_particles &964 + SIZE(grid_particles(kp,jp,ip)%particles)965 959 number_of_particles = number_of_particles & 966 960 + prt_count(kp,jp,ip)
Note: See TracChangeset
for help on using the changeset viewer.