Home

Context Navigation

write_restart_data_mod.f90 @ 4118

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

timestep.f90:

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

other:

remove dt_old

Property svn:keywords set to Id

File size: 34.1 KB

Rev	Line
[2894]	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
	6	! terms of the GNU General Public License as published by the Free Software
	7	! Foundation, either version 3 of the License, or (at your option) any later
	8	! version.
	9	!
	10	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
	11	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
	12	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
	13	!
	14	! You should have received a copy of the GNU General Public License along with
	15	! PALM. If not, see <http://www.gnu.org/licenses/>.
	16	!
[3655]	17	! Copyright 1997-2019 Leibniz Universitaet Hannover
[2894]	18	!------------------------------------------------------------------------------!
	19	!
	20	! Current revisions:
	21	! -----------------
	22	!
	23	!
	24	! Former revisions:
	25	! -----------------
	26	! $Id: write_restart_data_mod.f90 4101 2019-07-17 15:14:26Z suehring $
[4101]	27	! remove old_dt
	28	!
	29	! 4039 2019-06-18 10:32:41Z suehring
[4039]	30	! Modularize diagnostic output
	31	!
	32	! 4017 2019-06-06 12:16:46Z schwenkel
[3994]	33	! output of turbulence intensity added
	34	!
	35	! 3988 2019-05-22 11:32:37Z kanani
[3988]	36	! + time_virtual_measurement (to enable steering of output interval)
	37	!
	38	! 3936 2019-04-26 15:38:02Z kanani
[3936]	39	! Enable time-averaged output of theta_2m* with restarts
	40	!
	41	! 3668 2019-01-14 12:49:24Z maronga
[3637]	42	! Implementation of the PALM module interface
	43	!
	44	! 3589 2018-11-30 15:09:51Z suehring
[3589]	45	! Move the control parameter "salsa" from salsa_mod to control_parameters
	46	! (M. Kurppa)
	47	!
	48	!
	49	! 3582 2018-11-29 19:16:36Z suehring
[3467]	50	! Implementation of a new aerosol module salsa.
	51	!
	52	! 3355 2018-10-16 14:03:34Z knoop
[3294]	53	! changes concerning modularization of ocean option
	54	!
	55	! 3274 2018-09-24 15:42:55Z knoop
[3274]	56	! Modularization of all bulk cloud physics code components
	57	!
	58	! 3065 2018-06-12 07:03:02Z Giersch
[3065]	59	! New parameters concerning vertical grid stretching have been added
	60	!
	61	! 3004 2018-04-27 12:33:25Z Giersch
[3004]	62	! precipitation_rate_av removed
	63	!
	64	! 3003 2018-04-23 10:22:58Z Giersch
[3003]	65	! z_i is also written out to use the last known inversion height from the
	66	! initial run as the first inversion height which is written into the
	67	! run control file
	68	!
	69	! 2956 2018-04-10 11:01:03Z Giersch
[2956]	70	! spectrum_x and spectrum_y have been moved to global data
	71	!
	72	! 2921 2018-03-22 15:05:23Z Giersch
[2921]	73	! spinup_time, day_of_year_init and time_utc_init are also written out now
	74	!
	75	! 2912 2018-03-20 13:00:05Z knoop
[2912]	76	! Added gust module interface calls
	77	!
	78	! 2894 2018-03-15 09:17:58Z Giersch
[2894]	79	! Initial revision
	80	!
	81	!
	82	! Description:
	83	! ------------
	84	!> Writes restart data into binary file(s) for restart runs.
	85	!------------------------------------------------------------------------------!
	86	MODULE write_restart_data_mod
	87
	88
[3637]	89	USE arrays_3d, &
	90	ONLY: inflow_damping_factor, mean_inflow_profiles, pt_init, &
	91	q_init, ref_state, s_init, u_init, ug, v_init, vg, &
	92	e, kh, km, p, pt, q, ql, s, u, u_m_l, u_m_n, u_m_r, &
	93	u_m_s, v, v_m_l, v_m_n, v_m_r, v_m_s, vpt, w, w_m_l, w_m_n, &
	94	w_m_r, w_m_s
	95
	96	USE averaging
	97
[2894]	98	USE control_parameters
	99
[3637]	100	USE date_and_time_mod, &
	101	ONLY: day_of_year_init, time_utc_init
	102
	103	USE grid_variables, &
	104	ONLY: dx, dy
	105
	106	USE indices, &
	107	ONLY: nx, nxl, nxr, ny, nys, nyn, nz, nzb, nzt
	108
	109	USE kinds
	110
	111	USE model_1d_mod, &
	112	ONLY: damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
	113
	114	USE module_interface, &
	115	ONLY: module_interface_wrd_global, &
	116	module_interface_wrd_local
	117
	118	USE netcdf_interface, &
	119	ONLY: netcdf_precision, output_for_t0
	120
[2894]	121	USE pegrid, &
[3637]	122	ONLY: collective_wait, hor_index_bounds, myid, numprocs
[2894]	123
[3637]	124	USE radiation_model_mod, &
	125	ONLY: time_radiation
	126
	127	USE random_function_mod, &
	128	ONLY: random_iv, random_iy
	129
	130	USE random_generator_parallel, &
	131	ONLY: id_random_array, seq_random_array
	132
	133	USE spectra_mod, &
	134	ONLY: average_count_sp, spectrum_x, spectrum_y
	135
	136	USE statistics, &
	137	ONLY: statistic_regions, hom, hom_sum, u_max, u_max_ijk, v_max, &
	138	v_max_ijk, w_max, w_max_ijk, z_i
	139
	140	USE surface_mod, &
	141	ONLY : surface_wrd_local
	142
	143	USE vertical_nesting_mod, &
	144	ONLY: vnest_init
	145
[3988]	146	USE virtual_measurement_mod, &
	147	ONLY: time_virtual_measurement
