Gravity

Source terms for gravitational forces. More...

Modules
module	gravity_base_mod
	generic gravity terms module providing functionaly common to all gravity terms
module	gravity_binary_mod
	source terms module for gravitational acceleration due to two pointmasses
module	gravity_pointmass_mod
	source terms module for gravitational acceleration due to a point mass at the center of the coordinate system
module	gravity_sboxspectral_mod
	Poisson solver via spectral methods within the shearingsheet.
module	gravity_spectral_mod
	2D poisson solver using spectral methods for direct integration

Detailed Description

Source terms for gravitational forces.

• general parameters of gravity group as key-values

  INTEGER gtype of gravity source

  INTEGER output/accel (=1) output of acceleration

  INTEGER output/height (=1) output of disc height

• parameters of gravity_binary as key-values

  REAL mass1 = 1.0 of primary component

  REAL mass2 = 1.0 of secondary component

  REAL excentricity = 0.0

  REAL semimayoraxis = 1.0 major axis

  REAL softening1 = 0.0 parameter of primary component

  REAL softening2 = 0.0 parameter of secondary component

  REAL switchon1 = -1.0 switch on

  REAL omega_rot = 0.0 binary positions to a rotating reference frame

  INTEGER mesh_rot = 0 binary positions to a rotating mesh (0|1)

  with angular velocity omega

• parameters of gravity_pointmass as key-values

  INTEGER potential of the potential

  REAL mass = 1.0 of the point mass

  REAL x = 0.0 x-position of the point mass

  REAL y = 0.0 y-position of the point mass

  REAL z = 0.0 z-position of the point mass

  REAL softening = 0.0 (e.g. for planets inside the computational domain)

  REAL switchon = -1.0 switch on

  REAL acclimit = -1.0 accretion rate / mass [1/unit time] (if negative => disabled)

  INTEGER outbound = depends on mesh mass accretion by setting "outbound" to one of the four boundaries

• parameters of gravity_sboxspectral as key-values

  INTEGER output/phi (enable=1) output of gravitational potential

  INTEGER output/accel_x (enable=1) output of accel in x-direction

  INTEGER output/accel_y (enable=1) output of accel in y-direction

  INTEGER output/mass2D (enable=1) output of density that will be

  fourier transformed

  INTEGER output/Fmass2D (enable=1) output of density in spectral

  space

• parameters of gravity_spectral as key-values

  INTEGER green = 1 type of Green-function

  REAL sigma = 0.05 deviation

  INTEGER output/potential (=1) output of grav. potential