Home

Context Navigation

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

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