[3637]	148
[3988]	149
[2894]	150	IMPLICIT NONE
	151
	152
	153	INTERFACE wrd_global
	154	MODULE PROCEDURE wrd_global
	155	END INTERFACE wrd_global
	156
	157	INTERFACE wrd_local
	158	MODULE PROCEDURE wrd_local
	159	END INTERFACE wrd_local
	160
	161
	162	PUBLIC wrd_local, wrd_global
	163
	164
	165	CONTAINS
	166
	167
	168	! Description:
	169	! ------------
	170	!> Global data of control variables and arrays is written out for
	171	!> restarts (binary format).
	172	!> This information is only written to the file opened by PE0.
	173	!------------------------------------------------------------------------------!
	174	SUBROUTINE wrd_global
	175
	176
	177	CHARACTER (LEN=10) :: binary_version_global !<
	178
	179
[3668]	180	binary_version_global = '4.8'
[2894]	181
	182	CALL wrd_write_string( 'binary_version_global' )
	183	WRITE ( 14 ) binary_version_global
	184
	185	CALL wrd_write_string( 'numprocs' )
	186	WRITE ( 14 ) numprocs
	187
	188	CALL wrd_write_string( 'hor_index_bounds' )
	189	WRITE ( 14 ) hor_index_bounds
	190
	191	CALL wrd_write_string( 'nz' )
	192	WRITE ( 14 ) nz
	193
	194	CALL wrd_write_string( 'max_pr_user' )
	195	WRITE ( 14 ) max_pr_user
	196
	197	CALL wrd_write_string( 'statistic_regions' )
	198	WRITE ( 14 ) statistic_regions
	199
	200	!
	201	!-- Caution: After changes in the following parameter-list, the
	202	!-- ------- version number stored in the variable binary_version_global has to
	203	!-- be increased. The same changes must also be done in the parameter-
	204	!-- list in rrd_global.
	205
	206	CALL wrd_write_string( 'advected_distance_x' )
	207	WRITE ( 14 ) advected_distance_x
	208
	209	CALL wrd_write_string( 'advected_distance_y' )
	210	WRITE ( 14 ) advected_distance_y
	211
	212	CALL wrd_write_string( 'alpha_surface' )
	213	WRITE ( 14 ) alpha_surface
	214
	215	CALL wrd_write_string( 'average_count_pr' )
	216	WRITE ( 14 ) average_count_pr
	217
	218	CALL wrd_write_string( 'average_count_sp' )
	219	WRITE ( 14 ) average_count_sp
	220
	221	CALL wrd_write_string( 'average_count_3d' )
	222	WRITE ( 14 ) average_count_3d
	223
	224	CALL wrd_write_string( 'bc_e_b' )
	225	WRITE ( 14 ) bc_e_b
	226
	227	CALL wrd_write_string( 'bc_lr' )
	228	WRITE ( 14 ) bc_lr
	229
	230	CALL wrd_write_string( 'bc_ns' )
	231	WRITE ( 14 ) bc_ns
	232
	233	CALL wrd_write_string( 'bc_p_b' )
	234	WRITE ( 14 ) bc_p_b
	235
	236	CALL wrd_write_string( 'bc_p_t' )
	237	WRITE ( 14 ) bc_p_t
	238
	239	CALL wrd_write_string( 'bc_pt_b' )
	240	WRITE ( 14 ) bc_pt_b
	241
	242	CALL wrd_write_string( 'bc_pt_t' )
	243	WRITE ( 14 ) bc_pt_t
	244
	245	CALL wrd_write_string( 'bc_pt_t_val' )
	246	WRITE ( 14 ) bc_pt_t_val
	247
	248	CALL wrd_write_string( 'bc_q_b' )
	249	WRITE ( 14 ) bc_q_b
	250
	251	CALL wrd_write_string( 'bc_q_t' )
	252	WRITE ( 14 ) bc_q_t
	253
	254	CALL wrd_write_string( 'bc_q_t_val' )
	255	WRITE ( 14 ) bc_q_t_val
	256
	257	CALL wrd_write_string( 'bc_s_b' )
	258	WRITE ( 14 ) bc_s_b
	259
	260	CALL wrd_write_string( 'bc_s_t' )
	261	WRITE ( 14 ) bc_s_t
	262
	263	CALL wrd_write_string( 'bc_uv_b' )
	264	WRITE ( 14 ) bc_uv_b
	265
	266	CALL wrd_write_string( 'bc_uv_t' )
	267	WRITE ( 14 ) bc_uv_t
	268
	269	CALL wrd_write_string( 'building_height' )
	270	WRITE ( 14 ) building_height
	271
	272	CALL wrd_write_string( 'building_length_x' )
	273	WRITE ( 14 ) building_length_x
	274
	275	CALL wrd_write_string( 'building_length_y' )
	276	WRITE ( 14 ) building_length_y
	277
	278	CALL wrd_write_string( 'building_wall_left' )
	279	WRITE ( 14 ) building_wall_left
	280
	281	CALL wrd_write_string( 'building_wall_south' )
	282	WRITE ( 14 ) building_wall_south
	283
	284	CALL wrd_write_string( 'call_psolver_at_all_substeps' )
	285	WRITE ( 14 ) call_psolver_at_all_substeps
	286
	287	CALL wrd_write_string( 'canyon_height' )
	288	WRITE ( 14 ) canyon_height
	289
	290	CALL wrd_write_string( 'canyon_wall_left' )
	291	WRITE ( 14 ) canyon_wall_left
	292
	293	CALL wrd_write_string( 'canyon_wall_south' )
	294	WRITE ( 14 ) canyon_wall_south
	295
	296	CALL wrd_write_string( 'canyon_width_x' )
	297	WRITE ( 14 ) canyon_width_x
	298
	299	CALL wrd_write_string( 'canyon_width_y' )
	300	WRITE ( 14 ) canyon_width_y
	301
	302	CALL wrd_write_string( 'cfl_factor' )
	303	WRITE ( 14 ) cfl_factor
	304
	305	CALL wrd_write_string( 'cloud_droplets' )
	306	WRITE ( 14 ) cloud_droplets
	307
	308	CALL wrd_write_string( 'collective_wait' )
	309	WRITE ( 14 ) collective_wait
	310
	311	CALL wrd_write_string( 'conserve_volume_flow' )
	312	WRITE ( 14 ) conserve_volume_flow
	313
	314	CALL wrd_write_string( 'conserve_volume_flow_mode' )
	315	WRITE ( 14 ) conserve_volume_flow_mode
	316
	317	CALL wrd_write_string( 'constant_flux_layer' )
	318	WRITE ( 14 ) constant_flux_layer
	319
	320	CALL wrd_write_string( 'coupling_start_time' )
	321	WRITE ( 14 ) coupling_start_time
	322
	323	CALL wrd_write_string( 'current_timestep_number' )
	324	WRITE ( 14 ) current_timestep_number
	325
	326	CALL wrd_write_string( 'cycle_mg' )
	327	WRITE ( 14 ) cycle_mg
	328
	329	CALL wrd_write_string( 'damp_level_1d' )
	330	WRITE ( 14 ) damp_level_1d
	331
