Home

Context Navigation

read_restart_data_mod.f90 @ 4112

Last change on this file since 4112 was 4101, checked in by gronemeier, 5 years ago

timestep.f90:

consider 2*Km within diffusion criterion as Km is considered twice within the diffusion of e,
in RANS mode, instead of considering each wind component individually use the wind speed of 3d wind vector in CFL criterion
do not limit the increase of dt based on its previous value in RANS mode

other:

remove dt_old

Property svn:keywords set to Id

File size: 78.7 KB

Rev	Line
[2894]	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	!
[3655]	17	! Copyright 1997-2019 Leibniz Universitaet Hannover
[2894]	18	!------------------------------------------------------------------------------!
	19	!
	20	! Current revisions:
	21	! -----------------
	22	!
[3589]	23	!
[2894]	24	! Former revisions:
	25	! -----------------
	26	! $Id: read_restart_data_mod.f90 4101 2019-07-17 15:14:26Z suehring $
[4101]	27	! remove old_dt
	28	!
	29	! 4039 2019-06-18 10:32:41Z suehring
[4039]	30	! input of uu_av, vv_av, ww_av added
	31	!
	32	! 4017 2019-06-06 12:16:46Z schwenkel
[4009]	33	! bugfix for r3998, allocation of 3d temporary arrays of various dimensions revised
	34	!
	35	! 3998 2019-05-23 13:38:11Z suehring
[3998]	36	! Formatting adjustment
	37	!
	38	! 3994 2019-05-22 18:08:09Z suehring
[3994]	39	! output of turbulence intensity added
	40	!
	41	! 3988 2019-05-22 11:32:37Z kanani
[3988]	42	! + time_virtual_measurement (to enable steering of output interval)
	43	!
	44	! 3936 2019-04-26 15:38:02Z kanani
[3936]	45	! Enable time-averaged output of theta_2m* with restarts
	46	!
	47	! 3767 2019-02-27 08:18:02Z raasch
[3767]	48	! unused variables removed from rrd-subroutines parameter list
	49	!
	50	! 3766 2019-02-26 16:23:41Z raasch
[3766]	51	! first argument removed from module_interface_rrd_*
	52	!
	53	! 3668 2019-01-14 12:49:24Z maronga
[3668]	54	! Removed most_method and increased binary version
	55	!
	56	! 3655 2019-01-07 16:51:22Z knoop
[3637]	57	! Implementation of the PALM module interface
	58	!
	59	! 3589 2018-11-30 15:09:51Z suehring
[3589]	60	! Move the control parameter "salsa" from salsa_mod to control_parameters
	61	! (M. Kurppa)
	62	!
	63	! 3582 2018-11-29 19:16:36Z suehring
[3298]	64	! Added chemistry profiles for restart run (basit)
	65	!
	66	! 3294 2018-10-01 02:37:10Z raasch
[3294]	67	! changes concerning modularization of ocean option
	68	!
	69	! 3289 2018-09-28 10:23:58Z suehring
[3289]	70	! Introduce module parameter for number of inflow profiles
	71	!
	72	! 3288 2018-09-28 10:23:08Z suehring
[3274]	73	! Modularization of all bulk cloud physics code components
	74	!
	75	! 3241 2018-09-12 15:02:00Z raasch
[3241]	76	! unused variables removed
	77	!
	78	! 3183 2018-07-27 14:25:55Z suehring
[3183]	79	! Rename variables for boundary flags
	80	!
	81	! 3182 2018-07-27 13:36:03Z suehring
[3065]	82	! New parameters concerning vertical grid stretching have been added
	83	!
	84	! 3056 2018-06-04 07:49:35Z Giersch
[3056]	85	! found variable has to be set to false inside overlap loop
	86	!
	87	! 3049 2018-05-29 13:52:36Z Giersch
[3045]	88	! Error messages revised
	89	!
[3049]	90	! 3045 2018-05-28 07:55:41Z Giersch
	91	! Error messages revised
	92	!
[3045]	93	! 3004 2018-04-27 12:33:25Z Giersch
[3004]	94	! precipitation_rate_av removed
	95	!
	96	! 3003 2018-04-23 10:22:58Z Giersch
[3003]	97	! z_i is also read to use the last known inversion height from the
	98	! initial run as the first inversion height which is written into the
	99	! run control file
	100	!
	101	! 2956 2018-04-10 11:01:03Z Giersch
[2956]	102	! spectrum_x and spectrum_y have been moved to global data
	103	!
	104	! 2921 2018-03-22 15:05:23Z Giersch
[2921]	105	! spinup_time, day_of_year_init and time_utc_init are also read now
	106	!
	107	! 2912 2018-03-20 13:00:05Z knoop
[2912]	108	! Added gust module interface calls
	109	!
	110	! 2894 2018-03-15 09:17:58Z Giersch
[2894]	111	! Initial revision
	112	!
	113	!
	114	! Description:
	115	! ------------
	116	!> Reads restart data from restart-file(s) (binary format).
[3298]	117	!>
	118	!> @todo: Revise max_pr_cs (profiles for chemistry)
[4039]	119	!> @todo: Modularize reading of restart data for diagnostic quantities, which
	120	!> is not possible with the current module-interface structure
[2894]	121	!------------------------------------------------------------------------------!
	122	MODULE read_restart_data_mod
	123
	124
[3637]	125	USE arrays_3d, &
	126	ONLY: inflow_damping_factor, mean_inflow_profiles, pt_init, &
	127	q_init, ref_state, sa_init, s_init, u_init, ug, v_init, vg, &
	128	e, kh, km, p, pt, q, ql, s, u, u_m_l, u_m_n, u_m_r, u_m_s, &
	129	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
[2894]	130
[3637]	131	USE averaging
	132
[3298]	133	USE chem_modules, &
	134	ONLY: max_pr_cs
	135
[3637]	136	USE control_parameters
[3298]	137
[3637]	138	USE cpulog, &
	139	ONLY: cpu_log, log_point_s
	140
	141	USE date_and_time_mod, &
	142	ONLY: day_of_year_init, time_utc_init
	143
[3994]	144	USE diagnostic_output_quantities_mod, &
[4039]	145	ONLY: ti_av, uu_av, vv_av, ww_av
[3994]	146
[3637]	147	USE grid_variables, &
	148	ONLY: dx, dy
	149
	150	USE indices, &
	151	ONLY: nbgp, nx, nxl, nxlg, nxr, nxrg, nx_on_file, ny, nys, nysg, nyn, &
	152	nyng, ny_on_file, nz, nzb, nzt
	153
	154	USE kinds
	155
	156	USE model_1d_mod, &
	157	ONLY: damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
	158
	159	USE module_interface, &
	160	ONLY: module_interface_rrd_global, &
	161	module_interface_rrd_local
	162
	163	USE netcdf_interface, &
	164	ONLY: netcdf_precision, output_for_t0
	165
	166	USE pegrid
	167
	168	USE radiation_model_mod, &
	169	ONLY: time_radiation
	170
	171	USE random_function_mod, &
	172	ONLY: random_iv, random_iy
	173
	174	USE random_generator_parallel, &
	175	ONLY: id_random_array, seq_random_array
	176
	177	USE spectra_mod, &
	178	ONLY: average_count_sp, spectrum_x, spectrum_y
	179
	180	USE surface_mod, &
	181	ONLY : surface_rrd_local
	182
	183	USE statistics, &
	184	ONLY: statistic_regions, hom, hom_sum, pr_palm, u_max, u_max_ijk, &
	185	v_max, v_max_ijk, w_max, w_max_ijk, z_i
	186
	187	USE vertical_nesting_mod, &
	188	ONLY: vnest_init
	189
[3988]	190	USE virtual_measurement_mod, &
	191	ONLY: time_virtual_measurement
[3637]	192
[3988]	193
[2894]	194	IMPLICIT NONE
	195
	196
	197	INTERFACE rrd_global
	198	MODULE PROCEDURE rrd_global
	199	END INTERFACE rrd_global
	200
	201	INTERFACE rrd_read_parts_of_global
	202	MODULE PROCEDURE rrd_read_parts_of_global
	203	END INTERFACE rrd_read_parts_of_global
	204
	205	INTERFACE rrd_local
	206	MODULE PROCEDURE rrd_local
	207	END INTERFACE rrd_local
	208
	209	INTERFACE rrd_skip_global
	210	MODULE PROCEDURE rrd_skip_global
	211	END INTERFACE rrd_skip_global
	212
	213
[3294]	214	PUBLIC rrd_global, rrd_read_parts_of_global, rrd_local, rrd_skip_global
[2894]	215
	216
	217	CONTAINS
	218
[3294]	219	!------------------------------------------------------------------------------!
[2894]	220	! Description:
	221	! ------------
	222	!> Reads values of global control variables from restart-file (binary format)
	223	!> created by PE0 of the previous run
	224	!------------------------------------------------------------------------------!
	225	SUBROUTINE rrd_global
	226
	227
	228	CHARACTER (LEN=10) :: binary_version_global, version_on_file
	229
	230	LOGICAL :: found
	231
	232
	233	CALL check_open( 13 )
	234	!
[3294]	235	!-- Make version number check first
[2894]	236	READ ( 13 ) length
	237	READ ( 13 ) restart_string(1:length)
	238	READ ( 13 ) version_on_file
	239
[3668]	240	binary_version_global = '4.8'
[2894]	241	IF ( TRIM( version_on_file ) /= TRIM( binary_version_global ) ) THEN
[3045]	242	WRITE( message_string, * ) 'version mismatch concerning ', &
	243	'binary_version_global:', &
[3046]	244	'&version on file = "', &
[2894]	245	TRIM( version_on_file ), '"', &
[3046]	246	'&version on program = "', &
[2894]	247	TRIM( binary_version_global ), '"'
	248	CALL message( 'rrd_global', 'PA0296', 1, 2, 0, 6, 0 )
	249	ENDIF
	250
	251	!
[3294]	252	!-- Read number of PEs and horizontal index bounds of all PEs used in the
	253	!-- previous run
