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

Tested and works. We can now continue simulation from the chose snapshot number.

Browse Source
This commit is contained in:
Miikka Vaisala
2019-10-02 14:09:47 +08:00
parent 1b0e9803b0
commit 0dbbcd22d5
1 changed files with 10 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
src/standalone/simulation.cc
View File

@@ -168,7 +168,10 @@ read_mesh(AcMesh& read_mesh, const int step, AcReal* t_step)
for (size_t i = 0; i < n; ++i) { for (size_t i = 0; i < n; ++i) {
result = fread(&read_buf, sizeof(AcReal), 1, read_ptr); result = fread(&read_buf, sizeof(AcReal), 1, read_ptr);
read_mesh.vertex_buffer[VertexBufferHandle(w)][i] = read_buf; read_mesh.vertex_buffer[VertexBufferHandle(w)][i] = read_buf;
if (result != sizeof(AcReal)) {fprintf(stderr, "Reading error in %s, element %i\n", vtxbuf_names[w], int(i));} if (int(result) != 1) {
fprintf(stderr, "Reading error in %s, element %i\n", vtxbuf_names[w], int(i));
fprintf(stderr, "Result = %i, \n", int(result));
}
} }
fclose(read_ptr); fclose(read_ptr);
} }
@@ -252,11 +255,13 @@ run_simulation(const char* config_path)
diag_file = fopen("timeseries.ts", "a"); diag_file = fopen("timeseries.ts", "a");
// Generate the title row. // Generate the title row.
if (start_step == 0) {
fprintf(diag_file, "step t_step dt uu_total_min uu_total_rms uu_total_max "); fprintf(diag_file, "step t_step dt uu_total_min uu_total_rms uu_total_max ");
for (int i = 0; i < NUM_VTXBUF_HANDLES; ++i) { for (int i = 0; i < NUM_VTXBUF_HANDLES; ++i) {
fprintf(diag_file, "%s_min %s_rms %s_max ", vtxbuf_names[i], vtxbuf_names[i], fprintf(diag_file, "%s_min %s_rms %s_max ", vtxbuf_names[i], vtxbuf_names[i],
vtxbuf_names[i]); vtxbuf_names[i]);
} }
}
#if LSINK #if LSINK
fprintf(diag_file, "sink_mass accreted_mass "); fprintf(diag_file, "sink_mass accreted_mass ");
#endif #endif