|
Last change
on this file since 3325 was
3040,
checked in by jbclement, 15 months ago
|
|
Mars PCM:
The variable 'timeperi' (defined in "planete_h.F90" and computed in "iniorbit.F") is renamed into 'lsperi' and thus slightly changed to be coherent to the solar longitude of perihelion in radian. It can now be used out of the box by other subroutines/programs like the PEM.
JBC
|
|
File size:
937 bytes
|
| Line | |
|---|
| 1 | SUBROUTINE orbite(pls,pdist_sol,pdecli) |
|---|
| 2 | USE planete_h, ONLY: e_elips, p_elips, obliquit, lsperi |
|---|
| 3 | USE comcstfi_h, ONLY: pi |
|---|
| 4 | IMPLICIT NONE |
|---|
| 5 | |
|---|
| 6 | c======================================================================= |
|---|
| 7 | c |
|---|
| 8 | c Compute distance to Sun and declination as a function of the solar |
|---|
| 9 | c longitude Ls |
|---|
| 10 | c |
|---|
| 11 | c Arguments: |
|---|
| 12 | c ---------- |
|---|
| 13 | c |
|---|
| 14 | c Input: |
|---|
| 15 | c ------ |
|---|
| 16 | c pls Ls (radians) |
|---|
| 17 | c |
|---|
| 18 | c Output: |
|---|
| 19 | c ------- |
|---|
| 20 | c pdist_sol Distance Sun-Planet in UA |
|---|
| 21 | c pdecli Declinaison ( in radians ) |
|---|
| 22 | c |
|---|
| 23 | c======================================================================= |
|---|
| 24 | |
|---|
| 25 | c arguments: |
|---|
| 26 | c ---------- |
|---|
| 27 | |
|---|
| 28 | REAL,INTENT(IN) :: pls |
|---|
| 29 | REAL,INTENT(OUT) :: pdist_sol, pdecli |
|---|
| 30 | |
|---|
| 31 | c----------------------------------------------------------------------- |
|---|
| 32 | |
|---|
| 33 | c Distance Sun-Planet |
|---|
| 34 | pdist_sol = p_elips/(1.+e_elips*cos(pls-lsperi)) |
|---|
| 35 | |
|---|
| 36 | c Solar declination |
|---|
| 37 | pdecli = asin(sin(pls)*sin(obliquit*pi/180.)) |
|---|
| 38 | |
|---|
| 39 | END |
|---|
Note: See
TracBrowser
for help on using the repository browser.
