ref: c00c60d327168c25e56077c8310cf2380833e9be
dir: /sys/src/cmd/astro/pdate.c/
#include "astro.h"
char* month[] =
{
"January",
"February",
"March",
"April",
"May",
"June",
"July",
"August",
"September",
"October",
"November",
"December",
};
double
dsrc(double d, Tim *t, int i)
{
double y;
do {
t->ifa[i] += 1.;
y = convdate(t);
} while(d >= y);
do {
t->ifa[i] -= 1.;
y = convdate(t);
} while(d < y);
return d - y;
}
void
dtsetup(double d, Tim *t)
{
double v;
t->ifa[0] = floor(1900 + d/365.24220);
t->ifa[1] = 1;
t->ifa[2] = 1;
t->ifa[3] = 0;
t->ifa[4] = 0;
t->ifa[1] = floor(1 + dsrc(d, t, 0)/30);
t->ifa[2] = floor(1 + dsrc(d, t, 1));
dsrc(d, t, 2);
v = (d - convdate(t)) * 24;
t->ifa[3] = floor(v);
t->ifa[4] = (v - t->ifa[3]) * 60;
convdate(t); /* to set timezone */
}
void
pdate(double d)
{
int i;
Tim t;
dtsetup(d, &t);
if(flags['s']) {
i = t.ifa[1];
print("%s ", month[i-1]);
i = t.ifa[2];
numb(i);
print("...");
return;
}
/* year month day */
print("%4d %2d %2d",
(int)t.ifa[0],
(int)t.ifa[1],
(int)t.ifa[2]);
}
void
ptime(double d)
{
int h, m, s;
char *mer;
Tim t;
if(flags['s']) {
/* hour minute */
dtsetup(d + .5/(24*60), &t);
h = t.ifa[3];
m = floor(t.ifa[4]);
mer = "AM";
if(h >= 12) {
mer = "PM";
h -= 12;
}
if(h == 0)
h = 12;
numb(h);
if(m < 10) {
if(m == 0) {
print("%s exactly ...", mer);
return;
}
print("O ");
}
numb(m);
print("%s ...", mer);
return;
}
/* hour minute second */
dtsetup(d, &t);
h = t.ifa[3];
m = floor(t.ifa[4]);
s = floor((t.ifa[4]-m) * 60);
print("%.2d:%.2d:%.2d %.*s", h, m, s, utfnlen(t.tz, 3), t.tz);
}
char* unit[] =
{
"zero",
"one",
"two",
"three",
"four",
"five",
"six",
"seven",
"eight",
"nine",
"ten",
"eleven",
"twelve",
"thirteen",
"fourteen",
"fifteen",
"sixteen",
"seventeen",
"eighteen",
"nineteen"
};
char* decade[] =
{
"twenty",
"thirty",
"forty",
"fifty",
"sixty",
"seventy",
"eighty",
"ninety"
};
void
pstime(double d)
{
setime(d);
semi = 0;
motion = 0;
rad = 1.e9;
lambda = 0;
beta = 0;
// uses lambda, beta, rad, motion
// sets alpha, delta, rp
helio();
// uses alpha, delta, rp
// sets ra, decl, lha, decl2, az, el
geo();
print(" %R %D %D %4.0f", lha, nlat, awlong, elev/3.28084);
}
void
numb(int n)
{
if(n >= 100) {
print("%d ", n);
return;
}
if(n >= 20) {
print("%s ", decade[n/10 - 2]);
n %= 10;
if(n == 0)
return;
}
print("%s ", unit[n]);
}
double
tzone(double y, Tim *z)
{
double t, l1, l2;
Tm t1, t2;
/*
* get a rough approximation to unix mean time
*/
t = (y - 25567.5) * 86400;
/*
* if outside unix conversions,
* just call it GMT
*/
if(t < 0 || t > 2.1e9)
return y;
/*
* convert by both local and gmt
*/
t1 = *localtime((long)t);
t2 = *gmtime((long)t);
/*
* pick up year crossings
*/
if(t1.yday == 0 && t2.yday > 1)
t1.yday = t2.yday+1;
if(t2.yday == 0 && t1.yday > 1)
t2.yday = t1.yday+1;
/*
* convert times to days
*/
l1 = t1.yday + t1.hour/24. + t1.min/1440. + t1.sec/86400.;
l2 = t2.yday + t2.hour/24. + t2.min/1440. + t2.sec/86400.;
/*
* return difference
*/
strncpy(z->tz, t1.zone, sizeof(z->tz));
return y + (l2 - l1);
}
int dmo[12] =
{
0,
31,
59,
90,
120,
151,
181,
212,
243,
273,
304,
334
};
/*
* input date conversion
* output is done by zero crossing
* on this input conversion.
*/
double
convdate(Tim *t)
{
double y, d;
int m;
y = t->ifa[0];
m = t->ifa[1];
d = t->ifa[2];
/*
* normalize the month
*/
while(m < 1) {
m += 12;
y -= 1;
}
while(m > 12) {
m -= 12;
y += 1;
}
/*
* bc correction
*/
if(y < 0)
y += 1;
/*
* normal conversion
*/
y += 4712;
if(fmod(y, 4) == 0 && m > 2)
d += 1;
y = y*365 + floor((y+3)/4) + dmo[m-1] + d - 1;
/*
* gregorian change
*/
if(y > 2361232)
y -= floor((y-1794167)/36524.220) -
floor((y-1721117)/146100);
y += t->ifa[3]/24 + t->ifa[4]/1440 - 2415020.5;
/*
* kitchen clock correction
*/
strncpy(t->tz, "GMT", sizeof(t->tz));
if(flags['k'])
y = tzone(y, t);
return y;
}