[2894]	254	READ ( 13 ) length
	255	READ ( 13 ) restart_string(1:length)
	256
	257	IF ( TRIM( restart_string(1:length) ) /= 'numprocs' ) THEN
	258	WRITE( message_string, * ) 'numprocs not found in data from prior ', &
	259	'run on PE ', myid
	260	CALL message( 'rrd_global', 'PA0297', 1, 2, 0, 6, 0 )
	261	ENDIF
	262	READ ( 13 ) numprocs_previous_run
	263
	264	IF ( .NOT. ALLOCATED( hor_index_bounds_previous_run ) ) THEN
	265	ALLOCATE( hor_index_bounds_previous_run(4,0:numprocs_previous_run-1) )
	266	ENDIF
	267
	268	READ ( 13 ) length
	269	READ ( 13 ) restart_string(1:length)
	270
	271	IF ( restart_string(1:length) /= 'hor_index_bounds' ) THEN
	272	WRITE( message_string, * ) 'hor_index_bounds not found in data ', &
	273	'from prior run on PE ', myid
	274	CALL message( 'rrd_global', 'PA0298', 1, 2, 0, 6, 0 )
	275	ENDIF
	276	READ ( 13 ) hor_index_bounds_previous_run
	277
	278	!
[3294]	279	!-- Read vertical number of gridpoints and number of different areas used
	280	!-- for computing statistics. Allocate arrays depending on these values,
	281	!-- which are needed for the following read instructions.
[2894]	282	READ ( 13 ) length
	283	READ ( 13 ) restart_string(1:length)
	284
	285	IF ( restart_string(1:length) /= 'nz' ) THEN
	286	WRITE( message_string, * ) 'nz not found in data from prior run ', &
	287	'on PE ', myid
	288	CALL message( 'rrd_global', 'PA0299', 1, 2, 0, 6, 0 )
	289	ENDIF
	290	READ ( 13 ) nz
	291
	292	READ ( 13 ) length
	293	READ ( 13 ) restart_string(1:length)
	294
	295	IF ( restart_string(1:length) /= 'max_pr_user' ) THEN
	296	WRITE( message_string, * ) 'max_pr_user not found in data from ', &
	297	'prior run on PE ', myid
	298	CALL message( 'rrd_global', 'PA0300', 1, 2, 0, 6, 0 )
	299	ENDIF
	300	READ ( 13 ) max_pr_user ! This value is checked against the number of
	301	! user profiles given for the current run
	302	! in routine user_parin (it has to match)
	303
	304	READ ( 13 ) length
	305	READ ( 13 ) restart_string(1:length)
	306
	307	IF ( restart_string(1:length) /= 'statistic_regions' ) THEN
	308	WRITE( message_string, * ) 'statistic_regions not found in data ', &
	309	'from prior run on PE ', myid
	310	CALL message( 'rrd_global', 'PA0301', 1, 2, 0, 6, 0 )
	311	ENDIF
	312	READ ( 13 ) statistic_regions
[3294]	313
[2894]	314	IF ( .NOT. ALLOCATED( ug ) ) THEN
	315	ALLOCATE( ug(0:nz+1), u_init(0:nz+1), vg(0:nz+1), &
	316	v_init(0:nz+1), pt_init(0:nz+1), q_init(0:nz+1), &
	317	ref_state(0:nz+1), s_init(0:nz+1), sa_init(0:nz+1), &
[3298]	318	hom(0:nz+1,2,pr_palm+max_pr_user+max_pr_cs,0:statistic_regions), &
	319	hom_sum(0:nz+1,pr_palm+max_pr_user+max_pr_cs,0:statistic_regions) )
[2894]	320	ENDIF
	321
	322	!
[3294]	323	!-- Now read all control parameters:
	324	!-- Caution: When the following read instructions have been changed, the
	325	!-- ------- version number stored in the variable binary_version_global has
	326	!-- to be increased. The same changes must also be done in
	327	!-- wrd_write_global.
[2894]	328	READ ( 13 ) length
	329	READ ( 13 ) restart_string(1:length)
	330
	331	DO WHILE ( restart_string(1:length) /= 'binary_version_local' )
	332
	333	found = .FALSE.
	334
	335	SELECT CASE ( restart_string(1:length) )
	336
	337	CASE ( 'advected_distance_x' )
	338	READ ( 13 ) advected_distance_x
	339	CASE ( 'advected_distance_y' )
	340	READ ( 13 ) advected_distance_y
	341	CASE ( 'alpha_surface' )
	342	READ ( 13 ) alpha_surface
	343	CASE ( 'average_count_pr' )
	344	READ ( 13 ) average_count_pr
	345	CASE ( 'average_count_sp' )
	346	READ ( 13 ) average_count_sp
	347	CASE ( 'average_count_3d' )
	348	READ ( 13 ) average_count_3d
	349	CASE ( 'bc_e_b' )
	350	READ ( 13 ) bc_e_b
	351	CASE ( 'bc_lr' )
	352	READ ( 13 ) bc_lr
	353	CASE ( 'bc_ns' )
	354	READ ( 13 ) bc_ns
	355	CASE ( 'bc_p_b' )
	356	READ ( 13 ) bc_p_b
	357	CASE ( 'bc_p_t' )
	358	READ ( 13 ) bc_p_t
	359	CASE ( 'bc_pt_b' )
	360	READ ( 13 ) bc_pt_b
	361	CASE ( 'bc_pt_t' )
	362	READ ( 13 ) bc_pt_t
	363	CASE ( 'bc_pt_t_val' )
	364	READ ( 13 ) bc_pt_t_val
	365	CASE ( 'bc_q_b' )
	366	READ ( 13 ) bc_q_b
	367	CASE ( 'bc_q_t' )
	368	READ ( 13 ) bc_q_t
	369	CASE ( 'bc_q_t_val' )
	370	READ ( 13 ) bc_q_t_val
	371	CASE ( 'bc_s_b' )
	372	READ ( 13 ) bc_s_b
	373	CASE ( 'bc_s_t' )
	374	READ ( 13 ) bc_s_t
	375	CASE ( 'bc_uv_b' )
	376	READ ( 13 ) bc_uv_b
	377	CASE ( 'bc_uv_t' )
	378	READ ( 13 ) bc_uv_t
	379	CASE ( 'building_height' )
	380	READ ( 13 ) building_height
	381	CASE ( 'building_length_x' )
	382	READ ( 13 ) building_length_x
	383	CASE ( 'building_length_y' )
	384	READ ( 13 ) building_length_y
	385	CASE ( 'building_wall_left' )
	386	READ ( 13 ) building_wall_left
	387	CASE ( 'building_wall_south' )
	388	READ ( 13 ) building_wall_south
	389	CASE ( 'call_psolver_at_all_substeps' )
	390	READ ( 13 ) call_psolver_at_all_substeps
	391	CASE ( 'canyon_height' )
	392	READ ( 13 ) canyon_height
	393	CASE ( 'canyon_wall_left' )
	394	READ ( 13 ) canyon_wall_left
	395	CASE ( 'canyon_wall_south' )
	396	READ ( 13 ) canyon_wall_south
	397	CASE ( 'canyon_width_x' )
	398	READ ( 13 ) canyon_width_x
	399	CASE ( 'canyon_width_y' )
	400	READ ( 13 ) canyon_width_y
	401	CASE ( 'cfl_factor' )
	402	READ ( 13 ) cfl_factor
	403	CASE ( 'cloud_droplets' )
	404	READ ( 13 ) cloud_droplets
	405	CASE ( 'collective_wait' )
	406	READ ( 13 ) collective_wait
	407	CASE ( 'conserve_volume_flow' )
	408	READ ( 13 ) conserve_volume_flow
	409	CASE ( 'conserve_volume_flow_mode' )
	410	READ ( 13 ) conserve_volume_flow_mode
	411	CASE ( 'constant_flux_layer' )
	412	READ ( 13 ) constant_flux_layer
	413	CASE ( 'coupling_start_time' )
	414	READ ( 13 ) coupling_start_time
	415	CASE ( 'current_timestep_number' )
	416	READ ( 13 ) current_timestep_number
	417	CASE ( 'cycle_mg' )
	418	READ ( 13 ) cycle_mg
	419	CASE ( 'damp_level_1d' )
	420	READ ( 13 ) damp_level_1d
[2921]	421	CASE ( 'day_of_year_init' )
	422	READ ( 13 ) day_of_year_init
[2894]	423	CASE ( 'dissipation_1d' )
	424	READ ( 13 ) dissipation_1d
	425	CASE ( 'do2d_xy_time_count' )
	426	READ ( 13 ) do2d_xy_time_count
	427	CASE ( 'do2d_xz_time_count' )
	428	READ ( 13 ) do2d_xz_time_count
	429	CASE ( 'do2d_yz_time_count' )
	430	READ ( 13 ) do2d_yz_time_count
	431	CASE ( 'do3d_time_count' )
	432	READ ( 13 ) do3d_time_count
	433	CASE ( 'dp_external' )
	434	READ ( 13 ) dp_external
	435	CASE ( 'dp_level_b' )
	436	READ ( 13 ) dp_level_b
	437	CASE ( 'dp_smooth' )
	438	READ ( 13 ) dp_smooth
	439	CASE ( 'dpdxy' )
	440	READ ( 13 ) dpdxy
	441	CASE ( 'dt_3d' )
	442	READ ( 13 ) dt_3d
	443	CASE ( 'dt_pr_1d' )
	444	READ ( 13 ) dt_pr_1d
	445	CASE ( 'dt_run_control_1d' )
	446	READ ( 13 ) dt_run_control_1d
	447	CASE ( 'dx' )
	448	READ ( 13 ) dx
	449	CASE ( 'dy' )
	450	READ ( 13 ) dy
	451	CASE ( 'dz' )
	452	READ ( 13 ) dz
	453	CASE ( 'dz_max' )
	454	READ ( 13 ) dz_max
	455	CASE ( 'dz_stretch_factor' )
	456	READ ( 13 ) dz_stretch_factor
[3065]	457	CASE ( 'dz_stretch_factor_array' )
	458	READ ( 13 ) dz_stretch_factor_array
[2894]	459	CASE ( 'dz_stretch_level' )
	460	READ ( 13 ) dz_stretch_level
[3065]	461	CASE ( 'dz_stretch_level_end' )
	462	READ ( 13 ) dz_stretch_level_end
	463	CASE ( 'dz_stretch_level_start' )
	464	READ ( 13 ) dz_stretch_level_start
