| 1 | ! ----------------------------------------------------------- |
|---|
| 2 | |
|---|
| 3 | ! - Astronomical functions |
|---|
| 4 | ! you will find the description in the annex 1 of the documentation |
|---|
| 5 | ! RRS is the distance Sun-Earth |
|---|
| 6 | ! RDS is the declination of the Earth |
|---|
| 7 | ! RET is the equation of time |
|---|
| 8 | |
|---|
| 9 | ! Orbit of the earth |
|---|
| 10 | |
|---|
| 11 | REAL(KIND=JPRB) :: RTETA,REL,REM,RRS,RLLS,RLLLS,RDS,RET |
|---|
| 12 | REAL(KIND=JPRB) :: PTIME,PTETA |
|---|
| 13 | |
|---|
| 14 | RTETA(PTIME)=PTIME/(RDAY*365.25_JPRB) |
|---|
| 15 | REL(PTETA)=1.7535_JPRB+6.283076_JPRB*PTETA |
|---|
| 16 | REM(PTETA)=6.240075_JPRB+6.283020_JPRB*PTETA |
|---|
| 17 | RRS(PTETA)=REA*(1.0001_JPRB-0.0163_JPRB*SIN(REL(PTETA))& |
|---|
| 18 | &+0.0037_JPRB*COS(REL(PTETA))) |
|---|
| 19 | ! Relative movement Sun/Earth |
|---|
| 20 | RLLS(PTETA)=4.8951_JPRB+6.283076_JPRB*PTETA |
|---|
| 21 | RLLLS(PTETA)=4.8952_JPRB+6.283320_JPRB*PTETA-0.0075_JPRB*SIN(REL(PTETA))& |
|---|
| 22 | &-0.0326_JPRB*COS(REL(PTETA))-0.0003_JPRB*SIN(2.0_JPRB*REL(PTETA))& |
|---|
| 23 | &+0.0002_JPRB*COS(2.0_JPRB*REL(PTETA)) |
|---|
| 24 | RDS(PTETA)=ASIN(SIN(REPSM)*SIN(RLLLS(PTETA))) |
|---|
| 25 | RET(PTETA)=591.8_JPRB*SIN(2.0_JPRB*RLLS(PTETA))-459.4_JPRB*SIN(REM(PTETA))& |
|---|
| 26 | &+39.5_JPRB*SIN(REM(PTETA))*COS(2.0_JPRB*RLLS(PTETA))& |
|---|
| 27 | &-12.7_JPRB*SIN(4._JPRB*RLLS(PTETA))-4.8_JPRB*SIN(2.0_JPRB*REM(PTETA)) |
|---|
| 28 | ! ------------------------------------------------------------- |
|---|
| 29 | |
|---|
