shijianzhuanhuan.f90

计算卫星位置的fortran程序

SUBROUTINE G_GPS(Y,M,D,H,MIN,S,WN,TOW)    !格里高利历转化为儒略日
IMPLICIT NONE
REAL*8 TOW,S
INTEGER Y,M,D,MIN,H,SUM,WN
DOUBLE PRECISION JD,H1
    Y=Y+2000
!    WRITE(10,*)Y,"年",M,"月",D,"日",H,"时",MIN,"分",S,"秒"
    IF(M<=2)THEN
       Y=Y-1
       M=M+12
	ENDIF
	   JD=INT(365.25*Y)+INT(30.6001*(M+1))+D+1720981.5
!	   WRITE(*,*)"JD=",JD
	   WN=INT((JD-2444244.5)/7)
       TOW=((JD-2444244.5)-WN*7)*86400.0+H*3600+MIN*60+S !+WN*604800
!	   WRITE(*,*)"WN=",WN,"TOW=",TOW
	RETURN
END