source: palm/trunk/SOURCE/package_parin.f90 @ 828

Last change on this file since 828 was 828, checked in by raasch, 10 years ago

New:
---

Changed:


Optimization of collision kernels. Collision tables can be calculated once at
simulation start for defined radius (and dissipation) classes instead of
re-calculating them at every timestep and for the particle ensemble in
every gridbox.
For this purpose the particle feature color is renamed class.
New parpar parameters radius_classes and dissipation_classes.
(Makefile, advec_particles, check_parameters, data_output_dvrp, header, init_particles, lpm_collision_kernels, modules, package_parin, set_particle_attributes)

Lower limit for droplet radius changed from 1E-7 to 1E-8.
(advec_particles)

Complete re-formatting of collision code (including changes in names of
variables, modules and subroutines).
(advec_particles, lpm_collision_kernels)

Errors:


Bugfix: transformation factor for dissipation changed from 1E5 to 1E4
(advec_particles, lpm_collision_kernels)

  • Property svn:keywords set to Id
File size: 6.1 KB
Line 
1 SUBROUTINE package_parin
2
3!------------------------------------------------------------------------------!
4! Current revisions:
5! -----------------
6! +dissipation_classes, radius_classes in parpar
7!
8! Former revisions:
9! -----------------
10! $Id: package_parin.f90 828 2012-02-21 12:00:36Z raasch $
11!
12! 825 2012-02-19 03:03:44Z raasch
13! wang_collision_kernel and turbulence_effects_on_collision in particles_par
14! replaced by collision_kernel
15!
16! 790 2011-11-29 03:11:20Z raasch
17! +turbulence_effects_on_collision, wang_collision_kernel in particles_par
18!
19! 336 2009-06-10 11:19:35Z raasch
20! +clip_dvrp_*, cluster_size, color_interval, dvrpsize_interval,
21! groundplate_color, isosurface_color, particle_color, particle_dvrpsize
22! topography_color, in dvrp_graphics_par,
23! parameter dvrp_psize moved from particles_par to dvrp_graphics_par
24! Variables for dvrp-mode pathlines added
25!
26! 210 2008-11-06 08:54:02Z raasch
27! Variables for dvrp-mode pathlines added
28!
29! 116 2007-10-11 02:30:27Z raasch
30! +dt_sort_particles in package_parin
31!
32! 60 2007-03-11 11:50:04Z raasch
33! Particles-package is now part of the default code
34!
35! RCS Log replace by Id keyword, revision history cleaned up
36!
37! Revision 1.18  2006/08/04 14:52:23  raasch
38! +dt_dopts, dt_min_part, end_time_prel, particles_per_point,
39! use_sgs_for_particles in particles_par
40!
41! Revision 1.1  2000/12/28 13:21:57  raasch
42! Initial revision
43!
44!
45! Description:
46! ------------
47! This subroutine reads from the NAMELIST file variables controling model
48! software packages which are used optionally in the run.
49!------------------------------------------------------------------------------!
50
51    USE control_parameters
52    USE dvrp_variables
53    USE particle_attributes
54    USE spectrum
55
56    IMPLICIT NONE
57
58    CHARACTER (LEN=80) ::  zeile
59
60    NAMELIST /dvrp_graphics_par/  clip_dvrp_l, clip_dvrp_n, clip_dvrp_r,       &
61                                  clip_dvrp_s, cluster_size, color_interval,   &
62                                  dt_dvrp, dvrpsize_interval, dvrp_directory,  &
63                                  dvrp_file, dvrp_host, dvrp_output,           &
64                                  dvrp_password, dvrp_psize, dvrp_username,    &
65                                  groundplate_color, isosurface_color,         &
66                                  mode_dvrp, particle_color, particle_dvrpsize,&
67                                  pathlines_fadeintime, pathlines_fadeouttime, &
68                                  pathlines_linecount, pathlines_maxhistory,   &
69                                  pathlines_wavecount, pathlines_wavetime,     &
70                                  slicer_range_limits_dvrp, superelevation,    &
71                                  superelevation_x, superelevation_y,          &
72                                  threshold, topography_color, vc_alpha,       &
73                                  vc_gradient_normals, vc_mode, vc_size_x,     &
74                                  vc_size_y, vc_size_z
75
76    NAMELIST /particles_par/      bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,    &
77                                  collision_kernel, density_ratio,             &
78                                  dissipation_classes, dt_dopts,               &
79                                  dt_min_part, dt_prel, dt_sort_particles,     &
80                                  dt_write_particle_data,                      &
81                                  end_time_prel, initial_weighting_factor,     &
82                                  maximum_number_of_particles,                 &
83                                  maximum_number_of_tailpoints,                &
84                                  maximum_tailpoint_age,                       &
85                                  minimum_tailpoint_distance,                  &
86                                  number_of_particle_groups,                   &
87                                  particles_per_point,                         &
88                                  particle_advection_start,                    &
89                                  particle_maximum_age, pdx, pdy, pdz, psb,    &
90                                  psl, psn, psr, pss, pst, radius,             &
91                                  radius_classes, random_start_position,       &
92                                  read_particles_from_restartfile,             &
93                                  skip_particles_for_tail,                     &
94                                  use_particle_tails, use_sgs_for_particles,   &
95                                  vertical_particle_advection,                 &
96                                  write_particle_statistics
97    NAMELIST /spectra_par/        averaging_interval_sp, comp_spectra_level,   &
98                                  data_output_sp, dt_dosp, plot_spectra_level, &
99                                  skip_time_dosp, spectra_direction
100
101!
102!-- Position the namelist-file at the beginning (it was already opened in
103!-- parin), search for the namelist-group of the package and position the
104!-- file at this line. Do the same for each optionally used package.
105    zeile = ' '
106
107#if defined( __dvrp_graphics )
108    REWIND ( 11 )
109    zeile = ' '
110    DO   WHILE ( INDEX( zeile, '&dvrp_graphics_par' ) == 0 )
111       READ ( 11, '(A)', END=10 )  zeile
112    ENDDO
113    BACKSPACE ( 11 )
114
115!
116!-- Read user-defined namelist
117    READ ( 11, dvrp_graphics_par )
118
119 10 CONTINUE
120#endif
121
122!
123!-- Try to find particles package
124    REWIND ( 11 )
125    zeile = ' '
126    DO   WHILE ( INDEX( zeile, '&particles_par' ) == 0 )
127       READ ( 11, '(A)', END=20 )  zeile
128    ENDDO
129    BACKSPACE ( 11 )
130
131!
132!-- Read user-defined namelist
133    READ ( 11, particles_par )
134
135!
136!-- Set flag that indicates that particles are switched on
137    particle_advection = .TRUE.
138
139 20 CONTINUE
140
141
142#if defined( __spectra )
143    REWIND ( 11 )
144    zeile = ' '
145    DO   WHILE ( INDEX( zeile, '&spectra_par' ) == 0 )
146       READ ( 11, '(A)', END=30 )  zeile
147    ENDDO
148    BACKSPACE ( 11 )
149
150!
151!-- Read user-defined namelist
152    READ ( 11, spectra_par )
153
154!
155!-- Default setting of dt_dosp here (instead of check_parameters), because its
156!-- current value is needed in init_pegrid
157    IF ( dt_dosp == 9999999.9 )  dt_dosp = dt_data_output
158
159 30 CONTINUE
160#endif
161
162 END SUBROUTINE package_parin
Note: See TracBrowser for help on using the repository browser.