Commit 7f238e06 authored by David Deepwell's avatar David Deepwell
Browse files

Explicitly set dt_max value

If dt_max is not given in the spins.conf, it is now explicitly
set to 0 so as to be adjusted later to use the buoyancy frequency.
The dt_max is also printed in the initialization for clarity.
parent 99219705
......@@ -446,7 +446,7 @@ int main(int argc, char ** argv) {
option_category("Temporal options");
add_option("final_time",&final_time,"Final time");
add_option("plot_interval",&plot_interval,"Time between writes");
add_option("dt_max",&dt_max,"Maximum time step");
add_option("dt_max",&dt_max,0,"Maximum time step. Zero value results in the default");
option_category("Restart options");
add_option("restart",&restarting,false,"Restart from prior output time.");
......@@ -523,7 +523,8 @@ int main(int argc, char ** argv) {
fprintf(stdout,"Time between plots: %g s\n",plot_interval);
fprintf(stdout,"Initial velocity perturbation: %g\n",perturb);
fprintf(stdout,"Filter cutoff = %f, order = %f, strength = %f\n",f_cutoff,f_order,f_strength);
fprintf(stdout,"Buoyancy frequency squared %g\n",N2_max);
fprintf(stdout,"Approx. max buoyancy frequency squared: %g\n",N2_max);
fprintf(stdout,"Max time step: %g\n",dt_max);
/* ------------------ Do stuff --------------------- */
Supports Markdown
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