Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/380/1599 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-13
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/MNRAS/380/1599 Ephemeris of V382 Cyg and TU Mus (Qian+, 2007)
*================================================================================
*Evolutionary states of the two shortest period O-type overcontact binaries
*V382 Cyg and TU Mus.
* Qian S.-B., Yuan J.-Z., Liu L., He J.-J., Fernandez Lajus E., Kreiner J.M.
* <Mon. Not. R. Astron. Soc., 380, 1599-1607 (2007)>
* =2007MNRAS.380.1599Q
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Times of light minimum of V382 Cyg
integer*4 nr__
parameter (nr__=97) ! Number of records
character*62 ar__ ! Full-size record
real*8 HJD (nr__) ! (d) Heliocentric Julian date
real*4 e_HJD (nr__) ! (d) ? rms uncertainty on HJD
character*2 Min (nr__) ! [I II] First (I) or second (II) minimum
character*3 Obs (nr__) ! Observation type: CCD, Pe (photoelectric),
* Pg (photographic) or Vis (visual)
real*8 E (nr__) ! Ephemeris E value
real*8 v_O_C_1 (nr__) ! (d) Observed minus computed value
* for linear ephemeris (1)
real*8 v_O_C_2 (nr__) ! (d) Observed minus computed value
* for quadratic ephemeris (2)
real*8 Res (nr__) ! (d) Residuals
*Note (1): For V382 Cyg, Min.I=2436814.802+1.8855109dxE
* For TU Mus, Min.I=2441669.3255+1.3872766dxE
*Note (2):
* For V382 Cyg:
* Min.I=2436814.7610(+/-0.0021)+1.88550741d(+/-0.00000027)xE
* +1.13(+/-0.05)x10^-9^xE^2^
* +0.0148(+/-0.0020)sin[(0.0390{deg}xE+71.08{deg}(+/-0.06{deg})]
* For TU Mus:
* Min.I=2441669.3099(+/-0.0073)+1.38727095d(+/-0.00000113)xE
* +5.96(+/-0.93)x10^-10^xE^2^
* +0.0100(+/-0.0054)sin[0.0286{deg}xE+346.1{deg}(+/-9.9{deg})]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Times of light minimum of TU Mus
integer*4 nr__1
parameter (nr__1=38) ! Number of records
character*62 ar__1 ! Full-size record
real*8 HJD_1 (nr__1) ! (d) Heliocentric Julian date
real*4 e_HJD_1 (nr__1) ! (d) ? rms uncertainty on HJD
character*2 Min_1 (nr__1) ! [I II] First (I) or second (II) minimum
character*3 Obs_1 (nr__1) ! Observation type: CCD, Pe (photoelectric),
* Pg (photographic) or Vis (visual)
real*8 E_1 (nr__1) ! Ephemeris E value
real*8 v_O_C_1_1 (nr__1) ! (d) Observed minus computed value
* for linear ephemeris (1)
real*8 v_O_C_2_1 (nr__1) ! (d) Observed minus computed value
* for quadratic ephemeris (2)
real*8 Res_1 (nr__1) ! (d) Residuals
*Note (1): For V382 Cyg, Min.I=2436814.802+1.8855109dxE
* For TU Mus, Min.I=2441669.3255+1.3872766dxE
*Note (2):
* For V382 Cyg:
* Min.I=2436814.7610(+/-0.0021)+1.88550741d(+/-0.00000027)xE
* +1.13(+/-0.05)x10^-9^xE^2^
* +0.0148(+/-0.0020)sin[(0.0390{deg}xE+71.08{deg}(+/-0.06{deg})]
* For TU Mus:
* Min.I=2441669.3099(+/-0.0073)+1.38727095d(+/-0.00000113)xE
* +5.96(+/-0.93)x10^-10^xE^2^
* +0.0100(+/-0.0054)sin[0.0286{deg}xE+346.1{deg}(+/-9.9{deg})]
C=============================================================================
C Loading file 'table4.dat' ! Times of light minimum of V382 Cyg
C Format for file interpretation
1 format(
+ F12.4,1X,F6.4,2X,A2,1X,A3,1X,F8.1,3X,F7.4,1X,F7.4,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,97
read(1,'(A62)')ar__
read(ar__,1)
+ HJD(i__),e_HJD(i__),Min(i__),Obs(i__),E(i__),v_O_C_1(i__),
+ v_O_C_2(i__),Res(i__)
if(ar__(14:19) .EQ. '') e_HJD(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ HJD(i__),e_HJD(i__),Min(i__),Obs(i__),E(i__),v_O_C_1(i__),
+ v_O_C_2(i__),Res(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Times of light minimum of TU Mus
C Format for file interpretation
2 format(
+ F12.4,1X,F6.4,2X,A2,1X,A3,1X,F8.1,3X,F7.4,1X,F7.4,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,38
read(1,'(A62)')ar__1
read(ar__1,2)
+ HJD_1(i__),e_HJD_1(i__),Min_1(i__),Obs_1(i__),E_1(i__),
+ v_O_C_1_1(i__),v_O_C_2_1(i__),Res_1(i__)
if(ar__1(14:19) .EQ. '') e_HJD_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ HJD_1(i__),e_HJD_1(i__),Min_1(i__),Obs_1(i__),E_1(i__),
+ v_O_C_1_1(i__),v_O_C_2_1(i__),Res_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end