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