mirror
/
XEphem
mirror of https://github.com/XEphem/XEphem.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity
XEphem/libastro/constel.c

1802 lines
50 KiB
C
Raw Permalink Blame History

/* code to handle constellation boundaries and figures */
/*
Constellation boundaries:
Primary reference:
METHOD TO DETERMINE THE CONSTELLATION IN WHICH A POSITION IS LOCATED
Recently, Mr. Barry N. Rappaport of New Mexico State University
transcribed the constellation boundaries as fixed by the IAU into
machine-readable form. These have been transcribed by Dr. Nancy G.
Roman to make it possible to determine by computer the constellation
in which a position is located.
NSSDC catalog description:
6042 AN Catalog of Constellation Boundary Data (Delporte, E. 1930,
Cambridge Univ. Press)
Comment(s): includes constellation identification software
(ADC 1987; see Roman, N.G. 1987, Publ. Astron. Soc. Pacific
99, 695); 23 description, 118 software, 358 data records.
3 files: 23x80, 118x80, 358x29
Further adapted for xephem by:
Craig Counterman: conversion from original F77 to C
Elwood Downey: incorporation into xephem
Ernie Wright: additional speed and time improvments
Constellation figures:
Chris Marriott: original figures
Elwood Downey: incorporated into xephem
*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include <math.h>
#include "astro.h"
/*
======================================================================
Ernie Wright 2 Mar 94
Find the constellation for a given position.
First C version by Craig Counterman and Elwood Downey. Based on a
FORTRAN program by Nancy G. Roman (Roman, N.G. 1987, Publ. Astron.
Soc. Pacific 99, 695). IAU constellation boundaries transcribed into
machine-readable form by Barry N. Rappaport, New Mexico State Univ.
======================================================================
*/
#define NBOUNDS 357
/* constellation ids */
#define And 0
#define Ant 1
#define Aps 2
#define Aql 3
#define Aqr 4
#define Ara 5
#define Ari 6
#define Aur 7
#define Boo 8
#define CMa 9
#define CMi 10
#define CVn 11
#define Cae 12
#define Cam 13
#define Cap 14
#define Car 15
#define Cas 16
#define Cen 17
#define Cep 18
#define Cet 19
#define Cha 20
#define Cir 21
#define Cnc 22
#define Col 23
#define Com 24
#define CrA 25
#define CrB 26
#define Crt 27
#define Cru 28
#define Crv 29
#define Cyg 30
#define Del 31
#define Dor 32
#define Dra 33
#define Equ 34
#define Eri 35
#define For 36
#define Gem 37
#define Gru 38
#define Her 39
#define Hor 40
#define Hya 41
#define Hyi 42
#define Ind 43
#define LMi 44
#define Lac 45
#define Leo 46
#define Lep 47
#define Lib 48
#define Lup 49
#define Lyn 50
#define Lyr 51
#define Men 52
#define Mic 53
#define Mon 54
#define Mus 55
#define Nor 56
#define Oct 57
#define Oph 58
#define Ori 59
#define Pav 60
#define Peg 61
#define Per 62
#define Phe 63
#define Pic 64
#define PsA 65
#define Psc 66
#define Pup 67
#define Pyx 68
#define Ret 69
#define Scl 70
#define Sco 71
#define Sct 72
#define Se1 73
#define Sex 74
#define Sge 75
#define Sgr 76
#define Tau 77
#define Tel 78
#define TrA 79
#define Tri 80
#define Tuc 81
#define UMa 82
#define UMi 83
#define Vel 84
#define Vir 85
#define Vol 86
#define Vul 87
#define Se2 88
static char *cns_namemap[ NCNS ] = {
/* 0 */ "And: Andromeda",
/* 1 */ "Ant: Antlia",
/* 2 */ "Aps: Apus",
/* 3 */ "Aql: Aquila",
/* 4 */ "Aqr: Aquarius",
/* 5 */ "Ara: Ara",
/* 6 */ "Ari: Aries",
/* 7 */ "Aur: Auriga",
/* 8 */ "Boo: Bootes",
/* 9 */ "CMa: Canis Major",
/* 10 */ "CMi: Canis Minor",
/* 11 */ "CVn: Canes Venatici",
/* 12 */ "Cae: Caelum",
/* 13 */ "Cam: Camelopardalis",
/* 14 */ "Cap: Capricornus",
/* 15 */ "Car: Carina",
/* 16 */ "Cas: Cassiopeia",
/* 17 */ "Cen: Centaurus",
/* 18 */ "Cep: Cepheus",
/* 19 */ "Cet: Cetus",
/* 20 */ "Cha: Chamaeleon",
/* 21 */ "Cir: Circinus",
/* 22 */ "Cnc: Cancer",
/* 23 */ "Col: Columba",
/* 24 */ "Com: Coma Berenices",
/* 25 */ "CrA: Corona Australis",
/* 26 */ "CrB: Corona Borealis",
/* 27 */ "Crt: Crater",
/* 28 */ "Cru: Crux",
/* 29 */ "Crv: Corvus",
/* 30 */ "Cyg: Cygnus",
/* 31 */ "Del: Delphinus",
/* 32 */ "Dor: Dorado",
/* 33 */ "Dra: Draco",
/* 34 */ "Equ: Equuleus",
/* 35 */ "Eri: Eridanus",
/* 36 */ "For: Fornax",
/* 37 */ "Gem: Gemini",
/* 38 */ "Gru: Grus",
/* 39 */ "Her: Hercules",
/* 40 */ "Hor: Horologium",
/* 41 */ "Hya: Hydra",
/* 42 */ "Hyi: Hydrus",
/* 43 */ "Ind: Indus",
/* 44 */ "LMi: Leo Minor",
/* 45 */ "Lac: Lacerta",
/* 46 */ "Leo: Leo",
/* 47 */ "Lep: Lepus",
/* 48 */ "Lib: Libra",
/* 49 */ "Lup: Lupus",
/* 50 */ "Lyn: Lynx",
/* 51 */ "Lyr: Lyra",
/* 52 */ "Men: Mensa",
/* 53 */ "Mic: Microscopium",
/* 54 */ "Mon: Monoceros",
/* 55 */ "Mus: Musca",
/* 56 */ "Nor: Norma",
/* 57 */ "Oct: Octans",
/* 58 */ "Oph: Ophiuchus",
/* 59 */ "Ori: Orion",
/* 60 */ "Pav: Pavo",
/* 61 */ "Peg: Pegasus",
/* 62 */ "Per: Perseus",
/* 63 */ "Phe: Phoenix",
/* 64 */ "Pic: Pictor",
/* 65 */ "PsA: Piscis Austrinus",
/* 66 */ "Psc: Pisces",
/* 67 */ "Pup: Puppis",
/* 68 */ "Pyx: Pyxis",
/* 69 */ "Ret: Reticulum",
/* 70 */ "Scl: Sculptor",
/* 71 */ "Sco: Scorpius",
/* 72 */ "Sct: Scutum",
/* 73 */ "Se1: Serpens Caput",
/* 74 */ "Sex: Sextans",
/* 75 */ "Sge: Sagitta",
/* 76 */ "Sgr: Sagittarius",
/* 77 */ "Tau: Taurus",
/* 78 */ "Tel: Telescopium",
/* 79 */ "TrA: Triangulum Australe",
/* 80 */ "Tri: Triangulum",
/* 81 */ "Tuc: Tucana",
/* 82 */ "UMa: Ursa Major",
/* 83 */ "UMi: Ursa Minor",
/* 84 */ "Vel: Vela",
/* 85 */ "Vir: Virgo",
/* 86 */ "Vol: Volans",
/* 87 */ "Vul: Vulpecula",
/* 88 */ "Se2: Serpens Cauda",
};
static struct {
unsigned short lower_ra; /* hours * 1800 */
unsigned short upper_ra; /* hours * 1800 */
short lower_dec; /* degrees * 60 */
short index;
} cbound[ NBOUNDS ] = {
{ 0, 43200, 5280, UMi },
{ 14400, 26100, 5190, UMi },
{ 37800, 41400, 5170, UMi },
{ 32400, 37800, 5160, UMi },
{ 0, 14400, 5100, Cep },
{ 16500, 19200, 4920, Cam },
{ 0, 9000, 4800, Cep },
{ 19200, 26100, 4800, Cam },
{ 31500, 32400, 4800, UMi },
{ 36300, 37800, 4800, Dra },
{ 0, 6315, 4620, Cep },
{ 20700, 24450, 4620, Cam },
{ 29760, 31500, 4500, UMi },
{ 36300, 37200, 4500, Cep },
{ 14340, 16500, 4410, Cam },
{ 16500, 20400, 4410, Dra },
{ 23400, 29760, 4200, UMi },
{ 5580, 6150, 4080, Cas },
{ 36750, 37200, 4020, Dra },
{ 20400, 21600, 3990, Dra },
{ 0, 600, 3960, Cep },
{ 25200, 28200, 3960, UMi },
{ 42450, 43200, 3960, Cep },
{ 21600, 24300, 3840, Dra },
{ 24300, 25950, 3780, Dra },
{ 41700, 42450, 3780, Cep },
{ 10980, 12600, 3720, Cam },
{ 36000, 36750, 3690, Dra },
{ 36966, 37080, 3655, Cep },
{ 12600, 14340, 3600, Cam },
{ 14340, 15150, 3600, UMa },
{ 35580, 36000, 3570, Dra },
{ 36000, 36966, 3570, Cep },
{ 41160, 41700, 3545, Cep },
{ 0, 4380, 3510, Cas },
{ 34950, 35580, 3480, Dra },
{ 3060, 3435, 3450, Cas },
{ 4380, 5580, 3420, Cas },
{ 5580, 5700, 3420, Cam },
{ 40170, 41160, 3375, Cep },
{ 9000, 10980, 3360, Cam },
{ 25260, 25950, 3330, UMa },
{ 25950, 34950, 3330, Dra },
{ 5700, 6000, 3300, Cam },
{ 39840, 40170, 3300, Cep },
{ 37080, 39540, 3290, Cep },
{ 0, 3060, 3240, Cas },
{ 10980, 11700, 3240, Lyn },
{ 21750, 24300, 3180, UMa },
{ 27450, 28350, 3180, Dra },
{ 39540, 39840, 3165, Cep },
{ 6000, 9000, 3150, Cam },
{ 41160, 42000, 3150, Cas },
{ 28350, 30600, 3090, Dra },
{ 3675, 4530, 3030, Per },
{ 30600, 32820, 3030, Dra },
{ 0, 2460, 3000, Cas },
{ 2460, 3000, 3000, Per },
{ 11700, 12240, 3000, Lyn },
{ 42000, 43200, 3000, Cas },
{ 24300, 25260, 2910, UMa },
{ 0, 2010, 2880, Cas },
{ 42450, 43200, 2880, Cas },
{ 32715, 32820, 2850, Her },
{ 32820, 34350, 2850, Dra },
{ 34350, 34500, 2850, Cyg },
{ 3000, 3675, 2820, Per },
{ 15150, 16500, 2820, UMa },
{ 300, 1560, 2760, Cas },
{ 21600, 21750, 2700, UMa },
{ 12240, 13260, 2670, Lyn },
{ 39435, 39540, 2640, Cyg },
{ 39375, 39435, 2625, Cyg },
{ 34500, 34920, 2610, Cyg },
{ 16500, 18300, 2520, UMa },
{ 18300, 19410, 2400, UMa },
{ 27780, 28350, 2400, Boo },
{ 28350, 29400, 2400, Her },
{ 16650, 17250, 2385, Lyn },
{ 0, 4530, 2205, And },
{ 4530, 4620, 2205, Per },
{ 34845, 34920, 2190, Lyr },
{ 8100, 8445, 2160, Per },
{ 39120, 39375, 2160, Cyg },
{ 39375, 39600, 2160, Lac },
{ 11760, 13260, 2130, Aur },
{ 13260, 13950, 2130, Lyn },
{ 0, 3600, 2100, And },
{ 39600, 41070, 2100, Lac },
{ 41070, 41160, 2070, Lac },
{ 41160, 42300, 2070, And },
{ 4620, 4890, 2040, Per },
{ 19410, 19800, 2040, UMa },
{ 21600, 22200, 2040, CVn },
{ 13950, 16650, 2010, Lyn },
{ 16650, 17790, 2010, LMi },
{ 1290, 2535, 1980, And },
{ 27330, 27780, 1980, Boo },
{ 42300, 42750, 1925, And },
{ 22200, 23850, 1920, CVn },
{ 42750, 43200, 1880, And },
{ 25125, 25260, 1845, CVn },
{ 4350, 4890, 1840, Tri },
{ 4890, 8100, 1840, Per },
{ 8100, 8550, 1800, Aur },
{ 32715, 34845, 1800, Lyr },
{ 19800, 21600, 1740, UMa },
{ 35400, 37650, 1740, Cyg },
{ 8550, 10590, 1710, Aur },
{ 17790, 18900, 1710, LMi },
{ 23850, 25125, 1710, CVn },
{ 0, 120, 1680, And },
{ 2535, 3000, 1680, Tri },
{ 10590, 11760, 1680, Aur },
{ 14190, 14400, 1680, Gem },
{ 37650, 39120, 1680, Cyg },
{ 34665, 35400, 1650, Cyg },
{ 3450, 4350, 1635, Tri },
{ 29100, 29400, 1620, CrB },
{ 27150, 27330, 1560, Boo },
{ 27330, 29100, 1560, CrB },
{ 33060, 33960, 1560, Lyr },
{ 19350, 19800, 1530, LMi },
{ 33960, 34665, 1530, Lyr },
{ 3000, 3450, 1500, Tri },
{ 1290, 1530, 1425, Psc },
{ 18900, 19350, 1410, LMi },
{ 38250, 38550, 1410, Vul },
{ 10260, 10590, 1370, Tau },
{ 120, 255, 1320, And },
{ 28650, 28860, 1320, Se1 },
{ 10590, 11190, 1290, Gem },
{ 35700, 36450, 1275, Vul },
{ 33960, 34650, 1265, Vul },
{ 255, 1530, 1260, And },
{ 36450, 37020, 1230, Vul },
{ 14055, 14190, 1200, Gem },
{ 37020, 38250, 1170, Vul },
{ 34650, 35700, 1150, Vul },
{ 5910, 6060, 1140, Ari },
{ 33960, 34200, 1110, Sge },
{ 10260, 10380, 1080, Ori },
{ 11190, 11355, 1050, Gem },
{ 34200, 35700, 970, Sge },
{ 8940, 9600, 960, Tau },
{ 28650, 28950, 960, Her },
{ 35700, 36450, 945, Sge },
{ 8310, 8940, 930, Tau },
{ 9600, 10080, 930, Tau },
{ 23100, 24300, 900, Com },
{ 31050, 32850, 860, Her },
{ 21360, 23100, 840, Com },
{ 13500, 14055, 810, Gem },
{ 30150, 31050, 770, Her },
{ 0, 255, 750, Peg },
{ 10080, 10380, 750, Tau },
{ 12600, 13500, 750, Gem },
{ 38010, 38400, 750, Peg },
{ 11355, 12480, 720, Gem },
{ 32850, 33960, 720, Her },
{ 37575, 37890, 710, Del },
{ 37890, 38010, 710, Peg },
{ 20730, 21360, 660, Leo },
{ 11235, 11355, 600, Ori },
{ 12480, 12600, 600, Gem },
{ 14055, 14265, 600, Cnc },
{ 42900, 43200, 600, Peg },
{ 3000, 5910, 595, Ari },
{ 36255, 36540, 510, Del },
{ 24300, 27150, 480, Boo },
{ 40950, 42900, 450, Peg },
{ 14265, 16650, 420, Cnc },
{ 16650, 19350, 420, Leo },
{ 32850, 33592, 375, Oph },
{ 33592, 33960, 375, Aql },
{ 37500, 37575, 360, Del },
{ 12600, 12630, 330, CMi },
{ 32850, 33165, 270, Se2 },
{ 28950, 30150, 240, Her },
{ 32850, 33165, 180, Oph },
{ 38640, 39000, 165, Peg },
{ 0, 3600, 120, Psc },
{ 33450, 33960, 120, Se2 },
{ 36540, 37500, 120, Del },
{ 37500, 38400, 120, Equ },
{ 38400, 38640, 120, Peg },
{ 39600, 40950, 120, Peg },
{ 39000, 39600, 105, Peg },
{ 12630, 12960, 90, CMi },
{ 6450, 8310, 0, Tau },
{ 8310, 8400, 0, Ori },
{ 12960, 14550, 0, CMi },
{ 26400, 27150, 0, Vir },
{ 32100, 32850, 0, Oph },
{ 4770, 5910, -105, Cet },
{ 5910, 6450, -105, Tau },
{ 27150, 29280, -195, Se1 },
{ 8400, 9150, -240, Ori },
{ 10500, 11235, -240, Ori },
{ 32100, 32340, -240, Se2 },
{ 32850, 33450, -240, Se2 },
{ 33450, 33960, -240, Aql },
{ 40950, 42900, -240, Psc },
{ 19350, 20730, -360, Leo },
{ 20730, 21300, -360, Vir },
{ 0, 600, -420, Psc },
{ 42900, 43200, -420, Psc },
{ 25650, 26400, -480, Vir },
{ 28650, 29280, -480, Oph },
{ 36000, 36960, -540, Aql },
{ 38400, 39360, -540, Aqr },
{ 30900, 32340, -600, Oph },
{ 10500, 14550, -660, Mon },
{ 8850, 9150, -660, Eri },
{ 9150, 10500, -660, Ori },
{ 14550, 15060, -660, Hya },
{ 17250, 19350, -660, Sex },
{ 21300, 23100, -660, Vir },
{ 31650, 31800, -700, Oph },
{ 33960, 36000, -722, Aql },
{ 8700, 8850, -870, Eri },
{ 36960, 38400, -900, Aqr },
{ 30900, 32850, -960, Se2 },
{ 32850, 33960, -960, Sct },
{ 15060, 15450, -1020, Hya },
{ 29280, 29475, -1095, Oph },
{ 15450, 16350, -1140, Hya },
{ 19350, 19500, -1140, Crt },
{ 29280, 29475, -1155, Sco },
{ 28200, 28650, -1200, Lib },
{ 22650, 23100, -1320, Crv },
{ 23100, 25650, -1320, Vir },
{ 16350, 17550, -1440, Hya },
{ 3000, 4770, -1463, Cet },
{ 4770, 6750, -1463, Eri },
{ 19500, 21300, -1470, Crt },
{ 21300, 22650, -1470, Crv },
{ 25650, 26850, -1470, Lib },
{ 29280, 30150, -1475, Oph },
{ 0, 3000, -1530, Cet },
{ 38400, 39360, -1530, Cap },
{ 39360, 42900, -1530, Aqr },
{ 42900, 43200, -1530, Cet },
{ 17550, 18450, -1590, Hya },
{ 8460, 8700, -1635, Eri },
{ 8700, 11010, -1635, Lep },
{ 36000, 38400, -1680, Cap },
{ 18450, 19050, -1750, Hya },
{ 22650, 26850, -1770, Hya },
{ 26850, 28200, -1770, Lib },
{ 28200, 28800, -1770, Sco },
{ 8250, 8460, -1800, Eri },
{ 30150, 31680, -1800, Oph },
{ 31680, 32100, -1800, Sgr },
{ 19050, 19500, -1870, Hya },
{ 11010, 13260, -1980, CMa },
{ 22050, 22650, -1980, Hya },
{ 19500, 22050, -2100, Hya },
{ 6300, 6750, -2160, For },
{ 15060, 16860, -2205, Pyx },
{ 7680, 8250, -2220, Eri },
{ 32100, 34500, -2220, Sgr },
{ 38400, 41400, -2220, PsA },
{ 41400, 42000, -2220, Scl },
{ 5400, 6300, -2375, For },
{ 16860, 19800, -2385, Ant },
{ 0, 3000, -2400, Scl },
{ 3000, 5400, -2400, For },
{ 6960, 7680, -2400, Eri },
{ 42000, 43200, -2400, Scl },
{ 25500, 26850, -2520, Cen },
{ 28200, 28800, -2520, Lup },
{ 28800, 29557, -2520, Sco },
{ 8700, 9000, -2580, Cae },
{ 9000, 11850, -2580, Col },
{ 14400, 15060, -2580, Pup },
{ 6150, 6960, -2640, Eri },
{ 29557, 32100, -2730, Sco },
{ 32100, 34500, -2730, CrA },
{ 34500, 36600, -2730, Sgr },
{ 36600, 38400, -2730, Mic },
{ 5400, 6150, -2760, Eri },
{ 8100, 8700, -2790, Cae },
{ 27600, 28200, -2880, Lup },
{ 0, 4200, -2890, Phe },
{ 4800, 5400, -2940, Eri },
{ 7350, 7680, -2940, Hor },
{ 7680, 8100, -2940, Cae },
{ 38400, 39600, -3000, Gru },
{ 10800, 14400, -3045, Pup },
{ 14400, 14700, -3045, Vel },
{ 4350, 4800, -3060, Eri },
{ 6900, 7350, -3060, Hor },
{ 0, 3300, -3090, Phe },
{ 10800, 11100, -3150, Car },
{ 14700, 15210, -3180, Vel },
{ 6300, 6900, -3190, Hor },
{ 6900, 7200, -3190, Dor },
{ 0, 2850, -3210, Phe },
{ 3900, 4350, -3240, Eri },
{ 8100, 9000, -3240, Pic },
{ 27090, 27600, -3240, Lup },
{ 15210, 15900, -3270, Vel },
{ 11100, 11700, -3300, Car },
{ 21300, 23100, -3300, Cen },
{ 25500, 27090, -3300, Lup },
{ 27090, 27600, -3300, Nor },
{ 7200, 7800, -3390, Dor },
{ 15900, 19800, -3390, Vel },
{ 19800, 20250, -3390, Cen },
{ 31500, 32400, -3420, Ara },
{ 32400, 36600, -3420, Tel },
{ 39600, 42000, -3420, Gru },
{ 5760, 6300, -3450, Hor },
{ 9000, 9900, -3450, Pic },
{ 11700, 12300, -3480, Car },
{ 0, 2400, -3510, Phe },
{ 2400, 3900, -3510, Eri },
{ 42000, 43200, -3510, Phe },
{ 7800, 8250, -3540, Dor },
{ 27600, 29557, -3600, Nor },
{ 36600, 38400, -3600, Ind },
{ 9900, 10800, -3660, Pic },
{ 27300, 27600, -3660, Cir },
{ 29557, 29850, -3660, Ara },
{ 26850, 27300, -3815, Cir },
{ 29850, 30150, -3815, Ara },
{ 10800, 12300, -3840, Pic },
{ 12300, 16260, -3840, Car },
{ 20250, 21300, -3840, Cen },
{ 21300, 23100, -3840, Cru },
{ 23100, 26160, -3840, Cen },
{ 24300, 24600, -3900, Cir },
{ 30150, 30300, -3900, Ara },
{ 3900, 5760, -4050, Hor },
{ 5760, 8250, -4050, Ret },
{ 26550, 26850, -4050, Cir },
{ 30300, 31500, -4050, Ara },
{ 31500, 32400, -4050, Pav },
{ 39600, 42000, -4050, Tuc },
{ 8250, 11850, -4200, Dor },
{ 24600, 26550, -4200, Cir },
{ 26550, 30600, -4200, TrA },
{ 0, 2400, -4500, Tuc },
{ 6300, 8250, -4500, Hyi },
{ 11850, 16260, -4500, Vol },
{ 16260, 20250, -4500, Car },
{ 20250, 24600, -4500, Mus },
{ 32400, 38400, -4500, Pav },
{ 38400, 42000, -4500, Ind },
{ 42000, 43200, -4500, Tuc },
{ 1350, 2400, -4560, Tuc },
{ 0, 6300, -4950, Hyi },
{ 13800, 24600, -4950, Cha },
{ 24600, 32400, -4950, Aps },
{ 6300, 13800, -5100, Men },
{ 0, 43200, -5400, Oct }
};
static short start[] = {
355, 352, 343, 340, 332, 320,
303, 288, 277, 266, 257, 251,
239, 229, 221, 211, 203, 189,
177, 163, 149, 136, 124, 104,
87, 75, 69, 54, 43, 29,
23, 16, 12, 6, 4, 0, 0
};
/*
======================================================================
constellation_pick()
Do a constellation pick from RA and Dec.
INPUTS
r right ascension, radians
d declination, radians
e epoch to which r and d precessed, as an mjd
RESULTS
Returns an index for the constellation region that the coordinates
belong to, or -1 if no constellation pick can be found.
The constellation is identified by linear search. We look for a
member of cbound[] whose lower_dec is less than the declination of the
pick point and whose lower_ra and upper_ra bracket the pick point's
right ascension. The data for that cbound[] member describes a sub-
region belonging to the picked constellation.
In geometric terms, the search amounts to starting at the north pole
and traveling south on a line of constant right ascension through the
pick point. At every declination where a cbound[] member lives, we
wake up, and if (1) the pick point is now behind us, and (2) the
cbound[] segment is crossing our path, we know we've hit the edge of
a constellation subregion that the pick point belongs to.
The cbound[] right ascension and declination values are scaled into
integers; this conserves storage and makes comparisons faster. The
start[] array, which gives starting points in cbound[] for different
declinations in 5-degree increments, further speeds the search by
skipping early parts of the list for which we know the search can't
succeed--geometrically, we start no more than 5 degrees north of the
pick point, rather than at the north pole.
The data in cbound[] are for epoch 1875.
======================================================================
*/
int
cns_pick(double r, double d, double e)
{
double Mjd;
unsigned short ra;
short de, i;
cal_mjd( 1, 1.0, 1875, &Mjd );
precess( e, Mjd, &r, &d );
ra = ( unsigned short )( radhr( r ) * 1800 );
de = ( short )( raddeg( d ) * 60 );
if (d < 0.0) --de;
i = ( de + 5400 ) / 300;
if ( i < 0 || i > 36 ) return -1;
i = start[ i ];
for ( ; i < NBOUNDS; i++ )
if ( cbound[ i ].lower_dec <= de &&
cbound[ i ].upper_ra > ra &&
cbound[ i ].lower_ra <= ra ) break;
return ( i == NBOUNDS ) ? -1 : ( int ) cbound[ i ].index;
}
/* given a constellation id (as from cns_pick()), return pointer to static
* storage containg its name in the form "AAA: Name".
* return "???: ???" if id is invalid.
*/
char *
cns_name (int id)
{
if (id < 0 || id >= NCNS)
return ("???: ???");
return (cns_namemap[id]);
}
/* return cns_namemap index matching first three chars in abbrev[], else -1.
*/
int
cns_id (char *abbrev)
{
int i;
for (i = 0; i < NCNS; i++)
if (strncmp (abbrev, cns_namemap[i], 3) == 0)
return (i);
return (-1);
}
/* edges of constant ra */
static struct {
unsigned short ra; /* hours * 1800 */
short dec0, dec1; /* degrees * 60 */
} ra_edges[] = {
{ 0, -4950, -4500 },
{ 0, 600, 750 },
{ 0, 1680, 1879 },
{ 120, 1320, 1680 },
{ 255, 750, 1260 },
{ 255, 1260, 1320 },
{ 300, 2760, 2880 },
{ 599, -420, 0 },
{ 599, 0, 120 },
{ 599, 3960, 4620 },
{ 1290, 1425, 1980 },
{ 1350, -4500, -4560 },
{ 1530, 1260, 1425 },
{ 1560, 2760, 2880 },
{ 2010, 2880, 3000 },
{ 2399, -4560, -3510 },
{ 2399, -3510, -3210 },
{ 2460, 3000, 3240 },
{ 2534, 1680, 1980 },
{ 2534, 1980, 2100 },
{ 2849, -3210, -3090 },
{ 3000, -2400, -1530 },
{ 3000, -1530, -1462 },
{ 3000, 595, 1500 },
{ 3000, 1500, 1680 },
{ 3000, 2820, 3000 },
{ 3060, 3240, 3450 },
{ 3299, -3090, -2890 },
{ 3434, 3450, 3510 },
{ 3450, 1500, 1635 },
{ 3600, 120, 595 },
{ 3600, 2100, 2205 },
{ 3675, 2820, 3030 },
{ 3900, -4050, -3510 },
{ 3900, -3510, -3240 },
{ 4199, -2890, -2400 },
{ 4350, -3240, -3060 },
{ 4350, 1635, 1840 },
{ 4379, 3420, 3510 },
{ 4530, 2205, 3030 },
{ 4620, 2040, 2205 },
{ 4770, -1462, -105 },
{ 4800, -3060, -2940 },
{ 4890, 1840, 2040 },
{ 5400, -2940, -2760 },
{ 5400, -2400, -2374 },
{ 5580, 3420, 4080 },
{ 5700, 3300, 3420 },
{ 5760, -3450, -4050 },
{ 5909, -105, 0 },
{ 5909, 0, 595 },
{ 5909, 595, 1140 },
{ 5999, 3150, 3300 },
{ 6060, 1140, 1840 },
{ 6150, -2760, -2640 },
{ 6150, 4080, 4620 },
{ 6300, -5100, -4950 },
{ 6300, -4950, -4500 },
{ 6300, -3190, -3450 },
{ 6300, -2374, -2160 },
{ 6314, 4620, 4800 },
{ 6449, -105, 0 },
{ 6750, -2160, -1462 },
{ 6899, -3060, -3190 },
{ 6960, -2640, -2400 },
{ 7200, -3390, -3190 },
{ 7349, -2940, -3060 },
{ 7680, -2400, -2940 },
{ 7680, -2220, -2400 },
{ 7799, -3540, -3390 },
{ 8100, -3240, -2940 },
{ 8100, -2940, -2790 },
{ 8100, 1800, 1840 },
{ 8100, 1840, 2160 },
{ 8249, -4500, -4200 },
{ 8249, -4050, -4200 },
{ 8249, -4050, -3540 },
{ 8249, -1800, -2220 },
{ 8310, 0, 930 },
{ 8400, -240, 0 },
{ 8445, 2160, 3150 },
{ 8460, -1635, -1800 },
{ 8550, 1710, 1800 },
{ 8699, -2790, -2580 },
{ 8699, -870, -1635 },
{ 8850, -660, -870 },
{ 8940, 930, 960 },
{ 9000, -3450, -3240 },
{ 9000, -2580, -1635 },
{ 9000, 3150, 3360 },
{ 9000, 4800, 5100 },
{ 9149, -660, -240 },
{ 9599, 930, 960 },
{ 9900, -3660, -3450 },
{ 10080, 750, 930 },
{ 10260, 1080, 1369 },
{ 10380, 750, 1080 },
{ 10499, -660, -240 },
{ 10589, 1290, 1369 },
{ 10589, 1369, 1680 },
{ 10589, 1680, 1710 },
{ 10800, -3840, -3660 },
{ 10800, -3150, -3045 },
{ 10800, -3045, -2580 },
{ 10980, 3240, 3360 },
{ 10980, 3360, 3720 },
{ 11010, -1980, -1635 },
{ 11010, -1635, -660 },
{ 11100, -3300, -3150 },
{ 11190, 1050, 1290 },
{ 11235, -240, 0 },
{ 11235, 0, 600 },
{ 11354, 600, 720 },
{ 11354, 720, 1050 },
{ 11700, -3480, -3300 },
{ 11700, 3000, 3240 },
{ 11759, 1680, 2130 },
{ 11849, -4200, -4500 },
{ 11849, -4200, -3840 },
{ 11849, -1980, -2580 },
{ 12240, 2670, 3000 },
{ 12299, -3840, -3480 },
{ 12479, 600, 720 },
{ 12600, 330, 600 },
{ 12600, 600, 750 },
{ 12600, 3600, 3720 },
{ 12630, 90, 330 },
{ 12960, 0, 90 },
{ 13260, -660, -1980 },
{ 13260, 2130, 2670 },
{ 13500, 750, 810 },
{ 13800, -5100, -4950 },
{ 13800, -4500, -4950 },
{ 13950, 2010, 2130 },
{ 14054, 600, 810 },
{ 14054, 810, 1200 },
{ 14189, 1200, 1680 },
{ 14265, 420, 600 },
{ 14340, 3600, 4410 },
{ 14400, -3045, -2580 },
{ 14400, 1680, 2010 },
{ 14400, 5100, 5190 },
{ 14400, 5280, 5190 },
{ 14549, -660, 0 },
{ 14549, 0, 420 },
{ 14700, -3180, -3045 },
{ 15060, -2580, -2205 },
{ 15060, -2205, -1020 },
{ 15060, -1020, -660 },
{ 15150, 2820, 3600 },
{ 15209, -3270, -3180 },
{ 15449, -1140, -1020 },
{ 15899, -3390, -3270 },
{ 16259, -4500, -3840 },
{ 16349, -1440, -1140 },
{ 16500, 2520, 2820 },
{ 16500, 4410, 4920 },
{ 16650, 420, 2010 },
{ 16650, 2010, 2385 },
{ 16860, -2385, -2205 },
{ 16860, -1440, -2205 },
{ 17249, -660, 0 },
{ 17249, 0, 420 },
{ 17249, 2385, 2520 },
{ 17550, -1590, -1440 },
{ 17789, 1710, 2010 },
{ 18300, 2400, 2520 },
{ 18450, -1750, -1590 },
{ 18900, 1410, 1710 },
{ 19049, -1870, -1750 },
{ 19200, 4800, 4920 },
{ 19350, -660, -1140 },
{ 19350, -360, -660 },
{ 19350, -360, 0 },
{ 19350, 0, 420 },
{ 19350, 1410, 1530 },
{ 19409, 2040, 2400 },
{ 19499, -2100, -1870 },
{ 19499, -1140, -1470 },
{ 19800, -3390, -2385 },
{ 19800, -2385, -2100 },
{ 19800, 1530, 1740 },
{ 19800, 1740, 2040 },
{ 20250, -4500, -3840 },
{ 20250, -3840, -3390 },
{ 20399, 3990, 4410 },
{ 20700, 4620, 4800 },
{ 20730, -360, 0 },
{ 20730, 0, 660 },
{ 21299, -3840, -3300 },
{ 21299, -1470, -660 },
{ 21299, -660, -360 },
{ 21360, 660, 840 },
{ 21360, 840, 1740 },
{ 21600, 1740, 2040 },
{ 21600, 2040, 2700 },
{ 21600, 3840, 3990 },
{ 21749, 2700, 3180 },
{ 22050, -1980, -2100 },
{ 22199, 1920, 2040 },
{ 22649, -1770, -1980 },
{ 22649, -1470, -1320 },
{ 23099, -3840, -3300 },
{ 23099, -1320, -660 },
{ 23099, 840, 900 },
{ 23400, 4620, 4200 },
{ 23850, 1710, 1920 },
{ 24300, -3900, -3840 },
{ 24300, 480, 900 },
{ 24300, 900, 1710 },
{ 24300, 2910, 3180 },
{ 24300, 3780, 3840 },
{ 24449, 4800, 4620 },
{ 24600, -4950, -4500 },
{ 24600, -4500, -4200 },
{ 24600, -4200, -3900 },
{ 25124, 1710, 1845 },
{ 25200, 4200, 3960 },
{ 25259, 1845, 2910 },
{ 25259, 2910, 3330 },
{ 25500, -3300, -2520 },
{ 25650, -1320, -1470 },
{ 25650, -480, -1320 },
{ 25950, 3330, 3780 },
{ 26100, 5190, 4800 },
{ 26159, -3840, -3300 },
{ 26400, 0, -480 },
{ 26550, -4200, -4050 },
{ 26850, -4050, -3814 },
{ 26850, -2520, -1770 },
{ 26850, -1470, -1770 },
{ 27090, -3300, -3240 },
{ 27149, -195, 0 },
{ 27149, 0, 480 },
{ 27149, 480, 1560 },
{ 27300, -3814, -3660 },
{ 27329, 1560, 1980 },
{ 27450, 3180, 3330 },
{ 27599, -3660, -3600 },
{ 27599, -3600, -3300 },
{ 27599, -3240, -2880 },
{ 27779, 1980, 2400 },
{ 28200, -2880, -2520 },
{ 28200, -1770, -1200 },
{ 28200, 3960, 4200 },
{ 28350, 2400, 3090 },
{ 28350, 3090, 3180 },
{ 28650, -1200, -480 },
{ 28650, -480, -195 },
{ 28650, 960, 1320 },
{ 28800, -2520, -1770 },
{ 28859, 1320, 1560 },
{ 28949, 240, 960 },
{ 29100, 1560, 1620 },
{ 29280, -1474, -1155 },
{ 29280, -1095, -480 },
{ 29280, -195, 0 },
{ 29280, 0, 240 },
{ 29399, 1620, 2400 },
{ 29475, -1155, -1095 },
{ 29557, -3660, -3600 },
{ 29557, -3600, -2730 },
{ 29557, -2730, -2520 },
{ 29759, 4200, 4500 },
{ 29849, -3814, -3660 },
{ 30150, -3900, -3814 },
{ 30150, -1800, -1474 },
{ 30150, 240, 769 },
{ 30299, -4050, -3900 },
{ 30600, -4050, -4200 },
{ 30600, 3030, 3090 },
{ 30900, -960, -600 },
{ 31050, 769, 859 },
{ 31500, -3420, -4050 },
{ 31500, 4500, 4800 },
{ 31649, -700, -600 },
{ 31680, -960, -1800 },
{ 31800, -700, -600 },
{ 32099, -2730, -2220 },
{ 32099, -1800, -2220 },
{ 32099, -240, 0 },
{ 32340, -600, -240 },
{ 32400, -4950, -4500 },
{ 32400, -4500, -4050 },
{ 32400, -2730, -3420 },
{ 32400, 4800, 5160 },
{ 32715, 1800, 2850 },
{ 32819, 2850, 3030 },
{ 32850, -240, -960 },
{ 32850, 0, 180 },
{ 32850, 270, 375 },
{ 32850, 720, 859 },
{ 33060, 1560, 1800 },
{ 33165, 180, 270 },
{ 33449, -240, 0 },
{ 33449, 0, 120 },
{ 33591, 375, 720 },
{ 33960, -960, -721 },
{ 33960, -721, -240 },
{ 33960, 120, 375 },
{ 33960, 720, 1110 },
{ 33960, 1110, 1264 },
{ 33960, 1264, 1530 },
{ 33960, 1530, 1560 },
{ 34200, 970, 1110 },
{ 34349, 2850, 3330 },
{ 34500, -2220, -2730 },
{ 34500, 2610, 2850 },
{ 34650, 1150, 1264 },
{ 34664, 1530, 1650 },
{ 34664, 1650, 1800 },
{ 34844, 1800, 2190 },
{ 34920, 2190, 2610 },
{ 34950, 3330, 3480 },
{ 35400, 1650, 1740 },
{ 35580, 3480, 3570 },
{ 35699, 945, 970 },
{ 35699, 1150, 1275 },
{ 36000, -1680, -721 },
{ 36000, -540, -721 },
{ 36000, 3570, 3690 },
{ 36255, 510, 945 },
{ 36300, 4500, 4800 },
{ 36450, 945, 1230 },
{ 36450, 1230, 1275 },
{ 36540, 120, 510 },
{ 36599, -3600, -3420 },
{ 36599, -3420, -2730 },
{ 36599, -2730, -1680 },
{ 36750, 3690, 4020 },
{ 36959, -900, -540 },
{ 36959, -540, 0 },
{ 36959, 0, 120 },
{ 36966, 3570, 3655 },
{ 37020, 1170, 1230 },
{ 37080, 3289, 3655 },
{ 37200, 4020, 4500 },
{ 37499, 120, 360 },
{ 37575, 360, 709 },
{ 37650, 1680, 1740 },
{ 37800, 4800, 5160 },
{ 37800, 5160, 5169 },
{ 37890, 709, 1170 },
{ 38010, 709, 750 },
{ 38250, 1170, 1410 },
{ 38399, -4500, -3600 },
{ 38399, -3000, -2730 },
{ 38399, -2220, -2730 },
{ 38399, -1680, -2220 },
{ 38399, -1680, -1530 },
{ 38399, -540, -900 },
{ 38399, 120, 750 },
{ 38550, 1410, 1680 },
{ 38640, 120, 165 },
{ 39000, 105, 165 },
{ 39119, 1680, 2160 },
{ 39360, -1530, -540 },
{ 39375, 2160, 2625 },
{ 39434, 2625, 2640 },
{ 39540, 2640, 3165 },
{ 39540, 3165, 3289 },
{ 39600, -4050, -3420 },
{ 39600, -3420, -3000 },
{ 39600, 105, 120 },
{ 39600, 2100, 2160 },
{ 39839, 3165, 3300 },
{ 40170, 3300, 3375 },
{ 40950, -240, 0 },
{ 40950, 0, 120 },
{ 40950, 120, 450 },
{ 41070, 2070, 2100 },
{ 41160, 2070, 3150 },
{ 41160, 3150, 3375 },
{ 41160, 3375, 3544 },
{ 41400, -2220, -1530 },
{ 41400, 5169, 5280 },
{ 41700, 3544, 3780 },
{ 41999, -4500, -4050 },
{ 41999, -3510, -3420 },
{ 41999, -3420, -2400 },
{ 41999, -2400, -2220 },
{ 41999, 3000, 3150 },
{ 42300, 1924, 2070 },
{ 42449, 2880, 3000 },
{ 42449, 3780, 3960 },
{ 42750, 1879, 1924 },
{ 42899, -1530, -420 },
{ 42899, -420, -240 },
{ 42899, 450, 600 },
};
#define NRA ((int)(sizeof(ra_edges)/sizeof(ra_edges[0])))
/* edges of constant dec */
static struct {
short dec; /* degrees * 60 */
unsigned short ra0, ra1; /* hours * 1800 */
} dec_edges[] = {
{ -5100, 6300, 13800 },
{ -4950, 0, 6300 },
{ -4950, 13800, 24600 },
{ -4950, 24600, 32400 },
{ -4560, 1350, 2399 },
{ -4500, 0, 1350 },
{ -4500, 6300, 8249 },
{ -4500, 11849, 13800 },
{ -4500, 13800, 16259 },
{ -4500, 16259, 20250 },
{ -4500, 20250, 24600 },
{ -4500, 32400, 38399 },
{ -4500, 38399, 41999 },
{ -4500, 41999, 0 },
{ -4200, 8249, 11849 },
{ -4200, 24600, 26550 },
{ -4200, 30600, 26550 },
{ -4050, 3900, 5760 },
{ -4050, 5760, 8249 },
{ -4050, 26550, 26850 },
{ -4050, 30299, 30600 },
{ -4050, 31500, 30600 },
{ -4050, 31500, 32400 },
{ -4050, 39600, 41999 },
{ -3900, 24300, 24600 },
{ -3900, 30150, 30299 },
{ -3840, 10800, 11849 },
{ -3840, 11849, 12299 },
{ -3840, 12299, 16259 },
{ -3840, 20250, 21299 },
{ -3840, 21299, 23099 },
{ -3840, 23099, 24300 },
{ -3840, 24300, 26159 },
{ -3814, 26850, 27300 },
{ -3814, 29849, 30150 },
{ -3660, 9900, 10800 },
{ -3660, 27300, 27599 },
{ -3660, 29557, 29849 },
{ -3600, 27599, 29557 },
{ -3600, 36599, 38399 },
{ -3540, 7799, 8249 },
{ -3510, 2399, 3900 },
{ -3510, 2399, 41999 },
{ -3510, 3900, 2399 },
{ -3510, 41999, 2399 },
{ -3480, 11700, 12299 },
{ -3450, 6300, 5760 },
{ -3450, 9000, 9900 },
{ -3420, 32400, 31500 },
{ -3420, 32400, 36599 },
{ -3420, 39600, 41999 },
{ -3390, 7200, 7799 },
{ -3390, 15899, 19800 },
{ -3390, 19800, 20250 },
{ -3300, 11100, 11700 },
{ -3300, 21299, 23099 },
{ -3300, 25500, 26159 },
{ -3300, 26159, 27090 },
{ -3300, 27090, 27599 },
{ -3270, 15209, 15899 },
{ -3240, 3900, 4350 },
{ -3240, 8100, 9000 },
{ -3240, 27090, 27599 },
{ -3210, 2399, 2849 },
{ -3190, 6899, 6300 },
{ -3190, 6899, 7200 },
{ -3180, 14700, 15209 },
{ -3150, 10800, 11100 },
{ -3090, 2849, 3299 },
{ -3060, 4350, 4800 },
{ -3060, 7349, 6899 },
{ -3045, 10800, 14400 },
{ -3045, 14400, 14700 },
{ -3000, 38399, 39600 },
{ -2940, 4800, 5400 },
{ -2940, 7680, 7349 },
{ -2940, 7680, 8100 },
{ -2890, 3299, 4199 },
{ -2880, 27599, 28200 },
{ -2790, 8100, 8699 },
{ -2760, 5400, 6150 },
{ -2730, 29557, 32099 },
{ -2730, 32099, 32400 },
{ -2730, 34500, 32400 },
{ -2730, 34500, 36599 },
{ -2730, 38399, 36599 },
{ -2640, 6150, 6960 },
{ -2580, 8699, 9000 },
{ -2580, 10800, 9000 },
{ -2580, 11849, 10800 },
{ -2580, 14400, 15060 },
{ -2520, 25500, 26850 },
{ -2520, 28200, 28800 },
{ -2520, 29557, 28800 },
{ -2400, 3000, 41999 },
{ -2400, 4199, 3000 },
{ -2400, 4199, 5400 },
{ -2400, 6960, 7680 },
{ -2385, 16860, 19800 },
{ -2374, 5400, 6300 },
{ -2220, 8249, 7680 },
{ -2220, 32099, 34500 },
{ -2220, 38399, 41400 },
{ -2220, 41999, 41400 },
{ -2205, 16860, 15060 },
{ -2160, 6300, 6750 },
{ -2100, 19800, 19499 },
{ -2100, 22050, 19800 },
{ -1980, 11010, 11849 },
{ -1980, 13260, 11849 },
{ -1980, 22649, 22050 },
{ -1870, 19499, 19049 },
{ -1800, 8460, 8249 },
{ -1800, 30150, 31680 },
{ -1800, 31680, 32099 },
{ -1770, 26850, 22649 },
{ -1770, 26850, 28200 },
{ -1770, 28800, 28200 },
{ -1750, 19049, 18450 },
{ -1680, 36000, 36599 },
{ -1680, 36599, 38399 },
{ -1635, 8699, 8460 },
{ -1635, 8699, 9000 },
{ -1635, 9000, 11010 },
{ -1590, 18450, 17550 },
{ -1530, 3000, 42899 },
{ -1530, 38399, 39360 },
{ -1530, 41400, 39360 },
{ -1530, 41400, 42899 },
{ -1530, 42899, 3000 },
{ -1530, 42899, 41400 },
{ -1474, 29280, 30150 },
{ -1470, 19499, 21299 },
{ -1470, 21299, 22649 },
{ -1470, 25650, 26850 },
{ -1462, 3000, 4770 },
{ -1462, 4770, 6750 },
{ -1440, 16349, 16860 },
{ -1440, 17550, 16860 },
{ -1320, 22649, 23099 },
{ -1320, 23099, 25650 },
{ -1200, 28200, 28650 },
{ -1155, 29280, 29475 },
{ -1140, 15449, 16349 },
{ -1140, 19350, 19499 },
{ -1095, 29280, 29475 },
{ -1020, 15060, 15449 },
{ -960, 30900, 31680 },
{ -960, 32850, 31680 },
{ -960, 32850, 33960 },
{ -900, 38399, 36959 },
{ -870, 8850, 8699 },
{ -721, 36000, 33960 },
{ -700, 31649, 31800 },
{ -660, 9149, 8850 },
{ -660, 10499, 9149 },
{ -660, 11010, 10499 },
{ -660, 11010, 13260 },
{ -660, 14549, 13260 },
{ -660, 15060, 14549 },
{ -660, 17249, 19350 },
{ -660, 23099, 21299 },
{ -600, 30900, 31649 },
{ -600, 31800, 32340 },
{ -540, 36959, 36000 },
{ -540, 39360, 38399 },
{ -480, 26400, 25650 },
{ -480, 28650, 29280 },
{ -420, 599, 42899 },
{ -360, 20730, 19350 },
{ -360, 21299, 20730 },
{ -240, 8400, 9149 },
{ -240, 10499, 11235 },
{ -240, 32099, 32340 },
{ -240, 33449, 32850 },
{ -240, 33960, 33449 },
{ -240, 40950, 42899 },
{ -195, 28650, 27149 },
{ -195, 28650, 29280 },
{ -105, 4770, 5909 },
{ -105, 5909, 6449 },
{ 0, 6449, 8310 },
{ 0, 8310, 8400 },
{ 0, 12960, 14549 },
{ 0, 27149, 26400 },
{ 0, 32099, 32850 },
{ 90, 12630, 12960 },
{ 105, 39000, 39600 },
{ 120, 599, 3600 },
{ 120, 33449, 33960 },
{ 120, 36540, 36959 },
{ 120, 36959, 37499 },
{ 120, 37499, 38399 },
{ 120, 38399, 38640 },
{ 120, 39600, 40950 },
{ 165, 38640, 39000 },
{ 180, 32850, 33165 },
{ 240, 28949, 29280 },
{ 240, 29280, 30150 },
{ 270, 33165, 32850 },
{ 330, 12600, 12630 },
{ 360, 37499, 37575 },
{ 375, 32850, 33591 },
{ 375, 33591, 33960 },
{ 420, 14265, 14549 },
{ 420, 14549, 16650 },
{ 420, 16650, 17249 },
{ 420, 17249, 19350 },
{ 450, 40950, 42899 },
{ 480, 24300, 27149 },
{ 510, 36255, 36540 },
{ 595, 3000, 3600 },
{ 595, 3600, 5909 },
{ 600, 0, 42899 },
{ 600, 11235, 11354 },
{ 600, 12479, 12600 },
{ 600, 14054, 14265 },
{ 600, 42899, 0 },
{ 660, 20730, 21360 },
{ 709, 37575, 37890 },
{ 709, 37890, 38010 },
{ 720, 11354, 12479 },
{ 720, 32850, 33591 },
{ 720, 33591, 33960 },
{ 750, 0, 255 },
{ 750, 10080, 10380 },
{ 750, 12600, 13500 },
{ 750, 38010, 38399 },
{ 769, 30150, 31050 },
{ 810, 13500, 14054 },
{ 840, 21360, 23099 },
{ 859, 31050, 32850 },
{ 900, 23099, 24300 },
{ 930, 8310, 8940 },
{ 930, 9599, 10080 },
{ 945, 35699, 36255 },
{ 945, 36255, 36450 },
{ 960, 8940, 9599 },
{ 960, 28650, 28949 },
{ 970, 34200, 35699 },
{ 1050, 11190, 11354 },
{ 1080, 10260, 10380 },
{ 1110, 33960, 34200 },
{ 1140, 5909, 6060 },
{ 1150, 34650, 35699 },
{ 1170, 37020, 37890 },
{ 1170, 37890, 38250 },
{ 1200, 14054, 14189 },
{ 1230, 36450, 37020 },
{ 1260, 255, 1530 },
{ 1264, 33960, 34650 },
{ 1275, 35699, 36450 },
{ 1290, 10589, 11190 },
{ 1320, 120, 255 },
{ 1320, 28650, 28859 },
{ 1369, 10260, 10589 },
{ 1410, 18900, 19350 },
{ 1410, 38250, 38550 },
{ 1425, 1290, 1530 },
{ 1500, 3000, 3450 },
{ 1530, 19350, 19800 },
{ 1530, 33960, 34664 },
{ 1560, 27149, 27329 },
{ 1560, 27329, 28859 },
{ 1560, 28859, 29100 },
{ 1560, 33060, 33960 },
{ 1620, 29100, 29399 },
{ 1635, 3450, 4350 },
{ 1650, 34664, 35400 },
{ 1680, 0, 120 },
{ 1680, 2534, 3000 },
{ 1680, 10589, 11759 },
{ 1680, 14189, 14400 },
{ 1680, 37650, 38550 },
{ 1680, 38550, 39119 },
{ 1710, 8550, 10589 },
{ 1710, 17789, 18900 },
{ 1710, 23850, 24300 },
{ 1710, 24300, 25124 },
{ 1740, 19800, 21360 },
{ 1740, 21360, 21600 },
{ 1740, 35400, 37650 },
{ 1800, 8100, 8550 },
{ 1800, 32715, 33060 },
{ 1800, 34664, 34844 },
{ 1840, 4350, 4890 },
{ 1840, 4890, 6060 },
{ 1840, 6060, 8100 },
{ 1845, 25124, 25259 },
{ 1879, 0, 42750 },
{ 1920, 22199, 23850 },
{ 1924, 42300, 42750 },
{ 1980, 1290, 2534 },
{ 1980, 27329, 27779 },
{ 2010, 13950, 14400 },
{ 2010, 14400, 16650 },
{ 2010, 16650, 17789 },
{ 2040, 4620, 4890 },
{ 2040, 19409, 19800 },
{ 2040, 21600, 22199 },
{ 2070, 41070, 41160 },
{ 2070, 41160, 42300 },
{ 2100, 2534, 3600 },
{ 2100, 39600, 41070 },
{ 2130, 11759, 13260 },
{ 2130, 13260, 13950 },
{ 2160, 8100, 8445 },
{ 2160, 39119, 39375 },
{ 2160, 39375, 39600 },
{ 2190, 34844, 34920 },
{ 2205, 3600, 4530 },
{ 2205, 4530, 4620 },
{ 2385, 16650, 17249 },
{ 2400, 18300, 19409 },
{ 2400, 27779, 28350 },
{ 2400, 28350, 29399 },
{ 2520, 16500, 17249 },
{ 2520, 17249, 18300 },
{ 2610, 34500, 34920 },
{ 2625, 39375, 39434 },
{ 2640, 39434, 39540 },
{ 2670, 12240, 13260 },
{ 2700, 21600, 21749 },
{ 2760, 300, 1560 },
{ 2820, 3000, 3675 },
{ 2820, 15150, 16500 },
{ 2850, 32715, 32819 },
{ 2850, 32819, 34349 },
{ 2850, 34349, 34500 },
{ 2880, 300, 42449 },
{ 2880, 1560, 2010 },
{ 2880, 42449, 300 },
{ 2910, 24300, 25259 },
{ 3000, 2010, 2460 },
{ 3000, 2460, 3000 },
{ 3000, 11700, 12240 },
{ 3000, 41999, 42449 },
{ 3030, 3675, 4530 },
{ 3030, 30600, 32819 },
{ 3090, 28350, 30600 },
{ 3150, 5999, 8445 },
{ 3150, 8445, 9000 },
{ 3150, 41160, 41999 },
{ 3165, 39540, 39839 },
{ 3180, 21749, 24300 },
{ 3180, 27450, 28350 },
{ 3240, 2460, 3060 },
{ 3240, 10980, 11700 },
{ 3289, 37080, 39540 },
{ 3300, 5700, 5999 },
{ 3300, 39839, 40170 },
{ 3330, 25259, 25950 },
{ 3330, 25950, 27450 },
{ 3330, 34349, 34950 },
{ 3360, 9000, 10980 },
{ 3375, 40170, 41160 },
{ 3420, 4379, 5580 },
{ 3420, 5580, 5700 },
{ 3450, 3060, 3434 },
{ 3480, 34950, 35580 },
{ 3510, 3434, 4379 },
{ 3544, 41160, 41700 },
{ 3570, 35580, 36000 },
{ 3570, 36000, 36966 },
{ 3600, 12600, 14340 },
{ 3600, 14340, 15150 },
{ 3655, 36966, 37080 },
{ 3690, 36000, 36750 },
{ 3720, 10980, 12600 },
{ 3780, 24300, 25950 },
{ 3780, 41700, 42449 },
{ 3840, 21600, 24300 },
{ 3960, 599, 42449 },
{ 3960, 25200, 28200 },
{ 3960, 42449, 599 },
{ 3990, 20399, 21600 },
{ 4020, 36750, 37200 },
{ 4080, 5580, 6150 },
{ 4200, 23400, 25200 },
{ 4200, 28200, 29759 },
{ 4410, 14340, 16500 },
{ 4410, 16500, 20399 },
{ 4500, 29759, 31500 },
{ 4500, 36300, 37200 },
{ 4620, 599, 6150 },
{ 4620, 6150, 6314 },
{ 4620, 20700, 23400 },
{ 4620, 24449, 23400 },
{ 4800, 6314, 9000 },
{ 4800, 19200, 20700 },
{ 4800, 26100, 24449 },
{ 4800, 31500, 32400 },
{ 4800, 36300, 37800 },
{ 4920, 16500, 19200 },
{ 5100, 9000, 14400 },
{ 5160, 32400, 37800 },
{ 5169, 37800, 41400 },
{ 5190, 14400, 26100 },
{ 5280, 6300, 14400 },
{ 5280, 41400, 6300 },
};
#define NDEC ((int)(sizeof(dec_edges)/sizeof(dec_edges[0])))
/* given an epoch, give caller a list of all constellation edges.
* return count if ok, else -1.
* N.B. caller should *not* free what we return because we cache it here.
*/
int
cns_edges (double e, double **ra0p, double **dec0p, double **ra1p,
double **dec1p)
{
#define NEDGES (NRA+NDEC)
static double *ra0, *dec0, *ra1, *dec1;
static double laste = -12345.6; /* any bogus value */
double mjd0;
int i, n;
/* if same epoch just return the same list */
if (e == laste) {
*ra0p = ra0;
*dec0p = dec0;
*ra1p = ra1;
*dec1p = dec1;
return (NEDGES);
}
/* get space for arrays, first time only */
if (!ra0) {
ra0 = (double *)malloc (NEDGES * sizeof(double));
if (!ra0)
return (-1);
dec0 = (double *)malloc (NEDGES * sizeof(double));
if (!dec0) {
free ((void *)ra0);
return (-1);
}
ra1 = (double *)malloc (NEDGES * sizeof(double));
if (!ra1) {
free ((void *)ra0);
free ((void *)dec0);
return (-1);
}
dec1 = (double *)malloc (NEDGES * sizeof(double));
if (!dec1) {
free ((void *)ra0);
free ((void *)dec0);
free ((void *)ra1);
return (-1);
}
}
/* prepare for precession from 1875 */
cal_mjd (1, 1.0, 1875, &mjd0);
/* build the constant-ra edge lists */
n = 0;
for (i = 0; i < NRA; i++) {
ra0[n] = ra1[n] = hrrad((double)ra_edges[i].ra/1800.0);
dec0[n] = degrad((double)ra_edges[i].dec0/60.0);
dec1[n] = degrad((double)ra_edges[i].dec1/60.0);
precess (mjd0, e, &ra0[n], &dec0[n]);
precess (mjd0, e, &ra1[n], &dec1[n]);
n++;
}
/* add the constant-dec edge lists */
for (i = 0; i < NDEC; i++) {
ra0[n] = hrrad((double)dec_edges[i].ra0/1800.0);
ra1[n] = hrrad((double)dec_edges[i].ra1/1800.0);
dec0[n] = dec1[n] = degrad((double)dec_edges[i].dec/60.0);
precess (mjd0, e, &ra0[n], &dec0[n]);
precess (mjd0, e, &ra1[n], &dec1[n]);
n++;
}
/* sanity check the count */
if (n != NEDGES) {
printf ("cns_edges(): n=%d NEDGES=%ld\n", n, (long)NEDGES);
abort();
}
/* ok */
*ra0p = ra0;
*dec0p = dec0;
*ra1p = ra1;
*dec1p = dec1;
laste = e;
return (NEDGES);
}
/* given an ra, dec and epoch return the list of constellation ids which
* *may* fall within the given radius of said location.
* return the number of ids.
* ids[] need be no larger than 89.
*/
/* ARGSUSED */
int
cns_list (double ra, double dec, double e, double rad, int ids[])
{
int i;
/* TODO: this! */
for (i = 0; i < NCNS; i++)
ids[i] = i;
return (NCNS);
}
/* epoch 2000 RA/Dec of constellation figure end-points.
* drawcodes: 0=move to; 1=draw to; 2=draw to dashed; -1=end
*/
typedef struct {
int drawcode; /* draw code */
float ra; /* rads */
float dec; /* rads */
} ConFig;
/* array of malloced lists of ConFigs, same order as cns_namemap[]
*/
static ConFig *figmap[NCNS];
/* add one entry to the drawing code lists */
static void
addFigList (ConFig **new, int *nused, int c, int drawcode, double ra, double dec)
{
ConFig *cp;
new[c]= (ConFig*) realloc (new[c], (nused[c]+1)*sizeof(ConFig));
cp = &new[c][nused[c]++];
cp->drawcode = drawcode;
cp->ra = (float)hrrad(ra);
cp->dec = (float)degrad(dec);
}
/* load the given constellation definition file.
* return 0 if ok else reason why not in msg[] and -1.
*/
int
cns_loadfigs (FILE *fp, char *msg)
{
char line[1024]; /* one line from the file */
char cname[1024]; /* constellation name */
ConFig **new; /* array of ConFig[] for each cnstn */
int *nused; /* number of ConFig[] for each cnstn */
int c = -1; /* index, same as cns_namemap[] */
int s = 0; /* status */
/* init the temp lists */
new = (ConFig **) calloc (NCNS, sizeof(ConFig*));
nused = (int *) calloc (NCNS, sizeof(int));
/* read the file */
while (fgets (line, sizeof(line), fp)) {
char rastr[64], decstr[64];
char *lp;
int code;
/* skip leading/trailing whitespace, blank lines and # lines */
for (lp = line+strlen(line)-1; lp>=line && isspace(*lp); --lp)
*lp = '\0';
for (lp = line; isspace(*lp); lp++)
continue;
if (*lp == '#' || *lp == '\0')
continue;
/* ok, line looks interesting, look more carefully */
if (sscanf (lp, "%d %63s %63s", &code, rastr, decstr) == 3) {
/* looks like a drawing line */
double ra, dec;
/* must be working on a current constellation */
if (c < 0) {
sprintf (msg,"Found coord line before first constellation");
s = -1;
break;
}
/* check draw code */
if (code < 0 || code > 2) {
sprintf (msg, "Bad draw code in %s: %d", cname, code);
s = -1;
break;
}
/* crack ra dec */
if (f_scansexa (rastr, &ra) < 0 || ra < 0 || ra >= 24) {
sprintf (msg, "Bad RA format in %s: %s", cname, rastr);
s = -1;
break;
}
if (f_scansexa (decstr, &dec) < 0 || dec < -90 || dec > 90) {
sprintf (msg, "Bad Dec format in %s: %s", cname, decstr);
s = -1;
break;
}
/* add to list */
addFigList (new, nused, c, code, ra, dec);
} else {
/* finish previous list, if any */
if (c >= 0)
addFigList (new, nused, c, -1, 0.0, 0.0);
/* see if it's a recognized constellation name */
for (c = 0; c < NCNS; c++)
if (strcmp (lp, cns_namemap[c]+5) == 0)
break;
if (c == NCNS) {
sprintf (msg, "Unknown constellation: %s", lp);
s = -1;
break;
}
if (new[c]) {
sprintf (msg, "Duplicate definition for %s", lp);
s = -1;
break;
}
/* init its list */
strcpy (cname, lp);
new[c] = (ConFig *) malloc (1); /* realloc seed */
}
}
/* even if ok check we found all definitions */
if (s == 0) {
int l = 0;
/* finish last list */
addFigList (new, nused, c, -1, 0.0, 0.0);
for (c = 0; c < NCNS; c++)
if (!new[c])
l += sprintf (msg+l, "%s ", cns_namemap[c]+5);
if (l > 0) {
strcat (msg, ": no definition found");
s = -1;
}
}
/* handle ok or error */
if (s < 0) {
/* trouble: free temp lists */
for (c = 0; c < NCNS; c++)
if (new[c])
free (new[c]);
} else {
/* make temp lists persistent */
for (c = 0; c < NCNS; c++) {
if (figmap[c])
free (figmap[c]);
figmap[c] = new[c];
}
}
/* done with lists themselves regardless */
free (new);
free (nused);
/* done */
return (s);
}
/* given a constellation id and epoch, return arrays of ra[] and dec[]
* end-points precessed to the desired epoch that, if connected, will form the
* given constellation figure.
* dcodes is 0 if the coord is a "move-to", 1 if a "draw-to" or 2 if a "draw-to
* as dotted-line".
* return the total number of tripples or -1 if id is bogus.
* the arrays need be no larger than 35 entries.
*/
int
cns_figure (int id, double e, double ra[], double dec[], int dcodes[])
{
ConFig *cfp;
if (id < 0 || id >= NCNS)
return (-1);
for (cfp = figmap[id]; cfp->drawcode >= 0; cfp++) {
*ra = (double)cfp->ra;
*dec = (double)cfp->dec;
precess (J2000, e, ra, dec);
ra++;
dec++;
*dcodes++ = cfp->drawcode;
}
return (cfp - figmap[id]);
}
Reference in New Issue View Git Blame Copy Permalink