Commit 11ffdf27 authored by Stuart McAlpine's avatar Stuart McAlpine Committed by Matthieu Schaller
Browse files

Special case for when you have an outputlist entry on the final step

parent e0d70e55
......@@ -22,7 +22,7 @@ TimeIntegration:
# Parameters for the self-gravity scheme
Gravity:
eta: 0.025
MAC: adpative
MAC: adaptive
theta_cr: 0.7
epsilon_fmm: 0.001
comoving_DM_softening: 0.0889 # 1/25th of the mean inter-particle separation: 88.9 kpc
......@@ -35,7 +35,7 @@ Snapshots:
delta_time: 1.0816 # Only every second VELOCIraptor invoke gets a full snapshot dump.
scale_factor_first: 0.1 # z = 9
compression: 4
# Parameters governing the conserved quantities statistics
Statistics:
delta_time: 1.02
......
......@@ -22,7 +22,7 @@ TimeIntegration:
# Parameters for the self-gravity scheme
Gravity:
eta: 0.025
MAC: adpative
MAC: adaptive
theta_cr: 0.7
epsilon_fmm: 0.001
comoving_DM_softening: 0.0889 # 1/25th of the mean inter-particle separation: 88.9 kpc
......
......@@ -22,7 +22,7 @@ TimeIntegration:
# Parameters for the self-gravity scheme
Gravity:
eta: 0.025
MAC: adpative
MAC: adaptive
theta_cr: 0.7
epsilon_fmm: 0.001
comoving_DM_softening: 0.0889 # 1/25th of the mean inter-particle separation: 88.9 kpc
......
......@@ -22,7 +22,7 @@ TimeIntegration:
# Parameters for the self-gravity scheme
Gravity:
eta: 0.025
MAC: adpative
MAC: adaptive
theta_cr: 0.7
epsilon_fmm: 0.001
comoving_DM_softening: 0.0889 # 1/25th of the mean inter-particle separation: 88.9 kpc
......
......@@ -22,7 +22,7 @@ TimeIntegration:
# Parameters for the self-gravity scheme
Gravity:
eta: 0.025
MAC: adpative
MAC: adaptive
theta_cr: 0.7
epsilon_fmm: 0.001
comoving_DM_softening: 0.0889 # 1/25th of the mean inter-particle separation: 88.9 kpc
......
......@@ -243,7 +243,8 @@ void output_list_read_next_time(struct output_list *t, const struct engine *e,
* or a=1 is found in output_list.txt set the flag `final_step_dump`
* to 1 - this is not special behaviour that is controlled by a
* parameter file flag. */
if (time == time_end) {
if (time == time_end ||
(time > time_end && time - time_end < OUTPUT_LIST_EPS_TIME_END)) {
t->final_step_dump = 1;
if (e->verbose) {
if (is_cosmo) {
......
......@@ -33,6 +33,9 @@ struct engine;
/*! Maximal number of output lists */
#define OUTPUT_LIST_MAX_NUM_OF_SELECT_OUTPUT_STYLES 8
/*! How far past time_end do we check for an output time? */
#define OUTPUT_LIST_EPS_TIME_END 1e-4
/**
* @brief the different output_list type
*/
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment