trunk/phylmd/test_orbite.f90

program test_orbite

  use orbite_m, only: orbite

  implicit none

  REAL xjour ! jour de l'ann\'ee \`a compter du premier janvier

  REAL longi
  ! longitude vraie de la Terre dans son orbite solaire, par rapport
  ! au point vernal (21 mars), en degr\'es

  REAL dist ! distance terre-soleil, en ua
  integer i
  integer, parameter:: n = 100
  real delta_x

  !------------------------------------------------------------------------

  xjour = 1.
  delta_x = 359. / (n - 1)
  print *, '"" degrees ua'
  print *, "xjour longi dist"

  do i = 1, n
     call orbite(xjour, longi, dist)
     print *, xjour, longi, dist
     xjour = xjour + delta_x
  end do

end program test_orbite
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