Example #1
0
    iers = GLOBALutils.JPLiers(baryc_dir, mjd - 999.0, mjd + 999.0)
    obsradius, R0 = GLOBALutils.JPLR0(latitude, altitude)
    obpos = GLOBALutils.obspos(longitude, obsradius, R0)

    jplephem.set_ephemeris_dir(baryc_dir, ephemeris)
    jplephem.set_observer_coordinates(float(obpos[0]), float(obpos[1]),
                                      float(obpos[2]))

    res = jplephem.doppler_fraction(float(ra / 15.0), float(dec), long(mjd),
                                    float(mjd % 1), 1, 0.0)
    lbary_ltopo = 1.0 + res['frac'][0]
    bcvel_baryc = (lbary_ltopo - 1.0) * 2.99792458E5

    print "\t\tBarycentric velocity:", bcvel_baryc

    res = jplephem.pulse_delay(ra / 15.0, dec, int(mjd), mjd % 1, 1, 0.0)
    mbjd = mjd + res['delay'][0] / (3600.0 * 24.0)

    # Moon Phase Calculations
    gobs = ephem.Observer()
    gobs.name = 'VLT'
    gobs.lat = rad(latitude)  # lat/long in decimal degrees
    gobs.long = rad(longitude)
    gobs.date = h[0].header['DATE-OBS'].replace('T', ' ')
    mephem = ephem.Moon()
    mephem.compute(gobs)
    Mcoo = jplephem.object_track("Moon", int(mjd), float(mjd % 1), 1, 0.0)
    Mp = jplephem.barycentric_object_track("Moon", int(mjd), float(mjd % 1), 1,
                                           0.0)
    Sp = jplephem.barycentric_object_track("Sun", int(mjd), float(mjd % 1), 1,
                                           0.0)
Example #2
0
        if ra2 !=0 and dec2 != 0:
            RA = ra2
            DEC = dec2
        else:
            print '\t\tUsing the coordinates found in the image header.'

        # set info for compute the baricentric correction
        iers          = GLOBALutils.JPLiers( baryc_dir, scmjd-999.0, scmjd+999.0 )
        obsradius, R0 = GLOBALutils.JPLR0( latitude, altitude)
        obpos         = GLOBALutils.obspos( longitude, obsradius, R0 )
        jplephem.set_ephemeris_dir( baryc_dir , ephemeris )
        jplephem.set_observer_coordinates( obpos[0], obpos[1], obpos[2] )
        res         = jplephem.doppler_fraction(RA/15.0, DEC, int(scmjd), scmjd%1, 1, 0.0)
        lbary_ltopo = 1.0 + res['frac'][0]
        bcvel_baryc = ( lbary_ltopo - 1.0 ) * 2.99792458E5  #This in the barycentric velocity
        res         = jplephem.pulse_delay(RA/15.0, DEC, int(scmjd), scmjd%1, 1, 0.0)
        scmbjd      = scmjd + res['delay'][0] / (3600.0 * 24.0) #This is the modified barycentric julian day of the observation

        # set observatory info to retrive info about the moon
        gobs = ephem.Observer()
        gobs.name = 'VBT'
        gobs.lat  = rad(latitude)
        gobs.long = rad(longitude)
        #gobs.date = hd['UT-DATE'] + ' ' + hd['UT-TIME'].replace(':','_')
        gobs.date = hd['DATE-OBS'].replace('T',' ')

        mephem = ephem.Moon()
        mephem.compute(gobs)
        Mcoo = jplephem.object_track("Moon", int(scmjd), float(scmjd%1), 1, 0.0)
        Mp   = jplephem.barycentric_object_track("Moon", int(scmjd), float(scmjd%1), 1, 0.0)
        Sp   = jplephem.barycentric_object_track("Sun", int(scmjd), float(scmjd%1), 1, 0.0)