Skip to content

Instantly share code, notes, and snippets.

from __future__ import (absolute_import, division, print_function,
unicode_literals)
from astroplan import (Observer, FixedTarget, observability_table,
AltitudeConstraint, AtNightConstraint)
from astropy.time import Time
import astropy.units as u
import numpy as np
from astropy.coordinates import SkyCoord
import matplotlib.pyplot as plt
@bmorris3
bmorris3 / spaceprobes.py
Created November 5, 2015 23:40
a .astronomy hack from sydney!
from astropy.utils.data import download_file
import matplotlib.pyplot as plt
import json
import numpy as np
import astropy.units as u
from astropy.time import Time
from astropy.coordinates import (SkyCoord, CartesianRepresentation,
SphericalRepresentation, AltAz, Angle,
EarthLocation, CylindricalRepresentation)
import ephem
@bmorris3
bmorris3 / creating.py
Last active December 17, 2015 20:37
Sample HDF5 experiments
import h5py
import numpy as np
# Hack for interactive use with iPython notebook:
try:
f.close()
except (NameError, ValueError):
pass
# Create archive
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@bmorris3
bmorris3 / XSEDE-starter.md
Last active February 22, 2016 17:53
Minimal XSEDE starter scripts

Getting Started with XSEDE/OSG

Logging in

Login to the XSEDE login node with:

ssh [email protected]

Login to the Open Science Grid:

gsissh submit-1.osg.xsede.org
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@bmorris3
bmorris3 / mpi_goose_chase.py
Last active April 13, 2016 18:48
How do you send numpy arrays with metadata between processes using mpi?
"""
Run with
mpiexec -n 11 mpi_goose_chase.py
"""
from __future__ import print_function
from mpi4py import MPI
import numpy as np
import time
setinst echelle
hedit *.fits dispaxis 1 add=yes verify=no show=yes update=yes
!ls *fits | grep -v 'bias' > inlist
!sed s/.fits/.c.fits/ inlist > outlist
cosmicrays (input="@inlist", output="@outlist", answer="yes", crmasks="", threshold=50.0, fluxratio=10.0, npasses=20, window="5", interactive=no, train=no,objects="", savefile="", plotfile="")
!ls bias*fits > biaslist
zerocombine (input="@biaslist", output="Zero.fits", combine="average", reject="avsigclip", ccdtype="", process=no, delete=no, clobber=no, scale="none", nlow=0, nhigh=1, nkeep=1, mclip=yes, lsigma=3.0, hsigma=3.0, rdnoise="7", gain="3.8", snoise="0.")
!cp outlist inlist
!sed s/.c.fits/.pc.fits/ inlist > outlist
ccdproc (images="@inlist", output="@outlist", ccdtype="", fixpix = yes, overscan=no, trim=yes, zerocor=yes, darkcor=no, flatcor=no, illumcor=no, fringecor=no, readcor=no, scancor=no, readaxis="line", fixfile = "database/badpix.txt", biassec="", trimsec="[200:1850,1:2048]", zero="Zero.fits", dark="", flat="", illum="", fringe="", minreplace