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

Flexible boundary conditions called from simulation.cc

Browse Source
This commit is contained in:
Miikka Vaisala
2020-11-20 16:31:13 +08:00
parent 4add619e2f
commit 288693fab5
1 changed files with 15 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
samples/standalone/simulation.cc
View File

@@ -367,7 +367,8 @@ run_simulation(const char* config_path)
#endif #endif
} }
acBoundcondStep(); //acBoundcondStep();
acBoundcondStepGBC(mesh_info);
acStore(mesh); acStore(mesh);
if (start_step == 0) { if (start_step == 0) {
save_mesh(*mesh, 0, t_step); save_mesh(*mesh, 0, t_step);
@@ -426,7 +427,11 @@ run_simulation(const char* config_path)
loadForcingParamsToDevice(forcing_params); loadForcingParamsToDevice(forcing_params);
#endif #endif
acIntegrate(dt); /* Uses now flexible bokundary conditions */
//acIntegrate(dt);
acIntegrateGBC(mesh_info, dt);
t_step += dt; t_step += dt;
@@ -472,7 +477,8 @@ run_simulation(const char* config_path)
acBoundcondStep(); acBoundcondStep();
acStore(mesh); acStore(mesh);
*/ */
acBoundcondStep(); //acBoundcondStep();
acBoundcondStepGBC(mesh_info);
acStore(mesh); acStore(mesh);
save_mesh(*mesh, i, t_step); save_mesh(*mesh, i, t_step);
@@ -493,7 +499,8 @@ run_simulation(const char* config_path)
if (dt < dt_typical/AcReal(1e5)) { if (dt < dt_typical/AcReal(1e5)) {
if (dtcounter > 10) { if (dtcounter > 10) {
printf("dt = %e TOO LOW! Ending run at t = %#e \n", double(dt), double(t_step)); printf("dt = %e TOO LOW! Ending run at t = %#e \n", double(dt), double(t_step));
acBoundcondStep(); //acBoundcondStep();
acBoundcondStepGBC(mesh_info);
acStore(mesh); acStore(mesh);
save_mesh(*mesh, i, t_step); save_mesh(*mesh, i, t_step);
break; break;
@@ -507,7 +514,8 @@ run_simulation(const char* config_path)
// End loop if nan is found // End loop if nan is found
if (found_nan > 0) { if (found_nan > 0) {
printf("Found nan at t = %e \n", double(t_step)); printf("Found nan at t = %e \n", double(t_step));
acBoundcondStep(); //acBoundcondStep();
acBoundcondStepGBC(mesh_info);
acStore(mesh); acStore(mesh);
save_mesh(*mesh, i, t_step); save_mesh(*mesh, i, t_step);
break; break;
@@ -522,7 +530,8 @@ run_simulation(const char* config_path)
if (found_stop == 1) { if (found_stop == 1) {
printf("Found STOP file at t = %e \n", double(t_step)); printf("Found STOP file at t = %e \n", double(t_step));
acBoundcondStep(); //acBoundcondStep();
acBoundcondStepGBC(mesh_info);
acStore(mesh); acStore(mesh);
save_mesh(*mesh, i, t_step); save_mesh(*mesh, i, t_step);
break; break;