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

Last change on this file since 4173 was 4173, checked in by gronemeier, 3 years ago

add vdi_internal_controls

  • Property svn:keywords set to Id
File size: 47.8 KB
Line 
1!> @file parin.f90
2!------------------------------------------------------------------------------!
3! This file is part of the PALM model system.
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-2019 Leibniz Universitaet Hannover
18!------------------------------------------------------------------------------!
19!
20! Current revisions:
21! -----------------
22!
23!
24! Former revisions:
25! -----------------
26! $Id: parin.f90 4173 2019-08-20 12:04:06Z gronemeier $
27! add vdi_internal_controls
28!
29! 4131 2019-08-02 11:06:18Z monakurppa
30! Allocate hom and hom_sum to allow profile output for salsa variables.
31!
32! 4079 2019-07-09 18:04:41Z suehring
33! +monotonic_limiter_z
34!
35! 4022 2019-06-12 11:52:39Z suehring
36! Change default top boundary condition for pressure to Neumann in offline
37! nesting case
38!
39! 4017 2019-06-06 12:16:46Z schwenkel
40! Introduce alternative switch for debug output during timestepping
41!
42! 3885 2019-04-11 11:29:34Z kanani
43! Changes related to global restructuring of location messages and introduction
44! of additional debug messages
45!
46! 3806 2019-03-21 12:45:50Z raasch
47! additional check for lateral boundary conditions added
48!
49! 3747 2019-02-16 15:15:23Z gronemeier
50! removed setting of parameter region
51!
52! 3746 2019-02-16 12:41:27Z gronemeier
53! Removed most_method
54!
55! 3649 2019-01-02 16:52:21Z suehring
56! Delete debug-print statements
57!
58! 3646 2018-12-28 17:58:49Z kanani
59! Bugfix: for restarts, reset data_output_during_spinup to FALSE to enable
60! correct calculation of ntdim in check_parameters
61!
62! 3637 2018-12-20 01:51:36Z knoop
63! Implementation of the PALM module interface
64!
65! 3569 2018-11-27 17:03:40Z kanani
66! dom_dwd_user, Schrempf:
67! Remove uv exposure model code, this is now part of biometeorology_mod.
68!
69! 3545 2018-11-21 11:19:41Z gronemeier
70! remove rans_mode from initialization_parameters
71!
72! 3525 2018-11-14 16:06:14Z kanani
73! Changes related to clean-up of biometeorology (dom_dwd_user)
74!
75! 3473 2018-10-30 20:50:15Z suehring
76! Add virtual measurement module
77!
78! 3472 2018-10-30 20:43:50Z suehring
79! Add indoor model (kanani, srissman, tlang),
80! minor formatting
81!
82! 3467 2018-10-30 19:05:21Z suehring
83! Implementation of a new aerosol module salsa.
84!
85! 3448 2018-10-29 18:14:31Z kanani
86! Add biometeorology
87!
88! 3435 2018-10-26 18:25:44Z gronemeier
89! Add mask_k_over_surface
90!
91! 3421 2018-10-24 18:39:32Z gronemeier
92! Added module for data output at surfaces
93!
94! 3355 2018-10-16 14:03:34Z knoop
95! - offline nesting separated from large-scale forcing module
96! - top boundary condition for pressure in offline nesting changed
97!
98! 3343 2018-10-15 10:38:52Z suehring
99! Introduced reading of date_init to inipar.(Russo)
100! Add extra profiles for chemistry (basit)
101!
102! 3294 2018-10-01 02:37:10Z raasch
103! changes concerning modularization of ocean option
104!
105! 3274 2018-09-24 15:42:55Z knoop
106! Modularization of all bulk cloud physics code components
107!
108! 3248 2018-09-14 09:42:06Z sward
109! Minor formating changes
110!
111! 3246 2018-09-13 15:14:50Z sward
112! Added error handling for input namelist via parin_fail_message
113!
114! 3240 2018-09-12 12:04:40Z Giersch
115! A check that controls the number of user-defined profiles on the restart file
116! with the one given for the current run has been added.
117!
118! 3204 2018-08-23 10:49:05Z raasch
119! additional check for nz
120!
121! 3184 2018-07-27 17:13:09Z suehring
122! Bugfix, comment setting of chemistry and passive scalar top boundary condition
123! in case of offline nesting
124!
125! 3183 2018-07-27 14:25:55Z suehring
126! Rename variables and boundary conditions in mesoscale-offline nesting mode
127!
128! 3182 2018-07-27 13:36:03Z suehring
129! Added multi agent system
130!
131! 3157 2018-07-19 21:08:49Z maronga
132! added use_free_convection_scaling
133!
134! 3083 2018-06-19 14:03:12Z gronemeier
135! Added rans_const_c and rans_const_sigma as input parameters (TG)
136!
137! 3065 2018-06-12 07:03:02Z Giersch
138! New initialization parameters added
139!
140! 3049 2018-05-29 13:52:36Z Giersch
141! Error messages revised
142!
143! 3045 2018-05-28 07:55:41Z Giersch
144! z_max_do2d removed, error messages revised
145!
146! 2995 2018-04-19 12:13:16Z Giersch
147! time_since_reference_point must be calculated/initialized before the first 
148! call of functions related to the radiation model which occur in
149! time_integration_spinup or time_integration
150!
151! 2980 2018-04-17 15:19:27Z suehring
152! Revise message call
153!
154! 2975 2018-04-16 15:22:20Z suehring
155! - Informative message when initializing_actions has been changed
156!   to set_constant_profile in child domain
157! - Change location in message call
158!
159! 2967 2018-04-13 11:22:08Z raasch
160! bugfix: missing parallel cpp-directives added
161!
162! 2941 2018-04-03 11:54:58Z kanani
163! Fix for spinup in case of restart run
164!
165! 2938 2018-03-27 15:52:42Z suehring
166! Change initialization in case child domain should be initialized with Inifor.
167!
168! 2936 2018-03-27 14:49:27Z suehring
169! inipar renamed to initialization_parameters.
170! d3par renamed to runtime_parameters.
171!
172! 2921 2018-03-22 15:05:23Z Giersch
173! Activation of spinup has been moved from lsm/usm_parin to parin itself
174!
175! 2906 2018-03-19 08:56:40Z Giersch
176! ENVIRONMENT variables read/write_svf has been added
177!
178! 2894 2018-03-15 09:17:58Z Giersch
179! read_var_list has been renamed to rrd_global, all module related _parin
180! routines are called before reading the global restart data to overwrite them
181! in case of restart runs
182!
183! 2881 2018-03-13 16:24:40Z suehring
184! Added flag for switching on/off calculation of soil moisture
185!
186! 2849 2018-03-05 10:49:33Z Giersch
187! Position of d3par namelist in parameter file is unimportant now
188!
189! 2826 2018-02-21 12:39:28Z Giersch
190! Bugfix in setting the default boundary conditions for nest domains
191!
192! 2817 2018-02-19 16:32:21Z knoop
193! Preliminary gust module interface implemented
194!
195! 2773 2018-01-30 14:12:54Z suehring
196! Nesting for chemical species implemented
197!
198! 2766 2018-01-22 17:17:47Z kanani
199! Removed preprocessor directive __chem
200!
201! 2718 2018-01-02 08:49:38Z maronga
202! Corrected "Former revisions" section
203!
204! 2696 2017-12-14 17:12:51Z kanani
205! Change in file header (GPL part)
206! Implementation of uv exposure model (FK)
207! Added rans_mode and turbulence_closure to inipar (TG)
208! Implementation of chemistry module
209! Sorting of USE list (FK)
210! Forcing implemented, and initialization with inifor (MS)
211!
212! 2600 2017-11-01 14:11:20Z raasch
213! some comments added and variables renamed concerning r2599
214!
215! 2599 2017-11-01 13:18:45Z hellstea
216! The i/o grouping is updated to work correctly also in nested runs.
217!
218! 2575 2017-10-24 09:57:58Z maronga
219! Renamed phi -> latitude, added longitude
220!
221! 2563 2017-10-19 15:36:10Z Giersch
222! Changed position where restart files are closed.
223!
224! 2550 2017-10-16 17:12:01Z boeske
225! Added complex_terrain
226!
227! 2544 2017-10-13 18:09:32Z maronga
228! Moved day_of_year_init and time_utc_init to inipar.
229!
230! 2397 2017-09-04 16:22:48Z suehring
231! Enable initialization of 3d model by user in the child domain.
232!
233! 2375 2017-08-29 14:10:28Z schwenkel
234! Added aerosol initialization for bulk microphysics
235!
236! 2372 2017-08-25 12:37:32Z sward
237! y_shift added to namelist
238!
239! 2365 2017-08-21 14:59:59Z kanani
240! Vertical grid nesting: add vnest_start_time to d3par (SadiqHuq)
241!
242! 2339 2017-08-07 13:55:26Z gronemeier
243! corrected timestamp in header
244!
245! 2338 2017-08-07 12:15:38Z gronemeier
246! Modularize 1D model
247!
248! 2310 2017-07-11 09:37:02Z gronemeier
249! Bugfix: re-arranged call for error messages for ENVPAR file
250!
251! 2304 2017-07-04 14:35:55Z suehring
252! Bugfix, enable restarts for child domain.
253!
254! 2298 2017-06-29 09:28:18Z raasch
255! -return_addres, return_username in ENVPAR, -cross_ts_uymax, cross_ts_uymin in
256! d3par
257!
258! 2296 2017-06-28 07:53:56Z maronga
259! Added parameters for model spinup
260!
261! 2292 2017-06-20 09:51:42Z schwenkel
262! Implementation of new microphysic scheme: cloud_scheme = 'morrison'
263! includes two more prognostic equations for cloud drop concentration (nc) 
264! and cloud water content (qc).
265!
266! 2267 2017-06-09 09:33:25Z gronemeier
267! Bugfix: removed skipping of reading namelists in case of omitted d3par
268!
269! 2259 2017-06-08 09:09:11Z gronemeier
270! Implemented synthetic turbulence generator
271!
272! 2233 2017-05-30 18:08:54Z suehring
273!
274! 2232 2017-05-30 17:47:52Z suehring
275! typo corrected
276! +wall_salinityflux
277! +tunnel_height, tunnel_lenght, tunnel_width_x, tunnel_width_y,
278!  tunnel_wall_depth
279!
280! 2118 2017-01-17 16:38:49Z raasch
281! -background_communication from inipar
282!
283! 2050 2016-11-08 15:00:55Z gronemeier
284! Implement turbulent outflow condition
285!
286! 2037 2016-10-26 11:15:40Z knoop
287! Anelastic approximation implemented
288!
289! 2035 2016-10-24 15:06:17Z suehring
290! Remove check for npex and npey in nesting case
291!
292! 2011 2016-09-19 17:29:57Z kanani
293! Added flag lsf_exception to allow explicit enabling of large scale forcing
294! together with buildings on flat terrain.
295!
296! 2007 2016-08-24 15:47:17Z kanani
297! Added call to urban surface model for reading of &urban_surface_par
298!
299! 2004 2016-08-24 10:25:59Z suehring
300! Humidity and passive scalar treated separately in nesting mode
301!
302! 2000 2016-08-20 18:09:15Z knoop
303! Forced header and separation lines into 80 columns
304!
305! 1992 2016-08-12 15:14:59Z suehring
306! +top_scalarflux
307!
308! 1960 2016-07-12 16:34:24Z suehring
309! Allocation of s_init
310!
311! 1957 2016-07-07 10:43:48Z suehring
312! flight module added
313!
314! 1955 2016-07-01 12:38:59Z hellstea
315! The parameter intializating_actions is set to 'set_constant_profiles for
316! all nest domains in order to make sure that diagnostic variables are properly
317! initialized for nest domains. Prognostic variables are later initialized by
318! interpolation from the parent domain.
319!
320! 1917 2016-05-27 14:28:12Z witha
321! Initial version of purely vertical nesting introduced.
322!
323! 1914 2016-05-26 14:44:07Z witha
324! Added call to wind turbine model for reading of &wind_turbine_par
325!
326! 1849 2016-04-08 11:33:18Z hoffmann
327! Adapted for modularization of microphysics
328!
329! 1833 2016-04-07 14:23:03Z raasch
330! call of spectra_parin
331!
332! 1831 2016-04-07 13:15:51Z hoffmann
333! turbulence renamed collision_turbulence, drizzle renamed
334! cloud_water_sedimentation
335! curvature_solution_effects removed
336!
337! 1826 2016-04-07 12:01:39Z maronga
338! Added call to radiation model for reading of &radiation_par.
339! Added call to plant canopy model for reading of &canopy_par.
340!
341! 1817 2016-04-06 15:44:20Z maronga
342! Added call to land surface model for reading of &lsm_par
343!
344! 1804 2016-04-05 16:30:18Z maronga
345! Removed code for parameter file check (__check)
346!
347! 1783 2016-03-06 18:36:17Z raasch
348! +netcdf_deflate in d3par, netcdf module and variable names changed
349!
350! 1764 2016-02-28 12:45:19Z raasch
351! cpp-statements for nesting removed, explicit settings of boundary conditions
352! in nest domains,
353! bugfix: npex/npey message moved from inipar to d3par
354! bugfix: check of lateral boundary conditions from check_parameters to here,
355! because they will be already used in init_pegrid and init_grid
356!
357! 1762 2016-02-25 12:31:13Z hellstea
358! Introduction of nested domain feature
359!
360! 1691 2015-10-26 16:17:44Z maronga
361! Added parameter most_method. Renamed prandtl_layer to constant_flux_layer.
362!
363! 1682 2015-10-07 23:56:08Z knoop
364! Code annotations made doxygen readable
365!
366! 1560 2015-03-06 10:48:54Z keck
367! +recycling_yshift
368!
369! 1496 2014-12-02 17:25:50Z maronga
370! Renamed: "radiation -> "cloud_top_radiation"
371!
372! 1484 2014-10-21 10:53:05Z kanani
373! Changes due to new module structure of the plant canopy model:
374!   canopy-model related parameters moved to new package canopy_par in
375!   subroutine package_parin
376!
377! 1429 2014-07-15 12:53:45Z knoop
378! +ensemble_member_nr to prepare the random_generator for ensemble runs
379!
380! 1402 2014-05-09 14:25:13Z raasch
381! location messages modified, progress_bar_disabled included in envpar-NAMELIST
382!
383! 1384 2014-05-02 14:31:06Z raasch
384! location messages added
385!
386! 1365 2014-04-22 15:03:56Z boeske
387! Usage of large scale forcing enabled:
388! +use_subsidence_tendencies
389!
390! 1361 2014-04-16 15:17:48Z hoffmann
391! +call_microphysics_at_all_substeps
392!
393! 1359 2014-04-11 17:15:14Z hoffmann
394! REAL constants provided with KIND-attribute
395!
396! 1353 2014-04-08 15:21:23Z heinze
397! REAL constants provided with KIND-attribute
398!
399! 1327 2014-03-21 11:00:16Z raasch
400! -data_output_format, do3d_compress, do3d_comp_prec
401!
402! 1320 2014-03-20 08:40:49Z raasch
403! ONLY-attribute added to USE-statements,
404! kind-parameters added to all INTEGER and REAL declaration statements,
405! kinds are defined in new module kinds,
406! old module precision_kind is removed,
407! revision history before 2012 removed,
408! comment fields (!:) to be used for variable explanations added to
409! all variable declaration statements
410!
411! 1318 2014-03-17 13:35:16Z raasch
412! +cpu_log_barrierwait in d3par
413!
414! 1301 2014-03-06 13:29:46Z heinze
415! +large_scale_subsidence
416!
417! 1241 2013-10-30 11:36:58Z heinze
418! +nudging
419! +large_scale_forcing
420!
421! 1216 2013-08-26 09:31:42Z raasch
422! +transpose_compute_overlap in inipar
423!
424! 1195 2013-07-01 12:27:57Z heinze
425! Bugfix: allocate ref_state
426!
427! 1179 2013-06-14 05:57:58Z raasch
428! +reference_state in inipar
429!
430! 1159 2013-05-21 11:58:22Z fricke
431! +use_cmax
432!
433! 1128 2013-04-12 06:19:32Z raasch
434! +background_communication in inipar
435!
436! 1115 2013-03-26 18:16:16Z hoffmann
437! unused variables removed
438!
439! 1092 2013-02-02 11:24:22Z raasch
440! unused variables removed
441!
442! 1065 2012-11-22 17:42:36Z hoffmann
443! +nc, c_sedimentation, limiter_sedimentation, turbulence
444! -mu_constant, mu_constant_value
445!
446! 1053 2012-11-13 17:11:03Z hoffmann
447! necessary expansions according to the two new prognostic equations (nr, qr)
448! of the two-moment cloud physics scheme and steering parameters:
449! +*_init, *_surface, *_surface_initial_change, *_vertical_gradient,
450! +*_vertical_gradient_level, surface_waterflux_*,
451! +cloud_scheme, drizzle, mu_constant, mu_constant_value, ventilation_effect
452!
453! 1036 2012-10-22 13:43:42Z raasch
454! code put under GPL (PALM 3.9)
455!
456! 1015 2012-09-27 09:23:24Z raasch
457! -adjust_mixing_length
458!
459! 1003 2012-09-14 14:35:53Z raasch
460! -grid_matching
461!
462! 1001 2012-09-13 14:08:46Z raasch
463! -cut_spline_overshoot, long_filter_factor, overshoot_limit_*, ups_limit_*
464!
465! 996 2012-09-07 10:41:47Z raasch
466! -use_prior_plot1d_parameters
467!
468! 978 2012-08-09 08:28:32Z fricke
469! -km_damp_max, outflow_damping_width
470! +pt_damping_factor, pt_damping_width
471! +z0h_factor
472!
473! 964 2012-07-26 09:14:24Z raasch
474! -cross_normalized_x, cross_normalized_y, cross_xtext, z_max_do1d,
475! z_max_do1d_normalized
476!
477! 940 2012-07-09 14:31:00Z raasch
478! +neutral in inipar
479!
480! 927 2012-06-06 19:15:04Z raasch
481! +masking_method in inipar
482!
483! 824 2012-02-17 09:09:57Z raasch
484! +curvature_solution_effects in inipar
485!
486! 809 2012-01-30 13:32:58Z maronga
487! Bugfix: replaced .AND. and .NOT. with && and ! in the preprocessor directives
488!
489! 807 2012-01-25 11:53:51Z maronga
490! New cpp directive "__check" implemented which is used by check_namelist_files
491!
492! Revision 1.1  1997/07/24 11:22:50  raasch
493! Initial revision
494!
495!
496! Description:
497! ------------
498!> This subroutine reads variables controling the run from the NAMELIST files
499!>
500!> @todo: Revise max_pr_cs (profiles for chemistry)
501!------------------------------------------------------------------------------!
502 SUBROUTINE parin
503 
504
505    USE arrays_3d,                                                             &
506        ONLY:  pt_init, q_init, ref_state, s_init, sa_init,                    &
507               ug, u_init, v_init, vg
508
509    USE chem_modules
510
511    USE control_parameters
512
513    USE cpulog,                                                                &
514        ONLY:  cpu_log_barrierwait
515
516    USE date_and_time_mod,                                                     &
517        ONLY:  date_init, day_of_year_init, time_utc_init
518
519    USE grid_variables,                                                        &
520        ONLY:  dx, dy
521
522    USE indices,                                                               &
523        ONLY:  nx, ny, nz
524
525    USE kinds
526
527    USE model_1d_mod,                                                          &
528        ONLY:  damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
529
530    USE module_interface,                                                      &
531        ONLY:  module_interface_parin
532
533    USE netcdf_interface,                                                      &
534        ONLY:  netcdf_data_format, netcdf_deflate, netcdf_precision
535
536    USE pegrid
537
538    USE pmc_interface,                                                         &
539        ONLY:  nested_run, nesting_mode
540
541    USE profil_parameter,                                                      &
542        ONLY:  cross_profiles, profile_columns, profile_rows
543
544    USE progress_bar,                                                          &
545        ONLY :  progress_bar_disabled
546
547    USE read_restart_data_mod,                                                 &
548        ONLY:  rrd_global
549
550    USE statistics,                                                            &
551        ONLY:  hom, hom_sum, pr_palm, statistic_regions
552
553    USE turbulence_closure_mod,                                                &
554        ONLY:  rans_const_c, rans_const_sigma
555
556    USE vertical_nesting_mod,                                                  &
557        ONLY:  vnest_start_time
558
559
560    IMPLICIT NONE
561
562    CHARACTER (LEN=80) ::  line  !< dummy string that contains the current line of the parameter file
563
564    INTEGER(iwp) ::  global_id      !< process id with respect to MPI_COMM_WORLD
565    INTEGER(iwp) ::  global_procs   !< # of procs with respect to MPI_COMM_WORLD
566    INTEGER(iwp) ::  i              !<
567    INTEGER(iwp) ::  ioerr          !< error flag for open/read/write
568
569    NAMELIST /inipar/  alpha_surface, approximation, bc_e_b,     &
570                       bc_lr, bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b, &
571             bc_q_t,bc_s_b, bc_s_t, bc_uv_b, bc_uv_t,                 &
572             building_height, building_length_x,          &
573             building_length_y, building_wall_left, building_wall_south,       &
574             calc_soil_moisture_during_spinup,                                 &
575             call_psolver_at_all_substeps,  &
576             canyon_height,                                                    &
577             canyon_width_x, canyon_width_y, canyon_wall_left,                 &
578             canyon_wall_south, cfl_factor, cloud_droplets,   &
579             collective_wait, complex_terrain,           &
580             conserve_volume_flow,                                             &
581             conserve_volume_flow_mode, constant_flux_layer,                   &
582             coupling_start_time,             &
583             cycle_mg, damp_level_1d,                                          &
584             data_output_during_spinup,                                        &
585             date_init,                                                        &
586             day_of_year_init,                                                 &
587             dissipation_1d,                                                   &
588             dp_external, dp_level_b, dp_smooth, dpdxy,    &
589             dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max,   &
590             dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
591             dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
592             e_min, fft_method, flux_input_mode, flux_output_mode,             &
593             galilei_transformation, humidity,                                 &
594             inflow_damping_height, inflow_damping_width,                      &
595             inflow_disturbance_begin, inflow_disturbance_end,                 &
596             initializing_actions, km_constant,                                &
597             large_scale_forcing, large_scale_subsidence, latitude,            &
598             longitude,                                 &
599             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
600             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
601             monotonic_limiter_z,                                              &
602             netcdf_precision, neutral, ngsrb,                                 &
603             nsor, nsor_ini, nudging, nx, ny, nz, ocean_mode, omega,           &
604             omega_sor, outflow_source_plane, passive_scalar,                  &
605             prandtl_number, psolver, pt_damping_factor,        &
606             pt_damping_width, pt_reference, pt_surface,                       &
607             pt_surface_initial_change, pt_vertical_gradient,                  &
608             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
609             q_vertical_gradient, q_vertical_gradient_level,                   &
610             random_generator, random_heatflux, rans_const_c, rans_const_sigma,&
611             rayleigh_damping_factor, rayleigh_damping_height,                 &
612             recycling_width, recycling_yshift,                                &
613             reference_state, residual_limit,                                  &
614             roughness_length,                                                 &
615             scalar_advec,   &
616             scalar_rayleigh_damping,                              &
617             spinup_time, spinup_pt_amplitude, spinup_pt_mean,                 &
618             statistic_regions, subs_vertical_gradient,                        &
619             subs_vertical_gradient_level, surface_heatflux, surface_pressure, &
620             surface_scalarflux, surface_waterflux,                            &
621             s_surface, s_surface_initial_change, s_vertical_gradient,         &
622             s_vertical_gradient_level, time_utc_init, timestep_scheme,        &
623             topography, topography_grid_convention, top_heatflux,             &
624             top_momentumflux_u, top_momentumflux_v,                           &
625             top_scalarflux, transpose_compute_overlap,                        &
626             tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,     &
627             tunnel_wall_depth, turbulence_closure,                            &
628             turbulent_inflow, turbulent_outflow,                              &
629             use_subsidence_tendencies, ug_surface, ug_vertical_gradient,      &
630             use_free_convection_scaling,                                      &
631             ug_vertical_gradient_level, use_surface_fluxes, use_cmax,         &
632             use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke,      &
633             uv_heights, u_bulk, u_profile, vdi_checks, vg_surface,            &
634             vg_vertical_gradient,  &
635             vg_vertical_gradient_level, v_bulk, v_profile,&
636             wall_adjustment, wall_heatflux, wall_humidityflux,                &
637             wall_scalarflux, y_shift, zeta_max, zeta_min,  &
638             z0h_factor
639
640    NAMELIST /initialization_parameters/  alpha_surface,                       &
641             approximation, bc_e_b,                                            &
642             bc_lr, bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b,           &
643             bc_q_t,bc_s_b, bc_s_t, bc_uv_b, bc_uv_t,                          &
644             building_height, building_length_x,                               &
645             building_length_y, building_wall_left, building_wall_south,       &
646             calc_soil_moisture_during_spinup,                                 &
647             call_psolver_at_all_substeps,                                     &
648             canyon_height,                                                    &
649             canyon_width_x, canyon_width_y, canyon_wall_left,                 &
650             canyon_wall_south, cfl_factor, cloud_droplets,                    &
651             collective_wait, complex_terrain,                                 &
652             conserve_volume_flow,                                             &
653             conserve_volume_flow_mode, constant_flux_layer,                   &
654             coupling_start_time,                                              &
655             cycle_mg, damp_level_1d,                                          &
656             data_output_during_spinup,                                        &
657             date_init,                                                        &
658             day_of_year_init,                                                 &
659             dissipation_1d,                                                   &
660             dp_external, dp_level_b, dp_smooth, dpdxy,                        &
661             dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max,   &
662             dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
663             dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
664             e_min, fft_method, flux_input_mode, flux_output_mode,             &
665             galilei_transformation, humidity,                                 &
666             inflow_damping_height, inflow_damping_width,                      &
667             inflow_disturbance_begin, inflow_disturbance_end,                 &
668             initializing_actions, km_constant,                                &
669             large_scale_forcing, large_scale_subsidence, latitude,            &
670             longitude,                                                        &
671             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
672             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
673             monotonic_limiter_z,                                              &
674             netcdf_precision, neutral, ngsrb,                                 &
675             nsor, nsor_ini, nudging, nx, ny, nz, ocean_mode, omega,           &
676             omega_sor, outflow_source_plane, passive_scalar,                  &
677             prandtl_number, psolver, pt_damping_factor,                       &
678             pt_damping_width, pt_reference, pt_surface,                       &
679             pt_surface_initial_change, pt_vertical_gradient,                  &
680             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
681             q_vertical_gradient, q_vertical_gradient_level,                   &
682             random_generator, random_heatflux, rans_const_c, rans_const_sigma,&
683             rayleigh_damping_factor, rayleigh_damping_height,                 &
684             recycling_width, recycling_yshift,                                &
685             reference_state, residual_limit,                                  &
686             roughness_length, scalar_advec,                                   &
687             scalar_rayleigh_damping,                                          &
688             spinup_time, spinup_pt_amplitude, spinup_pt_mean,                 &
689             statistic_regions, subs_vertical_gradient,                        &
690             subs_vertical_gradient_level, surface_heatflux, surface_pressure, &
691             surface_scalarflux, surface_waterflux,                            &
692             s_surface, s_surface_initial_change, s_vertical_gradient,         &
693             s_vertical_gradient_level, time_utc_init, timestep_scheme,        &
694             topography, topography_grid_convention, top_heatflux,             &
695             top_momentumflux_u, top_momentumflux_v,                           &
696             top_scalarflux, transpose_compute_overlap,                        &
697             tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,     &
698             tunnel_wall_depth, turbulence_closure,                            &
699             turbulent_inflow, turbulent_outflow,                              &
700             use_subsidence_tendencies, ug_surface, ug_vertical_gradient,      &
701             ug_vertical_gradient_level, use_surface_fluxes, use_cmax,         &
702             use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke,      &
703             use_free_convection_scaling,                                      &
704             uv_heights, u_bulk, u_profile, vdi_checks,                        &
705             vg_surface, vg_vertical_gradient,  &
706             vg_vertical_gradient_level, v_bulk, v_profile,                    &
707             wall_adjustment, wall_heatflux, wall_humidityflux,                &
708             wall_scalarflux, y_shift, zeta_max, zeta_min, z0h_factor
709             
710    NAMELIST /d3par/  averaging_interval, averaging_interval_pr,               &
711             cpu_log_barrierwait, create_disturbances,                         &
712             cross_profiles, data_output, data_output_masks,                   &
713             data_output_pr, data_output_2d_on_each_pe,                        &
714             debug_output,                                                     &
715             debug_output_timestep,                                            &
716             disturbance_amplitude,                                            &
717             disturbance_energy_limit, disturbance_level_b,                    &
718             disturbance_level_t, do2d_at_begin, do3d_at_begin,                &
719             dt, dt_averaging_input, dt_averaging_input_pr,                    &
720             dt_coupling, dt_data_output, dt_data_output_av, dt_disturb,       &
721             dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy,         &
722             dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart,              &
723             dt_run_control,end_time, force_print_header, mask_k_over_surface, &
724             mask_scale_x,                                                     &
725             mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop,  &
726             mask_y_loop, mask_z_loop, netcdf_data_format, netcdf_deflate,     &
727             normalizing_region, npex, npey, nz_do3d,                          &
728             profile_columns, profile_rows,     &
729             restart_time, section_xy, section_xz, section_yz,                 &
730             skip_time_data_output, skip_time_data_output_av, skip_time_dopr,  &
731             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
732             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
733             termination_time_needed, vnest_start_time
734
735    NAMELIST /runtime_parameters/  averaging_interval, averaging_interval_pr,  &
736             cpu_log_barrierwait, create_disturbances,                         &
737             cross_profiles, data_output, data_output_masks,                   &
738             data_output_pr, data_output_2d_on_each_pe,                        &
739             debug_output,                                                     &
740             debug_output_timestep,                                            &
741             disturbance_amplitude,                                            &
742             disturbance_energy_limit, disturbance_level_b,                    &
743             disturbance_level_t, do2d_at_begin, do3d_at_begin,                &
744             dt, dt_averaging_input, dt_averaging_input_pr,                    &
745             dt_coupling, dt_data_output, dt_data_output_av, dt_disturb,       &
746             dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy,         &
747             dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart,              &
748             dt_run_control,end_time, force_print_header, mask_k_over_surface, &
749             mask_scale_x,                                                     &
750             mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop,  &
751             mask_y_loop, mask_z_loop, netcdf_data_format, netcdf_deflate,     &
752             normalizing_region, npex, npey, nz_do3d,                          &
753             profile_columns, profile_rows,     &
754             restart_time, section_xy, section_xz, section_yz,                 &
755             skip_time_data_output, skip_time_data_output_av, skip_time_dopr,  &
756             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
757             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
758             termination_time_needed, vnest_start_time
759
760    NAMELIST /envpar/  progress_bar_disabled, host,                            &
761                       maximum_cpu_time_allowed, maximum_parallel_io_streams,  &
762                       read_svf, revision, run_identifier, tasks_per_node,     &
763                       write_binary, write_svf
764
765!
766!-- First read values of environment variables (this NAMELIST file is
767!-- generated by palmrun)
768    CALL location_message( 'reading environment parameters from ENVPAR', 'start' )
769
770    OPEN ( 90, FILE='ENVPAR', STATUS='OLD', FORM='FORMATTED', IOSTAT=ioerr )
771
772    IF ( ioerr /= 0 )  THEN
773       message_string = 'local file ENVPAR not found' //                       &
774                        '&some variables for steering may not be properly set'
775       CALL message( 'parin', 'PA0276', 0, 1, 0, 6, 0 )
776    ELSE
777       READ ( 90, envpar, IOSTAT=ioerr )
778       IF ( ioerr < 0 )  THEN
779          message_string = 'no envpar-NAMELIST found in local file '  //       &
780                           'ENVPAR& or some variables for steering may '  //   &
781                           'not be properly set'
782          CALL message( 'parin', 'PA0278', 0, 1, 0, 6, 0 )
783       ELSEIF ( ioerr > 0 )  THEN
784          message_string = 'errors in local file ENVPAR' //                    &
785                           '&some variables for steering may not be properly set'
786          CALL message( 'parin', 'PA0277', 0, 1, 0, 6, 0 )
787       ENDIF
788       CLOSE ( 90 )
789    ENDIF
790
791    CALL location_message( 'reading environment parameters from ENVPAR', 'finished' )
792!
793!-- Calculate the number of groups into which parallel I/O is split.
794!-- The default for files which are opened by all PEs (or where each
795!-- PE opens his own independent file) is, that all PEs are doing input/output
796!-- in parallel at the same time. This might cause performance or even more
797!-- severe problems depending on the configuration of the underlying file
798!-- system.
799!-- Calculation of the number of blocks and the I/O group must be based on all
800!-- PEs involved in this run. Since myid and numprocs are related to the
801!-- comm2d communicator, which gives only a subset of all PEs in case of
802!-- nested runs, that information must be inquired again from the global
803!-- communicator.
804!-- First, set the default:
805#if defined( __parallel )
806    CALL MPI_COMM_RANK( MPI_COMM_WORLD, global_id, ierr )
807    CALL MPI_COMM_SIZE( MPI_COMM_WORLD, global_procs, ierr )
808#else
809    global_id    = 0
810    global_procs = 1
811#endif
812    IF ( maximum_parallel_io_streams == -1  .OR.                               &
813         maximum_parallel_io_streams > global_procs )  THEN
814       maximum_parallel_io_streams = global_procs
815    ENDIF
816!
817!-- Now calculate the number of io_blocks and the io_group to which the
818!-- respective PE belongs. I/O of the groups is done in serial, but in parallel
819!-- for all PEs belonging to the same group.
820    io_blocks = global_procs / maximum_parallel_io_streams
821    io_group  = MOD( global_id+1, io_blocks )
822   
823    CALL location_message( 'reading NAMELIST parameters from PARIN', 'start' )
824!
825!-- Data is read in parallel by groups of PEs
826    DO  i = 0, io_blocks-1
827       IF ( i == io_group )  THEN
828
829!
830!--       Open the NAMELIST-file which is send with this job
831          CALL check_open( 11 )
832
833!
834!--       Read the control parameters for initialization.
835!--       The namelist "inipar" must be provided in the NAMELIST-file.
836          READ ( 11, initialization_parameters, ERR=10, END=11 )
837          GOTO 14
838         
839 10       BACKSPACE( 11 )
840          READ( 11 , '(A)') line
841          CALL parin_fail_message( 'initialization_parameters', line )
842
843 11       REWIND ( 11 )
844          READ ( 11, inipar, ERR=12, END=13 )
845 
846          message_string = 'namelist inipar is deprecated and will be ' //    &
847                          'removed in near future. & Please use namelist ' // &
848                          'initialization_parameters instead'
849          CALL message( 'parin', 'PA0017', 0, 1, 0, 6, 0 )
850 
851          GOTO 14
852 
853 12       BACKSPACE( 11 )
854          READ( 11 , '(A)') line
855          CALL parin_fail_message( 'inipar', line )
856
857 13       message_string = 'no initialization_parameters-namelist found'
858          CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 )
859
860!
861!--       Try to read runtime parameters given by the user for this run
862!--       (namelist "runtime_parameters"). The namelist "runtime_parmeters"   
863!--       can be omitted. In that case default values are used for the         
864!--       parameters.
865 14       line = ' '
866
867          REWIND ( 11 )
868          line = ' '
869          DO WHILE ( INDEX( line, '&runtime_parameters' ) == 0 )
870             READ ( 11, '(A)', END=16 )  line
871          ENDDO
872          BACKSPACE ( 11 )
873
874!
875!--       Read namelist
876          READ ( 11, runtime_parameters, ERR = 15 )
877          GOTO 18
878
879 15       BACKSPACE( 11 )
880          READ( 11 , '(A)') line
881          CALL parin_fail_message( 'runtime_parameters', line )
882
883 16       REWIND ( 11 )
884          line = ' '
885          DO WHILE ( INDEX( line, '&d3par' ) == 0 )
886             READ ( 11, '(A)', END=18 )  line
887          ENDDO
888          BACKSPACE ( 11 )
889
890!
891!--       Read namelist
892          READ ( 11, d3par, ERR = 17, END = 18 )
893
894          message_string = 'namelist d3par is deprecated and will be ' //      &
895                          'removed in near future. &Please use namelist ' //   &
896                          'runtime_parameters instead'
897          CALL message( 'parin', 'PA0487', 0, 1, 0, 6, 0 )
898
899          GOTO 18
900
901 17       BACKSPACE( 11 )
902          READ( 11 , '(A)') line
903          CALL parin_fail_message( 'd3par', line )
904
905 18       CONTINUE
906
907!
908!--       Check for module namelists and read them
909          CALL module_interface_parin
910
911!
912!--       If required, read control parameters from restart file (produced by
913!--       a prior run). All PEs are reading from file created by PE0 (see
914!--       check_open)
915          IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
916
917             CALL rrd_global
918!
919!--          Increment the run count
920             runnr = runnr + 1
921!
922!--          In case of a restart run, the number of user-defined profiles on
923!--          the restart file (already stored in max_pr_user) has to match the
924!--          one given for the current run. max_pr_user_tmp is calculated in
925!--          user_parin and max_pr_user is read in via rrd_global.
926             IF ( max_pr_user /= max_pr_user_tmp )  THEN
927                WRITE( message_string, * ) 'the number of user-defined ',      &
928                      'profiles given in data_output_pr (', max_pr_user_tmp,   &
929                      ') does not match the one ',                             &
930                      'found in the restart file (', max_pr_user, ')'
931                CALL message( 'user_parin', 'UI0009', 1, 2, 0, 6, 0 )
932             ENDIF
933          ELSE
934             max_pr_user = max_pr_user_tmp
935          ENDIF
936!
937!--       Activate spinup
938          IF ( land_surface .OR. urban_surface )  THEN
939             IF ( spinup_time > 0.0_wp )  THEN
940                coupling_start_time = spinup_time
941                time_since_reference_point = simulated_time - coupling_start_time
942                IF ( spinup_pt_mean == 9999999.9_wp )  THEN
943                   spinup_pt_mean = pt_surface
944                ENDIF
945                end_time = end_time + spinup_time
946                IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
947                   spinup = .TRUE.
948                ELSEIF (        TRIM( initializing_actions ) == 'read_restart_data'      &
949                         .AND.  time_since_reference_point > 0.0_wp )  THEN
950                   data_output_during_spinup = .FALSE.  !< required for correct ntdim calculation
951                                                        !< in check_parameters for restart run
952                ENDIF
953             ENDIF
954          ENDIF
955
956!
957!--       In case of nested runs, explicitly set nesting boundary conditions.
958!--       This will overwrite the user settings and basic defaults.
959!--       bc_lr and bc_ns always need to be cyclic for vertical nesting.
960          IF ( nested_run )  THEN
961             IF ( nesting_mode == 'vertical' )  THEN
962                IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' )  THEN
963                   WRITE ( message_string, *) 'bc_lr and bc_ns were set to ,', &
964                        'cyclic for vertical nesting'
965                   CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 )
966                   bc_lr   = 'cyclic'
967                   bc_ns   = 'cyclic'
968                ENDIF
969                IF ( child_domain )  THEN
970                   bc_uv_t  = 'nested'
971                   bc_pt_t  = 'nested'
972                   bc_q_t   = 'nested'
973                   bc_s_t   = 'nested'
974                   bc_cs_t  = 'nested'
975                   bc_p_t   = 'neumann' 
976                ENDIF
977!
978!--          For other nesting modes only set boundary conditions for
979!--          nested domains.
980             ELSE
981                IF ( child_domain )  THEN
982                   bc_lr    = 'nested'
983                   bc_ns    = 'nested'
984                   bc_uv_t  = 'nested'
985                   bc_pt_t  = 'nested'
986                   bc_q_t   = 'nested'
987                   bc_s_t   = 'nested'
988                   bc_cs_t  = 'nested'
989                   bc_p_t   = 'neumann'
990                ENDIF
991             ENDIF
992          ENDIF
993!
994!--       Set boundary conditions also in case the model is offline-nested in
995!--       larger-scale models.
996          IF ( nesting_offline )  THEN
997             bc_lr    = 'nesting_offline'
998             bc_ns    = 'nesting_offline'
999             bc_uv_t  = 'nesting_offline'
1000             bc_pt_t  = 'nesting_offline'
1001             bc_q_t   = 'nesting_offline'
1002           !  bc_s_t   = 'nesting_offline'  ! scalar boundary condition is not clear yet
1003           !  bc_cs_t  = 'nesting_offline'  ! same for chemical species
1004             bc_p_t   = 'neumann'
1005          ENDIF
1006
1007!         
1008!--       In case of nested runs, make sure that initializing_actions =
1009!--       'set_constant_profiles' even though the constant-profiles
1010!--       initializations for the prognostic variables will be overwritten
1011!--       by pmci_child_initialize and pmci_parent_initialize. This is,
1012!--       however, important e.g. to make sure that diagnostic variables
1013!--       are set properly. An exception is made in case of restart runs and
1014!--       if user decides to do everything by its own.
1015          IF ( child_domain  .AND.  .NOT. (                                    &
1016               TRIM( initializing_actions ) == 'read_restart_data'      .OR.   &
1017               TRIM( initializing_actions ) == 'set_constant_profiles'  .OR.   &
1018               TRIM( initializing_actions ) == 'by_user' ) )  THEN
1019             message_string = 'initializing_actions = ' //                     &
1020                              TRIM( initializing_actions ) // ' has been ' //  &
1021                              'changed to set_constant_profiles in child ' //  &
1022                              'domain.' 
1023             CALL message( 'parin', 'PA0492', 0, 0, 0, 6, 0 )
1024
1025             initializing_actions = 'set_constant_profiles'
1026          ENDIF           
1027!
1028!--       Check validity of lateral boundary conditions. This has to be done
1029!--       here because they are already used in init_pegrid and init_grid and
1030!--       therefore cannot be check in check_parameters
1031          IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
1032               bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested'  .AND. &
1033               bc_lr /= 'nesting_offline' )  THEN
1034             message_string = 'unknown boundary condition: bc_lr = "' // &
1035                              TRIM( bc_lr ) // '"'
1036             CALL message( 'parin', 'PA0049', 1, 2, 0, 6, 0 )
1037          ENDIF
1038          IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
1039               bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested'  .AND. &
1040               bc_ns /= 'nesting_offline' )  THEN
1041             message_string = 'unknown boundary condition: bc_ns = "' // &
1042                              TRIM( bc_ns ) // '"'
1043             CALL message( 'parin', 'PA0050', 1, 2, 0, 6, 0 )
1044          ENDIF
1045!
1046!--       Set internal variables used for speed optimization in if clauses
1047          IF ( bc_lr /= 'cyclic' )               bc_lr_cyc    = .FALSE.
1048          IF ( bc_lr == 'dirichlet/radiation' )  bc_lr_dirrad = .TRUE.
1049          IF ( bc_lr == 'radiation/dirichlet' )  bc_lr_raddir = .TRUE.
1050          IF ( bc_ns /= 'cyclic' )               bc_ns_cyc    = .FALSE.
1051          IF ( bc_ns == 'dirichlet/radiation' )  bc_ns_dirrad = .TRUE.
1052          IF ( bc_ns == 'radiation/dirichlet' )  bc_ns_raddir = .TRUE.
1053!
1054!--       Radiation-Dirichlet conditions are allowed along one of the horizontal directions only.
1055!--       In general, such conditions along x and y may work, but require a) some code changes
1056!--       (e.g. concerning allocation of c_u, c_v ... arrays), and b) a careful model setup by the
1057!--       user, in order to guarantee that there is a clearly defined outflow at two sides.
1058!--       Otherwise, the radiation condition may produce wrong results.
1059          IF ( ( bc_lr_dirrad .OR. bc_lr_raddir )  .AND.  ( bc_ns_dirrad .OR. bc_ns_raddir ) )  THEN
1060             message_string = 'bc_lr = "' // TRIM( bc_lr ) // '" and bc_ns = "' //                 &
1061                              TRIM( bc_ns ) // '" are not allowed to be set at the same time'
1062             CALL message( 'parin', 'PA0589', 1, 2, 0, 6, 0 )
1063          ENDIF
1064!
1065!--       Check in case of initial run, if the grid point numbers are well
1066!--       defined and allocate some arrays which are already needed in
1067!--       init_pegrid or check_parameters. During restart jobs, these arrays
1068!--       will be allocated in rrd_global. All other arrays are allocated
1069!--       in init_3d_model.
1070          IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
1071
1072             IF ( nx <= 0 )  THEN
1073                WRITE( message_string, * ) 'no value or wrong value given',    &
1074                                           ' for nx: nx=', nx
1075                CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 )
1076             ENDIF
1077             IF ( ny <= 0 )  THEN
1078                WRITE( message_string, * ) 'no value or wrong value given',    &
1079                                           ' for ny: ny=', ny
1080                CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 )
1081             ENDIF
1082             IF ( nz <= 0 )  THEN
1083                WRITE( message_string, * ) 'no value or wrong value given',    &
1084                                           ' for nz: nz=', nz
1085                CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 )
1086             ENDIF
1087
1088!
1089!--          As a side condition, routine flow_statistics require at least 14
1090!--          vertical grid levels (they are used to store time-series data)
1091!>           @todo   Remove this restriction
1092             IF ( nz < 14 )  THEN
1093                WRITE( message_string, * ) 'nz >= 14 is required'
1094                CALL message( 'parin', 'PA0362', 1, 2, 0, 6, 0 )
1095             ENDIF
1096
1097!
1098!--          ATTENTION: in case of changes to the following statement please
1099!--                  also check the allocate statement in routine rrd_global
1100             ALLOCATE( pt_init(0:nz+1), q_init(0:nz+1), s_init(0:nz+1),        &
1101                       ref_state(0:nz+1), sa_init(0:nz+1), ug(0:nz+1),         &
1102                       u_init(0:nz+1), v_init(0:nz+1), vg(0:nz+1),             &
1103                    hom(0:nz+1,2,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions),  &
1104                    hom_sum(0:nz+1,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions) )
1105
1106             hom = 0.0_wp
1107
1108          ENDIF
1109
1110!
1111!--       NAMELIST-file is not needed anymore
1112          CALL close_file( 11 )
1113
1114       ENDIF
1115#if defined( __parallel )
1116       CALL MPI_BARRIER( MPI_COMM_WORLD, ierr )
1117#endif
1118    ENDDO
1119
1120    CALL location_message( 'reading NAMELIST parameters from PARIN', 'finished' )
1121
1122 END SUBROUTINE parin
Note: See TracBrowser for help on using the repository browser.