空间大地测量与gps导航定位时间系统相互转换,格里高利时通用时儒略日,gps时,年积日相互转换的源代码程序

cond);Р pdoy->tod.tos=pct->second-int(pct->second); /*pct->hour*3600Р +pct->minute*60+pct->second-pdoy->tod.sn;*/Р}Рvoid monTime (PDOY pdoy, PCOMMONTIME pct)Р{Р?PCOMMONTIME pcto=MONTIME;Р?pcto->year=pdoy->year;Р?pcto->month=1;Р?pcto->day=1;Р?pcto->hour=0;Р?pcto->minute=0;Р?pcto->second=0;Р?PJULIANDAY pjdo=new JULIANDAY;Р double JD,JDO;Р?CommonTimeToJulianDay(pcto,pjdo);Р?JDO=pjdo->day+(pjdo->tod.sn+pjdo->tod.tos)/86400;Р?JD=JDO+pdoy->day+(pdoy->tod.sn+pdoy->tod.tos)/86400-1;Р?long a,b,c,d,e;Р?a=(long)(JD+0.5);Р?b=a+1537;Р?c=(long)((b-122.1)/365.25);Р?d=(long)(365.25*c);Р?e=(long)((b-d)/30.6001);Р?pct->day=short(b-d-(long)(30.6001*e)+FRAC(JD+0.5));Р?pct->month=short(e-1-12*(long)(e/14));Р?pct->year=short(c-4715-(long)((7+pct->month)/10));Р?pct->hour=short((pdoy->tod.sn+pdoy->tod.tos)/3600);