Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/498/313 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-19
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/498/313 Photometry of 3 main belt asteroids (Marciniak+, 2009)
*================================================================================
*Photometry and models of selected main belt asteroids.
*VI. 160 Una, 747 Winchester, and 849 Ara.
* Marciniak A., Michalowski T., Hirsch R., Polinska M.,
* Kaminski K., Kwiatkowski T., Kryszczynska A., Behrend R.,
* Bernasconi L., Michalowski J., Starczewski S., Fagas M., Sobkowiak K.
* <Astron. Astrophys. 498, 313 (2009)>
* =2009A&A...498..313M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'objects.dat' ! Asteroid fundamental parameters
integer*4 nr__
parameter (nr__=3) ! Number of records
character*68 ar__ ! Full-size record
integer*4 Ast ! Asteroid number
character*10 Name ! Asteroid name
real*4 H ! (mag) Absolute magnitude H parameter (1)
real*4 Diam ! (km) IRAS diameter (1)
real*8 i ! (deg) Inclination (1)
real*8 e ! Eccentricity (1)
real*8 a ! (AU) Semimajor axis (1)
character*10 FileName ! Name of the file containing the light curve
*Note (1): Data from The Asteroid Orbital Elements Database, Cat. <B/astorb>.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'aspdat.dat' ! Aspect data of 160 Una, 747 Winchester and
849 Ara
integer*4 nr__1
parameter (nr__1=61) ! Number of records
character*67 ar__1 ! Full-size record
C J2000 position composed of: Elon Elat
integer*4 Ast_1 ! Asteroid number
character*10 Name_1 ! Asteroid name
integer*4 Obs_Y ! (yr) Year of an observation
integer*4 Obs_M ! (m) Month of an observation
real*4 Obs_D ! (d) Mid-time of an observation in days
real*4 r ! (AU) Distance from object to the Sun
real*4 delta ! (AU) Distance from object to the Earth
real*4 phAngle ! (deg) Sun-asteroid-Earth phase angle
real*4 Elon ! (deg) J2000 ecliptic longitude of object ({lambda})
real*4 Elat ! (deg) J2000 ecliptic latitude of object ({beta})
character*4 Site ! Observation site (1)
*Note (1): Observatory codes as follows:
* Bor = Borowiec
* EnO = Les Engarouines Observatory
* SAAO = South African Astronomical Observatory
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '160lcs.dat' ! Asteroid 160 Una individual lightcurves
integer*4 nr__2
parameter (nr__2=933) ! Number of records
character*111 ar__2 ! Full-size record
real*8 JD ! (d) JD epoch corrected for the light-time
* corresponding to the Earth-asteroid distance
real*4 br ! Relative brightness in intensity units,
* mean brightness of each lightcurve is unity
real*4 Sx ! (AU) x component of the vector from the asteroid to
* the Sun in J2000 ecliptic Cartesian coordinates
real*4 Sy ! (AU) y component of the vector from the asteroid to
* the Sun in J2000 ecliptic Cartesian coordinates
real*4 Sz ! (AU) z component of the vector from the asteroid to
* the Sun in J2000 ecliptic Cartesian coordinates
real*4 Ex ! (AU) x component of the vector from the asteroid to
* the Earth in J2000 ecliptic Cartesian coordinates
real*4 Ey ! (AU) y component of the vector from the asteroid to
* the Earth in J2000 ecliptic Cartesian coordinates
real*4 Ez ! (AU) z component of the vector from the asteroid to
* the Earth in J2000 ecliptic Cartesian coordinates
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '747lcs.dat' ! Asteroid 747 Winchester individual lightcurves
integer*4 nr__3
parameter (nr__3=2173) ! Number of records
character*111 ar__3 ! Full-size record
real*8 JD_1 ! (d) JD epoch corrected for the light-time
* corresponding to the Earth-asteroid distance
real*4 br_1 ! Relative brightness in intensity units,
* mean brightness of each lightcurve is unity
real*4 Sx_1 ! (AU) x component of the vector from the asteroid to
* the Sun in J2000 ecliptic Cartesian coordinates
real*4 Sy_1 ! (AU) y component of the vector from the asteroid to
* the Sun in J2000 ecliptic Cartesian coordinates
real*4 Sz_1 ! (AU) z component of the vector from the asteroid to
* the Sun in J2000 ecliptic Cartesian coordinates
real*4 Ex_1 ! (AU) x component of the vector from the asteroid to
* the Earth in J2000 ecliptic Cartesian coordinates
real*4 Ey_1 ! (AU) y component of the vector from the asteroid to
* the Earth in J2000 ecliptic Cartesian coordinates
real*4 Ez_1 ! (AU) z component of the vector from the asteroid to
* the Earth in J2000 ecliptic Cartesian coordinates
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '849lcs.dat' ! Asteroid 849 Ara individual lightcurves
integer*4 nr__4
parameter (nr__4=480) ! Number of records
character*111 ar__4 ! Full-size record
real*8 JD_2 ! (d) JD epoch corrected for the light-time
* corresponding to the Earth-asteroid distance
real*4 br_2 ! Relative brightness in intensity units,
* mean brightness of each lightcurve is unity
real*4 Sx_2 ! (AU) x component of the vector from the asteroid to
* the Sun in J2000 ecliptic Cartesian coordinates
real*4 Sy_2 ! (AU) y component of the vector from the asteroid to
* the Sun in J2000 ecliptic Cartesian coordinates
real*4 Sz_2 ! (AU) z component of the vector from the asteroid to
* the Sun in J2000 ecliptic Cartesian coordinates
real*4 Ex_2 ! (AU) x component of the vector from the asteroid to
* the Earth in J2000 ecliptic Cartesian coordinates
real*4 Ey_2 ! (AU) y component of the vector from the asteroid to
* the Earth in J2000 ecliptic Cartesian coordinates
real*4 Ez_2 ! (AU) z component of the vector from the asteroid to
* the Earth in J2000 ecliptic Cartesian coordinates
C=============================================================================
C Loading file 'objects.dat' ! Asteroid fundamental parameters
C Format for file interpretation
1 format(
+ I3,1X,A10,1X,F4.2,1X,F5.1,1X,F9.6,1X,F10.8,1X,F10.8,1X,A10)
C Effective file loading
open(unit=1,status='old',file=
+'objects.dat')
write(6,*) '....Loading file: objects.dat'
do i__=1,3
read(1,'(A68)')ar__
read(ar__,1)Ast,Name,H,Diam,i,e,a,FileName
c ..............Just test output...........
write(6,1)Ast,Name,H,Diam,i,e,a,FileName
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'aspdat.dat' ! Aspect data of 160 Una, 747 Winchester and
* 849 Ara
C Format for file interpretation
2 format(
+ I3,1X,A10,1X,I4,1X,I2,1X,F4.1,2X,F6.4,1X,F6.4,1X,F5.2,1X,F6.2,
+ 1X,F6.2,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'aspdat.dat')
write(6,*) '....Loading file: aspdat.dat'
do i__=1,61
read(1,'(A67)')ar__1
read(ar__1,2)
+ Ast_1,Name_1,Obs_Y,Obs_M,Obs_D,r,delta,phAngle,Elon,Elat,Site
c ..............Just test output...........
write(6,2)
+ Ast_1,Name_1,Obs_Y,Obs_M,Obs_D,r,delta,phAngle,Elon,Elat,Site
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file '160lcs.dat' ! Asteroid 160 Una individual lightcurves
C Format for file interpretation
3 format(
+ F14.6,1X,E12.7,1X,E13.8,1X,E13.8,1X,E13.6,1X,E13.8,1X,E13.8,
+ 1X,E13.8)
C Effective file loading
open(unit=1,status='old',file=
+'160lcs.dat')
write(6,*) '....Loading file: 160lcs.dat'
do i__=1,933
read(1,'(A111)')ar__2
read(ar__2,3)JD,br,Sx,Sy,Sz,Ex,Ey,Ez
c ..............Just test output...........
write(6,3)JD,br,Sx,Sy,Sz,Ex,Ey,Ez
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file '747lcs.dat' ! Asteroid 747 Winchester individual lightcurves
C Format for file interpretation
4 format(
+ F14.6,1X,E12.7,1X,E13.8,1X,E13.8,1X,E13.6,1X,E13.8,1X,E13.8,
+ 1X,E13.8)
C Effective file loading
open(unit=1,status='old',file=
+'747lcs.dat')
write(6,*) '....Loading file: 747lcs.dat'
do i__=1,2173
read(1,'(A111)')ar__3
read(ar__3,4)JD_1,br_1,Sx_1,Sy_1,Sz_1,Ex_1,Ey_1,Ez_1
c ..............Just test output...........
write(6,4)JD_1,br_1,Sx_1,Sy_1,Sz_1,Ex_1,Ey_1,Ez_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file '849lcs.dat' ! Asteroid 849 Ara individual lightcurves
C Format for file interpretation
5 format(
+ F14.6,1X,E12.7,1X,E13.8,1X,E13.8,1X,E13.6,1X,E13.8,1X,E13.8,
+ 1X,E13.8)
C Effective file loading
open(unit=1,status='old',file=
+'849lcs.dat')
write(6,*) '....Loading file: 849lcs.dat'
do i__=1,480
read(1,'(A111)')ar__4
read(ar__4,5)JD_2,br_2,Sx_2,Sy_2,Sz_2,Ex_2,Ey_2,Ez_2
c ..............Just test output...........
write(6,5)JD_2,br_2,Sx_2,Sy_2,Sz_2,Ex_2,Ey_2,Ez_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end