Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/601/A79 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-18
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/601/A79 VLTS. Properties of O dwarf in 30 Dor (Sabin-Sanjulian+, 2017)
*================================================================================
*The VLT-FLAMES Tarantula Survey.
*XXVI: Properties of the O-dwarf population in 30 Doradus.
* Sabin-Sanjulian C., Simon-Diaz S., Herrero A., Puls J., Schneider F.R.N.,
* Evans C.J., Garcia M., Najarro F., Brott I., Castro N., Crowther P.A.,
* de Koter A., de Mink S.E., Grafener G., Grin N.J., Holgado G., Langer N.,
* Lennon D.J., Maiz Apellaniz J., Ramirez-Agudelo O.H., Sana H., Taylor W.D.,
* Vink J.S., Walborn N.R.
* <Astron. Astrophys. 601, A79 (2017)>
* =2017A&A...601A..79S (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Stellar and wind parameters for stars analyzed
with HHe diagnostics
integer*4 nr__
parameter (nr__=85) ! Number of records
character*88 ar__ ! Full-size record
integer*4 VFTS ! VLT-FLAMES Tarantula Survey identifier
* (2011A&A...530A.108E)
character*17 SpType ! Spectral classification
integer*4 vsini ! (km/s) Projected rotational velocity
integer*4 Teff ! (K) Effective temperature
integer*4 e_Teff ! (K) Mean error on effective temperature
real*4 logg ! ([cm/s2]) Surface gravity
character*1 n_logg ! [*] Note on logg (1)
real*4 loggc ! ([cm/s2]) ? Rotation-corrected gravity (G1)
real*4 e_loggc ! ([cm/s2]) Mean error on gravity (minimum of 0.10)
character*1 l_Y_He ! Limit flag on Y(He)
real*4 Y_He ! Helium abundance N(He)/N(H) (error=0.02)
character*1 l_logQ ! Limit flag on logQ
real*4 logQ ! Wind-strength parameter (2)
real*4 e_logQ ! ? Mean error on logQ
character*1 l_logDmom ! Limit flag on logDmom
real*4 logDmom ! ? Wind momentum (error=0.4)
character*8 Comm ! Comments regarding possible binarity or
* multiplicity (G2)
*Note (1): * for stars with too high gravity (larger than 4.2 dex).
*Note (2): Q=(dM/dt)/(v_inf_R*)^3/2, where dM/dt is the mass-loss rate and
* v_inf_ the terminal wind velocity and R the stellar radius.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! *First estimates of stellar and wind parameters
for stars analyzed with HHeN diagnostics
integer*4 nr__1
parameter (nr__1=20) ! Number of records
character*64 ar__1 ! Full-size record
integer*4 VFTS_1 ! VLT-FLAMES Tarantula Survey identifier
* (2011A&A...530A.108E)
character*18 SpType_1 ! Spectral classification
integer*4 vsini_1 ! (km/s) Projected rotational velocity
integer*4 Teff_1 ! (K) Effective temperature
* (adopted error of 1500K)
real*4 logg_1 ! ([cm/s2]) Surface gravity (adopted error of 0.1)
real*4 loggc_1 ! ([cm/s2]) Rotation-corrected gravity (G1)
character*1 l_logQ_1 ! Limit flag on logQ
real*4 logQ_1 ! Wind-strength parameter
* (adopted error of 0.20)
character*1 l_logDmom_1 ! Limit flag on logDmom
real*4 logDmom_1 ! Wind momentum (error=0.4)
character*6 Comm_1 ! Comments regarding possible binarity or
* multiplicity (G2)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! Radii, luminosities and masses for stars
analyzed with HHe diagnostics
integer*4 nr__2
parameter (nr__2=85) ! Number of records
character*107 ar__2 ! Full-size record
integer*4 VFTS_2 ! VLT-FLAMES Tarantula Survey identifier
* (2011A&A...530A.108E)
character*17 SpType_2 ! Spectral classification
integer*4 Teff_2 ! (K) Effective temperature
real*4 loggc_2 ! ([cm/s2]) ? Rotation-corrected gravity (G1)
real*4 VMag ! (mag) ? Absolute magnitude in V (G3)
real*4 R ! (Rsun) ? Stellar radius
real*4 e_R ! (Rsun) ? Mean error in R
real*4 logL ! ([Lsun]) ? Stellar luminosity
real*4 e_logL ! ([Lsun]) ? Mean error in logL
real*4 Msp ! (Msun) ? Spectroscopic mass
real*4 e_Msp ! (Msun) ? Mean error in Msp
real*4 Mev_GT ! (Msun) ? Evolutionary mass from Kiel diagram
real*4 e_Mev_GT ! (Msun) ? Mean error on Mev(GT)
real*4 Mev_LT ! (Msun) ? Evolutionary mass from HR diagram
real*4 e_Mev_LT ! (Msun) ? Mean error on Mev(LT)
real*4 Mev_B ! (Msun) ? Evolutionary mass from BONNSAI
* (Schneider et al., 2014A&A...570A..66S)
real*4 e_Mev_B ! (Msun) ? Mean error in Mev(B)
character*8 Comm_2 ! Comments regarding possible binarity or
* multiplicity (G2)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb2.dat' ! First estimates of radii, luminosities and
masses for stars analyzed with HHeN diagnostics
integer*4 nr__3
parameter (nr__3=20) ! Number of records
character*84 ar__3 ! Full-size record
integer*4 VFTS_3 ! VLT-FLAMES Tarantula Survey identifier
* (2011A&A...530A.108E)
character*18 SpType_3 ! Spectral classification
integer*4 Teff_3 ! (K) Effective temperature
real*4 loggc_3 ! ([cm/s2]) Rotation-corrected gravity (G1)
real*4 VMag_1 ! (mag) Absolute magnitude in V (G3)
real*4 R_1 ! (Rsun) Stellar radius
real*4 logL_1 ! ([Lsun]) Stellar luminosity
real*4 Msp_1 ! (Msun) Spectroscopic mass
character*1 l_Mev_GT ! ? Limit flag on Mev(GT)
integer*4 Mev_GT_1 ! (Msun) Evolutionary mass from Kiel diagram
integer*4 Mev_LT_1 ! (Msun) Evolutionary mass from HR diagram
real*4 Mev_B_1 ! (Msun) Evolutionary mass from BONNSAI
* ( Schneider et al., 2014A&A...570A..66S)
real*4 e_Mev_B_1 ! (Msun) Mean error in Mev(B)
character*6 Comm_3 ! Comments regarding possible binarity or
* multiplicity (G2)
C=============================================================================
C Loading file 'tablea1.dat' ! Stellar and wind parameters for stars analyzed
* with HHe diagnostics
C Format for file interpretation
1 format(
+ I3,1X,A17,1X,I3,1X,I5,1X,I4,1X,F4.2,1X,A1,1X,F4.2,1X,F4.2,1X,
+ A1,1X,F4.2,1X,A1,1X,F5.1,1X,F3.1,1X,A1,1X,F4.1,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,85
read(1,'(A88)')ar__
read(ar__,1)
+ VFTS,SpType,vsini,Teff,e_Teff,logg,n_logg,loggc,e_loggc,
+ l_Y_He,Y_He,l_logQ,logQ,e_logQ,l_logDmom,logDmom,Comm
if(ar__(45:48) .EQ. '') loggc = rNULL__
if(ar__(70:72) .EQ. '') e_logQ = rNULL__
if(ar__(76:79) .EQ. '') logDmom = rNULL__
c ..............Just test output...........
write(6,1)
+ VFTS,SpType,vsini,Teff,e_Teff,logg,n_logg,loggc,e_loggc,
+ l_Y_He,Y_He,l_logQ,logQ,e_logQ,l_logDmom,logDmom,Comm
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! *First estimates of stellar and wind parameters
* for stars analyzed with HHeN diagnostics
C Format for file interpretation
2 format(
+ I3,1X,A18,1X,I3,1X,I5,1X,F4.2,1X,F4.2,1X,A1,1X,F5.1,1X,A1,1X,
+ F4.1,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,20
read(1,'(A64)')ar__1
read(ar__1,2)
+ VFTS_1,SpType_1,vsini_1,Teff_1,logg_1,loggc_1,l_logQ_1,logQ_1,
+ l_logDmom_1,logDmom_1,Comm_1
c ..............Just test output...........
write(6,2)
+ VFTS_1,SpType_1,vsini_1,Teff_1,logg_1,loggc_1,l_logQ_1,logQ_1,
+ l_logDmom_1,logDmom_1,Comm_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb1.dat' ! Radii, luminosities and masses for stars
* analyzed with HHe diagnostics
C Format for file interpretation
3 format(
+ I3,1X,A17,1X,I5,1X,F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,
+ 1X,F3.1,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,85
read(1,'(A107)')ar__2
read(ar__2,3)
+ VFTS_2,SpType_2,Teff_2,loggc_2,VMag,R,e_R,logL,e_logL,Msp,
+ e_Msp,Mev_GT,e_Mev_GT,Mev_LT,e_Mev_LT,Mev_B,e_Mev_B,Comm_2
if(ar__2(29:32) .EQ. '') loggc_2 = rNULL__
if(ar__2(34:38) .EQ. '') VMag = rNULL__
if(ar__2(40:44) .EQ. '') R = rNULL__
if(ar__2(46:49) .EQ. '') e_R = rNULL__
if(ar__2(51:54) .EQ. '') logL = rNULL__
if(ar__2(56:59) .EQ. '') e_logL = rNULL__
if(ar__2(61:64) .EQ. '') Msp = rNULL__
if(ar__2(66:69) .EQ. '') e_Msp = rNULL__
if(ar__2(71:74) .EQ. '') Mev_GT = rNULL__
if(ar__2(76:79) .EQ. '') e_Mev_GT = rNULL__
if(ar__2(81:84) .EQ. '') Mev_LT = rNULL__
if(ar__2(86:89) .EQ. '') e_Mev_LT = rNULL__
if(ar__2(91:94) .EQ. '') Mev_B = rNULL__
if(ar__2(96:98) .EQ. '') e_Mev_B = rNULL__
c ..............Just test output...........
write(6,3)
+ VFTS_2,SpType_2,Teff_2,loggc_2,VMag,R,e_R,logL,e_logL,Msp,
+ e_Msp,Mev_GT,e_Mev_GT,Mev_LT,e_Mev_LT,Mev_B,e_Mev_B,Comm_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb2.dat' ! First estimates of radii, luminosities and
* masses for stars analyzed with HHeN diagnostics
C Format for file interpretation
4 format(
+ I3,1X,A18,1X,I5,1X,F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.1,1X,A1,
+ 1X,I3,1X,I3,1X,F5.1,1X,F4.1,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'tableb2.dat')
write(6,*) '....Loading file: tableb2.dat'
do i__=1,20
read(1,'(A84)')ar__3
read(ar__3,4)
+ VFTS_3,SpType_3,Teff_3,loggc_3,VMag_1,R_1,logL_1,Msp_1,
+ l_Mev_GT,Mev_GT_1,Mev_LT_1,Mev_B_1,e_Mev_B_1,Comm_3
c ..............Just test output...........
write(6,4)
+ VFTS_3,SpType_3,Teff_3,loggc_3,VMag_1,R_1,logL_1,Msp_1,
+ l_Mev_GT,Mev_GT_1,Mev_LT_1,Mev_B_1,e_Mev_B_1,Comm_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end