calclos.f90 File Reference

Go to the source code of this file.

Modules

module  calclos_mod

Data Types

type  calclos_mod::param_type
 type to hold global information about the lines of sight More...
type  calclos_mod::pix_type
 line of sight pixel type More...
type  calclos_mod::work_particle_type
 local work particle type More...

Functions/Subroutines

subroutine calclos_mod::init_linesofsight (OmegaM, OmegaL, OmegaB, XH, a, h, npix, boxlen_kpch, origin, xyfile)
 subroutine that initializes the kernel table, physical parameters and reads the x-y coords file
subroutine calclos_mod::make_linetable ()
 Creates the linetable by doing a series of 1D line integrals through the kernel at different impact parameters.
real calclos_mod::interp_linetable (b)
 function that interpolates from the kernel line integral table.
subroutine calclos_mod::calc_linesofsight (pos, vel, mass, id, T, xHI, hsml)
 subroutine that calculates quantities along the lines of sight
subroutine calclos_mod::normalize_linesofsight ()
 normalizes the mass weighted quantities, converts rhob and rhoHI to critical density units and sets nHI and taureal
subroutine calclos_mod::convolve_linesofsight ()
 convolves the line of sight quantities with the velocity field and adds thermal broadening
subroutine calclos_mod::output_linesofsight (outfile)
 subroutine to output line of sight information ot file

Variables

real(8), parameter calclos_mod::zero = 0.0d0
 symbolic 0.0
real(8), parameter calclos_mod::one = 1.0d0
 symbolic 1.0
real(8), parameter calclos_mod::two = 2.0d0
 symbolic 2.0
real(8), parameter calclos_mod::GLen = 3.085678d21
 Gadget length unit [cm/h].
real(8), parameter calclos_mod::GMass = 1.989d43
 Gadget mass unit [g/h].
real(8), parameter calclos_mod::GVel = 1.0d5
 Gadget velocity unit [cm/s].
real(8), parameter calclos_mod::GTime = GLen / GVel
 Gadget time unit [s/h].
real(8), parameter calclos_mod::GRho = GMass / GLen**3
 Gadget density unit [cgs*h^2].
real(8), parameter calclos_mod::GEnrg = GMass * GVel**2
 Gadget energy unit [ergs/h].
real(8), parameter calclos_mod::Mpc2km = 3.08568025d19
 convert Mpc -> km
real(8), parameter calclos_mod::c = 2.99792458d10
 speed of light [cm/s]
real(8), parameter calclos_mod::protonmass = 1.67262158d-24
 protonmass [g]
real(8), parameter calclos_mod::lyac = 1.34d-12
 LyA cross-section [cm^2].
real(8), parameter calclos_mod::G_cgs = 6.67259d-8
 grav. constant [cgs]
real(8), parameter calclos_mod::pi = 3.1415926536d0
 pi
type(param_type) calclos_mod::params
 holds line of sight parameters
real(8) calclos_mod::z
 red shift
real(8) calclos_mod::OmegaFac
 sqrt(OmegaM * a^-3 + OmegaL)
real(8) calclos_mod::Mpch2kms
 hubble flow Mpc/h -> km/s
real(8) calclos_mod::kpch2kms
 hubble flow kpc/h -> km/s
real(8) calclos_mod::H0_cgs
 hubble parameter now in cgs [1/s]
real(8) calclos_mod::rhocrit_cgs
 critical density now [cgs]
real(8) calclos_mod::gascrit_cgs
 critical density * OmegaB [cgs]
real(8) calclos_mod::rhocrit_gad
 critical density now [gadget]
real(8) calclos_mod::gascrit_gad
 critical density * OmegaB [gadget]
real(8) calclos_mod::boxlen_kms
 box length [km/s]
real(8) calclos_mod::pixlen_kpch
 pixel length [kpc/h]
real(8) calclos_mod::pixlen_kms
 pixel length [km/s]
real(8) calclos_mod::pixlen_cm
 pixel length [cm]
type(pix_type), dimension(:,:),
allocatable 
calclos_mod::los
 lines of sight [npix,nlos]
real, dimension(:,:), allocatable calclos_mod::xycoords
 x-y coords [2,nlos]
real, dimension(:), allocatable calclos_mod::addbit
 tracks kernel quant. for pixels
real, dimension(:), allocatable calclos_mod::addtau
 tracks tau quant. for pixels
real, dimension(:), allocatable calclos_mod::zaxis_kms
 z values for pixels [km/s]
real, dimension(:), allocatable calclos_mod::zaxis_kpch
 z values for pixels [kpc/h]
integer, parameter calclos_mod::kbins = 1000
 kernel table bins
real(8), parameter calclos_mod::db = two / (kbins-1)
 kernel pixel spacing
real(8), dimension(kbins) calclos_mod::linetable
 line integral table
real(8), dimension(kbins) calclos_mod::btable
 impact parameters for table
real(8) calclos_mod::maxhsml
 maximum smoothing length
integer calclos_mod::ncells
 number of chaing mesh cells in 1D
real calclos_mod::cellsize
 linear extent of chaining mesh cell
integer, dimension(:,:),
allocatable 
calclos_mod::hoc
 head of cell array
integer, dimension(:), allocatable calclos_mod::ll
 linked list


Detailed Description

Definition in file calclos.f90.


Generated on Wed Mar 18 12:10:44 2009 for Fortran 90 SPH Analysis Tools by  doxygen 1.5.5