casa6
Public
Commits
| 6 6 | from math import sqrt, sin, cos, log |
| 7 7 | |
| 8 8 | from casatools import table, quanta, measures |
| 9 9 | from casatasks import casalog |
| 10 10 | |
| 11 11 | _tb = table() |
| 12 12 | _qa = quanta() |
| 13 13 | _me = measures() |
| 14 14 | |
| 15 15 | # debug |
| 16 - | debug = False |
| 16 + | _debug = False |
| 17 17 | |
| 18 18 | def gethorizonsephem(objectname, starttime, stoptime, incr, outtable, asis=False, rawdatafile=''): |
| 19 19 | """ |
| 20 20 | Main driver function for ephemeris data query from JPL-Horizons |
| 21 21 | |
| 22 22 | arguments: |
| 23 23 | objectname: ephemeris object name (case insensitive). Try to convert |
| 24 24 | a common name to ID if asis = False. |
| 25 25 | starttime: start time of the ephemeris data (expects YYYY/MM/DD/HH:MM |
| 26 26 | format |
| 85 85 | start_time = _qa.time(starttime, form='ymd') |
| 86 86 | stop_time = _qa.time(stoptime, form='ymd') |
| 87 87 | except ValueError as e: |
| 88 88 | casalog.post(e) |
| 89 89 | |
| 90 90 | try: |
| 91 91 | step_size = incr.replace(' ', '') |
| 92 92 | except ValueError as e: |
| 93 93 | casalog.post(e) |
| 94 94 | |
| 95 - | if debug: |
| 95 + | if _debug: |
| 96 96 | print("target=",target) |
| 97 97 | print(f"start_time={start_time}, stop_time={stop_time}, step_size={step_size}") |
| 98 98 | print(f"quantities={quantities}, ang_format={ang_format}") |
| 99 99 | ephemdata = queryhorizons(target, start_time, stop_time, step_size, quantities, ang_format, rawdatafile) |
| 100 100 | # return ephemdata |
| 101 101 | if ephemdata and 'result' in ephemdata: |
| 102 102 | casalog.post('converting ephemeris data to a CASA table') |
| 103 103 | tocasatb(ephemdata, outtable) |
| 104 104 | |
| 105 105 | def queryhorizons(target, starttime, stoptime, stepsize, quantities, ang_format, rawdatafile=''): |