def main(): site = sites.load(args.site) if args.date is None: date = datetime.date.today() else: date = args.date run(site, args.sources, args.start_utc, args.end_utc, date, show_extra_panels=args.extra_info)
def main(): site = sites.load(args.site) if (args.utc is not None): args.lst = site.utc_to_lst(args.utc, args.date) if (args.lst is not None) and (args.date is None): # A fixed LST is used, but no date is provided. # Fix date to today's date. args.date = datetime.date.today() modes.run(args.modename, site, args.lst, args.date, args.targets, \ args.testsources, args.calibrators, args)
def main(): site = sites.load(args.site) if args.lst is None: lst = site.lstnow() else: lst = args.lst alt_deg, az_deg = args.coords ra_deg, decl_deg = site.get_skyposn(alt_deg, az_deg, lst=lst) lst_hms = utils.deg_to_hmsstr(lst*15)[0] print "%s oriented at (alt=%g deg, az=%g deg) at LST=%s is pointed towards" \ "\n RA=%s\n Dec=%s" % (site.name, alt_deg, az_deg, lst_hms, \ utils.deg_to_hmsstr(ra_deg)[0], utils.deg_to_dmsstr(decl_deg)[0])
def main(): site = sites.load(args.site) if args.date is None: date = datetime.date.today() else: date = args.date print "Converting LST to UTC for %s on %s" % \ (site.name, date.strftime("%b %d, %Y")) for lst in args.lst: lst_hours = utils.parse_timestr(lst) lst_hms = utils.deg_to_hmsstr(lst_hours*15) utc_hours = site.lst_to_utc(lst_hours, date) utc_hms = utils.deg_to_hmsstr(utc_hours*15) print " LST: %s = UTC: %s" % (lst_hms[0], utc_hms[0])
def main(): site = sites.load(args.site) if args.date is None: date = datetime.date.today() else: date = args.date print "Converting UTC to LST for %s on %s" % \ (site.name, date.strftime("%b %d, %Y")) for utc in args.utc: utc_hours = utils.parse_timestr(utc) utc_hms = utils.deg_to_hmsstr(utc_hours*15) lst_hours = site.utc_to_lst(utc_hours, date) lst_hms = utils.deg_to_hmsstr(lst_hours*15) print " UTC: %s = LST: %s" % (utc_hms[0], lst_hms[0])
def get_riseset_info(site, src, utc=None, lst=None, date=None): """Convenience function to get rise/set info for a site/source. Inputs: site: The observing site to get info for. src: The source to generate info for. utc: The UTC to get info for. (Default: now) lst: The LST to get info for. (Default: now) date: The date to get info for. (Default: now) Outputs: riseset_info: A dictionary of rise/set info. """ if (utc is not None) and (lst is not None): raise ValueError("Only one of 'utc' and 'lst' can be provided.") if not isinstance(site, sites.base.BaseSite): site = sites.load(site) if not isinstance(src, sources.Source): src = sources.Source.from_string(src) if utc is not None: lst = site.utc_to_lst(utc, date) if (lst is not None) and (date is None): # A fixed LST is used, but no date is provided. # Fix date to today's date. date = datetime.date.today() if (lst is None) and (date is None): date = datetime.date.today() lst = site.lstnow() if (lst is None) and (date is not None): # Want current LST for a specified date. This doesn't make sense. raise ValueError("Incompatible inputs: 'lst' is None, but 'date' is not None!") rs_info = {'msg':"", 'circumpolar':False, 'neverrises':False} ra_deg, dec_deg = src.get_posn(lst, date) rs_info['ra_deg'] = ra_deg rs_info['dec_deg'] = dec_deg rs_info['ra_hmsstr'] = utils.deg_to_hmsstr(ra_deg, 2)[0] rs_info['dec_dmsstr'] = utils.deg_to_dmsstr(dec_deg, 2)[0] alt, az = src.get_altaz(site, lst, date) rs_info['alt_deg'] = alt[0] rs_info['az_deg'] = az[0] try: risetime, settime = src.get_rise_set_times(site, date) except errors.SourceIsCircumpolar: rs_info['msg'] = "Source is circumpolar." rs_info['circumpolar'] = True rs_info['is_visible'] = True rs_info['next_riseset'] = None except errors.SourceNeverRises: rs_info['msg'] = "Source never rises." rs_info['neverrises'] = True rs_info['is_visible'] = False rs_info['next_riseset'] = None except errors.MultipleRiseSets: rs_info['msg'] = "Multiple rise/set times?!" except: # Any other error raise else: if src.is_visible(site, lst, date): rs_info['is_visible'] = True rs_info['next_riseset'] = \ utils.deg_to_hmsstr(((settime-lst)%24)*15)[0] else: rs_info['is_visible'] = False rs_info['next_riseset'] = \ utils.deg_to_hmsstr(((risetime-lst)%24)*15)[0] rs_info['uptime'] = \ utils.deg_to_hmsstr(((settime-risetime)%24)*15)[0] rs_info['rise_lst'] = \ utils.deg_to_hmsstr((risetime%24)*15)[0] rs_info['rise_utc'] = \ utils.deg_to_hmsstr((site.lst_to_utc(risetime, \ date)%24)*15)[0] rs_info['set_lst'] = \ utils.deg_to_hmsstr((settime%24)*15)[0] rs_info['set_utc'] = \ utils.deg_to_hmsstr((site.lst_to_utc(settime, \ date)%24)*15)[0] return rs_info