def plot_airmass_moon():
# coordinates to center of OGLE field BLG502 (2017 objects)
ra = "18:00:00"
dec = "-30:00:00"
months = np.array([8, 9, 9])
days = np.array([15, 1, 15])
# outdir = '/Users/jlu/doc/proposals/keck/uc/18A/'
outdir = '/u/casey/scratch/code/JLU-python-code/jlu/papers/'
# Keck 2
skycalc.plot_airmass(ra,
dec,
2020,
months,
days,
'keck2',
outfile=outdir + 'microlens_airmass_keck2_20B.png',
date_idx=-1)
skycalc.plot_moon(ra,
dec,
2020,
np.array([8, 9]),
outfile=outdir + 'microlens_moon_20B.png')
# Keck 1
skycalc.plot_airmass(ra,
dec,
2020,
months,
days,
'keck1',
outfile=outdir + 'microlens_airmass_keck1_20B.png',
date_idx=-1)
return
def plot_airmass_moon():
# coordinates to center of OGLE field BLG502 (2017 objects)
ra = "18:00:00"
dec = "-30:00:00"
months = np.array([8, 8])
days = np.array([1, 30])
outdir = '/Users/fatima/Desktop/keck_prop_18B_mlens/'
# Keck 2
skycalc.plot_airmass(ra,
dec,
2018,
months,
days,
'keck2',
outfile=outdir + 'microlens_airmass_keck2.png',
date_idx=-1)
skycalc.plot_moon(ra,
dec,
2018,
months,
outfile=outdir + 'microlens_moon.png')
# Keck 1
#skycalc.plot_airmass(ra, dec, 2018, months, days, 'keck1', outfile=outdir + 'microlens_airmass_keck1.png', date_idx=-1)
return
def plot_airmass_moon():
# coordinates to center of OGLE field BLG502 (2017 objects)
ra = "18:00:00"
dec = "-30:00:00"
months = np.array([4, 5, 6, 7])
days = np.array([10, 10, 10, 30])
outdir = '/Users/jlu/doc/proposals/keck/uc/18A/'
# outdir = '/Users/fatima/Desktop/'
# Keck 2
skycalc.plot_airmass(ra,
dec,
2018,
months,
days,
'keck2',
outfile=outdir + 'microlens_airmass_keck2.png',
date_idx=-1)
skycalc.plot_moon(ra,
dec,
2018,
months,
outfile=outdir + 'microlens_moon.png')
# Keck 1
skycalc.plot_airmass(ra,
dec,
2018,
months,
days,
'keck1',
outfile=outdir + 'microlens_airmass_keck1.png',
date_idx=-1)
return
def gc_plot_airmass_moon():
# coordinates to GC
ra = "17:45:40.04"
dec = "-29:00:28.12"
months = np.array([4, 5, 6, 7])
days = np.array([10, 10, 10, 30])
outdir = '/Users/jlu/doc/proposals/ifa/16A/'
skycalc.plot_airmass(ra, dec, 2016, months, days, outfile=outdir + 'gc_airmass.png')
skycalc.plot_moon(ra, dec, 2016, months, outfile=outdir + 'gc_moon.png')
return
def gc_plot_airmass_moon():
# coordinates to GC
ra = "05:35:17.30"
dec = "-05:23:28.0"
months = np.array([10, 11, 12])
days = np.array([10, 10, 10])
outdir = '/Users/jlu/doc/proposals/ifa/16B/'
skycalc.plot_airmass(ra, dec, 2016, months, days, outfile=outdir + 'orion_airmass.png')
skycalc.plot_moon(ra, dec, 2016, months, outfile=outdir + 'orion_moon.png')
return
def gc_plot_airmass_moon():
# coordinates to GC
ra = "17:45:40.04"
dec = "-29:00:28.12"
months = np.array([4, 5, 6, 7])
days = np.array([10, 10, 10, 30])
outdir = '/Users/jlu/doc/proposals/keck/uc/17A/'
# Keck 2
skycalc.plot_airmass(ra, dec, 2017, months, days, 'keck2', outfile=outdir + 'gc_airmass_keck2.png')
skycalc.plot_moon(ra, dec, 2017, months, outfile=outdir + 'gc_moon.png')
# Keck 1
skycalc.plot_airmass(ra, dec, 2017, months, days, 'keck1', outfile=outdir + 'gc_airmass_keck1.png')
return
def plot_airmass_moon():
# coordinates to center of OGLE field BLG502 (2017 objects)
ra = "18:00:00"
dec = "-30:00:00"
months = np.array([8, 8])
days = np.array([1, 30])
outdir = '/Users/fatima/Desktop/keck_prop_18B_mlens/'
# Keck 2
skycalc.plot_airmass(ra, dec, 2018, months, days, 'keck2', outfile=outdir + 'microlens_airmass_keck2.png', date_idx=-1)
skycalc.plot_moon(ra, dec, 2018, months, outfile=outdir + 'microlens_moon.png')
# Keck 1
#skycalc.plot_airmass(ra, dec, 2018, months, days, 'keck1', outfile=outdir + 'microlens_airmass_keck1.png', date_idx=-1)
return
def gc_plot_airmass_moon():
# coordinates to GC
ra = "17:45:40.04"
dec = "-29:00:28.12"
months = np.array([4, 5, 6, 7])
days = np.array([10, 10, 10, 30])
outdir = '/Users/jlu/doc/proposals/ifa/16A/'
skycalc.plot_airmass(ra,
dec,
2016,
months,
days,
outfile=outdir + 'gc_airmass.png')
skycalc.plot_moon(ra, dec, 2016, months, outfile=outdir + 'gc_moon.png')
return
def gc_plot_airmass_moon():
# coordinates to GC
ra = "05:35:17.30"
dec = "-05:23:28.0"
months = np.array([10, 11, 12])
days = np.array([10, 10, 10])
outdir = '/Users/jlu/doc/proposals/ifa/16B/'
skycalc.plot_airmass(ra,
dec,
2016,
months,
days,
outfile=outdir + 'orion_airmass.png')
skycalc.plot_moon(ra, dec, 2016, months, outfile=outdir + 'orion_moon.png')
return
def plot_airmass_moon():
# coordinates to center of OGLE field BLG502 (2017 objects)
ra = "18:00:00"
dec = "-30:00:00"
months = np.array([4, 5, 6, 7])
days = np.array([10, 15, 15, 30])
# outdir = '/Users/jlu/doc/proposals/keck/uc/18A/'
# outdir = '/Users/fatima/Desktop/'
outdir = '/Users/casey/scratch/'
# Keck 2
skycalc.plot_airmass(ra, dec, 2019, months, days, 'keck2', outfile=outdir + 'microlens_airmass_keck2_19A.png', date_idx=-1)
skycalc.plot_moon(ra, dec, 2019, months, outfile=outdir + 'microlens_moon_19A.png')
# Keck 1
skycalc.plot_airmass(ra, dec, 2019, months, days, 'keck1', outfile=outdir + 'microlens_airmass_keck1_19A.png', date_idx=-1)
return