Module Temperature

Define a temperature power law.

The functional form is: Tg(r) = T0 * (R0/r)^epsilon

Parameters:

• r (array) - The radial grid in cm
• T0 (float) - The initial temperature at Ri in K
• r0 (float) - The inner radius in cm
• epsilon (float) - The exponent of the temp law

Returns: array

The temperature grid (K)

Return a dust temperature power law of the form as suggested by observational evidence.

See Thesis p32, where power is -2/(4+s) in T(r) = T_eff*(2*r/R_STAR)**(-2/(4+s)) and s is typically 1 in the Rayleigh limit and optically thin case, with T_eff = T0, R_STAR = r0.

Parameters:

• r (array) - The radial grid for the power law in Rstar
• T0 (float) - The stellar effective temperature
• r0 - The initial value for the radius, typically the stellar radius. @type r0
• s (int) - The s parameter in the power law T(r) given above.

  (default: 1)

Returns: array

The temperature grid (K)