1 |
guez |
319 |
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 |