https://github.com/pseudotensor/grmonty
grmonty: relativistic Monte Carlo code (from Dolence et al. 2009)
https://github.com/pseudotensor/grmonty
Last synced: 8 months ago
JSON representation
grmonty: relativistic Monte Carlo code (from Dolence et al. 2009)
- Host: GitHub
- URL: https://github.com/pseudotensor/grmonty
- Owner: pseudotensor
- License: gpl-2.0
- Created: 2016-05-19T03:07:18.000Z (about 10 years ago)
- Default Branch: master
- Last Pushed: 2016-05-19T03:08:15.000Z (about 10 years ago)
- Last Synced: 2025-01-12T23:10:21.732Z (over 1 year ago)
- Language: C
- Size: 21.1 MB
- Stars: 3
- Watchers: 3
- Forks: 13
- Open Issues: 0
-
Metadata Files:
- Readme: README
- License: COPYING
Awesome Lists containing this project
README
GRMONTY is a relativistic Monte Carlo code.
GRMONTY is parallelized using openMP.
This version is configured to use input from harm2d, which is also available on this site.
quick start:
# unpack the tarball
tar -xzvf grmonty.tgz
# make (requires openmp enabled gcc)
make
# set number of threads (example is for csh and 8 threads)
setenv OMP_NUM_THREADS 8
# run the code on the supplied harm output file
# arguments are:
# -estimate of photon number (actual number is probabilistic due to scattering)
# -harm dump file for model
# -mass unit (few x 10^19 is appropriate for Sgr A*)
./grmonty 5000000 dump019 4.e19 >&! grmonty.err
# will output spectrum to grmonty.spec
# which should be identical to grmonty_spec_verify
Then use SM script plspec to plot up broad-band spectral energy distribution.
To calculate spectra from other sources, replace harm_model.c
with your own source model. Begin by modifying harm_model.c.
You must supply
init_model
make_super_photon
bias_func
get_fluid_params
report_spectrum
stop_criterion
record_criterion
gcon_func
gcov_func
get_connection
in the model file.