Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/682/A29 into FORTRAN code for reading data files line by line.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Oct-08
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/A+A/682/A29 Detection of the Yarkovsky effect on NEAs (Fenucci+, 2024)
*================================================================================
*An automated procedure for the detection of the Yarkovsky effect
*and results from the ESA NEO Coordination Centre.
* Fenucci M., Micheli M., Gianotto F., Faggioli L., Oliviero D., Porru A.,
* Rudawska R., Cano J.L., Conversi L., Moissl R.
* <Astron. Astrophys. 682, A29 (2024)>
* =2024A&A...682A..29F (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! Determinations of A2 with S/N>3
integer*4 nr__
parameter (nr__=463) ! Number of records
character*185 ar__ ! Full-size record
character*10 Asteroid ! Name of the asteroid
real*4 H ! (mag) Absolute magnitude
real*8 RMS ! RMS of normalized residuals, 7D OD
real*8 RMS6D ! RMS of normalized residuals, 6D OD
real*4 A2 ! (au/d2) Transversal acceleration component
real*4 e_A2 ! (au/d2) Error in A2
real*4 da_dt ! (au/Myr) Semi-major axis drift
real*4 e_da_dt ! (au/Myr) Error in da/dt
real*4 max_da_dt ! (au/Myr) Maximum da/dt from Monte Carlo model
real*8 SNR ! Signal-to-noise of A2 detection
character*4 FAccept ! [1 Rej.] Flag for acceptance of the
* detection
integer*4 NOptObs ! Number of optical observations
integer*4 NRejOpt ! Number of rejected optical observations
* in 7D OD
integer*4 NRej6D ! Number of rejected optical observations
* in 6D OD
integer*4 NRadObs ! Number of radar observations
integer*4 NRejRad ! Number of rejected radar obs in 7D OD
integer*4 NRejRad6D ! Number of rejected radar obs in 6D OD
integer*4 NOptOld ! Number of old observations
real*8 Dlow ! (m) 15-th percentile of diameter
real*8 Dmed ! (m) 50-th percentile of diameter
real*8 Dhigh ! (m) 85-th percentile of diameter
integer*4 ModFlag ! [0/1] Flag for model used in Monte Carlo
real*8 Prot ! (h) ?=-1 Rotation period of the asteroid
character*1 Tax ! Taxonomic complex of the asteroid
real*4 deltat ! (yr) Length of observational arc
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb2.dat' ! Comparisons with JPL data
integer*4 nr__1
parameter (nr__1=350) ! Number of records
character*119 ar__1 ! Full-size record
character*10 Asteroid_1 ! Name of the asteroid
real*4 A2_1 ! (au/d2) Transversal acceleration from NEOCC
real*4 e_A2_1 ! (au/d2) Error in A2 from NEOCC
real*8 SNR_1 ! Signal-to-noise of A2 from NEOCC
real*4 A2j ! (au/d2) ?=- Transversal acceleration from JPL
real*4 e_A2j ! (au/d2) ?=- Error in A2 from JPL
real*8 SNRj ! ?=- Signal-to-noise of A2 from JPL
integer*4 Model ! [1/3]?=- Dynamical model of JPL
real*8 Err1 ! ?=- Relative error of JPL determination
real*8 Err2 ! ?=- Relative error of NEOCC determination
real*8 chiA2 ! ?=- Chi value of determination
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb3.dat' ! Comparisons with JPL data not accepted or
not found by NEOCC
integer*4 nr__2
parameter (nr__2=64) ! Number of records
character*120 ar__2 ! Full-size record
character*10 Asteroid_2 ! Name of the asteroid
real*4 da_dt_NEOCC ! (au/Myr) ?=- Semi-major axis drift from NEOCC
real*4 e_da_dt_NEOCC ! (au/Myr) ?=- Error of NEOCC determination
real*8 SNR_NEOCC ! ?=- Signal-to-noise of NEOCC determination
real*8 YM ! ?=- Maximum da/dt from Monte Carlo model
real*4 da_dt_JPL ! (au/Myr) Semi-major axis drift from JPL
real*4 e_da_dt_JPL ! (au/Myr) Error of JPL determination
real*8 SNR_JPL ! Signal-to-noise of JPL determination
integer*4 Model_1 ! [1/3] Dynamical model used by JPL
C=============================================================================
C Loading file 'tableb1.dat' ! Determinations of A2 with S/N>3
C Format for file interpretation
1 format(
+ A10,1X,F6.3,1X,F8.6,1X,F8.6,1X,E12.6,1X,E11.6,1X,E12.6,1X,
+ E11.6,1X,E11.9,1X,F9.5,1X,A4,1X,I5,1X,I2,1X,I2,1X,I2,1X,I1,1X,
+ I2,1X,I4,1X,F9.3,1X,F9.3,F9.3,1X,I1,1X,F7.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,463
read(1,'(A185)')ar__
read(ar__,1)
+ Asteroid,H,RMS,RMS6D,A2,e_A2,da_dt,e_da_dt,max_da_dt,SNR,
+ FAccept,NOptObs,NRejOpt,NRej6D,NRadObs,NRejRad,NRejRad6D,
+ NOptOld,Dlow,Dmed,Dhigh,ModFlag,Prot,Tax,deltat
c ..............Just test output...........
write(6,1)
+ Asteroid,H,RMS,RMS6D,A2,e_A2,da_dt,e_da_dt,max_da_dt,SNR,
+ FAccept,NOptObs,NRejOpt,NRej6D,NRadObs,NRejRad,NRejRad6D,
+ NOptOld,Dlow,Dmed,Dhigh,ModFlag,Prot,Tax,deltat
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb2.dat' ! Comparisons with JPL data
C Format for file interpretation
2 format(
+ A10,1X,E12.5,1X,E11.5,1X,F9.5,1X,E13.5,1X,E9.5,1X,F9.5,1X,I1,
+ 1X,F13.10,1X,F11.9,1X,F11.9)
C Effective file loading
open(unit=1,status='old',file=
+'tableb2.dat')
write(6,*) '....Loading file: tableb2.dat'
do i__=1,350
read(1,'(A119)')ar__1
read(ar__1,2)
+ Asteroid_1,A2_1,e_A2_1,SNR_1,A2j,e_A2j,SNRj,Model,Err1,Err2,
+ chiA2
if (idig(ar__1(47:59)).EQ.0) A2j = rNULL__
if (idig(ar__1(61:69)).EQ.0) e_A2j = rNULL__
if (idig(ar__1(71:79)).EQ.0) SNRj = rNULL__
if (Model .EQ. 45) Model = iNULL__
if (idig(ar__1(83:95)).EQ.0) Err1 = rNULL__
if (idig(ar__1(97:107)).EQ.0) Err2 = rNULL__
if (idig(ar__1(109:119)).EQ.0) chiA2 = rNULL__
c ..............Just test output...........
write(6,2)
+ Asteroid_1,A2_1,e_A2_1,SNR_1,A2j,e_A2j,SNRj,Model,Err1,Err2,
+ chiA2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb3.dat' ! Comparisons with JPL data not accepted or
* not found by NEOCC
C Format for file interpretation
3 format(
+ A10,2X,E13.5,2X,E19.5,2X,F10.7,4X,F11.9,2X,E12.9,1X,E11.5,3X,
+ F9.6,8X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'tableb3.dat')
write(6,*) '....Loading file: tableb3.dat'
do i__=1,64
read(1,'(A120)')ar__2
read(ar__2,3)
+ Asteroid_2,da_dt_NEOCC,e_da_dt_NEOCC,SNR_NEOCC,YM,da_dt_JPL,
+ e_da_dt_JPL,SNR_JPL,Model_1
if (idig(ar__2(13:25)).EQ.0) da_dt_NEOCC = rNULL__
if (idig(ar__2(28:46)).EQ.0) e_da_dt_NEOCC = rNULL__
if (idig(ar__2(49:58)).EQ.0) SNR_NEOCC = rNULL__
if (idig(ar__2(63:73)).EQ.0) YM = rNULL__
c ..............Just test output...........
write(6,3)
+ Asteroid_2,da_dt_NEOCC,e_da_dt_NEOCC,SNR_NEOCC,YM,da_dt_JPL,
+ e_da_dt_JPL,SNR_JPL,Model_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end