pointmasses

Module

Description

$Id: pointmasses.f90,v 1.1 2019/02/02 03:54:41 wlyra Exp $

This module takes care of direct N-body gravity between point masses.

** AUTOMATIC CPARAM.INC GENERATION ************************

Declare (for generation of cparam.inc) the number of f array variables and auxiliary variables added by this module

MQVAR CONTRIBUTION 7 MQAUX CONTRIBUTION 0 MPVAR CONTRIBUTION 3 CPARAM logical, parameter :: lpointmasses=.true.

---

Quick access

Types:

indexdustparticles

Variables:

accrete_hills_frac, bcqx, bcqy, bcqz, cdtq, dfq, dfq_cart, eccentricity, final_ramped_mass, fq, frac_smooth, gnewton1, hills_tempering_fraction, imass, indexdustparticles, initvvq, initxxq, ipotential_pointmass, iprimary, isecondary, ivpx_cart, ivpy_cart, ivpz_cart, ivxq, ivyq, ivzq, ixq, iyq, izq, l2d, l3d, laccretion, ladd_dragforce, lcentrifugal_force_pointmass, lcoriolis_force_pointmass, lcorrect_gasgravity_lstart, lcylindrical_gravity_nbody, ldt_pointmasses, ldust, ldust_gravity, lexclude_frozen, lexclude_hills, lfollow_particle, lgas_gravity, lgas_removal, llinear_drag, llive_secondary, lmomentum_removal, lnogravz_star, lnorm, lnoselfgrav_primary, lquadratic_drag, lramp, lreset_cm, lretrograde, lshear_accel_pointmass, nqvar, pmass, pmass1, qvarname, r1_smooth, r_smooth, ramp_orbits, semimajor_axis, st1, stokesnumber, tau_accretion, totmass, totmass1, ugas, vxq0, vyq0, vzq0, xq0, yq0, zq0

Routines:

boundconds_pointmasses(), calc_diagnostics_pointmasses(), calc_pencils_pointmasses(), init_pointmasses(), initialize_pointmasses(), pencil_criteria_pointmasses(), pencil_interdep_pointmasses(), pointmasses_pde(), pointmasses_pde_pencil(), pointmasses_read_snapshot(), pointmasses_timestep_first(), pointmasses_timestep_second(), pointmasses_write_qdim(), pointmasses_write_snapshot(), read_pointmasses_init_pars(), read_pointmasses_run_pars(), register_pointmasses(), rprint_pointmasses(), write_pointmasses_init_pars(), write_pointmasses_run_pars()

Needed modules

• cdata: $Id$

• messages: $Id$

• general (keep_compiler_quiet(), find_proc()): $Id$

Types

• type  pointmasses/indexdustparticles

  Type fields:

  • % acc_grav (nqpar,3) [real]

  • % idiag_mdot_pt [integer,optional/default=0]

  • % idiag_period (nqpar) [integer,optional/default=0]

  • % idiag_torqext (nqpar) [integer,optional/default=0]

  • % idiag_torqext_gas (nqpar) [integer,optional/default=0]

  • % idiag_torqext_par (nqpar) [integer,optional/default=0]

  • % idiag_torqint (nqpar) [integer,optional/default=0]

  • % idiag_torqint_gas (nqpar) [integer,optional/default=0]

  • % idiag_torqint_par (nqpar) [integer,optional/default=0]

  • % idiag_torque (nqpar) [integer,optional/default=0]

  • % idiag_totenergy [integer,optional/default=0]

  • % idiag_vvq (nqpar,3) [integer,optional/default=0]

  • % idiag_xxq (nqpar,3) [integer,optional/default=0]

  • % index [indexdustparticles,save]

  • % ivxw [integer,optional/default=0]

  • % ivyw [integer,optional/default=0]

  • % ivzw [integer,optional/default=0]

  • % ixw [integer,optional/default=0]

  • % iyw [integer,optional/default=0]

  • % izw [integer,optional/default=0]

Variables

• pointmasses/accrete_hills_frac (nqpar) [real,private/optional/default=0.2]

• pointmasses/bcqx [character,private/optional/default='p']

• pointmasses/bcqy [character,private/optional/default='p']

• pointmasses/bcqz [character,private/optional/default='p']

• pointmasses/cdtq [real,private/optional/default=0.1]

• pointmasses/density_scale [real,private/optional/default=0.001]

• pointmasses/dfq (nqpar,mqvar) [real,private/optional/default=0.0]

• pointmasses/dfq_cart (nqpar,3) [real,private/optional/default=0.0]

• pointmasses/eccentricity [real,private/optional/default=0.0]

• pointmasses/final_ramped_mass (nqpar) [real,private/optional/default=0.0]

• pointmasses/fq (nqpar,mqarray) [real,private/optional/default=0.0]

• pointmasses/frac_smooth (nqpar) [real,private/optional/default=0.4]

