disk_ICs.cpp File Reference

Definitions of initial conditions for hydrostatic disks. Note that the grid is mapped to 1D as i + (x_dim)*j + (x_dim*y_dim)*k. More...

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <cmath>
#include <string.h>
#include <time.h>
#include "global.h"
#include "grid3D.h"
#include "mpi_routines.h"
#include "io.h"
#include "error_handling.h"

Functions
Real	log_func (Real y)
Real	gz_halo_D3D (Real R, Real z, Real *hdp)
Real	gr_halo_D3D (Real R, Real z, Real *hdp)
Real	Sigma_disk_D3D (Real r, Real *hdp)
Real	gz_disk_D3D (Real R, Real z, Real *hdp)
Real	gr_disk_D3D (Real R, Real z, Real *hdp)
Real	phi_halo_D3D (Real R, Real z, Real *hdp)
Real	phi_disk_D3D (Real R, Real z, Real *hdp)
Real	phi_total_D3D (Real R, Real z, Real *hdp)
Real	phi_hot_halo_D3D (Real r, Real *hdp)
Real	z_hc_D3D (int k, Real dz, int nz, int ng)
Real	r_hc_D3D (int i, Real dr)
void	hydrostatic_ray_analytical_D3D (Real *rho, Real *r, Real *hdp, Real dr, int nr)
	Calculate the density at spherical radius r due to a hydrostatic halo. Uses an analytic expression normalized by the value of the potential at the cooling radius.
void	hydrostatic_column_isothermal_D3D (Real *rho, Real R, Real *hdp, Real dz, int nz, int ng)
	Calculate the 1D density distribution in a hydrostatic column, assuming an isothermal gas. Uses an iterative to scheme to determine the density at (R, z=0) relative to (R=0,z=0), then sets the densities according to an analytic expression.
void	hydrostatic_column_analytical_D3D (Real *rho, Real R, Real *hdp, Real dz, int nz, int ng)
	Calculate the 1D density distribution in a hydrostatic column. Uses an iterative to scheme to determine the density at (R, z=0) relative to (R=0,z=0), then sets the densities according to an analytic expression.
Real	determine_rho_eos_D3D (Real cs, Real Sigma_0, Real *hdp)
Real	halo_density_D3D (Real r, Real *r_halo, Real *rho_halo, Real dr, int nr)

Detailed Description

Definitions of initial conditions for hydrostatic disks. Note that the grid is mapped to 1D as i + (x_dim)*j + (x_dim*y_dim)*k.