Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/509/A3 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/509/A3 Earth Orientation Catalog 4 (EOC-4) (Vondrak+, 2010)
*================================================================================
*The Earth Orientation Catalog 4.
*An optical reference frame for monitoring Earth orientation in 20th century.
* Vondrak J., Stefka V.
* <Astron. Astrophys. 509, A3 (2010)>
* =2010A&A...509A...3V
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'eoc4.dat' ! IV version of Earth Orientation Catalogue
integer*4 nr__
parameter (nr__=4418) ! Number of records
character*287 ar__ ! Full-size record
C Position composed of: RAsec DEsec (Epoch=J2000)
integer*4 Seq ! Star number
real*8 RAsec ! (s) [0/86400] Right ascension (ICRS, Ep=J2000)
real*8 DEsec ! (arcsec) [-324000/324000] Declination (ICRS, Ep=J2000)
real*8 pmRA ! (arcsec/yr) Proper motion in right ascension multiplied
* by cos(declination) (pmRA*cos(DE))
real*8 pmDE ! (arcsec/yr) Proper motion in declination
real*4 EpRA ! (yr) Epoch of RA (offset 1900)
real*4 e_RAsec ! (mas) rms uncertainty in RA at EpRA multiplied
* by cos(declination)
real*4 e_pmRA ! (mas/yr) rms uncertainty in pmRA
real*4 EpDE ! (yr) Epoch of DE (offset 1900)
real*4 e_DEsec ! (mas) rms uncertainty in DEsec at EpDE
real*4 e_pmDE ! (mas/yr) rms uncertainty in pmDE
real*8 plx ! (arcsec) Trigonometric parallax
real*4 RV ! (km/s) Radial velocity
real*4 mag ! (mag) Apparent magnitude
real*4 P1 ! (yr) Period P1 (1/f1) (1)
real*4 P2 ! (yr) Period P2 (1/f2) (1)
real*8 c1R ! (arcsec) c1R coefficient (1)
real*8 s1R ! (arcsec) s1R coefficient (1)
real*8 c2R ! (arcsec) c2R coefficient (1)
real*8 s2R ! (arcsec) s2R coefficient (1)
real*8 c3R ! (arcsec) c3R coefficient (1)
real*8 s3R ! (arcsec) s3R coefficient (1)
real*8 c4R ! (arcsec) c4R coefficient (1)
real*8 s4R ! (arcsec) s4R coefficient (1)
real*8 c1D ! (arcsec) c1D coefficient (1)
real*8 s1D ! (arcsec) s1D coefficient (1)
real*8 c2D ! (arcsec) c2D coefficient (1)
real*8 s2D ! (arcsec) s2D coefficient (1)
real*8 c3D ! (arcsec) c3D coefficient (1)
real*8 s3D ! (arcsec) s3D coefficient (1)
real*8 c4D ! (arcsec) c4D coefficient (1)
real*8 s4D ! (arcsec) s4D coefficient (1)
real*4 e_ARA ! (mas) rms uncertainty in amplitudes of RA
real*4 e_ADE ! (mas) rms uncertainty in amplitudes of DE
integer*4 AE ! [0/3] Astrometric excellency (2)
character*3 CAT ! Astrometric Catalogue (3)
character*1 BS ! [ABPQ] Binary system (4)
integer*4 PF ! [1/3]? How the periods were found out (5)
*Note (1): Equations:
* RA(t) = RA+[pmRA*t+c1R*cos(arg1)+s1R*sin(arg1)+c2R*cos(2*arg1)
* +s2R*sin(2*arg1)+c3R*cos(arg2)+s3R*sin(arg2)
* +c4R*cos(2*arg2)+s4R*sin(2*arg2)]/15cos(DE)
* DE(t) = DE+pmDE*t+c1D*cos(arg1)+s1D*sin(arg1)+c2D*cos(2*arg1)+s2D*sin(2*arg1)
* +c3D*cos(arg2)+s3D*sin(arg2)+c4D*cos(2*arg2)+s4D*sin(2*arg2)
* where
* arg1 = 2*{pi}*t/P1,
* arg2 = 2*{pi}*t/P2,
* t = time in years from J2000.0
*Note (2): Astrometric excellency from 3 (highest rank) to 0 (lowest rank)
*Note (3): Catalogue that was combined with Earth Orientation observations
* as follows:
* F61 = Sixth Catalogue of Fundamental Stars I.
* F63 = Sixth Catalogue of Fundamental Stars III.
* GCH = Combination of Hipparcos with Boss' General Catalogue
* HIP = Hipparcos Catalogue
* PPM = Position and Proper Motion catalogue
* TYC = Tycho Catalogue
* T2H = Tycho-2
* XXX = 'local' catalogue
*Note (4): Binary system as follows:
* A = component A
* B = component B
* P = photocenter
* Q = other (404214 is the photocenter of stars 104214 and 104217)
*Note (5): How the periods were found out as follows:
* 1 = Period was found by spectral analysis only.
* 2 = Period was found both by spectral analysis and in the USNO Sixth
* Catalog of Orbits of Visual Binary Stars, their difference was
* smaller than a 20% relative tolerance, catalog's period was accepted.
* 3 = The difference between the two values for the period exceeded the
* tolerance, period found by spectral analysis was not
* statistically significant, but for the period taken over from
* the Sixth Catalog the condition of significance was met.
C=============================================================================
C Loading file 'eoc4.dat' ! IV version of Earth Orientation Catalogue
C Format for file interpretation
1 format(
+ I6,1X,F11.5,1X,F12.4,1X,F8.5,1X,F8.5,1X,F5.2,F6.2,1X,F5.2,1X,
+ F5.2,F6.2,1X,F5.2,1X,F7.4,1X,F6.1,1X,F5.2,1X,F6.2,F6.2,2X,
+ F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,
+ 2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,
+ F7.4,1X,F6.2,1X,F6.2,1X,I1,1X,A3,1X,A1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'eoc4.dat')
write(6,*) '....Loading file: eoc4.dat'
do i__=1,4418
read(1,'(A287)')ar__
read(ar__,1)
+ Seq,RAsec,DEsec,pmRA,pmDE,EpRA,e_RAsec,e_pmRA,EpDE,e_DEsec,
+ e_pmDE,plx,RV,mag,P1,P2,c1R,s1R,c2R,s2R,c3R,s3R,c4R,s4R,c1D,
+ s1D,c2D,s2D,c3D,s3D,c4D,s4D,e_ARA,e_ADE,AE,CAT,BS,PF
if(ar__(287:287) .EQ. '') PF = iNULL__
c ..............Just test output...........
write(6,1)
+ Seq,RAsec,DEsec,pmRA,pmDE,EpRA,e_RAsec,e_pmRA,EpDE,e_DEsec,
+ e_pmDE,plx,RV,mag,P1,P2,c1R,s1R,c2R,s2R,c3R,s3R,c4R,s4R,c1D,
+ s1D,c2D,s2D,c3D,s3D,c4D,s4D,e_ARA,e_ADE,AE,CAT,BS,PF
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end