[2921]	332	CALL wrd_write_string( 'day_of_year_init' )
	333	WRITE ( 14 ) day_of_year_init
	334
[2894]	335	CALL wrd_write_string( 'dissipation_1d' )
	336	WRITE ( 14 ) dissipation_1d
	337
	338	CALL wrd_write_string( 'do2d_xy_time_count' )
	339	WRITE ( 14 ) do2d_xy_time_count
	340
	341	CALL wrd_write_string( 'do2d_xz_time_count' )
	342	WRITE ( 14 ) do2d_xz_time_count
	343
	344	CALL wrd_write_string( 'do2d_yz_time_count' )
	345	WRITE ( 14 ) do2d_yz_time_count
	346
	347	CALL wrd_write_string( 'do3d_time_count' )
	348	WRITE ( 14 ) do3d_time_count
	349
	350	CALL wrd_write_string( 'dp_external' )
	351	WRITE ( 14 ) dp_external
	352
	353	CALL wrd_write_string( 'dp_level_b' )
	354	WRITE ( 14 ) dp_level_b
	355
	356	CALL wrd_write_string( 'dp_smooth' )
	357	WRITE ( 14 ) dp_smooth
	358
	359	CALL wrd_write_string( 'dpdxy' )
	360	WRITE ( 14 ) dpdxy
	361
	362	CALL wrd_write_string( 'dt_3d' )
	363	WRITE ( 14 ) dt_3d
	364
	365	CALL wrd_write_string( 'dt_pr_1d' )
	366	WRITE ( 14 ) dt_pr_1d
	367
	368	CALL wrd_write_string( 'dt_run_control_1d' )
	369	WRITE ( 14 ) dt_run_control_1d
	370
	371	CALL wrd_write_string( 'dx' )
	372	WRITE ( 14 ) dx
	373
	374	CALL wrd_write_string( 'dy' )
	375	WRITE ( 14 ) dy
	376
	377	CALL wrd_write_string( 'dz' )
	378	WRITE ( 14 ) dz
[3065]	379
[2894]	380	CALL wrd_write_string( 'dz_max' )
	381	WRITE ( 14 ) dz_max
	382
	383	CALL wrd_write_string( 'dz_stretch_factor' )
	384	WRITE ( 14 ) dz_stretch_factor
[3065]	385
	386	CALL wrd_write_string( 'dz_stretch_factor_array' )
	387	WRITE ( 14 ) dz_stretch_factor_array
	388
[2894]	389	CALL wrd_write_string( 'dz_stretch_level' )
	390	WRITE ( 14 ) dz_stretch_level
	391
[3065]	392	CALL wrd_write_string( 'dz_stretch_level_end' )
	393	WRITE ( 14 ) dz_stretch_level_end
	394
	395	CALL wrd_write_string( 'dz_stretch_level_start' )
	396	WRITE ( 14 ) dz_stretch_level_start
	397
[2894]	398	CALL wrd_write_string( 'e_min' )
	399	WRITE ( 14 ) e_min
	400
	401	CALL wrd_write_string( 'end_time_1d' )
	402	WRITE ( 14 ) end_time_1d
	403
	404	CALL wrd_write_string( 'fft_method' )
	405	WRITE ( 14 ) fft_method
	406
	407	CALL wrd_write_string( 'first_call_lpm' )
	408	WRITE ( 14 ) first_call_lpm
	409
	410	CALL wrd_write_string( 'galilei_transformation' )
	411	WRITE ( 14 ) galilei_transformation
	412
	413	CALL wrd_write_string( 'hom' )
	414	WRITE ( 14 ) hom
	415
	416	CALL wrd_write_string( 'hom_sum' )
	417	WRITE ( 14 ) hom_sum
	418
	419	CALL wrd_write_string( 'humidity' )
	420	WRITE ( 14 ) humidity
	421
	422	IF ( ALLOCATED( inflow_damping_factor ) ) THEN
	423	CALL wrd_write_string( 'inflow_damping_factor' )
	424	WRITE ( 14 ) inflow_damping_factor
	425	ENDIF
	426
	427	CALL wrd_write_string( 'inflow_damping_height' )
	428	WRITE ( 14 ) inflow_damping_height
	429
	430	CALL wrd_write_string( 'inflow_damping_width' )
	431	WRITE ( 14 ) inflow_damping_width
	432
	433	CALL wrd_write_string( 'inflow_disturbance_begin' )
	434	WRITE ( 14 ) inflow_disturbance_begin
	435
	436	CALL wrd_write_string( 'inflow_disturbance_end' )
	437	WRITE ( 14 ) inflow_disturbance_end
	438
	439	CALL wrd_write_string( 'km_constant' )
	440	WRITE ( 14 ) km_constant
	441
	442	CALL wrd_write_string( 'large_scale_forcing' )
	443	WRITE ( 14 ) large_scale_forcing
	444
	445	CALL wrd_write_string( 'large_scale_subsidence' )
	446	WRITE ( 14 ) large_scale_subsidence
	447
	448	CALL wrd_write_string( 'latitude' )
	449	WRITE ( 14 ) latitude
	450
	451	CALL wrd_write_string( 'longitude' )
	452	WRITE ( 14 ) longitude
	453
	454	CALL wrd_write_string( 'loop_optimization' )
	455	WRITE ( 14 ) loop_optimization
	456
	457	CALL wrd_write_string( 'masking_method' )
	458	WRITE ( 14 ) masking_method
	459
	460	IF ( ALLOCATED( mean_inflow_profiles ) ) THEN
	461	CALL wrd_write_string( 'mean_inflow_profiles' )
	462	WRITE ( 14 ) mean_inflow_profiles
	463	ENDIF
	464
	465	CALL wrd_write_string( 'mg_cycles' )
	466	WRITE ( 14 ) mg_cycles
	467
	468	CALL wrd_write_string( 'mg_switch_to_pe0_level' )
	469	WRITE ( 14 ) mg_switch_to_pe0_level
	470
	471	CALL wrd_write_string( 'mixing_length_1d' )
	472	WRITE ( 14 ) mixing_length_1d
	473
	474	CALL wrd_write_string( 'momentum_advec' )
	475	WRITE ( 14 ) momentum_advec
	476
	477	CALL wrd_write_string( 'netcdf_precision' )
	478	WRITE ( 14 ) netcdf_precision
	479
	480	CALL wrd_write_string( 'neutral' )
	481	WRITE ( 14 ) neutral
	482
	483	CALL wrd_write_string( 'ngsrb' )
	484	WRITE ( 14 ) ngsrb
	485
	486	CALL wrd_write_string( 'nsor' )
	487	WRITE ( 14 ) nsor
	488
	489	CALL wrd_write_string( 'nsor_ini' )
	490	WRITE ( 14 ) nsor_ini
	491
	492	CALL wrd_write_string( 'nudging' )
	493	WRITE ( 14 ) nudging
	494
	495	CALL wrd_write_string( 'num_leg' )
	496	WRITE ( 14 ) num_leg
	497
	498	CALL wrd_write_string( 'nx' )
	499	WRITE ( 14 ) nx
	500
	501	CALL wrd_write_string( 'ny' )
	502	WRITE ( 14 ) ny
	503