[2894]	465	CASE ( 'e_min' )
	466	READ ( 13 ) e_min
	467	CASE ( 'end_time_1d' )
	468	READ ( 13 ) end_time_1d
	469	CASE ( 'fft_method' )
	470	READ ( 13 ) fft_method
	471	CASE ( 'first_call_lpm' )
	472	READ ( 13 ) first_call_lpm
	473	CASE ( 'galilei_transformation' )
	474	READ ( 13 ) galilei_transformation
	475	CASE ( 'hom' )
	476	READ ( 13 ) hom
	477	CASE ( 'hom_sum' )
	478	READ ( 13 ) hom_sum
	479	CASE ( 'humidity' )
	480	READ ( 13 ) humidity
	481	CASE ( 'inflow_damping_factor' )
	482	IF ( .NOT. ALLOCATED( inflow_damping_factor ) ) THEN
	483	ALLOCATE( inflow_damping_factor(0:nz+1) )
	484	ENDIF
	485	READ ( 13 ) inflow_damping_factor
	486	CASE ( 'inflow_damping_height' )
	487	READ ( 13 ) inflow_damping_height
	488	CASE ( 'inflow_damping_width' )
	489	READ ( 13 ) inflow_damping_width
	490	CASE ( 'inflow_disturbance_begin' )
	491	READ ( 13 ) inflow_disturbance_begin
	492	CASE ( 'inflow_disturbance_end' )
	493	READ ( 13 ) inflow_disturbance_end
	494	CASE ( 'km_constant' )
	495	READ ( 13 ) km_constant
	496	CASE ( 'large_scale_forcing' )
	497	READ ( 13 ) large_scale_forcing
	498	CASE ( 'large_scale_subsidence' )
	499	READ ( 13 ) large_scale_subsidence
	500	CASE ( 'latitude' )
	501	READ ( 13 ) latitude
	502	CASE ( 'longitude' )
	503	READ ( 13 ) longitude
	504	CASE ( 'loop_optimization' )
	505	READ ( 13 ) loop_optimization
	506	CASE ( 'masking_method' )
	507	READ ( 13 ) masking_method
	508	CASE ( 'mean_inflow_profiles' )
	509	IF ( .NOT. ALLOCATED( mean_inflow_profiles ) ) THEN
[3288]	510	ALLOCATE( mean_inflow_profiles(0:nz+1,1:num_mean_inflow_profiles) )
[2894]	511	ENDIF
	512	READ ( 13 ) mean_inflow_profiles
	513	CASE ( 'mg_cycles' )
	514	READ ( 13 ) mg_cycles
	515	CASE ( 'mg_switch_to_pe0_level' )
	516	READ ( 13 ) mg_switch_to_pe0_level
	517	CASE ( 'mixing_length_1d' )
	518	READ ( 13 ) mixing_length_1d
	519	CASE ( 'momentum_advec' )
	520	READ ( 13 ) momentum_advec
	521	CASE ( 'netcdf_precision' )
	522	READ ( 13 ) netcdf_precision
	523	CASE ( 'neutral' )
	524	READ ( 13 ) neutral
	525	CASE ( 'ngsrb' )
	526	READ ( 13 ) ngsrb
	527	CASE ( 'nsor' )
	528	READ ( 13 ) nsor
	529	CASE ( 'nsor_ini' )
	530	READ ( 13 ) nsor_ini
	531	CASE ( 'nudging' )
	532	READ ( 13 ) nudging
	533	CASE ( 'num_leg' )
	534	READ ( 13 ) num_leg
	535	CASE ( 'nx' )
	536	READ ( 13 ) nx
	537	nx_on_file = nx
	538	CASE ( 'ny' )
	539	READ ( 13 ) ny
	540	ny_on_file = ny
[3294]	541	CASE ( 'ocean_mode' )
	542	READ ( 13 ) ocean_mode
[2894]	543	CASE ( 'omega' )
	544	READ ( 13 ) omega
	545	CASE ( 'omega_sor' )
	546	READ ( 13 ) omega_sor
	547	CASE ( 'output_for_t0' )
	548	READ (13) output_for_t0
	549	CASE ( 'passive_scalar' )
	550	READ ( 13 ) passive_scalar
	551	CASE ( 'prandtl_number' )
	552	READ ( 13 ) prandtl_number
	553	CASE ( 'psolver' )
	554	READ ( 13 ) psolver
	555	CASE ( 'pt_damping_factor' )
	556	READ ( 13 ) pt_damping_factor
	557	CASE ( 'pt_damping_width' )
	558	READ ( 13 ) pt_damping_width
	559	CASE ( 'pt_init' )
	560	READ ( 13 ) pt_init
	561	CASE ( 'pt_reference' )
	562	READ ( 13 ) pt_reference
	563	CASE ( 'pt_surface' )
	564	READ ( 13 ) pt_surface
	565	CASE ( 'pt_surface_initial_change' )
	566	READ ( 13 ) pt_surface_initial_change
	567	CASE ( 'pt_vertical_gradient' )
	568	READ ( 13 ) pt_vertical_gradient
	569	CASE ( 'pt_vertical_gradient_level' )
	570	READ ( 13 ) pt_vertical_gradient_level
	571	CASE ( 'pt_vertical_gradient_level_ind' )
	572	READ ( 13 ) pt_vertical_gradient_level_ind
	573	CASE ( 'q_init' )
	574	READ ( 13 ) q_init
	575	CASE ( 'q_surface' )
	576	READ ( 13 ) q_surface
	577	CASE ( 'q_surface_initial_change' )
	578	READ ( 13 ) q_surface_initial_change
	579	CASE ( 'q_vertical_gradient' )
	580	READ ( 13 ) q_vertical_gradient
	581	CASE ( 'q_vertical_gradient_level' )
	582	READ ( 13 ) q_vertical_gradient_level
	583	CASE ( 'q_vertical_gradient_level_ind' )
	584	READ ( 13 ) q_vertical_gradient_level_ind
	585	CASE ( 'random_generator' )
	586	READ ( 13 ) random_generator
	587	CASE ( 'random_heatflux' )
	588	READ ( 13 ) random_heatflux
	589	CASE ( 'rans_mode' )
	590	READ ( 13 ) rans_mode
	591	CASE ( 'rayleigh_damping_factor' )
	592	READ ( 13 ) rayleigh_damping_factor
	593	CASE ( 'rayleigh_damping_height' )
	594	READ ( 13 ) rayleigh_damping_height
	595	CASE ( 'recycling_width' )
	596	READ ( 13 ) recycling_width
	597	CASE ( 'recycling_yshift' )
	598	READ ( 13 ) recycling_yshift
	599	CASE ( 'ref_state' )
	600	READ ( 13 ) ref_state
	601	CASE ( 'reference_state' )
	602	READ ( 13 ) reference_state
	603	CASE ( 'residual_limit' )
	604	READ ( 13 ) residual_limit
	605	CASE ( 'roughness_length' )
	606	READ ( 13 ) roughness_length
	607	CASE ( 'run_coupled' )
	608	READ ( 13 ) run_coupled
	609	CASE ( 'runnr' )
	610	READ ( 13 ) runnr
	611	CASE ( 's_init' )
	612	READ ( 13 ) s_init
	613	CASE ( 's_surface' )
	614	READ ( 13 ) s_surface
	615	CASE ( 's_surface_initial_change' )
	616	READ ( 13 ) s_surface_initial_change
	617	CASE ( 's_vertical_gradient' )
	618	READ ( 13 ) s_vertical_gradient
	619	CASE ( 's_vertical_gradient_level' )
	620	READ ( 13 ) s_vertical_gradient_level
	621	CASE ( 's_vertical_gradient_level_ind' )
	622	READ ( 13 ) s_vertical_gradient_level_ind
	623	CASE ( 'scalar_advec' )
	624	READ ( 13 ) scalar_advec
	625	CASE ( 'simulated_time' )
	626	READ ( 13 ) simulated_time
[2956]	627	CASE ( 'spectrum_x' )
	628	IF ( .NOT. ALLOCATED( spectrum_x ) ) THEN
	629	ALLOCATE( spectrum_x( 1:nx/2, 1:100, 1:10 ) )
	630	ENDIF
	631	READ ( 13 ) spectrum_x
	632	CASE ( 'spectrum_y' )
	633	IF ( .NOT. ALLOCATED( spectrum_y ) ) THEN
	634	ALLOCATE( spectrum_y( 1:ny/2, 1:100, 1:10 ) )
	635	ENDIF
	636	READ ( 13 ) spectrum_y
[2921]	637	CASE ( 'spinup_time' )
	638	READ ( 13 ) spinup_time
[2894]	639	CASE ( 'surface_heatflux' )
	640	READ ( 13 ) surface_heatflux
	641	CASE ( 'surface_pressure' )
	642	READ ( 13 ) surface_pressure
	643	CASE ( 'surface_scalarflux' )
	644	READ ( 13 ) surface_scalarflux
	645	CASE ( 'surface_waterflux' )
	646	READ ( 13 ) surface_waterflux
	647	CASE ( 'time_coupling' )
	648	READ ( 13 ) time_coupling
	649	CASE ( 'time_disturb' )
	650	READ ( 13 ) time_disturb
	651	CASE ( 'time_do2d_xy' )
	652	READ ( 13 ) time_do2d_xy
	653	CASE ( 'time_do2d_xz' )
	654	READ ( 13 ) time_do2d_xz
	655	CASE ( 'time_do2d_yz' )
	656	READ ( 13 ) time_do2d_yz
	657	CASE ( 'time_do3d' )
	658	READ ( 13 ) time_do3d
	659	CASE ( 'time_do_av' )
	660	READ ( 13 ) time_do_av
	661	CASE ( 'time_do_sla' )
	662	READ ( 13 ) time_do_sla
	663	CASE ( 'time_domask' )
	664	READ ( 13 ) time_domask
	665	CASE ( 'time_dopr' )
	666	READ ( 13 ) time_dopr
	667	CASE ( 'time_dopr_av' )
	668	READ ( 13 ) time_dopr_av
	669	CASE ( 'time_dopr_listing' )
	670	READ ( 13 ) time_dopr_listing
	671	CASE ( 'time_dopts' )
	672	READ ( 13 ) time_dopts
	673	CASE ( 'time_dosp' )
	674	READ ( 13 ) time_dosp
	675	CASE ( 'time_dots' )
	676	READ ( 13 ) time_dots
	677	CASE ( 'time_radiation' )
	678	READ ( 13 ) time_radiation
	679	CASE ( 'time_restart' )
	680	READ ( 13 ) time_restart
	681	CASE ( 'time_run_control' )
	682	READ ( 13 ) time_run_control
	683	CASE ( 'time_since_reference_point' )
	684	READ ( 13 ) time_since_reference_point
