case the output was not acceptable, we restarted the calculation using more iterations. In case of a run with unacceptable output, we selected a starting model that had a different log $g$ from the initial starting model and used it to restart the calculation. Then the previously described convergence test was performed and more restarts were done if it was necessary. If the output remained unacceptable, the effective temperature of the output model was changed by 10, 50, and 100K. If the output of any of these runs was acceptable, then it was used as an input model to calculate the atmosphere with the original effective temperature.

To test our model atmospheres, we replicated the Rosseland opacity calculations by @castelli01. For this we used Castelli’s scripts without any modification to calculate the ODF, and Rosseland opacities for the @grevesse01 abundances. These calculations concluded in perfect agreement within numerical precision. We then attempted to reproduce model atmospheres with the same parameters found on Castelli’s website with our scripts using the ODFs and Rosseland opacity files generated with @grevesse01 abundances. This test also concluded with near perfect agreement with only 0.1$-$0.2K maximum differences coming from the different version of compilers used.

MARCS 