[3294]	504	CALL wrd_write_string( 'ocean_mode' )
	505	WRITE ( 14 ) ocean_mode
[2894]	506
	507	CALL wrd_write_string( 'omega' )
	508	WRITE ( 14 ) omega
	509
	510	CALL wrd_write_string( 'omega_sor' )
	511	WRITE ( 14 ) omega_sor
	512
	513	CALL wrd_write_string( 'output_for_t0' )
	514	WRITE ( 14 ) output_for_t0
	515
	516	CALL wrd_write_string( 'passive_scalar' )
	517	WRITE ( 14 ) passive_scalar
	518
	519	CALL wrd_write_string( 'prandtl_number' )
	520	WRITE ( 14 ) prandtl_number
	521
	522	CALL wrd_write_string( 'psolver' )
	523	WRITE ( 14 ) psolver
	524
	525	CALL wrd_write_string( 'pt_damping_factor' )
	526	WRITE ( 14 ) pt_damping_factor
	527
	528	CALL wrd_write_string( 'pt_damping_width' )
	529	WRITE ( 14 ) pt_damping_width
	530
	531	CALL wrd_write_string( 'pt_init' )
	532	WRITE ( 14 ) pt_init
	533
	534	CALL wrd_write_string( 'pt_reference' )
	535	WRITE ( 14 ) pt_reference
	536
	537	CALL wrd_write_string( 'pt_surface' )
	538	WRITE ( 14 ) pt_surface
	539
	540	CALL wrd_write_string( 'pt_surface_initial_change' )
	541	WRITE ( 14 ) pt_surface_initial_change
	542
	543	CALL wrd_write_string( 'pt_vertical_gradient' )
	544	WRITE ( 14 ) pt_vertical_gradient
	545
	546	CALL wrd_write_string( 'pt_vertical_gradient_level' )
	547	WRITE ( 14 ) pt_vertical_gradient_level
	548
	549	CALL wrd_write_string( 'pt_vertical_gradient_level_ind' )
	550	WRITE ( 14 ) pt_vertical_gradient_level_ind
	551
	552	CALL wrd_write_string( 'q_init' )
	553	WRITE ( 14 ) q_init
	554
	555	CALL wrd_write_string( 'q_surface' )
	556	WRITE ( 14 ) q_surface
	557
	558	CALL wrd_write_string( 'q_surface_initial_change' )
	559	WRITE ( 14 ) q_surface_initial_change
	560
	561	CALL wrd_write_string( 'q_vertical_gradient' )
	562	WRITE ( 14 ) q_vertical_gradient
	563
	564	CALL wrd_write_string( 'q_vertical_gradient_level' )
	565	WRITE ( 14 ) q_vertical_gradient_level
	566
	567	CALL wrd_write_string( 'q_vertical_gradient_level_ind' )
	568	WRITE ( 14 ) q_vertical_gradient_level_ind
	569
	570	CALL wrd_write_string( 'random_generator' )
	571	WRITE ( 14 ) random_generator
	572
	573	CALL wrd_write_string( 'random_heatflux' )
	574	WRITE ( 14 ) random_heatflux
	575
	576	CALL wrd_write_string( 'rans_mode' )
	577	WRITE ( 14 ) rans_mode
	578
	579	CALL wrd_write_string( 'rayleigh_damping_factor' )
	580	WRITE ( 14 ) rayleigh_damping_factor
	581
	582	CALL wrd_write_string( 'rayleigh_damping_height' )
	583	WRITE ( 14 ) rayleigh_damping_height
	584
	585	CALL wrd_write_string( 'recycling_width' )
	586	WRITE ( 14 ) recycling_width
	587
	588	CALL wrd_write_string( 'recycling_yshift' )
	589	WRITE ( 14 ) recycling_yshift
	590
	591	CALL wrd_write_string( 'ref_state' )
	592	WRITE ( 14 ) ref_state
	593
	594	CALL wrd_write_string( 'reference_state' )
	595	WRITE ( 14 ) reference_state
	596
	597	CALL wrd_write_string( 'residual_limit' )
	598	WRITE ( 14 ) residual_limit
	599
	600	CALL wrd_write_string( 'roughness_length' )
	601	WRITE ( 14 ) roughness_length
	602
	603	CALL wrd_write_string( 'run_coupled' )
	604	WRITE ( 14 ) run_coupled
	605
	606	CALL wrd_write_string( 'runnr' )
	607	WRITE ( 14 ) runnr
	608
	609	CALL wrd_write_string( 's_init' )
	610	WRITE ( 14 ) s_init
	611
	612	CALL wrd_write_string( 's_surface' )
	613	WRITE ( 14 ) s_surface
	614
	615	CALL wrd_write_string( 's_surface_initial_change' )
	616	WRITE ( 14 ) s_surface_initial_change
	617
	618	CALL wrd_write_string( 's_vertical_gradient' )
	619	WRITE ( 14 ) s_vertical_gradient
	620
	621	CALL wrd_write_string( 's_vertical_gradient_level' )
	622	WRITE ( 14 ) s_vertical_gradient_level
	623
	624	CALL wrd_write_string( 's_vertical_gradient_level_ind' )
	625	WRITE ( 14 ) s_vertical_gradient_level_ind
	626
	627	CALL wrd_write_string( 'scalar_advec' )
	628	WRITE ( 14 ) scalar_advec
	629
	630	CALL wrd_write_string( 'simulated_time' )
	631	WRITE ( 14 ) simulated_time
	632
