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

Changed to step function

Browse Source
This commit is contained in:
Miikka Vaisala
2019-08-16 11:44:41 +08:00
parent a8fbced89d
commit f921840505
1 changed files with 11 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
acc/mhd_solver/stencil_process.sps
View File

@@ -105,9 +105,18 @@ sink_accretion(int3 globalVertexIdx, in ScalarField lnrho, Scalar dt){
// return accretion_density;
// multiplying the truelove density by a wave function to avoid step-function like accretion profile.
const Scalar weight = exp(-(accretion_distance/profile_range));
Scalar weight;
// const Scalar weight = exp(-(accretion_distance/profile_range));
// Step function weighting
if ((accretion_distance) <= profile_range){
weight = Scalar(1.0);
} else {
weight = Scalar(0.0);
}
const Scalar sink_mass = DCONST_REAL(AC_M_sink);
const Scalar lnrho_min = Scalar(-10.0);
const Scalar lnrho_min = Scalar(-10.0); //TODO Define from astaroth.conf
// const Scalar rate = truelove_density(lnrho);
// const Scalar B = Scalar(0.5);
// const Scalar k = Scalar(1.5);