rometsch / fargocpt

Fluid dynamics code to study the interaction between protoplanetary disks and planets.
GNU Affero General Public License v3.0
6 stars 3 forks source link

Sg scale height #58

Closed rometsch closed 10 months ago

rometsch commented 10 months ago

When using self-gravity, the scaleheight should be adjusted. See 3.7 of the code paper, especially Eqs. (23)-(24). This is implemented here.

In case of SG, the Toomre parameter is now calculated every time the scale height is calculated. Watch out if this is a problem in your usecase.

Included is also a fix of the calculation of the smoothing length for direct summation in the sg test.