[2921]	685	CASE ( 'time_utc_init' )
	686	READ ( 13 ) time_utc_init
[3988]	687	CASE ( 'time_virtual_measurement' )
	688	READ ( 13 ) time_virtual_measurement
[2894]	689	CASE ( 'timestep_scheme' )
	690	READ ( 13 ) timestep_scheme
	691	CASE ( 'top_heatflux' )
	692	READ ( 13 ) top_heatflux
	693	CASE ( 'top_momentumflux_u' )
	694	READ ( 13 ) top_momentumflux_u
	695	CASE ( 'top_momentumflux_v' )
	696	READ ( 13 ) top_momentumflux_v
	697	CASE ( 'top_scalarflux' )
	698	READ ( 13 ) top_scalarflux
	699	CASE ( 'topography' )
	700	READ ( 13 ) topography
	701	CASE ( 'topography_grid_convention' )
	702	READ ( 13 ) topography_grid_convention
	703	CASE ( 'tsc' )
	704	READ ( 13 ) tsc
	705	CASE ( 'tunnel_height' )
	706	READ ( 13 ) tunnel_height
	707	CASE ( 'tunnel_length' )
	708	READ ( 13 ) tunnel_length
	709	CASE ( 'tunnel_wall_depth' )
	710	READ ( 13 ) tunnel_wall_depth
	711	CASE ( 'tunnel_width_x' )
	712	READ ( 13 ) tunnel_width_x
	713	CASE ( 'tunnel_width_y' )
	714	READ ( 13 ) tunnel_width_y
	715	CASE ( 'turbulence_closure' )
	716	READ ( 13 ) turbulence_closure
	717	CASE ( 'turbulent_inflow' )
	718	READ ( 13 ) turbulent_inflow
	719	CASE ( 'u_bulk' )
	720	READ ( 13 ) u_bulk
	721	CASE ( 'u_init' )
	722	READ ( 13 ) u_init
	723	CASE ( 'u_max' )
	724	READ ( 13 ) u_max
	725	CASE ( 'u_max_ijk' )
	726	READ ( 13 ) u_max_ijk
	727	CASE ( 'ug' )
	728	READ ( 13 ) ug
	729	CASE ( 'ug_surface' )
	730	READ ( 13 ) ug_surface
	731	CASE ( 'ug_vertical_gradient' )
	732	READ ( 13 ) ug_vertical_gradient
	733	CASE ( 'ug_vertical_gradient_level' )
	734	READ ( 13 ) ug_vertical_gradient_level
	735	CASE ( 'ug_vertical_gradient_level_ind' )
	736	READ ( 13 ) ug_vertical_gradient_level_ind
	737	CASE ( 'use_surface_fluxes' )
	738	READ ( 13 ) use_surface_fluxes
	739	CASE ( 'use_top_fluxes' )
	740	READ ( 13 ) use_top_fluxes
	741	CASE ( 'use_ug_for_galilei_tr' )
	742	READ ( 13 ) use_ug_for_galilei_tr
	743	CASE ( 'use_upstream_for_tke' )
	744	READ ( 13 ) use_upstream_for_tke
	745	CASE ( 'v_bulk' )
	746	READ ( 13 ) v_bulk
	747	CASE ( 'v_init' )
	748	READ ( 13 ) v_init
	749	CASE ( 'v_max' )
	750	READ ( 13 ) v_max
	751	CASE ( 'v_max_ijk' )
	752	READ ( 13 ) v_max_ijk
	753	CASE ( 'vg' )
	754	READ ( 13 ) vg
	755	CASE ( 'vg_surface' )
	756	READ ( 13 ) vg_surface
	757	CASE ( 'vg_vertical_gradient' )
	758	READ ( 13 ) vg_vertical_gradient
	759	CASE ( 'vg_vertical_gradient_level' )
	760	READ ( 13 ) vg_vertical_gradient_level
	761	CASE ( 'vg_vertical_gradient_level_ind' )
	762	READ ( 13 ) vg_vertical_gradient_level_ind
	763	CASE ( 'virtual_flight' )
	764	READ ( 13 ) virtual_flight
	765	CASE ( 'vnest_init' )
	766	READ ( 13 ) vnest_init
	767	CASE ( 'volume_flow_area' )
	768	READ ( 13 ) volume_flow_area
	769	CASE ( 'volume_flow_initial' )
	770	READ ( 13 ) volume_flow_initial
	771	CASE ( 'subs_vertical_gradient' )
	772	READ ( 13 ) subs_vertical_gradient
	773	CASE ( 'subs_vertical_gradient_level' )
	774	READ ( 13 ) subs_vertical_gradient_level
	775	CASE ( 'subs_vertical_gradient_level_i' )
	776	READ ( 13 ) subs_vertical_gradient_level_i
	777	CASE ( 'w_max' )
	778	READ ( 13 ) w_max
	779	CASE ( 'w_max_ijk' )
	780	READ ( 13 ) w_max_ijk
	781	CASE ( 'wall_adjustment' )
	782	READ ( 13 ) wall_adjustment
	783	CASE ( 'wall_heatflux' )
	784	READ ( 13 ) wall_heatflux
	785	CASE ( 'wall_humidityflux' )
	786	READ ( 13 ) wall_humidityflux
	787	CASE ( 'wall_scalarflux' )
	788	READ ( 13 ) wall_scalarflux
	789	CASE ( 'y_shift' )
	790	READ ( 13 ) y_shift
	791	CASE ( 'z0h_factor' )
	792	READ ( 13 ) z0h_factor
	793	CASE ( 'zeta_max' )
	794	READ ( 13 ) zeta_max
	795	CASE ( 'zeta_min' )
	796	READ ( 13 ) zeta_min
[3003]	797	CASE ( 'z_i' )
	798	READ ( 13 ) z_i
[2894]	799
	800	CASE DEFAULT
[3294]	801	!
	802	!-- Read global variables from of other modules
[3766]	803	CALL module_interface_rrd_global( found )
[2894]	804
	805	IF ( .NOT. found ) THEN
	806	WRITE( message_string, * ) 'unknown variable named "', &
	807	restart_string(1:length), &
	808	'" found in global data from ', &
	809	'prior run on PE ', myid
	810	CALL message( 'rrd_global', 'PA0302', 1, 2, 0, 6, 0 )
	811
	812	ENDIF
	813
	814	END SELECT
	815	!
	816	!-- Read next string
	817	READ ( 13 ) length
	818	READ ( 13 ) restart_string(1:length)
	819
	820	ENDDO
	821
	822
	823	CALL close_file( 13 )
	824
	825
	826	END SUBROUTINE rrd_global
	827
	828
	829
	830	!------------------------------------------------------------------------------!
	831	! Description:
	832	! ------------
	833	!> Skipping the global control variables from restart-file (binary format)
	834	!> except some information needed when reading restart data from a previous
	835	!> run which used a smaller total domain or/and a different domain decomposition
	836	!> (initializing_actions == 'cyclic_fill').
	837	!------------------------------------------------------------------------------!
	838
	839	SUBROUTINE rrd_read_parts_of_global
	840
	841
	842	CHARACTER (LEN=10) :: version_on_file
	843	CHARACTER (LEN=20) :: momentum_advec_check
	844	CHARACTER (LEN=20) :: scalar_advec_check
	845	CHARACTER (LEN=1) :: cdum
	846
	847	INTEGER(iwp) :: max_pr_user_on_file
	848	INTEGER(iwp) :: nz_on_file
	849	INTEGER(iwp) :: statistic_regions_on_file
	850	INTEGER(iwp) :: tmp_mpru
	851	INTEGER(iwp) :: tmp_sr
	852
	853	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: hom_sum_on_file
	854	REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: hom_on_file
	855
	856
	857	CALL check_open( 13 )
	858
	859	READ ( 13 ) length
	860	READ ( 13 ) restart_string(1:length)
	861	READ ( 13 ) version_on_file
	862
	863	!
	864	!-- Read number of PEs and horizontal index bounds of all PEs used in previous
	865	!-- run
	866	READ ( 13 ) length
	867	READ ( 13 ) restart_string(1:length)
	868
	869	IF ( restart_string(1:length) /= 'numprocs' ) THEN
	870	WRITE( message_string, * ) 'numprocs not found in data from prior ', &
	871	'run on PE ', myid
	872	CALL message( 'rrd_read_parts_of_global', 'PA0297', 1, 2, 0, 6, 0 )
	873	ENDIF
	874	READ ( 13 ) numprocs_previous_run
	875
	876	IF ( .NOT. ALLOCATED( hor_index_bounds_previous_run ) ) THEN
	877	ALLOCATE( hor_index_bounds_previous_run(4,0:numprocs_previous_run-1) )
	878	ENDIF
	879
	880	READ ( 13 ) length
	881	READ ( 13 ) restart_string(1:length)
	882
	883	IF ( restart_string(1:length) /= 'hor_index_bounds' ) THEN
	884	WRITE( message_string, * ) 'hor_index_bounds not found in data ', &
	885	'from prior run on PE ', myid
	886	CALL message( 'rrd_read_parts_of_global', 'PA0298', 1, 2, 0, 6, 0 )
	887	ENDIF
	888	READ ( 13 ) hor_index_bounds_previous_run
	889
	890	!
	891	!-- Read vertical number of gridpoints and number of different areas used
	892	!-- for computing statistics. Allocate arrays depending on these values,
	893	!-- which are needed for the following read instructions.
	894	READ ( 13 ) length
	895	READ ( 13 ) restart_string(1:length)
	896
	897	IF ( restart_string(1:length) /= 'nz' ) THEN
	898	message_string = 'nz not found in restart data file'
	899	CALL message( 'rrd_read_parts_of_global', 'PA0303', 1, 2, 0, 6, 0 )
	900	ENDIF
	901	READ ( 13 ) nz_on_file
	902	IF ( nz_on_file /= nz ) THEN
	903	WRITE( message_string, * ) 'mismatch concerning number of ', &
