source: LMDZ5/branches/testing/libf/phylmd/rrtm/fctast.h @ 5407

Last change on this file since 5407 was 1999, checked in by Laurent Fairhead, 11 years ago

Merged trunk changes r1920:1997 into testing branch
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 1.2 KB
RevLine 
[1989]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(KIND=JPRB) :: RTETA,REL,REM,RRS,RLLS,RLLLS,RDS,RET
12REAL(KIND=JPRB) :: 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(2.0_JPRB*REL(PTETA))&
23         &+0.0002_JPRB*COS(2.0_JPRB*REL(PTETA))
24RDS(PTETA)=ASIN(SIN(REPSM)*SIN(RLLLS(PTETA)))
25RET(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
Note: See TracBrowser for help on using the repository browser.