Exemplo n.º 1
0
def earth_rotation(epoch, jd):
    ERA = coords.earth_rotation_angle(jd, False)
    print 'ERA', ERA
    GAST = coords.greenwich_sidereal_time(jd, True) * pi / 12
    print '2000B GAST', GAST
    GMST = coords.greenwich_sidereal_time(jd, False) * pi / 12
    print 'GMST', GMST

    def cmpfunc(cstdout, ERA, GAST, GMST):
        from warnings import warn

        ls = [
            l.split(':')[-1].strip() for l in cstdout.split('\n')
            if l.strip() != ''
        ]
        dERA = ERA - float(ls[0])
        dGAST = GAST - float(ls[1])
        dGMST = GMST - float(ls[2])

        for s in ('dERA', 'dGAST', 'dGMST'):
            val = locals()[s]
            if abs(val) > 1e-8:
                warn('%s too large:%g !!!' % (s, val))

    return cmpfunc, ERA, GAST, GMST
Exemplo n.º 2
0
def greenwich_sidereal_time_in_radians():
    """Calculate using GMT (according to client's time settings)."""
    #Function astropysics.obstools.epoch_to_jd wants a decimal year as input
    #Function astropysics.obstools.calendar_to_jd can take a datetime object
    gmt_jd = obstools.calendar_to_jd(site_time_gmt_as_datetime())
    #Convert from hours to radians... 24hr = 2*pi
    return coords.greenwich_sidereal_time(gmt_jd) * pi / 12
Exemplo n.º 3
0
def greenwich_sidereal_time_in_radians():
    """Calculate using GMT (according to client's time settings)."""
    # Function astropysics.obstools.epoch_to_jd wants a decimal year as input
    # Function astropysics.obstools.calendar_to_jd can take a datetime object
    gmt_jd = obstools.calendar_to_jd(site_time_gmt_as_datetime())
    # Convert from hours to radians... 24hr = 2*pi
    return coords.greenwich_sidereal_time(gmt_jd) * pi / 12
Exemplo n.º 4
0
def earth_rotation(epoch,jd):
    ERA = coords.earth_rotation_angle(jd,False)
    print 'ERA',ERA
    GAST = coords.greenwich_sidereal_time(jd,True)*pi/12
    print '2000B GAST',GAST
    GMST = coords.greenwich_sidereal_time(jd,False)*pi/12
    print 'GMST', GMST
    
    def cmpfunc(cstdout,ERA,GAST,GMST):
        from warnings import warn
        
        ls = [l.split(':')[-1].strip() for l in cstdout.split('\n') if l.strip()!='']
        dERA = ERA-float(ls[0])
        dGAST = GAST-float(ls[1])
        dGMST = GMST-float(ls[2])
        
        for s in ('dERA','dGAST','dGMST'):
            val = locals()[s]
            if abs(val)>1e-8:
                warn('%s too large:%g !!!'%(s,val))
        
    return cmpfunc,ERA,GAST,GMST