Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/522/A93 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/522/A93 Short-term variability of 29 minor planets (Thirouin+, 2010)
*================================================================================
*Short-term variability of a sample of 29 trans-Neptunian objects and Centaurs.
* Thirouin A., Ortiz J.L., Duffard R., Santos-Sanz P., Aceituno F.J.,
* Morales N.
* <Astron. Astrophys. 522, A93 (2010)>
* =2010A&A...522A..93T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'list.dat' ! List of studied objects
integer*4 nr__
parameter (nr__=29) ! Number of records
character*63 ar__ ! Full-size record
integer*4 Planet ! ? Planet number, if assigned (G1)
character*10 Name ! Planet name (G1)
real*4 H ! (mag) ? Absolute magnitude H parameter
real*8 i ! (deg) ? Orbit inclination
real*8 e ! ? Orbit eccentricity
real*8 a ! (AU) ? Orbit semimajor axis
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Dates, geometric and photometric data of the
observations
integer*4 nr__1
parameter (nr__1=239) ! Number of records
character*69 ar__1 ! Full-size record
integer*4 Planet_1 ! ? Planet number, if assigned (G1)
character*10 Name_1 ! Planet name (G1)
character*10 Date ! ("DD/MM/YYYY") Observation date
integer*4 Nimg ! Number of images
real*4 rh ! (AU) Topocentric distance
real*4 delta ! (AU) Heliocentric distance
real*4 phi ! (deg) Phase angle
character*1 Filt ! [R] Filter used
character*10 Tel ! Telescope: Calar Alto, INT or OSN (1)
*Note (1): Telescopes are:
* Calar Alto = 2.2m at Calar Alto Obs. (Almeria, Spain)
* INT = 2.5m Isaac Newton Telescope at El Roque de Los Muchachos (Spain)
* OSN = 1.5m telescope at Sierra Nevada Obs. (Granada, Spain)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Photometric results
integer*4 nr__2
parameter (nr__2=5377) ! Number of records
character*87 ar__2 ! Full-size record
integer*4 Planet_2 ! ? Planet number, if assigned (G1)
character*10 Name_2 ! Planet name (G1)
character*1 n_JD ! [*] * for the date of the image in which we
* performed a crude absolute calibration
real*8 JD ! (d) Julian date, not corrected for travel time
real*4 relmag ! (mag) Relative magnitude
real*4 e_relmag ! (mag) 1-{sigma} error associated to relmag
real*4 Rmag ! (mag) R magnitude
character*5 Filter ! Filter used during observational runs
* (Clear or R)
real*4 phi_1 ! (deg) Phase angle
real*4 rh_1 ! (AU) Topocentric distance
real*4 delta_1 ! (AU) Heliocentric distance
real*4 mR_1_1 ! (mag) Magnitude at 1AU from the Earth and at 1AU
* from the Sun
C=============================================================================
C Loading file 'list.dat' ! List of studied objects
C Format for file interpretation
1 format(1X,I6,2X,A10,5X,F5.2,2X,F9.6,1X,F10.8,1X,F11.8)
C Effective file loading
open(unit=1,status='old',file=
+'list.dat')
write(6,*) '....Loading file: list.dat'
do i__=1,29
read(1,'(A63)')ar__
read(ar__,1)Planet,Name,H,i,e,a
if(ar__(2:7) .EQ. '') Planet = iNULL__
if(ar__(25:29) .EQ. '') H = rNULL__
if(ar__(32:40) .EQ. '') i = rNULL__
if(ar__(42:51) .EQ. '') e = rNULL__
if(ar__(53:63) .EQ. '') a = rNULL__
c ..............Just test output...........
write(6,1)Planet,Name,H,i,e,a
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table1.dat' ! Dates, geometric and photometric data of the
* observations
C Format for file interpretation
2 format(
+ 1X,I6,2X,A10,4X,A10,1X,I3,1X,F6.3,1X,F6.3,1X,F4.2,1X,A1,1X,A10)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,239
read(1,'(A69)')ar__1
read(ar__1,2)Planet_1,Name_1,Date,Nimg,rh,delta,phi,Filt,Tel
if(ar__1(2:7) .EQ. '') Planet_1 = iNULL__
c ..............Just test output...........
write(6,2)Planet_1,Name_1,Date,Nimg,rh,delta,phi,Filt,Tel
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Photometric results
C Format for file interpretation
3 format(
+ 1X,I6,2X,A10,3X,A1,1X,F13.5,1X,F6.3,1X,F5.3,1X,F5.2,1X,A5,1X,
+ F4.2,1X,F6.3,1X,F6.3,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,5377
read(1,'(A87)')ar__2
read(ar__2,3)
+ Planet_2,Name_2,n_JD,JD,relmag,e_relmag,Rmag,Filter,phi_1,
+ rh_1,delta_1,mR_1_1
if(ar__2(2:7) .EQ. '') Planet_2 = iNULL__
c ..............Just test output...........
write(6,3)
+ Planet_2,Name_2,n_JD,JD,relmag,e_relmag,Rmag,Filter,phi_1,
+ rh_1,delta_1,mR_1_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end