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
RevLine 
[253]1      subroutine orbite(pls,pdist_star,pdecli)
[1308]2     
3      use planete_mod, only: p_elips, e_elips, timeperi, obliquit
[253]4      implicit none
[135]5
[253]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
[135]24c   Declarations:
25c   -------------
26
[1308]27!#include "planete.h"
[135]28#include "comcstfi.h"
29
30c arguments:
31c ----------
32
[253]33      REAL pday,pdist_star,pdecli,pls,i
[135]34
35c-----------------------------------------------------------------------
36
[253]37c Star-Planet Distance
[135]38
[253]39      pdist_star = p_elips/(1.+e_elips*cos(pls+timeperi))
[135]40
[253]41c Stellar declination
[135]42
43c ********************* version before 01/01/2000 *******
44
[253]45      pdecli = asin (sin(pls)*sin(obliquit*pi/180.))
[135]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.