Ejemplos de next_sunrise en Python

Ejemplo n.º 1

Archivo: lc2_019-ncp.py Proyecto: AdriaanRenting/lofar-obs-xml

#!/usr/bin/env python2

import momxml
import ephem
import sys

# 17:29 - 06:10

cal_duration_s   = 5*60
source_catalogue = momxml.SourceCatalogue()

mid_day    = momxml.ephem.Date(sys.argv[1])
start_date = momxml.ephem.Date(momxml.next_sunset(mid_day)  + 20*momxml.ephem.minute)
end_date   = momxml.ephem.Date(momxml.next_sunrise(start_date) - 20*momxml.ephem.minute)

total_duration_s = (end_date - start_date)*24*3600.0
target_duration_s = total_duration_s - 2*cal_duration_s - 2*61.0

target       = momxml.simbad('NCP')
target.name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]

pre_cal      = source_catalogue.cal_source(start_date, 'HBA')
post_cal     = source_catalogue.cal_source(start_date+(target_duration_s+2*cal_duration_s)*ephem.second, 'HBA')
cal_fields   = [momxml.simbad('3C61.1')]
ra_3c61_1    = cal_fields[0].ra_angle
dec_3c61_1   = cal_fields[0].dec_angle
ra_inc       = momxml.Angle(deg = 60.0)

antenna_set     = 'HBA_DUAL_INNER'
band            = 'HBA_LOW'
stations        = momxml.station_list('all', exclude = [])

Ejemplo n.º 2

Archivo: eor.py Proyecto: brentjens/lofar-obs-xml

import momxml
import ephem
import sys

# 17:29 - 06:10

cal_duration_s = 5 * 60
source_catalogue = momxml.SourceCatalogue()

mid_day = momxml.ephem.Date(sys.argv[1])
start_date = momxml.ephem.Date(
    momxml.next_sunset(mid_day) + 20 * momxml.ephem.minute)
end_date = momxml.ephem.Date(
    momxml.next_sunrise(mid_day) - 20 * momxml.ephem.minute)

total_duration_s = (end_date - start_date) * 24 * 3600.0
target_duration_s = total_duration_s - 2 * cal_duration_s - 2 * 61.0

target = momxml.simbad('NCP')
target.name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]

pre_cal = source_catalogue.cal_source(start_date, 'HBA')
post_cal = source_catalogue.cal_source(
    start_date + (target_duration_s + 2 * cal_duration_s) * ephem.second,
    'HBA')
cal_fields = [momxml.simbad('3C61.1')]
ra_3c61_1 = cal_fields[0].ra_angle
dec_3c61_1 = cal_fields[0].dec_angle
ra_inc = momxml.Angle(deg=60.0)

antenna_set = 'HBA_DUAL_INNER'

Ejemplo n.º 3

#!/usr/bin/env python2

import momxml
import ephem
import sys

# 17:29 - 06:10

cal_duration_s = 5 * 60
source_catalogue = momxml.SourceCatalogue()

mid_day = momxml.ephem.Date(sys.argv[1])
start_date = momxml.ephem.Date(
    momxml.next_sunset(mid_day) + 20 * momxml.ephem.minute)
end_date = momxml.ephem.Date(
    momxml.next_sunrise(start_date) - 20 * momxml.ephem.minute)

total_duration_s = (end_date - start_date) * 24 * 3600.0
target_duration_s = total_duration_s - 2 * cal_duration_s - 2 * 61.0

target = momxml.simbad('NCP')
target.name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]

pre_cal = source_catalogue.cal_source(start_date, 'HBA')
post_cal = source_catalogue.cal_source(
    start_date + (target_duration_s + 2 * cal_duration_s) * ephem.second,
    'HBA')
cal_fields = [momxml.simbad('3C61.1')]
ra_3c61_1 = cal_fields[0].ra_angle
dec_3c61_1 = cal_fields[0].dec_angle
ra_inc = momxml.Angle(deg=60.0)

Ejemplo n.º 4

Archivo: eor.py Proyecto: AdriaanRenting/lofar-obs-xml

import momxml
import ephem
import sys

# 17:29 - 06:10

cal_duration_s   = 5*60
source_catalogue = momxml.SourceCatalogue()

mid_day    = momxml.ephem.Date(sys.argv[1])
start_date = momxml.ephem.Date(momxml.next_sunset(mid_day)  + 20*momxml.ephem.minute)
end_date   = momxml.ephem.Date(momxml.next_sunrise(mid_day) - 20*momxml.ephem.minute)

total_duration_s = (end_date - start_date)*24*3600.0
target_duration_s = total_duration_s - 2*cal_duration_s - 2*61.0

target       = momxml.simbad('NCP')
target.name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]

pre_cal      = source_catalogue.cal_source(start_date, 'HBA')
post_cal     = source_catalogue.cal_source(start_date+(target_duration_s+2*cal_duration_s)*ephem.second, 'HBA')
cal_fields   = [momxml.simbad('3C61.1')]
ra_3c61_1    = cal_fields[0].ra_angle
dec_3c61_1   = cal_fields[0].dec_angle
ra_inc       = momxml.Angle(deg = 60.0)

antenna_set     = 'HBA_DUAL_INNER'
band            = 'HBA_LOW'
stations        = momxml.station_list('nl', exclude = [])
int_s           = 2.0
chan            = 64

Ejemplo n.º 5

Archivo: lt5_009-ncp.py Proyecto: SterVeen/lofar-obs-xml

#!/usr/bin/env python

import momxml
import ephem
import sys
from numpy import array

# 17:29 - 06:10

cal_duration_s   = 5*60
source_catalogue = momxml.SourceCatalogue()

mid_day    = momxml.ephem.Date(sys.argv[1])
start_date = momxml.ephem.Date(momxml.next_sunset(mid_day)  + 20*momxml.ephem.minute)
end_date   = momxml.ephem.Date(momxml.next_sunrise(start_date) - 20*momxml.ephem.minute)

total_duration_s = (end_date - start_date)*24*3600.0
target_duration_s = total_duration_s - 2*cal_duration_s - 2*61.0

target       = momxml.simbad('NCP')
target.name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]

pre_cal      = source_catalogue.cal_source(start_date, 'HBA')
post_cal     = source_catalogue.cal_source(start_date+(target_duration_s+2*cal_duration_s)*ephem.second, 'HBA')
cal_fields   = [momxml.simbad('3C61.1')]
ra_3c61_1    = cal_fields[0].ra_angle
dec_3c61_1   = cal_fields[0].dec_angle
ra_inc       = momxml.Angle(deg = 60.0)
ra_offset    = momxml.Angle(deg=30.0)

antenna_set     = 'HBA_DUAL_INNER'