Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/610/A21 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 2024-Apr-20
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/610/A21 KiDS Survey for solar system objects mining (Mahlke+, 2018)
*================================================================================
*Mining the Kilo-Degree Survey for solar system objects.
* Mahlke M., Bouy H., Altieri B., Verdoes Kleijn G., Carry B., Bertin E.,
* de Jong J.T.A., Kuijken K., McFarland J., Valentijn E.
* <Astron. Astrophys. 610, A21 (2018)>
* =2018A&A...610A..21M (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'large.dat' ! Sample of 20221 SSOs with false- positive
content below 0.05%
integer*4 nr__
parameter (nr__=91454) ! Number of records
character*139 ar__ ! Full-size record
C J2000.0 position composed of: RAdeg DEdeg
integer*4 Number ! ID number of single observation
integer*4 Source ! ID number of SSO candidate
real*8 RAdeg ! (deg) Right ascension (J2000.0)
real*8 e_RAdeg ! (arcsec) Uncertainty of right ascension (J2000.0)
real*8 DEdeg ! (deg) Declination (J2000.0)
real*8 e_DEdeg ! (arcsec) Uncertainty of declination (J2000.0)
real*8 Epoch ! (yr) Observation epoch
character*1 Band ! [ugri] Observation band
real*8 mag ! (mag) Observed SSO magnitude (1)
real*8 e_mag ! (mag) Uncertainty of observed SSO magnitude
real*8 pmRA ! (arcsec/h) Proper motion in RA
real*8 pmDE ! (arcsec/h) Proper motion in DE
real*8 e_pm ! (arcsec/h) Uncertainty on proper motion
character*15 SkybotName ! SkyBoT match Name
real*4 Skybotmag ! ? SkyBoT match visual magnitude
*Note (1): The method is susceptible to artifacts in the images and the small.dat
* sample contains about 24% of false-positive.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'small.dat' ! Sample of 8069 SSOs with false-positive
content of about 24%
integer*4 nr__1
parameter (nr__1=25448) ! Number of records
character*139 ar__1 ! Full-size record
C J2000.0 position composed of: RAdeg DEdeg
integer*4 Number_1 ! ID number of single observation
integer*4 Source_1 ! ID number of SSO candidate
real*8 RAdeg_1 ! (deg) Right ascension (J2000.0)
real*8 e_RAdeg_1 ! (arcsec) Uncertainty of right ascension (J2000.0)
real*8 DEdeg_1 ! (deg) Declination (J2000.0)
real*8 e_DEdeg_1 ! (arcsec) Uncertainty of declination (J2000.0)
real*8 Epoch_1 ! (yr) Observation epoch
character*1 Band_1 ! [ugri] Observation band
real*8 mag_1 ! (mag) Observed SSO magnitude (1)
real*8 e_mag_1 ! (mag) Uncertainty of observed SSO magnitude
real*8 pmRA_1 ! (arcsec/h) Proper motion in RA
real*8 pmDE_1 ! (arcsec/h) Proper motion in DE
real*8 e_pm_1 ! (arcsec/h) Uncertainty on proper motion
character*15 SkybotName_1 ! SkyBoT match Name
real*4 Skybotmag_1 ! ? SkyBoT match visual magnitude
*Note (1): The method is susceptible to artifacts in the images and the small.dat
* sample contains about 24% of false-positive.
C=============================================================================
C Loading file 'large.dat' ! Sample of 20221 SSOs with false- positive
* content below 0.05%
C Format for file interpretation
1 format(
+ I6,1X,I5,1X,F10.6,1X,F8.6,1X,F10.6,1X,F8.6,1X,F11.6,1X,A1,
+ F10.6,1X,F8.6,1X,F11.6,1X,F11.6,1X,F8.6,1X,A15,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'large.dat')
write(6,*) '....Loading file: large.dat'
do i__=1,91454
read(1,'(A139)')ar__
read(ar__,1)
+ Number,Source,RAdeg,e_RAdeg,DEdeg,e_DEdeg,Epoch,Band,mag,
+ e_mag,pmRA,pmDE,e_pm,SkybotName,Skybotmag
if(ar__(136:139) .EQ. '') Skybotmag = rNULL__
c ..............Just test output...........
write(6,1)
+ Number,Source,RAdeg,e_RAdeg,DEdeg,e_DEdeg,Epoch,Band,mag,
+ e_mag,pmRA,pmDE,e_pm,SkybotName,Skybotmag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'small.dat' ! Sample of 8069 SSOs with false-positive
* content of about 24%
C Format for file interpretation
2 format(
+ I6,1X,I5,1X,F10.6,1X,F8.6,1X,F10.6,1X,F8.6,1X,F11.6,1X,A1,
+ F10.6,1X,F8.6,1X,F11.6,1X,F11.6,1X,F8.6,1X,A15,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'small.dat')
write(6,*) '....Loading file: small.dat'
do i__=1,25448
read(1,'(A139)')ar__1
read(ar__1,2)
+ Number_1,Source_1,RAdeg_1,e_RAdeg_1,DEdeg_1,e_DEdeg_1,Epoch_1,
+ Band_1,mag_1,e_mag_1,pmRA_1,pmDE_1,e_pm_1,SkybotName_1,
+ Skybotmag_1
if(ar__1(136:139) .EQ. '') Skybotmag_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Number_1,Source_1,RAdeg_1,e_RAdeg_1,DEdeg_1,e_DEdeg_1,Epoch_1,
+ Band_1,mag_1,e_mag_1,pmRA_1,pmDE_1,e_pm_1,SkybotName_1,
+ Skybotmag_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end