|
Last change
on this file since 3572 was
3175,
checked in by emillour, 18 months 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
|
| Rev | Line | |
|---|
| [3175] | 1 | c ******************************************************** |
|---|
| 2 | subroutine call_dayperi(Lsperi,e_elips,dayperi,year_day) |
|---|
| 3 | IMPLICIT NONE |
|---|
| 4 | |
|---|
| 5 | c======================================================================= |
|---|
| 6 | c |
|---|
| 7 | c Objet: |
|---|
| 8 | c ------ |
|---|
| 9 | c Computing the Plutonian date (number of sols since Ls=0 = spring |
|---|
| 10 | c equinox) at perihelion. Called by testphys1D only |
|---|
| 11 | c |
|---|
| 12 | c Arguments: |
|---|
| 13 | c ---------- |
|---|
| 14 | c |
|---|
| 15 | c Input: |
|---|
| 16 | c ------ |
|---|
| 17 | c Lsperi solar longitude (Ls) of perohelion (rad) |
|---|
| 18 | c e_elips Excentricity |
|---|
| 19 | c |
|---|
| 20 | c output |
|---|
| 21 | c ------ |
|---|
| 22 | c dayperi Plutonian date at perihelion (sol) |
|---|
| 23 | c real year_day ! number of sols per Pluto year |
|---|
| 24 | c----------------------------------------------------------------------- |
|---|
| 25 | |
|---|
| 26 | c arguments: |
|---|
| 27 | c ---------- |
|---|
| 28 | real Lsperi,e_elips,dayperi |
|---|
| 29 | real year_day ! number of sols per Pluto year ~ 14164.46 |
|---|
| 30 | |
|---|
| 31 | |
|---|
| 32 | c 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.