[2956]	633	IF ( ALLOCATED( spectrum_x ) ) THEN
	634	CALL wrd_write_string( 'spectrum_x' )
	635	WRITE ( 14 ) spectrum_x
	636	CALL wrd_write_string( 'spectrum_y' )
	637	WRITE ( 14 ) spectrum_y
	638	ENDIF
	639
[2921]	640	CALL wrd_write_string( 'spinup_time ' )
	641	WRITE ( 14 ) spinup_time
	642
[2894]	643	CALL wrd_write_string( 'surface_heatflux' )
	644	WRITE ( 14 ) surface_heatflux
	645
	646	CALL wrd_write_string( 'surface_pressure' )
	647	WRITE ( 14 ) surface_pressure
	648
	649	CALL wrd_write_string( 'surface_scalarflux' )
	650	WRITE ( 14 ) surface_scalarflux
	651
	652	CALL wrd_write_string( 'surface_waterflux' )
	653	WRITE ( 14 ) surface_waterflux
	654
	655	CALL wrd_write_string( 'time_coupling' )
	656	WRITE ( 14 ) time_coupling
	657
	658	CALL wrd_write_string( 'time_disturb' )
	659	WRITE ( 14 ) time_disturb
	660
	661	CALL wrd_write_string( 'time_do2d_xy' )
	662	WRITE ( 14 ) time_do2d_xy
	663
	664	CALL wrd_write_string( 'time_do2d_xz' )
	665	WRITE ( 14 ) time_do2d_xz
	666
	667	CALL wrd_write_string( 'time_do2d_yz' )
	668	WRITE ( 14 ) time_do2d_yz
	669
	670	CALL wrd_write_string( 'time_do3d' )
	671	WRITE ( 14 ) time_do3d
	672
	673	CALL wrd_write_string( 'time_do_av' )
	674	WRITE ( 14 ) time_do_av
	675
	676	CALL wrd_write_string( 'time_do_sla' )
	677	WRITE ( 14 ) time_do_sla
	678
	679	CALL wrd_write_string( 'time_domask' )
	680	WRITE ( 14 ) time_domask
	681
	682	CALL wrd_write_string( 'time_dopr' )
	683	WRITE ( 14 ) time_dopr
	684
	685	CALL wrd_write_string( 'time_dopr_av' )
	686	WRITE ( 14 ) time_dopr_av
	687
	688	CALL wrd_write_string( 'time_dopr_listing' )
	689	WRITE ( 14 ) time_dopr_listing
	690
	691	CALL wrd_write_string( 'time_dopts' )
	692	WRITE ( 14 ) time_dopts
	693
	694	CALL wrd_write_string( 'time_dosp' )
	695	WRITE ( 14 ) time_dosp
	696
	697	CALL wrd_write_string( 'time_dots' )
	698	WRITE ( 14 ) time_dots
	699
	700	CALL wrd_write_string( 'time_radiation' )
	701	WRITE ( 14 ) time_radiation
	702
	703	CALL wrd_write_string( 'time_restart' )
	704	WRITE ( 14 ) time_restart
	705
	706	CALL wrd_write_string( 'time_run_control' )
	707	WRITE ( 14 ) time_run_control
	708
	709	CALL wrd_write_string( 'time_since_reference_point' )
	710	WRITE ( 14 ) time_since_reference_point
	711
[2921]	712	CALL wrd_write_string( 'time_utc_init' )
	713	WRITE ( 14 ) time_utc_init
	714
[3988]	715	CALL wrd_write_string( 'time_virtual_measurement' )
	716	WRITE ( 14 ) time_virtual_measurement
	717
