I've been involved with this project for 25 years. I wrote a couple of applications in the early 1980s that have changed little over that time
Also, for the code below I have hardcopy from 1978 and it's identical line-for-line, including the error 2nd line.
FUNCTION RACRAE(IYRDY,IHMS,RAC) C RACRAE PHILLI 0174 NAV: CONVERTS CELESTIAL ONG TO EARTH LON C $ FUNCTION RACRAE(IYRDY, IHMS, RAC) (DAS) C $ CONVERT CONVERT CELESTIAL LONGITUDE TO EARTH LONGITUDE. FN VAL IS C $ IN REAL*4 DEGREES. C $ IYRDY = (I) INPUT YEAR AND JULIAN DAY (YYDDD) C $ IHMS = (I) INPUT TIME (HHMMSS) C $ RAE = (R) INPUT CELESTIAL LONGITUDE (DEGREES) C $$ RACRAE = NAVIGATION, CONVERT, LONGITUDE C C
DOUBLE PRECISION TIMDIF,RAHA,SOLSID,SHA
SHA=100.26467D0
IRAYD=74001
IRAHMS=0
SOLSID=1.00273791D0
RAHA=RAC-SHA+TIMDIF(IYRDY,IHMS,IRAYD,IRAHMS)*SOLSI D/4.0D0
RAE=DMOD(RAHA,360.0D0)
IF(RAE.LT.0.0)RAE=RAE+360.0
RACRAE=RAE
RETURN
END
I've been involved with this project for 25 years. I wrote a couple of applications in the early 1980s that have changed little over that time
I D/4.0D0
Also, for the code below I have hardcopy from 1978 and it's identical line-for-line, including the error 2nd line.
FUNCTION RACRAE(IYRDY,IHMS,RAC)
C RACRAE PHILLI 0174 NAV: CONVERTS CELESTIAL ONG TO EARTH LON
C $ FUNCTION RACRAE(IYRDY, IHMS, RAC) (DAS)
C $ CONVERT CONVERT CELESTIAL LONGITUDE TO EARTH LONGITUDE. FN VAL IS
C $ IN REAL*4 DEGREES.
C $ IYRDY = (I) INPUT YEAR AND JULIAN DAY (YYDDD)
C $ IHMS = (I) INPUT TIME (HHMMSS)
C $ RAE = (R) INPUT CELESTIAL LONGITUDE (DEGREES)
C $$ RACRAE = NAVIGATION, CONVERT, LONGITUDE
C
C
DOUBLE PRECISION TIMDIF,RAHA,SOLSID,SHA
SHA=100.26467D0
IRAYD=74001
IRAHMS=0
SOLSID=1.00273791D0
RAHA=RAC-SHA+TIMDIF(IYRDY,IHMS,IRAYD,IRAHMS)*SOLS
RAE=DMOD(RAHA,360.0D0)
IF(RAE.LT.0.0)RAE=RAE+360.0
RACRAE=RAE
RETURN
END