1 |
program test_orbite |
2 |
|
3 |
use orbite_m, only: orbite |
4 |
|
5 |
implicit none |
6 |
|
7 |
REAL xjour ! jour de l'ann\'ee \`a compter du premier janvier |
8 |
|
9 |
REAL longi |
10 |
! longitude vraie de la Terre dans son orbite solaire, par rapport |
11 |
! au point vernal (21 mars), en degr\'es |
12 |
|
13 |
REAL dist ! distance terre-soleil, en ua |
14 |
integer i |
15 |
integer, parameter:: n = 100 |
16 |
real delta_x |
17 |
|
18 |
!------------------------------------------------------------------------ |
19 |
|
20 |
xjour = 1. |
21 |
delta_x = 359. / (n - 1) |
22 |
print *, '"" degrees ua' |
23 |
print *, "xjour longi dist" |
24 |
|
25 |
do i = 1, n |
26 |
call orbite(xjour, longi, dist) |
27 |
print *, xjour, longi, dist |
28 |
xjour = xjour + delta_x |
29 |
end do |
30 |
|
31 |
end program test_orbite |