GitHub - hummel/gadget

Branches Tags

Name		Name	Last commit message	Last commit date
Latest commit History 100 Commits
H2-cooling-ratios.dat		H2-cooling-ratios.dat
README		README
accel.c		accel.c
accrete.c		accrete.c
allocate.c		allocate.c
allvars.c		allvars.c
allvars.h		allvars.h
begrun.c		begrun.c
calc_photo.F		calc_photo.F
chebyshev.F		chebyshev.F
chemcool.c		chemcool.c
chemcool_consts.h		chemcool_consts.h
cheminmo.F		cheminmo.F
compute_gamma.F		compute_gamma.F
compute_gamma_ARS.c		compute_gamma_ARS.c
const_rates.F		const_rates.F
cool.h		cool.h
cool_func.F		cool_func.F
cool_util.F		cool_util.F
coolinmo.F		coolinmo.F
density.c		density.c
density_ARS.c		density_ARS.c
domain.c		domain.c
driftfac.c		driftfac.c
dvode.F		dvode.F
endrun.c		endrun.c
evolve_abundances.F		evolve_abundances.F
f2c.h		f2c.h
forcetree.c		forcetree.c
fs_data.h		fs_data.h
global.c		global.c
gravtree.c		gravtree.c
gravtree_forcetest.c		gravtree_forcetest.c
heat_ion_rates.c		heat_ion_rates.c
highn_gamma.h		highn_gamma.h
hydra.c		hydra.c
init.c		init.c
io.c		io.c
iohr.c		iohr.c
jac.F		jac.F
jxr.dat		jxr.dat
khcr.dat		khcr.dat
khxr.dat		khxr.dat
longrange.c		longrange.c
main.c		main.c
makefile		makefile
mol_data.h		mol_data.h
ngb.c		ngb.c
non_eq.h		non_eq.h
outputs.txt		outputs.txt
peano.c		peano.c
photoinit_ism.F		photoinit_ism.F
photoinit_lowZ.F		photoinit_lowZ.F
pm_nonperiodic.c		pm_nonperiodic.c
pm_periodic.c		pm_periodic.c
potential.c		potential.c
predict.c		predict.c
proto.h		proto.h
rate_eq.F		rate_eq.F
rate_eq_GMC.F		rate_eq_GMC.F
rate_eq_highn.F		rate_eq_highn.F
rate_eq_mol.F		rate_eq_mol.F
rate_eq_primordial.F		rate_eq_primordial.F
rate_eq_simple.F		rate_eq_simple.F
ray.c		ray.c
ray_ARS.c		ray_ARS.c
raytrace.c		raytrace.c
read_ic.c		read_ic.c
read_ichr.c		read_ichr.c
restart.c		restart.c
rsk_turbdriving_NGP.c		rsk_turbdriving_NGP.c
rsk_turbdriving_field.c		rsk_turbdriving_field.c
run.c		run.c
sd_metal.h		sd_metal.h
shield_data.h		shield_data.h
sink.c		sink.c
spline.F		spline.F
system.c		system.c
tags.h		tags.h
timestep.c		timestep.c
turbulence.h		turbulence.h
ucr.dat		ucr.dat
validate_rates.F		validate_rates.F

Repository files navigation

/* 
   The numbers in the "sinkmasses" file  are printed (e.g. line 112 in sink.c) 
   at regular time intervals and for every sink in the simulation:
*/
fprintf(sinkmasses,"%17.13g %8d %15.6g %15.6d %15.11g %15.11g %15.11g\n", 
	All.Time, P[i].ID, P[i].Mass, All.NumCurrentTiStep,
	P[i].Pos[0], P[i].Pos[1], P[i].Pos[2]);

/* 
   So that corresponds to the simulation time (in units of 1/(1+z)), the sink
   ID, sink Mass, simulation timestep, and the x, y, and z position of the 
   sink(s).
*/


/* 
   Numbers in "sinkdat" are printed whenever a particle has been accreted 
   (print command is in both sink.c and accrete.c): 
*/
fprintf(sinkdat,"%17.13g %8d %15.6g %15.6g %15.6g %10d %10d %15.6g\n", 
	All.Time, n_acc, r_sink_phys, 
	(SphP[k].Entropy/(SphP[k].Gamma-1.0))
	*pow(SphP[k].Density*a3inv,(SphP[k].Gamma-1.0)), 
	SphP[num_sink_part].Entropy, P[k].ID, 
	sink_ID, 
	SphP[num_sink_part].Pressure);

/*
  That corresponds to simulation time, number of particles the sink has accreted
  up to that point in that timestep, the physical accretion radius in kpc, the 
  internal energy of the accreted particle, the entropy of the accreting sink,
  ID of the accreting sink, and pressure of the accreting sink.
*/

/*
  Numbers in "sinktemp" are the same as in sinkmasses, but only for the main, 
  first-formed and most massive sink particle.  It's just one line that contains
  the most recent info on the main sink particle, as opposed to sinkmasses which
  contains records of the sink characteristics at regular time intervals and for
  every sink.
*/

/*
  Numbers in "sinkmergers" are printed only when two sinks have merged (print 
  command is in both sink.c and accrete.c):
*/

fprintf(sinkmergers,"%17.13g %8d %8d %15.6g %15.6g \n", 
	All.Time, P[k].ID, P[num_sink_part].ID, P[k].Mass, 
	P[num_sink_part].Mass);

/*
  This corresponds to simulation time, ID of the first of the two merged sinks,
  ID of the second of the two merged sinks, mass of the first sink, and mass of
  the second sink.  After the merger, the 'new' sink is given the ID of the 
  larger of the two original sinks.  In other words, sink mergers are tracked 
  such that it's the bigger sink eating the smaller sink.
*/

/*
  Numbers in "sinkangmom" are also printed every time a particle has been 
  accreted (print command is in both sink.c and accrete.c):
*/
fprintf(sinkangmom,"%17.13g %15.6g %15.6g %15.6g %15.6g %15.6g %15.6g %15.6g %15.6g %15.6g %10d\n", 
	All.Time, d_sink_phys, vx_phys, vy_phys, vz_phys, jx, jy, jz, 
	jtot, sink_Mass, sink_ID);
/*
  This is the simulation time,  physical distance in cm between the sink and the
  accreted particle, velocity in x, y, and z directions, specific angular 
  momentum of the accreted particle wrt the sink in x, y, and z directions, 
  total specific angular momentum of the accreted particle, mass of the 
  accreting sink, and ID of the accreting sink.
*/

/*
  As you can tell, a lot of the info printed may be unnecessary or repetitious, 
  so feel free to modify as you see fit.  Note that there are multiple places in
  sink.c and accrete.c where it will print to these files.  More specifically, 
  sinkmerger, sinkangmom, and sinkdat each have three spots in both accrete.c 
  and sink.c where there is a command to print to these files.
*/