set_binary

Module

Description

This module is contains utilities for setting up binaries

Our conventions for binary orbital parameters are consistent with those produced by the imorbel code (Pearce, Wyatt & Kennedy 2015) which can be used to produce orbits matching observed orbital arcs of binary companions on the sky

References:

Eggleton (1983) ApJ 268, 368-369 (ref:eggleton83) Lucy (2014), A&A 563, A126 Pearce, Wyatt & Kennedy (2015), MNRAS 448, 3679 https://en.wikipedia.org/wiki/Orbital_elements

Owner:

Daniel Price

Runtime parameters:

None

Dependencies:

orbits

Quick access

Variables:

get_mean_angmom_vector, ierr_ecc, ierr_hier1, ierr_hier2, ierr_m1, ierr_m2, ierr_missstar, ierr_omegasubst, ierr_semi, ierr_subststar, l1_point, rochelobe_estimate, rotate

Routines:

set_binary()

Needed modules

  • orbits (rochelobe_estimate(), l1_point(), get_mean_angmom_vector(), get_inclination(), get_e(), get_e_from_mean_anomaly(), get_e_from_true_anomaly(), pi(), deg_to_rad(), rad_to_deg(), orbit_is_parabolic(), get_specific_energy())

Variables

  • setbinary/get_mean_angmom_vector [public]
  • setbinary/ierr_ecc [integer,private/parameter/optional/default=3]
  • setbinary/ierr_hier1 [integer,private/parameter/optional/default=5]
  • setbinary/ierr_hier2 [integer,private/parameter/optional/default=6]
  • setbinary/ierr_m1 [integer,private/parameter/optional/default=1]
  • setbinary/ierr_m2 [integer,private/parameter/optional/default=2]
  • setbinary/ierr_missstar [integer,private/parameter/optional/default=9]
  • setbinary/ierr_omegasubst [integer,private/parameter/optional/default=8]
  • setbinary/ierr_semi [integer,private/parameter/optional/default=4]
  • setbinary/ierr_subststar [integer,private/parameter/optional/default=7]
  • setbinary/l1_point [public]
  • setbinary/rochelobe_estimate [public]
  • setbinary/rotate [private]

Subroutines and functions

subroutine  setbinary/set_binary(m1, m2, semimajoraxis, eccentricity, accretion_radius1, accretion_radius2, xyzmh_ptmass, vxyz_ptmass, nptmass, ierr[, omega_corotate[, posang_ascnode[, arg_peri[, incl[, f[, mean_anomaly[, verbose]]]]]]])
Parameters:

  • m1 [real,in]

  • m2 [real,in]

  • semimajoraxis [real,in]

  • eccentricity [real,in]

  • accretion_radius1 [real,in]

  • accretion_radius2 [real,in]

  • xyzmh_ptmass (,) [real,inout]

  • vxyz_ptmass (,) [real,inout]

  • nptmass [integer,inout]

  • ierr [integer,out]

  • omega_corotate [real,out,]

  • posang_ascnode [real,in,]

  • arg_peri [real,in,]

  • incl [real,in,]

  • f [real,in,]

  • mean_anomaly [real,in,]

  • verbose [logical,in,]