MOOG17 with scattering
Switch branches/tags
Nothing to show
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Permalink
Failed to load latest commit information.
.gitignore
Abfind.f
Abpop.f
Abunplot.f
AngWeight.f
Atmos.com
Barklem.dat
BarklemUV.dat
Batom.f
Begin.f
Binary.f
Binplot.f
Binplotprep.f
Blankstring.f
Blends.f
Bmolec.f
Boxit.f
Calmod.f
Cdcalc.f
Cdcalc_JS.f
Chabund.f
Cog.f
Cogplot.f
Cogsyn.f
Correl.f
Crosscorr.f
Curve.f
Dampdat.com
Damping.f
Defcolor.f
Discov.f
Doflux.f
Drawcurs.f
Dummy.com
Eqlib.f
Equivs.com
Estim.f
Ewfind.f
Ewweighted.f
Factor.com
Fakeline.f
Findtic.f
Finish.f
Fluxplot.f
Gammabark.f
Getasci.f
Getcount.f
Getnum.f
Getsyns.f
Gridplo.f
Gridsyn.f
Infile.f
Inlines.f
Inmodel.f
Invert.f
Jexpint.f
Kappa.com
LICENSE
Lineabund.f
Lineinfo.f
Linex.com
Linlimit.f
MOOGSILENT
Makefile.macdesk
Makefile.macdesksilent
Makefile.maclap
Makefile.maclapsilent
Makefile.rh
Makefile.rh64
Makefile.rh64silent
Makefile.rhsilent
Makefile.sobeck
Makeplot.f
Minimax.f
Mol.com
Molquery.f
Moog.f
Moogsilent.f
Multimod.com
Multistar.com
Mydriver.f
Nansi.f
Nearly.f
Number.f
Obshead.f
Obspars.com
Oneline.f
OpacHelium.f
OpacHydrogen.f
Opaccouls.f
Opacit.f
Opacmetals.f
Opacscat.f
Params.f
Partfn.f
Partnew.f
Plotit.f
Plotremember.f
Plotval.com
Pltabun.f
Pltcog.f
Pltflux.f
Pltspec.f
Pointcurs.f
Prinfo.f
Pstuff.com
Putasci.f
Quants.com
README.md
Readobs.f
Rinteg.f
Smooth.f
Source.com
Sourcefunc_scat_cont.f
Sourcefunc_scat_line.f
Specplot.f
Spectemp.f
Stats.f
Sunder.f
Synpop.f
Synspec.f
Synth.f
Tablepop.f
Taukap.f
Total.f
Trudamp.f
Ucalc.f
Vargauss.f
Vmacro.f
Voigt.f
Wavecalc.f
Weedout.f
Writenumber.f

README.md

moog17scat

MOOG 2017 version with scattering

IMPORTANT UPDATE: Sep 13, 2018.

A bug was introduced with the "trudamp" parameter in the 2014 refactoring. It was turned on by default, rather than kept off. This is now fixed to be off by default. If you would like to turn trudamp on, use "trudamp 1" in the parameters.

There is also a typo bug in Trudamp.f that has been fixed for Ca3933 (thanks to Rana Ezzeddine for identifying the bug).

The bug only affected these specific lines (see Damping.f and Trudamp.f):

  • Ca II (3933, 8498, 8542, 8662)
  • CH (3693)
  • Ca I (6717, 6318, 6343, 6361)
  • Ca I autoionization (6318, 6343, 6361)

IMPORTANT UPDATE: Aug 7, 2017.

I mistakenly introduced a bug causing problems with abfind. It is now fixed.

MOOG

MOOG is a radiative transfer code for stellar abundances written by Chris Sneden. http://www.as.utexas.edu/~chris/moog.html

The publicly available version of MOOG does not treat scattering properly for blue, metal-poor stars. This was fixed in 2011 by Jennifer Sobeck (Sobeck et al. 2011, AJ 141, 175), and her version was widely distributed. A proper treatment of scattering is important to accurately compute abundances of blue lines in metal-poor stars. However, MOOG has since experienced several upgrades, with a big update coming in 2014. Scattering was not propagated to this version.

Here, I have taken the publicly available version of MOOG (from Feb 2017) and patched in the code written by Jennifer Sobeck. It can be turned on by setting "scat 1" in the configuration (batch.par) files.

WARNING: Blends.f and Ewfind.f were not modified in the 2011 version I have. I may be mistaken, but it does not appear that they use Sobeck's scattering. I will eventually fix this, but this likely affected anything that used the blends driver for bluer lines (<4500A) in metal-poor stars. UPDATE: I have learned that some versions of the scattering code did have Blends and Ewfind fixed. If you would like to check your version, you can do so by looking to see if "cdcalc_JS" is in Blends.f or Ewfind.f. If so, then you are probably good!

Citation

If you use this code in a paper, please cite Sneden 1973 and Sobeck et al. 2011, AJ 141, 175. I would also appreciate if you include a link to this page in a footnote (www.github.com/alexji/moog17scat).

Installation Instructions

Instructions copied from http://www.as.utexas.edu/~chris/moog.html

  • Download/extract or clone the repository (green button in top right)
  • Edit the moogpath variable in Moog.f and Moogsilent.f to suit your download location
  • Compile MOOG with make -f Makefile.xxx where xxx is the setup that you want.
    • You will have to have a fortran compiler. g77 is the one assumed by MOOG by default.
    • For Mac OS X, and old version of g77 can be downloaded here: http://hpc.sourceforge.net/
    • However, gfortran is really the only fortran compiler supported on newer Mac laptops. So I recommend you use this instead.
    • I have edited Makefile.maclapsilent to work with gfortran. See that file if you want to see what compiler flags you should set to use gfortran.
  • Ensure this new MOOG or MOOGSILENT is your default in your path.
    • You can edit your .cshrc, .bashrc, and/or .bash_profile to add this new directory.
    • For example, my .bash_profile has export PATH="/Users/alexji/lib/moog17scat:$PATH"
    • Check this by opening a new terminal and typing which MOOGSILENT. It should point to the correct place.