[2894]	718	CALL wrd_write_string( 'timestep_scheme' )
	719	WRITE ( 14 ) timestep_scheme
	720
	721	CALL wrd_write_string( 'top_heatflux' )
	722	WRITE ( 14 ) top_heatflux
	723
	724	CALL wrd_write_string( 'top_momentumflux_u' )
	725	WRITE ( 14 ) top_momentumflux_u
	726
	727	CALL wrd_write_string( 'top_momentumflux_v' )
	728	WRITE ( 14 ) top_momentumflux_v
	729
	730	CALL wrd_write_string( 'top_scalarflux' )
	731	WRITE ( 14 ) top_scalarflux
	732
	733	CALL wrd_write_string( 'topography' )
	734	WRITE ( 14 ) topography
	735
	736	CALL wrd_write_string( 'topography_grid_convention' )
	737	WRITE ( 14 ) topography_grid_convention
	738
	739	CALL wrd_write_string( 'tsc' )
	740	WRITE ( 14 ) tsc
	741
	742	CALL wrd_write_string( 'tunnel_height' )
	743	WRITE ( 14 ) tunnel_height
	744
	745	CALL wrd_write_string( 'tunnel_length' )
	746	WRITE ( 14 ) tunnel_length
	747
	748	CALL wrd_write_string( 'tunnel_wall_depth' )
	749	WRITE ( 14 ) tunnel_wall_depth
	750
	751	CALL wrd_write_string( 'tunnel_width_x' )
	752	WRITE ( 14 ) tunnel_width_x
	753
	754	CALL wrd_write_string( 'tunnel_width_y' )
	755	WRITE ( 14 ) tunnel_width_y
	756
	757	CALL wrd_write_string( 'turbulence_closure' )
	758	WRITE ( 14 ) turbulence_closure
	759
	760	CALL wrd_write_string( 'turbulent_inflow' )
	761	WRITE ( 14 ) turbulent_inflow
	762
	763	CALL wrd_write_string( 'u_bulk' )
	764	WRITE ( 14 ) u_bulk
	765
	766	CALL wrd_write_string( 'u_init' )
	767	WRITE ( 14 ) u_init
	768
	769	CALL wrd_write_string( 'u_max' )
	770	WRITE ( 14 ) u_max
	771
	772	CALL wrd_write_string( 'u_max_ijk' )
	773	WRITE ( 14 ) u_max_ijk
	774
	775	CALL wrd_write_string( 'ug' )
	776	WRITE ( 14 ) ug
	777
	778	CALL wrd_write_string( 'ug_surface' )
	779	WRITE ( 14 ) ug_surface
	780
	781	CALL wrd_write_string( 'ug_vertical_gradient' )
	782	WRITE ( 14 ) ug_vertical_gradient
	783
	784	CALL wrd_write_string( 'ug_vertical_gradient_level' )
	785	WRITE ( 14 ) ug_vertical_gradient_level
	786
	787	CALL wrd_write_string( 'ug_vertical_gradient_level_ind' )
	788	WRITE ( 14 ) ug_vertical_gradient_level_ind
	789
	790	CALL wrd_write_string( 'use_surface_fluxes' )
	791	WRITE ( 14 ) use_surface_fluxes
	792
	793	CALL wrd_write_string( 'use_top_fluxes' )
	794	WRITE ( 14 ) use_top_fluxes
	795
	796	CALL wrd_write_string( 'use_ug_for_galilei_tr' )
	797	WRITE ( 14 ) use_ug_for_galilei_tr
	798
	799	CALL wrd_write_string( 'use_upstream_for_tke' )
	800	WRITE ( 14 ) use_upstream_for_tke
	801
	802	CALL wrd_write_string( 'v_bulk' )
	803	WRITE ( 14 ) v_bulk
	804
	805	CALL wrd_write_string( 'v_init' )
	806	WRITE ( 14 ) v_init
	807
	808	CALL wrd_write_string( 'v_max' )
	809	WRITE ( 14 ) v_max
	810
	811	CALL wrd_write_string( 'v_max_ijk' )
	812	WRITE ( 14 ) v_max_ijk
	813
	814	CALL wrd_write_string( 'vg' )
	815	WRITE ( 14 ) vg
	816
	817	CALL wrd_write_string( 'vg_surface' )
	818	WRITE ( 14 ) vg_surface
	819
	820	CALL wrd_write_string( 'vg_vertical_gradient' )
	821	WRITE ( 14 ) vg_vertical_gradient
	822
	823	CALL wrd_write_string( 'vg_vertical_gradient_level' )
	824	WRITE ( 14 ) vg_vertical_gradient_level
	825
	826	CALL wrd_write_string( 'vg_vertical_gradient_level_ind' )
	827	WRITE ( 14 ) vg_vertical_gradient_level_ind
	828
	829	CALL wrd_write_string( 'virtual_flight' )
	830	WRITE ( 14 ) virtual_flight
	831
	832	CALL wrd_write_string( 'vnest_init' )
	833	WRITE ( 14 ) vnest_init
	834
	835	CALL wrd_write_string( 'volume_flow_area' )
	836	WRITE ( 14 ) volume_flow_area
	837
	838	CALL wrd_write_string( 'volume_flow_initial' )
	839	WRITE ( 14 ) volume_flow_initial
	840
	841	CALL wrd_write_string( 'subs_vertical_gradient' )
	842	WRITE ( 14 ) subs_vertical_gradient
	843
	844	CALL wrd_write_string( 'subs_vertical_gradient_level' )
	845	WRITE ( 14 ) subs_vertical_gradient_level
	846
	847	CALL wrd_write_string( 'subs_vertical_gradient_level_i' )
	848	WRITE ( 14 ) subs_vertical_gradient_level_i
	849
	850	CALL wrd_write_string( 'w_max' )
	851	WRITE ( 14 ) w_max
	852
	853	CALL wrd_write_string( 'w_max_ijk' )
	854	WRITE ( 14 ) w_max_ijk
	855
	856	CALL wrd_write_string( 'wall_adjustment' )
	857	WRITE ( 14 ) wall_adjustment
	858
	859	CALL wrd_write_string( 'wall_heatflux' )
	860	WRITE ( 14 ) wall_heatflux
	861
	862	CALL wrd_write_string( 'wall_humidityflux' )
	863	WRITE ( 14 ) wall_humidityflux
	864
	865	CALL wrd_write_string( 'wall_scalarflux' )
	866	WRITE ( 14 ) wall_scalarflux
	867
	868	CALL wrd_write_string( 'y_shift' )
	869	WRITE ( 14 ) y_shift
	870
	871	CALL wrd_write_string( 'z0h_factor' )
	872	WRITE ( 14 ) z0h_factor
	873
	874	CALL wrd_write_string( 'zeta_max' )
	875	WRITE ( 14 ) zeta_max
	876
	877	CALL wrd_write_string( 'zeta_min' )
	878	WRITE ( 14 ) zeta_min
	879
[3003]	880	CALL wrd_write_string( 'z_i' )
	881	WRITE ( 14 ) z_i
[2894]	882
	883	!