[3045]	904	'gridpoints along z:', &
[3046]	905	'&nz on file = "', nz_on_file, '"', &
	906	'&nz from run = "', nz, '"'
[2894]	907	CALL message( 'rrd_read_parts_of_global', 'PA0304', 1, 2, 0, 6, 0 )
	908	ENDIF
	909
	910	READ ( 13 ) length
	911	READ ( 13 ) restart_string(1:length)
	912
	913	IF ( restart_string(1:length) /= 'max_pr_user' ) THEN
	914	message_string = 'max_pr_user not found in restart data file'
	915	CALL message( 'rrd_read_parts_of_global', 'PA0305', 1, 2, 0, 6, 0 )
	916	ENDIF
	917	READ ( 13 ) max_pr_user_on_file
	918	IF ( max_pr_user_on_file /= max_pr_user ) THEN
	919	WRITE( message_string, * ) 'number of user profiles on res', &
	920	'tart data file differs from the ', &
[3046]	921	'current run:&max_pr_user on file = "',&
[2894]	922	max_pr_user_on_file, '"', &
[3046]	923	'&max_pr_user from run = "', &
[2894]	924	max_pr_user, '"'
	925	CALL message( 'rrd_read_parts_of_global', 'PA0306', 0, 0, 0, 6, 0 )
	926	tmp_mpru = MIN( max_pr_user_on_file, max_pr_user )
	927	ELSE
	928	tmp_mpru = max_pr_user
	929	ENDIF
	930
	931	READ ( 13 ) length
	932	READ ( 13 ) restart_string(1:length)
	933
	934	IF ( restart_string(1:length) /= 'statistic_regions' ) THEN
	935	message_string = 'statistic_regions not found in restart data file'
	936	CALL message( 'rrd_read_parts_of_global', 'PA0307', 1, 2, 0, 6, 0 )
	937	ENDIF
	938	READ ( 13 ) statistic_regions_on_file
	939	IF ( statistic_regions_on_file /= statistic_regions ) THEN
	940	WRITE( message_string, * ) 'statistic regions on restart data file ',&
[3045]	941	'differ from the current run:', &
[3046]	942	'&statistic regions on file = "', &
[2894]	943	statistic_regions_on_file, '"', &
[3046]	944	'&statistic regions from run = "', &
[2894]	945	statistic_regions, '"', &
[3046]	946	'&statistic data may be lost!'
[2894]	947	CALL message( 'rrd_read_parts_of_global', 'PA0308', 0, 1, 0, 6, 0 )
	948	tmp_sr = MIN( statistic_regions_on_file, statistic_regions )
	949	ELSE
	950	tmp_sr = statistic_regions
	951	ENDIF
	952
	953	!
	954	!-- Now read and check some control parameters and skip the rest
	955	READ ( 13 ) length
	956	READ ( 13 ) restart_string(1:length)
	957
	958	DO WHILE ( restart_string(1:length) /= 'binary_version_local' )
	959
	960	SELECT CASE ( restart_string(1:length) )
	961
	962	CASE ( 'average_count_pr' )
	963	READ ( 13 ) average_count_pr
	964	IF ( average_count_pr /= 0 ) THEN
	965	WRITE( message_string, * ) 'inflow profiles not ', &
[3046]	966	'temporally averaged. &Averaging will be ', &
[3637]	967	'done now using', average_count_pr, &
[2894]	968	' samples.'
	969	CALL message( 'rrd_read_parts_of_global', 'PA0309', &
	970	0, 1, 0, 6, 0 )
	971	ENDIF
	972
	973	CASE ( 'hom' )
	974	ALLOCATE( hom_on_file(0:nz+1,2,pr_palm+max_pr_user_on_file, &
	975	0:statistic_regions_on_file) )
	976	READ ( 13 ) hom_on_file
	977	hom(:,:,1:pr_palm+tmp_mpru,0:tmp_sr) = &
	978	hom_on_file(:,:,1:pr_palm+tmp_mpru,0:tmp_sr)
	979	DEALLOCATE( hom_on_file )
	980
	981	CASE ( 'hom_sum' )
	982	ALLOCATE( hom_sum_on_file(0:nz+1,pr_palm+max_pr_user_on_file, &
	983	0:statistic_regions_on_file) )
	984	READ ( 13 ) hom_sum_on_file
	985	hom_sum(:,1:pr_palm+tmp_mpru,0:tmp_sr) = &
	986	hom_sum_on_file(:,1:pr_palm+tmp_mpru,0:tmp_sr)
	987	DEALLOCATE( hom_sum_on_file )
	988
	989	CASE ( 'momentum_advec' )
	990	momentum_advec_check = momentum_advec
	991	READ ( 13 ) momentum_advec
	992	IF ( TRIM( momentum_advec_check ) /= TRIM( momentum_advec ) ) &
	993	THEN
	994	WRITE( message_string, * ) 'momentum_advec of the restart ',&
	995	'run differs from momentum_advec of the ', &
	996	'initial run.'
	997	CALL message( 'rrd_read_parts_of_global', 'PA0100', &
	998	1, 2, 0, 6, 0 )
	999	ENDIF
	1000
	1001	CASE ( 'nx' )
	1002	READ ( 13 ) nx_on_file
	1003
	1004	CASE ( 'ny' )
	1005	READ ( 13 ) ny_on_file
	1006
	1007	CASE ( 'ref_state' )
	1008	READ ( 13 ) ref_state
	1009
	1010	CASE ( 'scalar_advec' )
	1011	scalar_advec_check = scalar_advec
	1012	READ ( 13 ) scalar_advec
	1013	IF ( TRIM( scalar_advec_check ) /= TRIM( scalar_advec ) ) &
	1014	THEN
	1015	WRITE( message_string, * ) 'scalar_advec of the restart ', &
	1016	'run differs from scalar_advec of the ', &
	1017	'initial run.'
	1018	CALL message( 'rrd_read_parts_of_global', 'PA0101', &
	1019	1, 2, 0, 6, 0 )
	1020	ENDIF
	1021
	1022	CASE DEFAULT
	1023
	1024	READ ( 13 ) cdum
	1025
	1026	END SELECT
	1027
	1028	READ ( 13 ) length
	1029	READ ( 13 ) restart_string(1:length)
	1030
	1031	ENDDO
	1032
	1033	!
	1034	!-- Calculate the temporal average of vertical profiles, if neccessary
	1035	IF ( average_count_pr /= 0 ) THEN
	1036	hom_sum = hom_sum / REAL( average_count_pr, KIND=wp )
	1037	ENDIF
	1038
	1039
	1040	CALL close_file( 13 )
	1041
	1042
	1043	END SUBROUTINE rrd_read_parts_of_global
	1044
	1045
	1046	! Description:
	1047	! ------------
	1048	!> Reads processor specific data of variables and arrays from restart file
	1049	!> (binary format).
	1050	!------------------------------------------------------------------------------!
[3637]	1051	SUBROUTINE rrd_local
[2894]	1052
	1053
	1054	CHARACTER (LEN=7) :: myid_char_save
	1055	CHARACTER (LEN=10) :: binary_version_local
	1056	CHARACTER (LEN=10) :: version_on_file
	1057
	1058	INTEGER(iwp) :: files_to_be_opened !<
	1059	INTEGER(iwp) :: i !<
	1060	INTEGER(iwp) :: j !<
	1061	INTEGER(iwp) :: k !<
	1062	INTEGER(iwp) :: myid_on_file !<
	1063	INTEGER(iwp) :: numprocs_on_file !<
	1064	INTEGER(iwp) :: nxlc !<
	1065	INTEGER(iwp) :: nxlf !<
	1066	INTEGER(iwp) :: nxlpr !<
	1067	INTEGER(iwp) :: nxl_on_file !<
	1068	INTEGER(iwp) :: nxrc !<
	1069	INTEGER(iwp) :: nxrf !<
	1070	INTEGER(iwp) :: nxrpr !<
	1071	INTEGER(iwp) :: nxr_on_file !<
	1072	INTEGER(iwp) :: nync !<
	1073	INTEGER(iwp) :: nynf !<
	1074	INTEGER(iwp) :: nynpr !<
	1075	INTEGER(iwp) :: nyn_on_file !<
	1076	INTEGER(iwp) :: nysc !<
	1077	INTEGER(iwp) :: nysf !<
	1078	INTEGER(iwp) :: nyspr !<
	1079	INTEGER(iwp) :: nys_on_file !<
	1080	INTEGER(iwp) :: nzb_on_file !<
	1081	INTEGER(iwp) :: nzt_on_file !<
	1082	INTEGER(iwp) :: offset_x !<
	1083	INTEGER(iwp) :: offset_y !<
	1084	INTEGER(iwp) :: shift_x !<
	1085	INTEGER(iwp) :: shift_y !<
	1086
	1087	INTEGER(iwp), DIMENSION(numprocs_previous_run) :: file_list !<
	1088	INTEGER(iwp), DIMENSION(numprocs_previous_run) :: overlap_count !<
	1089
	1090	INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) :: nxlfa !<
	1091	INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) :: nxrfa !<
	1092	INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) :: nynfa !<
	1093	INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) :: nysfa !<
	1094	INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) :: offset_xa !<
	1095	INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) :: offset_ya !<
	1096
	1097	INTEGER(isp), DIMENSION(:,:), ALLOCATABLE :: tmp_2d_id_random !< temporary array for storing random generator data
	1098	INTEGER(isp), DIMENSION(:,:,:), ALLOCATABLE :: tmp_2d_seq_random !< temporary array for storing random generator data
	1099
	1100	LOGICAL :: found
	1101
[3998]	1102	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tmp_2d !< temporary array for storing 2D data
	1103	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: tmp_3d !< temporary array for storing 3D data
[4009]	1104	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: tmp_3d_non_standard !< temporary array for storing 3D data
	1105	!< with non standard dimensions
