Rosenfeld2d

class sf3dmodels.model.disc2d.Rosenfeld2d(grid)

Bases: sf3dmodels.model.disc2d.Velocity, sf3dmodels.model.disc2d.Intensity, sf3dmodels.model.disc2d.Linewidth, sf3dmodels.model.disc2d.Tools

Host class for the Rosenfeld+2013 model which describes the velocity field of a flared disc in 2D. This model assumes a (Keplerian) double cone to account for the near and far sides of the disc and solves analytical equations to find the line-of-sight velocity v_obs projected on the sky-plane from both sides.

Parameters:

grid : array_like, shape (nrows, ncols)

(x’, y’) map of the sky-plane onto which the disc velocity field will be projected.

Attributes:

velocity_func : function(coord, **kwargs)

Velocity function describing the kinematics of the disc. The argument coord is a dictionary of coordinates (e.g. ‘x’, ‘y’, ‘z’, ‘r’, ‘R’, ‘theta’, ‘phi’) where the function will be evaluated. Additional arguments are optional and depend upon the function definition, e.g. Mstar=1.0*sfu.Msun

Methods Summary

make_model(self, incl, psi[, PA, …])

Executes the Rosenfeld+2013 model.

Methods Documentation

make_model(self, incl, psi, PA=0.0, int_kwargs={}, vel_kwargs={}, lw_kwargs=None, ls_kwargs=None)

Executes the Rosenfeld+2013 model. The sum of incl+psi must be < 90, otherwise the quadratic equation will have imaginary roots as some portions of the cone (which has finite extent) do not intersect with the sky plane.

Parameters:

incl : scalar

Inclination of the disc midplane with respect to the x’y’ plane; pi/2 radians is edge-on.

psi : scalar

Opening angle of the cone describing the velocity field of the gas emitting layer in the disc; 0 radians returns the projected velocity field of the disc midplane (i.e no conic emission).

PA : scalar, optional

Position angle in radians. Measured from North (+y) to East (-x).

Attributes:

velocity : array_like, size (n,)

Velocity field computed using the Rosenfeld+2013 model.

velocity2d : array_like, size (nx, ny)

If set get_2d=True: Velocity field computed using the Rosenfeld+2013 model, reshaped to 2D to facilitate plotting.