- s int j;
double jd0=2451419; double dt1=0.5E0; t_calcephbin *peph; double G[6];
/* open the ephemeris file */ peph = calceph_open("example2_rotangmom.dat"); if (peph) {
calceph_prefetch(peph);
- calceph_rotangmom_unit(peph, jd0, dt1, NAIFID_EARTH,
CALCEPH_USE_NAIFID+CALCEPH_UNIT_SEC, G);
for(j=0; j<6; j++) printf("%23.16En", G[j]);
/* close the ephemeris file */ calceph_close(peph);
}
integer*8 peph
integer res
double precision jd0
double precision dt1
double precision G(6)
jd0 = 2451419
dt1 = 0.5D0
res = f90calceph_open(peph, "example2_rotangmom.dat")
if (res.eq.1) then
res = f90calceph_rotangmom_unit(peph,jd0, dt1, NAIFID_EARTH,
& CALCEPH_USE_NAIFID+CALCEPH_UNIT_SEC,
& G)
write(*,*) G
call f90calceph_close(peph)
endif