[2894]	1106
	1107	!
	1108	!-- Read data from previous model run.
	1109	CALL cpu_log( log_point_s(14), 'rrd_local', 'start' )
	1110
	1111	!
	1112	!-- Check which of the restart files contain data needed for the subdomain
	1113	!-- of this PE
	1114	files_to_be_opened = 0
	1115
	1116	DO i = 1, numprocs_previous_run
	1117	!
	1118	!-- Store array bounds of the previous run ("pr") in temporary scalars
	1119	nxlpr = hor_index_bounds_previous_run(1,i-1)
	1120	nxrpr = hor_index_bounds_previous_run(2,i-1)
	1121	nyspr = hor_index_bounds_previous_run(3,i-1)
	1122	nynpr = hor_index_bounds_previous_run(4,i-1)
	1123
	1124	!
	1125	!-- Determine the offsets. They may be non-zero in case that the total domain
	1126	!-- on file is smaller than the current total domain.
	1127	offset_x = ( nxl / ( nx_on_file + 1 ) ) * ( nx_on_file + 1 )
	1128	offset_y = ( nys / ( ny_on_file + 1 ) ) * ( ny_on_file + 1 )
	1129
	1130	!
	1131	!-- Start with this offset and then check, if the subdomain on file
	1132	!-- matches another time(s) in the current subdomain by shifting it
	1133	!-- for nx_on_file+1, ny_on_file+1 respectively
	1134
	1135	shift_y = 0
	1136	j = 0
	1137	DO WHILE ( nyspr+shift_y <= nyn-offset_y )
	1138
	1139	IF ( nynpr+shift_y >= nys-offset_y ) THEN
	1140
	1141	shift_x = 0
	1142	DO WHILE ( nxlpr+shift_x <= nxr-offset_x )
	1143
	1144	IF ( nxrpr+shift_x >= nxl-offset_x ) THEN
	1145	j = j +1
	1146	IF ( j > 1000 ) THEN
	1147	!
	1148	!-- Array bound exceeded
[4009]	1149	message_string = 'data from subdomain of previous' // &
[2894]	1150	' run mapped more than 1000 times'
[4009]	1151	CALL message( 'rrd_local', 'PA0284', 2, 2, -1, 6, 1 )
[2894]	1152	ENDIF
	1153
	1154	IF ( j == 1 ) THEN
	1155	files_to_be_opened = files_to_be_opened + 1
	1156	file_list(files_to_be_opened) = i-1
	1157	ENDIF
	1158
	1159	offset_xa(files_to_be_opened,j) = offset_x + shift_x
	1160	offset_ya(files_to_be_opened,j) = offset_y + shift_y
	1161	!
	1162	!-- Index bounds of overlapping data
[4009]	1163	nxlfa(files_to_be_opened,j) = MAX( nxl-offset_x-shift_x, nxlpr )
	1164	nxrfa(files_to_be_opened,j) = MIN( nxr-offset_x-shift_x, nxrpr )
	1165	nysfa(files_to_be_opened,j) = MAX( nys-offset_y-shift_y, nyspr )
	1166	nynfa(files_to_be_opened,j) = MIN( nyn-offset_y-shift_y, nynpr )
[2894]	1167
	1168	ENDIF
	1169
	1170	shift_x = shift_x + ( nx_on_file + 1 )
	1171	ENDDO
	1172
	1173	ENDIF
	1174
	1175	shift_y = shift_y + ( ny_on_file + 1 )
	1176	ENDDO
	1177
	1178	IF ( j > 0 ) overlap_count(files_to_be_opened) = j
	1179
	1180	ENDDO
	1181
	1182	!
	1183	!-- Save the id-string of the current process, since myid_char may now be used
	1184	!-- to open files created by PEs with other id.
	1185	myid_char_save = myid_char
	1186
[4009]	1187	IF ( files_to_be_opened /= 1 .OR. numprocs /= numprocs_previous_run ) THEN
	1188	WRITE( message_string, * ) 'number of PEs or virtual PE-grid changed ', &
	1189	'in restart run & PE', myid, ' will read from files ', &
	1190	file_list(1:files_to_be_opened)
[2894]	1191	CALL message( 'rrd_local', 'PA0285', 0, 0, 0, 6, 0 )
	1192	ENDIF
	1193
	1194	!
	1195	!-- Read data from all restart files determined above
	1196	DO i = 1, files_to_be_opened
	1197
	1198	j = file_list(i)
	1199	!
	1200	!-- Set the filename (underscore followed by four digit processor id)
	1201	WRITE (myid_char,'(''_'',I6.6)') j
	1202
	1203	!
	1204	!-- Open the restart file. If this file has been created by PE0 (_000000),
	1205	!-- the global variables at the beginning of the file have to be skipped
	1206	!-- first.
	1207	CALL check_open( 13 )
	1208	IF ( j == 0 ) CALL rrd_skip_global
	1209
	1210	!
	1211	!-- First compare the version numbers
	1212	READ ( 13 ) length
	1213	READ ( 13 ) restart_string(1:length)
	1214	READ ( 13 ) version_on_file
	1215
	1216	binary_version_local = '4.7'
	1217	IF ( TRIM( version_on_file ) /= TRIM( binary_version_local ) ) THEN
[4009]	1218	WRITE( message_string, * ) 'version mismatch concerning ', &
	1219	'binary_version_local:', &
	1220	'&version on file = "', TRIM( version_on_file ), '"', &
	1221	'&version in program = "', TRIM( binary_version_local ), '"'
[2894]	1222	CALL message( 'rrd_local', 'PA0286', 1, 2, 0, 6, 0 )
	1223	ENDIF
	1224
	1225	!
	1226	!-- Read number of processors, processor-id, and array ranges.
	1227	!-- Compare the array ranges with those stored in the index bound array.
[4009]	1228	READ ( 13 ) numprocs_on_file, myid_on_file, nxl_on_file, nxr_on_file, nys_on_file, &
	1229	nyn_on_file, nzb_on_file, nzt_on_file
[2894]	1230
	1231	IF ( nxl_on_file /= hor_index_bounds_previous_run(1,j) ) THEN
[4009]	1232	WRITE( message_string, * ) 'problem with index bound nxl on ', &
	1233	'restart file "', myid_char, '"', &
	1234	'&nxl = ', nxl_on_file, ' but it should be', &
	1235	'&= ', hor_index_bounds_previous_run(1,j), &
	1236	'&from the index bound information array'
[2894]	1237	CALL message( 'rrd_local', 'PA0287', 2, 2, -1, 6, 1 )
	1238	ENDIF
	1239
	1240	IF ( nxr_on_file /= hor_index_bounds_previous_run(2,j) ) THEN
[4009]	1241	WRITE( message_string, * ) 'problem with index bound nxr on ', &
	1242	'restart file "', myid_char, '"' , &
	1243	' nxr = ', nxr_on_file, ' but it should be', &
	1244	' = ', hor_index_bounds_previous_run(2,j), &
	1245	' from the index bound information array'
[2894]	1246	CALL message( 'rrd_local', 'PA0288', 2, 2, -1, 6, 1 )
	1247
	1248	ENDIF
	1249
	1250	IF ( nys_on_file /= hor_index_bounds_previous_run(3,j) ) THEN
[4009]	1251	WRITE( message_string, * ) 'problem with index bound nys on ', &
	1252	'restart file "', myid_char, '"', &
	1253	'&nys = ', nys_on_file, ' but it should be', &
	1254	'&= ', hor_index_bounds_previous_run(3,j), &
	1255	'&from the index bound information array'
[2894]	1256	CALL message( 'rrd_local', 'PA0289', 2, 2, -1, 6, 1 )
	1257	ENDIF
	1258
	1259	IF ( nyn_on_file /= hor_index_bounds_previous_run(4,j) ) THEN
[4009]	1260	WRITE( message_string, * ) 'problem with index bound nyn on ', &
	1261	'restart file "', myid_char, '"', &
	1262	'&nyn = ', nyn_on_file, ' but it should be', &
	1263	'&= ', hor_index_bounds_previous_run(4,j), &
	1264	'&from the index bound information array'
[2894]	1265	CALL message( 'rrd_local', 'PA0290', 2, 2, -1, 6, 1 )
	1266	ENDIF
	1267
	1268	IF ( nzb_on_file /= nzb ) THEN
[4009]	1269	WRITE( message_string, * ) 'mismatch between actual data and data ', &
	1270	'from prior run on PE ', myid, &
	1271	'&nzb on file = ', nzb_on_file, &
[3046]	1272	'&nzb = ', nzb
[2894]	1273	CALL message( 'rrd_local', 'PA0291', 1, 2, 0, 6, 0 )
	1274	ENDIF
	1275
	1276	IF ( nzt_on_file /= nzt ) THEN
[4009]	1277	WRITE( message_string, * ) 'mismatch between actual data and data ', &
	1278	'from prior run on PE ', myid, &
	1279	'&nzt on file = ', nzt_on_file, &
[3046]	1280	'&nzt = ', nzt
[2894]	1281	CALL message( 'rrd_local', 'PA0292', 1, 2, 0, 6, 0 )
	1282	ENDIF
	1283
	1284	!
	1285	!-- Allocate temporary arrays sized as the arrays on the restart file
[4009]	1286	ALLOCATE( tmp_2d(nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp), &
	1287	tmp_3d(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &
	1288	nxl_on_file-nbgp:nxr_on_file+nbgp) )
[2894]	1289
	1290	!
	1291	!-- Read arrays
	1292	!-- ATTENTION: If the following read commands have been altered, the
	1293	!-- ---------- version number of the variable binary_version_local must
	1294	!-- be altered, too. Furthermore, the output list of arrays in
	1295	!-- wrd_write_local must also be altered
	1296	!-- accordingly.
	1297	READ ( 13 ) length
	1298	READ ( 13 ) restart_string(1:length)
	1299
	1300
	1301	!
	1302	!-- Loop over processor specific field data
	1303	DO WHILE ( restart_string(1:length) /= '* end *' )
	1304
	1305	!
	1306	!-- Map data on file as often as needed (data are read only for k=1)
	1307	DO k = 1, overlap_count(i)
	1308
[3056]	1309	found = .FALSE.
	1310
