1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
!/*-----------------------------------------------------------------*/
!/*!
! \file f77getpositionrecordindex.f
! \brief Check if calceph_getpositionrecord works with fortran 77 compiler.
!
! \author M. Gastineau
! Astronomie et Systemes Dynamiques, IMCCE, CNRS, Observatoire de Paris.
!
! Copyright, 2008-2021, CNRS
! email of the author : Mickael.Gastineau@obspm.fr
!
!*/
!/*-----------------------------------------------------------------*/
!/*-----------------------------------------------------------------*/
!/* License of this file :
! This file is "triple-licensed", you have to choose one of the three licenses
! below to apply on this file.
!
! CeCILL-C
! The CeCILL-C license is close to the GNU LGPL.
! ( http://www.cecill.info/licences/Licence_CeCILL-C_V1-en.html )
!
! or CeCILL-B
! The CeCILL-B license is close to the BSD.
! ( http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.txt)
!
! or CeCILL v2.1
! The CeCILL license is compatible with the GNU GPL.
! ( http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.html )
!
!
! This library is governed by the CeCILL-C, CeCILL-B or the CeCILL license under
! French law and abiding by the rules of distribution of free software.
! You can use, modify and/ or redistribute the software under the terms
! of the CeCILL-C,CeCILL-B or CeCILL license as circulated by CEA, CNRS and INRIA
! at the following URL "http://www.cecill.info".
!
! As a counterpart to the access to the source code and rights to copy,
! modify and redistribute granted by the license, users are provided only
! with a limited warranty and the software's author, the holder of the
! economic rights, and the successive licensors have only limited
! liability.
!
! In this respect, the user's attention is drawn to the risks associated
! with loading, using, modifying and/or developing or reproducing the
! software by the user in light of its specific status of free software,
! that may mean that it is complicated to manipulate, and that also
! therefore means that it is reserved for developers and experienced
! professionals having in-depth computer knowledge. Users are therefore
! encouraged to load and test the software's suitability as regards their
! requirements in conditions enabling the security of their systems and/or
! data to be ensured and, more generally, to use and operate it in the
! same conditions as regards security.
!
! The fact that you are presently reading this means that you have had
! knowledge of the CeCILL-C,CeCILL-B or CeCILL license and that you accept its terms.
!*/
!/*-----------------------------------------------------------------*/
!/*-----------------------------------------------------------------*/
!/* main program */
!/*-----------------------------------------------------------------*/
program f77getpositionrecordindex
implicit none
include 'f90calceph.h'
integer*8 peph
integer res
integer j
double precision firsttime, lasttime
integer countconst, frame, target1, center1,segid
include 'fopenfiles.h'
res=f90calceph_open(peph, trim(TOPSRCDIR) &
& //"../examples/example1.dat")
if (res.eq.1) then
countconst = f90calceph_getpositionrecordcount(peph)
do j=1, countconst
res = f90calceph_getpositionrecordindex(peph,j, &
& target1, center1, firsttime, lasttime, frame)
if (res.eq.0) then
write(*,*) 'res=0 for j=',j
stop 2
elseif (frame.ne.1) then
write(*,*)'not a valid frame',frame
stop 2
endif
res = f90calceph_getpositionrecordindex2(peph,j, &
& target1, center1, firsttime, lasttime, frame, &
& segid)
if (res.eq.0) then
stop 2
elseif (frame.ne.1) then
write(*,*)'bad frame',frame
stop 2
elseif (segid.ne.CALCEPH_SEGTYPE_ORIG_0) then
write(*,*)'bad segment type',segid
stop 2
endif
enddo
if (countconst.ne.12) then
write(*,*) 'countconst=',countconst
res = 0
endif
if (f90calceph_getpositionrecordindex(peph, 0, target1, &
& center1, firsttime, lasttime, frame).ne.0) then
write(*,*) 'invalid index 0'
res = 0
endif
if (f90calceph_getpositionrecordindex(peph, countconst+1, &
& target1, center1, firsttime, lasttime, frame).ne.0 &
& ) then
write(*,*) 'invalid index over', countconst+1
res = 0
endif
endif
! stop on sucess
if (res.eq.1) then
stop
endif
! stop on error
stop 2
end
|
