Code-equations/Gravitational-Potential-including-indirect-terms.md

@@ -41,8 +41,10 @@ $`d -> 0`$. The choice of the smoothing length depends mainly on the resolution:
 - In the 2 Dimensional case 
 we use  a potential of the form:
 $`\Phi _p^{\epsilon} = -{m_pG\over {\sqrt {d^2+(\epsilon R_H)^2}}}`$
-In 2 dimensional models we do not only have the problem of avoiding the singularity but we choose a potential that allows to mimic the average influence  that the planet would have  on the vertical gas column.
-A value often used for the smoothing is $`0.7H`$ where $`H`$ is the disc  scale height. We remark that the parameter $`\epsilon`$ is considered in the code in units of the Hill radius and not in disc scale height.
+
+In 2 dimensional models we certainly have the problem of avoiding the singularity as well as in the 3D case; however the main constraint is to choose a potential that allows to mimic the average influence  that the planet would have  on the vertical gas column.
+
+A value often used for the smoothing lenght is $`0.7H`$ where $`H`$ is the disc  scale height. We remark that the smoothing parameter $`\epsilon`$ is considered in the code in units of the Hill radius and not in disc scale height.
 To select this choice for a simulation in the config file:
 ```
     Smoothing