source: trunk/LMDZ.PLUTO.old/libf/phypluto/call_dayperi.F @ 3670

Last change on this file since 3670 was 3175, checked in by emillour, 2 years ago

Pluto PCM:
Add the old Pluto LMDZ for reference (required prior step to making
an LMDZ.PLUTO using the same framework as the other physics packages).
TB+EM
File size: 1.2 KB
Line 
1c ********************************************************
2      subroutine call_dayperi(Lsperi,e_elips,dayperi,year_day)
3      IMPLICIT NONE
4
5c=======================================================================
6c
7c   Objet:
8c   ------
9c   Computing the Plutonian date (number of sols since Ls=0 = spring
10c   equinox) at perihelion. Called by testphys1D only
11c
12c   Arguments:
13c   ----------
14c
15c   Input:
16c   ------
17c   Lsperi      solar longitude (Ls) of perohelion (rad)
18c   e_elips       Excentricity
19c
20c   output
21c   ------
22c   dayperi       Plutonian date at perihelion (sol)
23c   real year_day      ! number of sols per Pluto year
24c-----------------------------------------------------------------------
25
26c arguments:
27c ----------
28       real Lsperi,e_elips,dayperi
29       real year_day      ! number of sols per Pluto year ~ 14164.46
30
31
32c      Local
33       real x1, x2,pi
34
35
36       pi=2.*asin(1.)
37     
38       x1 = sqrt((1.-e_elips)/(1.+e_elips))
39       x2 = e_elips*sqrt(1.-e_elips**2.)
40
41       dayperi = 0.5*(year_day/pi)*
42     &   ( 2.*atan(x1*tan(0.5*Lsperi))
43     &     -x2*sin(Lsperi)/(1.+e_elips*cos(Lsperi)) )
44       if(dayperi.lt.0) dayperi=dayperi+year_day
45       return
46       end
Note: See TracBrowser for help on using the repository browser.