integer*8 peph
integer res, nvalue
double precision, allocatable :: radii
double precision svalue

res = f90calceph_open(peph, "example1.dat")
if (res.eq.1) then
     ! get the number of values
     nvalue = calceph_getconstantsd(peph, "BODY399_RADII", svalue)
     ! fill the array
     allocate(radii(1:nvalue))
     res = calceph_getconstantvd(peph, "BODY399_RADII", radii, nvalue)
     write(*,*) radii

  call f90calceph_close(peph)
endif