Home

Context Navigation

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

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