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

Last change on this file since 2287 was 2263, checked in by schwenkel, 7 years ago

Implemented splitting and merging algorithm

  • Property svn:keywords set to Id
File size: 11.7 KB
Line 
1!> @file package_parin.f90
2!------------------------------------------------------------------------------!
3! This file is part of PALM.
4!
5! PALM is free software: you can redistribute it and/or modify it under the
6! terms of the GNU General Public License as published by the Free Software
7! Foundation, either version 3 of the License, or (at your option) any later
8! version.
9!
10! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
13!
14! You should have received a copy of the GNU General Public License along with
15! PALM. If not, see <http://www.gnu.org/licenses/>.
16!
17! Copyright 1997-2017 Leibniz Universitaet Hannover
18!------------------------------------------------------------------------------!
19!
20! Current revisions:
21! -----------------
22!
23!
24! Former revisions:
25! -----------------
26! $Id: package_parin.f90 2263 2017-06-08 14:59:01Z suehring $
27! Implemented splitting and merging algorithm
28!
29! 2183 2017-03-17 14:29:15Z schwenkel
30!
31! 2182 2017-03-17 14:27:40Z schwenkel
32! Added parameters for simplified particle initialization.
33!
34! 2000 2016-08-20 18:09:15Z knoop
35! Forced header and separation lines into 80 columns
36!
37! 1936 2016-06-13 13:37:44Z suehring
38! +deallocate_memory, step_dealloc
39!
40! 1871 2016-04-15 11:46:09Z hoffmann
41! Initialization of aerosols added.
42!
43! 1833 2016-04-07 14:23:03Z raasch
44! reading of spectra_par moved to spectra_mod
45!
46! 1831 2016-04-07 13:15:51Z hoffmann
47! curvature_solution_effects added
48!
49! 1826 2016-04-07 12:01:39Z maronga
50! Reading of &radiation_par moved to radiation_model_mod.
51! Reading of &canopy_par moved to plant_canopy_model_mod.
52!
53! 822 2016-04-07 07:49:42Z hoffmann
54! +collision_algorithm
55! Tails removed.
56!
57! 1817 2016-04-06 15:44:20Z maronga
58! Reading of &lsm_par moved to land_surface_model_mod.
59!
60! 1788 2016-03-10 11:01:04Z maronga
61! Parameter dewfall removed.
62!
63! 1786 2016-03-08 05:49:27Z raasch
64! cpp-direktives for spectra removed
65!
66! 1757 2016-02-22 15:49:32Z maronga
67! Added parameter unscheduled_radiation_calls
68!
69! 1691 2015-10-26 16:17:44Z maronga
70! Added skip_time_do_lsm, skip_time_do_radiation, and emissivity
71!
72! 1682 2015-10-07 23:56:08Z knoop
73! Code annotations made doxygen readable
74!
75! 1585 2015-04-30 07:05:52Z maronga
76! Added several radiation_par parameters
77!
78! 1575 2015-03-27 09:56:27Z raasch
79! +seed_follows_topography in particles_par
80!
81! 1553 2015-03-03 17:33:54Z maronga
82! Resorting of lsm_par
83!
84! 1551 2015-03-03 14:18:16Z maronga
85! Several changes in the liste for land surface model and radiation model
86!
87! 1496 2014-12-02 17:25:50Z maronga
88! Added support for the land surface model and radiation scheme
89!
90! 1484 2014-10-21 10:53:05Z kanani
91! Changes due to new module structure of the plant canopy model:
92!   module plant_canopy_model_mod added,
93!   new package/namelist canopy_par added, i.e. the canopy model is no longer
94!   steered over the inipar-namelist,
95!   drag_coefficient, leaf_surface_concentration and scalar_exchange_coefficient
96!   renamed to canopy_drag_coeff, leaf_surface_conc and leaf_scalar_exch_coeff.
97! Changed statement tags in CONTINUE-statement
98!
99! 1367 2014-04-23 15:18:30Z witha
100! Bugfix: module kinds must be used
101!
102! 1359 2014-04-11 17:15:14Z hoffmann
103! +alloc_factor, + min_nr_particle
104! -dt_sort_particles, -maximum_number_of_particles
105!
106! 1340 2014-03-25 19:45:13Z kanani
107! REAL constants defined as wp-kinds
108!
109! 1324 2014-03-21 09:13:16Z suehring
110! Bugfix: Missing variable dt_data_output output added to ONLY statement
111!
112! 1320 2014-03-20 08:40:49Z raasch
113! ONLY-attribute added to USE-statements,
114! kind-parameters added to all INTEGER and REAL declaration statements,
115! kinds are defined in new module kinds,
116! old module precision_kind is removed,
117! revision history before 2012 removed,
118! comment fields (!:) to be used for variable explanations added to
119! all variable declaration statements
120!
121! 1036 2012-10-22 13:43:42Z raasch
122! code put under GPL (PALM 3.9)
123!
124! 828 2012-02-21 12:00:36Z raasch
125! +dissipation_classes, radius_classes in parpar
126!
127! 825 2012-02-19 03:03:44Z raasch
128! wang_collision_kernel and turbulence_effects_on_collision in particles_par
129! replaced by collision_kernel
130!
131! Revision 1.1  2000/12/28 13:21:57  raasch
132! Initial revision
133!
134!
135! Description:
136! ------------
137!> This subroutine reads from the NAMELIST file variables controling model
138!> software packages which are used optionally in the run.
139!>
140!> @todo Perform all actions in the respective submodules and remove
141!>       package_parin
142!------------------------------------------------------------------------------!
143 SUBROUTINE package_parin
144 
145
146    USE control_parameters,                                                    &
147        ONLY:  dt_data_output, dt_dopts, dt_dvrp, particle_maximum_age,        &
148               threshold
149
150    USE dvrp_variables,                                                        &
151        ONLY:  clip_dvrp_l, clip_dvrp_n, clip_dvrp_r, clip_dvrp_s,             &
152               cluster_size, color_interval, dvrpsize_interval,                &
153               dvrp_directory, dvrp_file, dvrp_host, dvrp_output,              &
154               dvrp_password, dvrp_username, groundplate_color,                &
155               isosurface_color, mode_dvrp, particle_color,                    &
156               particle_dvrpsize, pathlines_fadeintime,                        &
157               pathlines_fadeouttime, pathlines_linecount,                     &
158               pathlines_maxhistory, pathlines_wavecount,                      &
159               pathlines_wavetime, slicer_range_limits_dvrp, superelevation,   &
160               superelevation_x, superelevation_y, topography_color,           &
161               vc_alpha, vc_gradient_normals, vc_mode, vc_size_x, vc_size_y,   &
162               vc_size_z
163
164    USE kinds
165
166    USE particle_attributes,                                                   &
167        ONLY:  alloc_factor, bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,         &
168               collision_algorithm, collision_kernel,                          &
169               curvature_solution_effects, deallocate_memory, density_ratio,   &
170               dissipation_classes, dt_min_part, dt_prel,                      &
171               dt_write_particle_data, end_time_prel, initial_weighting_factor,&
172               init_aerosol_probabilistic, max_number_particles_per_gridbox,   & 
173               merging, min_nr_particle, monodisperse_aerosols, n1, n2, n3,    &
174               number_concentration, number_of_particle_groups,                &
175               number_particles_per_gridbox, particles_per_point,              &
176               particle_advection, particle_advection_start, pdx, pdy, pdz,    &
177               psb, psl, psn, psr, pss, pst, radius, radius_classes,           &
178               radius_merge, radius_split, random_start_position,              &
179               read_particles_from_restartfile, rm1, rm2, rm3,                 &
180               seed_follows_topography, splitting, splitting_factor,           & 
181               splitting_factor_max, splitting_function, splitting_mode,       &
182               step_dealloc, s1, s2, s3, use_sgs_for_particles,                &
183               vertical_particle_advection, weight_factor_merge,               &
184               weight_factor_split, write_particle_statistics               
185
186    IMPLICIT NONE
187
188    CHARACTER (LEN=80) ::  line  !<
189
190
191
192    NAMELIST /dvrp_graphics_par/  clip_dvrp_l, clip_dvrp_n, clip_dvrp_r,       &
193                                  clip_dvrp_s, cluster_size, color_interval,   &
194                                  dt_dvrp, dvrpsize_interval, dvrp_directory,  &
195                                  dvrp_file, dvrp_host, dvrp_output,           &
196                                  dvrp_password, dvrp_username,                &
197                                  groundplate_color, isosurface_color,         &
198                                  mode_dvrp, particle_color, particle_dvrpsize,&
199                                  pathlines_fadeintime, pathlines_fadeouttime, &
200                                  pathlines_linecount, pathlines_maxhistory,   &
201                                  pathlines_wavecount, pathlines_wavetime,     &
202                                  slicer_range_limits_dvrp, superelevation,    &
203                                  superelevation_x, superelevation_y,          &
204                                  threshold, topography_color, vc_alpha,       &
205                                  vc_gradient_normals, vc_mode, vc_size_x,     &
206                                  vc_size_y, vc_size_z
207
208    NAMELIST /particles_par/      alloc_factor, bc_par_b, bc_par_lr,           &
209                                  bc_par_ns, bc_par_t, collision_algorithm,    &
210                                  collision_kernel, curvature_solution_effects,&
211                                  deallocate_memory, density_ratio,            &
212                                  dissipation_classes, dt_dopts,               &
213                                  dt_min_part, dt_prel,                        &
214                                  dt_write_particle_data,                      &
215                                  end_time_prel, initial_weighting_factor,     &
216                                  init_aerosol_probabilistic,                  &
217                                  max_number_particles_per_gridbox, merging,   &
218                                  min_nr_particle, monodisperse_aerosols,      &
219                                  n1, n2, n3, number_concentration,            &
220                                  number_of_particle_groups,                   &
221                                  number_particles_per_gridbox,                &
222                                  particles_per_point,                         &
223                                  particle_advection_start,                    &
224                                  particle_maximum_age, pdx, pdy, pdz, psb,    &
225                                  psl, psn, psr, pss, pst, radius,             &
226                                  radius_classes, radius_merge, radius_split,  &   
227                                  random_start_position,                       &
228                                  read_particles_from_restartfile,             &
229                                  rm1, rm2, rm3,                               &
230                                  seed_follows_topography,                     &
231                                  splitting, splitting_factor,                 &
232                                  splitting_factor_max, splitting_function,    &
233                                  splitting_mode, step_dealloc, s1, s2, s3,    &
234                                  use_sgs_for_particles,                       &
235                                  vertical_particle_advection,                 &
236                                  weight_factor_merge, weight_factor_split,    &
237                                  write_particle_statistics
238
239!
240!-- Position the namelist-file at the beginning (it was already opened in
241!-- parin), search for the namelist-group of the package and position the
242!-- file at this line. Do the same for each optionally used package.
243    line = ' '
244
245
246#if defined( __dvrp_graphics )
247    REWIND ( 11 )
248    line = ' '
249    DO   WHILE ( INDEX( line, '&dvrp_graphics_par' ) == 0 )
250       READ ( 11, '(A)', END=20 )  line
251    ENDDO
252    BACKSPACE ( 11 )
253
254!
255!-- Read user-defined namelist
256    READ ( 11, dvrp_graphics_par )
257
258 20 CONTINUE
259#endif
260
261!
262!-- Try to find particles package
263    REWIND ( 11 )
264    line = ' '
265    DO   WHILE ( INDEX( line, '&particles_par' ) == 0 )
266       READ ( 11, '(A)', END=30 )  line
267    ENDDO
268    BACKSPACE ( 11 )
269
270!
271!-- Read user-defined namelist
272    READ ( 11, particles_par )
273
274!
275!-- Set flag that indicates that particles are switched on
276    particle_advection = .TRUE.
277
278 30 CONTINUE
279
280 END SUBROUTINE package_parin
Note: See TracBrowser for help on using the repository browser.