cwpearson/astaroth
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Tested and works. Now it is possible to continue the simulatiom with a specific file number.

Browse Source 
Next task: move I/O to src/utils/
This commit is contained in:
Miikka Vaisala
2019-10-02 14:30:13 +08:00
parent 0dbbcd22d5
commit 79fe634a84
1 changed files with 8 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/standalone/simulation.cc
View File

@@ -268,15 +268,20 @@ run_simulation(const char* config_path)
fprintf(diag_file, "\n"); fprintf(diag_file, "\n");
write_mesh_info(&mesh_info); write_mesh_info(&mesh_info);
if (start_step == 0) {
#if LSINK #if LSINK
print_diagnostics(0, AcReal(.0), t_step, diag_file, mesh_info.real_params[AC_M_sink_init], 0.0); print_diagnostics(0, AcReal(.0), t_step, diag_file, mesh_info.real_params[AC_M_sink_init], 0.0);
#else #else
print_diagnostics(0, AcReal(.0), t_step, diag_file, -1.0, -1.0); print_diagnostics(0, AcReal(.0), t_step, diag_file, -1.0, -1.0);
#endif #endif
}
acBoundcondStep(); acBoundcondStep();
acStore(mesh); acStore(mesh);
if (start_step == 0) {
save_mesh(*mesh, 0, t_step); save_mesh(*mesh, 0, t_step);
}
const int max_steps = mesh_info.int_params[AC_max_steps]; const int max_steps = mesh_info.int_params[AC_max_steps];
const int save_steps = mesh_info.int_params[AC_save_steps]; const int save_steps = mesh_info.int_params[AC_save_steps];