[3294]	884	!-- Write restart data of the other modules
[3637]	885	CALL module_interface_wrd_global
[2894]	886
	887
	888	END SUBROUTINE wrd_global
	889
	890
	891	!------------------------------------------------------------------------------!
	892	! Description:
	893	! ------------
	894	!> Processor specific data of variables and arrays is written out for
	895	!> restarts (binary format).
	896	!> This information is written to the file opened by each PE.
	897	!------------------------------------------------------------------------------!
	898	SUBROUTINE wrd_local
	899
	900
	901	CHARACTER (LEN=10) :: binary_version_local !<
	902
	903	!
	904	!-- Write arrays.
	905	binary_version_local = '4.7'
	906
	907	CALL wrd_write_string( 'binary_version_local' )
	908	WRITE ( 14 ) binary_version_local
	909
	910	WRITE ( 14 ) numprocs, myid, nxl, nxr, nys, nyn, nzb, nzt
	911
	912	!
	913	!-- Attention: After changes to the following output commands the version number
	914	!-- --------- of the variable binary_version_local must be changed!
	915	!-- Also, the list of arrays to be read in rrd_local must be
	916	!-- adjusted accordingly.
	917	CALL wrd_write_string( 'e' )
	918	WRITE ( 14 ) e
	919
	920	IF ( ALLOCATED( e_av ) ) THEN
	921	CALL wrd_write_string( 'e_av' )
	922	WRITE ( 14 ) e_av
	923	ENDIF
	924
	925	CALL wrd_write_string( 'kh' )
	926	WRITE ( 14 ) kh
	927
	928
	929	IF ( ALLOCATED( kh_av ) ) THEN
	930	CALL wrd_write_string( 'kh_av' )
	931	WRITE ( 14 ) kh_av
	932	ENDIF
	933
	934	CALL wrd_write_string( 'km' )
	935	WRITE ( 14 ) km
	936
	937	IF ( ALLOCATED( km_av ) ) THEN
	938	CALL wrd_write_string( 'km_av' )
	939	WRITE ( 14 ) km_av
	940	ENDIF
	941
	942	IF ( ALLOCATED( lpt_av ) ) THEN
	943	CALL wrd_write_string( 'lpt_av' )
	944	WRITE ( 14 ) lpt_av
	945	ENDIF
	946
	947	IF ( ALLOCATED( lwp_av ) ) THEN
	948	CALL wrd_write_string( 'lwp_av' )
	949	WRITE ( 14 ) lwp_av
	950	ENDIF
	951
	952	CALL wrd_write_string( 'p' )
	953	WRITE ( 14 ) p
	954
	955	IF ( ALLOCATED( p_av ) ) THEN
	956	CALL wrd_write_string( 'p_av' )
	957	WRITE ( 14 ) p_av
	958	ENDIF
	959
	960	IF ( ALLOCATED( pc_av ) ) THEN
	961	CALL wrd_write_string( 'pc_av' )
	962	WRITE ( 14 ) pc_av
	963	ENDIF
	964
	965	IF ( ALLOCATED( pr_av ) ) THEN
	966	CALL wrd_write_string( 'pr_av' )
	967	WRITE ( 14 ) pr_av
	968	ENDIF
	969
	970	CALL wrd_write_string( 'pt' )
	971	WRITE ( 14 ) pt
	972
	973	IF ( ALLOCATED( pt_av ) ) THEN
	974	CALL wrd_write_string( 'pt_av' )
	975	WRITE ( 14 ) pt_av
	976	ENDIF
	977
[3936]	978	IF ( ALLOCATED( pt_2m_av ) ) THEN
	979	CALL wrd_write_string( 'pt_2m_av' )
	980	WRITE ( 14 ) pt_2m_av
	981	ENDIF
	982
[2894]	983	IF ( humidity ) THEN
	984
	985	CALL wrd_write_string( 'q' )
	986	WRITE ( 14 ) q
	987
	988	IF ( ALLOCATED( q_av ) ) THEN
	989	CALL wrd_write_string( 'q_av' )
	990	WRITE ( 14 ) q_av
	991	ENDIF
	992
[3274]	993	IF ( cloud_droplets ) THEN
[2894]	994
[3274]	995	CALL wrd_write_string( 'ql' )
[2894]	996	WRITE ( 14 ) ql
	997
	998	IF ( ALLOCATED( ql_av ) ) THEN
