duffell2020

dustpylib.substructures.gaps.duffell2020(r, a, q, h, alpha0)

Function calculates the planetary gap profile according Duffell (2020).

Parameters:

• r (array-like, (Nr,)) – Radial grid

• a (float) – Semi-major axis of planet

• q (float) – Planet-star mass ratio

• h (float) – Aspect ratio at planet location

• alpha0 (float) – Unperturbed alpha viscosity parameter

Returns:

f – Pertubation of surface density due to planet

Return type:

array-like, (Nr,)