modified the sink_mass initial value to make it stable.
This commit is contained in:
@@ -54,7 +54,7 @@ AC_sink_pos_x = 3.14
|
|||||||
AC_sink_pos_y = 3.14
|
AC_sink_pos_y = 3.14
|
||||||
AC_sink_pos_z = 3.14
|
AC_sink_pos_z = 3.14
|
||||||
//AC_M_sink_Msun = 0.005
|
//AC_M_sink_Msun = 0.005
|
||||||
AC_M_sink_Msun = 1e-8
|
AC_M_sink_Msun = 1e-10
|
||||||
AC_soft = 0.36
|
AC_soft = 0.36
|
||||||
|
|
||||||
// Accretion Parameters
|
// Accretion Parameters
|
||||||
|
@@ -255,7 +255,7 @@ run_simulation(void)
|
|||||||
accreted_mass = accreted_mass + sum_mass;
|
accreted_mass = accreted_mass + sum_mass;
|
||||||
AcReal sink_mass = 0.0;
|
AcReal sink_mass = 0.0;
|
||||||
//if (i > 1000 ) {
|
//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("sink mass is: %e \n", sink_mass);
|
||||||
printf("accreted mass is: %e \n", accreted_mass);
|
printf("accreted mass is: %e \n", accreted_mass);
|
||||||
|
Reference in New Issue
Block a user