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

Helical forcing funtion works.

Browse Source 
But we will need a wavenumber generator to add stochasticity.
This commit is contained in:
Miikka Vaisala
2019-06-26 13:15:28 +08:00
parent 8191c47fa0
commit 5cacda2850
1 changed files with 1 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
acc/mhd_solver/stencil_process.sps
View File

@@ -263,7 +263,7 @@ forcing(int3 globalVertexIdx)
//Determine that forcing funtion type at this point. //Determine that forcing funtion type at this point.
//Vector force = simple_vortex_forcing(a, xx, magnitude); //Vector force = simple_vortex_forcing(a, xx, magnitude);
//Vector force = simple_outward_flow_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; } if (is_valid(force)) { return force; }
else { return (Vector){0, 0, 0}; } else { return (Vector){0, 0, 0}; }