[3274]	999	CALL wrd_write_string( 'ql_av' )
[2894]	1000	WRITE ( 14 ) ql_av
	1001	ENDIF
	1002
	1003	ENDIF
	1004
	1005	IF ( ALLOCATED( qsws_av ) ) THEN
	1006	CALL wrd_write_string( 'qsws_av' )
	1007	WRITE ( 14 ) qsws_av
	1008	ENDIF
	1009
	1010	ENDIF
	1011
	1012	IF ( passive_scalar ) THEN
	1013
	1014	CALL wrd_write_string( 's' )
	1015	WRITE ( 14 ) s
	1016
	1017	IF ( ALLOCATED( s_av ) ) THEN
	1018	CALL wrd_write_string( 's_av' )
	1019	WRITE ( 14 ) s_av
	1020	ENDIF
	1021
	1022	IF ( ALLOCATED( ssws_av ) ) THEN
	1023	CALL wrd_write_string( 'ssws_av' )
	1024	WRITE ( 14 ) ssws_av
	1025	ENDIF
	1026
	1027	ENDIF
	1028
	1029	IF ( ALLOCATED( ql_c_av ) ) THEN
	1030	CALL wrd_write_string( 'ql_c_av' )
	1031	WRITE ( 14 ) ql_c_av
	1032	ENDIF
	1033
	1034	IF ( ALLOCATED( ql_v_av ) ) THEN
	1035	CALL wrd_write_string( 'ql_v_av' )
	1036	WRITE ( 14 ) ql_v_av
	1037	ENDIF
	1038
	1039	IF ( ALLOCATED( ql_vp_av ) ) THEN
	1040	CALL wrd_write_string( 'ql_vp_av' )
	1041	WRITE ( 14 ) ql_vp_av
	1042	ENDIF
	1043
	1044	IF ( ALLOCATED( qv_av ) ) THEN
	1045	CALL wrd_write_string( 'qv_av' )
	1046	WRITE ( 14 ) qv_av
	1047	ENDIF
	1048
	1049	CALL wrd_write_string( 'random_iv' )
	1050	WRITE ( 14 ) random_iv
	1051	WRITE ( 14 ) random_iy
	1052
	1053	IF ( ALLOCATED( seq_random_array ) ) THEN
	1054	CALL wrd_write_string( 'seq_random_array' )
	1055	WRITE ( 14 ) id_random_array
	1056	WRITE ( 14 ) seq_random_array
	1057	ENDIF
	1058
	1059	IF ( ALLOCATED( s_av ) ) THEN
	1060	CALL wrd_write_string( 's_av' )
	1061	WRITE ( 14 ) s_av
	1062	ENDIF
	1063
	1064	IF ( ALLOCATED( shf_av ) ) THEN
	1065	CALL wrd_write_string( 'shf_av' )
	1066	WRITE ( 14 ) shf_av
	1067	ENDIF
	1068
	1069	IF ( ALLOCATED( ts_av ) ) THEN
	1070	CALL wrd_write_string( 'ts_av' )
	1071	WRITE ( 14 ) ts_av
	1072	ENDIF
	1073
	1074	CALL wrd_write_string( 'u' )
	1075	WRITE ( 14 ) u
	1076
	1077	IF ( ALLOCATED( u_av ) ) THEN
	1078	CALL wrd_write_string( 'u_av' )
	1079	WRITE ( 14 ) u_av
	1080	ENDIF
	1081
	1082	IF ( ALLOCATED( u_m_l ) ) THEN
	1083	CALL wrd_write_string( 'u_m_l' )
	1084	WRITE ( 14 ) u_m_l
	1085	ENDIF
	1086
	1087	IF ( ALLOCATED( u_m_n ) ) THEN
	1088	CALL wrd_write_string( 'u_m_n' )
	1089	WRITE ( 14 ) u_m_n
	1090	ENDIF
	1091
	1092	IF ( ALLOCATED( u_m_r ) ) THEN
	1093	CALL wrd_write_string( 'u_m_r' )
	1094	WRITE ( 14 ) u_m_r
	1095	ENDIF
	1096
	1097	IF ( ALLOCATED( u_m_s ) ) THEN
	1098	CALL wrd_write_string( 'u_m_s' )
	1099	WRITE ( 14 ) u_m_s
	1100	ENDIF
	1101
	1102	IF ( ALLOCATED( us_av ) ) THEN
	1103	CALL wrd_write_string( 'us_av' )
	1104	WRITE ( 14 ) us_av
	1105	ENDIF
	1106
	1107	CALL wrd_write_string( 'v' )
	1108	WRITE ( 14 ) v
	1109
	1110	IF ( ALLOCATED( v_av ) ) THEN
	1111	CALL wrd_write_string( 'v_av' )
	1112	WRITE ( 14 ) v_av
	1113	ENDIF
	1114
	1115	IF ( ALLOCATED( v_m_l ) ) THEN
	1116	CALL wrd_write_string( 'v_m_l' )
	1117	WRITE ( 14 ) v_m_l
	1118	ENDIF
	1119
	1120	IF ( ALLOCATED( v_m_n ) ) THEN
	1121	CALL wrd_write_string( 'v_m_n' )
	1122	WRITE ( 14 ) v_m_n
	1123	ENDIF
	1124
	1125	IF ( ALLOCATED( v_m_r ) ) THEN
	1126	CALL wrd_write_string( 'v_m_r' )
	1127	WRITE ( 14 ) v_m_r
	1128	ENDIF
	1129
	1130	IF ( ALLOCATED( v_m_s ) ) THEN
	1131	CALL wrd_write_string( 'v_m_s' )
	1132	WRITE ( 14 ) v_m_s
	1133	ENDIF
	1134
	1135	IF ( humidity ) THEN
	1136
	1137	CALL wrd_write_string( 'vpt' )
	1138	WRITE ( 14 ) vpt
	1139
	1140	IF ( ALLOCATED( vpt_av ) ) THEN
	1141	CALL wrd_write_string( 'vpt_av' )
	1142	WRITE ( 14 ) vpt_av
	1143	ENDIF
	1144
	1145	ENDIF
	1146
	1147	CALL wrd_write_string( 'w' )
	1148	WRITE ( 14 ) w
	1149
	1150	IF ( ALLOCATED( w_av ) ) THEN
	1151	CALL wrd_write_string( 'w_av' )
	1152	WRITE ( 14 ) w_av
	1153	ENDIF
	1154
	1155	IF ( ALLOCATED( w_m_l ) ) THEN
	1156	CALL wrd_write_string( 'w_m_l' )
	1157	WRITE ( 14 ) w_m_l
	1158	ENDIF
	1159
	1160	IF ( ALLOCATED( w_m_n ) ) THEN
	1161	CALL wrd_write_string( 'w_m_n' )
	1162	WRITE ( 14 ) w_m_n
	1163	ENDIF
	1164
	1165	IF ( ALLOCATED( w_m_r ) ) THEN
	1166	CALL wrd_write_string( 'w_m_r' )
	1167	WRITE ( 14 ) w_m_r
	1168	ENDIF
	1169
	1170	IF ( ALLOCATED( w_m_s ) ) THEN
	1171	CALL wrd_write_string( 'w_m_s' )
	1172	WRITE ( 14 ) w_m_s
	1173	ENDIF
	1174
	1175	IF ( ALLOCATED( z0_av ) ) THEN
	1176	CALL wrd_write_string( 'z0_av' )
	1177	WRITE ( 14 ) z0_av
	1178	ENDIF
	1179
	1180	IF ( ALLOCATED( z0h_av ) ) THEN
	1181	CALL wrd_write_string( 'z0h_av' )
	1182	WRITE ( 14 ) z0h_av
	1183	ENDIF
	1184
	1185	IF ( ALLOCATED( z0q_av ) ) THEN
	1186	CALL wrd_write_string( 'z0q_av' )
	1187	WRITE ( 14 ) z0q_av
	1188	ENDIF
	1189
	1190	IF ( land_surface .OR. urban_surface ) THEN
	1191
	1192	IF ( ALLOCATED( ghf_av ) ) THEN
	1193	CALL wrd_write_string( 'ghf_av' )
	1194	WRITE ( 14 ) ghf_av
	1195	ENDIF
	1196
	1197	IF ( ALLOCATED( r_a_av ) ) THEN
	1198	CALL wrd_write_string( 'r_a_av' )
	1199	WRITE ( 14 ) r_a_av
	1200	ENDIF
	1201
	1202	ENDIF
	1203
	1204	IF ( ALLOCATED( tsurf_av ) ) THEN
	1205	CALL wrd_write_string( 'tsurf_av' )
	1206	WRITE ( 14 ) tsurf_av
	1207	ENDIF
	1208
[3637]	1209
[3582]	1210	CALL surface_wrd_local
[2894]	1211
	1212	!
[3637]	1213	!-- Write restart data of other modules
	1214	CALL module_interface_wrd_local
[2894]	1215
	1216	!
[3294]	1217	!-- Write end label
[2894]	1218	CALL wrd_write_string( '* end *' )
	1219
	1220	END SUBROUTINE wrd_local
	1221
	1222
	1223	END MODULE write_restart_data_mod

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

Context Navigation

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

Download in other formats: