UNIX time_t

/*3:*/

#line 80 "timedate.w"

#ifndef TIMEDATE_HEADER_DEFINES

#define TIMEDATE_HEADER_DEFINES

#include <stdio.h> 

#include <math.h>        

#include <time.h> 

#include <assert.h> 

#include <string> 

#include <stdexcept> 

using namespace std;

/*6:*/

#line 117 "timedate.w"

class angle{

#define PI 3.14159265358979323846          

public:

static double Pi(void){

return PI;

}

static double Asec(double x){

return x/3600.0;

}

static double dtr(double d){

return d*(PI/180);

}

static double rtd(double r){

return r/(PI/180);

}

static double fixangle(double d){

return d-360.0*(floor(d/360.0));

}

static double fixangr(double r){

return r-(PI*2)*(floor(r/(PI*2)));

}

static double d_m_s_to_decimal(int d,int m,double s){

return d+(m/60.0)+(s/(60.0*60.0));

}

static string degrees_to_d_m_s(double a);

#undef PI

};

/*:6*//*8:*/

#line 209 "timedate.w"

class systemtime:protected angle{

private:

time_t t;

public:

static const int SecondsPerDay= 24*60*60,

SecondsPerHour= 60*60,

SecondsPerMinute= 60;

systemtime(time_t it= 0){

set_time(it);

}

systemtime(string dateTime){

fromString(dateTime);

}

time_t get_time(void){

return t;

}

void set_time(const time_t nt){

t= nt;

}

void now(void){

t= time(NULL);

}

time_t midnight(void){

return t-(t%SecondsPerDay);

}

void dayStep(int days){

t= ((t/SecondsPerDay)+days)*SecondsPerDay;

}

void nextDay(void){

dayStep(1);

}

void previousDay(void){

dayStep(-1);

}

void toUTC(int*year= NULL,int*month= NULL,int*mday= NULL,

int*hour= NULL,int*min= NULL,int*sec= NULL);

void fromUTC(int year= 1970,int month= 1,int mday= 1,

int hour= 0,int min= 0,int sec= 0);

string dateToString(void);

string timeToString(void);

void fromString(string s);

double toJulian(void);

void fromJulian(double jd);

int weekday(void){

return jwday(toJulian());

}

double siderealTime(void){

return gmst(toJulian());

}

static double utctoj(long year,int mon,int mday,

int hour= 0,int min= 0,int sec= 0);

static void jyear(double td,

int*yy= NULL,int*mm= NULL,int*dd= NULL);

static void jhms(double j,int*h= NULL,int*m= NULL,int*s= NULL);

static int jwday(double j){

return((int)(j+1.5))%7;

}

static double gmst(double jd);

static double meanObliquityOfEcliptic(double jd);

static void nutation(double jd,double*deltaPsi,double*deltaEpsilon);

static double equationTime(double jd);

static void sunpos(double jd,bool apparent,

double*ra,double*dec,double*rv,double*slong);

static void moon_and_sun(double pdate,double*phaseang,

double*pphase,double*mage,double*dist,

double*angdia,double*sudist,double*suangdia);

static double kepler(double m,double ecc);

static void moonphases(double sdate,double phases[5]);

protected:

static void store_if(int*p,const int v){

if(p!=NULL){

*p= v;

}

}

static void dstore_if(double*p,const double v){

if(p!=NULL){

*p= v;

}

}

static double meanphase(double sdate,double k);

static double truephase(double k,double phase);

};

/*:8*/

#line 92 "timedate.w"

#endif

/*:3*/

////////////////////////////////////////////////

/*2:*/

#line 67 "timedate.w"

#include "config.h"                   

#define JulianCentury 36525.0

#define JulianMillennium (JulianCentury*10) 

#define J2000 2451545.0

#define SynMonth 29.53058868 \

 \

#define NUTERMS 63 \

#line 70 "timedate.w"

/*4:*/

#line 99 "timedate.w"

#include "timedate.h"       

/*:4*/

#line 72 "timedate.w"

/*5:*/

#line 106 "timedate.w"

/*7:*/

#line 168 "timedate.w"

string angle::degrees_to_d_m_s(double a){

int dd,mm;

double ss;

char result[16];

char*sign= "";

if(a<0){

sign= "-";

a= -a;

}

dd= (int)a;

mm= (int)((a-dd)*60);

ss= (a-(dd+(mm/60.0)))*3600;

sprintf(result,"%s%d°%d'%.3f\"",sign,dd,mm,ss);

return string(result);

}

/*:7*/

#line 107 "timedate.w"

/*9:*/

#line 374 "timedate.w"

void systemtime::toUTC(int*year,int*month,int*mday,

int*hour,int*min,int*sec){

struct tm*gt= gmtime(&t);

store_if(year,gt->tm_year+1900);

store_if(month,gt->tm_mon+1);

store_if(mday,gt->tm_mday);

store_if(hour,gt->tm_hour);

store_if(min,gt->tm_min);

store_if(sec,gt->tm_sec);

}

/*:9*//*10:*/

#line 394 "timedate.w"

void systemtime::fromUTC(int year,int month,int mday,

int hour,int min,int sec){

double jd;

const double J1970= 2440587.5;

jd= utctoj(year,month,mday,hour,min,sec);

t= (time_t)(((jd-J1970)*(60*60*24))+0.5);

}

/*:10*//*11:*/

#line 409 "timedate.w"

string systemtime::dateToString(void){

int year,month,day;

char edate[16];

toUTC(&year,&month,&day);

sprintf(edate,"%04d-%02d-%02d",year,month,day);

return string(edate);

}

/*:11*//*12:*/

#line 423 "timedate.w"

string systemtime::timeToString(void){

int hour,min,sec;

char etime[16];

toUTC(NULL,NULL,NULL,&hour,&min,&sec);

sprintf(etime,"%02d:%02d:%02d",hour,min,sec);

return string(etime);

}

/*:12*//*13:*/

#line 439 "timedate.w"

void systemtime::fromString(string s)

{

int year= 1970,month= 1,mday= 1,

hour= 0,min= 0,sec= 0,n;

char d1,d2,d3,d4,d5;

n= sscanf(s.c_str(),"%d%c%d%c%d%c%d%c%d%c%d",

&year,&d1,&month,&d2,&mday,&d3,

&hour,&d4,&min,&d5,&sec);

if(n<5){

throw(invalid_argument("Invalid date format"));

}

fromUTC(year,month,mday,hour,min,sec);

}

/*:13*//*14:*/

#line 466 "timedate.w"

double systemtime::utctoj(long year,int mon,int mday,

int hour,int min,int sec)

{

int a,b,m;

long y;

mon--;

assert(mon>=0&&mon<12);

assert(mday> 0&&mday<32);

assert(hour>=0&&hour<24);

assert(min>=0&&min<60);

assert(sec>=0&&sec<60);

m= mon+1;

y= year;

if(m<=2){

y--;

m+= 12;

}

if((year<1582)||((year==1582)&&((mon<9)||(mon==9&&mday<5)))){

b= 0;

}else{

a= ((int)(y/100));

b= 2-a+(a/4);

}

return(((long)(365.25*(y+4716)))+((int)(30.6001*(m+1)))+

mday+b-1524.5)+

((sec+60L*(min+60L*hour))/86400.0);

}

/*:14*//*15:*/

#line 513 "timedate.w"

void systemtime::jyear(double td,int*yy,int*mm,int*dd)

{

double z,f,a,alpha,b,c,d,e;

td+= 0.5;

z= floor(td);

f= td-z;

if(z<2299161.0){

a= z;

}else{

alpha= floor((z-1867216.25)/36524.25);

a= z+1+alpha-floor(alpha/4);

}

b= a+1524;

c= floor((b-122.1)/365.25);

d= floor(365.25*c);

e= floor((b-d)/30.6001);

store_if(dd,(int)(b-d-floor(30.6001*e)+f));

store_if(mm,(int)((e<14)?(e-1):(e-13)));

store_if(yy,(int)((*mm> 2)?(c-4716):(c-4715)));

}

/*:15*//*16:*/

#line 544 "timedate.w"

void systemtime::jhms(double j,int*h,int*m,int*s)

{

long ij;

j+= 0.5;

ij= (long)((j-floor(j))*86400.0);

store_if(h,(int)(ij/3600L));

store_if(m,(int)((ij/60L)%60L));

store_if(s,(int)(ij%60L));

}

/*:16*//*17:*/

#line 561 "timedate.w"

double systemtime::toJulian(void){

int year,month,mday,hour,min,sec;

toUTC(&year,&month,&mday,&hour,&min,&sec);

return utctoj(year,month,mday,hour,min,sec);

}

void systemtime::fromJulian(double jd){

int year,month,mday,hour,min,sec;

jyear(jd,&year,&month,&mday);

jhms(jd,&hour,&min,&sec);

fromUTC(year,month,mday,hour,min,sec);

}

/*:17*//*18:*/

#line 586 "timedate.w"

double systemtime::gmst(double jd){

double t,theta0;

t= ((floor(jd+0.5)-0.5)-J2000)/JulianCentury;

theta0= d_m_s_to_decimal(6,41,50.54841)+

(d_m_s_to_decimal(0,0,8640184.812866)*t)+

(d_m_s_to_decimal(0,0,0.093104)*t*t)+

(d_m_s_to_decimal(0,0,0.0000062)*t*t*t);

t= (jd+0.5)-(floor(jd+0.5));

theta0+= (t*24)*1.00273790935;

theta0= (theta0-24.0*(floor(theta0/24.0)));

return theta0;

}

/*:18*//*19:*/

#line 633 "timedate.w"

double systemtime::equationTime(double jd)

{

double T,tau,eps,deltaEps,y,L0,e,M,E;

T= (jd-J2000)/JulianCentury;

tau= T/10;

eps= meanObliquityOfEcliptic(jd);

nutation(jd,NULL,&deltaEps);

y= eps+deltaEps;

y= tan(y/2);

y*= y;

L0= 280.4664567+(360007.6982779*tau)+

(0.03032028*tau*tau)+

((tau*tau*tau)/49931)-

((tau*tau*tau*tau)/15300)-

((tau*tau*tau*tau*tau)/2000000);

L0= fixangle(L0);

e= 0.016708634-(0.000042037*T)-

(0.0000001267*T*T);

M= 357.52911+(35999.05029*T)-(0.0001537*T*T);

E= y*sin(dtr(2*L0))-(2*e*sin(dtr(M)))+

(4*e*y*sin(dtr(M))*cos(dtr(2*L0)))-

(((y*y)*sin(dtr(4*L0)))/2)-

((5*(e*e)*sin(dtr(2*M)))/4);

#if 0

cout<<"E = "<<E<<"  degrees = "<<rtd(E)<<

"  min = "<<((rtd(E)/15)*60)<<"\n";

#endif

return rtd(E)*(3600/15);

}

/*:19*//*20:*/

#line 684 "timedate.w"

double systemtime::meanObliquityOfEcliptic(double jd)

{

/*29:*/

#line 1197 "timedate.w"

static double oterms[10]= {

Asec(-4680.93),

Asec(-1.55),

Asec(1999.25),

Asec(-51.38),

Asec(-249.67),

Asec(-39.05),

Asec(7.12),

Asec(27.87),

Asec(5.79),

Asec(2.45)

};

/*:29*/

#line 687 "timedate.w"

;

double eps= 23+(26/60.0)+(21.448/3600.0),u,v;

int i;

v= u= (jd-J2000)/(JulianCentury*100);

if(fabs(u)<1.0){

for(i= 0;i<10;i++){

eps+= oterms[i]*v;

v*= u;

}

}

return dtr(eps);

}

/*:20*//*21:*/

#line 725 "timedate.w"

void systemtime::nutation(double jd,double*deltaPsi,double*deltaEpsilon)

{

int i,j;

double t= (jd-J2000)/JulianCentury,t2,t3,to10;

double ta[5];

double dp= 0,de= 0,ang;

/*30:*/

#line 1222 "timedate.w"

static signed char nutArgMult[NUTERMS][5]= {

{0,0,0,0,1},

{-2,0,0,2,2},

{0,0,0,2,2},

{0,0,0,0,2},

{0,1,0,0,0},

{0,0,1,0,0},

{-2,1,0,2,2},

{0,0,0,2,1},

{0,0,1,2,2},

{-2,-1,0,2,2},

{-2,0,1,0,0},

{-2,0,0,2,1},

{0,0,-1,2,2},

{2,0,0,0,0},

{0,0,1,0,1},

{2,0,-1,2,2},

{0,0,-1,0,1},

{0,0,1,2,1},

{-2,0,2,0,0},

{0,0,-2,2,1},

{2,0,0,2,2},

{0,0,2,2,2},

{0,0,2,0,0},

{-2,0,1,2,2},

{0,0,0,2,0},

{-2,0,0,2,0},

{0,0,-1,2,1},

{0,2,0,0,0},

{2,0,-1,0,1},

{-2,2,0,2,2},

{0,1,0,0,1},

{-2,0,1,0,1},

{0,-1,0,0,1},

{0,0,2,-2,0},

{2,0,-1,2,1},

{2,0,1,2,2},

{0,1,0,2,2},

{-2,1,1,0,0},

{0,-1,0,2,2},

{2,0,0,2,1},

{2,0,1,0,0},

{-2,0,2,2,2},

{-2,0,1,2,1},

{2,0,-2,0,1},

{2,0,0,0,1},

{0,-1,1,0,0},

{-2,-1,0,2,1},

{-2,0,0,0,1},

{0,0,2,2,1},

{-2,0,2,0,1},

{-2,1,0,2,1},

{0,0,1,-2,0},

{-1,0,1,0,0},

{-2,1,0,0,0},

{1,0,0,0,0},

{0,0,1,2,0},

{-1,-1,1,0,0},

{0,1,1,0,0},

{0,-1,1,2,2},

{2,-1,-1,2,2},

{0,0,-2,2,2},

{0,0,3,2,2},

{2,-1,0,2,2}

};

/*:30*/

#line 734 "timedate.w"

;

/*31:*/

#line 1301 "timedate.w"

static long nutArgCoeff[NUTERMS][4]= {

{-171996,-1742,92095,89},

{-13187,-16,5736,-31},

{-2274,-2,977,-5},

{2062,2,-895,5},

{1426,-34,54,-1},

{712,1,-7,0},

{-517,12,224,-6},

{-386,-4,200,0},

{-301,0,129,-1},

{217,-5,-95,3},

{-158,0,0,0},

{129,1,-70,0},

{123,0,-53,0},

{63,0,0,0},

{63,1,-33,0},

{-59,0,26,0},

{-58,-1,32,0},

{-51,0,27,0},

{48,0,0,0},

{46,0,-24,0},

{-38,0,16,0},

{-31,0,13,0},

{29,0,0,0},

{29,0,-12,0},

{26,0,0,0},

{-22,0,0,0},

{21,0,-10,0},

{17,-1,0,0},

{16,0,-8,0},

{-16,1,7,0},

{-15,0,9,0},

{-13,0,7,0},

{-12,0,6,0},

{11,0,0,0},

{-10,0,5,0},

{-8,0,3,0},

{7,0,-3,0},

{-7,0,0,0},

{-7,0,3,0},

{-7,0,3,0},

{6,0,0,0},

{6,0,-3,0},

{6,0,-3,0},

{-6,0,3,0},

{-6,0,3,0},

{5,0,0,0},

{-5,0,3,0},

{-5,0,3,0},

{-5,0,3,0},

{4,0,0,0},

{4,0,0,0},

{4,0,0,0},

{-4,0,0,0},

{-4,0,0,0},

{-4,0,0,0},

{3,0,0,0},

{-3,0,0,0},

{-3,0,0,0},

{-3,0,0,0},

{-3,0,0,0},

{-3,0,0,0},

{-3,0,0,0},

{-3,0,0,0}

};

/*:31*/

#line 735 "timedate.w"

;

t3= t*(t2= t*t);

ta[0]= dtr(297.85036+445267.11148*t-0.0019142*t2+

t3/189474.0);

ta[1]= dtr(357.52772+35999.05034*t-0.0001603*t2-

t3/300000.0);

ta[2]= dtr(134.96298+477198.867398*t+0.0086972*t2+

t3/56250.0);

ta[3]= dtr(93.27191+483202.017538*t-0.0036825*t2+

t3/327270.0);

ta[4]= dtr(125.04452-1934.136261*t+0.0020708*t2+

t3/450000.0);

for(i= 0;i<5;i++){

ta[i]= fixangr(ta[i]);

}

to10= t/10.0;

for(i= 0;i<NUTERMS;i++){

ang= 0;

for(j= 0;j<5;j++){

if(nutArgMult[i][j]!=0){

ang+= nutArgMult[i][j]*ta[j];

}

}

dp+= (nutArgCoeff[i][0]+nutArgCoeff[i][1]*to10)*sin(ang);

de+= (nutArgCoeff[i][2]+nutArgCoeff[i][3]*to10)*cos(ang);

}

dstore_if(deltaPsi,dtr(dp/(3600.0*10000.0)));

dstore_if(deltaEpsilon,dtr(de/(3600.0*10000.0)));

}

/*:21*//*22:*/

#line 788 "timedate.w"

void systemtime::sunpos(double jd,bool apparent,

double*ra,double*dec,double*rv,double*slong)

{

double t,t2,C,l,m,e,v,theta,omega,eps;

t= (jd-J2000)/JulianCentury;

t2= t*t;

l= fixangle(280.46646+(36000.76983*t)+(0.0003032*t2));

m= fixangle(357.52911+(35999.05029*t)-(0.0001537*t2));

e= 0.016708634-(0.000042037*t)-(0.0000001267*t2);

C= ((1.914602-(0.004817*t)-(0.000014*t2))*sin(dtr(m)))+

((0.019993-(0.000101*t))*sin(2*dtr(m)))+

(0.000289*sin(3*dtr(m)));

theta= l+C;

v= m+C;

eps= rtd(meanObliquityOfEcliptic(jd));

if(apparent){

omega= fixangle(125.04-1934.136*t);

theta= theta-0.00569-0.00478*sin(dtr(omega));

eps+= 0.00256*cos(dtr(omega));

}

dstore_if(slong,dtr(theta));

dstore_if(rv,(1.000001018*(1-(e*e)))/(1+e*cos(dtr(v))));

dstore_if(ra,fixangr(atan2(cos(dtr(eps))*sin(dtr(theta)),cos(dtr(theta)))));

dstore_if(dec,asin(sin(dtr(eps))*sin(dtr(theta))));

}

/*:22*//*23:*/

#line 868 "timedate.w"

void systemtime::moon_and_sun(double jd,double*phaseang,

double*pphase,double*mage,double*dist,

double*angdia,double*sudist,

double*suangdia)

{

double Day,N,M,Ec,Lambdasun,ml,MM,MN,Ev,Ae,A3,MmP,

mEc,A4,lP,Varia,lPP,NP,y,x,Lambdamoon,

MoonAge,MoonPhase,

MoonDist,MoonDFrac,MoonAng,

F,SunDist,SunAng;

static const double

AstronomicalUnit= 149597870.0,

SunSMAX= (AstronomicalUnit*1.000001018),

epoch= 2444238.5,

elonge= 278.833540,

elongp= 282.596403,

eccent= 0.016718,

sunangsiz= 0.533128,

mmlong= 64.975464,

mmlongp= 349.383063,

mlnode= 151.950429,

minc= 5.145396,

mecc= 0.054900,

mangsiz= 0.5181,

msmax= 384401.0;

Day= jd-epoch;

N= fixangle((360/365.2422)*Day);

M= fixangle(N+elonge-elongp);

Ec= kepler(M,eccent);

Ec= sqrt((1+eccent)/(1-eccent))*tan(Ec/2);

Ec= 2*rtd(atan(Ec));

Lambdasun= fixangle(Ec+elongp);

F= ((1+eccent*cos(dtr(Ec)))/(1-eccent*eccent));

SunDist= SunSMAX/F;

SunAng= F*sunangsiz;

ml= fixangle(13.1763966*Day+mmlong);

MM= fixangle(ml-0.1114041*Day-mmlongp);

MN= fixangle(mlnode-0.0529539*Day);

Ev= 1.2739*sin(dtr(2*(ml-Lambdasun)-MM));

Ae= 0.1858*sin(dtr(M));

A3= 0.37*sin(dtr(M));

MmP= MM+Ev-Ae-A3;

mEc= 6.2886*sin(dtr(MmP));

A4= 0.214*sin(dtr(2*MmP));

lP= ml+Ev+mEc-Ae+A4;

Varia= 0.6583*sin(dtr(2*(lP-Lambdasun)));

lPP= lP+Varia;

NP= MN-0.16*sin(dtr(M));

y= sin(dtr(lPP-NP))*cos(dtr(minc));

x= cos(dtr(lPP-NP));

Lambdamoon= rtd(atan2(y,x));

Lambdamoon+= NP;

MoonAge= lPP-Lambdasun;

MoonPhase= (1-cos(dtr(MoonAge)))/2;

MoonDist= (msmax*(1-mecc*mecc))/

(1+mecc*cos(dtr(MmP+mEc)));

MoonDFrac= MoonDist/msmax;

MoonAng= mangsiz/MoonDFrac;

dstore_if(pphase,MoonPhase);

dstore_if(mage,SynMonth*(fixangle(MoonAge)/360.0));

dstore_if(dist,MoonDist);

dstore_if(angdia,dtr(MoonAng));

dstore_if(sudist,SunDist);

dstore_if(suangdia,dtr(SunAng));

dstore_if(phaseang,fixangle(MoonAge)/360.0);

}

/*:23*//*24:*/

#line 1024 "timedate.w"

double systemtime::kepler(double m,double ecc)

{

double E,delta;

static const double EPSILON= 1E-6;

E= m= dtr(m);

do{

delta= E-ecc*sin(E)-m;

E-= delta/(1-ecc*cos(E));

}while(fabs(delta)> EPSILON);

return E;

}

/*:24*//*25:*/

#line 1045 "timedate.w"

double systemtime::meanphase(double sdate,double k)

{

double t,t2,t3,nt1;

t= (sdate-2415020.0)/JulianCentury;

t2= t*t;

t3= t2*t;

nt1= 2415020.75933+SynMonth*k

+0.0001178*t2

-0.000000155*t3

+0.00033*sin(dtr(166.56+132.87*t-0.009173*t2));

return nt1;

}

/*:25*//*26:*/

#line 1068 "timedate.w"

double systemtime::truephase(double k,double phase)

{

double t,t2,t3,pt,m,mprime,f;

k+= phase;

t= k/1236.85;

t2= t*t;

t3= t2*t;

pt= 2415020.75933

+SynMonth*k

+0.0001178*t2

-0.000000155*t3

+0.00033*sin(dtr(166.56+132.87*t-0.009173*t2));

m= 359.2242

+29.10535608*k

-0.0000333*t2

-0.00000347*t3;

mprime= 306.0253

+385.81691806*k

+0.0107306*t2

+0.00001236*t3;

f= 21.2964

+390.67050646*k

-0.0016528*t2

-0.00000239*t3;

if((phase<0.01)||(fabs(phase-0.5)<0.01)){

pt+= (0.1734-0.000393*t)*sin(dtr(m))

+0.0021*sin(dtr(2*m))

-0.4068*sin(dtr(mprime))

+0.0161*sin(dtr(2*mprime))

-0.0004*sin(dtr(3*mprime))

+0.0104*sin(dtr(2*f))

-0.0051*sin(dtr(m+mprime))

-0.0074*sin(dtr(m-mprime))

+0.0004*sin(dtr(2*f+m))

-0.0004*sin(dtr(2*f-m))

-0.0006*sin(dtr(2*f+mprime))

+0.0010*sin(dtr(2*f-mprime))

+0.0005*sin(dtr(m+2*mprime));

}else if((fabs(phase-0.25)<0.01||(fabs(phase-0.75)<0.01))){

pt+= (0.1721-0.0004*t)*sin(dtr(m))

+0.0021*sin(dtr(2*m))

-0.6280*sin(dtr(mprime))

+0.0089*sin(dtr(2*mprime))

-0.0004*sin(dtr(3*mprime))

+0.0079*sin(dtr(2*f))

-0.0119*sin(dtr(m+mprime))

-0.0047*sin(dtr(m-mprime))

+0.0003*sin(dtr(2*f+m))

-0.0004*sin(dtr(2*f-m))

-0.0006*sin(dtr(2*f+mprime))

+0.0021*sin(dtr(2*f-mprime))

+0.0003*sin(dtr(m+2*mprime))

+0.0004*sin(dtr(m-2*mprime))

-0.0003*sin(dtr(2*m+mprime));

if(phase<0.5)

pt+= 0.0028-0.0004*cos(dtr(m))+0.0003*cos(dtr(mprime));

else

pt+= -0.0028+0.0004*cos(dtr(m))-0.0003*cos(dtr(mprime));

}

return pt;

}

