source: LMDZ5/branches/IPSLCM6.0.11/libf/phymar/fctast.h @ 3110

Last change on this file since 3110 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR


Integration of MAR physics

File size: 1.1 KB
Line 
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
11REAL_B :: RTETA,REL,REM,RRS,RLLS,RLLLS,RDS,RET
12REAL_B :: PTIME,PTETA
13
14RTETA(PTIME)=PTIME/(RDAY*365.25_JPRB)
15REL(PTETA)=1.7535_JPRB+6.283076_JPRB*PTETA
16REM(PTETA)=6.240075_JPRB+6.283020_JPRB*PTETA
17RRS(PTETA)=REA*(1.0001_JPRB-0.0163_JPRB*SIN(REL(PTETA))&
18          &+0.0037_JPRB*COS(REL(PTETA)))
19! Relative movement Sun/Earth
20RLLS(PTETA)=4.8951_JPRB+6.283076_JPRB*PTETA
21RLLLS(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(_TWO_*REL(PTETA))&
23         &+0.0002_JPRB*COS(_TWO_*REL(PTETA))
24RDS(PTETA)=ASIN(SIN(REPSM)*SIN(RLLLS(PTETA)))
25RET(PTETA)=591.8_JPRB*SIN(_TWO_*RLLS(PTETA))-459.4_JPRB*SIN(REM(PTETA))&
26  &+39.5_JPRB*SIN(REM(PTETA))*COS(_TWO_*RLLS(PTETA))&
27  &-12.7_JPRB*SIN(4._JPRB*RLLS(PTETA))-4.8_JPRB*SIN(_TWO_*REM(PTETA))
28!    -------------------------------------------------------------
29
Note: See TracBrowser for help on using the repository browser.