davidgardenier

from bokeh.io import save, export_svgs
from bokeh.plotting import figure, output_file, show

filename = "test.svg"

p = figure(plot_width=971, plot_height=600)
p.circle([1,2,3,4,5], [2,5,8,2,7], size=10)

# Set axis titles
p.xaxis.axis_label = "Important"

davidgardenier

import subprocess
import numpy as np

template = 'python ../code/mySpitzerTransit/4_finalPipelineMCMC.py '
template += '{0} {1} {2} {3} {4} {5} {6} > outfile_{7}.txt'

planets = ['Hatp32b',
           'XO1b',
           'Hatp1b',
           'Wasp17b',

davidgardenier

for f in `find ./ -name "*.zip"` ; do p=`pwd`; d=`dirname $f`; cd $d; b=`basename $f`; unzip -o $b; cd $p; done;
find . -name *.txt -exec rm -rf {} \;
for f in `find ./ -name "*.py"` ; do p=`pwd`; d=`dirname $f`; cd $d; b=`basename $f`; code2pdf -l -s a4 $b; cd $p; done
for f in */; do p=`pwd`; d=`dirname $f`; cd $d; pdftk "$f"*.pdf output "$p"/"${f::-1}".pdf; cd $p; done

davidgardenier

import sys
import time

# main
parms = sys.argv
for p in parms:
  time.sleep(5)

davidgardenier

import subprocess

template = 'python script.py {0} {1} {2} {3} > file_{0}_{1}_{2}_{3}.txt'

args = [[1, 2, 3, 4], [5, 6, 7, 8]]

# Run commands in parallel
processes = []

for arg in args:

davidgardenier

import subprocess

template = 'python script.py {} {} {} {}'

args = [[1, 2, 3, 4], [5, 6, 7, 8]]

# Run commands in parallel
processes = []

for arg in args:

davidgardenier

def load_T_sky():
    """
    Read the Haslam sky temperature map into a list from which temperatures can
    be retrieved. The temperature sky map is given in the weird units of
    HealPix, and despite looking up info on this coordinate system, I don't
    have the foggiest idea of how to transform these to galactic coordinates. I
    have therefore directly copied the following code from psrpoppy in the
    assumption Sam Bates managed to figure it out.
    Returns:
        t_sky_list (list): List of sky temperatures in HealPix? coordinates?

davidgardenier

import logging
import logging.config
import os
import sys

class Log(object):
    """Setting default logging styles. Only really helpful when invoking as

    >>> logger = Log().logger()

davidgardenier

def sky_frac(beam_size):
    """
    Calculate which fraction of the sky a beam covers

    Args:
        beam_size (float): In square degrees
    Returns:
        frac (float): Fraction of the sky that the beam covers
    """

davidgardenier

def test_sky_frac(self):

    frac = go.sky_frac(-90, 90, 0, 360)
    self.assertTrue(frac > 0.99999999999)

    frac = go.sky_frac(0, 90, 0, 360)
    self.assertTrue(0.49 < frac < 0.51)

    frac = go.sky_frac(0, 90, 0, 90)
    self.assertTrue(0.124 < frac < 0.126)