/*:26*//*27:*/

#line 1155 "timedate.w"

void systemtime::moonphases(double sdate,double phases[5])

{

double adate,k1,k2;

int yy,mm,dd;

adate= sdate-45;

jyear(adate,&yy,&mm,&dd);

k1= floor((yy+((mm-1)*(1.0/12.0))-1900)*12.3685);

adate= meanphase(adate,k1);

while(true){

adate+= SynMonth;

k2= k1+1;

if(truephase(k1,0.0)<=sdate&&truephase(k2,0.0)> sdate){

break;

}

k1= k2;

}

phases[0]= truephase(k1,0.0);

phases[1]= truephase(k1,0.25);

phases[2]= truephase(k1,0.5);

phases[3]= truephase(k1,0.75);

phases[4]= truephase(k2,0.0);

}

/*:27*/

#line 108 "timedate.w"

/*:5*/

#line 73 "timedate.w"

/*:2*/

///////////////////////////////////////////////////////////////////////

/*1:*/

#line 58 "timedate.w"

#define REVDATE "2nd February 2002"

/*:1*//*32:*/

#line 1377 "timedate.w"

/*35:*/

#line 1541 "timedate.w"

#include "config.h"         

#include <iostream> 

#include <exception> 

#include <stdexcept> 

#include <string> 

using namespace std;

#include <stdio.h> 

#include <stdlib.h> 

#include <ctype.h> 

#ifdef HAVE_GETOPT

#ifdef HAVE_UNISTD_H

#include <unistd.h> 

#endif

#else

#include "getopt.h"     

#endif

#include "timedate.h"        

/*:35*/

#line 1379 "timedate.w"

;

/*34:*/

#line 1523 "timedate.w"

static void usage(void)

{

cout<<PRODUCT<<"  --  Analyse eggsummary files.  Call:\n";

cout<<"              "<<PRODUCT<<" [options] [infile] [outfile]\n";

cout<<"\n";

cout<<"Options:\n";

cout<<"           --copyright       Print copyright information\n";

cout<<"           -u, --help        Print this message\n";

cout<<"           --version         Print version number\n";

cout<<"\n";

cout<<"by John Walker\n";

cout<<"http://www.fourmilab.ch/\n";

}

/*:34*/

#line 1381 "timedate.w"

;

int main(int argc,char*argv[])

