Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |