Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/587/A48 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-Jun-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/587/A48 Lowell Photometric Database asteroid models (Durech+, 2016)
*================================================================================
*Asteroid models from the Lowell photometric database.
* Durech J., Hanus J, Oszkiewicz D., Vanco R.
* <Astron. Astrophys. 587, A48 (2016)>
* =2016A&A...587A..48D (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! List of new asteroid models derived from the
full period interval 2-100 hours
integer*4 nr__
parameter (nr__=259) ! Number of records
character*85 ar__ ! Full-size record
integer*4 Number (nr__) ! Asteroid number
character*21 Name (nr__) ! Asteroid name or designation
character*1 Ast (nr__) ! [*] '*' if published by Hanus et al.
* (2016A&A...586A.108H)
real*8 lambda1 (nr__) ! (deg) Ecliptic pole longitude (J2000.0) for model 1
real*8 beta1 (nr__) ! (deg) Ecliptic pole latitude (J2000.0) for model 1
real*8 lambda2 (nr__) ! (deg) ? Ecliptic pole longitude (J2000.0) for model 2
real*8 beta2 (nr__) ! (deg) ? Ecliptic pole latitude (J2000.0) for model 2
real*8 P (nr__) ! (h) Sidereal period of rotation
real*8 PLCDB (nr__) ! (h) ? Rotation period in the LCDB
real*4 Amin (nr__) ! (mag) ? Minimum lightcurve amplitude in LCDB
real*4 Amax (nr__) ! (mag) ? Maximum lightcurve amplitude in LCDB
character*2 U (nr__) ! Uncertainty code according to LCDB
integer*4 N (nr__) ! Number of photometric points
character*2 Method (nr__) ! Method used for period determination (1)
*Note (1): The method which was used to derive the unique rotation period:
* C - convex inversion, E - ellipsoids, CE - both methods gave the same
* unique period.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! List of new asteroid models derived from a
restricted period interval
integer*4 nr__1
parameter (nr__1=69) ! Number of records
character*85 ar__1 ! Full-size record
integer*4 Number_1 (nr__1) ! Asteroid number
character*21 Name_1 (nr__1) ! Asteroid name or designation
character*1 Ast_1 (nr__1) ! [*] '*' if published by Hanus et al.
* (2016A&A...586A.108H)
real*8 lambda1_1 (nr__1) ! (deg) Ecliptic pole longitude (J2000.0) for model 1
real*8 beta1_1 (nr__1) ! (deg) Ecliptic pole latitude (J2000.0) for model 1
real*8 lambda2_1 (nr__1) ! (deg) ? Ecliptic pole longitude (J2000.0) for model 2
real*8 beta2_1 (nr__1) ! (deg) ? Ecliptic pole latitude (J2000.0) for model 2
real*8 P_1 (nr__1) ! (h) Sidereal period of rotation
real*8 PLCDB_1 (nr__1) ! (h) ? Rotation period in the LCDB
real*4 Amin_1 (nr__1) ! (mag) ? Minimum lightcurve amplitude in LCDB
real*4 Amax_1 (nr__1) ! (mag) ? Maximum lightcurve amplitude in LCDB
character*2 U_1 (nr__1) ! Uncertainty code according to LCDB
integer*4 N_1 (nr__1) ! Number of photometric points
character*2 Method_1 (nr__1) ! Method used for period determination (1)
*Note (1): The method which was used to derive the unique rotation period:
* C - convex inversion, E - ellipsoids, CE - both methods gave the same
* unique period.
C=============================================================================
C Loading file 'table1.dat' ! List of new asteroid models derived from the
* full period interval 2-100 hours
C Format for file interpretation
1 format(
+ I5,2X,A21,1X,A1,1X,D3.2,1X,D3.2,1X,D3.2,1X,D3.2,1X,F9.6,1X,
+ F8.5,1X,F4.2,1X,F4.2,1X,A2,1X,I3,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,259
read(1,'(A85)')ar__
read(ar__,1)
+ Number(i__),Name(i__),Ast(i__),lambda1(i__),beta1(i__),
+ lambda2(i__),beta2(i__),P(i__),PLCDB(i__),Amin(i__),Amax(i__),
+ U(i__),N(i__),Method(i__)
if(ar__(40:42) .EQ. '') lambda2(i__) = rNULL__
if(ar__(44:46) .EQ. '') beta2(i__) = rNULL__
if(ar__(58:65) .EQ. '') PLCDB(i__) = rNULL__
if(ar__(67:70) .EQ. '') Amin(i__) = rNULL__
if(ar__(72:75) .EQ. '') Amax(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Number(i__),Name(i__),Ast(i__),lambda1(i__),beta1(i__),
+ lambda2(i__),beta2(i__),P(i__),PLCDB(i__),Amin(i__),Amax(i__),
+ U(i__),N(i__),Method(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! List of new asteroid models derived from a
* restricted period interval
C Format for file interpretation
2 format(
+ I5,2X,A21,1X,A1,1X,D3.2,1X,D3.2,1X,D3.2,1X,D3.2,1X,F9.6,1X,
+ F8.5,1X,F4.2,1X,F4.2,1X,A2,1X,I3,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,69
read(1,'(A85)')ar__1
read(ar__1,2)
+ Number_1(i__),Name_1(i__),Ast_1(i__),lambda1_1(i__),
+ beta1_1(i__),lambda2_1(i__),beta2_1(i__),P_1(i__),
+ PLCDB_1(i__),Amin_1(i__),Amax_1(i__),U_1(i__),N_1(i__),
+ Method_1(i__)
if(ar__1(40:42) .EQ. '') lambda2_1(i__) = rNULL__
if(ar__1(44:46) .EQ. '') beta2_1(i__) = rNULL__
if(ar__1(58:65) .EQ. '') PLCDB_1(i__) = rNULL__
if(ar__1(67:70) .EQ. '') Amin_1(i__) = rNULL__
if(ar__1(72:75) .EQ. '') Amax_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Number_1(i__),Name_1(i__),Ast_1(i__),lambda1_1(i__),
+ beta1_1(i__),lambda2_1(i__),beta2_1(i__),P_1(i__),
+ PLCDB_1(i__),Amin_1(i__),Amax_1(i__),U_1(i__),N_1(i__),
+ Method_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end