{

extern char*optarg;

extern int optind;

int opt;

/*33:*/

#line 1491 "timedate.w"

while((opt= getopt(argc,argv,"nu-:"))!=-1){

switch(opt){

case'u':

case'?':

usage();

return 0;

case'-':

switch(optarg[0]){

case'c':

cout<<"This program is in the public domain.\n";

return 0;

case'h':

usage();

return 0;

case'v':

cout<<PRODUCT<<" "<<VERSION<<"\n";

cout<<"Last revised: "<<REVDATE<<"\n";

cout<<"The latest version is always available\n";

cout<<"at http://www.fourmilab.ch/eggtools/eggshell\n";

return 0;

}

}

}

/*:33*/

#line 1390 "timedate.w"

;

string s;

systemtime td,td1;

td.now();

cout<<td.dateToString()<<" "<<td.timeToString()<<"\n";

td1= td;

td1.nextDay();

cout<<"Tomorrow: "<<td1.dateToString()<<" "<<td1.timeToString()<<"\n";

td1= td;

td1.previousDay();

cout<<"Yesterday: "<<td1.dateToString()<<" "<<td1.timeToString()<<"\n";

td1= td;

td1.dayStep(7);

cout<<"Next week: "<<td1.dateToString()<<" "<<td1.timeToString()<<"\n";

td.fromString("1981-04-12 13:51:18");

cout<<td.dateToString()<<" "<<td.timeToString()<<"\n";

td.fromString("2021-04-12");

cout<<td.dateToString()<<" "<<td.timeToString()<<"\n";

cout<<"Weekday = "<<systemtime::jwday(systemtime::utctoj(2001,

5,12,23,59,59))<<"\n";

td.fromJulian(2436116.31);

cout<<td.dateToString()<<" "<<td.timeToString()<<"\n";

cout<<"To Julian = "<<td.toJulian()<<" Weekday = "<<td.weekday()<<"\n";

td.fromString("1987-04-10 19:21:00");

cout<<"Sidereal time = "<<angle::degrees_to_d_m_s(td.siderealTime())<<"\n";

cout.flush();

int t1,t2,t3;

systemtime::jyear(2436116.31,&t1,&t2,&t3);

printf("%04d-%02d-%02d\n",t1,t2,t3);

systemtime::jhms(systemtime::utctoj(1989,11,14,16,23,18),&t1,&t2,&t3);

printf("%02d:%02d:%02d\n",t1,t2,t3);

{

double j= 2446895.5;

double o= angle::rtd(systemtime::meanObliquityOfEcliptic(j));

double delpsi,deleps;

cout<<"meanObliquityOfEcliptic("<<((long)j)<<") = "<<

angle::degrees_to_d_m_s(o)<<"\n";

systemtime::nutation(j,&delpsi,&deleps);

cout<<"  delpsi = "<<

angle::degrees_to_d_m_s(angle::rtd(delpsi))<<

" deleps = "<<

angle::degrees_to_d_m_s(angle::rtd(deleps))<<"\n";

double eOt= systemtime::equationTime(2448908.5);

cout<<"Equation of time = "<<((int)(eOt/60))<<" minutes "<<

(eOt-(60*((int)(eOt/60))))<<" seconds\n";

j= 2448908.5;

double sura,sudec,surv,suslong;

systemtime::sunpos(j,true,

&sura,&sudec,&surv,&suslong);

cout<<"Sun position:  RA = "<<angle::rtd(sura)<<

"  Dec = "<<angle::rtd(sudec)<<"  Rv = "<<surv<<

"  Long = "<<angle::rtd(suslong)<<"\n";

}

{

double phaseang,pphase,mage,dist,angdia,

sudist,suangdia;

systemtime::moon_and_sun(2452060.13874,

&phaseang,&pphase,&mage,&dist,&angdia,

&sudist,&suangdia);

cout<<"Moon:  Phaseang = "<<(phaseang*100)<<

"  Pphase = "<<((int)(pphase*100))<<"%\n       "<<

"  Moon age = "<<((int)mage)<<" days, "<<

angle::degrees_to_d_m_s((mage-((int)mage))*24)<<

"  Dist = "<<dist<<" km\n"<<

"Sun:   Angdia = "<<angle::rtd(angdia)<<

"  Sudist = "<<((int)sudist)<<" km"<<

"  Suangdia = "<<angle::rtd(suangdia)<<"\n";

}

{

double phases[5];

systemtime pt[5];

int i;

systemtime::moonphases(2452060.13874,phases);

for(i= 0;i<5;i++){

pt[i].fromJulian(phases[i]);

cout<<"Phase "<<i<<": "<<pt[i].dateToString()<<

" "<<pt[i].timeToString()<<"\n";

}

}

return 0;

}

/*:32*/