• pointmasses/gnewton [real,private/optional/default=impossible]

• pointmasses/gnewton1 [real,private]

• pointmasses/hills_tempering_fraction [real,private/optional/default=0.8]

• pointmasses/imass [integer,private/optional/default=0]

• pointmasses/indexdustparticles [private]

• pointmasses/initvvq [character,private/optional/default='nothing']

• pointmasses/initxxq [character,private/optional/default='random']

• pointmasses/ipotential_pointmass (nqpar) [character,private/optional/default='newton']

• pointmasses/iprimary [integer,private/optional/default=1]

• pointmasses/isecondary [integer,private/optional/default=2]

• pointmasses/ivpx_cart [integer,private/optional/default=0]

• pointmasses/ivpy_cart [integer,private/optional/default=0]

• pointmasses/ivpz_cart [integer,private/optional/default=0]

• pointmasses/ivxq [integer,private/optional/default=0]

• pointmasses/ivyq [integer,private/optional/default=0]

• pointmasses/ivzq [integer,private/optional/default=0]

• pointmasses/ixq [integer,private/optional/default=0]

• pointmasses/iyq [integer,private/optional/default=0]

• pointmasses/izq [integer,private/optional/default=0]

• pointmasses/l2d [logical,private]

• pointmasses/l3d [logical,private]

• pointmasses/laccretion (nqpar) [logical,private/optional/default=.false.]

• pointmasses/ladd_dragforce [logical,private/optional/default=.false.]

• pointmasses/lcentrifugal_force_pointmass [logical,private/optional/default=.false.]

• pointmasses/lcoriolis_force_pointmass [logical,private/optional/default=.true.]

• pointmasses/lcorrect_gasgravity_lstart [logical,private/optional/default=.false.]

• pointmasses/lcylindrical_gravity_nbody (nqpar) [logical,private/optional/default=.false.]

• pointmasses/ldt_pointmasses [logical,private/optional/default=.true.]

• pointmasses/ldust [logical,private/optional/default=.false.]

• pointmasses/ldust_gravity [logical,private/optional/default=.false.]

• pointmasses/lexclude_frozen [logical,private/optional/default=.false.]

• pointmasses/lexclude_hills [logical,private/optional/default=.false.]

• pointmasses/lfollow_particle (nqpar) [logical,private/optional/default=.false.]

• pointmasses/lgas_gravity [logical,private/optional/default=.true.]

• pointmasses/lgas_removal [logical,private/optional/default=.false.]

• pointmasses/llinear_drag [logical,private/optional/default=.true.]

• pointmasses/llive_secondary [logical,private/optional/default=.false.]

• pointmasses/lmomentum_removal [logical,private/optional/default=.false.]

• pointmasses/lnogravz_star [logical,private/optional/default=.false.]

• pointmasses/lnorm [logical,private/optional/default=.true.]

• pointmasses/lnoselfgrav_primary [logical,private/optional/default=.true.]

• pointmasses/lquadratic_drag [logical,private/optional/default=.false.]

• pointmasses/lramp [logical,private/optional/default=.false.]

• pointmasses/lreset_cm [logical,private/optional/default=.true.]

• pointmasses/lretrograde [logical,private/optional/default=.false.]

• pointmasses/lshear_accel_pointmass [logical,private/optional/default=.true.]

• pointmasses/nqvar [integer,private/optional/default=0]

• pointmasses/pmass (nqpar) [real,private/optional/default=0.0]

• pointmasses/pmass1 (nqpar) [real,private]

• pointmasses/qvarname (mqarray) [character,private]

• pointmasses/r1_smooth (nqpar) [real,private]

• pointmasses/r_smooth (nqpar) [real,private/optional/default=impossible]

• pointmasses/ramp_orbits [integer,private/optional/default=5]

• pointmasses/rhs_poisson_const [real,private/pointer]

• pointmasses/semimajor_axis [real,private/optional/default=1.0]

• pointmasses/st1 (nqpar) [real,private]

• pointmasses/stokesnumber (nqpar) [real,private/optional/default=1.]

• pointmasses/tau_accretion [real,private/optional/default=1.0]

• pointmasses/totmass [real,private]

• pointmasses/totmass1 [real,private]

• pointmasses/tstart_selfgrav [real,private/pointer]

• pointmasses/ugas [real,private/optional/default=0.0]

• pointmasses/vxq0 (nqpar) [real,private/optional/default=0.0]

• pointmasses/vyq0 (nqpar) [real,private/optional/default=0.0]

• pointmasses/vzq0 (nqpar) [real,private/optional/default=0.0]

• pointmasses/xq0 (nqpar) [real,private/optional/default=0.0]

• pointmasses/yq0 (nqpar) [real,private/optional/default=0.0]

• pointmasses/zq0 (nqpar) [real,private/optional/default=0.0]