From 5cacda285076d3e41f59d83ceff42535b87a9eb1 Mon Sep 17 00:00:00 2001 From: Miikka Vaisala Date: Wed, 26 Jun 2019 13:15:28 +0800 Subject: [PATCH] Helical forcing funtion works. But we will need a wavenumber generator to add stochasticity. --- acc/mhd_solver/stencil_process.sps | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/acc/mhd_solver/stencil_process.sps b/acc/mhd_solver/stencil_process.sps index 9b6780d..714e742 100644 --- a/acc/mhd_solver/stencil_process.sps +++ b/acc/mhd_solver/stencil_process.sps @@ -263,7 +263,7 @@ forcing(int3 globalVertexIdx) //Determine that forcing funtion type at this point. //Vector force = simple_vortex_forcing(a, xx, magnitude); //Vector force = simple_outward_flow_forcing(a, xx, magnitude); - Vector force = helical_forcing(magnitude, k_force, xx, ff_re,ff_im, phi); + Vector force = helical_forcing(magnitude, k_force, xx, ff_re,ff_im, phase); if (is_valid(force)) { return force; } else { return (Vector){0, 0, 0}; }