source: palm/trunk/SOURCE/read_restart_data_mod.f90 @ 4600

Last change on this file since 4600 was 4590, checked in by suehring, 4 years ago

Biomet: Enable restart via mpi-IO; Restart: Bugfix in allocation of hom and hom_sum in case of mpi-io restart when chemistry or salsa are employed

  • Property svn:keywords set to Id
File size: 115.4 KB
Line 
1!> @file read_restart_data_mod.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-2020 Leibniz Universitaet Hannover
18!------------------------------------------------------------------------------!
19!
20! Current revisions:
21! -----------------
22!
23!
24! Former revisions:
25! -----------------
26! $Id: read_restart_data_mod.f90 4590 2020-07-06 14:34:59Z suehring $
27! Bugfix in allocation of hom and hom_sum in case of mpi-io restart when
28! chemistry or salsa are employed
29!
30! 4580 2020-06-29 07:54:21Z raasch
31! data handling with MPI-IO for cyclic-fill added (so far only for global data)
32!
33! 4564 2020-06-12 14:03:36Z raasch
34! Vertical nesting method of Huq et al. (2019) removed
35!
36! 4539 2020-05-18 14:05:17Z raasch
37! location message added
38!
39! 4536 2020-05-17 17:24:13Z raasch
40! binary version incremented
41!
42! 4534 2020-05-14 18:35:22Z raasch
43! adjustments for I/O on reduced number of cores using shared memory MPI
44!
45! 4518 2020-05-04 15:44:28Z suehring
46! Move input of diagnostic output quantities to doq_rrd_local
47!
48! 4517 2020-05-03 14:29:30Z raasch
49! qsurf and ssurf added
50!
51! 4498 2020-04-15 14:26:31Z raasch
52! argument removed from rd_mpi_io_open
53!
54! 4496 2020-04-15 08:37:26Z raasch
55! bugfix: MPI barrier removed, coupling character added to input filename
56!
57! 4495 2020-04-13 20:11:20Z raasch
58! restart data handling with MPI-IO added
59!
60! 4435 2020-03-03 10:38:41Z raasch
61! bugfix for message that reports about files that are read from in case that the virtual PE grid
62! has chenged (in case of large number of files format was exceeded), detailed messages about the
63! files are now output to the debug file
64!
65! 4431 2020-02-27 23:23:01Z gronemeier
66! added u_center_av, v_center_av, wspeed_av
67!
68! 4360 2020-01-07 11:25:50Z suehring
69! Change automatic arrays to allocatable ones in rrd_local, in order to avoid
70! memory problems due to too small stack size for large jobs with intel
71! compiler. (J.Resler)
72!
73! 4331 2019-12-10 18:25:02Z suehring
74! Enable restart data for 2-m potential temperature output
75!
76! 4301 2019-11-22 12:09:09Z oliver.maas
77! removed recycling_yshift
78!
79! 4227 2019-09-10 18:04:34Z gronemeier
80! implement new palm_date_time_mod and increased binary version
81!
82! 4146 2019-08-07 07:47:36Z gronemeier
83! Corrected "Former revisions" section
84!
85! 4131 2019-08-02 11:06:18Z monakurppa
86! Allocate hom and hom_sum to allow profile output for salsa variables.
87!
88! 4101 2019-07-17 15:14:26Z gronemeier
89! remove old_dt
90!
91! 4039 2019-06-18 10:32:41Z suehring
92! input of uu_av, vv_av, ww_av added
93!
94! 4017 2019-06-06 12:16:46Z schwenkel
95! bugfix for r3998, allocation of 3d temporary arrays of various dimensions revised
96!
97! 3998 2019-05-23 13:38:11Z suehring
98! Formatting adjustment
99!
100! 3994 2019-05-22 18:08:09Z suehring
101! output of turbulence intensity added
102!
103! 3988 2019-05-22 11:32:37Z kanani
104! + time_virtual_measurement (to enable steering of output interval)
105!
106! 3936 2019-04-26 15:38:02Z kanani
107! Enable time-averaged output of theta_2m* with restarts
108!
109! 3767 2019-02-27 08:18:02Z raasch
110! unused variables removed from rrd-subroutines parameter list
111!
112! 3766 2019-02-26 16:23:41Z raasch
113! first argument removed from module_interface_rrd_*
114!
115! 3668 2019-01-14 12:49:24Z maronga
116! Removed most_method and increased binary version
117!
118! 3655 2019-01-07 16:51:22Z knoop
119! Implementation of the PALM module interface
120!
121! 2894 2018-03-15 09:17:58Z Giersch
122! Initial revision
123!
124!
125! Description:
126! ------------
127!> Reads restart data from restart-file(s) (binary format).
128!>
129!> @todo: Revise max_pr_cs (profiles for chemistry)
130!> @todo: Modularize reading of restart data for diagnostic quantities, which
131!>        is not possible with the current module-interface structure
132!------------------------------------------------------------------------------!
133 MODULE read_restart_data_mod
134
135
136    USE arrays_3d,                                                             &
137        ONLY:  inflow_damping_factor, mean_inflow_profiles, pt_init,           &
138               q_init, ref_state, sa_init, s_init, u_init, ug, v_init, vg,     &
139               e, kh, km, p, pt, q, ql, s, u, u_m_l, u_m_n, u_m_r, u_m_s,      &
140               v, v_m_l, v_m_n, v_m_r, v_m_s, vpt, w, w_m_l, w_m_n, w_m_r, w_m_s
141
142    USE averaging
143
144    USE bulk_cloud_model_mod,                                                                      &
145        ONLY:  bulk_cloud_model
146
147    USE chem_modules,                                                                              &
148       ONLY: max_pr_cs
149
150    USE control_parameters
151
152    USE cpulog,                                                                &
153        ONLY:  cpu_log, log_point_s
154
155    USE grid_variables,                                                        &
156        ONLY:  dx, dy
157
158    USE gust_mod,                                                                                  &
159        ONLY:  gust_module_enabled
160
161    USE indices,                                                               &
162        ONLY:  nbgp, nx, nxl, nxlg, nxr, nxrg, nx_on_file, ny, nys, nysg, nyn, &
163               nyng, ny_on_file, nz, nzb, nzt
164
165    USE kinds
166
167    USE model_1d_mod,                                                          &
168        ONLY:  damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
169
170    USE module_interface,                                                      &
171        ONLY:  module_interface_rrd_global,                                    &
172               module_interface_rrd_local
173
174    USE netcdf_interface,                                                      &
175        ONLY:  netcdf_precision, output_for_t0
176
177    USE particle_attributes,                                                                       &
178        ONLY:  particle_advection
179
180    USE pegrid
181
182    USE radiation_model_mod,                                                   &
183        ONLY:  time_radiation
184
185    USE random_function_mod,                                                   &
186        ONLY:  random_iv, random_iy
187
188    USE random_generator_parallel,                                             &
189        ONLY:  id_random_array, seq_random_array
190
191    USE restart_data_mpi_io_mod,                                                                   &
192        ONLY:  rd_mpi_io_check_array, rd_mpi_io_close, rd_mpi_io_open, rrd_mpi_io,                 &
193               rrd_mpi_io_global_array
194
195    USE spectra_mod,                                                           &
196        ONLY:  average_count_sp, spectrum_x, spectrum_y
197
198    USE surface_mod,                                                           &
199        ONLY :  surface_rrd_local
200
201    USE statistics,                                                            &
202        ONLY:  statistic_regions, hom, hom_sum, pr_palm, u_max, u_max_ijk,     &
203               v_max, v_max_ijk, w_max, w_max_ijk, z_i
204
205    USE user,                                                                                      &
206        ONLY:  user_module_enabled
207
208    USE virtual_measurement_mod,                                               &
209        ONLY:  time_virtual_measurement
210
211
212    IMPLICIT NONE
213
214
215    INTERFACE rrd_global
216       MODULE PROCEDURE rrd_global
217    END INTERFACE rrd_global
218
219    INTERFACE rrd_read_parts_of_global
220       MODULE PROCEDURE rrd_read_parts_of_global
221    END INTERFACE rrd_read_parts_of_global
222
223    INTERFACE rrd_local
224       MODULE PROCEDURE rrd_local
225    END INTERFACE rrd_local
226
227    INTERFACE rrd_skip_global
228       MODULE PROCEDURE rrd_skip_global
229    END INTERFACE rrd_skip_global
230
231
232    PUBLIC rrd_global, rrd_read_parts_of_global, rrd_local, rrd_skip_global
233
234
235 CONTAINS
236
237!------------------------------------------------------------------------------!
238! Description:
239! ------------
240!> Reads values of global control variables from restart-file (binary format)
241!> created by PE0 of the previous run
242!------------------------------------------------------------------------------!
243 SUBROUTINE rrd_global
244
245
246    CHARACTER(LEN=10) ::  binary_version_global  !<
247    CHARACTER(LEN=10) ::  version_on_file        !<
248    CHARACTER(LEN=20) ::  tmp_name               !< temporary variable
249
250    INTEGER ::  i                                !< loop index
251
252    LOGICAL ::  array_found                      !<
253    LOGICAL ::  found                            !<
254
255
256    CALL location_message( 'read global restart data', 'start' )
257
258    IF ( TRIM( restart_data_format_input ) == 'fortran_binary' )  THEN
259!
260!--    Input in Fortran binary format
261       CALL check_open( 13 )
262!
263!--    Make version number check first
264       READ ( 13 )  length
265       READ ( 13 )  restart_string(1:length)
266       READ ( 13 )  version_on_file
267
268       binary_version_global = '5.0'
269       IF ( TRIM( version_on_file ) /= TRIM( binary_version_global ) )  THEN
270          WRITE( message_string, * ) 'version mismatch concerning ',           &
271                                     'binary_version_global:',                 &
272                                     '&version on file    = "',                &
273                                     TRIM( version_on_file ), '"',             &
274                                     '&version on program = "',                &
275                                     TRIM( binary_version_global ), '"'
276          CALL message( 'rrd_global', 'PA0296', 1, 2, 0, 6, 0 )
277       ENDIF
278
279!
280!--    Read number of PEs and horizontal index bounds of all PEs used in the
281!--    previous run
282       READ ( 13 )  length
283       READ ( 13 )  restart_string(1:length)
284
285       IF ( TRIM( restart_string(1:length) ) /= 'numprocs' )  THEN
286          WRITE( message_string, * ) 'numprocs not found in data from prior ', &
287                                     'run on PE ', myid
288          CALL message( 'rrd_global', 'PA0297', 1, 2, 0, 6, 0 )
289       ENDIF
290       READ ( 13 )  numprocs_previous_run
291
292       IF ( .NOT. ALLOCATED( hor_index_bounds_previous_run ) )  THEN
293          ALLOCATE( hor_index_bounds_previous_run(4,0:numprocs_previous_run-1) )
294       ENDIF
295
296       READ ( 13 )  length
297       READ ( 13 )  restart_string(1:length)
298
299       IF ( restart_string(1:length) /= 'hor_index_bounds' )  THEN
300          WRITE( message_string, * ) 'hor_index_bounds not found in data ',    &
301                                     'from prior run on PE ', myid
302          CALL message( 'rrd_global', 'PA0298', 1, 2, 0, 6, 0 )
303       ENDIF
304       READ ( 13 )  hor_index_bounds_previous_run
305
306!
307!--    Read vertical number of gridpoints and number of different areas used
308!--    for computing statistics. Allocate arrays depending on these values,
309!--    which are needed for the following read instructions.
310       READ ( 13 )  length
311       READ ( 13 )  restart_string(1:length)
312
313       IF ( restart_string(1:length) /= 'nz' )  THEN
314          WRITE( message_string, * ) 'nz not found in data from prior run ',   &
315                                     'on PE ', myid
316          CALL message( 'rrd_global', 'PA0299', 1, 2, 0, 6, 0 )
317       ENDIF
318       READ ( 13 )  nz
319
320       READ ( 13 )  length
321       READ ( 13 )  restart_string(1:length)
322
323       IF ( restart_string(1:length) /= 'max_pr_user' )  THEN
324          WRITE( message_string, * ) 'max_pr_user not found in data from ',    &
325                                     'prior run on PE ', myid
326          CALL message( 'rrd_global', 'PA0300', 1, 2, 0, 6, 0 )
327       ENDIF
328       READ ( 13 )  max_pr_user    ! This value is checked against the number of
329                                   ! user profiles given for the current run
330                                   ! in routine user_parin (it has to match)
331
332       READ ( 13 )  length
333       READ ( 13 )  restart_string(1:length)
334
335       IF ( restart_string(1:length) /= 'statistic_regions' )  THEN
336          WRITE( message_string, * ) 'statistic_regions not found in data ',   &
337                                     'from prior run on PE ', myid
338          CALL message( 'rrd_global', 'PA0301', 1, 2, 0, 6, 0 )
339       ENDIF
340       READ ( 13 )  statistic_regions
341
342!
343!--    The following global arrays (better to say, they have the same size and values on each
344!--    subdomain) are by default allocated in routine parin, but not in case of restarts!
345       IF ( .NOT. ALLOCATED( ug ) )  THEN
346          ALLOCATE( ug(0:nz+1), u_init(0:nz+1), vg(0:nz+1),                                        &
347                    v_init(0:nz+1), pt_init(0:nz+1), q_init(0:nz+1),                               &
348                    ref_state(0:nz+1), s_init(0:nz+1), sa_init(0:nz+1),                            &
349                    hom(0:nz+1,2,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions),  &
350                    hom_sum(0:nz+1,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions) )
351       ENDIF
352
353!
354!--    Now read all control parameters:
355!--    Caution: When the following read instructions have been changed, the
356!--    -------  version number stored in the variable binary_version_global has
357!--             to be increased. The same changes must also be done in
358!--             wrd_write_global.
359       READ ( 13 )  length
360       READ ( 13 )  restart_string(1:length)
361
362       DO WHILE ( restart_string(1:length) /= 'binary_version_local' )
363
364          found = .FALSE.
365
366          SELECT CASE ( restart_string(1:length) )
367
368             CASE ( 'advected_distance_x' )
369                READ ( 13 )  advected_distance_x
370             CASE ( 'advected_distance_y' )
371                READ ( 13 )  advected_distance_y
372             CASE ( 'alpha_surface' )
373                READ ( 13 )  alpha_surface
374             CASE ( 'average_count_pr' )
375                READ ( 13 )  average_count_pr
376             CASE ( 'average_count_sp' )
377                READ ( 13 )  average_count_sp
378             CASE ( 'average_count_3d' )
379                READ ( 13 )  average_count_3d
380             CASE ( 'bc_e_b' )
381                READ ( 13 )  bc_e_b
382             CASE ( 'bc_lr' )
383                READ ( 13 )  bc_lr
384             CASE ( 'bc_ns' )
385                READ ( 13 )  bc_ns
386             CASE ( 'bc_p_b' )
387                READ ( 13 )  bc_p_b
388             CASE ( 'bc_p_t' )
389                READ ( 13 )  bc_p_t
390             CASE ( 'bc_pt_b' )
391                READ ( 13 )  bc_pt_b
392             CASE ( 'bc_pt_t' )
393                READ ( 13 )  bc_pt_t
394             CASE ( 'bc_pt_t_val' )
395                READ ( 13 )  bc_pt_t_val
396             CASE ( 'bc_q_b' )
397                READ ( 13 )  bc_q_b
398             CASE ( 'bc_q_t' )
399                READ ( 13 )  bc_q_t
400             CASE ( 'bc_q_t_val' )
401                READ ( 13 )  bc_q_t_val
402             CASE ( 'bc_s_b' )
403                READ ( 13 )  bc_s_b
404             CASE ( 'bc_s_t' )
405                READ ( 13 )  bc_s_t
406             CASE ( 'bc_uv_b' )
407                READ ( 13 )  bc_uv_b
408             CASE ( 'bc_uv_t' )
409                READ ( 13 )  bc_uv_t
410             CASE ( 'building_height' )
411                READ ( 13 )  building_height
412             CASE ( 'building_length_x' )
413                READ ( 13 )  building_length_x
414             CASE ( 'building_length_y' )
415                READ ( 13 )  building_length_y
416             CASE ( 'building_wall_left' )
417                READ ( 13 )  building_wall_left
418             CASE ( 'building_wall_south' )
419                READ ( 13 )  building_wall_south
420             CASE ( 'call_psolver_at_all_substeps' )
421                READ ( 13 )  call_psolver_at_all_substeps
422             CASE ( 'canyon_height' )
423                READ ( 13 )  canyon_height
424             CASE ( 'canyon_wall_left' )
425                READ ( 13 )  canyon_wall_left
426             CASE ( 'canyon_wall_south' )
427                READ ( 13 )  canyon_wall_south
428             CASE ( 'canyon_width_x' )
429                READ ( 13 )  canyon_width_x
430             CASE ( 'canyon_width_y' )
431                READ ( 13 )  canyon_width_y
432             CASE ( 'cfl_factor' )
433                READ ( 13 )  cfl_factor
434             CASE ( 'cloud_droplets' )
435                READ ( 13 )  cloud_droplets
436             CASE ( 'collective_wait' )
437                READ ( 13 )  collective_wait
438             CASE ( 'conserve_volume_flow' )
439                READ ( 13 )  conserve_volume_flow
440             CASE ( 'conserve_volume_flow_mode' )
441                READ ( 13 )  conserve_volume_flow_mode
442             CASE ( 'constant_flux_layer' )
443                READ ( 13 )  constant_flux_layer
444             CASE ( 'coupling_start_time' )
445                READ ( 13 )  coupling_start_time
446             CASE ( 'current_timestep_number' )
447                READ ( 13 )  current_timestep_number
448             CASE ( 'cycle_mg' )
449                READ ( 13 )  cycle_mg
450             CASE ( 'damp_level_1d' )
451                READ ( 13 )  damp_level_1d
452             CASE ( 'origin_date_time' )
453                READ ( 13 )  origin_date_time
454             CASE ( 'dissipation_1d' )
455                READ ( 13 )  dissipation_1d
456             CASE ( 'do2d_xy_time_count' )
457                READ ( 13 )  do2d_xy_time_count
458             CASE ( 'do2d_xz_time_count' )
459                READ ( 13 )  do2d_xz_time_count
460             CASE ( 'do2d_yz_time_count' )
461                READ ( 13 )  do2d_yz_time_count
462             CASE ( 'do3d_time_count' )
463                READ ( 13 )  do3d_time_count
464             CASE ( 'dp_external' )
465                READ ( 13 )  dp_external
466             CASE ( 'dp_level_b' )
467                READ ( 13 )  dp_level_b
468             CASE ( 'dp_smooth' )
469                READ ( 13 )  dp_smooth
470             CASE ( 'dpdxy' )
471                READ ( 13 )  dpdxy
472             CASE ( 'dt_3d' )
473                READ ( 13 )  dt_3d
474             CASE ( 'dt_pr_1d' )
475                READ ( 13 )  dt_pr_1d
476             CASE ( 'dt_run_control_1d' )
477                READ ( 13 )  dt_run_control_1d
478             CASE ( 'dx' )
479                READ ( 13 )  dx
480             CASE ( 'dy' )
481                READ ( 13 )  dy
482             CASE ( 'dz' )
483                READ ( 13 )  dz
484             CASE ( 'dz_max' )
485                READ ( 13 )  dz_max
486             CASE ( 'dz_stretch_factor' )
487                READ ( 13 )  dz_stretch_factor
488             CASE ( 'dz_stretch_factor_array' )
489                READ ( 13 )  dz_stretch_factor_array
490             CASE ( 'dz_stretch_level' )
491                READ ( 13 )  dz_stretch_level
492             CASE ( 'dz_stretch_level_end' )
493                READ ( 13 )  dz_stretch_level_end
494             CASE ( 'dz_stretch_level_start' )
495                READ ( 13 )  dz_stretch_level_start
496             CASE ( 'e_min' )
497                READ ( 13 )  e_min
498             CASE ( 'end_time_1d' )
499                READ ( 13 )  end_time_1d
500             CASE ( 'fft_method' )
501                READ ( 13 )  fft_method
502             CASE ( 'first_call_lpm' )
503                READ ( 13 )  first_call_lpm
504             CASE ( 'galilei_transformation' )
505                READ ( 13 )  galilei_transformation
506             CASE ( 'hom' )
507                READ ( 13 )  hom
508             CASE ( 'hom_sum' )
509                READ ( 13 )  hom_sum
510             CASE ( 'humidity' )
511                READ ( 13 )  humidity
512             CASE ( 'inflow_damping_factor' )
513                IF ( .NOT. ALLOCATED( inflow_damping_factor ) )  THEN
514                   ALLOCATE( inflow_damping_factor(0:nz+1) )
515                ENDIF
516                READ ( 13 )  inflow_damping_factor
517             CASE ( 'inflow_damping_height' )
518                READ ( 13 )  inflow_damping_height
519             CASE ( 'inflow_damping_width' )
520                READ ( 13 )  inflow_damping_width
521             CASE ( 'inflow_disturbance_begin' )
522                READ ( 13 )  inflow_disturbance_begin
523             CASE ( 'inflow_disturbance_end' )
524                READ ( 13 )  inflow_disturbance_end
525             CASE ( 'km_constant' )
526                READ ( 13 )  km_constant
527             CASE ( 'large_scale_forcing' )
528                READ ( 13 )  large_scale_forcing
529             CASE ( 'large_scale_subsidence' )
530                READ ( 13 )  large_scale_subsidence
531             CASE ( 'latitude' )
532                READ ( 13 )  latitude
533             CASE ( 'longitude' )
534                READ ( 13 )  longitude
535             CASE ( 'loop_optimization' )
536                READ ( 13 )  loop_optimization
537             CASE ( 'masking_method' )
538                READ ( 13 )  masking_method
539             CASE ( 'mean_inflow_profiles' )
540                IF ( .NOT. ALLOCATED( mean_inflow_profiles ) )  THEN
541                   ALLOCATE( mean_inflow_profiles(0:nz+1,1:num_mean_inflow_profiles) )
542                ENDIF
543                READ ( 13 )  mean_inflow_profiles
544             CASE ( 'mg_cycles' )
545                READ ( 13 )  mg_cycles
546             CASE ( 'mg_switch_to_pe0_level' )
547                READ ( 13 )  mg_switch_to_pe0_level
548             CASE ( 'mixing_length_1d' )
549                READ ( 13 )  mixing_length_1d
550             CASE ( 'momentum_advec' )
551                READ ( 13 )  momentum_advec
552             CASE ( 'netcdf_precision' )
553                READ ( 13 )  netcdf_precision
554             CASE ( 'neutral' )
555                READ ( 13 )  neutral
556             CASE ( 'ngsrb' )
557                READ ( 13 )  ngsrb
558             CASE ( 'nsor' )
559                READ ( 13 )  nsor
560             CASE ( 'nsor_ini' )
561                READ ( 13 )  nsor_ini
562             CASE ( 'nudging' )
563                READ ( 13 )  nudging
564             CASE ( 'num_leg' )
565                READ ( 13 )  num_leg
566             CASE ( 'nx' )
567                READ ( 13 )  nx
568                nx_on_file = nx
569             CASE ( 'ny' )
570                READ ( 13 )  ny
571                ny_on_file = ny
572             CASE ( 'ocean_mode' )
573                READ ( 13 )  ocean_mode
574             CASE ( 'omega' )
575                READ ( 13 )  omega
576             CASE ( 'omega_sor' )
577                READ ( 13 )  omega_sor
578             CASE ( 'output_for_t0' )
579                READ (13)    output_for_t0
580             CASE ( 'passive_scalar' )
581                READ ( 13 )  passive_scalar
582             CASE ( 'prandtl_number' )
583                READ ( 13 )  prandtl_number
584             CASE ( 'psolver' )
585                READ ( 13 )  psolver
586             CASE ( 'pt_damping_factor' )
587                READ ( 13 )  pt_damping_factor
588             CASE ( 'pt_damping_width' )
589                READ ( 13 )  pt_damping_width
590             CASE ( 'pt_init' )
591                READ ( 13 )  pt_init
592             CASE ( 'pt_reference' )
593                READ ( 13 )  pt_reference
594             CASE ( 'pt_surface' )
595                READ ( 13 )  pt_surface
596             CASE ( 'pt_surface_initial_change' )
597                READ ( 13 )  pt_surface_initial_change
598             CASE ( 'pt_vertical_gradient' )
599                READ ( 13 )  pt_vertical_gradient
600             CASE ( 'pt_vertical_gradient_level' )
601                READ ( 13 )  pt_vertical_gradient_level
602             CASE ( 'pt_vertical_gradient_level_ind' )
603                READ ( 13 )  pt_vertical_gradient_level_ind
604             CASE ( 'q_init' )
605                READ ( 13 )  q_init
606             CASE ( 'q_surface' )
607                READ ( 13 )  q_surface
608             CASE ( 'q_surface_initial_change' )
609                READ ( 13 )  q_surface_initial_change
610             CASE ( 'q_vertical_gradient' )
611                READ ( 13 )  q_vertical_gradient
612             CASE ( 'q_vertical_gradient_level' )
613                READ ( 13 )  q_vertical_gradient_level
614             CASE ( 'q_vertical_gradient_level_ind' )
615                READ ( 13 )  q_vertical_gradient_level_ind
616             CASE ( 'random_generator' )
617                READ ( 13 )  random_generator
618             CASE ( 'random_heatflux' )
619                READ ( 13 )  random_heatflux
620             CASE ( 'rans_mode' )
621                READ ( 13 )  rans_mode
622             CASE ( 'rayleigh_damping_factor' )
623                READ ( 13 )  rayleigh_damping_factor
624             CASE ( 'rayleigh_damping_height' )
625                READ ( 13 )  rayleigh_damping_height
626             CASE ( 'recycling_width' )
627                READ ( 13 )  recycling_width
628             CASE ( 'ref_state' )
629                READ ( 13 )  ref_state
630             CASE ( 'reference_state' )
631                READ ( 13 )  reference_state
632             CASE ( 'residual_limit' )
633                READ ( 13 )  residual_limit
634             CASE ( 'roughness_length' )
635                READ ( 13 )  roughness_length
636             CASE ( 'run_coupled' )
637                READ ( 13 )  run_coupled
638             CASE ( 'runnr' )
639                READ ( 13 )  runnr
640             CASE ( 's_init' )
641                READ ( 13 )  s_init
642             CASE ( 's_surface' )
643                READ ( 13 )  s_surface
644             CASE ( 's_surface_initial_change' )
645                READ ( 13 )  s_surface_initial_change
646             CASE ( 's_vertical_gradient' )
647                READ ( 13 )  s_vertical_gradient
648             CASE ( 's_vertical_gradient_level' )
649                READ ( 13 )  s_vertical_gradient_level
650             CASE ( 's_vertical_gradient_level_ind' )
651                READ ( 13 )  s_vertical_gradient_level_ind
652             CASE ( 'scalar_advec' )
653                READ ( 13 )  scalar_advec
654             CASE ( 'simulated_time' )
655                READ ( 13 )  simulated_time
656             CASE ( 'spectrum_x' )
657                IF ( .NOT. ALLOCATED( spectrum_x ) )  THEN
658                   ALLOCATE( spectrum_x( 1:nx/2, 1:100, 1:10 ) )
659                ENDIF
660                READ ( 13 )  spectrum_x
661             CASE ( 'spectrum_y' )
662                IF ( .NOT. ALLOCATED( spectrum_y ) )  THEN
663                   ALLOCATE( spectrum_y( 1:ny/2, 1:100, 1:10 ) )
664                ENDIF
665                READ ( 13 )  spectrum_y
666             CASE ( 'spinup_time' )
667                READ ( 13 )  spinup_time
668             CASE ( 'subs_vertical_gradient' )
669                READ ( 13 )  subs_vertical_gradient
670             CASE ( 'subs_vertical_gradient_level' )
671                READ ( 13 )  subs_vertical_gradient_level
672             CASE ( 'subs_vertical_gradient_level_i' )
673                READ ( 13 )  subs_vertical_gradient_level_i
674             CASE ( 'surface_heatflux' )
675                READ ( 13 )  surface_heatflux
676             CASE ( 'surface_pressure' )
677                READ ( 13 )  surface_pressure
678             CASE ( 'surface_scalarflux' )
679                READ ( 13 )  surface_scalarflux
680             CASE ( 'surface_waterflux' )
681                READ ( 13 )  surface_waterflux
682             CASE ( 'time_coupling' )
683                READ ( 13 )  time_coupling
684             CASE ( 'time_disturb' )
685                READ ( 13 )  time_disturb
686             CASE ( 'time_do2d_xy' )
687                READ ( 13 )  time_do2d_xy
688             CASE ( 'time_do2d_xz' )
689                READ ( 13 )  time_do2d_xz
690             CASE ( 'time_do2d_yz' )
691                READ ( 13 )  time_do2d_yz
692             CASE ( 'time_do3d' )
693                READ ( 13 )  time_do3d
694             CASE ( 'time_do_av' )
695                READ ( 13 )  time_do_av
696             CASE ( 'time_do_sla' )
697                READ ( 13 )  time_do_sla
698             CASE ( 'time_domask' )
699                READ ( 13 )  time_domask
700             CASE ( 'time_dopr' )
701                READ ( 13 )  time_dopr
702             CASE ( 'time_dopr_av' )
703                READ ( 13 )  time_dopr_av
704             CASE ( 'time_dopr_listing' )
705                READ ( 13 )  time_dopr_listing
706             CASE ( 'time_dopts' )
707                READ ( 13 )  time_dopts
708             CASE ( 'time_dosp' )
709                READ ( 13 )  time_dosp
710             CASE ( 'time_dots' )
711                READ ( 13 )  time_dots
712             CASE ( 'time_radiation' )
713                READ ( 13 )  time_radiation
714             CASE ( 'time_restart' )
715                READ ( 13 )  time_restart
716             CASE ( 'time_run_control' )
717                READ ( 13 )  time_run_control
718             CASE ( 'time_since_reference_point' )
719                READ ( 13 )  time_since_reference_point
720             CASE ( 'time_virtual_measurement' )
721                READ ( 13 )  time_virtual_measurement
722             CASE ( 'timestep_scheme' )
723                READ ( 13 )  timestep_scheme
724             CASE ( 'top_heatflux' )
725                READ ( 13 )  top_heatflux
726             CASE ( 'top_momentumflux_u' )
727                READ ( 13 )  top_momentumflux_u
728             CASE ( 'top_momentumflux_v' )
729                READ ( 13 )  top_momentumflux_v
730             CASE ( 'top_scalarflux' )
731                READ ( 13 )  top_scalarflux
732             CASE ( 'topography' )
733                READ ( 13 )  topography
734             CASE ( 'topography_grid_convention' )
735                READ ( 13 )  topography_grid_convention
736             CASE ( 'tsc' )
737                READ ( 13 )  tsc
738             CASE ( 'tunnel_height' )
739                READ ( 13 )  tunnel_height
740             CASE ( 'tunnel_length' )
741                READ ( 13 )  tunnel_length
742             CASE ( 'tunnel_wall_depth' )
743                READ ( 13 )  tunnel_wall_depth
744             CASE ( 'tunnel_width_x' )
745                READ ( 13 )  tunnel_width_x
746             CASE ( 'tunnel_width_y' )
747                READ ( 13 )  tunnel_width_y
748             CASE ( 'turbulence_closure' )
749                READ ( 13 )  turbulence_closure
750             CASE ( 'turbulent_inflow' )
751                READ ( 13 )  turbulent_inflow
752             CASE ( 'u_bulk' )
753                READ ( 13 )  u_bulk
754             CASE ( 'u_init' )
755                READ ( 13 )  u_init
756             CASE ( 'u_max' )
757                READ ( 13 )  u_max
758             CASE ( 'u_max_ijk' )
759                READ ( 13 )  u_max_ijk
760             CASE ( 'ug' )
761                READ ( 13 )  ug
762             CASE ( 'ug_surface' )
763                READ ( 13 )  ug_surface
764             CASE ( 'ug_vertical_gradient' )
765                READ ( 13 )  ug_vertical_gradient
766             CASE ( 'ug_vertical_gradient_level' )
767                READ ( 13 )  ug_vertical_gradient_level
768             CASE ( 'ug_vertical_gradient_level_ind' )
769                READ ( 13 )  ug_vertical_gradient_level_ind
770             CASE ( 'use_surface_fluxes' )
771                READ ( 13 )  use_surface_fluxes
772             CASE ( 'use_top_fluxes' )
773                READ ( 13 )  use_top_fluxes
774             CASE ( 'use_ug_for_galilei_tr' )
775                READ ( 13 )  use_ug_for_galilei_tr
776             CASE ( 'use_upstream_for_tke' )
777                READ ( 13 )  use_upstream_for_tke
778             CASE ( 'v_bulk' )
779                READ ( 13 )  v_bulk
780             CASE ( 'v_init' )
781                READ ( 13 )  v_init
782             CASE ( 'v_max' )
783                READ ( 13 )  v_max
784             CASE ( 'v_max_ijk' )
785                READ ( 13 )  v_max_ijk
786             CASE ( 'vg' )
787                READ ( 13 )  vg
788             CASE ( 'vg_surface' )
789                READ ( 13 )  vg_surface
790             CASE ( 'vg_vertical_gradient' )
791                READ ( 13 )  vg_vertical_gradient
792             CASE ( 'vg_vertical_gradient_level' )
793                READ ( 13 )  vg_vertical_gradient_level
794             CASE ( 'vg_vertical_gradient_level_ind' )
795                READ ( 13 )  vg_vertical_gradient_level_ind
796             CASE ( 'virtual_flight' )
797                READ ( 13 )  virtual_flight
798             CASE ( 'volume_flow_area' )
799                READ ( 13 )  volume_flow_area
800             CASE ( 'volume_flow_initial' )
801                READ ( 13 )  volume_flow_initial
802             CASE ( 'w_max' )
803                READ ( 13 )  w_max
804             CASE ( 'w_max_ijk' )
805                READ ( 13 )  w_max_ijk
806             CASE ( 'wall_adjustment' )
807                READ ( 13 )  wall_adjustment
808             CASE ( 'wall_heatflux' )
809                READ ( 13 )  wall_heatflux
810             CASE ( 'wall_humidityflux' )
811                READ ( 13 )  wall_humidityflux
812             CASE ( 'wall_scalarflux' )
813                READ ( 13 )  wall_scalarflux
814             CASE ( 'y_shift' )
815                READ ( 13 )  y_shift
816             CASE ( 'z0h_factor' )
817                READ ( 13 )  z0h_factor
818             CASE ( 'zeta_max' )
819                READ ( 13 )  zeta_max
820             CASE ( 'zeta_min' )
821                READ ( 13 )  zeta_min
822             CASE ( 'z_i' )
823                READ ( 13 )  z_i
824
825             CASE DEFAULT
826!
827!--             Read global variables from of other modules
828                CALL module_interface_rrd_global( found )
829
830                IF ( .NOT. found )  THEN
831                   WRITE( message_string, * ) 'unknown variable named "',      &
832                                           restart_string(1:length),           &
833                                          '" found in global data from ',      &
834                                          'prior run on PE ', myid
835                CALL message( 'rrd_global', 'PA0302', 1, 2, 0, 6, 0 )
836
837                ENDIF
838
839          END SELECT
840!
841!--       Read next string
842          READ ( 13 )  length
843          READ ( 13 )  restart_string(1:length)
844
845       ENDDO  ! End of loop for reading the restart string
846
847       CALL close_file( 13 )
848
849    ELSEIF ( restart_data_format_input(1:3) == 'mpi' )  THEN
850!
851!--    Read global restart data using MPI-IO
852!--    ATTENTION: Arrays need to be read with routine rrd_mpi_io_global_array!
853!--    Caution: When any of the following read instructions have been changed, the
854!--    -------  version number stored in the variable binary_version_global has
855!--             to be increased. The same changes must also be done in
856!--             wrd_write_global.
857!
858!--    Open the MPI-IO restart file.
859       CALL rd_mpi_io_open( 'read', 'BININ' // TRIM( coupling_char ),                              &
860                            open_for_global_io_only = .TRUE. )
861
862!
863!--    Make version number check first
864       CALL rrd_mpi_io( 'binary_version_global',  version_on_file )
865
866       binary_version_global = '5.0'
867       IF ( TRIM( version_on_file ) /= TRIM( binary_version_global ) )  THEN
868          WRITE( message_string, * ) 'version mismatch concerning binary_version_global:',         &
869                                     '&version on file    = "', TRIM( version_on_file ), '"',      &
870                                     '&version in program = "', TRIM( binary_version_global ), '"'
871          CALL message( 'rrd_global', 'PA0296', 1, 2, 0, 6, 0 )
872       ENDIF
873
874       CALL rrd_mpi_io( 'numprocs',  numprocs_previous_run )
875       CALL rrd_mpi_io( 'nz' , nz )
876       CALL rrd_mpi_io( 'max_pr_user',  max_pr_user )
877       CALL rrd_mpi_io( 'statistic_regions', statistic_regions )
878
879!
880!--    The following global arrays (better to say, they have the same size and values on each
881!--    subdomain) are by default allocated in routine parin, but not in case of restarts!
882       IF ( .NOT. ALLOCATED( ug ) )  THEN
883           ALLOCATE( ug(0:nz+1), u_init(0:nz+1), vg(0:nz+1),                                       &
884                     v_init(0:nz+1), pt_init(0:nz+1), q_init(0:nz+1),                              &
885                     ref_state(0:nz+1), s_init(0:nz+1), sa_init(0:nz+1),                           &
886                     hom(0:nz+1,2,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions), &
887                     hom_sum(0:nz+1,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions) )
888       ENDIF
889
890       CALL rrd_mpi_io( 'advected_distance_x',  advected_distance_x )
891       CALL rrd_mpi_io( 'advected_distance_y', advected_distance_y )
892       CALL rrd_mpi_io( 'alpha_surface', alpha_surface )
893       CALL rrd_mpi_io( 'average_count_pr', average_count_pr )
894       CALL rrd_mpi_io( 'average_count_sp', average_count_sp )
895       CALL rrd_mpi_io( 'average_count_3d', average_count_3d )
896       CALL rrd_mpi_io( 'bc_e_b', bc_e_b )
897       CALL rrd_mpi_io( 'bc_lr', bc_lr )
898       CALL rrd_mpi_io( 'bc_ns', bc_ns )
899       CALL rrd_mpi_io( 'bc_p_b', bc_p_b )
900       CALL rrd_mpi_io( 'bc_p_t', bc_p_t )
901       CALL rrd_mpi_io( 'bc_pt_b', bc_pt_b )
902       CALL rrd_mpi_io( 'bc_pt_t', bc_pt_t )
903       CALL rrd_mpi_io( 'bc_pt_t_val', bc_pt_t_val )
904       CALL rrd_mpi_io( 'bc_q_b', bc_q_b )
905       CALL rrd_mpi_io( 'bc_q_t', bc_q_t )
906       CALL rrd_mpi_io( 'bc_q_t_val', bc_q_t_val )
907       CALL rrd_mpi_io( 'bc_s_b', bc_s_b )
908       CALL rrd_mpi_io( 'bc_s_t', bc_s_t )
909       CALL rrd_mpi_io( 'bc_uv_b', bc_uv_b )
910       CALL rrd_mpi_io( 'bc_uv_t', bc_uv_t )
911       CALL rrd_mpi_io( 'biometeorology', biometeorology )
912       CALL rrd_mpi_io( 'building_height', building_height )
913       CALL rrd_mpi_io( 'building_length_x', building_length_x )
914       CALL rrd_mpi_io( 'building_length_y', building_length_y )
915       CALL rrd_mpi_io( 'building_wall_left', building_wall_left )
916       CALL rrd_mpi_io( 'building_wall_south', building_wall_south )
917       CALL rrd_mpi_io( 'bulk_cloud_model', bulk_cloud_model )
918       CALL rrd_mpi_io( 'call_psolver_at_all_substeps', call_psolver_at_all_substeps )
919       CALL rrd_mpi_io( 'canyon_height', canyon_height )
920       CALL rrd_mpi_io( 'canyon_wall_left', canyon_wall_left )
921       CALL rrd_mpi_io( 'canyon_wall_south', canyon_wall_south )
922       CALL rrd_mpi_io( 'canyon_width_x',  canyon_width_x )
923       CALL rrd_mpi_io( 'canyon_width_y', canyon_width_y )
924       CALL rrd_mpi_io( 'cfl_factor', cfl_factor )
925       CALL rrd_mpi_io( 'cloud_droplets',  cloud_droplets )
926       CALL rrd_mpi_io( 'collective_wait', collective_wait )
927       CALL rrd_mpi_io( 'conserve_volume_flow', conserve_volume_flow )
928       CALL rrd_mpi_io( 'conserve_volume_flow_mode', conserve_volume_flow_mode )
929       CALL rrd_mpi_io( 'constant_flux_layer', constant_flux_layer )
930       CALL rrd_mpi_io( 'coupling_start_time', coupling_start_time )
931       CALL rrd_mpi_io( 'current_timestep_number', current_timestep_number )
932       CALL rrd_mpi_io( 'cycle_mg', cycle_mg )
933       CALL rrd_mpi_io( 'damp_level_1d', damp_level_1d )
934       CALL rrd_mpi_io( 'dissipation_1d', dissipation_1d )
935       CALL rrd_mpi_io_global_array( 'do2d_xy_time_count', do2d_xy_time_count )
936       CALL rrd_mpi_io_global_array( 'do2d_xz_time_count', do2d_xz_time_count )
937       CALL rrd_mpi_io_global_array( 'do2d_yz_time_count', do2d_yz_time_count )
938       CALL rrd_mpi_io_global_array( 'do3d_time_count', do3d_time_count )
939       CALL rrd_mpi_io( 'dp_external', dp_external )
940       CALL rrd_mpi_io( 'dp_level_b', dp_level_b )
941       CALL rrd_mpi_io( 'dp_smooth', dp_smooth )
942       CALL rrd_mpi_io_global_array( 'dpdxy', dpdxy )
943       CALL rrd_mpi_io( 'dt_3d', dt_3d )
944       CALL rrd_mpi_io( 'dt_pr_1d', dt_pr_1d )
945       CALL rrd_mpi_io( 'dt_run_control_1d', dt_run_control_1d )
946       CALL rrd_mpi_io( 'dx', dx )
947       CALL rrd_mpi_io( 'dy', dy )
948       CALL rrd_mpi_io_global_array( 'dz', dz )
949       CALL rrd_mpi_io( 'dz_max', dz_max )
950       CALL rrd_mpi_io( 'dz_stretch_factor', dz_stretch_factor )
951       CALL rrd_mpi_io_global_array( 'dz_stretch_factor_array', dz_stretch_factor_array )
952       CALL rrd_mpi_io( 'dz_stretch_level', dz_stretch_level )
953       CALL rrd_mpi_io_global_array( 'dz_stretch_level_end', dz_stretch_level_end )
954       CALL rrd_mpi_io_global_array( 'dz_stretch_level_start', dz_stretch_level_start )
955       CALL rrd_mpi_io( 'e_min', e_min )
956       CALL rrd_mpi_io( 'end_time_1d', end_time_1d )
957       CALL rrd_mpi_io( 'fft_method', fft_method )
958       CALL rrd_mpi_io( 'first_call_lpm', first_call_lpm )
959       CALL rrd_mpi_io( 'galilei_transformation', galilei_transformation )
960       CALL rrd_mpi_io( 'gust_module_enabled', gust_module_enabled )
961       CALL rrd_mpi_io_global_array( 'hom', hom )
962       CALL rrd_mpi_io_global_array( 'hom_sum', hom_sum )
963       CALL rrd_mpi_io( 'humidity', humidity )
964       CALL rd_mpi_io_check_array( 'inflow_damping_factor', found = array_found )
965       IF ( array_found )  THEN
966           IF ( .NOT. ALLOCATED( inflow_damping_factor ) )  THEN
967               ALLOCATE( inflow_damping_factor(0:nz+1) )
968           ENDIF
969           CALL rrd_mpi_io_global_array( 'inflow_damping_factor', inflow_damping_factor )
970       ENDIF
971       CALL rrd_mpi_io( 'inflow_damping_height', inflow_damping_height )
972       CALL rrd_mpi_io( 'inflow_damping_width', inflow_damping_width )
973       CALL rrd_mpi_io( 'inflow_disturbance_begin', inflow_disturbance_begin )
974       CALL rrd_mpi_io( 'inflow_disturbance_end', inflow_disturbance_end )
975       CALL rrd_mpi_io( 'km_constant', km_constant )
976       CALL rrd_mpi_io( 'large_scale_forcing', large_scale_forcing )
977       CALL rrd_mpi_io( 'large_scale_subsidence', large_scale_subsidence )
978       CALL rrd_mpi_io( 'latitude', latitude )
979       CALL rrd_mpi_io( 'longitude', longitude )
980       CALL rrd_mpi_io( 'loop_optimization', loop_optimization )
981       CALL rrd_mpi_io( 'masking_method', masking_method )
982       CALL rd_mpi_io_check_array( 'mean_inflow_profiles', found = array_found )
983       IF ( array_found)  THEN
984          IF ( .NOT. ALLOCATED( mean_inflow_profiles ) )  THEN
985             ALLOCATE( mean_inflow_profiles(0:nz+1,7) )
986          ENDIF
987          CALL rrd_mpi_io_global_array( 'mean_inflow_profiles', mean_inflow_profiles )
988       ENDIF
989       CALL rrd_mpi_io( 'mg_cycles', mg_cycles )
990       CALL rrd_mpi_io( 'mg_switch_to_pe0_level', mg_switch_to_pe0_level )
991       CALL rrd_mpi_io( 'mixing_length_1d', mixing_length_1d )
992       CALL rrd_mpi_io( 'momentum_advec', momentum_advec )
993!
994!--    There is no module procedure for CHARACTER arrays
995       DO  i = 1, SIZE( netcdf_precision , 1 )
996          WRITE( tmp_name, '(A,I2.2)' )  'netcdf_precision_', i
997          CALL rrd_mpi_io( TRIM( tmp_name ), netcdf_precision(i) )
998       ENDDO
999       CALL rrd_mpi_io( 'neutral', neutral )
1000       CALL rrd_mpi_io( 'ngsrb', ngsrb )
1001       CALL rrd_mpi_io( 'nsor', nsor )
1002       CALL rrd_mpi_io( 'nsor_ini', nsor_ini )
1003       CALL rrd_mpi_io( 'nudging', nudging )
1004       CALL rrd_mpi_io( 'num_leg', num_leg )
1005       CALL rrd_mpi_io( 'nx', nx )
1006       nx_on_file = nx
1007       CALL rrd_mpi_io( 'ny', ny )
1008       ny_on_file = ny
1009       CALL rrd_mpi_io( 'ocean_mode', ocean_mode )
1010       CALL rrd_mpi_io( 'omega', omega )
1011       CALL rrd_mpi_io( 'omega_sor', omega_sor )
1012       CALL rrd_mpi_io( 'origin_date_time', origin_date_time )
1013       CALL rrd_mpi_io( 'output_for_t0', output_for_t0 )
1014       CALL rrd_mpi_io( 'particle_advection', particle_advection )
1015       CALL rrd_mpi_io( 'passive_scalar', passive_scalar )
1016       CALL rrd_mpi_io( 'prandtl_number', prandtl_number )
1017       CALL rrd_mpi_io( 'psolver', psolver )
1018       CALL rrd_mpi_io( 'pt_damping_factor', pt_damping_factor )
1019       CALL rrd_mpi_io( 'pt_damping_width', pt_damping_width )
1020       CALL rrd_mpi_io_global_array( 'pt_init', pt_init )
1021       CALL rrd_mpi_io( 'pt_reference', pt_reference )
1022       CALL rrd_mpi_io( 'pt_surface', pt_surface )
1023       CALL rrd_mpi_io( 'pt_surface_initial_change', pt_surface_initial_change )
1024       CALL rrd_mpi_io_global_array( 'pt_vertical_gradient', pt_vertical_gradient )
1025       CALL rrd_mpi_io_global_array( 'pt_vertical_gradient_level', pt_vertical_gradient_level )
1026       CALL rrd_mpi_io_global_array( 'pt_vertical_gradient_level_ind', pt_vertical_gradient_level_ind )
1027       CALL rrd_mpi_io_global_array( 'q_init', q_init )
1028       CALL rrd_mpi_io( 'q_surface', q_surface )
1029       CALL rrd_mpi_io( 'q_surface_initial_change', q_surface_initial_change )
1030       CALL rrd_mpi_io_global_array( 'q_vertical_gradient', q_vertical_gradient )
1031       CALL rrd_mpi_io_global_array( 'q_vertical_gradient_level', q_vertical_gradient_level )
1032       CALL rrd_mpi_io_global_array( 'q_vertical_gradient_level_ind', q_vertical_gradient_level_ind )
1033       CALL rrd_mpi_io( 'random_generator', random_generator )
1034       CALL rrd_mpi_io( 'random_heatflux', random_heatflux )
1035       CALL rrd_mpi_io( 'rans_mode', rans_mode )
1036       CALL rrd_mpi_io( 'rayleigh_damping_factor', rayleigh_damping_factor )
1037       CALL rrd_mpi_io( 'rayleigh_damping_height', rayleigh_damping_height )
1038       CALL rrd_mpi_io( 'recycling_width', recycling_width )
1039       CALL rrd_mpi_io_global_array( 'ref_state', ref_state )
1040       CALL rrd_mpi_io( 'reference_state', reference_state )
1041       CALL rrd_mpi_io( 'residual_limit', residual_limit )
1042       CALL rrd_mpi_io( 'roughness_length', roughness_length )
1043       CALL rrd_mpi_io( 'run_coupled', run_coupled )
1044       CALL rrd_mpi_io( 'runnr', runnr )
1045       CALL rrd_mpi_io_global_array( 's_init', s_init )
1046       CALL rrd_mpi_io( 's_surface', s_surface )
1047       CALL rrd_mpi_io( 's_surface_initial_change', s_surface_initial_change )
1048       CALL rrd_mpi_io_global_array( 's_vertical_gradient', s_vertical_gradient )
1049       CALL rrd_mpi_io_global_array( 's_vertical_gradient_level', s_vertical_gradient_level )
1050       CALL rrd_mpi_io_global_array( 's_vertical_gradient_level_ind', s_vertical_gradient_level_ind )
1051       CALL rrd_mpi_io( 'scalar_advec', scalar_advec )
1052       CALL rrd_mpi_io( 'simulated_time', simulated_time )
1053       CALL rd_mpi_io_check_array( 'spectrum_x', found = array_found )
1054       IF (array_found )  THEN
1055           IF ( .NOT. ALLOCATED( spectrum_x ) )  THEN
1056              ALLOCATE( spectrum_x( 1:nx/2, 1:100, 1:10 ) )
1057           ENDIF
1058           CALL rrd_mpi_io( 'spectrum_y', spectrum_y )
1059       ENDIF
1060       CALL rd_mpi_io_check_array( 'spectrum_y', found = array_found )
1061       IF ( array_found )  THEN
1062           IF ( .NOT. ALLOCATED( spectrum_y ) )  THEN
1063              ALLOCATE( spectrum_y( 1:ny/2, 1:100, 1:10 ) )
1064           ENDIF
1065           CALL rrd_mpi_io( 'spectrum_y', spectrum_y )
1066       ENDIF
1067       CALL rrd_mpi_io( 'spinup_time ', spinup_time )
1068       CALL rrd_mpi_io_global_array( 'subs_vertical_gradient', subs_vertical_gradient )
1069       CALL rrd_mpi_io_global_array( 'subs_vertical_gradient_level', subs_vertical_gradient_level )
1070       CALL rrd_mpi_io_global_array( 'subs_vertical_gradient_level_i', subs_vertical_gradient_level_i )
1071       CALL rrd_mpi_io( 'surface_heatflux', surface_heatflux )
1072       CALL rrd_mpi_io( 'surface_pressure', surface_pressure )
1073       CALL rrd_mpi_io( 'surface_output', surface_output )
1074       CALL rrd_mpi_io( 'surface_scalarflux', surface_scalarflux )
1075       CALL rrd_mpi_io( 'surface_waterflux', surface_waterflux )
1076       CALL rrd_mpi_io( 'syn_turb_gen', syn_turb_gen )
1077       CALL rrd_mpi_io( 'time_coupling', time_coupling )
1078       CALL rrd_mpi_io( 'time_disturb', time_disturb )
1079       CALL rrd_mpi_io( 'time_do2d_xy', time_do2d_xy )
1080       CALL rrd_mpi_io( 'time_do2d_xz', time_do2d_xz )
1081       CALL rrd_mpi_io( 'time_do2d_yz', time_do2d_yz )
1082       CALL rrd_mpi_io( 'time_do3d', time_do3d )
1083       CALL rrd_mpi_io( 'time_do_av', time_do_av )
1084       CALL rrd_mpi_io( 'time_do_sla', time_do_sla )
1085       CALL rrd_mpi_io_global_array( 'time_domask', time_domask )
1086       CALL rrd_mpi_io( 'time_dopr', time_dopr )
1087       CALL rrd_mpi_io( 'time_dopr_av', time_dopr_av )
1088       CALL rrd_mpi_io( 'time_dopr_listing', time_dopr_listing )
1089       CALL rrd_mpi_io( 'time_dopts', time_dopts )
1090       CALL rrd_mpi_io( 'time_dosp', time_dosp )
1091       CALL rrd_mpi_io( 'time_dots', time_dots )
1092       CALL rrd_mpi_io( 'time_radiation', time_radiation )
1093       CALL rrd_mpi_io( 'time_restart', time_restart )
1094       CALL rrd_mpi_io( 'time_run_control', time_run_control )
1095       CALL rrd_mpi_io( 'time_since_reference_point', time_since_reference_point )
1096       CALL rrd_mpi_io( 'time_virtual_measurement', time_virtual_measurement )
1097       CALL rrd_mpi_io( 'timestep_scheme', timestep_scheme )
1098       CALL rrd_mpi_io( 'top_heatflux', top_heatflux )
1099       CALL rrd_mpi_io( 'top_momentumflux_u', top_momentumflux_u )
1100       CALL rrd_mpi_io( 'top_momentumflux_v', top_momentumflux_v )
1101       CALL rrd_mpi_io( 'top_scalarflux', top_scalarflux )
1102       CALL rrd_mpi_io( 'topography', topography )
1103       CALL rrd_mpi_io( 'topography_grid_convention', topography_grid_convention )
1104       CALL rrd_mpi_io_global_array( 'tsc', tsc )
1105       CALL rrd_mpi_io( 'tunnel_height', tunnel_height )
1106       CALL rrd_mpi_io( 'tunnel_length', tunnel_length )
1107       CALL rrd_mpi_io( 'tunnel_wall_depth', tunnel_wall_depth )
1108       CALL rrd_mpi_io( 'tunnel_width_x', tunnel_width_x )
1109       CALL rrd_mpi_io( 'tunnel_width_y', tunnel_width_y )
1110       CALL rrd_mpi_io( 'turbulence_closure', turbulence_closure )
1111       CALL rrd_mpi_io( 'turbulent_inflow', turbulent_inflow )
1112       CALL rrd_mpi_io( 'u_bulk', u_bulk )
1113       CALL rrd_mpi_io_global_array( 'u_init', u_init )
1114       CALL rrd_mpi_io( 'u_max', u_max )
1115       CALL rrd_mpi_io_global_array( 'u_max_ijk', u_max_ijk )
1116       CALL rrd_mpi_io_global_array( 'ug', ug )
1117       CALL rrd_mpi_io( 'ug_surface', ug_surface )
1118       CALL rrd_mpi_io_global_array( 'ug_vertical_gradient', ug_vertical_gradient )
1119       CALL rrd_mpi_io_global_array( 'ug_vertical_gradient_level', ug_vertical_gradient_level )
1120       CALL rrd_mpi_io_global_array( 'ug_vertical_gradient_level_ind', ug_vertical_gradient_level_ind )
1121       CALL rrd_mpi_io( 'use_surface_fluxes', use_surface_fluxes )
1122       CALL rrd_mpi_io( 'use_top_fluxes', use_top_fluxes )
1123       CALL rrd_mpi_io( 'use_ug_for_galilei_tr', use_ug_for_galilei_tr )
1124       CALL rrd_mpi_io( 'use_upstream_for_tke', use_upstream_for_tke )
1125       CALL rrd_mpi_io( 'user_module_enabled', user_module_enabled )
1126       CALL rrd_mpi_io( 'v_bulk', v_bulk )
1127       CALL rrd_mpi_io_global_array( 'v_init', v_init )
1128       CALL rrd_mpi_io( 'v_max', v_max )
1129       CALL rrd_mpi_io_global_array( 'v_max_ijk', v_max_ijk )
1130       CALL rrd_mpi_io_global_array( 'vg', vg )
1131       CALL rrd_mpi_io( 'vg_surface', vg_surface )
1132       CALL rrd_mpi_io_global_array( 'vg_vertical_gradient', vg_vertical_gradient )
1133       CALL rrd_mpi_io_global_array( 'vg_vertical_gradient_level', vg_vertical_gradient_level )
1134       CALL rrd_mpi_io_global_array( 'vg_vertical_gradient_level_ind', vg_vertical_gradient_level_ind )
1135       CALL rrd_mpi_io( 'virtual_flight', virtual_flight )
1136       CALL rrd_mpi_io_global_array( 'volume_flow_area', volume_flow_area )
1137       CALL rrd_mpi_io_global_array( 'volume_flow_initial', volume_flow_initial )
1138       CALL rrd_mpi_io( 'w_max', w_max )
1139       CALL rrd_mpi_io_global_array( 'w_max_ijk', w_max_ijk )
1140       CALL rrd_mpi_io( 'wall_adjustment', wall_adjustment )
1141       CALL rrd_mpi_io_global_array( 'wall_heatflux', wall_heatflux )
1142       CALL rrd_mpi_io_global_array( 'wall_humidityflux', wall_humidityflux )
1143       CALL rrd_mpi_io_global_array( 'wall_scalarflux', wall_scalarflux )
1144       CALL rrd_mpi_io( 'y_shift', y_shift )
1145       CALL rrd_mpi_io( 'z0h_factor', z0h_factor )
1146       CALL rrd_mpi_io( 'zeta_max', zeta_max )
1147       CALL rrd_mpi_io( 'zeta_min', zeta_min )
1148       CALL rrd_mpi_io_global_array( 'z_i', z_i )
1149
1150!
1151!--    Read global variables from of other modules
1152       CALL module_interface_rrd_global
1153
1154!
1155!--    Close restart file
1156       CALL rd_mpi_io_close
1157
1158    ENDIF
1159
1160    CALL location_message( 'read global restart data', 'finished' )
1161
1162 END SUBROUTINE rrd_global
1163
1164
1165
1166!------------------------------------------------------------------------------!
1167! Description:
1168! ------------
1169!> Skipping the global control variables from restart-file (binary format)
1170!> except some information needed when reading restart data from a previous
1171!> run which used a smaller total domain or/and a different domain decomposition
1172!> (initializing_actions  == 'cyclic_fill').
1173!------------------------------------------------------------------------------!
1174 SUBROUTINE rrd_read_parts_of_global
1175
1176
1177    CHARACTER (LEN=10) ::  version_on_file
1178    CHARACTER (LEN=20) ::  momentum_advec_check
1179    CHARACTER (LEN=20) ::  scalar_advec_check
1180    CHARACTER (LEN=1)  ::  cdum
1181
1182    INTEGER(iwp) ::  max_pr_user_on_file
1183    INTEGER(iwp) ::  nz_on_file
1184    INTEGER(iwp) ::  statistic_regions_on_file
1185    INTEGER(iwp) ::  tmp_mpru
1186    INTEGER(iwp) ::  tmp_sr
1187
1188    REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE ::  hom_sum_on_file
1189    REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE ::  hom_on_file
1190
1191
1192    IF ( TRIM( restart_data_format_input ) == 'fortran_binary' )  THEN
1193!
1194!--    Input in Fortran binary format
1195       CALL check_open( 13 )
1196
1197       READ ( 13 )  length
1198       READ ( 13 )  restart_string(1:length)
1199       READ ( 13 )  version_on_file
1200
1201!
1202!--    Read number of PEs and horizontal index bounds of all PEs used in previous run
1203       READ ( 13 )  length
1204       READ ( 13 )  restart_string(1:length)
1205
1206       IF ( restart_string(1:length) /= 'numprocs' )  THEN
1207          WRITE( message_string, * ) 'numprocs not found in data from prior ', &
1208                                     'run on PE ', myid
1209          CALL message( 'rrd_read_parts_of_global', 'PA0297', 1, 2, 0, 6, 0 )
1210       ENDIF
1211       READ ( 13 )  numprocs_previous_run
1212
1213       IF ( .NOT. ALLOCATED( hor_index_bounds_previous_run ) )  THEN
1214          ALLOCATE( hor_index_bounds_previous_run(4,0:numprocs_previous_run-1) )
1215       ENDIF
1216
1217       READ ( 13 )  length
1218       READ ( 13 )  restart_string(1:length)
1219
1220       IF ( restart_string(1:length) /= 'hor_index_bounds' )  THEN
1221          WRITE( message_string, * ) 'hor_index_bounds not found in data ',    &
1222                                     'from prior run on PE ', myid
1223          CALL message( 'rrd_read_parts_of_global', 'PA0298', 1, 2, 0, 6, 0 )
1224       ENDIF
1225       READ ( 13 )  hor_index_bounds_previous_run
1226
1227!
1228!--    Read vertical number of gridpoints and number of different areas used for computing
1229!--    statistics. Allocate arrays depending on these values, which are needed for the following
1230!--    read instructions.
1231       READ ( 13 )  length
1232       READ ( 13 )  restart_string(1:length)
1233
1234       IF ( restart_string(1:length) /= 'nz' )  THEN
1235          message_string = 'nz not found in restart data file'
1236          CALL message( 'rrd_read_parts_of_global', 'PA0303', 1, 2, 0, 6, 0 )
1237       ENDIF
1238       READ ( 13 )  nz_on_file
1239       IF ( nz_on_file /= nz )  THEN
1240          WRITE( message_string, * ) 'mismatch concerning number of ',         &
1241                                     'gridpoints along z:',                    &
1242                                     '&nz on file    = "', nz_on_file, '"',    &
1243                                     '&nz from run   = "', nz, '"'
1244          CALL message( 'rrd_read_parts_of_global', 'PA0304', 1, 2, 0, 6, 0 )
1245       ENDIF
1246
1247       READ ( 13 )  length
1248       READ ( 13 )  restart_string(1:length)
1249
1250       IF ( restart_string(1:length) /= 'max_pr_user' )  THEN
1251          message_string = 'max_pr_user not found in restart data file'
1252          CALL message( 'rrd_read_parts_of_global', 'PA0305', 1, 2, 0, 6, 0 )
1253       ENDIF
1254       READ ( 13 )  max_pr_user_on_file
1255       IF ( max_pr_user_on_file /= max_pr_user )  THEN
1256          WRITE( message_string, * ) 'number of user profiles on res',         &
1257                                     'tart data file differs from the ',       &
1258                                     'current run:&max_pr_user on file    = "',&
1259                                     max_pr_user_on_file, '"',                 &
1260                                     '&max_pr_user from run   = "',            &
1261                                     max_pr_user, '"'
1262          CALL message( 'rrd_read_parts_of_global', 'PA0306', 0, 0, 0, 6, 0 )
1263          tmp_mpru = MIN( max_pr_user_on_file, max_pr_user )
1264       ELSE
1265          tmp_mpru = max_pr_user
1266       ENDIF
1267
1268       READ ( 13 )  length
1269       READ ( 13 )  restart_string(1:length)
1270
1271       IF ( restart_string(1:length) /= 'statistic_regions' )  THEN
1272          message_string = 'statistic_regions not found in restart data file'
1273          CALL message( 'rrd_read_parts_of_global', 'PA0307', 1, 2, 0, 6, 0 )
1274       ENDIF
1275       READ ( 13 )  statistic_regions_on_file
1276       IF ( statistic_regions_on_file /= statistic_regions )  THEN
1277          WRITE( message_string, * ) 'statistic regions on restart data file ',&
1278                                     'differ from the current run:',           &
1279                                     '&statistic regions on file    = "',      &
1280                                     statistic_regions_on_file, '"',           &
1281                                     '&statistic regions from run   = "',      &
1282                                      statistic_regions, '"',                  &
1283                                     '&statistic data may be lost!'
1284          CALL message( 'rrd_read_parts_of_global', 'PA0308', 0, 1, 0, 6, 0 )
1285          tmp_sr = MIN( statistic_regions_on_file, statistic_regions )
1286       ELSE
1287          tmp_sr = statistic_regions
1288       ENDIF
1289
1290!
1291!--    Now read and check some control parameters and skip the rest
1292       READ ( 13 )  length
1293       READ ( 13 )  restart_string(1:length)
1294
1295       DO  WHILE ( restart_string(1:length) /= 'binary_version_local' )
1296
1297          SELECT CASE ( restart_string(1:length) )
1298
1299             CASE ( 'average_count_pr' )
1300                READ ( 13 )  average_count_pr
1301                IF ( average_count_pr /= 0 )  THEN
1302                   WRITE( message_string, * ) 'inflow profiles not ',          &
1303                                  'temporally averaged. &Averaging will be ',  &
1304                                  'done now using', average_count_pr,          &
1305                                  ' samples.'
1306                   CALL message( 'rrd_read_parts_of_global', 'PA0309',         &
1307                                 0, 1, 0, 6, 0 )
1308                ENDIF
1309
1310             CASE ( 'hom' )
1311                ALLOCATE( hom_on_file(0:nz+1,2,pr_palm+max_pr_user_on_file,    &
1312                          0:statistic_regions_on_file) )
1313                READ ( 13 )  hom_on_file
1314                hom(:,:,1:pr_palm+tmp_mpru,0:tmp_sr) =                         &
1315                             hom_on_file(:,:,1:pr_palm+tmp_mpru,0:tmp_sr)
1316                DEALLOCATE( hom_on_file )
1317
1318             CASE ( 'hom_sum' )
1319                ALLOCATE( hom_sum_on_file(0:nz+1,pr_palm+max_pr_user_on_file,  &
1320                          0:statistic_regions_on_file) )
1321                READ ( 13 )  hom_sum_on_file
1322                hom_sum(:,1:pr_palm+tmp_mpru,0:tmp_sr) =                       &
1323                             hom_sum_on_file(:,1:pr_palm+tmp_mpru,0:tmp_sr)
1324                DEALLOCATE( hom_sum_on_file )
1325
1326             CASE ( 'momentum_advec' )
1327                momentum_advec_check = momentum_advec
1328                READ ( 13 )  momentum_advec
1329                IF ( TRIM( momentum_advec_check ) /= TRIM( momentum_advec ) )  &
1330                THEN
1331                   WRITE( message_string, * ) 'momentum_advec of the restart ',&
1332                                  'run differs from momentum_advec of the ',   &
1333                                  'initial run.'
1334                   CALL message( 'rrd_read_parts_of_global', 'PA0100',         &
1335                                 1, 2, 0, 6, 0 )
1336                ENDIF
1337
1338             CASE ( 'nx' )
1339                READ ( 13 )  nx_on_file
1340
1341             CASE ( 'ny' )
1342                READ ( 13 )  ny_on_file
1343
1344             CASE ( 'ref_state' )
1345                READ ( 13 )  ref_state
1346
1347             CASE ( 'scalar_advec' )
1348                scalar_advec_check = scalar_advec
1349                READ ( 13 )  scalar_advec
1350                IF ( TRIM( scalar_advec_check ) /= TRIM( scalar_advec ) )      &
1351                THEN
1352                   WRITE( message_string, * ) 'scalar_advec of the restart ',  &
1353                                  'run differs from scalar_advec of the ',     &
1354                                  'initial run.'
1355                   CALL message( 'rrd_read_parts_of_global', 'PA0101',         &
1356                                 1, 2, 0, 6, 0 )
1357                ENDIF
1358
1359             CASE DEFAULT
1360
1361                READ ( 13 )  cdum
1362
1363          END SELECT
1364
1365          READ ( 13 )  length
1366          READ ( 13 )  restart_string(1:length)
1367
1368       ENDDO
1369
1370       CALL close_file( 13 )
1371
1372    ELSEIF ( restart_data_format_input(1:3) == 'mpi' )  THEN
1373!
1374!--    Open the MPI-IO restart file.
1375       CALL rd_mpi_io_open( 'read', 'BININ' // TRIM( coupling_char ),                              &
1376                            open_for_global_io_only = .TRUE. )
1377
1378!
1379!--    Read vertical number of gridpoints and number of different areas used for computing
1380!--    statistics. Allocate arrays depending on these values, which are required for the following
1381!--    read instructions.
1382       CALL rrd_mpi_io( 'nz', nz_on_file )
1383       IF ( nz_on_file /= nz )  THEN
1384          WRITE( message_string, * ) 'mismatch concerning number of gridpoints along z:',          &
1385                                     '&nz on file    = "', nz_on_file, '"',                        &
1386                                     '&nz from run   = "', nz, '"'
1387          CALL message( 'rrd_read_parts_of_global', 'PA0304', 1, 2, 0, 6, 0 )
1388       ENDIF
1389
1390       CALL rrd_mpi_io( 'max_pr_user', max_pr_user_on_file )
1391       IF ( max_pr_user_on_file /= max_pr_user )  THEN
1392          WRITE( message_string, * ) 'number of user profiles on restart data file differs from ', &
1393                                     'the current run:&max_pr_user on file    = "',                &
1394                                     max_pr_user_on_file, '" &max_pr_user from run   = "',         &
1395                                     max_pr_user, '"'
1396          CALL message( 'rrd_read_parts_of_global', 'PA0306', 0, 0, 0, 6, 0 )
1397          tmp_mpru = MIN( max_pr_user_on_file, max_pr_user )
1398       ELSE
1399          tmp_mpru = max_pr_user
1400       ENDIF
1401
1402       CALL rrd_mpi_io( 'statistic_regions', statistic_regions_on_file )
1403       IF ( statistic_regions_on_file /= statistic_regions )  THEN
1404          WRITE( message_string, * ) 'statistic regions on restart data file ',&
1405                                     'differ from the current run:',           &
1406                                     '&statistic regions on file    = "',      &
1407                                     statistic_regions_on_file, '"',           &
1408                                     '&statistic regions from run   = "',      &
1409                                      statistic_regions, '"',                  &
1410                                     '&statistic data may be lost!'
1411          CALL message( 'rrd_read_parts_of_global', 'PA0308', 0, 1, 0, 6, 0 )
1412          tmp_sr = MIN( statistic_regions_on_file, statistic_regions )
1413       ELSE
1414          tmp_sr = statistic_regions
1415       ENDIF
1416
1417!
1418!--    Now read and check some control parameters and skip the rest.
1419       CALL rrd_mpi_io( 'average_count_pr', average_count_pr )
1420       IF ( average_count_pr /= 0 )  THEN
1421          WRITE( message_string, * ) 'inflow profiles not ',          &
1422                         'temporally averaged. &Averaging will be ',  &
1423                         'done now using', average_count_pr,          &
1424                         ' samples.'
1425          CALL message( 'rrd_read_parts_of_global', 'PA0309', 0, 1, 0, 6, 0 )
1426       ENDIF
1427
1428       ALLOCATE( hom_on_file(0:nz+1,2,pr_palm+max_pr_user_on_file,0:statistic_regions_on_file) )
1429       CALL rrd_mpi_io_global_array( 'hom', hom_on_file )
1430       hom(:,:,1:pr_palm+tmp_mpru,0:tmp_sr) = hom_on_file(:,:,1:pr_palm+tmp_mpru,0:tmp_sr)
1431       DEALLOCATE( hom_on_file )
1432
1433       ALLOCATE( hom_sum_on_file(0:nz+1,pr_palm+max_pr_user_on_file, 0:statistic_regions_on_file) )
1434       CALL rrd_mpi_io_global_array( 'hom_sum', hom_sum_on_file )
1435       hom_sum(:,1:pr_palm+tmp_mpru,0:tmp_sr) = hom_sum_on_file(:,1:pr_palm+tmp_mpru,0:tmp_sr)
1436       DEALLOCATE( hom_sum_on_file )
1437
1438       momentum_advec_check = momentum_advec
1439       CALL rrd_mpi_io( 'momentum_advec', momentum_advec )
1440       IF ( TRIM( momentum_advec_check ) /= TRIM( momentum_advec ) )  THEN
1441          WRITE( message_string, * ) 'momentum_advec of the restart ',&
1442                                  'run differs from momentum_advec of the ',   &
1443                                  'initial run.'
1444          CALL message( 'rrd_read_parts_of_global', 'PA0100', 1, 2, 0, 6, 0 )
1445       ENDIF
1446
1447       CALL rrd_mpi_io( 'nx', nx_on_file )
1448       CALL rrd_mpi_io( 'ny', ny_on_file )
1449       CALL rrd_mpi_io_global_array( 'ref_state', ref_state )
1450
1451       scalar_advec_check = scalar_advec
1452       CALL rrd_mpi_io( 'scalar_advec', scalar_advec )
1453       IF ( TRIM( scalar_advec_check ) /= TRIM( scalar_advec ) )  THEN
1454          WRITE( message_string, * ) 'scalar_advec of the restart ',  &
1455                                  'run differs from scalar_advec of the ',     &
1456                                  'initial run.'
1457          CALL message( 'rrd_read_parts_of_global', 'PA0101', 1, 2, 0, 6, 0 )
1458       ENDIF
1459
1460!
1461!--    Close restart file
1462       CALL rd_mpi_io_close
1463
1464    ENDIF
1465
1466!
1467!-- Calculate the temporal average of vertical profiles, if neccessary
1468    IF ( average_count_pr /= 0 )  THEN
1469       hom_sum = hom_sum / REAL( average_count_pr, KIND=wp )
1470    ENDIF
1471
1472 END SUBROUTINE rrd_read_parts_of_global
1473
1474
1475! Description:
1476! ------------
1477!> Reads processor (subdomain) specific data of variables and arrays from restart file
1478!> (binary format).
1479!------------------------------------------------------------------------------!
1480 SUBROUTINE rrd_local
1481
1482
1483    CHARACTER (LEN=7)  ::  myid_char_save
1484    CHARACTER (LEN=10) ::  binary_version_local
1485    CHARACTER (LEN=10) ::  version_on_file
1486    CHARACTER (LEN=20) ::  tmp_name               !< temporary variable
1487
1488    INTEGER(iwp) ::  files_to_be_opened  !<
1489    INTEGER(iwp) ::  i                   !<
1490    INTEGER(iwp) ::  j                   !<
1491    INTEGER(iwp) ::  k                   !<
1492    INTEGER(iwp) ::  myid_on_file        !<
1493    INTEGER(iwp) ::  numprocs_on_file    !<
1494    INTEGER(iwp) ::  nxlc                !<
1495    INTEGER(iwp) ::  nxlf                !<
1496    INTEGER(iwp) ::  nxlpr               !<
1497    INTEGER(iwp) ::  nxl_on_file         !<
1498    INTEGER(iwp) ::  nxrc                !<
1499    INTEGER(iwp) ::  nxrf                !<
1500    INTEGER(iwp) ::  nxrpr               !<
1501    INTEGER(iwp) ::  nxr_on_file         !<
1502    INTEGER(iwp) ::  nync                !<
1503    INTEGER(iwp) ::  nynf                !<
1504    INTEGER(iwp) ::  nynpr               !<
1505    INTEGER(iwp) ::  nyn_on_file         !<
1506    INTEGER(iwp) ::  nysc                !<
1507    INTEGER(iwp) ::  nysf                !<
1508    INTEGER(iwp) ::  nyspr               !<
1509    INTEGER(iwp) ::  nys_on_file         !<
1510    INTEGER(iwp) ::  nzb_on_file         !<
1511    INTEGER(iwp) ::  nzt_on_file         !<
1512    INTEGER(iwp) ::  offset_x            !<
1513    INTEGER(iwp) ::  offset_y            !<
1514    INTEGER(iwp) ::  shift_x             !<
1515    INTEGER(iwp) ::  shift_y             !<
1516
1517    INTEGER(iwp), DIMENSION(numprocs_previous_run) ::  file_list       !<
1518    INTEGER(iwp), DIMENSION(numprocs_previous_run) ::  overlap_count   !<
1519
1520    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nxlfa      !<
1521    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nxrfa      !<
1522    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nynfa      !<
1523    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nysfa      !<
1524    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  offset_xa  !<
1525    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  offset_ya  !<
1526
1527    INTEGER(isp), DIMENSION(:,:),   ALLOCATABLE ::  tmp_2d_id_random   !< temporary array for storing random generator data
1528    INTEGER(isp), DIMENSION(:,:,:), ALLOCATABLE ::  tmp_2d_seq_random  !< temporary array for storing random generator data
1529
1530    LOGICAL ::  array_found                      !<
1531    LOGICAL ::  found                            !<
1532
1533    REAL(wp), DIMENSION(:,:),   ALLOCATABLE   ::  tmp_2d         !< temporary array for storing 2D data
1534    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE   ::  tmp_3d         !< temporary array for storing 3D data
1535    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE   ::  tmp_3d_non_standard !< temporary array for storing 3D data
1536                                                                 !< with non standard dimensions
1537
1538!
1539!-- Read data from previous model run.
1540    CALL cpu_log( log_point_s(14), 'read-restart-data-local', 'start' )
1541
1542    CALL location_message( 'reading local restart data', 'start' )
1543
1544    IF ( TRIM( restart_data_format_input ) == 'fortran_binary' )  THEN
1545!
1546!--    Input in Fortran binary format
1547
1548!
1549!--    Allocate temporary buffer arrays. In previous versions, they were
1550!--    declared as automated arrays, causing memory problems when these
1551!--    were allocate on stack.
1552       ALLOCATE( nxlfa(numprocs_previous_run,1000) )
1553       ALLOCATE( nxrfa(numprocs_previous_run,1000) )
1554       ALLOCATE( nynfa(numprocs_previous_run,1000) )
1555       ALLOCATE( nysfa(numprocs_previous_run,1000) )
1556       ALLOCATE( offset_xa(numprocs_previous_run,1000) )
1557       ALLOCATE( offset_ya(numprocs_previous_run,1000) )
1558
1559!
1560!--    Check which of the restart files contain data needed for the subdomain
1561!--    of this PE
1562       files_to_be_opened = 0
1563
1564       DO  i = 1, numprocs_previous_run
1565!
1566!--       Store array bounds of the previous run ("pr") in temporary scalars
1567          nxlpr = hor_index_bounds_previous_run(1,i-1)
1568          nxrpr = hor_index_bounds_previous_run(2,i-1)
1569          nyspr = hor_index_bounds_previous_run(3,i-1)
1570          nynpr = hor_index_bounds_previous_run(4,i-1)
1571
1572!
1573!--       Determine the offsets. They may be non-zero in case that the total domain
1574!--       on file is smaller than the current total domain.
1575          offset_x = ( nxl / ( nx_on_file + 1 ) ) * ( nx_on_file + 1 )
1576          offset_y = ( nys / ( ny_on_file + 1 ) ) * ( ny_on_file + 1 )
1577
1578!
1579!--       Start with this offset and then check, if the subdomain on file
1580!--       matches another time(s) in the current subdomain by shifting it
1581!--       for nx_on_file+1, ny_on_file+1 respectively
1582
1583          shift_y = 0
1584          j       = 0
1585          DO WHILE (  nyspr+shift_y <= nyn-offset_y )
1586
1587             IF ( nynpr+shift_y >= nys-offset_y ) THEN
1588
1589                shift_x = 0
1590                DO WHILE ( nxlpr+shift_x <= nxr-offset_x )
1591
1592                   IF ( nxrpr+shift_x >= nxl-offset_x ) THEN
1593                      j = j +1
1594                      IF ( j > 1000 )  THEN
1595!
1596!--                      Array bound exceeded
1597                         message_string = 'data from subdomain of previous' //                        &
1598                                          ' run mapped more than 1000 times'
1599                         CALL message( 'rrd_local', 'PA0284', 2, 2, -1, 6, 1 )
1600                      ENDIF
1601
1602                      IF ( j == 1 )  THEN
1603                         files_to_be_opened = files_to_be_opened + 1
1604                         file_list(files_to_be_opened) = i-1
1605                      ENDIF
1606
1607                      offset_xa(files_to_be_opened,j) = offset_x + shift_x
1608                      offset_ya(files_to_be_opened,j) = offset_y + shift_y
1609!
1610!--                   Index bounds of overlapping data
1611                      nxlfa(files_to_be_opened,j) = MAX( nxl-offset_x-shift_x, nxlpr )
1612                      nxrfa(files_to_be_opened,j) = MIN( nxr-offset_x-shift_x, nxrpr )
1613                      nysfa(files_to_be_opened,j) = MAX( nys-offset_y-shift_y, nyspr )
1614                      nynfa(files_to_be_opened,j) = MIN( nyn-offset_y-shift_y, nynpr )
1615
1616                   ENDIF
1617
1618                   shift_x = shift_x + ( nx_on_file + 1 )
1619                ENDDO
1620
1621             ENDIF
1622
1623             shift_y = shift_y + ( ny_on_file + 1 )
1624          ENDDO
1625
1626          IF ( j > 0 )  overlap_count(files_to_be_opened) = j
1627
1628       ENDDO
1629
1630!
1631!--    Save the id-string of the current process, since myid_char may now be used
1632!--    to open files created by PEs with other id.
1633       myid_char_save = myid_char
1634
1635       IF ( files_to_be_opened /= 1  .OR.  numprocs /= numprocs_previous_run )  THEN
1636          WRITE( message_string, * ) 'number of PEs or virtual PE-grid changed in restart run. & ',   &
1637                                     'Set debug_output =.T. to get a list of files from which the & ',&
1638                                     'single PEs will read respectively'
1639          CALL message( 'rrd_local', 'PA0285', 0, 0, 0, 6, 0 )
1640          IF ( debug_output )  THEN
1641             IF ( files_to_be_opened <= 120 )  THEN
1642                WRITE( debug_string, '(2A,1X,120(I6.6,1X))' )                                         &
1643                     'number of PEs or virtual PE-grid changed in restart run.  PE will read from ',  &
1644                     'files ', file_list(1:files_to_be_opened)
1645             ELSE
1646                WRITE( debug_string, '(3A,1X,120(I6.6,1X),A)' )                                      &
1647                     'number of PEs or virtual PE-grid changed in restart run.  PE will read from ',  &
1648                     'files ', file_list(1:120), '... and more'
1649             ENDIF
1650             CALL debug_message( 'rrd_local', 'info' )
1651          ENDIF
1652       ENDIF
1653
1654!
1655!--    Read data from all restart files determined above
1656       DO  i = 1, files_to_be_opened
1657
1658          j = file_list(i)
1659!
1660!--       Set the filename (underscore followed by four digit processor id)
1661          WRITE (myid_char,'(''_'',I6.6)')  j
1662
1663!
1664!--       Open the restart file. If this file has been created by PE0 (_000000),
1665!--       the global variables at the beginning of the file have to be skipped
1666!--       first.
1667          CALL check_open( 13 )
1668          IF ( j == 0 )  CALL rrd_skip_global
1669
1670!
1671!--       First compare the version numbers
1672          READ ( 13 )  length
1673          READ ( 13 )  restart_string(1:length)
1674          READ ( 13 )  version_on_file
1675
1676          binary_version_local = '5.0'
1677          IF ( TRIM( version_on_file ) /= TRIM( binary_version_local ) )  THEN
1678             WRITE( message_string, * ) 'version mismatch concerning ',                               &
1679                                        'binary_version_local:',                                      &
1680                                        '&version on file    = "', TRIM( version_on_file ), '"',      &
1681                                        '&version in program = "', TRIM( binary_version_local ), '"'
1682             CALL message( 'rrd_local', 'PA0286', 1, 2, 0, 6, 0 )
1683          ENDIF
1684
1685!
1686!--       Read number of processors, processor-id, and array ranges.
1687!--       Compare the array ranges with those stored in the index bound array.
1688          READ ( 13 )  numprocs_on_file, myid_on_file, nxl_on_file, nxr_on_file, nys_on_file,         &
1689                       nyn_on_file, nzb_on_file, nzt_on_file
1690
1691          IF ( nxl_on_file /= hor_index_bounds_previous_run(1,j) )  THEN
1692             WRITE( message_string, * ) 'problem with index bound nxl on ',                           &
1693                                        'restart file "', myid_char, '"',                             &
1694                                        '&nxl = ', nxl_on_file, ' but it should be',                  &
1695                                        '&= ', hor_index_bounds_previous_run(1,j),                    &
1696                                        '&from the index bound information array'
1697             CALL message( 'rrd_local', 'PA0287', 2, 2, -1, 6, 1 )
1698          ENDIF
1699
1700          IF ( nxr_on_file /= hor_index_bounds_previous_run(2,j) )  THEN
1701              WRITE( message_string, * ) 'problem with index bound nxr on ',                          &
1702                                         'restart file "', myid_char, '"'  ,                          &
1703                                         ' nxr = ', nxr_on_file, ' but it should be',                 &
1704                                         ' = ', hor_index_bounds_previous_run(2,j),                   &
1705                                         ' from the index bound information array'
1706             CALL message( 'rrd_local', 'PA0288', 2, 2, -1, 6, 1 )
1707
1708          ENDIF
1709
1710          IF ( nys_on_file /= hor_index_bounds_previous_run(3,j) )  THEN
1711             WRITE( message_string, * ) 'problem with index bound nys on ',                           &
1712                                        'restart file "', myid_char, '"',                             &
1713                                        '&nys = ', nys_on_file, ' but it should be',                  &
1714                                        '&= ', hor_index_bounds_previous_run(3,j),                    &
1715                                        '&from the index bound information array'
1716             CALL message( 'rrd_local', 'PA0289', 2, 2, -1, 6, 1 )
1717          ENDIF
1718
1719          IF ( nyn_on_file /= hor_index_bounds_previous_run(4,j) )  THEN
1720             WRITE( message_string, * ) 'problem with index bound nyn on ',                           &
1721                                        'restart file "', myid_char, '"',                             &
1722                                        '&nyn = ', nyn_on_file, ' but it should be',                  &
1723                                        '&= ', hor_index_bounds_previous_run(4,j),                    &
1724                                        '&from the index bound information array'
1725             CALL message( 'rrd_local', 'PA0290', 2, 2, -1, 6, 1 )
1726          ENDIF
1727
1728          IF ( nzb_on_file /= nzb )  THEN
1729             WRITE( message_string, * ) 'mismatch between actual data and data ',                     &
1730                                        'from prior run on PE ', myid,                                &
1731                                        '&nzb on file = ', nzb_on_file,                               &
1732                                        '&nzb         = ', nzb
1733             CALL message( 'rrd_local', 'PA0291', 1, 2, 0, 6, 0 )
1734          ENDIF
1735
1736          IF ( nzt_on_file /= nzt )  THEN
1737             WRITE( message_string, * ) 'mismatch between actual data and data ',                     &
1738                                        'from prior run on PE ', myid,                                &
1739                                        '&nzt on file = ', nzt_on_file,                               &
1740                                        '&nzt         = ', nzt
1741             CALL message( 'rrd_local', 'PA0292', 1, 2, 0, 6, 0 )
1742          ENDIF
1743
1744!
1745!--       Allocate temporary arrays sized as the arrays on the restart file
1746          ALLOCATE( tmp_2d(nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp),      &
1747                    tmp_3d(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,                               &
1748                           nxl_on_file-nbgp:nxr_on_file+nbgp) )
1749
1750!
1751!--       Read arrays
1752!--       ATTENTION: If the following read commands have been altered, the
1753!--       ---------- version number of the variable binary_version_local must
1754!--                  be altered, too. Furthermore, the output list of arrays in
1755!--                  wrd_write_local must also be altered
1756!--                  accordingly.
1757          READ ( 13 )  length
1758          READ ( 13 )  restart_string(1:length)
1759
1760
1761!
1762!--       Loop over processor specific field data
1763          DO  WHILE ( restart_string(1:length) /= '*** end ***' )
1764
1765!
1766!--          Map data on file as often as needed (data are read only for k=1)
1767             DO  k = 1, overlap_count(i)
1768
1769                found = .FALSE.
1770
1771!
1772!--             Get the index range of the subdomain on file which overlap with
1773!--             the current subdomain
1774                nxlf = nxlfa(i,k)
1775                nxlc = nxlfa(i,k) + offset_xa(i,k)
1776                nxrf = nxrfa(i,k)
1777                nxrc = nxrfa(i,k) + offset_xa(i,k)
1778                nysf = nysfa(i,k)
1779                nysc = nysfa(i,k) + offset_ya(i,k)
1780                nynf = nynfa(i,k)
1781                nync = nynfa(i,k) + offset_ya(i,k)
1782
1783
1784                SELECT CASE ( restart_string(1:length) )
1785
1786                   CASE ( 'ghf_av' )
1787                      IF ( .NOT. ALLOCATED( ghf_av ) )  THEN
1788                         ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
1789                      ENDIF
1790                      IF ( k == 1 )  READ ( 13 )  tmp_2d
1791                      ghf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
1792                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1793
1794                   CASE ( 'e' )
1795                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1796                      e(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                  &
1797                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1798
1799                   CASE ( 'e_av' )
1800                      IF ( .NOT. ALLOCATED( e_av ) )  THEN
1801                         ALLOCATE( e_av(nzb:nzt+1,nys-nbgp:nyn+nbgp,                                  &
1802                                        nxl-nbgp:nxr+nbgp) )
1803                      ENDIF
1804                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1805                      e_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
1806                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1807
1808                   CASE ( 'kh' )
1809                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1810                      kh(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                 &
1811                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1812
1813                   CASE ( 'kh_av' )
1814                      IF ( .NOT. ALLOCATED( kh_av ) )  THEN
1815                         ALLOCATE( kh_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ))
1816                      ENDIF
1817                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1818                      kh_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                              &
1819                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1820
1821                   CASE ( 'km' )
1822                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1823                      km(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                 &
1824                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1825
1826                   CASE ( 'km_av' )
1827                      IF ( .NOT. ALLOCATED( km_av ) )  THEN
1828                         ALLOCATE( km_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ))
1829                      ENDIF
1830                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1831                      km_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                              &
1832                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1833
1834                   CASE ( 'lpt_av' )
1835                      IF ( .NOT. ALLOCATED( lpt_av ) )  THEN
1836                         ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ))
1837                      ENDIF
1838                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1839                      lpt_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                             &
1840                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1841
1842                   CASE ( 'lwp_av' )
1843                      IF ( .NOT. ALLOCATED( lwp_av ) )  THEN
1844                         ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
1845                      ENDIF
1846                      IF ( k == 1 )  READ ( 13 )  tmp_2d
1847                      lwp_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
1848                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1849
1850                   CASE ( 'p' )
1851                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1852                      p(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                  &
1853                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1854
1855                   CASE ( 'p_av' )
1856                      IF ( .NOT. ALLOCATED( p_av ) )  THEN
1857                         ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
1858                      ENDIF
1859                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1860                      p_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
1861                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1862
1863                   CASE ( 'pt' )
1864                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1865                      pt(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                 &
1866                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1867
1868                   CASE ( 'pt_av' )
1869                      IF ( .NOT. ALLOCATED( pt_av ) )  THEN
1870                         ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
1871                      ENDIF
1872                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1873                      pt_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                              &
1874                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1875
1876                   CASE ( 'q' )
1877                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1878                      q(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                  &
1879                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1880
1881                   CASE ( 'q_av' )
1882                      IF ( .NOT. ALLOCATED( q_av ) )  THEN
1883                         ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ))
1884                      ENDIF
1885                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1886                      q_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
1887                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1888
1889                   CASE ( 'ql' )
1890                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1891                      ql(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                 &
1892                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1893
1894                   CASE ( 'ql_av' )
1895                      IF ( .NOT. ALLOCATED( ql_av ) )  THEN
1896                         ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
1897                      ENDIF
1898                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1899                      ql_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                              &
1900                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1901
1902                   CASE ( 'qsurf_av' )
1903                      IF ( .NOT. ALLOCATED( qsurf_av ) )  THEN
1904                         ALLOCATE( qsurf_av(nysg:nyng,nxlg:nxrg) )
1905                      ENDIF
1906                      IF ( k == 1 )  READ ( 13 )  tmp_2d
1907                      qsurf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                             &
1908                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1909
1910                   CASE ( 'qsws_av' )
1911                      IF ( .NOT. ALLOCATED( qsws_av ) )  THEN
1912                         ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
1913                      ENDIF
1914                      IF ( k == 1 )  READ ( 13 )  tmp_2d
1915                      qsws_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  =                             &
1916                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1917
1918                   CASE ( 'qv_av' )
1919                      IF ( .NOT. ALLOCATED( qv_av ) )  THEN
1920                         ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
1921                      ENDIF
1922                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1923                      qv_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                              &
1924                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1925
1926                   CASE ( 'r_a_av' )
1927                      IF ( .NOT. ALLOCATED( r_a_av ) )  THEN
1928                         ALLOCATE( r_a_av(nysg:nyng,nxlg:nxrg) )
1929                      ENDIF
1930                      IF ( k == 1 )  READ ( 13 )  tmp_2d
1931                      r_a_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
1932                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1933
1934                   CASE ( 'random_iv' )  ! still unresolved issue
1935                      IF ( k == 1 )  READ ( 13 )  random_iv
1936                      IF ( k == 1 )  READ ( 13 )  random_iy
1937
1938                   CASE ( 'seq_random_array' )
1939                      ALLOCATE( tmp_2d_id_random(nys_on_file:nyn_on_file,nxl_on_file:nxr_on_file) )
1940                      ALLOCATE( tmp_2d_seq_random(5,nys_on_file:nyn_on_file,nxl_on_file:nxr_on_file) )
1941                      IF ( .NOT. ALLOCATED( id_random_array ) )  THEN
1942                         ALLOCATE( id_random_array(nys:nyn,nxl:nxr) )
1943                      ENDIF
1944                      IF ( .NOT. ALLOCATED( seq_random_array ) )  THEN
1945                         ALLOCATE( seq_random_array(5,nys:nyn,nxl:nxr) )
1946                      ENDIF
1947                      IF ( k == 1 )  READ ( 13 )  tmp_2d_id_random
1948                      IF ( k == 1 )  READ ( 13 )  tmp_2d_seq_random
1949                      id_random_array(nysc:nync,nxlc:nxrc) = tmp_2d_id_random(nysf:nynf,nxlf:nxrf)
1950                      seq_random_array(:,nysc:nync,nxlc:nxrc) = tmp_2d_seq_random(:,nysf:nynf,nxlf:nxrf)
1951                      DEALLOCATE( tmp_2d_id_random, tmp_2d_seq_random )
1952
1953                   CASE ( 's' )
1954                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1955                      s(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                  &
1956                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1957
1958                   CASE ( 's_av' )
1959                      IF ( .NOT. ALLOCATED( s_av ) )  THEN
1960                         ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg))
1961                      ENDIF
1962                      IF ( k == 1 )  READ ( 13 )  tmp_3d
1963                      s_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
1964                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1965
1966                   CASE ( 'shf_av' )
1967                      IF ( .NOT. ALLOCATED( shf_av ) )  THEN
1968                         ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
1969                      ENDIF
1970                      IF ( k == 1 )  READ ( 13 )  tmp_2d
1971                      shf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  =                              &
1972                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1973
1974                   CASE ( 'ssurf_av' )
1975                      IF ( .NOT. ALLOCATED( ssurf_av ) )  THEN
1976                         ALLOCATE( ssurf_av(nysg:nyng,nxlg:nxrg) )
1977                      ENDIF
1978                      IF ( k == 1 )  READ ( 13 )  tmp_2d
1979                      ssurf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                             &
1980                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1981
1982                   CASE ( 'ssws_av' )
1983                      IF ( .NOT. ALLOCATED( ssws_av ) )  THEN
1984                         ALLOCATE( ssws_av(nysg:nyng,nxlg:nxrg) )
1985                      ENDIF
1986                      IF ( k == 1 )  READ ( 13 )  tmp_2d
1987                      ssws_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  =                             &
1988                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1989
1990                   CASE ( 'ts_av' )
1991                      IF ( .NOT. ALLOCATED( ts_av ) )  THEN
1992                         ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
1993                      ENDIF
1994                      IF ( k == 1 )  READ ( 13 )  tmp_2d
1995                      ts_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  =                               &
1996                           tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
1997
1998                   CASE ( 'tsurf_av' )
1999                      IF ( .NOT. ALLOCATED( tsurf_av ) )  THEN
2000                         ALLOCATE( tsurf_av(nysg:nyng,nxlg:nxrg) )
2001                      ENDIF
2002                      IF ( k == 1 )  READ ( 13 )  tmp_2d
2003                      tsurf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                             &
2004                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2005
2006                   CASE ( 'u' )
2007                      IF ( k == 1 )  READ ( 13 )  tmp_3d
2008                      u(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                  &
2009                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2010
2011                   CASE ( 'u_av' )
2012                      IF ( .NOT. ALLOCATED( u_av ) )  THEN
2013                         ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2014                      ENDIF
2015                      IF ( k == 1 )  READ ( 13 )  tmp_3d
2016                      u_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
2017                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2018
2019                   CASE ( 'u_m_l' )
2020                      IF ( k == 1 )  THEN
2021                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,   &
2022                                                       1:2) )
2023                         READ ( 13 )  tmp_3d_non_standard
2024                      ENDIF
2025                      IF ( bc_radiation_l )  THEN
2026                         u_m_l(:,nysc-nbgp:nync+nbgp,:) =  tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
2027                      ENDIF
2028
2029                   CASE ( 'u_m_n' )
2030                      IF ( k == 1 )  THEN
2031                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,ny-1:ny,                             &
2032                                                       nxl_on_file-nbgp:nxr_on_file+nbgp) )
2033                         READ ( 13 )  tmp_3d_non_standard
2034                      ENDIF
2035                      IF ( bc_radiation_n )  THEN
2036                         u_m_n(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
2037                      ENDIF
2038
2039                   CASE ( 'u_m_r' )
2040                      IF ( k == 1 )  THEN
2041                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,   &
2042                                                       nx-1:nx) )
2043                         READ ( 13 )  tmp_3d_non_standard
2044                      ENDIF
2045                      IF ( bc_radiation_r )  THEN
2046                         u_m_r(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
2047                      ENDIF
2048
2049                   CASE ( 'u_m_s' )
2050                      IF ( k == 1 )  THEN
2051                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,0:1,                                 &
2052                                                       nxl_on_file-nbgp:nxr_on_file+nbgp) )
2053                         READ ( 13 )  tmp_3d_non_standard
2054                      ENDIF
2055                      IF ( bc_radiation_s )  THEN
2056                         u_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
2057                      ENDIF
2058
2059                   CASE ( 'us_av' )
2060                      IF ( .NOT. ALLOCATED( us_av ) )  THEN
2061                         ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
2062                      ENDIF
2063                      IF ( k == 1 )  READ ( 13 )  tmp_2d
2064                      us_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  =                               &
2065                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2066
2067                   CASE ( 'v' )
2068                      IF ( k == 1 )  READ ( 13 )  tmp_3d
2069                      v(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                  &
2070                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2071
2072                   CASE ( 'v_av' )
2073                      IF ( .NOT. ALLOCATED( v_av ) )  THEN
2074                         ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2075                      ENDIF
2076                      IF ( k == 1 )  READ ( 13 )  tmp_3d
2077                      v_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
2078                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2079
2080                   CASE ( 'v_m_l' )
2081                      IF ( k == 1 )  THEN
2082                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,   &
2083                                                       0:1) )
2084                         READ ( 13 )  tmp_3d_non_standard
2085                      ENDIF
2086                      IF ( bc_radiation_l )  THEN
2087                         v_m_l(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
2088                      ENDIF
2089
2090                   CASE ( 'v_m_n' )
2091                      IF ( k == 1 )  THEN
2092                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,ny-1:ny,                             &
2093                                                       nxl_on_file-nbgp:nxr_on_file+nbgp) )
2094                         READ ( 13 )  tmp_3d_non_standard
2095                      ENDIF
2096                      IF ( bc_radiation_n )  THEN
2097                         v_m_n(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
2098                      ENDIF
2099
2100                   CASE ( 'v_m_r' )
2101                      IF ( k == 1 )  THEN
2102                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,   &
2103                                                       nx-1:nx) )
2104                         READ ( 13 )  tmp_3d_non_standard
2105                      ENDIF
2106                      IF ( bc_radiation_r )  THEN
2107                         v_m_r(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
2108                      ENDIF
2109
2110                   CASE ( 'v_m_s' )
2111                      IF ( k == 1 )  THEN
2112                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,1:2,                                 &
2113                                                       nxl_on_file-nbgp:nxr_on_file+nbgp) )
2114                         READ ( 13 )  tmp_3d_non_standard
2115                      ENDIF
2116                      IF ( bc_radiation_s )  THEN
2117                         v_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
2118                      ENDIF
2119
2120                   CASE ( 'vpt' )
2121                      IF ( k == 1 )  READ ( 13 )  tmp_3d
2122                      vpt(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                &
2123                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2124
2125                   CASE ( 'vpt_av' )
2126                      IF ( .NOT. ALLOCATED( vpt_av ) )  THEN
2127                         ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2128                      ENDIF
2129                      IF ( k == 1 )  READ ( 13 )  tmp_3d
2130                      vpt_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                             &
2131                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2132
2133                   CASE ( 'w' )
2134                      IF ( k == 1 )  READ ( 13 )  tmp_3d
2135                      w(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                  &
2136                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2137
2138                   CASE ( 'w_av' )
2139                      IF ( .NOT. ALLOCATED( w_av ) )  THEN
2140                         ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2141                      ENDIF
2142                      IF ( k == 1 )  READ ( 13 )  tmp_3d
2143                      w_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                               &
2144                         tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2145
2146                   CASE ( 'w_m_l' )
2147                      IF ( k == 1 )  THEN
2148                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,   &
2149                                                       0:1) )
2150                         READ ( 13 )  tmp_3d_non_standard
2151                      ENDIF
2152                      IF ( bc_radiation_l )  THEN
2153                         w_m_l(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
2154                      ENDIF
2155
2156                   CASE ( 'w_m_n' )
2157                      IF ( k == 1 )  THEN
2158                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,ny-1:ny,                             &
2159                                                       nxl_on_file-nbgp:nxr_on_file+nbgp) )
2160                         READ ( 13 )  tmp_3d_non_standard
2161                      ENDIF
2162                      IF ( bc_radiation_n )  THEN
2163                         w_m_n(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
2164                      ENDIF
2165
2166                   CASE ( 'w_m_r' )
2167                      IF ( k == 1 )  THEN
2168                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,   &
2169                                                       nx-1:nx) )
2170                         READ ( 13 )  tmp_3d_non_standard
2171                      ENDIF
2172                      IF ( bc_radiation_r )  THEN
2173                         w_m_r(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
2174                      ENDIF
2175
2176                   CASE ( 'w_m_s' )
2177                      IF ( k == 1 )  THEN
2178                         ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,0:1,                                 &
2179                                                       nxl_on_file-nbgp:nxr_on_file+nbgp) )
2180                         READ ( 13 )  tmp_3d_non_standard
2181                      ENDIF
2182                      IF ( bc_radiation_s )  THEN
2183                         w_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
2184                      ENDIF
2185
2186                   CASE ( 'z0_av' )
2187                      IF ( .NOT. ALLOCATED( z0_av ) )  THEN
2188                         ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
2189                      ENDIF
2190                      IF ( k == 1 )  READ ( 13 )  tmp_2d
2191                      z0_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  =                               &
2192                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2193
2194                   CASE ( 'z0h_av' )
2195                      IF ( .NOT. ALLOCATED( z0h_av ) )  THEN
2196                         ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
2197                      ENDIF
2198                      IF ( k == 1 )  READ ( 13 )  tmp_2d
2199                      z0h_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  =                              &
2200                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2201
2202                   CASE ( 'z0q_av' )
2203                      IF ( .NOT. ALLOCATED( z0q_av ) )  THEN
2204                         ALLOCATE( z0q_av(nysg:nyng,nxlg:nxrg) )
2205                      ENDIF
2206                      IF ( k == 1 )  READ ( 13 )  tmp_2d
2207                      z0q_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  =                              &
2208                      tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
2209
2210                   CASE DEFAULT
2211
2212!
2213!--                   Read restart data of surfaces
2214                      IF ( .NOT. found )  CALL surface_rrd_local( k, nxlf, nxlc, nxl_on_file, nxrf,   &
2215                                                                  nxr_on_file, nynf, nyn_on_file,     &
2216                                                                  nysf, nysc, nys_on_file, found )
2217!
2218!--                   Read restart data of other modules
2219                      IF ( .NOT. found ) CALL module_interface_rrd_local(                             &
2220                                                                  k, nxlf, nxlc, nxl_on_file, nxrf,   &
2221                                                                  nxrc, nxr_on_file, nynf, nync,      &
2222                                                                  nyn_on_file, nysf, nysc,            &
2223                                                                  nys_on_file, tmp_2d, tmp_3d, found )
2224
2225
2226                      IF ( .NOT. found )  THEN
2227                         WRITE( message_string, * ) 'unknown variable named "',                       &
2228                                                    restart_string(1:length),                         &
2229                                                   '" found in subdomain data ',                      &
2230                                                   'from prior run on PE ', myid
2231                         CALL message( 'rrd_local', 'PA0302', 1, 2, 0, 6, 0 )
2232
2233                      ENDIF
2234
2235                END SELECT
2236
2237             ENDDO ! overlaploop
2238
2239!
2240!--          Deallocate non standard array needed for specific variables only
2241             IF ( ALLOCATED( tmp_3d_non_standard ) )  DEALLOCATE( tmp_3d_non_standard )
2242
2243!
2244!--          Read next character string
2245             READ ( 13 )  length
2246             READ ( 13 )  restart_string(1:length)
2247
2248          ENDDO ! dataloop
2249!
2250!--       Close the restart file
2251          CALL close_file( 13 )
2252
2253          DEALLOCATE( tmp_2d, tmp_3d )
2254
2255       ENDDO  ! loop over restart files
2256!
2257!--    Deallocate temporary buffer arrays
2258       DEALLOCATE( nxlfa )
2259       DEALLOCATE( nxrfa )
2260       DEALLOCATE( nynfa )
2261       DEALLOCATE( nysfa )
2262       DEALLOCATE( offset_xa )
2263       DEALLOCATE( offset_ya )
2264!
2265!--    Restore the original filename for the restart file to be written
2266       myid_char = myid_char_save
2267
2268
2269    ELSEIF ( restart_data_format_input(1:3) == 'mpi' )  THEN
2270
2271!
2272!--    Read global restart data using MPI-IO
2273!
2274!--    Open the MPI-IO restart file.
2275       CALL rd_mpi_io_open( 'read', 'BININ' // TRIM( coupling_char ) )
2276
2277
2278       CALL rd_mpi_io_check_array( 'ghf_av' , found = array_found )
2279       IF ( array_found )  THEN
2280          IF (.NOT. ALLOCATED( ghf_av ) )  ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
2281          CALL rrd_mpi_io( 'ghf_av', ghf_av )
2282       ENDIF
2283
2284       CALL rrd_mpi_io( 'e', e )
2285
2286       CALL rd_mpi_io_check_array( 'e_av' , found = array_found )
2287       IF ( array_found  )  THEN
2288          IF ( .NOT. ALLOCATED( e_av ) )  ALLOCATE( e_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2289          CALL rrd_mpi_io( 'e_av', e_av )
2290       ENDIF
2291
2292       CALL rrd_mpi_io( 'kh', kh )
2293
2294       CALL rd_mpi_io_check_array( 'kh_av' , found = array_found )
2295       IF ( array_found )  THEN
2296          IF ( .NOT. ALLOCATED( kh_av ) )  ALLOCATE( kh_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2297          CALL rrd_mpi_io( 'kh_av', kh_av )
2298       ENDIF
2299
2300       CALL rrd_mpi_io( 'km' , km)
2301
2302       CALL rd_mpi_io_check_array( 'km_av' , found = array_found )
2303       IF ( array_found )  THEN
2304          IF ( .NOT. ALLOCATED( km_av ) )  ALLOCATE( km_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2305          CALL rrd_mpi_io( 'km_av', km_av )
2306       ENDIF
2307
2308       CALL rd_mpi_io_check_array( 'lpt_av' , found = array_found )
2309       IF ( array_found )  THEN
2310          IF ( .NOT. ALLOCATED( lpt_av ) )  ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2311          CALL rrd_mpi_io( 'lpt_av', lpt_av )
2312       ENDIF
2313
2314       CALL rd_mpi_io_check_array( 'lwp_av' , found = array_found )
2315       IF ( array_found )  THEN
2316          IF ( .NOT. ALLOCATED( lwp_av ) )  ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
2317          CALL rrd_mpi_io( 'lwp_av', lwp_av )
2318       ENDIF
2319
2320       CALL rrd_mpi_io( 'p', p)
2321
2322       CALL rd_mpi_io_check_array( 'p_av' , found = array_found )
2323       IF ( array_found )  THEN
2324          IF ( .NOT. ALLOCATED( p_av ) )  ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2325          CALL rrd_mpi_io( 'p_av', p_av )
2326       ENDIF
2327
2328       CALL rrd_mpi_io( 'pt', pt)
2329
2330       CALL rd_mpi_io_check_array( 'pt_av' , found = array_found )
2331       IF ( array_found )  THEN
2332          IF ( .NOT. ALLOCATED( pt_av ) )  ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2333          CALL rrd_mpi_io( 'pt_av', pt_av )
2334       ENDIF
2335
2336       CALL rd_mpi_io_check_array( 'q' , found = array_found )
2337       IF ( array_found )  THEN
2338          CALL rrd_mpi_io( 'q', q )
2339       ENDIF
2340
2341       CALL rd_mpi_io_check_array( 'q_av' , found = array_found )
2342       IF ( array_found )  THEN
2343          IF ( .NOT. ALLOCATED( q_av ) )  ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2344          CALL rrd_mpi_io( 'q_av', q_av )
2345       ENDIF
2346
2347       CALL rd_mpi_io_check_array( 'ql' , found = array_found )
2348       IF ( array_found )  THEN
2349          CALL rrd_mpi_io( 'ql', ql )
2350       ENDIF
2351
2352       CALL rd_mpi_io_check_array( 'ql_av' , found = array_found )
2353       IF ( array_found )  THEN
2354          IF ( .NOT. ALLOCATED( ql_av ) )  ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2355          CALL rrd_mpi_io( 'ql_av', ql_av )
2356       ENDIF
2357
2358       CALL rd_mpi_io_check_array( 'qsurf_av' , found = array_found )
2359       IF ( array_found )  THEN
2360          IF ( .NOT. ALLOCATED( qsurf_av ) )  ALLOCATE( qsurf_av(nysg:nyng,nxlg:nxrg) )
2361          CALL rrd_mpi_io( 'qsurf_av', qsurf_av )
2362       ENDIF
2363
2364       CALL rd_mpi_io_check_array( 'qsws_av' , found = array_found )
2365       IF ( array_found )  THEN
2366          IF ( .NOT. ALLOCATED( qsws_av ) )  ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
2367          CALL rrd_mpi_io( 'qsws_av', qsws_av )
2368       ENDIF
2369
2370       CALL rd_mpi_io_check_array( 'qv_av' , found = array_found )
2371       IF ( array_found )  THEN
2372          IF ( .NOT. ALLOCATED( qv_av ) )  ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2373          CALL rrd_mpi_io( 'qv_av', qv_av )
2374       ENDIF
2375
2376       CALL rd_mpi_io_check_array( 'r_a_av' , found = array_found )
2377       IF ( array_found )  THEN
2378          IF ( .NOT. ALLOCATED( r_a_av ) )  ALLOCATE( r_a_av(nysg:nyng,nxlg:nxrg) )
2379          CALL rrd_mpi_io( 'r_a_av', r_a_av )
2380       ENDIF
2381
2382!
2383!--    ATTENTION: The random seeds are global data! If independent values for every PE are required,
2384!--    the general approach of PE indendent restart will be lost. That means that in general the
2385!--    parallel random number generator in random_generator_parallel_mod should be used!
2386       CALL rrd_mpi_io_global_array( 'random_iv', random_iv )
2387       CALL rrd_mpi_io( 'random_iy', random_iy )
2388
2389       CALL rd_mpi_io_check_array( 'id_random_array' , found = array_found )
2390       IF ( array_found )  THEN
2391          IF ( .NOT. ALLOCATED( id_random_array ) )  ALLOCATE( id_random_array(nys:nyn,nxl:nxr) )
2392          IF ( .NOT. ALLOCATED( seq_random_array ) )  ALLOCATE( seq_random_array(5,nys:nyn,nxl:nxr) )
2393          CALL rrd_mpi_io( 'id_random_array', id_random_array)
2394          DO  i = 1, SIZE( seq_random_array, 1 )
2395             WRITE( tmp_name, '(A,I2.2)' )  'seq_random_array', i
2396             CALL rrd_mpi_io( TRIM(tmp_name), seq_random_array(i,:,:) )
2397          ENDDO
2398       ENDIF
2399
2400       CALL rd_mpi_io_check_array( 's' , found = array_found )
2401       IF ( array_found )  THEN
2402          CALL rrd_mpi_io( 's', s )
2403       ENDIF
2404
2405       CALL rd_mpi_io_check_array( 's_av' , found = array_found )
2406       IF ( array_found )  THEN
2407          IF ( .NOT. ALLOCATED( s_av ) )  ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2408          CALL rrd_mpi_io( 's_av', s_av )
2409       ENDIF
2410
2411       CALL rd_mpi_io_check_array( 'shf_av' , found = array_found )
2412       IF ( array_found )  THEN
2413          IF ( .NOT. ALLOCATED( shf_av ) )  ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
2414          CALL rrd_mpi_io( 'shf_av', shf_av )
2415       ENDIF
2416
2417       CALL rd_mpi_io_check_array( 'ssurf_av' , found = array_found )
2418       IF ( array_found )  THEN
2419          IF ( .NOT. ALLOCATED( ssurf_av ) )  ALLOCATE( ssurf_av(nysg:nyng,nxlg:nxrg) )
2420          CALL rrd_mpi_io( 'ssurf_av', ssurf_av )
2421       ENDIF
2422
2423       CALL rd_mpi_io_check_array( 'ssws_av' , found = array_found )
2424       IF ( array_found )  THEN
2425          IF ( .NOT. ALLOCATED( ssws_av ) )  ALLOCATE( ssws_av(nysg:nyng,nxlg:nxrg) )
2426          CALL rrd_mpi_io( 'ssws_av', ssws_av )
2427       ENDIF
2428
2429       CALL rd_mpi_io_check_array( 'ts_av' , found = array_found )
2430       IF ( array_found )  THEN
2431          IF ( .NOT. ALLOCATED( ts_av ) )  ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
2432          CALL rrd_mpi_io( 'ts_av', ts_av )
2433       ENDIF
2434
2435       CALL rd_mpi_io_check_array( 'tsurf_av' , found = array_found )
2436       IF ( array_found )  THEN
2437          IF ( .NOT. ALLOCATED( tsurf_av ) )  ALLOCATE( tsurf_av(nysg:nyng,nxlg:nxrg) )
2438          CALL rrd_mpi_io( 'tsurf_av', tsurf_av )
2439       ENDIF
2440
2441       CALL rrd_mpi_io( 'u', u)
2442
2443       CALL rd_mpi_io_check_array( 'u_av' , found = array_found )
2444       IF ( array_found )  THEN
2445          IF ( .NOT. ALLOCATED( u_av ) )  ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2446          CALL rrd_mpi_io( 'u_av', u_av )
2447       ENDIF
2448
2449       CALL rd_mpi_io_check_array( 'u_m_l' , found = array_found )
2450       IF ( array_found )  THEN
2451          IF ( .NOT. ALLOCATED( u_m_l ) )  ALLOCATE( u_m_l(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2452          CALL rrd_mpi_io( 'u_m_l', u_m_l )
2453       ENDIF
2454
2455       CALL rd_mpi_io_check_array( 'u_m_n' , found = array_found )
2456       IF ( array_found )  THEN
2457          IF ( .NOT. ALLOCATED( u_m_n ) )  ALLOCATE( u_m_n(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2458          CALL rrd_mpi_io( 'u_m_n', u_m_n )
2459       ENDIF
2460
2461       CALL rd_mpi_io_check_array( 'u_m_r' , found = array_found )
2462       IF ( array_found )  THEN
2463          IF ( .NOT. ALLOCATED( u_m_r ) )  ALLOCATE( u_m_r(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2464          CALL rrd_mpi_io( 'u_m_r', u_m_r )
2465       ENDIF
2466
2467       CALL rd_mpi_io_check_array( 'u_m_s' , found = array_found )
2468       IF ( array_found )  THEN
2469          IF ( .NOT. ALLOCATED( u_m_s ) )  ALLOCATE( u_m_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2470          CALL rrd_mpi_io( 'u_m_s', u_m_s )
2471       ENDIF
2472
2473       CALL rd_mpi_io_check_array( 'us_av' , found = array_found )
2474       IF ( array_found )  THEN
2475          IF ( .NOT. ALLOCATED( us_av ) )  ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
2476          CALL rrd_mpi_io( 'us_av', us_av )
2477       ENDIF
2478
2479       CALL rrd_mpi_io( 'v', v )
2480
2481       CALL rd_mpi_io_check_array( 'v_av' , found = array_found )
2482       IF ( array_found )  THEN
2483          IF ( .NOT. ALLOCATED( v_av ) )  ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2484          CALL rrd_mpi_io( 'v_av', v_av )
2485       ENDIF
2486
2487       CALL rd_mpi_io_check_array( 'v_m_l' , found = array_found )
2488       IF ( array_found )  THEN
2489          IF ( .NOT. ALLOCATED( v_m_l ) )  ALLOCATE( v_m_l(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2490          CALL rrd_mpi_io( 'v_m_l', v_m_l )
2491       ENDIF
2492
2493       CALL rd_mpi_io_check_array( 'v_m_n' , found = array_found )
2494       IF ( array_found )  THEN
2495          IF ( .NOT. ALLOCATED( v_m_n ) )  ALLOCATE( v_m_n(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2496          CALL rrd_mpi_io( 'v_m_n', v_m_n )
2497       ENDIF
2498
2499       CALL rd_mpi_io_check_array( 'v_m_r' , found = array_found )
2500       IF ( array_found )  THEN
2501          IF ( .NOT. ALLOCATED( v_m_r ) )  ALLOCATE( v_m_r(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2502          CALL rrd_mpi_io( 'v_m_r', v_m_r )
2503       ENDIF
2504
2505       CALL rd_mpi_io_check_array( 'v_m_s' , found = array_found )
2506       IF ( array_found )  THEN
2507          IF ( .NOT. ALLOCATED( v_m_s ) )  ALLOCATE( v_m_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2508          CALL rrd_mpi_io( 'v_m_s', v_m_s )
2509       ENDIF
2510
2511       CALL rd_mpi_io_check_array( 'vpt' , found = array_found )
2512       IF ( array_found )  THEN
2513          CALL rrd_mpi_io( 'vpt',  vpt)
2514       ENDIF
2515
2516       CALL rd_mpi_io_check_array( 'vpt_av' , found = array_found )
2517       IF ( array_found )  THEN
2518          IF ( .NOT. ALLOCATED( vpt_av ) )  ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2519          CALL rrd_mpi_io( 'vpt_av', vpt_av )
2520       ENDIF
2521
2522       CALL rrd_mpi_io( 'w', w)
2523
2524       CALL rd_mpi_io_check_array( 'w_av' , found = array_found )
2525       IF ( array_found )  THEN
2526          IF ( .NOT. ALLOCATED( w_av ) )  ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2527          CALL rrd_mpi_io( 'w_av', w_av )
2528       ENDIF
2529
2530       CALL rd_mpi_io_check_array( 'w_m_l' , found = array_found )
2531       IF ( array_found )  THEN
2532          IF ( .NOT. ALLOCATED( w_m_l ) )  ALLOCATE( w_m_l(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2533          CALL rrd_mpi_io( 'w_m_l', w_m_l )
2534       ENDIF
2535
2536       CALL rd_mpi_io_check_array( 'w_m_n' , found = array_found )
2537       IF ( array_found )  THEN
2538          IF ( .NOT. ALLOCATED( w_m_n ) )  ALLOCATE( w_m_n(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2539          CALL rrd_mpi_io( 'w_m_n', w_m_n )
2540       ENDIF
2541
2542       CALL rd_mpi_io_check_array( 'w_m_r' , found = array_found )
2543       IF ( array_found )  THEN
2544          IF ( .NOT. ALLOCATED( w_m_r ) )  ALLOCATE( w_m_r(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2545          CALL rrd_mpi_io( 'w_m_r', w_m_r )
2546       ENDIF
2547
2548       CALL rd_mpi_io_check_array( 'w_m_s' , found = array_found )
2549       IF ( array_found )  THEN
2550          IF ( .NOT. ALLOCATED( w_m_s ) )  ALLOCATE( w_m_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
2551          CALL rrd_mpi_io( 'w_m_s', w_m_s )
2552       ENDIF
2553
2554       CALL rd_mpi_io_check_array( 'z0_av' , found = array_found )
2555       IF ( array_found )  THEN
2556          IF ( .NOT. ALLOCATED( z0_av ) )  ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
2557          CALL rrd_mpi_io( 'z0_av', z0_av )
2558       ENDIF
2559
2560       CALL rd_mpi_io_check_array( 'z0h_av' , found = array_found )
2561       IF ( array_found )  THEN
2562          IF ( .NOT. ALLOCATED( z0h_av ) )  ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
2563          CALL rrd_mpi_io( 'z0h_av', z0h_av )
2564       ENDIF
2565
2566       CALL rd_mpi_io_check_array( 'z0q_av' , found = array_found )
2567       IF ( array_found )  THEN
2568          IF ( .NOT. ALLOCATED( z0q_av ) )  ALLOCATE( z0q_av(nysg:nyng,nxlg:nxrg) )
2569          CALL rrd_mpi_io( 'z0q_av', z0q_av )
2570       ENDIF
2571
2572!
2573!--    Read restart data of surfaces
2574       CALL surface_rrd_local
2575
2576!
2577!--    Read restart data of other modules
2578       CALL module_interface_rrd_local
2579
2580!
2581!--    Close restart file
2582       CALL rd_mpi_io_close
2583
2584    ENDIF
2585
2586    CALL location_message( 'reading local restart data', 'finished' )
2587!
2588!-- End of time measuring for reading binary data
2589    CALL cpu_log( log_point_s(14), 'read-restart-data-local', 'stop' )
2590
2591 END SUBROUTINE rrd_local
2592
2593
2594!------------------------------------------------------------------------------!
2595! Description:
2596! ------------
2597!> Skipping the global control variables from restart-file (binary format)
2598!------------------------------------------------------------------------------!
2599 SUBROUTINE rrd_skip_global
2600
2601
2602    CHARACTER (LEN=1) ::  cdum
2603
2604
2605    READ ( 13 )  length
2606    READ ( 13 )  restart_string(1:length)
2607
2608    DO  WHILE ( restart_string(1:length) /= 'binary_version_local' )
2609
2610       READ ( 13 )  cdum
2611       READ ( 13 )  length
2612       READ ( 13 )  restart_string(1:length)
2613
2614    ENDDO
2615
2616    BACKSPACE ( 13 )
2617    BACKSPACE ( 13 )
2618
2619
2620 END SUBROUTINE rrd_skip_global
2621
2622
2623 END MODULE read_restart_data_mod
Note: See TracBrowser for help on using the repository browser.