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	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