Enforce domain radius = L1 when roche lobe overflow is activated

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

Enforce domain radius = L1 when roche lobe overflow is activated #52

Closed Lumajord closed 10 months ago

Lumajord commented 10 months ago

Roche lobe overflow is intended to be used when the outer domain radius is the L1 point. We should enforce this criteria to make the option more fail safe.