source: palm/trunk/SOURCE/parin.f90 @ 2292

Last change on this file since 2292 was 2292, checked in by schwenkel, 4 years ago

implementation of new bulk microphysics scheme

  • Property svn:keywords set to Id
File size: 27.8 KB
Line 
1!> @file 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: parin.f90 2292 2017-06-20 09:51:42Z schwenkel $
27! Implementation of new microphysic scheme: cloud_scheme = 'morrison'
28! includes two more prognostic equations for cloud drop concentration (nc) 
29! and cloud water content (qc).
30!
31! 2267 2017-06-09 09:33:25Z gronemeier
32! Bugfix: removed skipping of reading namelists in case of omitted d3par
33!
34! 2259 2017-06-08 09:09:11Z gronemeier
35! Implemented synthetic turbulence generator
36!
37! 2233 2017-05-30 18:08:54Z suehring
38!
39! 2232 2017-05-30 17:47:52Z suehring
40! typo corrected
41! +wall_salinityflux
42! +tunnel_height, tunnel_lenght, tunnel_width_x, tunnel_width_y,
43!  tunnel_wall_depth
44!
45! 2118 2017-01-17 16:38:49Z raasch
46! -background_communication from inipar
47!
48! 2050 2016-11-08 15:00:55Z gronemeier
49! Implement turbulent outflow condition
50!
51! 2037 2016-10-26 11:15:40Z knoop
52! Anelastic approximation implemented
53!
54! 2035 2016-10-24 15:06:17Z suehring
55! Remove check for npex and npey in nesting case
56!
57! 2011 2016-09-19 17:29:57Z kanani
58! Added flag lsf_exception to allow explicit enabling of large scale forcing
59! together with buildings on flat terrain.
60!
61! 2007 2016-08-24 15:47:17Z kanani
62! Added call to urban surface model for reading of &urban_surface_par
63!
64! 2004 2016-08-24 10:25:59Z suehring
65! Humidity and passive scalar treated separately in nesting mode
66!
67! 2000 2016-08-20 18:09:15Z knoop
68! Forced header and separation lines into 80 columns
69!
70! 1992 2016-08-12 15:14:59Z suehring
71! +top_scalarflux
72!
73! 1960 2016-07-12 16:34:24Z suehring
74! Allocation of s_init
75!
76! 1957 2016-07-07 10:43:48Z suehring
77! flight module added
78!
79! 1955 2016-07-01 12:38:59Z hellstea
80! The parameter intializating_actions is set to 'set_constant_profiles for
81! all nest domains in order to make sure that diagnostic variables are properly
82! initialized for nest domains. Prognostic variables are later initialized by
83! interpolation from the parent domain.
84!
85! 1917 2016-05-27 14:28:12Z witha
86! Initial version of purely vertical nesting introduced.
87!
88! 1914 2016-05-26 14:44:07Z witha
89! Added call to wind turbine model for reading of &wind_turbine_par
90!
91! 1849 2016-04-08 11:33:18Z hoffmann
92! Adapted for modularization of microphysics
93!
94! 1833 2016-04-07 14:23:03Z raasch
95! call of spectra_parin
96!
97! 1831 2016-04-07 13:15:51Z hoffmann
98! turbulence renamed collision_turbulence, drizzle renamed
99! cloud_water_sedimentation
100! curvature_solution_effects removed
101!
102! 1826 2016-04-07 12:01:39Z maronga
103! Added call to radiation model for reading of &radiation_par.
104! Added call to plant canopy model for reading of &canopy_par.
105!
106! 1817 2016-04-06 15:44:20Z maronga
107! Added call to land surface model for reading of &lsm_par
108!
109! 1804 2016-04-05 16:30:18Z maronga
110! Removed code for parameter file check (__check)
111!
112! 1783 2016-03-06 18:36:17Z raasch
113! +netcdf_deflate in d3par, netcdf module and variable names changed
114!
115! 1764 2016-02-28 12:45:19Z raasch
116! cpp-statements for nesting removed, explicit settings of boundary conditions
117! in nest domains,
118! bugfix: npex/npey message moved from inipar to d3par
119! bugfix: check of lateral boundary conditions from check_parameters to here,
120! because they will be already used in init_pegrid and init_grid
121!
122! 1762 2016-02-25 12:31:13Z hellstea
123! Introduction of nested domain feature
124!
125! 1691 2015-10-26 16:17:44Z maronga
126! Added parameter most_method. Renamed prandtl_layer to constant_flux_layer.
127!
128! 1682 2015-10-07 23:56:08Z knoop
129! Code annotations made doxygen readable
130!
131! 1560 2015-03-06 10:48:54Z keck
132! +recycling_yshift
133!
134! 1496 2014-12-02 17:25:50Z maronga
135! Renamed: "radiation -> "cloud_top_radiation"
136!
137! 1484 2014-10-21 10:53:05Z kanani
138! Changes due to new module structure of the plant canopy model:
139!   canopy-model related parameters moved to new package canopy_par in
140!   subroutine package_parin
141!
142! 1429 2014-07-15 12:53:45Z knoop
143! +ensemble_member_nr to prepare the random_generator for ensemble runs
144!
145! 1402 2014-05-09 14:25:13Z raasch
146! location messages modified, batch_job included in envpar-NAMELIST
147!
148! 1384 2014-05-02 14:31:06Z raasch
149! location messages added
150!
151! 1365 2014-04-22 15:03:56Z boeske
152! Usage of large scale forcing enabled:
153! +use_subsidence_tendencies
154!
155! 1361 2014-04-16 15:17:48Z hoffmann
156! +call_microphysics_at_all_substeps
157!
158! 1359 2014-04-11 17:15:14Z hoffmann
159! REAL constants provided with KIND-attribute
160!
161! 1353 2014-04-08 15:21:23Z heinze
162! REAL constants provided with KIND-attribute
163!
164! 1327 2014-03-21 11:00:16Z raasch
165! -data_output_format, do3d_compress, do3d_comp_prec
166!
167! 1320 2014-03-20 08:40:49Z raasch
168! ONLY-attribute added to USE-statements,
169! kind-parameters added to all INTEGER and REAL declaration statements,
170! kinds are defined in new module kinds,
171! old module precision_kind is removed,
172! revision history before 2012 removed,
173! comment fields (!:) to be used for variable explanations added to
174! all variable declaration statements
175!
176! 1318 2014-03-17 13:35:16Z raasch
177! +cpu_log_barrierwait in d3par
178!
179! 1301 2014-03-06 13:29:46Z heinze
180! +large_scale_subsidence
181!
182! 1241 2013-10-30 11:36:58Z heinze
183! +nudging
184! +large_scale_forcing
185!
186! 1216 2013-08-26 09:31:42Z raasch
187! +transpose_compute_overlap in inipar
188!
189! 1195 2013-07-01 12:27:57Z heinze
190! Bugfix: allocate ref_state
191!
192! 1179 2013-06-14 05:57:58Z raasch
193! +reference_state in inipar
194!
195! 1159 2013-05-21 11:58:22Z fricke
196! +use_cmax
197!
198! 1128 2013-04-12 06:19:32Z raasch
199! +background_communication in inipar
200!
201! 1115 2013-03-26 18:16:16Z hoffmann
202! unused variables removed
203!
204! 1092 2013-02-02 11:24:22Z raasch
205! unused variables removed
206!
207! 1065 2012-11-22 17:42:36Z hoffmann
208! +nc, c_sedimentation, limiter_sedimentation, turbulence
209! -mu_constant, mu_constant_value
210!
211! 1053 2012-11-13 17:11:03Z hoffmann
212! necessary expansions according to the two new prognostic equations (nr, qr)
213! of the two-moment cloud physics scheme and steering parameters:
214! +*_init, *_surface, *_surface_initial_change, *_vertical_gradient,
215! +*_vertical_gradient_level, surface_waterflux_*,
216! +cloud_scheme, drizzle, mu_constant, mu_constant_value, ventilation_effect
217!
218! 1036 2012-10-22 13:43:42Z raasch
219! code put under GPL (PALM 3.9)
220!
221! 1015 2012-09-27 09:23:24Z raasch
222! -adjust_mixing_length
223!
224! 1003 2012-09-14 14:35:53Z raasch
225! -grid_matching
226!
227! 1001 2012-09-13 14:08:46Z raasch
228! -cut_spline_overshoot, long_filter_factor, overshoot_limit_*, ups_limit_*
229!
230! 996 2012-09-07 10:41:47Z raasch
231! -use_prior_plot1d_parameters
232!
233! 978 2012-08-09 08:28:32Z fricke
234! -km_damp_max, outflow_damping_width
235! +pt_damping_factor, pt_damping_width
236! +z0h_factor
237!
238! 964 2012-07-26 09:14:24Z raasch
239! -cross_normalized_x, cross_normalized_y, cross_xtext, z_max_do1d,
240! z_max_do1d_normalized
241!
242! 940 2012-07-09 14:31:00Z raasch
243! +neutral in inipar
244!
245! 927 2012-06-06 19:15:04Z raasch
246! +masking_method in inipar
247!
248! 824 2012-02-17 09:09:57Z raasch
249! +curvature_solution_effects in inipar
250!
251! 809 2012-01-30 13:32:58Z maronga
252! Bugfix: replaced .AND. and .NOT. with && and ! in the preprocessor directives
253!
254! 807 2012-01-25 11:53:51Z maronga
255! New cpp directive "__check" implemented which is used by check_namelist_files
256!
257! Revision 1.1  1997/07/24 11:22:50  raasch
258! Initial revision
259!
260!
261! Description:
262! ------------
263!> This subroutine reads variables controling the run from the NAMELIST files
264!------------------------------------------------------------------------------!
265 SUBROUTINE parin
266 
267
268    USE arrays_3d,                                                             &
269        ONLY:  pt_init, q_init, ref_state, s_init, sa_init, ug, u_init, v_init,&
270               vg
271
272    USE plant_canopy_model_mod,                                                &
273         ONLY: pcm_parin 
274
275    USE control_parameters
276
277    USE cpulog,                                                                &
278        ONLY:  cpu_log_barrierwait
279
280    USE dvrp_variables,                                                        &
281        ONLY:  local_dvrserver_running
282
283    USE flight_mod,                                                            &
284        ONLY:  flight_parin
285
286    USE grid_variables,                                                        &
287        ONLY:  dx, dy
288
289    USE indices,                                                               &
290        ONLY:  nx, ny, nz
291
292    USE kinds
293
294    USE land_surface_model_mod,                                                &
295        ONLY: lsm_parin
296
297    USE microphysics_mod,                                                      &
298        ONLY:  c_sedimentation, cloud_water_sedimentation,                     &
299               collision_turbulence, curvature_solution_effects_bulk,          &
300               limiter_sedimentation ,na_init, nc_const, ventilation_effect
301
302    USE model_1d,                                                              &
303        ONLY:  damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
304
305    USE pegrid
306
307    USE netcdf_interface,                                                      &
308        ONLY:  netcdf_data_format, netcdf_deflate, netcdf_precision
309
310    USE pmc_interface,                                                         &
311        ONLY:  nested_run, nesting_mode
312
313    USE profil_parameter,                                                      &
314        ONLY:  cross_profiles, cross_ts_uymax, cross_ts_uymin,                 &
315               profile_columns, profile_rows
316
317    USE progress_bar,                                                          &
318        ONLY :  batch_job
319
320    USE radiation_model_mod,                                                   &
321        ONLY: radiation_parin 
322
323    USE spectra_mod,                                                           &
324        ONLY :  spectra_parin
325
326    USE statistics,                                                            &
327        ONLY:  hom, hom_sum, pr_palm, region, statistic_regions
328
329    USE synthetic_turbulence_generator_mod,                                    &
330        ONLY:  stg_parin
331
332    USE urban_surface_mod,                                                     &
333        ONLY: usm_parin
334
335    USE wind_turbine_model_mod,                                                &
336        ONLY:  wtm_parin
337
338
339    IMPLICIT NONE
340
341    INTEGER(iwp) ::  i   !<
342
343
344    NAMELIST /inipar/  alpha_surface, approximation, bc_e_b, bc_lr,            &
345                       bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b,        &
346             bc_q_t,bc_s_b, bc_s_t, bc_sa_t, bc_uv_b, bc_uv_t,                 &
347             bottom_salinityflux, building_height, building_length_x,          &
348             building_length_y, building_wall_left, building_wall_south,       &
349             call_psolver_at_all_substeps, call_microphysics_at_all_substeps,  &
350             canyon_height,                                                    &
351             canyon_width_x, canyon_width_y, canyon_wall_left,                 &
352             canyon_wall_south, c_sedimentation, cfl_factor, cloud_droplets,   &
353             cloud_physics, cloud_scheme, cloud_top_radiation,                 &
354             cloud_water_sedimentation,                                        &
355             collective_wait, collision_turbulence, conserve_volume_flow,      &
356             conserve_volume_flow_mode, constant_flux_layer,                   &
357             coupling_start_time, curvature_solution_effects_bulk,             &
358             cycle_mg, damp_level_1d,                                          &
359             dissipation_1d,                                                   &
360             dp_external, dp_level_b, dp_smooth, dpdxy,                        &
361             dt, dt_pr_1d, dt_run_control_1d, dx, dy, dz, dz_max,              &
362             dz_stretch_factor, dz_stretch_level, end_time_1d,                 &
363             ensemble_member_nr, e_init, e_min, fft_method,                    &
364             flux_input_mode, flux_output_mode,                                &
365             galilei_transformation, humidity,                                 &
366             inflow_damping_height, inflow_damping_width,                      &
367             inflow_disturbance_begin, inflow_disturbance_end,                 &
368             initializing_actions, km_constant,                                &
369             large_scale_forcing, large_scale_subsidence,                      &
370             limiter_sedimentation,                                            &
371             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
372             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
373             most_method, na_init, nc_const, netcdf_precision, neutral, ngsrb, &
374             nsor, nsor_ini, nudging, nx, ny, nz, ocean, omega, omega_sor,     &
375             outflow_source_plane, passive_scalar, phi,                        &
376             prandtl_number, precipitation, psolver, pt_damping_factor,        &
377             pt_damping_width, pt_reference, pt_surface,                       &
378             pt_surface_initial_change, pt_vertical_gradient,                  &
379             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
380             q_vertical_gradient, q_vertical_gradient_level,                   &
381             random_generator, random_heatflux,                                &
382             rayleigh_damping_factor, rayleigh_damping_height,                 &
383             recycling_width, recycling_yshift,                                &
384             reference_state, residual_limit,                                  &
385             roughness_length, sa_surface,                                     &
386             sa_vertical_gradient, sa_vertical_gradient_level, scalar_advec,   &
387             scalar_rayleigh_damping,                                          &
388             statistic_regions, subs_vertical_gradient,                        &
389             subs_vertical_gradient_level, surface_heatflux, surface_pressure, &
390             surface_scalarflux, surface_waterflux,                            &
391             s_surface, s_surface_initial_change, s_vertical_gradient,         &
392             s_vertical_gradient_level, timestep_scheme,                       &
393             topography, topography_grid_convention, top_heatflux,             &
394             top_momentumflux_u, top_momentumflux_v, top_salinityflux,         &
395             top_scalarflux, transpose_compute_overlap,                        &
396             tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,     &
397             tunnel_wall_depth,                                                &
398             turbulent_inflow, turbulent_outflow,                              &
399             use_subsidence_tendencies, ug_surface, ug_vertical_gradient,      &
400             ug_vertical_gradient_level, use_surface_fluxes, use_cmax,         &
401             use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke,      &
402             uv_heights, u_bulk, u_profile, vg_surface, vg_vertical_gradient,  &
403             vg_vertical_gradient_level, v_bulk, v_profile, ventilation_effect,&
404             wall_adjustment, wall_heatflux, wall_humidityflux,                &
405             wall_salinityflux, wall_scalarflux, zeta_max, zeta_min, z0h_factor
406     
407    NAMELIST /d3par/  averaging_interval, averaging_interval_pr,               &
408             cpu_log_barrierwait, create_disturbances,                         &
409             cross_profiles, cross_ts_uymax, cross_ts_uymin,                   &
410             data_output, data_output_masks,                                   &
411             data_output_pr, data_output_2d_on_each_pe, disturbance_amplitude, &
412             disturbance_energy_limit, disturbance_level_b,                    &
413             disturbance_level_t, do2d_at_begin, do3d_at_begin,                &
414             dt, dt_averaging_input, dt_averaging_input_pr,                    &
415             dt_coupling, dt_data_output, dt_data_output_av, dt_disturb,       &
416             dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy,         &
417             dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart,              &
418             dt_run_control,end_time, force_print_header, mask_scale_x,        &
419             mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop,  &
420             mask_y_loop, mask_z_loop, netcdf_data_format, netcdf_deflate,     &
421             normalizing_region, npex, npey, nz_do3d,                          &
422             precipitation_amount_interval, profile_columns, profile_rows,     &
423             restart_time, section_xy, section_xz, section_yz,                 &
424             skip_time_data_output, skip_time_data_output_av, skip_time_dopr,  &
425             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
426             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
427             termination_time_needed, z_max_do2d
428
429
430    NAMELIST /envpar/  batch_job, host, local_dvrserver_running,               &
431                       maximum_cpu_time_allowed, maximum_parallel_io_streams,  &
432                       revision, return_addres, return_username,               &
433                       run_identifier, tasks_per_node, write_binary
434
435!
436!-- First read values of environment variables (this NAMELIST file is
437!-- generated by mrun)
438    CALL location_message( 'reading environment parameters from ENVPAR', .FALSE. )
439    OPEN ( 90, FILE='ENVPAR', STATUS='OLD', FORM='FORMATTED', ERR=30 )
440    READ ( 90, envpar, ERR=31, END=32 )
441    CLOSE ( 90 )
442    CALL location_message( 'finished', .TRUE. )
443
444!
445!-- Calculate the number of groups into which parallel I/O is split.
446!-- The default for files which are opened by all PEs (or where each
447!-- PE opens his own independent file) is, that all PEs are doing input/output
448!-- in parallel at the same time. This might cause performance or even more
449!-- severe problems depending on the configuration of the underlying file
450!-- system.
451!-- First, set the default:
452    IF ( maximum_parallel_io_streams == -1  .OR.                               &
453         maximum_parallel_io_streams > numprocs )  THEN
454       maximum_parallel_io_streams = numprocs
455    ENDIF
456!
457!-- Now calculate the number of io_blocks and the io_group to which the
458!-- respective PE belongs. I/O of the groups is done in serial, but in parallel
459!-- for all PEs belonging to the same group.
460!-- These settings are repeated in init_pegrid for the communicator comm2d,
461!-- which is not available here
462    io_blocks = numprocs / maximum_parallel_io_streams
463    io_group  = MOD( myid+1, io_blocks )
464
465    CALL location_message( 'reading NAMELIST parameters from PARIN', .FALSE. )
466!
467!-- Data is read in parallel by groups of PEs
468    DO  i = 0, io_blocks-1
469       IF ( i == io_group )  THEN
470
471!
472!--       Open the NAMELIST-file which is send with this job
473          CALL check_open( 11 )
474
475!
476!--       Read the control parameters for initialization.
477!--       The namelist "inipar" must be provided in the NAMELIST-file.
478          READ ( 11, inipar, ERR=10, END=11 )
479
480          GOTO 12
481
482 10       message_string = 'errors in \$inipar &or no \$inipar-namelist ' //   &
483                           'found (CRAY-machines only)'
484          CALL message( 'parin', 'PA0271', 1, 2, 0, 6, 0 )
485
486 11       message_string = 'no \$inipar-namelist found'
487          CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 )
488
489!
490!--       If required, read control parameters from restart file (produced by
491!--       a prior run). All PEs are reading from file created by PE0 (see
492!--       check_open)
493 12       IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
494             CALL read_var_list
495!
496!--          The restart file will be reopened when reading the subdomain data
497             CALL close_file( 13 )
498
499!
500!--          Increment the run count
501             runnr = runnr + 1
502          ENDIF
503
504!
505!--       In case of nested runs, explicitly set nesting boundary conditions.
506!--       This will overwrite the user settings. bc_lr and bc_ns always need
507!--       to be cyclic for vertical nesting.
508          IF ( nesting_mode == 'vertical' )  THEN
509             IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' )  THEN
510                WRITE ( message_string, *)  'bc_lr and bc_ns were set to ,',   &
511                                      'cyclic for vertical nesting'
512                CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 )
513                bc_lr   = 'cyclic'
514                bc_ns   = 'cyclic'
515             ENDIF
516             IF ( nest_domain )  THEN
517                bc_uv_t = 'nested'
518                bc_pt_t = 'nested'
519                bc_q_t  = 'nested'
520                bc_s_t  = 'nested'
521                bc_p_t  = 'neumann'
522             ENDIF
523          ELSE
524         
525!
526!--       For other nesting modes only set boundary conditions for
527!--       nested domains.
528             IF ( nest_domain )  THEN
529                bc_lr   = 'nested'
530                bc_ns   = 'nested'
531                bc_uv_t = 'nested'
532                bc_pt_t = 'nested'
533                bc_q_t  = 'nested'
534                bc_s_t  = 'nested'
535                bc_p_t  = 'neumann'
536             ENDIF
537          ENDIF
538
539!         
540!--       In case of nested runs, make sure that initializing_actions =
541!--       'set_constant_profiles' even though the constant-profiles
542!--       initializations for the prognostic variables will be overwritten
543!--       by pmci_child_initialize and pmci_parent_initialize. This is,
544!--       however, important e.g. to make sure that diagnostic variables
545!--       are set properly.
546          IF ( nest_domain )  THEN
547             initializing_actions = 'set_constant_profiles'
548          ENDIF
549
550           
551!
552!--       Check validity of lateral boundary conditions. This has to be done
553!--       here because they are already used in init_pegrid and init_grid and
554!--       therefore cannot be check in check_parameters
555          IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
556               bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested' )  THEN
557             message_string = 'unknown boundary condition: bc_lr = "' // &
558                              TRIM( bc_lr ) // '"'
559             CALL message( 'check_parameters', 'PA0049', 1, 2, 0, 6, 0 )
560          ENDIF
561          IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
562               bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested' )  THEN
563             message_string = 'unknown boundary condition: bc_ns = "' // &
564                              TRIM( bc_ns ) // '"'
565             CALL message( 'check_parameters', 'PA0050', 1, 2, 0, 6, 0 )
566          ENDIF
567
568!
569!--       Set internal variables used for speed optimization in if clauses
570          IF ( bc_lr /= 'cyclic' )               bc_lr_cyc    = .FALSE.
571          IF ( bc_lr == 'dirichlet/radiation' )  bc_lr_dirrad = .TRUE.
572          IF ( bc_lr == 'radiation/dirichlet' )  bc_lr_raddir = .TRUE.
573          IF ( bc_ns /= 'cyclic' )               bc_ns_cyc    = .FALSE.
574          IF ( bc_ns == 'dirichlet/radiation' )  bc_ns_dirrad = .TRUE.
575          IF ( bc_ns == 'radiation/dirichlet' )  bc_ns_raddir = .TRUE.
576
577!
578!--       Definition of names of areas used for computing statistics. They must
579!--       be defined at this place, because they are allowed to be redefined by
580!--       the user in user_parin.
581          region = 'total domain'
582
583!
584!--       Read runtime parameters given by the user for this run (namelist
585!--       "d3par"). The namelist "d3par" can be omitted. In that case, default
586!--       values are used for the parameters.
587          READ ( 11, d3par, END=20 )
588 20       CONTINUE
589
590!
591!--       Check if land surface model is used and read &lsm_par if required
592          CALL lsm_parin
593
594!
595!--       Check if urban surface model is used and read &urban_surface_par if required
596          CALL usm_parin
597
598!
599!--       Check if spectra shall be calculated and read spectra_par if required
600          CALL spectra_parin
601
602!
603!--       Check if radiation model is used and read &radiation_par if required
604          CALL radiation_parin
605 
606 
607!--       Check if plant canopy model is used and read &canopy_par if required
608          CALL pcm_parin
609 
610!
611!--       Read control parameters for optionally used model software packages
612          CALL package_parin
613
614!
615!--       Check if wind turbine model is used and read &wind_turbine_par if
616!--       required
617          CALL wtm_parin
618!
619!--       Check if virtual flights should be carried out and read &flight_par
620!--       if required
621          CALL flight_parin
622!
623!--       Check if synthetic turbulence generator is used and read stg_par if
624!--       required
625          CALL stg_parin
626!
627!--       Read user-defined variables
628          CALL user_parin
629
630!
631!--       Check in case of initial run, if the grid point numbers are well
632!--       defined and allocate some arrays which are already needed in
633!--       init_pegrid or check_parameters. During restart jobs, these arrays
634!--       will be allocated in read_var_list. All other arrays are allocated
635!--       in init_3d_model.
636          IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
637
638             IF ( nx <= 0 )  THEN
639                WRITE( message_string, * ) 'no value or wrong value given',    &
640                                           ' for nx: nx=', nx
641                CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 )
642             ENDIF
643             IF ( ny <= 0 )  THEN
644                WRITE( message_string, * ) 'no value or wrong value given',    &
645                                           ' for ny: ny=', ny
646                CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 )
647             ENDIF
648             IF ( nz <= 0 )  THEN
649                WRITE( message_string, * ) 'no value or wrong value given',    &
650                                           ' for nz: nz=', nz
651                CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 )
652             ENDIF
653!
654!--          ATTENTION: in case of changes to the following statement please
655!--                  also check the allocate statement in routine read_var_list
656             ALLOCATE( pt_init(0:nz+1), q_init(0:nz+1), s_init(0:nz+1),        &
657                       ref_state(0:nz+1), sa_init(0:nz+1), ug(0:nz+1),         &
658                       u_init(0:nz+1), v_init(0:nz+1), vg(0:nz+1),             &
659                       hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions),  &
660                       hom_sum(0:nz+1,pr_palm+max_pr_user,0:statistic_regions) )
661
662             hom = 0.0_wp
663
664          ENDIF
665
666!
667!--       NAMELIST-file is not needed anymore
668          CALL close_file( 11 )
669
670       ENDIF
671#if defined( __parallel )
672       CALL MPI_BARRIER( MPI_COMM_WORLD, ierr )
673#endif
674    ENDDO
675
676    CALL location_message( 'finished', .TRUE. )
677
678    RETURN
679
680 30 message_string = 'local file ENVPAR not found' //                          &
681                     '&some variables for steering may not be properly set'
682    CALL message( 'parin', 'PA0276', 0, 1, 0, 6, 0 )
683    RETURN
684
685 31 message_string = 'errors in local file ENVPAR' //                          &
686                     '&some variables for steering may not be properly set'
687    CALL message( 'parin', 'PA0277', 0, 1, 0, 6, 0 )
688    RETURN
689
690 32 message_string = 'no envpar-NAMELIST found in local file ENVPAR'  //       &
691                     '&some variables for steering may not be properly set'
692    CALL message( 'parin', 'PA0278', 0, 1, 0, 6, 0 )
693
694 END SUBROUTINE parin
Note: See TracBrowser for help on using the repository browser.