[2894]	1311	!
	1312	!-- Get the index range of the subdomain on file which overlap with
	1313	!-- the current subdomain
	1314	nxlf = nxlfa(i,k)
	1315	nxlc = nxlfa(i,k) + offset_xa(i,k)
	1316	nxrf = nxrfa(i,k)
	1317	nxrc = nxrfa(i,k) + offset_xa(i,k)
	1318	nysf = nysfa(i,k)
	1319	nysc = nysfa(i,k) + offset_ya(i,k)
	1320	nynf = nynfa(i,k)
	1321	nync = nynfa(i,k) + offset_ya(i,k)
	1322
	1323
	1324	SELECT CASE ( restart_string(1:length) )
	1325
	1326	CASE ( 'ghf_av' )
	1327	IF ( .NOT. ALLOCATED( ghf_av ) ) THEN
	1328	ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
	1329	ENDIF
	1330	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1331	ghf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1332	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1333
	1334	CASE ( 'e' )
	1335	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1336	e(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1337	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1338
	1339	CASE ( 'e_av' )
	1340	IF ( .NOT. ALLOCATED( e_av ) ) THEN
[4009]	1341	ALLOCATE( e_av(nzb:nzt+1,nys-nbgp:nyn+nbgp, &
[2894]	1342	nxl-nbgp:nxr+nbgp) )
	1343	ENDIF
	1344	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1345	e_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1346	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1347
	1348	CASE ( 'kh' )
	1349	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1350	kh(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1351	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1352
	1353	CASE ( 'kh_av' )
	1354	IF ( .NOT. ALLOCATED( kh_av ) ) THEN
	1355	ALLOCATE( kh_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ))
	1356	ENDIF
	1357	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1358	kh_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1359	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1360
	1361	CASE ( 'km' )
	1362	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1363	km(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1364	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1365
	1366	CASE ( 'km_av' )
	1367	IF ( .NOT. ALLOCATED( km_av ) ) THEN
	1368	ALLOCATE( km_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ))
	1369	ENDIF
	1370	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1371	km_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1372	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1373
	1374	CASE ( 'lpt_av' )
	1375	IF ( .NOT. ALLOCATED( lpt_av ) ) THEN
	1376	ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ))
	1377	ENDIF
	1378	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1379	lpt_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1380	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1381
	1382	CASE ( 'lwp_av' )
	1383	IF ( .NOT. ALLOCATED( lwp_av ) ) THEN
	1384	ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
	1385	ENDIF
	1386	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1387	lwp_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1388	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1389
	1390	CASE ( 'p' )
	1391	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1392	p(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1393	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1394
	1395	CASE ( 'p_av' )
	1396	IF ( .NOT. ALLOCATED( p_av ) ) THEN
	1397	ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	1398	ENDIF
	1399	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1400	p_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1401	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1402
	1403	CASE ( 'pt' )
	1404	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1405	pt(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1406	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1407
	1408	CASE ( 'pt_av' )
	1409	IF ( .NOT. ALLOCATED( pt_av ) ) THEN
	1410	ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	1411	ENDIF
	1412	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1413	pt_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1414	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1415
[3936]	1416	CASE ( 'pt_2m_av' )
	1417	IF ( .NOT. ALLOCATED( pt_2m_av ) ) THEN
	1418	ALLOCATE( pt_2m_av(nysg:nyng,nxlg:nxrg) )
	1419	ENDIF
	1420	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1421	pt_2m_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[3936]	1422	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1423
[2894]	1424	CASE ( 'q' )
	1425	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1426	q(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1427	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1428
	1429	CASE ( 'q_av' )
	1430	IF ( .NOT. ALLOCATED( q_av ) ) THEN
	1431	ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ))
	1432	ENDIF
	1433	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1434	q_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1435	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1436
	1437	CASE ( 'ql' )
	1438	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1439	ql(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1440	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1441
	1442	CASE ( 'ql_av' )
	1443	IF ( .NOT. ALLOCATED( ql_av ) ) THEN
	1444	ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	1445	ENDIF
	1446	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1447	ql_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1448	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1449
	1450	CASE ( 'qsws_av' )
	1451	IF ( .NOT. ALLOCATED( qsws_av ) ) THEN
	1452	ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
	1453	ENDIF
	1454	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1455	qsws_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1456	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1457
	1458	CASE ( 'qv_av' )
	1459	IF ( .NOT. ALLOCATED( qv_av ) ) THEN
	1460	ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	1461	ENDIF
	1462	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1463	qv_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1464	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1465
	1466	CASE ( 'r_a_av' )
	1467	IF ( .NOT. ALLOCATED( r_a_av ) ) THEN
	1468	ALLOCATE( r_a_av(nysg:nyng,nxlg:nxrg) )
	1469	ENDIF
	1470	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1471	r_a_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1472	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1473
	1474	CASE ( 'random_iv' ) ! still unresolved issue
	1475	IF ( k == 1 ) READ ( 13 ) random_iv
	1476	IF ( k == 1 ) READ ( 13 ) random_iy
	1477
	1478	CASE ( 'seq_random_array' )
[4009]	1479	ALLOCATE( tmp_2d_id_random(nys_on_file:nyn_on_file,nxl_on_file:nxr_on_file) )
	1480	ALLOCATE( tmp_2d_seq_random(5,nys_on_file:nyn_on_file,nxl_on_file:nxr_on_file) )
[2894]	1481	IF ( .NOT. ALLOCATED( id_random_array ) ) THEN
	1482	ALLOCATE( id_random_array(nys:nyn,nxl:nxr) )
	1483	ENDIF
	1484	IF ( .NOT. ALLOCATED( seq_random_array ) ) THEN
	1485	ALLOCATE( seq_random_array(5,nys:nyn,nxl:nxr) )
	1486	ENDIF
	1487	IF ( k == 1 ) READ ( 13 ) tmp_2d_id_random
	1488	IF ( k == 1 ) READ ( 13 ) tmp_2d_seq_random
[4009]	1489	id_random_array(nysc:nync,nxlc:nxrc) = tmp_2d_id_random(nysf:nynf,nxlf:nxrf)
	1490	seq_random_array(:,nysc:nync,nxlc:nxrc) = tmp_2d_seq_random(:,nysf:nynf,nxlf:nxrf)
[2894]	1491	DEALLOCATE( tmp_2d_id_random, tmp_2d_seq_random )
	1492
	1493	CASE ( 's' )
	1494	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1495	s(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1496	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1497
	1498	CASE ( 's_av' )
	1499	IF ( .NOT. ALLOCATED( s_av ) ) THEN
	1500	ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg))
	1501	ENDIF
	1502	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1503	s_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1504	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1505
	1506	CASE ( 'shf_av' )
	1507	IF ( .NOT. ALLOCATED( shf_av ) ) THEN
	1508	ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
	1509	ENDIF
	1510	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1511	shf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1512	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1513
	1514	CASE ( 'ssws_av' )
	1515	IF ( .NOT. ALLOCATED( ssws_av ) ) THEN
	1516	ALLOCATE( ssws_av(nysg:nyng,nxlg:nxrg) )
	1517	ENDIF
	1518	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1519	ssws_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1520	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
[3994]	1521
	1522	CASE ( 'ti_av' )
	1523	IF ( .NOT. ALLOCATED( ti_av ) ) THEN
	1524	ALLOCATE( ti_av(nzb:nzt+1,nys:nyn,nxl:nxr) )
	1525	ENDIF
[4009]	1526	IF ( k == 1 ) THEN
	1527	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &
	1528	nxl_on_file:nxr_on_file) )
	1529	READ ( 13 ) tmp_3d_non_standard
	1530	ENDIF
	1531	ti_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)
[3994]	1532
[2894]	1533	CASE ( 'ts_av' )
	1534	IF ( .NOT. ALLOCATED( ts_av ) ) THEN
	1535	ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
	1536	ENDIF
	1537	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1538	ts_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1539	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1540
	1541	CASE ( 'tsurf_av' )
	1542	IF ( .NOT. ALLOCATED( tsurf_av ) ) THEN
	1543	ALLOCATE( tsurf_av(nysg:nyng,nxlg:nxrg) )
	1544	ENDIF
	1545	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1546	tsurf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1547	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1548
	1549	CASE ( 'u' )
	1550	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1551	u(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1552	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1553
	1554	CASE ( 'u_av' )
	1555	IF ( .NOT. ALLOCATED( u_av ) ) THEN
	1556	ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	1557	ENDIF
	1558	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1559	u_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1560	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
[4039]	1561
	1562	CASE ( 'uu_av' )
	1563	IF ( .NOT. ALLOCATED( uu_av ) ) THEN
	1564	ALLOCATE( uu_av(nzb:nzt+1,nys:nyn,nxl:nxr) )
	1565	ENDIF
	1566	IF ( k == 1 ) THEN
	1567	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &
	1568	nxl_on_file:nxr_on_file) )
	1569	READ ( 13 ) tmp_3d_non_standard
	1570	ENDIF
	1571	uu_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)
[2894]	1572
	1573	CASE ( 'u_m_l' )
	1574	IF ( k == 1 ) THEN
[4009]	1575	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &
	1576	1:2) )
	1577	READ ( 13 ) tmp_3d_non_standard
[2894]	1578	ENDIF
[3182]	1579	IF ( bc_radiation_l ) THEN
[4009]	1580	u_m_l(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
[2894]	1581	ENDIF
	1582
	1583	CASE ( 'u_m_n' )
	1584	IF ( k == 1 ) THEN
[4009]	1585	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,ny-1:ny, &
	1586	nxl_on_file-nbgp:nxr_on_file+nbgp) )
	1587	READ ( 13 ) tmp_3d_non_standard
[2894]	1588	ENDIF
[3182]	1589	IF ( bc_radiation_n ) THEN
[4009]	1590	u_m_n(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
[2894]	1591	ENDIF
	1592
	1593	CASE ( 'u_m_r' )
	1594	IF ( k == 1 ) THEN
[4009]	1595	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &
	1596	nx-1:nx) )
	1597	READ ( 13 ) tmp_3d_non_standard
[2894]	1598	ENDIF
[3182]	1599	IF ( bc_radiation_r ) THEN
[4009]	1600	u_m_r(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
[2894]	1601	ENDIF
	1602
	1603	CASE ( 'u_m_s' )
	1604	IF ( k == 1 ) THEN
[4009]	1605	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,0:1, &
	1606	nxl_on_file-nbgp:nxr_on_file+nbgp) )
	1607	READ ( 13 ) tmp_3d_non_standard
