Home

Context Navigation

source: palm/trunk/SOURCE/write_restart_data_mod.f90 @ 4802

Last change on this file since 4802 was 4802, checked in by suehring, 4 years ago
Bugfix in time-control of indoor model in case of restarts and at the very first time step
Property svn:keywords set to `Id`
File size: 55.0 KB

Line
1	!> @file write_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 terms of the GNU General
6	! Public License as published by the Free Software Foundation, either version 3 of the License, or
7	! (at your option) any later version.
8	!
9	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the
10	! implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
11	! Public License for more details.
12	!
13	! You should have received a copy of the GNU General Public License along with PALM. If not, see
14	! <http://www.gnu.org/licenses/>.
15	!
16	! Copyright 1997-2020 Leibniz Universitaet Hannover
17	!--------------------------------------------------------------------------------------------------!
18	!
19	! Current revisions:
20	! -----------------
21	!
22	!
23	! Former revisions:
24	! -----------------
25	! $Id: write_restart_data_mod.f90 4802 2020-11-30 14:16:31Z suehring $
26	! include time_indoor into restart mechanism
27	!
28	! 4777 2020-11-06 14:50:49Z raasch
29	! bugfix for writing spectra data with MPI-I/O (they are global arrays!)
30	!
31	! 4671 2020-09-09 20:27:58Z pavelkrc
32	! Implementation of downward facing USM and LSM surfaces
33	!
34	! 4564 2020-06-12 14:03:36Z raasch
35	! Vertical nesting method of Huq et al. (2019) removed
36	!
37	! 4536 2020-05-17 17:24:13Z raasch
38	! binary version incremented
39	!
40	! 4535 2020-05-15 12:07:23Z raasch
41	! bugfix for restart data format query
42	!
43	! 4514 2020-04-30 16:29:59Z suehring
44	! qsurf and ssurf added
45	!
46	! 4495 2020-04-13 20:11:20Z raasch
47	! restart data handling with MPI-IO added
48	!
49	! 4489 2020-04-04 10:54:08Z raasch
50	! file re-formatted to follow the PALM coding standard
51	!
52	! 4360 2020-01-07 11:25:50Z suehring
53	! - Move 2-m potential temperature to diagnostic_output_quantities.
54	!
55	! 4301 2019-11-22 12:09:09Z oliver.maas
56	! removed recycling_yshift
57	!
58	! 4227 2019-09-10 18:04:34Z gronemeier
59	! implement new palm_date_time_mod and increased binary version
60	!
61	! 4146 2019-08-07 07:47:36Z gronemeier
62	! Corrected "Former revisions" section
63	!
64	! 4101 2019-07-17 15:14:26Z gronemeier
65	! remove old_dt
66	!
67	! 4039 2019-06-18 10:32:41Z suehring
68	! Modularize diagnostic output
69	!
70	! 4017 2019-06-06 12:16:46Z schwenkel
71	! output of turbulence intensity added
72	!
73	! 3988 2019-05-22 11:32:37Z kanani
74	! + time_virtual_measurement (to enable steering of output interval)
75	!
76	! 3936 2019-04-26 15:38:02Z kanani
77	! Enable time-averaged output of theta_2m* with restarts
78	!
79	! 3668 2019-01-14 12:49:24Z maronga
80	! Implementation of the PALM module interface
81	!
82	! 2894 2018-03-15 09:17:58Z Giersch
83	! Initial revision
84	!
85	!
86	! Description:
87	! ------------
88	!> Writes restart data into binary file(s) for restart runs.
89	!--------------------------------------------------------------------------------------------------!
90	MODULE write_restart_data_mod
91
92
93	USE arrays_3d, &
94	ONLY: e, inflow_damping_factor, kh, km, mean_inflow_profiles, p, pt, pt_init, q, q_init, &
95	ql, ref_state, s, s_init, u, u_init, ug, u_m_l, u_m_n, u_m_r, u_m_s, v, v_init, vg, &
96	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
97
98	USE averaging
99
100	USE bulk_cloud_model_mod, &
101	ONLY: bulk_cloud_model
102
103	USE control_parameters
104
105	USE grid_variables, &
106	ONLY: dx, dy
107
108	USE gust_mod, &
109	ONLY: gust_module_enabled
110
111	USE indices, &
112	ONLY: nx, nxl, nxr, ny, nys, nyn, nz, nzb, nzt
113
114	USE indoor_model_mod, &
115	ONLY: time_indoor
116
117	USE kinds
118
119	USE model_1d_mod, &
120	ONLY: damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
121
122	USE module_interface, &
123	ONLY: module_interface_wrd_global, &
124	module_interface_wrd_local
125
126	USE netcdf_interface, &
127	ONLY: netcdf_precision, output_for_t0
128
129	USE particle_attributes, &
130	ONLY: particle_advection
131
132	USE pegrid, &
133	ONLY: collective_wait, hor_index_bounds, myid, numprocs
134
135	USE radiation_model_mod, &
136	ONLY: time_radiation
137
138	USE random_function_mod, &
139	ONLY: random_iv, random_iy
140
141	USE random_generator_parallel, &
142	ONLY: id_random_array, seq_random_array
143
144	USE restart_data_mpi_io_mod, &
145	ONLY: wrd_mpi_io, wrd_mpi_io_global_array
146
147	USE spectra_mod, &
148	ONLY: average_count_sp, spectrum_x, spectrum_y
149
150	USE statistics, &
151	ONLY: hom, hom_sum, statistic_regions, u_max, u_max_ijk, v_max, v_max_ijk, w_max, &
152	w_max_ijk, z_i
153
154	USE surface_mod, &
155	ONLY: surface_wrd_local
156
157	USE user, &
158	ONLY: user_module_enabled
159
160	USE virtual_measurement_mod, &
161	ONLY: time_virtual_measurement
162
163
164	IMPLICIT NONE
165
166
167	INTERFACE wrd_global
168	MODULE PROCEDURE wrd_global
169	END INTERFACE wrd_global
170
171	INTERFACE wrd_local
172	MODULE PROCEDURE wrd_local
173	END INTERFACE wrd_local
174
175
176	PUBLIC wrd_local, wrd_global
177
178
179	CONTAINS
180
181
182	! Description:
183	! ------------
184	!> Global data of control variables and arrays is written out for restarts.
185	!> In case of output in Fortran binary format, this information is only written to the file
186	!> opened by PE0.
187	!--------------------------------------------------------------------------------------------------!
188
189
190	SUBROUTINE wrd_global
191
192
193	CHARACTER(LEN=10) :: binary_version_global !<
194	CHARACTER(LEN=20) :: tmp_name !< temporary variable
195
196	INTEGER :: i !< loop index
197
198	binary_version_global = '5.1'
199
200	IF ( restart_data_format_output == 'fortran_binary' ) THEN
201	!
202	!-- Output in Fortran binary format
203	CALL wrd_write_string( 'binary_version_global' )
204	WRITE ( 14 ) binary_version_global
205
206	CALL wrd_write_string( 'numprocs' )
207	WRITE ( 14 ) numprocs
208
209	CALL wrd_write_string( 'hor_index_bounds' )
210	WRITE ( 14 ) hor_index_bounds
211
212	CALL wrd_write_string( 'nz' )
213	WRITE ( 14 ) nz
214
215	CALL wrd_write_string( 'max_pr_user' )
216	WRITE ( 14 ) max_pr_user
217
218	CALL wrd_write_string( 'statistic_regions' )
219	WRITE ( 14 ) statistic_regions
220
221	!
222	!-- Caution: After changes in the following parameter-list, the version number stored in the
223	!-- -------- variable binary_version_global has to be increased. The same changes must also be
224	!-- done in the parameter-list in rrd_global.
225	CALL wrd_write_string( 'advected_distance_x' )
226	WRITE ( 14 ) advected_distance_x
227
228	CALL wrd_write_string( 'advected_distance_y' )
229	WRITE ( 14 ) advected_distance_y
230
231	CALL wrd_write_string( 'alpha_surface' )
232	WRITE ( 14 ) alpha_surface
233
234	CALL wrd_write_string( 'average_count_pr' )
235	WRITE ( 14 ) average_count_pr
236
237	CALL wrd_write_string( 'average_count_sp' )
238	WRITE ( 14 ) average_count_sp
239
240	CALL wrd_write_string( 'average_count_3d' )
241	WRITE ( 14 ) average_count_3d
242
243	CALL wrd_write_string( 'bc_e_b' )
244	WRITE ( 14 ) bc_e_b
245
246	CALL wrd_write_string( 'bc_lr' )
247	WRITE ( 14 ) bc_lr
248
249	CALL wrd_write_string( 'bc_ns' )
250	WRITE ( 14 ) bc_ns
251
252	CALL wrd_write_string( 'bc_p_b' )
253	WRITE ( 14 ) bc_p_b
254
255	CALL wrd_write_string( 'bc_p_t' )
256	WRITE ( 14 ) bc_p_t
257
258	CALL wrd_write_string( 'bc_pt_b' )
259	WRITE ( 14 ) bc_pt_b
260
261	CALL wrd_write_string( 'bc_pt_t' )
262	WRITE ( 14 ) bc_pt_t
263
264	CALL wrd_write_string( 'bc_pt_t_val' )
265	WRITE ( 14 ) bc_pt_t_val
266
267	CALL wrd_write_string( 'bc_q_b' )
268	WRITE ( 14 ) bc_q_b
269
270	CALL wrd_write_string( 'bc_q_t' )
271	WRITE ( 14 ) bc_q_t
272
273	CALL wrd_write_string( 'bc_q_t_val' )
274	WRITE ( 14 ) bc_q_t_val
275
276	CALL wrd_write_string( 'bc_s_b' )
277	WRITE ( 14 ) bc_s_b
278
279	CALL wrd_write_string( 'bc_s_t' )
280	WRITE ( 14 ) bc_s_t
281
282	CALL wrd_write_string( 'bc_uv_b' )
283	WRITE ( 14 ) bc_uv_b
284
285	CALL wrd_write_string( 'bc_uv_t' )
286	WRITE ( 14 ) bc_uv_t
287
288	CALL wrd_write_string( 'building_height' )
289	WRITE ( 14 ) building_height
290
291	CALL wrd_write_string( 'building_length_x' )
292	WRITE ( 14 ) building_length_x
293
294	CALL wrd_write_string( 'building_length_y' )
295	WRITE ( 14 ) building_length_y
296
297	CALL wrd_write_string( 'building_wall_left' )
298	WRITE ( 14 ) building_wall_left
299
300	CALL wrd_write_string( 'building_wall_south' )
301	WRITE ( 14 ) building_wall_south
302
303	CALL wrd_write_string( 'call_psolver_at_all_substeps' )
304	WRITE ( 14 ) call_psolver_at_all_substeps
305
306	CALL wrd_write_string( 'canyon_height' )
307	WRITE ( 14 ) canyon_height
308
309	CALL wrd_write_string( 'canyon_wall_left' )
310	WRITE ( 14 ) canyon_wall_left
311
312	CALL wrd_write_string( 'canyon_wall_south' )
313	WRITE ( 14 ) canyon_wall_south
314
315	CALL wrd_write_string( 'canyon_width_x' )
316	WRITE ( 14 ) canyon_width_x
317
318	CALL wrd_write_string( 'canyon_width_y' )
319	WRITE ( 14 ) canyon_width_y
320
321	CALL wrd_write_string( 'cfl_factor' )
322	WRITE ( 14 ) cfl_factor
323
324	CALL wrd_write_string( 'cloud_droplets' )
325	WRITE ( 14 ) cloud_droplets
326
327	CALL wrd_write_string( 'collective_wait' )
328	WRITE ( 14 ) collective_wait
329
330	CALL wrd_write_string( 'conserve_volume_flow' )
331	WRITE ( 14 ) conserve_volume_flow
332
333	CALL wrd_write_string( 'conserve_volume_flow_mode' )
334	WRITE ( 14 ) conserve_volume_flow_mode
335
336	CALL wrd_write_string( 'constant_flux_layer' )
337	WRITE ( 14 ) constant_flux_layer
338
339	CALL wrd_write_string( 'coupling_start_time' )
340	WRITE ( 14 ) coupling_start_time
341
342	CALL wrd_write_string( 'current_timestep_number' )
343	WRITE ( 14 ) current_timestep_number
344
345	CALL wrd_write_string( 'cycle_mg' )
346	WRITE ( 14 ) cycle_mg
347
348	CALL wrd_write_string( 'damp_level_1d' )
349	WRITE ( 14 ) damp_level_1d
350
351	CALL wrd_write_string( 'origin_date_time' )
352	WRITE ( 14 ) origin_date_time
353
354	CALL wrd_write_string( 'dissipation_1d' )
355	WRITE ( 14 ) dissipation_1d
356
357	CALL wrd_write_string( 'do2d_xy_time_count' )
358	WRITE ( 14 ) do2d_xy_time_count
359
360	CALL wrd_write_string( 'do2d_xz_time_count' )
361	WRITE ( 14 ) do2d_xz_time_count
362
363	CALL wrd_write_string( 'do2d_yz_time_count' )
364	WRITE ( 14 ) do2d_yz_time_count
365
366	CALL wrd_write_string( 'do3d_time_count' )
367	WRITE ( 14 ) do3d_time_count
368
369	CALL wrd_write_string( 'dp_external' )
370	WRITE ( 14 ) dp_external
371
372	CALL wrd_write_string( 'dp_level_b' )
373	WRITE ( 14 ) dp_level_b
374
375	CALL wrd_write_string( 'dp_smooth' )
376	WRITE ( 14 ) dp_smooth
377
378	CALL wrd_write_string( 'dpdxy' )
379	WRITE ( 14 ) dpdxy
380
381	CALL wrd_write_string( 'dt_3d' )
382	WRITE ( 14 ) dt_3d
383
384	CALL wrd_write_string( 'dt_pr_1d' )
385	WRITE ( 14 ) dt_pr_1d
386
387	CALL wrd_write_string( 'dt_run_control_1d' )
388	WRITE ( 14 ) dt_run_control_1d
389
390	CALL wrd_write_string( 'dx' )
391	WRITE ( 14 ) dx
392
393	CALL wrd_write_string( 'dy' )
394	WRITE ( 14 ) dy
395
396	CALL wrd_write_string( 'dz' )
397	WRITE ( 14 ) dz
398
399	CALL wrd_write_string( 'dz_max' )
400	WRITE ( 14 ) dz_max
401
402	CALL wrd_write_string( 'dz_stretch_factor' )
403	WRITE ( 14 ) dz_stretch_factor
404
405	CALL wrd_write_string( 'dz_stretch_factor_array' )
406	WRITE ( 14 ) dz_stretch_factor_array
407
408	CALL wrd_write_string( 'dz_stretch_level' )
409	WRITE ( 14 ) dz_stretch_level
410
411	CALL wrd_write_string( 'dz_stretch_level_end' )
412	WRITE ( 14 ) dz_stretch_level_end
413
414	CALL wrd_write_string( 'dz_stretch_level_start' )
415	WRITE ( 14 ) dz_stretch_level_start
416
417	CALL wrd_write_string( 'e_min' )
418	WRITE ( 14 ) e_min
419
420	CALL wrd_write_string( 'end_time_1d' )
421	WRITE ( 14 ) end_time_1d
422
423	CALL wrd_write_string( 'fft_method' )
424	WRITE ( 14 ) fft_method
425
426	CALL wrd_write_string( 'first_call_lpm' )
427	WRITE ( 14 ) first_call_lpm
428
429	CALL wrd_write_string( 'galilei_transformation' )
430	WRITE ( 14 ) galilei_transformation
431
432	CALL wrd_write_string( 'hom' )
433	WRITE ( 14 ) hom
434
435	CALL wrd_write_string( 'hom_sum' )
436	WRITE ( 14 ) hom_sum
437
438	CALL wrd_write_string( 'humidity' )
439	WRITE ( 14 ) humidity
440
441	IF ( ALLOCATED( inflow_damping_factor ) ) THEN
442	CALL wrd_write_string( 'inflow_damping_factor' )
443	WRITE ( 14 ) inflow_damping_factor
444	ENDIF
445
446	CALL wrd_write_string( 'inflow_damping_height' )
447	WRITE ( 14 ) inflow_damping_height
448
449	CALL wrd_write_string( 'inflow_damping_width' )
450	WRITE ( 14 ) inflow_damping_width
451
452	CALL wrd_write_string( 'inflow_disturbance_begin' )
453	WRITE ( 14 ) inflow_disturbance_begin
454
455	CALL wrd_write_string( 'inflow_disturbance_end' )
456	WRITE ( 14 ) inflow_disturbance_end
457
458	CALL wrd_write_string( 'km_constant' )
459	WRITE ( 14 ) km_constant
460
461	CALL wrd_write_string( 'large_scale_forcing' )
462	WRITE ( 14 ) large_scale_forcing
463
464	CALL wrd_write_string( 'large_scale_subsidence' )
465	WRITE ( 14 ) large_scale_subsidence
466
467	CALL wrd_write_string( 'latitude' )
468	WRITE ( 14 ) latitude
469
470	CALL wrd_write_string( 'longitude' )
471	WRITE ( 14 ) longitude
472
473	CALL wrd_write_string( 'loop_optimization' )
474	WRITE ( 14 ) loop_optimization
475
476	CALL wrd_write_string( 'masking_method' )
477	WRITE ( 14 ) masking_method
478
479	IF ( ALLOCATED( mean_inflow_profiles ) ) THEN
480	CALL wrd_write_string( 'mean_inflow_profiles' )
481	WRITE ( 14 ) mean_inflow_profiles
482	ENDIF
483
484	CALL wrd_write_string( 'mg_cycles' )
485	WRITE ( 14 ) mg_cycles
486
487	CALL wrd_write_string( 'mg_switch_to_pe0_level' )
488	WRITE ( 14 ) mg_switch_to_pe0_level
489
490	CALL wrd_write_string( 'mixing_length_1d' )
491	WRITE ( 14 ) mixing_length_1d
492
493	CALL wrd_write_string( 'momentum_advec' )
494	WRITE ( 14 ) momentum_advec
495
496	CALL wrd_write_string( 'netcdf_precision' )
497	WRITE ( 14 ) netcdf_precision
498
499	CALL wrd_write_string( 'neutral' )
500	WRITE ( 14 ) neutral
501
502	CALL wrd_write_string( 'ngsrb' )
503	WRITE ( 14 ) ngsrb
504
505	CALL wrd_write_string( 'nsor' )
506	WRITE ( 14 ) nsor
507
508	CALL wrd_write_string( 'nsor_ini' )
509	WRITE ( 14 ) nsor_ini
510
511	CALL wrd_write_string( 'nudging' )
512	WRITE ( 14 ) nudging
513
514	CALL wrd_write_string( 'num_leg' )
515	WRITE ( 14 ) num_leg
516
517	CALL wrd_write_string( 'nx' )
518	WRITE ( 14 ) nx
519
520	CALL wrd_write_string( 'ny' )
521	WRITE ( 14 ) ny
522
523	CALL wrd_write_string( 'ocean_mode' )
524	WRITE ( 14 ) ocean_mode
525
526	CALL wrd_write_string( 'omega' )
527	WRITE ( 14 ) omega
528
529	CALL wrd_write_string( 'omega_sor' )
530	WRITE ( 14 ) omega_sor
531
532	CALL wrd_write_string( 'output_for_t0' )
533	WRITE ( 14 ) output_for_t0
534
535	CALL wrd_write_string( 'passive_scalar' )
536	WRITE ( 14 ) passive_scalar
537
538	CALL wrd_write_string( 'prandtl_number' )
539	WRITE ( 14 ) prandtl_number
540
541	CALL wrd_write_string( 'psolver' )
542	WRITE ( 14 ) psolver
543
544	CALL wrd_write_string( 'pt_damping_factor' )
545	WRITE ( 14 ) pt_damping_factor
546
547	CALL wrd_write_string( 'pt_damping_width' )
548	WRITE ( 14 ) pt_damping_width
549
550	CALL wrd_write_string( 'pt_init' )
551	WRITE ( 14 ) pt_init
552
553	CALL wrd_write_string( 'pt_reference' )
554	WRITE ( 14 ) pt_reference
555
556	CALL wrd_write_string( 'pt_surface' )
557	WRITE ( 14 ) pt_surface
558
559	CALL wrd_write_string( 'pt_surface_initial_change' )
560	WRITE ( 14 ) pt_surface_initial_change
561
562	CALL wrd_write_string( 'pt_vertical_gradient' )
563	WRITE ( 14 ) pt_vertical_gradient
564
565	CALL wrd_write_string( 'pt_vertical_gradient_level' )
566	WRITE ( 14 ) pt_vertical_gradient_level
567
568	CALL wrd_write_string( 'pt_vertical_gradient_level_ind' )
569	WRITE ( 14 ) pt_vertical_gradient_level_ind
570
571	CALL wrd_write_string( 'q_init' )
572	WRITE ( 14 ) q_init
573
574	CALL wrd_write_string( 'q_surface' )
575	WRITE ( 14 ) q_surface
576
577	CALL wrd_write_string( 'q_surface_initial_change' )
578	WRITE ( 14 ) q_surface_initial_change
579
580	CALL wrd_write_string( 'q_vertical_gradient' )
581	WRITE ( 14 ) q_vertical_gradient
582
583	CALL wrd_write_string( 'q_vertical_gradient_level' )
584	WRITE ( 14 ) q_vertical_gradient_level
585
586	CALL wrd_write_string( 'q_vertical_gradient_level_ind' )
587	WRITE ( 14 ) q_vertical_gradient_level_ind
588
589	CALL wrd_write_string( 'random_generator' )
590	WRITE ( 14 ) random_generator
591
592	CALL wrd_write_string( 'random_heatflux' )
593	WRITE ( 14 ) random_heatflux
594
595	CALL wrd_write_string( 'rans_mode' )
596	WRITE ( 14 ) rans_mode
597
598	CALL wrd_write_string( 'rayleigh_damping_factor' )
599	WRITE ( 14 ) rayleigh_damping_factor
600
601	CALL wrd_write_string( 'rayleigh_damping_height' )
602	WRITE ( 14 ) rayleigh_damping_height
603
604	CALL wrd_write_string( 'recycling_width' )
605	WRITE ( 14 ) recycling_width
606
607	CALL wrd_write_string( 'ref_state' )
608	WRITE ( 14 ) ref_state
609
610	CALL wrd_write_string( 'reference_state' )
611	WRITE ( 14 ) reference_state
612
613	CALL wrd_write_string( 'residual_limit' )
614	WRITE ( 14 ) residual_limit
615
616	CALL wrd_write_string( 'roughness_length' )
617	WRITE ( 14 ) roughness_length
618
619	CALL wrd_write_string( 'run_coupled' )
620	WRITE ( 14 ) run_coupled
621
622	CALL wrd_write_string( 'runnr' )
623	WRITE ( 14 ) runnr
624
625	CALL wrd_write_string( 's_init' )
626	WRITE ( 14 ) s_init
627
628	CALL wrd_write_string( 's_surface' )
629	WRITE ( 14 ) s_surface
630
631	CALL wrd_write_string( 's_surface_initial_change' )
632	WRITE ( 14 ) s_surface_initial_change
633
634	CALL wrd_write_string( 's_vertical_gradient' )
635	WRITE ( 14 ) s_vertical_gradient
636
637	CALL wrd_write_string( 's_vertical_gradient_level' )
638	WRITE ( 14 ) s_vertical_gradient_level
639
640	CALL wrd_write_string( 's_vertical_gradient_level_ind' )
641	WRITE ( 14 ) s_vertical_gradient_level_ind
642
643	CALL wrd_write_string( 'scalar_advec' )
644	WRITE ( 14 ) scalar_advec
645
646	CALL wrd_write_string( 'simulated_time' )
647	WRITE ( 14 ) simulated_time
648
649	IF ( ALLOCATED( spectrum_x ) ) THEN
650	CALL wrd_write_string( 'spectrum_x' )
651	WRITE ( 14 ) spectrum_x
652	CALL wrd_write_string( 'spectrum_y' )
653	WRITE ( 14 ) spectrum_y
654	ENDIF
655
656	CALL wrd_write_string( 'spinup_time ' )
657	WRITE ( 14 ) spinup_time
658
659	CALL wrd_write_string( 'subs_vertical_gradient' )
660	WRITE ( 14 ) subs_vertical_gradient
661
662	CALL wrd_write_string( 'subs_vertical_gradient_level' )
663	WRITE ( 14 ) subs_vertical_gradient_level
664
665	CALL wrd_write_string( 'subs_vertical_gradient_level_i' )
666	WRITE ( 14 ) subs_vertical_gradient_level_i
667
668	CALL wrd_write_string( 'surface_heatflux' )
669	WRITE ( 14 ) surface_heatflux
670
671	CALL wrd_write_string( 'surface_pressure' )
672	WRITE ( 14 ) surface_pressure
673
674	CALL wrd_write_string( 'surface_scalarflux' )
675	WRITE ( 14 ) surface_scalarflux
676
677	CALL wrd_write_string( 'surface_waterflux' )
678	WRITE ( 14 ) surface_waterflux
679
680	CALL wrd_write_string( 'time_coupling' )
681	WRITE ( 14 ) time_coupling
682
683	CALL wrd_write_string( 'time_disturb' )
684	WRITE ( 14 ) time_disturb
685
686	CALL wrd_write_string( 'time_do2d_xy' )
687	WRITE ( 14 ) time_do2d_xy
688
689	CALL wrd_write_string( 'time_do2d_xz' )
690	WRITE ( 14 ) time_do2d_xz
691
692	CALL wrd_write_string( 'time_do2d_yz' )
693	WRITE ( 14 ) time_do2d_yz
694
695	CALL wrd_write_string( 'time_do3d' )
696	WRITE ( 14 ) time_do3d
697
698	CALL wrd_write_string( 'time_do_av' )
699	WRITE ( 14 ) time_do_av
700
701	CALL wrd_write_string( 'time_do_sla' )
702	WRITE ( 14 ) time_do_sla
703
704	CALL wrd_write_string( 'time_domask' )
705	WRITE ( 14 ) time_domask
706
707	CALL wrd_write_string( 'time_dopr' )
708	WRITE ( 14 ) time_dopr
709
710	CALL wrd_write_string( 'time_dopr_av' )
711	WRITE ( 14 ) time_dopr_av
712
713	CALL wrd_write_string( 'time_dopr_listing' )
714	WRITE ( 14 ) time_dopr_listing
715
716	CALL wrd_write_string( 'time_dopts' )
717	WRITE ( 14 ) time_dopts
718
719	CALL wrd_write_string( 'time_dosp' )
720	WRITE ( 14 ) time_dosp
721
722	CALL wrd_write_string( 'time_dots' )
723	WRITE ( 14 ) time_dots
724
725	CALL wrd_write_string( 'time_indoor' )
726	WRITE ( 14 ) time_indoor
727
728	CALL wrd_write_string( 'time_radiation' )
729	WRITE ( 14 ) time_radiation
730
731	CALL wrd_write_string( 'time_restart' )
732	WRITE ( 14 ) time_restart
733
734	CALL wrd_write_string( 'time_run_control' )
735	WRITE ( 14 ) time_run_control
736
737	CALL wrd_write_string( 'time_since_reference_point' )
738	WRITE ( 14 ) time_since_reference_point
739
740	CALL wrd_write_string( 'time_virtual_measurement' )
741	WRITE ( 14 ) time_virtual_measurement
742
743	CALL wrd_write_string( 'timestep_scheme' )
744	WRITE ( 14 ) timestep_scheme
745
746	CALL wrd_write_string( 'top_heatflux' )
747	WRITE ( 14 ) top_heatflux
748
749	CALL wrd_write_string( 'top_momentumflux_u' )
750	WRITE ( 14 ) top_momentumflux_u
751
752	CALL wrd_write_string( 'top_momentumflux_v' )
753	WRITE ( 14 ) top_momentumflux_v
754
755	CALL wrd_write_string( 'top_scalarflux' )
756	WRITE ( 14 ) top_scalarflux
757
758	CALL wrd_write_string( 'topography' )
759	WRITE ( 14 ) topography
760
761	CALL wrd_write_string( 'topography_grid_convention' )
762	WRITE ( 14 ) topography_grid_convention
763
764	CALL wrd_write_string( 'tsc' )
765	WRITE ( 14 ) tsc
766
767	CALL wrd_write_string( 'tunnel_height' )
768	WRITE ( 14 ) tunnel_height
769
770	CALL wrd_write_string( 'tunnel_length' )
771	WRITE ( 14 ) tunnel_length
772
773	CALL wrd_write_string( 'tunnel_wall_depth' )
774	WRITE ( 14 ) tunnel_wall_depth
775
776	CALL wrd_write_string( 'tunnel_width_x' )
777	WRITE ( 14 ) tunnel_width_x
778
779	CALL wrd_write_string( 'tunnel_width_y' )
780	WRITE ( 14 ) tunnel_width_y
781
782	CALL wrd_write_string( 'turbulence_closure' )
783	WRITE ( 14 ) turbulence_closure
784
785	CALL wrd_write_string( 'turbulent_inflow' )
786	WRITE ( 14 ) turbulent_inflow
787
788	CALL wrd_write_string( 'u_bulk' )
789	WRITE ( 14 ) u_bulk
790
791	CALL wrd_write_string( 'u_init' )
792	WRITE ( 14 ) u_init
793
794	CALL wrd_write_string( 'u_max' )
795	WRITE ( 14 ) u_max
796
797	CALL wrd_write_string( 'u_max_ijk' )
798	WRITE ( 14 ) u_max_ijk
799
800	CALL wrd_write_string( 'ug' )
801	WRITE ( 14 ) ug
802
803	CALL wrd_write_string( 'ug_surface' )
804	WRITE ( 14 ) ug_surface
805
806	CALL wrd_write_string( 'ug_vertical_gradient' )
807	WRITE ( 14 ) ug_vertical_gradient
808
809	CALL wrd_write_string( 'ug_vertical_gradient_level' )
810	WRITE ( 14 ) ug_vertical_gradient_level
811
812	CALL wrd_write_string( 'ug_vertical_gradient_level_ind' )
813	WRITE ( 14 ) ug_vertical_gradient_level_ind
814
815	CALL wrd_write_string( 'use_surface_fluxes' )
816	WRITE ( 14 ) use_surface_fluxes
817
818	CALL wrd_write_string( 'use_top_fluxes' )
819	WRITE ( 14 ) use_top_fluxes
820
821	CALL wrd_write_string( 'use_ug_for_galilei_tr' )
822	WRITE ( 14 ) use_ug_for_galilei_tr
823
824	CALL wrd_write_string( 'use_upstream_for_tke' )
825	WRITE ( 14 ) use_upstream_for_tke
826
827	CALL wrd_write_string( 'v_bulk' )
828	WRITE ( 14 ) v_bulk
829
830	CALL wrd_write_string( 'v_init' )
831	WRITE ( 14 ) v_init
832
833	CALL wrd_write_string( 'v_max' )
834	WRITE ( 14 ) v_max
835
836	CALL wrd_write_string( 'v_max_ijk' )
837	WRITE ( 14 ) v_max_ijk
838
839	CALL wrd_write_string( 'vg' )
840	WRITE ( 14 ) vg
841
842	CALL wrd_write_string( 'vg_surface' )
843	WRITE ( 14 ) vg_surface
844
845	CALL wrd_write_string( 'vg_vertical_gradient' )
846	WRITE ( 14 ) vg_vertical_gradient
847
848	CALL wrd_write_string( 'vg_vertical_gradient_level' )
849	WRITE ( 14 ) vg_vertical_gradient_level
850
851	CALL wrd_write_string( 'vg_vertical_gradient_level_ind' )
852	WRITE ( 14 ) vg_vertical_gradient_level_ind
853
854	CALL wrd_write_string( 'virtual_flight' )
855	WRITE ( 14 ) virtual_flight
856
857	CALL wrd_write_string( 'volume_flow_area' )
858	WRITE ( 14 ) volume_flow_area
859
860	CALL wrd_write_string( 'volume_flow_initial' )
861	WRITE ( 14 ) volume_flow_initial
862
863	CALL wrd_write_string( 'w_max' )
864	WRITE ( 14 ) w_max
865
866	CALL wrd_write_string( 'w_max_ijk' )
867	WRITE ( 14 ) w_max_ijk
868
869	CALL wrd_write_string( 'wall_adjustment' )
870	WRITE ( 14 ) wall_adjustment
871
872	CALL wrd_write_string( 'wall_heatflux' )
873	WRITE ( 14 ) wall_heatflux
874
875	CALL wrd_write_string( 'wall_humidityflux' )
876	WRITE ( 14 ) wall_humidityflux
877
878	CALL wrd_write_string( 'wall_scalarflux' )
879	WRITE ( 14 ) wall_scalarflux
880
881	CALL wrd_write_string( 'y_shift' )
882	WRITE ( 14 ) y_shift
883
884	CALL wrd_write_string( 'z0h_factor' )
885	WRITE ( 14 ) z0h_factor
886
887	CALL wrd_write_string( 'zeta_max' )
888	WRITE ( 14 ) zeta_max
889
890	CALL wrd_write_string( 'zeta_min' )
891	WRITE ( 14 ) zeta_min
892
893	CALL wrd_write_string( 'z_i' )
894	WRITE ( 14 ) z_i
895
896	ELSEIF ( restart_data_format_output(1:3) == 'mpi' ) THEN
897	!
898	!-- Write global restart data using MPI-IO
899	!-- ATTENTION: Arrays need to be written with routine wrd_mpi_io_global_array!
900	CALL wrd_mpi_io( 'binary_version_global', binary_version_global )
901	CALL wrd_mpi_io( 'numprocs', numprocs )
902	CALL wrd_mpi_io( 'nz' , nz )
903	CALL wrd_mpi_io( 'max_pr_user', max_pr_user )
904	CALL wrd_mpi_io( 'statistic_regions', statistic_regions )
905
906	!
907	!-- Caution: After changes in the following parameter-list, the version number stored in the
908	!-- -------- variable binary_version_global has to be increased. The same changes must also be
909	!-- done in the parameter-list in rrd_global.
910	CALL wrd_mpi_io( 'advected_distance_x', advected_distance_x )
911	CALL wrd_mpi_io( 'advected_distance_y', advected_distance_y )
912	CALL wrd_mpi_io( 'alpha_surface', alpha_surface )
913	CALL wrd_mpi_io( 'average_count_pr', average_count_pr )
914	CALL wrd_mpi_io( 'average_count_sp', average_count_sp )
915	CALL wrd_mpi_io( 'average_count_3d', average_count_3d )
916	CALL wrd_mpi_io( 'bc_e_b', bc_e_b )
917	CALL wrd_mpi_io( 'bc_lr', bc_lr )
918	CALL wrd_mpi_io( 'bc_ns', bc_ns )
919	CALL wrd_mpi_io( 'bc_p_b', bc_p_b )
920	CALL wrd_mpi_io( 'bc_p_t', bc_p_t )
921	CALL wrd_mpi_io( 'bc_pt_b', bc_pt_b )
922	CALL wrd_mpi_io( 'bc_pt_t', bc_pt_t )
923	CALL wrd_mpi_io( 'bc_pt_t_val', bc_pt_t_val )
924	CALL wrd_mpi_io( 'bc_q_b', bc_q_b )
925	CALL wrd_mpi_io( 'bc_q_t', bc_q_t )
926	CALL wrd_mpi_io( 'bc_q_t_val', bc_q_t_val )
927	CALL wrd_mpi_io( 'bc_s_b', bc_s_b )
928	CALL wrd_mpi_io( 'bc_s_t', bc_s_t )
929	CALL wrd_mpi_io( 'bc_uv_b', bc_uv_b )
930	CALL wrd_mpi_io( 'bc_uv_t', bc_uv_t )
931	CALL wrd_mpi_io( 'biometeorology', biometeorology )
932	CALL wrd_mpi_io( 'building_height', building_height )
933	CALL wrd_mpi_io( 'building_length_x', building_length_x )
934	CALL wrd_mpi_io( 'building_length_y', building_length_y )
935	CALL wrd_mpi_io( 'building_wall_left', building_wall_left )
936	CALL wrd_mpi_io( 'building_wall_south', building_wall_south )
937	CALL wrd_mpi_io( 'bulk_cloud_model', bulk_cloud_model )
938	CALL wrd_mpi_io( 'call_psolver_at_all_substeps', call_psolver_at_all_substeps )
939	CALL wrd_mpi_io( 'canyon_height', canyon_height )
940	CALL wrd_mpi_io( 'canyon_wall_left', canyon_wall_left )
941	CALL wrd_mpi_io( 'canyon_wall_south', canyon_wall_south )
942	CALL wrd_mpi_io( 'canyon_width_x', canyon_width_x )
943	CALL wrd_mpi_io( 'canyon_width_y', canyon_width_y )
944	CALL wrd_mpi_io( 'cfl_factor', cfl_factor )
945	CALL wrd_mpi_io( 'cloud_droplets', cloud_droplets )
946	CALL wrd_mpi_io( 'collective_wait', collective_wait )
947	CALL wrd_mpi_io( 'conserve_volume_flow', conserve_volume_flow )
948	CALL wrd_mpi_io( 'conserve_volume_flow_mode', conserve_volume_flow_mode )
949	CALL wrd_mpi_io( 'constant_flux_layer', constant_flux_layer )
950	CALL wrd_mpi_io( 'coupling_start_time', coupling_start_time )
951	CALL wrd_mpi_io( 'current_timestep_number', current_timestep_number )
952	CALL wrd_mpi_io( 'cycle_mg', cycle_mg )
953	CALL wrd_mpi_io( 'damp_level_1d', damp_level_1d )
954	CALL wrd_mpi_io( 'dissipation_1d', dissipation_1d )
955	CALL wrd_mpi_io_global_array( 'do2d_xy_time_count', do2d_xy_time_count )
956	CALL wrd_mpi_io_global_array( 'do2d_xz_time_count', do2d_xz_time_count )
957	CALL wrd_mpi_io_global_array( 'do2d_yz_time_count', do2d_yz_time_count )
958	CALL wrd_mpi_io_global_array( 'do3d_time_count', do3d_time_count )
959	CALL wrd_mpi_io( 'dp_external', dp_external )
960	CALL wrd_mpi_io( 'dp_level_b', dp_level_b )
961	CALL wrd_mpi_io( 'dp_smooth', dp_smooth )
962	CALL wrd_mpi_io_global_array( 'dpdxy', dpdxy )
963	CALL wrd_mpi_io( 'dt_3d', dt_3d )
964	CALL wrd_mpi_io( 'dt_pr_1d', dt_pr_1d )
965	CALL wrd_mpi_io( 'dt_run_control_1d', dt_run_control_1d )
966	CALL wrd_mpi_io( 'dx', dx )
967	CALL wrd_mpi_io( 'dy', dy )
968	CALL wrd_mpi_io_global_array( 'dz', dz )
969	CALL wrd_mpi_io( 'dz_max', dz_max )
970	CALL wrd_mpi_io( 'dz_stretch_factor', dz_stretch_factor )
971	CALL wrd_mpi_io_global_array( 'dz_stretch_factor_array', dz_stretch_factor_array )
972	CALL wrd_mpi_io( 'dz_stretch_level', dz_stretch_level )
973	CALL wrd_mpi_io_global_array( 'dz_stretch_level_end', dz_stretch_level_end )
974	CALL wrd_mpi_io_global_array( 'dz_stretch_level_start', dz_stretch_level_start )
975	CALL wrd_mpi_io( 'e_min', e_min )
976	CALL wrd_mpi_io( 'end_time_1d', end_time_1d )
977	CALL wrd_mpi_io( 'fft_method', fft_method )
978	CALL wrd_mpi_io( 'first_call_lpm', first_call_lpm )
979	CALL wrd_mpi_io( 'galilei_transformation', galilei_transformation )
980	CALL wrd_mpi_io( 'gust_module_enabled', gust_module_enabled )
981	CALL wrd_mpi_io_global_array( 'hom', hom )
982	CALL wrd_mpi_io_global_array( 'hom_sum', hom_sum )
983	CALL wrd_mpi_io( 'humidity', humidity )
984	IF ( ALLOCATED( inflow_damping_factor ) ) THEN
985	CALL wrd_mpi_io_global_array( 'inflow_damping_factor', inflow_damping_factor )
986	ENDIF
987	CALL wrd_mpi_io( 'inflow_damping_height', inflow_damping_height )
988	CALL wrd_mpi_io( 'inflow_damping_width', inflow_damping_width )
989	CALL wrd_mpi_io( 'inflow_disturbance_begin', inflow_disturbance_begin )
990	CALL wrd_mpi_io( 'inflow_disturbance_end', inflow_disturbance_end )
991	CALL wrd_mpi_io( 'km_constant', km_constant )
992	CALL wrd_mpi_io( 'large_scale_forcing', large_scale_forcing )
993	CALL wrd_mpi_io( 'large_scale_subsidence', large_scale_subsidence )
994	CALL wrd_mpi_io( 'latitude', latitude )
995	CALL wrd_mpi_io( 'longitude', longitude )
996	CALL wrd_mpi_io( 'loop_optimization', loop_optimization )
997	CALL wrd_mpi_io( 'masking_method', masking_method )
998	IF ( ALLOCATED( mean_inflow_profiles ) ) THEN
999	CALL wrd_mpi_io_global_array( 'mean_inflow_profiles', mean_inflow_profiles )
1000	ENDIF
1001	CALL wrd_mpi_io( 'mg_cycles', mg_cycles )
1002	CALL wrd_mpi_io( 'mg_switch_to_pe0_level', mg_switch_to_pe0_level )
1003	CALL wrd_mpi_io( 'mixing_length_1d', mixing_length_1d )
1004	CALL wrd_mpi_io( 'momentum_advec', momentum_advec )
1005	!
1006	!-- There is no module procedure for CHARACTER arrays
1007	DO i = 1, SIZE( netcdf_precision, 1 )
1008	WRITE ( tmp_name, '(A,I2.2)' ) 'netcdf_precision_', i
1009	CALL wrd_mpi_io( TRIM( tmp_name ), netcdf_precision(i) )
1010	ENDDO
1011	CALL wrd_mpi_io( 'neutral', neutral )
1012	CALL wrd_mpi_io( 'ngsrb', ngsrb )
1013	CALL wrd_mpi_io( 'nsor', nsor )
1014	CALL wrd_mpi_io( 'nsor_ini', nsor_ini )
1015	CALL wrd_mpi_io( 'nudging', nudging )
1016	CALL wrd_mpi_io( 'num_leg', num_leg )
1017	CALL wrd_mpi_io( 'nx', nx )
1018	CALL wrd_mpi_io( 'ny', ny )
1019	CALL wrd_mpi_io( 'ocean_mode', ocean_mode )
1020	CALL wrd_mpi_io( 'omega', omega )
1021	CALL wrd_mpi_io( 'omega_sor', omega_sor )
1022	CALL wrd_mpi_io( 'origin_date_time', origin_date_time )
1023	CALL wrd_mpi_io( 'output_for_t0', output_for_t0 )
1024	CALL wrd_mpi_io( 'particle_advection', particle_advection )
1025	CALL wrd_mpi_io( 'passive_scalar', passive_scalar )
1026	CALL wrd_mpi_io( 'prandtl_number', prandtl_number )
1027	CALL wrd_mpi_io( 'psolver', psolver )
1028	CALL wrd_mpi_io( 'pt_damping_factor', pt_damping_factor )
1029	CALL wrd_mpi_io( 'pt_damping_width', pt_damping_width )
1030	CALL wrd_mpi_io_global_array( 'pt_init', pt_init )
1031	CALL wrd_mpi_io( 'pt_reference', pt_reference )
1032	CALL wrd_mpi_io( 'pt_surface', pt_surface )
1033	CALL wrd_mpi_io( 'pt_surface_initial_change', pt_surface_initial_change )
1034	CALL wrd_mpi_io_global_array( 'pt_vertical_gradient', pt_vertical_gradient )
1035	CALL wrd_mpi_io_global_array( 'pt_vertical_gradient_level', pt_vertical_gradient_level )
1036	CALL wrd_mpi_io_global_array( 'pt_vertical_gradient_level_ind', pt_vertical_gradient_level_ind )
1037	CALL wrd_mpi_io_global_array( 'q_init', q_init )
1038	CALL wrd_mpi_io( 'q_surface', q_surface )
1039	CALL wrd_mpi_io( 'q_surface_initial_change', q_surface_initial_change )
1040	CALL wrd_mpi_io_global_array( 'q_vertical_gradient', q_vertical_gradient )
1041	CALL wrd_mpi_io_global_array( 'q_vertical_gradient_level', q_vertical_gradient_level )
1042	CALL wrd_mpi_io_global_array( 'q_vertical_gradient_level_ind', q_vertical_gradient_level_ind )
1043	CALL wrd_mpi_io( 'random_generator', random_generator )
1044	CALL wrd_mpi_io( 'random_heatflux', random_heatflux )
1045	CALL wrd_mpi_io( 'rans_mode', rans_mode )
1046	CALL wrd_mpi_io( 'rayleigh_damping_factor', rayleigh_damping_factor )
1047	CALL wrd_mpi_io( 'rayleigh_damping_height', rayleigh_damping_height )
1048	CALL wrd_mpi_io( 'recycling_width', recycling_width )
1049	CALL wrd_mpi_io_global_array( 'ref_state', ref_state )
1050	CALL wrd_mpi_io( 'reference_state', reference_state )
1051	CALL wrd_mpi_io( 'residual_limit', residual_limit )
1052	CALL wrd_mpi_io( 'roughness_length', roughness_length )
1053	CALL wrd_mpi_io( 'run_coupled', run_coupled )
1054	CALL wrd_mpi_io( 'runnr', runnr )
1055	CALL wrd_mpi_io_global_array( 's_init', s_init )
1056	CALL wrd_mpi_io( 's_surface', s_surface )
1057	CALL wrd_mpi_io( 's_surface_initial_change', s_surface_initial_change )
1058	CALL wrd_mpi_io_global_array( 's_vertical_gradient', s_vertical_gradient )
1059	CALL wrd_mpi_io_global_array( 's_vertical_gradient_level', s_vertical_gradient_level )
1060	CALL wrd_mpi_io_global_array( 's_vertical_gradient_level_ind', s_vertical_gradient_level_ind )
1061	CALL wrd_mpi_io( 'scalar_advec', scalar_advec )
1062	CALL wrd_mpi_io( 'simulated_time', simulated_time )
1063	IF ( ALLOCATED( spectrum_x ) ) THEN
1064	CALL wrd_mpi_io_global_array( 'spectrum_x', spectrum_x )
1065	CALL wrd_mpi_io_global_array( 'spectrum_y', spectrum_y )
1066	ENDIF
1067	CALL wrd_mpi_io( 'spinup_time ', spinup_time )
1068	CALL wrd_mpi_io_global_array( 'subs_vertical_gradient', subs_vertical_gradient )
1069	CALL wrd_mpi_io_global_array( 'subs_vertical_gradient_level', subs_vertical_gradient_level )
1070	CALL wrd_mpi_io_global_array( 'subs_vertical_gradient_level_i', subs_vertical_gradient_level_i )
1071	CALL wrd_mpi_io( 'surface_heatflux', surface_heatflux )
1072	CALL wrd_mpi_io( 'surface_output', surface_output )
1073	CALL wrd_mpi_io( 'surface_pressure', surface_pressure )
1074	CALL wrd_mpi_io( 'surface_scalarflux', surface_scalarflux )
1075	CALL wrd_mpi_io( 'surface_waterflux', surface_waterflux )
1076	CALL wrd_mpi_io( 'syn_turb_gen', syn_turb_gen )
1077	CALL wrd_mpi_io( 'time_coupling', time_coupling )
1078	CALL wrd_mpi_io( 'time_disturb', time_disturb )
1079	CALL wrd_mpi_io( 'time_do2d_xy', time_do2d_xy )
1080	CALL wrd_mpi_io( 'time_do2d_xz', time_do2d_xz )
1081	CALL wrd_mpi_io( 'time_do2d_yz', time_do2d_yz )
1082	CALL wrd_mpi_io( 'time_do3d', time_do3d )
1083	CALL wrd_mpi_io( 'time_do_av', time_do_av )
1084	CALL wrd_mpi_io( 'time_do_sla', time_do_sla )
1085	CALL wrd_mpi_io_global_array( 'time_domask', time_domask )
1086	CALL wrd_mpi_io( 'time_dopr', time_dopr )
1087	CALL wrd_mpi_io( 'time_dopr_av', time_dopr_av )
1088	CALL wrd_mpi_io( 'time_dopr_listing', time_dopr_listing )
1089	CALL wrd_mpi_io( 'time_dopts', time_dopts )
1090	CALL wrd_mpi_io( 'time_dosp', time_dosp )
1091	CALL wrd_mpi_io( 'time_dots', time_dots )
1092	CALL wrd_mpi_io( 'time_indoor', time_indoor )
1093	CALL wrd_mpi_io( 'time_radiation', time_radiation )
1094	CALL wrd_mpi_io( 'time_restart', time_restart )
1095	CALL wrd_mpi_io( 'time_run_control', time_run_control )
1096	CALL wrd_mpi_io( 'time_since_reference_point', time_since_reference_point )
1097	CALL wrd_mpi_io( 'time_virtual_measurement', time_virtual_measurement )
1098	CALL wrd_mpi_io( 'timestep_scheme', timestep_scheme )
1099	CALL wrd_mpi_io( 'top_heatflux', top_heatflux )
1100	CALL wrd_mpi_io( 'top_momentumflux_u', top_momentumflux_u )
1101	CALL wrd_mpi_io( 'top_momentumflux_v', top_momentumflux_v )
1102	CALL wrd_mpi_io( 'top_scalarflux', top_scalarflux )
1103	CALL wrd_mpi_io( 'topography', topography )
1104	CALL wrd_mpi_io( 'topography_grid_convention', topography_grid_convention )
1105	CALL wrd_mpi_io_global_array( 'tsc', tsc )
1106	CALL wrd_mpi_io( 'tunnel_height', tunnel_height )
1107	CALL wrd_mpi_io( 'tunnel_length', tunnel_length )
1108	CALL wrd_mpi_io( 'tunnel_wall_depth', tunnel_wall_depth )
1109	CALL wrd_mpi_io( 'tunnel_width_x', tunnel_width_x )
1110	CALL wrd_mpi_io( 'tunnel_width_y', tunnel_width_y )
1111	CALL wrd_mpi_io( 'turbulence_closure', turbulence_closure )
1112	CALL wrd_mpi_io( 'turbulent_inflow', turbulent_inflow )
1113	CALL wrd_mpi_io( 'u_bulk', u_bulk )
1114	CALL wrd_mpi_io_global_array( 'u_init', u_init )
1115	CALL wrd_mpi_io( 'u_max', u_max )
1116	CALL wrd_mpi_io_global_array( 'u_max_ijk', u_max_ijk )
1117	CALL wrd_mpi_io_global_array( 'ug', ug )
1118	CALL wrd_mpi_io( 'ug_surface', ug_surface )
1119	CALL wrd_mpi_io_global_array( 'ug_vertical_gradient', ug_vertical_gradient )
1120	CALL wrd_mpi_io_global_array( 'ug_vertical_gradient_level', ug_vertical_gradient_level )
1121	CALL wrd_mpi_io_global_array( 'ug_vertical_gradient_level_ind', ug_vertical_gradient_level_ind )
1122	CALL wrd_mpi_io( 'use_surface_fluxes', use_surface_fluxes )
1123	CALL wrd_mpi_io( 'use_top_fluxes', use_top_fluxes )
1124	CALL wrd_mpi_io( 'use_ug_for_galilei_tr', use_ug_for_galilei_tr )
1125	CALL wrd_mpi_io( 'use_upstream_for_tke', use_upstream_for_tke )
1126	CALL wrd_mpi_io( 'user_module_enabled', user_module_enabled )
1127	CALL wrd_mpi_io( 'v_bulk', v_bulk )
1128	CALL wrd_mpi_io_global_array( 'v_init', v_init )
1129	CALL wrd_mpi_io( 'v_max', v_max )
1130	CALL wrd_mpi_io_global_array( 'v_max_ijk', v_max_ijk )
1131	CALL wrd_mpi_io_global_array( 'vg', vg )
1132	CALL wrd_mpi_io( 'vg_surface', vg_surface )
1133	CALL wrd_mpi_io_global_array( 'vg_vertical_gradient', vg_vertical_gradient )
1134	CALL wrd_mpi_io_global_array( 'vg_vertical_gradient_level', vg_vertical_gradient_level )
1135	CALL wrd_mpi_io_global_array( 'vg_vertical_gradient_level_ind', vg_vertical_gradient_level_ind )
1136	CALL wrd_mpi_io( 'virtual_flight', virtual_flight )
1137	CALL wrd_mpi_io_global_array( 'volume_flow_area', volume_flow_area )
1138	CALL wrd_mpi_io_global_array( 'volume_flow_initial', volume_flow_initial )
1139	CALL wrd_mpi_io( 'w_max', w_max )
1140	CALL wrd_mpi_io_global_array( 'w_max_ijk', w_max_ijk )
1141	CALL wrd_mpi_io( 'wall_adjustment', wall_adjustment )
1142	CALL wrd_mpi_io_global_array( 'wall_heatflux', wall_heatflux )
1143	CALL wrd_mpi_io_global_array( 'wall_humidityflux', wall_humidityflux )
1144	CALL wrd_mpi_io_global_array( 'wall_scalarflux', wall_scalarflux )
1145	CALL wrd_mpi_io( 'y_shift', y_shift )
1146	CALL wrd_mpi_io( 'z0h_factor', z0h_factor )
1147	CALL wrd_mpi_io( 'zeta_max', zeta_max )
1148	CALL wrd_mpi_io( 'zeta_min', zeta_min )
1149	CALL wrd_mpi_io_global_array( 'z_i', z_i )
1150
1151	ENDIF
1152
1153	!
1154	!-- Write restart data of the other modules
1155	CALL module_interface_wrd_global
1156
1157
1158	END SUBROUTINE wrd_global
1159
1160
1161	!--------------------------------------------------------------------------------------------------!
1162	! Description:
1163	! ------------
1164	!> Processor specific data of variables and arrays is written out for restarts (binary format).
1165	!> This information is written to the file opened by each PE.
1166	!--------------------------------------------------------------------------------------------------!
1167
1168
1169	SUBROUTINE wrd_local
1170
1171
1172	CHARACTER(LEN=10) :: binary_version_local !<
1173	CHARACTER(LEN=20) :: tmp_name !< temporary variable
1174
1175	INTEGER :: i !< loop index
1176
1177	!
1178	!-- Write arrays.
1179	binary_version_local = '5.1'
1180
1181	IF ( TRIM( restart_data_format_output ) == 'fortran_binary' ) THEN
1182	!
1183	!-- Output in Fortran binary format
1184
1185	CALL wrd_write_string( 'binary_version_local' )
1186	WRITE ( 14 ) binary_version_local
1187
1188	WRITE ( 14 ) numprocs, myid, nxl, nxr, nys, nyn, nzb, nzt
1189
1190	!
1191	!-- Attention: After changes to the following output commands the version number of the variable
1192	!-- ---------- binary_version_local must be changed!
1193	!-- Also, the list of arrays to be read in rrd_local must be adjusted accordingly.
1194	CALL wrd_write_string( 'e' )
1195	WRITE ( 14 ) e
1196
1197	IF ( ALLOCATED( e_av ) ) THEN
1198	CALL wrd_write_string( 'e_av' )
1199	WRITE ( 14 ) e_av
1200	ENDIF
1201
1202	CALL wrd_write_string( 'kh' )
1203	WRITE ( 14 ) kh
1204
1205
1206	IF ( ALLOCATED( kh_av ) ) THEN
1207	CALL wrd_write_string( 'kh_av' )
1208	WRITE ( 14 ) kh_av
1209	ENDIF
1210
1211	CALL wrd_write_string( 'km' )
1212	WRITE ( 14 ) km
1213
1214	IF ( ALLOCATED( km_av ) ) THEN
1215	CALL wrd_write_string( 'km_av' )
1216	WRITE ( 14 ) km_av
1217	ENDIF
1218
1219	IF ( ALLOCATED( lpt_av ) ) THEN
1220	CALL wrd_write_string( 'lpt_av' )
1221	WRITE ( 14 ) lpt_av
1222	ENDIF
1223
1224	IF ( ALLOCATED( lwp_av ) ) THEN
1225	CALL wrd_write_string( 'lwp_av' )
1226	WRITE ( 14 ) lwp_av
1227	ENDIF
1228
1229	CALL wrd_write_string( 'p' )
1230	WRITE ( 14 ) p
1231
1232	IF ( ALLOCATED( p_av ) ) THEN
1233	CALL wrd_write_string( 'p_av' )
1234	WRITE ( 14 ) p_av
1235	ENDIF
1236
1237	IF ( ALLOCATED( pc_av ) ) THEN
1238	CALL wrd_write_string( 'pc_av' )
1239	WRITE ( 14 ) pc_av
1240	ENDIF
1241
1242	IF ( ALLOCATED( pr_av ) ) THEN
1243	CALL wrd_write_string( 'pr_av' )
1244	WRITE ( 14 ) pr_av
1245	ENDIF
1246
1247	CALL wrd_write_string( 'pt' )
1248	WRITE ( 14 ) pt
1249
1250	IF ( ALLOCATED( pt_av ) ) THEN
1251	CALL wrd_write_string( 'pt_av' )
1252	WRITE ( 14 ) pt_av
1253	ENDIF
1254
1255	IF ( humidity ) THEN
1256
1257	CALL wrd_write_string( 'q' )
1258	WRITE ( 14 ) q
1259
1260	IF ( ALLOCATED( q_av ) ) THEN
1261	CALL wrd_write_string( 'q_av' )
1262	WRITE ( 14 ) q_av
1263	ENDIF
1264
1265	IF ( cloud_droplets ) THEN
1266
1267	CALL wrd_write_string( 'ql' )
1268	WRITE ( 14 ) ql
1269
1270	IF ( ALLOCATED( ql_av ) ) THEN
1271	CALL wrd_write_string( 'ql_av' )
1272	WRITE ( 14 ) ql_av
1273	ENDIF
1274
1275	ENDIF
1276
1277	IF ( ALLOCATED( qsurf_av ) ) THEN
1278	CALL wrd_write_string( 'qsurf_av' )
1279	WRITE ( 14 ) qsurf_av
1280	ENDIF
1281
1282	IF ( ALLOCATED( qsws_av ) ) THEN
1283	CALL wrd_write_string( 'qsws_av' )
1284	WRITE ( 14 ) qsws_av
1285	ENDIF
1286
1287	ENDIF
1288
1289	IF ( passive_scalar ) THEN
1290
1291	CALL wrd_write_string( 's' )
1292	WRITE ( 14 ) s
1293
1294	IF ( ALLOCATED( s_av ) ) THEN
1295	CALL wrd_write_string( 's_av' )
1296	WRITE ( 14 ) s_av
1297	ENDIF
1298
1299	IF ( ALLOCATED( ssurf_av ) ) THEN
1300	CALL wrd_write_string( 'ssurf_av' )
1301	WRITE ( 14 ) ssurf_av
1302	ENDIF
1303
1304	IF ( ALLOCATED( ssws_av ) ) THEN
1305	CALL wrd_write_string( 'ssws_av' )
1306	WRITE ( 14 ) ssws_av
1307	ENDIF
1308
1309	ENDIF
1310
1311	IF ( ALLOCATED( ql_c_av ) ) THEN
1312	CALL wrd_write_string( 'ql_c_av' )
1313	WRITE ( 14 ) ql_c_av
1314	ENDIF
1315
1316	IF ( ALLOCATED( ql_v_av ) ) THEN
1317	CALL wrd_write_string( 'ql_v_av' )
1318	WRITE ( 14 ) ql_v_av
1319	ENDIF
1320
1321	IF ( ALLOCATED( ql_vp_av ) ) THEN
1322	CALL wrd_write_string( 'ql_vp_av' )
1323	WRITE ( 14 ) ql_vp_av
1324	ENDIF
1325
1326	IF ( ALLOCATED( qv_av ) ) THEN
1327	CALL wrd_write_string( 'qv_av' )
1328	WRITE ( 14 ) qv_av
1329	ENDIF
1330
1331	CALL wrd_write_string( 'random_iv' )
1332	WRITE ( 14 ) random_iv
1333	WRITE ( 14 ) random_iy
1334
1335	IF ( ALLOCATED( seq_random_array ) ) THEN
1336	CALL wrd_write_string( 'seq_random_array' )
1337	WRITE ( 14 ) id_random_array
1338	WRITE ( 14 ) seq_random_array
1339	ENDIF
1340
1341	IF ( ALLOCATED( s_av ) ) THEN
1342	CALL wrd_write_string( 's_av' )
1343	WRITE ( 14 ) s_av
1344	ENDIF
1345
1346	IF ( ALLOCATED( shf_av ) ) THEN
1347	CALL wrd_write_string( 'shf_av' )
1348	WRITE ( 14 ) shf_av
1349	ENDIF
1350
1351	IF ( ALLOCATED( ts_av ) ) THEN
1352	CALL wrd_write_string( 'ts_av' )
1353	WRITE ( 14 ) ts_av
1354	ENDIF
1355
1356	CALL wrd_write_string( 'u' )
1357	WRITE ( 14 ) u
1358
1359	IF ( ALLOCATED( u_av ) ) THEN
1360	CALL wrd_write_string( 'u_av' )
1361	WRITE ( 14 ) u_av
1362	ENDIF
1363
1364	IF ( ALLOCATED( u_m_l ) ) THEN
1365	CALL wrd_write_string( 'u_m_l' )
1366	WRITE ( 14 ) u_m_l
1367	ENDIF
1368
1369	IF ( ALLOCATED( u_m_n ) ) THEN
1370	CALL wrd_write_string( 'u_m_n' )
1371	WRITE ( 14 ) u_m_n
1372	ENDIF
1373
1374	IF ( ALLOCATED( u_m_r ) ) THEN
1375	CALL wrd_write_string( 'u_m_r' )
1376	WRITE ( 14 ) u_m_r
1377	ENDIF
1378
1379	IF ( ALLOCATED( u_m_s ) ) THEN
1380	CALL wrd_write_string( 'u_m_s' )
1381	WRITE ( 14 ) u_m_s
1382	ENDIF
1383
1384	IF ( ALLOCATED( us_av ) ) THEN
1385	CALL wrd_write_string( 'us_av' )
1386	WRITE ( 14 ) us_av
1387	ENDIF
1388
1389	CALL wrd_write_string( 'v' )
1390	WRITE ( 14 ) v
1391
1392	IF ( ALLOCATED( v_av ) ) THEN
1393	CALL wrd_write_string( 'v_av' )
1394	WRITE ( 14 ) v_av
1395	ENDIF
1396
1397	IF ( ALLOCATED( v_m_l ) ) THEN
1398	CALL wrd_write_string( 'v_m_l' )
1399	WRITE ( 14 ) v_m_l
1400	ENDIF
1401
1402	IF ( ALLOCATED( v_m_n ) ) THEN
1403	CALL wrd_write_string( 'v_m_n' )
1404	WRITE ( 14 ) v_m_n
1405	ENDIF
1406
1407	IF ( ALLOCATED( v_m_r ) ) THEN
1408	CALL wrd_write_string( 'v_m_r' )
1409	WRITE ( 14 ) v_m_r
1410	ENDIF
1411
1412	IF ( ALLOCATED( v_m_s ) ) THEN
1413	CALL wrd_write_string( 'v_m_s' )
1414	WRITE ( 14 ) v_m_s
1415	ENDIF
1416
1417	IF ( humidity ) THEN
1418
1419	CALL wrd_write_string( 'vpt' )
1420	WRITE ( 14 ) vpt
1421
1422	IF ( ALLOCATED( vpt_av ) ) THEN
1423	CALL wrd_write_string( 'vpt_av' )
1424	WRITE ( 14 ) vpt_av
1425	ENDIF
1426
1427	ENDIF
1428
1429	CALL wrd_write_string( 'w' )
1430	WRITE ( 14 ) w
1431
1432	IF ( ALLOCATED( w_av ) ) THEN
1433	CALL wrd_write_string( 'w_av' )
1434	WRITE ( 14 ) w_av
1435	ENDIF
1436
1437	IF ( ALLOCATED( w_m_l ) ) THEN
1438	CALL wrd_write_string( 'w_m_l' )
1439	WRITE ( 14 ) w_m_l
1440	ENDIF
1441
1442	IF ( ALLOCATED( w_m_n ) ) THEN
1443	CALL wrd_write_string( 'w_m_n' )
1444	WRITE ( 14 ) w_m_n
1445	ENDIF
1446
1447	IF ( ALLOCATED( w_m_r ) ) THEN
1448	CALL wrd_write_string( 'w_m_r' )
1449	WRITE ( 14 ) w_m_r
1450	ENDIF
1451
1452	IF ( ALLOCATED( w_m_s ) ) THEN
1453	CALL wrd_write_string( 'w_m_s' )
1454	WRITE ( 14 ) w_m_s
1455	ENDIF
1456
1457	IF ( ALLOCATED( z0_av ) ) THEN
1458	CALL wrd_write_string( 'z0_av' )
1459	WRITE ( 14 ) z0_av
1460	ENDIF
1461
1462	IF ( ALLOCATED( z0h_av ) ) THEN
1463	CALL wrd_write_string( 'z0h_av' )
1464	WRITE ( 14 ) z0h_av
1465	ENDIF
1466
1467	IF ( ALLOCATED( z0q_av ) ) THEN
1468	CALL wrd_write_string( 'z0q_av' )
1469	WRITE ( 14 ) z0q_av
1470	ENDIF
1471
1472	IF ( land_surface .OR. urban_surface ) THEN
1473
1474	IF ( ALLOCATED( ghf_av ) ) THEN
1475	CALL wrd_write_string( 'ghf_av' )
1476	WRITE ( 14 ) ghf_av
1477	ENDIF
1478
1479	IF ( ALLOCATED( r_a_av ) ) THEN
1480	CALL wrd_write_string( 'r_a_av' )
1481	WRITE ( 14 ) r_a_av
1482	ENDIF
1483
1484	ENDIF
1485
1486	IF ( ALLOCATED( tsurf_av ) ) THEN
1487	CALL wrd_write_string( 'tsurf_av' )
1488	WRITE ( 14 ) tsurf_av
1489	ENDIF
1490
1491	ELSEIF ( restart_data_format_output(1:3) == 'mpi' ) THEN
1492	!
1493	!-- Write local restart data using MPI-IO
1494	!-- Attention: After changes to the following output commands the version number
1495	!-- --------- of the variable binary_version must be changed!
1496	!-- Also, the list of arrays to be read in rrd_local must be adjusted accordingly.
1497	CALL wrd_mpi_io( 'e', e)
1498	IF ( ALLOCATED( e_av ) ) CALL wrd_mpi_io( 'e_av', e_av)
1499	CALL wrd_mpi_io( 'kh', kh )
1500	IF ( ALLOCATED( kh_av ) ) CALL wrd_mpi_io( 'kh_av', kh_av )
1501	CALL wrd_mpi_io( 'km' , km )
1502	IF ( ALLOCATED( km_av ) ) CALL wrd_mpi_io( 'km_av', km_av )
1503	IF ( ALLOCATED( lpt_av ) ) CALL wrd_mpi_io( 'lpt_av', lpt_av )
1504	IF ( ALLOCATED( lwp_av ) ) CALL wrd_mpi_io( 'lwp_av', lwp_av )
1505	CALL wrd_mpi_io( 'p', p )
1506	IF ( ALLOCATED( p_av ) ) CALL wrd_mpi_io( 'p_av', p_av )
1507	IF ( ALLOCATED( pc_av ) ) CALL wrd_mpi_io( 'pc_av', pc_av )
1508	IF ( ALLOCATED( pr_av ) ) CALL wrd_mpi_io( 'pr_av', pr_av )
1509	CALL wrd_mpi_io( 'pt', pt )
1510	IF ( ALLOCATED( pt_av ) ) CALL wrd_mpi_io( 'pt_av', pt_av )
1511
1512	IF ( humidity ) THEN
1513
1514	CALL wrd_mpi_io( 'q', q )
1515	IF ( ALLOCATED( q_av ) ) CALL wrd_mpi_io( 'q_av', q_av )
1516
1517	IF ( cloud_droplets ) THEN
1518
1519	CALL wrd_mpi_io( 'ql', ql )
1520	IF ( ALLOCATED( ql_av ) ) CALL wrd_mpi_io( 'ql_av', ql_av )
1521
1522	ENDIF
1523
1524	IF ( ALLOCATED( qsurf_av ) ) CALL wrd_mpi_io( 'qsurf_av', qsurf_av )
1525	IF ( ALLOCATED( qsws_av ) ) CALL wrd_mpi_io( 'qsws_av', qsws_av )
1526
1527	ENDIF
1528
1529	IF ( passive_scalar ) THEN
1530
1531	CALL wrd_mpi_io( 's', s )
1532	IF ( ALLOCATED( s_av ) ) CALL wrd_mpi_io( 's_av', s_av )
1533	IF ( ALLOCATED( ssurf_av ) ) CALL wrd_mpi_io( 'ssurf_av', ssurf_av )
1534	IF ( ALLOCATED( ssws_av ) ) CALL wrd_mpi_io( 'ssws_av', ssws_av )
1535
1536	ENDIF
1537
1538	IF ( ALLOCATED( ql_c_av ) ) CALL wrd_mpi_io( 'ql_c_av', ql_c_av )
1539	IF ( ALLOCATED( ql_v_av ) ) CALL wrd_mpi_io( 'ql_v_av', ql_v_av )
1540	IF ( ALLOCATED( ql_vp_av ) ) CALL wrd_mpi_io( 'ql_vp_av', ql_vp_av )
1541	IF ( ALLOCATED( qv_av ) ) CALL wrd_mpi_io( 'qv_av', qv_av )
1542	!
1543	!-- Only PE0 writes random_iv and random_iy to restart file.
1544	!-- ATTENTION: If one value for every PE is required, the general approach of PE indendent
1545	!-- restart will be lost. That means that in general the parallel random number generator
1546	!-- in random_generatot_parallel_mod should be used!
1547	CALL wrd_mpi_io_global_array( 'random_iv', random_iv )
1548	CALL wrd_mpi_io( 'random_iy', random_iy )
1549
1550	IF ( ALLOCATED( seq_random_array ) ) THEN
1551	CALL wrd_mpi_io( 'id_random_array', id_random_array )
1552	DO i = 1, SIZE( seq_random_array, 1 )
1553	WRITE( tmp_name, '(A,I2.2)' ) 'seq_random_array', i
1554	CALL wrd_mpi_io( TRIM( tmp_name ), seq_random_array(i,:,:) )
1555	ENDDO
1556	ENDIF
1557	IF ( ALLOCATED( s_av ) ) CALL wrd_mpi_io( 's_av', s_av )
1558	IF ( ALLOCATED( shf_av ) ) CALL wrd_mpi_io( 'shf_av', shf_av )
1559	IF ( ALLOCATED( ts_av ) ) CALL wrd_mpi_io( 'ts_av', ts_av )
1560	CALL wrd_mpi_io( 'u', u)
1561	IF ( ALLOCATED( u_av ) ) CALL wrd_mpi_io( 'u_av', u_av )
1562	IF ( ALLOCATED( u_m_l ) ) CALL wrd_mpi_io( 'u_m_l', u_m_l )
1563	IF ( ALLOCATED( u_m_n ) ) CALL wrd_mpi_io( 'u_m_n', u_m_n )
1564	IF ( ALLOCATED( u_m_r ) ) CALL wrd_mpi_io( 'u_m_r', u_m_r )
1565	IF ( ALLOCATED( u_m_s ) ) CALL wrd_mpi_io( 'u_m_s', u_m_s )
1566	IF ( ALLOCATED( us_av ) ) CALL wrd_mpi_io( 'us_av', us_av )
1567	CALL wrd_mpi_io( 'v', v )
1568	IF ( ALLOCATED( v_av ) ) CALL wrd_mpi_io( 'v_av', v_av )
1569	IF ( ALLOCATED( v_m_l ) ) CALL wrd_mpi_io( 'v_m_l', v_m_l )
1570	IF ( ALLOCATED( v_m_n ) ) CALL wrd_mpi_io( 'v_m_n', v_m_n )
1571	IF ( ALLOCATED( v_m_r ) ) CALL wrd_mpi_io( 'v_m_r', v_m_r )
1572	IF ( ALLOCATED( v_m_s ) ) CALL wrd_mpi_io( 'v_m_s', v_m_s )
1573	IF ( humidity ) THEN
1574	CALL wrd_mpi_io( 'vpt', vpt )
1575	IF ( ALLOCATED( vpt_av ) ) CALL wrd_mpi_io( 'vpt_av', vpt_av )
1576	ENDIF
1577	CALL wrd_mpi_io( 'w', w)
1578	IF ( ALLOCATED( w_av ) ) CALL wrd_mpi_io( 'w_av', w_av )
1579	IF ( ALLOCATED( w_m_l ) ) CALL wrd_mpi_io( 'w_m_l', w_m_l )
1580	IF ( ALLOCATED( w_m_n ) ) CALL wrd_mpi_io( 'w_m_n', w_m_n )
1581	IF ( ALLOCATED( w_m_r ) ) CALL wrd_mpi_io( 'w_m_r', w_m_r )
1582	IF ( ALLOCATED( w_m_s ) ) CALL wrd_mpi_io( 'w_m_s', w_m_s )
1583	IF ( ALLOCATED( z0_av ) ) CALL wrd_mpi_io( 'z0_av', z0_av )
1584	IF ( ALLOCATED( z0h_av ) ) CALL wrd_mpi_io( 'z0h_av', z0h_av )
1585	IF ( ALLOCATED( z0q_av ) ) CALL wrd_mpi_io( 'z0q_av', z0q_av )
1586	IF ( land_surface .OR. urban_surface ) THEN
1587	IF ( ALLOCATED( ghf_av ) ) CALL wrd_mpi_io( 'ghf_av', ghf_av )
1588	IF ( ALLOCATED( r_a_av ) ) CALL wrd_mpi_io( 'r_a_av', r_a_av )
1589	ENDIF
1590	IF ( ALLOCATED( tsurf_av ) ) CALL wrd_mpi_io( 'tsurf_av', tsurf_av )
1591
1592	ENDIF
1593
1594	CALL surface_wrd_local
1595
1596	!
1597	!-- Write restart data of other modules
1598	CALL module_interface_wrd_local
1599
1600	!
1601	!-- Write end label
1602	IF ( restart_data_format_output == 'fortran_binary' ) THEN
1603	CALL wrd_write_string( '* end *' )
1604	ENDIF
1605
1606	END SUBROUTINE wrd_local
1607
1608
1609	END MODULE write_restart_data_mod

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |