source: trunk/LMDZ.GENERIC/libf/phystd/orbite.F @ 1316

Last change on this file since 1316 was 1308, checked in by emillour, 11 years ago

Generic GCM:
Some cleanup to simplify dynamics/physics interactions by getting rid
of dimphys.h (i.e. the nlayermx parameter) and minimizing use of
dimension.h in the physics.
EM
File size: 1.2 KB
Line 
1      subroutine orbite(pls,pdist_star,pdecli)
2     
3      use planete_mod, only: p_elips, e_elips, timeperi, obliquit
4      implicit none
5
6!==================================================================
7!     
8!     Purpose
9!     -------
10!     Distance from star and declination as a function of the stellar
11!     longitude Ls
12!     
13!     Inputs
14!     ------
15!     pls          Ls
16!
17!     Outputs
18!     -------
19!     pdist_star    Distance Star-Planet in UA
20!     pdecli        declinaison ( in radians )
21!
22!=======================================================================
23
24c   Declarations:
25c   -------------
26
27!#include "planete.h"
28#include "comcstfi.h"
29
30c arguments:
31c ----------
32
33      REAL pday,pdist_star,pdecli,pls,i
34
35c-----------------------------------------------------------------------
36
37c Star-Planet Distance
38
39      pdist_star = p_elips/(1.+e_elips*cos(pls+timeperi))
40
41c Stellar declination
42
43c ********************* version before 01/01/2000 *******
44
45      pdecli = asin (sin(pls)*sin(obliquit*pi/180.))
46
47c ********************* version after 01/01/2000 *******
48c     i=obliquit*pi/180.
49c     pdecli=asin(sin(pls)*sin(i)/sqrt(sin(pls)**2+
50c    & cos(pls)**2*cos(i)**2))
51c ******************************************************
52
53      RETURN
54      END
Note: See TracBrowser for help on using the repository browser.