[2894]	1608	ENDIF
[3182]	1609	IF ( bc_radiation_s ) THEN
[4009]	1610	u_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
[2894]	1611	ENDIF
	1612
	1613	CASE ( 'us_av' )
	1614	IF ( .NOT. ALLOCATED( us_av ) ) THEN
	1615	ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
	1616	ENDIF
	1617	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1618	us_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1619	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1620
	1621	CASE ( 'v' )
	1622	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1623	v(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1624	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1625
	1626	CASE ( 'v_av' )
	1627	IF ( .NOT. ALLOCATED( v_av ) ) THEN
	1628	ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	1629	ENDIF
	1630	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1631	v_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1632	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
[4039]	1633
	1634	CASE ( 'vv_av' )
	1635	IF ( .NOT. ALLOCATED( vv_av ) ) THEN
	1636	ALLOCATE( vv_av(nzb:nzt+1,nys:nyn,nxl:nxr) )
	1637	ENDIF
	1638	IF ( k == 1 ) THEN
	1639	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &
	1640	nxl_on_file:nxr_on_file) )
	1641	READ ( 13 ) tmp_3d_non_standard
	1642	ENDIF
	1643	vv_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)
[2894]	1644
	1645	CASE ( 'v_m_l' )
	1646	IF ( k == 1 ) THEN
[4009]	1647	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &
	1648	0:1) )
	1649	READ ( 13 ) tmp_3d_non_standard
[2894]	1650	ENDIF
[3182]	1651	IF ( bc_radiation_l ) THEN
[4009]	1652	v_m_l(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
[2894]	1653	ENDIF
	1654
	1655	CASE ( 'v_m_n' )
	1656	IF ( k == 1 ) THEN
[4009]	1657	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,ny-1:ny, &
	1658	nxl_on_file-nbgp:nxr_on_file+nbgp) )
	1659	READ ( 13 ) tmp_3d_non_standard
[2894]	1660	ENDIF
[3182]	1661	IF ( bc_radiation_n ) THEN
[4009]	1662	v_m_n(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
[2894]	1663	ENDIF
	1664
	1665	CASE ( 'v_m_r' )
	1666	IF ( k == 1 ) THEN
[4009]	1667	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &
	1668	nx-1:nx) )
	1669	READ ( 13 ) tmp_3d_non_standard
[2894]	1670	ENDIF
[3182]	1671	IF ( bc_radiation_r ) THEN
[4009]	1672	v_m_r(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
[2894]	1673	ENDIF
	1674
	1675	CASE ( 'v_m_s' )
	1676	IF ( k == 1 ) THEN
[4009]	1677	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,1:2, &
	1678	nxl_on_file-nbgp:nxr_on_file+nbgp) )
	1679	READ ( 13 ) tmp_3d_non_standard
[2894]	1680	ENDIF
[3182]	1681	IF ( bc_radiation_s ) THEN
[4009]	1682	v_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
[2894]	1683	ENDIF
	1684
	1685	CASE ( 'vpt' )
	1686	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1687	vpt(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1688	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1689
	1690	CASE ( 'vpt_av' )
	1691	IF ( .NOT. ALLOCATED( vpt_av ) ) THEN
	1692	ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	1693	ENDIF
	1694	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1695	vpt_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1696	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1697
	1698	CASE ( 'w' )
	1699	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1700	w(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1701	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1702
	1703	CASE ( 'w_av' )
	1704	IF ( .NOT. ALLOCATED( w_av ) ) THEN
	1705	ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
	1706	ENDIF
	1707	IF ( k == 1 ) READ ( 13 ) tmp_3d
[4009]	1708	w_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1709	tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
[4039]	1710
	1711	CASE ( 'ww_av' )
	1712	IF ( .NOT. ALLOCATED( ww_av ) ) THEN
	1713	ALLOCATE( ww_av(nzb:nzt+1,nys:nyn,nxl:nxr) )
	1714	ENDIF
	1715	IF ( k == 1 ) THEN
	1716	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &
	1717	nxl_on_file:nxr_on_file) )
	1718	READ ( 13 ) tmp_3d_non_standard
	1719	ENDIF
	1720	ww_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)
[2894]	1721
	1722	CASE ( 'w_m_l' )
	1723	IF ( k == 1 ) THEN
[4009]	1724	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &
	1725	0:1) )
	1726	READ ( 13 ) tmp_3d_non_standard
[2894]	1727	ENDIF
[3182]	1728	IF ( bc_radiation_l ) THEN
[4009]	1729	w_m_l(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
[2894]	1730	ENDIF
	1731
	1732	CASE ( 'w_m_n' )
	1733	IF ( k == 1 ) THEN
[4009]	1734	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,ny-1:ny, &
	1735	nxl_on_file-nbgp:nxr_on_file+nbgp) )
	1736	READ ( 13 ) tmp_3d_non_standard
[2894]	1737	ENDIF
[3182]	1738	IF ( bc_radiation_n ) THEN
[4009]	1739	w_m_n(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
[2894]	1740	ENDIF
	1741
	1742	CASE ( 'w_m_r' )
	1743	IF ( k == 1 ) THEN
[4009]	1744	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &
	1745	nx-1:nx) )
	1746	READ ( 13 ) tmp_3d_non_standard
[2894]	1747	ENDIF
[3182]	1748	IF ( bc_radiation_r ) THEN
[4009]	1749	w_m_r(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)
[2894]	1750	ENDIF
	1751
	1752	CASE ( 'w_m_s' )
	1753	IF ( k == 1 ) THEN
[4009]	1754	ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,0:1, &
	1755	nxl_on_file-nbgp:nxr_on_file+nbgp) )
	1756	READ ( 13 ) tmp_3d_non_standard
[2894]	1757	ENDIF
[3182]	1758	IF ( bc_radiation_s ) THEN
[4009]	1759	w_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)
[2894]	1760	ENDIF
	1761
	1762	CASE ( 'z0_av' )
	1763	IF ( .NOT. ALLOCATED( z0_av ) ) THEN
	1764	ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
	1765	ENDIF
	1766	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1767	z0_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1768	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1769
	1770	CASE ( 'z0h_av' )
	1771	IF ( .NOT. ALLOCATED( z0h_av ) ) THEN
	1772	ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
	1773	ENDIF
	1774	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1775	z0h_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1776	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1777
	1778	CASE ( 'z0q_av' )
	1779	IF ( .NOT. ALLOCATED( z0q_av ) ) THEN
	1780	ALLOCATE( z0q_av(nysg:nyng,nxlg:nxrg) )
	1781	ENDIF
	1782	IF ( k == 1 ) READ ( 13 ) tmp_2d
[4009]	1783	z0q_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
[2894]	1784	tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
	1785
	1786	CASE DEFAULT
	1787
	1788	!
[3637]	1789	!-- Read restart data of surfaces
[4009]	1790	IF ( .NOT. found ) CALL surface_rrd_local( k, nxlf, nxlc, nxl_on_file, nxrf, &
	1791	nxr_on_file, nynf, nyn_on_file, &
	1792	nysf, nysc, nys_on_file, found )
[2912]	1793	!
[3637]	1794	!-- Read restart data of other modules
[4009]	1795	IF ( .NOT. found ) CALL module_interface_rrd_local( &
	1796	k, nxlf, nxlc, nxl_on_file, nxrf, &
	1797	nxrc, nxr_on_file, nynf, nync, &
	1798	nyn_on_file, nysf, nysc, &
	1799	nys_on_file, tmp_2d, tmp_3d, found )
[2894]	1800
	1801
	1802	IF ( .NOT. found ) THEN
[4009]	1803	WRITE( message_string, * ) 'unknown variable named "', &
	1804	restart_string(1:length), &
	1805	'" found in subdomain data ', &
[2894]	1806	'from prior run on PE ', myid
	1807	CALL message( 'rrd_local', 'PA0302', 1, 2, 0, 6, 0 )
	1808
	1809	ENDIF
	1810
	1811	END SELECT
	1812
	1813	ENDDO ! overlaploop
	1814
	1815	!
[4009]	1816	!-- Deallocate non standard array needed for specific variables only
	1817	IF ( ALLOCATED( tmp_3d_non_standard ) ) DEALLOCATE( tmp_3d_non_standard )
[2894]	1818
	1819	!
	1820	!-- Read next character string
	1821	READ ( 13 ) length
	1822	READ ( 13 ) restart_string(1:length)
	1823
	1824	ENDDO ! dataloop
	1825
	1826	!
	1827	!-- Close the restart file
	1828	CALL close_file( 13 )
	1829
	1830	DEALLOCATE( tmp_2d, tmp_3d )
	1831
	1832	ENDDO ! loop over restart files
	1833
	1834	!
	1835	!-- Restore the original filename for the restart file to be written
	1836	myid_char = myid_char_save
	1837
	1838	!
	1839	!-- End of time measuring for reading binary data
	1840	CALL cpu_log( log_point_s(14), 'rrd_local', 'stop' )
	1841
	1842	END SUBROUTINE rrd_local
	1843
	1844
	1845	!------------------------------------------------------------------------------!
	1846	! Description:
	1847	! ------------
	1848	!> Skipping the global control variables from restart-file (binary format)
	1849	!------------------------------------------------------------------------------!
	1850
	1851	SUBROUTINE rrd_skip_global
	1852
	1853
	1854	CHARACTER (LEN=1) :: cdum
	1855
	1856
	1857	READ ( 13 ) length
	1858	READ ( 13 ) restart_string(1:length)
	1859
	1860	DO WHILE ( restart_string(1:length) /= 'binary_version_local' )
	1861
	1862	READ ( 13 ) cdum
[3637]	1863	READ ( 13 ) length
[2894]	1864	READ ( 13 ) restart_string(1:length)
	1865
	1866	ENDDO
	1867
	1868	BACKSPACE ( 13 )
	1869	BACKSPACE ( 13 )
	1870
	1871
	1872	END SUBROUTINE rrd_skip_global
	1873
	1874
	1875	END MODULE read_restart_data_mod

Note: See TracBrowser for help on using the repository browser.

Context Navigation

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

Download in other formats: