From 9a1a6fef17ae81b09e1c34441991425ec2e6fa52 Mon Sep 17 00:00:00 2001 From: JackHsu Date: Mon, 19 Aug 2019 15:38:35 +0800 Subject: [PATCH] modified the sink_mass initial value to make it stable. --- config/templates/sinktest.conf | 2 +- src/standalone/simulation.cc | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/config/templates/sinktest.conf b/config/templates/sinktest.conf index 258e0a3..1088da3 100644 --- a/config/templates/sinktest.conf +++ b/config/templates/sinktest.conf @@ -54,7 +54,7 @@ AC_sink_pos_x = 3.14 AC_sink_pos_y = 3.14 AC_sink_pos_z = 3.14 //AC_M_sink_Msun = 0.005 -AC_M_sink_Msun = 1e-8 +AC_M_sink_Msun = 1e-10 AC_soft = 0.36 // Accretion Parameters diff --git a/src/standalone/simulation.cc b/src/standalone/simulation.cc index c883086..c72f1d9 100644 --- a/src/standalone/simulation.cc +++ b/src/standalone/simulation.cc @@ -255,7 +255,7 @@ run_simulation(void) accreted_mass = accreted_mass + sum_mass; AcReal sink_mass = 0.0; //if (i > 1000 ) { - sink_mass = AC_M_sink_init + accreted_mass; + sink_mass = mesh_info.real_params[AC_M_sink_init] + accreted_mass; //} printf("sink mass is: %e \n", sink_mass); printf("accreted mass is: %e \n", accreted_mass);