From 6927ddff5ae77264de02d0de8eb8c0a4c9331b61 Mon Sep 17 00:00:00 2001 From: Miikka Vaisala Date: Thu, 1 Aug 2019 17:57:21 +0800 Subject: [PATCH] Correct an obscure bug. --- acc/mhd_solver/stencil_defines.h | 1 + src/standalone/config_loader.cc | 10 +++++----- 2 files changed, 6 insertions(+), 5 deletions(-) diff --git a/acc/mhd_solver/stencil_defines.h b/acc/mhd_solver/stencil_defines.h index ed67f5f..dbd632e 100644 --- a/acc/mhd_solver/stencil_defines.h +++ b/acc/mhd_solver/stencil_defines.h @@ -78,6 +78,7 @@ FUNC(AC_sink_pos_y),\ FUNC(AC_sink_pos_z),\ FUNC(AC_M_sink),\ + FUNC(AC_M_sink_Msun),\ FUNC(AC_soft),\ /* Run params */\ FUNC(AC_cdt), \ diff --git a/src/standalone/config_loader.cc b/src/standalone/config_loader.cc index f3d7aa3..abda35a 100644 --- a/src/standalone/config_loader.cc +++ b/src/standalone/config_loader.cc @@ -137,11 +137,11 @@ update_config(AcMeshInfo* config) 6.674e-8); // g/cm3/s GGS definition //TODO define in a separate module AcReal M_sun = AcReal(1.989e33); // g solar mass - config->real_params[AC_M_sink] = config->real_params[AC_M_sink] * M_sun / - ((config->real_params[AC_unit_length] * - config->real_params[AC_unit_length] * - config->real_params[AC_unit_length]) * - config->real_params[AC_unit_density]); + config->real_params[AC_M_sink_Msun] = config->real_params[AC_M_sink] * M_sun / + ((config->real_params[AC_unit_length] * + config->real_params[AC_unit_length] * + config->real_params[AC_unit_length]) * + config->real_params[AC_unit_density]); config->real_params[AC_G_const] = G_CONST_CGS / ((config->real_params[AC_unit_velocity] * config->real_params[AC_unit_velocity]) /