Runtime parameters ¶

Output steering ¶

Run steering ¶

Output steering: ¶

Parameter Name	FORTRAN Type	Default Value	Explanation
averaging_interval ¶	R	0.0	Averaging interval for all output of temporally averaged data (in s). This parameter defines the time interval length for temporally averaged data (vertical profiles, spectra, 2d cross-sections, 3d volume data). By default, data are not subject to temporal averaging. The interval length is limited by the parameter dt_data_output_av. In any case, averaging_interval <= dt_data_output_av must hold. If an interval is defined, then by default the average is calculated from the data values of all timesteps lying within this interval. The number of time levels entering into the average can be reduced with the parameter dt_averaging_input?. If an averaging interval can not be completed at the end of a run, it will be finished at the beginning of the next restart run. Thus for restart runs, averaging intervals do not necessarily begin at the beginning of the run. Parameters averaging_interval_pr and averaging_interval_sp? can be used to define different averaging intervals for vertical profile data and spectra, respectively.
averaging_interval_pr ¶	R	value of averaging_interval	Averaging interval for output of vertical profiles to local file DATA_1D_PR_NETCDF? (in s). If this parameter is given a non-zero value, temporally averaged vertical profile data are output. By default, profile data data are not subject to temporal averaging. The interval length is limited by the parameter dt_dopr?. In any case averaging_interval_pr <= dt_dopr must hold. If an interval is defined, then by default the average is calculated from the data values of all timesteps lying within this interval. The number of time levels entering into the average can be reduced with the parameter dt_averaging_input_pr?. If an averaging interval can not be completed at the end of a run, it will be finished at the beginning of the next restart run. Thus for restart runs, averaging intervals do not necessarily begin at the beginning of the run.
dt_data_output_av ¶	R	0.0

Run steering: ¶

Parameter Name	FORTRAN Type	Default Value	Explanation
call_psolver_at_all_substeps ¶	L	.T.	Switch to steer the call of the pressure solver. In order to speed-up performance, the Poisson equation for perturbation pressure (see psolver?) can be called only at the last substep of multistep Runge-Kutta timestep schemes (see timestep_scheme?) by setting call_psolver_at_all_substeps = .F.. In many cases, this sufficiently reduces the divergence of the velocity field. Nevertheless, small-scale ripples (2-delta-x) may occur. In this case and in case of non-cyclic lateral boundary conditions, call_psolver_at_all_timesteps = .T. should be used.
cfl_factor ¶	R	0.1, 0.8 or 0.9 (see right)	Time step limiting factor. In the model, the maximum allowed time step according to CFL and diffusion-criterion dt_max is reduced by dt? = dt_max * cfl_factor in order to avoid stability problems which may arise in the vicinity of the maximum allowed timestep. The condition 0.0 < cfl_factor < 1.0 applies. The default value of cfl_factor depends on the timestep_scheme? used: For the third order Runge-Kutta scheme it is cfl_factor = 0.9. In case of the leapfrog scheme a quite restrictive value of cfl_factor = 0.1 is used because for larger values the velocity divergence significantly effects the accuracy of the model results. Possibly larger values may be used with the leapfrog scheme but these are to be determined by appropriate test runs. The default value for the Euler scheme is cfl_factor = 0.8 .
create_disturbances ¶	L	.T.	Switch to impose random perturbations to the horizontal velocity field. With create_disturbances = .T., random perturbations can be imposed to the horizontal velocity field at certain times e.g. in order to trigger off the onset of convection, etc.. The temporal interval between these times can be steered with dt_disturb?, the vertical range of the perturbations with disturbance_level_b? and disturbance_level_t?, and the perturbation amplitude with disturbance_amplitude?. In case of non-cyclic lateral boundary conditions (see bc_lr? and bc_ns?), the horizontal range of the perturbations is determined by inflow_disturbance_begin? and inflow_disturbance_end?. A perturbation is added to each grid point with its individual value determined by multiplying the disturbance amplitude with a uniformly distributed random number. After this, the arrays of u and v are smoothed by applying a Shuman-filter twice and made divergence-free by applying the pressure solver. The random number generator to be used can be chosen with random_generator?. As soon as the desired flow features have developed (e.g. convection has started), further imposing of perturbations is not necessary and can be omitted (this does not hold for non-cyclic lateral boundaries!). This can be steered by assigning an upper limit value for the perturbation energy (the perturbation energy is defined by the deviation of the velocity from the mean flow) using the parameter disturbance_energy_limit?. As soon as the perturbation energy has exceeded this energy limit, no more random perturbations are assigned. Timesteps where a random perturbation has been imposed are marked in the local file RUN_CONTROL? by the character "D" appended to the values of the maximum horizontal velocities.
call_psolver_at_all_substeps ¶	L	.T.
call_psolver_at_all_substeps ¶	L	.T.
call_psolver_at_all_substeps ¶	L	.T.