from Starfish.grid_tools import download_PHOENIX_models
except ImportError:
raise ImportError("You need to install Starfish.")
def _parser():
"""Take care of all the argparse stuff."""
parser = argparse.ArgumentParser(
description=
"Download Phoenix ACES test spectra. A single spectra each for an M0, M3, M6, and M9 M-dwarf."
)
parser.add_argument(
"-o",
"--outdir",
help="Output directory. Default is configured by the config.yaml.",
type=str,
default=join(config.pathdir, config.paths["phoenix_raw"]),
)
return parser.parse_args()
if __name__ == "__main__":
args = _parser()
outdir = args.outdir
models = np.array(list(product((2600, 2800, 3500, 3900), (4.5, ),
(0.0, ))))
download_PHOENIX_models(models, outdir)
sys.exit(0)
def AlphaPHOENIXModels():
params = [(6100, 4.5, 0.0, -0.2)]
outdir = os.path.join(test_base, "data", "phoenix")
download_PHOENIX_models(outdir, parameters=params)
yield outdir
Fitting for best-fit parameter values of Teff, logg, and Av.
Created on Thu Apr 14 12:02:45 2020
@author: Christian
"""
import numpy as np
from Starfish.grid_tools import download_PHOENIX_models
#%%
# Download Pheonix models with specified ranges
ranges = [[5700, 8600], [4.0, 6.0], [-0.5, 0.5]]
download_PHOENIX_models(path="PHOENIX", ranges=ranges)
#%%
# Setup a grid interface with the models
from Starfish.grid_tools import PHOENIXGridInterfaceNoAlpha
grid = PHOENIXGridInterfaceNoAlpha(path="PHOENIX")
# Setup an HDF5 interface in order to allow much quicker reading and writing
# than compared to loading FITS files over and over again.
from Starfish.grid_tools.instruments import SPEX
from Starfish.grid_tools import HDF5Creator
creator = HDF5Creator(grid,
"F_SPEX_grid.hdf5",
instrument=SPEX(),
wl_range=(0.9e4, np.inf),
def PHOENIXModels(grid_points):
ranges = [[6000, 6200], [4.0, 5.0], [-1.0, 0.0]]
outdir = os.path.join(test_base, "data", "phoenix")
download_PHOENIX_models(outdir, ranges=ranges)
yield outdir