Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/682/A29 into FORTRAN code for loading all data files into arrays.
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 (nr__) ! Name of the asteroid
real*4 H (nr__) ! (mag) Absolute magnitude
real*8 RMS (nr__) ! RMS of normalized residuals, 7D OD
real*8 RMS6D (nr__) ! RMS of normalized residuals, 6D OD
real*4 A2 (nr__) ! (au/d2) Transversal acceleration component
real*4 e_A2 (nr__) ! (au/d2) Error in A2
real*4 da_dt (nr__) ! (au/Myr) Semi-major axis drift
real*4 e_da_dt (nr__) ! (au/Myr) Error in da/dt
real*4 max_da_dt (nr__) ! (au/Myr) Maximum da/dt from Monte Carlo model
real*8 SNR (nr__) ! Signal-to-noise of A2 detection
character*4 FAccept (nr__) ! [1 Rej.] Flag for acceptance of the
* detection
integer*4 NOptObs (nr__) ! Number of optical observations
integer*4 NRejOpt (nr__) ! Number of rejected optical observations
* in 7D OD
integer*4 NRej6D (nr__) ! Number of rejected optical observations
* in 6D OD
integer*4 NRadObs (nr__) ! Number of radar observations
integer*4 NRejRad (nr__) ! Number of rejected radar obs in 7D OD
integer*4 NRejRad6D (nr__) ! Number of rejected radar obs in 6D OD
integer*4 NOptOld (nr__) ! Number of old observations
real*8 Dlow (nr__) ! (m) 15-th percentile of diameter
real*8 Dmed (nr__) ! (m) 50-th percentile of diameter
real*8 Dhigh (nr__) ! (m) 85-th percentile of diameter
integer*4 ModFlag (nr__) ! [0/1] Flag for model used in Monte Carlo
real*8 Prot (nr__) ! (h) ?=-1 Rotation period of the asteroid
character*1 Tax (nr__) ! Taxonomic complex of the asteroid
real*4 deltat (nr__) ! (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 (nr__1) ! Name of the asteroid
real*4 A2_1 (nr__1) ! (au/d2) Transversal acceleration from NEOCC
real*4 e_A2_1 (nr__1) ! (au/d2) Error in A2 from NEOCC
real*8 SNR_1 (nr__1) ! Signal-to-noise of A2 from NEOCC
real*4 A2j (nr__1) ! (au/d2) ?=- Transversal acceleration from JPL
real*4 e_A2j (nr__1) ! (au/d2) ?=- Error in A2 from JPL
real*8 SNRj (nr__1) ! ?=- Signal-to-noise of A2 from JPL
integer*4 Model (nr__1) ! [1/3]?=- Dynamical model of JPL
real*8 Err1 (nr__1) ! ?=- Relative error of JPL determination
real*8 Err2 (nr__1) ! ?=- Relative error of NEOCC determination
real*8 chiA2 (nr__1) ! ?=- 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 (nr__2) ! Name of the asteroid
real*4 da_dt_NEOCC(nr__2) ! (au/Myr) ?=- Semi-major axis drift from NEOCC
real*4 e_da_dt_NEOCC(nr__2) ! (au/Myr) ?=- Error of NEOCC determination
real*8 SNR_NEOCC (nr__2) ! ?=- Signal-to-noise of NEOCC determination
real*8 YM (nr__2) ! ?=- Maximum da/dt from Monte Carlo model
real*4 da_dt_JPL (nr__2) ! (au/Myr) Semi-major axis drift from JPL
real*4 e_da_dt_JPL(nr__2) ! (au/Myr) Error of JPL determination
real*8 SNR_JPL (nr__2) ! Signal-to-noise of JPL determination
integer*4 Model_1 (nr__2) ! [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(i__),H(i__),RMS(i__),RMS6D(i__),A2(i__),e_A2(i__),
+ da_dt(i__),e_da_dt(i__),max_da_dt(i__),SNR(i__),FAccept(i__),
+ NOptObs(i__),NRejOpt(i__),NRej6D(i__),NRadObs(i__),
+ NRejRad(i__),NRejRad6D(i__),NOptOld(i__),Dlow(i__),Dmed(i__),
+ Dhigh(i__),ModFlag(i__),Prot(i__),Tax(i__),deltat(i__)
c ..............Just test output...........
write(6,1)
+ Asteroid(i__),H(i__),RMS(i__),RMS6D(i__),A2(i__),e_A2(i__),
+ da_dt(i__),e_da_dt(i__),max_da_dt(i__),SNR(i__),FAccept(i__),
+ NOptObs(i__),NRejOpt(i__),NRej6D(i__),NRadObs(i__),
+ NRejRad(i__),NRejRad6D(i__),NOptOld(i__),Dlow(i__),Dmed(i__),
+ Dhigh(i__),ModFlag(i__),Prot(i__),Tax(i__),deltat(i__)
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(i__),A2_1(i__),e_A2_1(i__),SNR_1(i__),A2j(i__),
+ e_A2j(i__),SNRj(i__),Model(i__),Err1(i__),Err2(i__),chiA2(i__)
if (idig(ar__1(47:59)).EQ.0) A2j(i__) = rNULL__
if (idig(ar__1(61:69)).EQ.0) e_A2j(i__) = rNULL__
if (idig(ar__1(71:79)).EQ.0) SNRj(i__) = rNULL__
if (Model(i__) .EQ. 45) Model(i__) = iNULL__
if (idig(ar__1(83:95)).EQ.0) Err1(i__) = rNULL__
if (idig(ar__1(97:107)).EQ.0) Err2(i__) = rNULL__
if (idig(ar__1(109:119)).EQ.0) chiA2(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Asteroid_1(i__),A2_1(i__),e_A2_1(i__),SNR_1(i__),A2j(i__),
+ e_A2j(i__),SNRj(i__),Model(i__),Err1(i__),Err2(i__),chiA2(i__)
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(i__),da_dt_NEOCC(i__),e_da_dt_NEOCC(i__),
+ SNR_NEOCC(i__),YM(i__),da_dt_JPL(i__),e_da_dt_JPL(i__),
+ SNR_JPL(i__),Model_1(i__)
if (idig(ar__2(13:25)).EQ.0) da_dt_NEOCC(i__) = rNULL__
if (idig(ar__2(28:46)).EQ.0) e_da_dt_NEOCC(i__) = rNULL__
if (idig(ar__2(49:58)).EQ.0) SNR_NEOCC(i__) = rNULL__
if (idig(ar__2(63:73)).EQ.0) YM(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Asteroid_2(i__),da_dt_NEOCC(i__),e_da_dt_NEOCC(i__),
+ SNR_NEOCC(i__),YM(i__),da_dt_JPL(i__),e_da_dt_JPL(i__),
+ SNR_JPL(i__),Model_1(i__)
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