{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "### Preparing a raw BOSZ Grid installation ###\n", "\n", "This is needed for making custom StarKit grids - not suitable for most users. \n", "\n", "You first need to navigate to the Phoenix folder that contains the grid ('bosz')" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "ExecuteTime": { "end_time": "2018-05-27T22:08:57.003460Z", "start_time": "2018-05-27T22:08:56.989575Z" }, "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/scigarfs/home/wkerzend/data/skgrid/bosz\n" ] } ], "source": [ "cd ~/data/skgrid/bosz/" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "ExecuteTime": { "end_time": "2018-05-27T22:09:13.686620Z", "start_time": "2018-05-27T22:08:57.659953Z" }, "collapsed": true }, "outputs": [], "source": [ "from starkit.gridkit.io.phoenix import PhoenixProcessGrid\n", "from starkit.gridkit.io.bosz.process import BOSZProcessGrid\n", "from starkit.gridkit.io.bosz.base import make_raw_index, make_grid_info, cache_bosz_grid\n", "from starkit.gridkit import load_grid\n", "import pandas as pd\n", "from astropy import units as u\n", "from astropy.io import fits\n", "import numpy as np\n", "import uuid" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "ExecuteTime": { "end_time": "2018-05-27T22:09:13.698972Z", "start_time": "2018-05-27T22:09:13.692794Z" }, "collapsed": true }, "outputs": [], "source": [ "#make_grid_info('bosz_grid_info.h5')" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "ExecuteTime": { "end_time": "2018-05-27T22:09:13.732207Z", "start_time": "2018-05-27T22:09:13.703529Z" }, "collapsed": true }, "outputs": [], "source": [ "#cache_bosz_grid()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "ExecuteTime": { "end_time": "2018-05-27T22:09:14.697612Z", "start_time": "2018-05-27T22:09:13.737263Z" }, "collapsed": true }, "outputs": [], "source": [ "meta = pd.read_hdf('bosz_grid_info.h5', 'meta')\n", "raw_index = pd.read_hdf('bosz_grid_info.h5', 'index')\n", "wavelength = pd.read_hdf('bosz_grid_info.h5', 'wavelength')[0].values * u.Unit(meta['wavelength_unit'])" ] }, { "cell_type": "code", "execution_count": 43, "metadata": { "ExecuteTime": { "end_time": "2018-05-27T22:20:46.934884Z", "start_time": "2018-05-27T22:20:46.911510Z" }, "collapsed": true }, "outputs": [], "source": [ "for col in raw_index.columns[:-1]:\n", " raw_index[col] = raw_index[col].astype(np.float64)\n", "raw_index.logg *= 0.1" ] }, { "cell_type": "code", "execution_count": 47, "metadata": { "ExecuteTime": { "end_time": "2018-05-27T22:21:08.437521Z", "start_time": "2018-05-27T22:21:08.417510Z" }, "collapsed": true }, "outputs": [], "source": [ "index_filter = (raw_index.teff.between(4000, 7000) &\n", " raw_index.logg.between(-1, 5) &\n", " raw_index.mh.between(-2.5, 0.5) &\n", " (raw_index.alpha == 0.0))\n", "\n", "new_index = raw_index.loc[index_filter]" ] }, { "cell_type": "code", "execution_count": 54, "metadata": { "ExecuteTime": { "end_time": "2018-05-27T22:24:03.079043Z", "start_time": "2018-05-27T22:24:03.031943Z" }, "collapsed": true }, "outputs": [], "source": [ "bgrid = BOSZProcessGrid(new_index, wavelength, meta, \n", " wavelength_start=2000*u.angstrom, \n", " wavelength_stop=25000*u.angstrom, R=20000.0)" ] }, { "cell_type": "code", "execution_count": 55, "metadata": { "ExecuteTime": { "end_time": "2018-05-27T22:32:23.589762Z", "start_time": "2018-05-27T22:24:05.015956Z" }, "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "100% (10530 of 10530) |###################| Elapsed Time: 0:07:47 Time: 0:07:47\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "done\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/lustre/home/wkerzend/miniconda3/envs/starkit/lib/python2.7/site-packages/pandas/core/generic.py:1299: PerformanceWarning: \n", "your performance may suffer as PyTables will pickle object types that it cannot\n", "map directly to c-types [inferred_type->mixed-integer,key->values] [items->None]\n", "\n", " return pytables.to_hdf(path_or_buf, key, self, **kwargs)\n" ] } ], "source": [ "bgrid.to_hdf('rcw86_fs1_bosz_grid.h5')" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "ExecuteTime": { "end_time": "2017-12-19T17:23:31.989687Z", "start_time": "2017-12-19T17:23:31.927872Z" }, "collapsed": false }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
teffloggmhalpha
0230000.50
123000.50.50
2230010.50
323001.50.50
4230020.50
523002.50.50
6230030.50
723003.50.50
8230040.50
923004.50.50
10230050.50
1123005.50.50
12230060.50
132400-0.50.50
14240000.50
1524000.50.50
16240010.50
1724001.50.50
18240020.50
1924002.50.50
20240030.50
2124003.50.50
22240040.50
2324004.50.50
24240050.50
2524005.50.50
26240060.50
272500-0.50.50
28250000.50
2925000.50.50
...............
47338140002-4-0.4
47339140002.5-4-0.4
47340140003-4-0.4
47341140003.5-4-0.4
47342140004-4-0.4
47343140004.5-4-0.4
47344140005-4-0.4
47345140005.5-4-0.4
47346140006-4-0.4
47347140006.5-4-0.4
47348145002-4-0.4
47349145002.5-4-0.4
47350145003-4-0.4
47351145003.5-4-0.4
47352145004-4-0.4
47353145004.5-4-0.4
47354145005-4-0.4
47355145005.5-4-0.4
47356145006-4-0.4
47357145006.5-4-0.4
47358150002-4-0.4
47359150002.5-4-0.4
47360150003-4-0.4
47361150003.5-4-0.4
47362150004-4-0.4
47363150004.5-4-0.4
47364150005-4-0.4
47365150005.5-4-0.4
47366150006-4-0.4
47367150006.5-4-0.4
\n", "

47368 rows × 4 columns

\n", "
" ], "text/plain": [ " teff logg mh alpha\n", "0 2300 0 0.5 0\n", "1 2300 0.5 0.5 0\n", "2 2300 1 0.5 0\n", "3 2300 1.5 0.5 0\n", "4 2300 2 0.5 0\n", "5 2300 2.5 0.5 0\n", "6 2300 3 0.5 0\n", "7 2300 3.5 0.5 0\n", "8 2300 4 0.5 0\n", "9 2300 4.5 0.5 0\n", "10 2300 5 0.5 0\n", "11 2300 5.5 0.5 0\n", "12 2300 6 0.5 0\n", "13 2400 -0.5 0.5 0\n", "14 2400 0 0.5 0\n", "15 2400 0.5 0.5 0\n", "16 2400 1 0.5 0\n", "17 2400 1.5 0.5 0\n", "18 2400 2 0.5 0\n", "19 2400 2.5 0.5 0\n", "20 2400 3 0.5 0\n", "21 2400 3.5 0.5 0\n", "22 2400 4 0.5 0\n", "23 2400 4.5 0.5 0\n", "24 2400 5 0.5 0\n", "25 2400 5.5 0.5 0\n", "26 2400 6 0.5 0\n", "27 2500 -0.5 0.5 0\n", "28 2500 0 0.5 0\n", "29 2500 0.5 0.5 0\n", "... ... ... ... ...\n", "47338 14000 2 -4 -0.4\n", "47339 14000 2.5 -4 -0.4\n", "47340 14000 3 -4 -0.4\n", "47341 14000 3.5 -4 -0.4\n", "47342 14000 4 -4 -0.4\n", "47343 14000 4.5 -4 -0.4\n", "47344 14000 5 -4 -0.4\n", "47345 14000 5.5 -4 -0.4\n", "47346 14000 6 -4 -0.4\n", "47347 14000 6.5 -4 -0.4\n", "47348 14500 2 -4 -0.4\n", "47349 14500 2.5 -4 -0.4\n", "47350 14500 3 -4 -0.4\n", "47351 14500 3.5 -4 -0.4\n", "47352 14500 4 -4 -0.4\n", "47353 14500 4.5 -4 -0.4\n", "47354 14500 5 -4 -0.4\n", "47355 14500 5.5 -4 -0.4\n", "47356 14500 6 -4 -0.4\n", "47357 14500 6.5 -4 -0.4\n", "47358 15000 2 -4 -0.4\n", "47359 15000 2.5 -4 -0.4\n", "47360 15000 3 -4 -0.4\n", "47361 15000 3.5 -4 -0.4\n", "47362 15000 4 -4 -0.4\n", "47363 15000 4.5 -4 -0.4\n", "47364 15000 5 -4 -0.4\n", "47365 15000 5.5 -4 -0.4\n", "47366 15000 6 -4 -0.4\n", "47367 15000 6.5 -4 -0.4\n", "\n", "[47368 rows x 4 columns]" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "raw_index.drop('filename', axis=1)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "ExecuteTime": { "end_time": "2017-12-19T16:32:50.976345Z", "start_time": "2017-12-19T16:32:50.971056Z" }, "collapsed": true }, "outputs": [], "source": [ "flux = fits.getdata('PHOENIX-ACES-AGSS-COND-2011/Z+0.5/lte02300-0.00+0.5.PHOENIX-ACES-AGSS-COND-2011-HiRes.fits')" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.11" }, "toc": { "nav_menu": {}, "number_sections": true, "sideBar": true, "skip_h1_title": false, "toc_cell": false, "toc_position": {}, "toc_section_display": "block", "toc_window_display": false } }, "nbformat": 4, "nbformat_minor": 1 }