From 0587d08b01f7e2ee6bfad544b361ca97bc35c96b Mon Sep 17 00:00:00 2001
From: JackHsu <ssp201903@asiaa.sinica.edu.tw>
Date: Tue, 13 Aug 2019 18:31:57 +0800
Subject: [PATCH] Had some rough draft of accretion criterion.

---
 acc/mhd_solver/stencil_process.sps | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/acc/mhd_solver/stencil_process.sps b/acc/mhd_solver/stencil_process.sps
index 2e0b2fe..7074a4b 100644
--- a/acc/mhd_solver/stencil_process.sps
+++ b/acc/mhd_solver/stencil_process.sps
@@ -118,7 +118,10 @@ accretion_profile(int3 globalVertexIdx, in Scalar lnrho){
      
     // multiplying the truelove density by a wave function to avoid step-function like accretion profile.   
     const Scalar weight = exp(-(accretion_distance/profile_range));    
-    const Scalar rate = truelove_density(lnrho);
+//    const Scalar rate = truelove_density(lnrho);
+    const Scalar B = Scalar(3.0);
+    const Scalar k = Scalar(2.0);
+    const Scalar rate = B * (pow(sink_mass, k) / (dsx * dsy * dsz));
     accretion_density = weight * rate;
     
     return accretion_density;