Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/603/A56 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-20
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/603/A56 Stellar parameters and assumed wind parameters (Cazorla+, 2017)
*================================================================================
*Chemical abundances of fast-rotating massive stars.
*I. Description of the methods and individual results.
* Cazorla C., Morel T., Naze Y., Rauw G., Semaan T., Daflon S., Oey S.
* <Astron. Astrophys. 603, A56 (2017)>
* =2017A&A...603A..56C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cool.dat' ! Cooler stars stellar parameters (table F2, part 1)
integer*4 nr__
parameter (nr__=17) ! Number of records
character*155 ar__ ! Full-size record
character*8 Name ! Star name
integer*4 n_Name ! [1/2]? Note on Name (1)
integer*4 vsini ! (km/s) Projected rotational velocity
integer*4 e_vsini ! (km/s) Error on vsini
character*1 l_vmac ! Upper limit sign for vmac
integer*4 vmac ! (km/s) Macroturbulence velocity
character*12 mult ! Multiplicity from spectroscopy
character*3 runStat ! [Yes No] Runaway?
integer*4 Teff ! (K) Effective temperature
integer*4 e_Teff ! (K) Error on Teff
real*4 logg ! ([cm/s2]) Surface gravity
real*4 e_logg ! ([cm/s2]) Error on logg
real*4 loggC ! ([cm/s2]) Surface gravity corrected for centrifuge
* forces
real*4 e_loggC ! ([cm/s2]) Error on loggC
real*4 Y ! Helium abundance
real*4 e_Y ! Error on Y
character*1 l_logE_C ! Lower/upper limit in logE(C)
real*4 logE_C ! ? Carbon abundance
real*4 e_logE_C ! Error on logE(C)
character*1 l_logE_N ! Lower/upper limit in logE(N)
real*4 logE_N ! Nitrogen abundance
real*4 e_logE_N ! Error on logE(N)
character*1 l_logE_O1 ! Lower/upper limit in logE(O) in 4060-4082{AA}
real*4 logE_O1 ! ? Oxygen abundance in 4060-4082{AA}
real*4 logE_O2 ! Oxygen abundance in 4691-4709{AA}
real*4 logE_O ! Adopted oxygen abundance
real*4 e_logE_O ! Error on logE(O)
character*1 l__N_C_ ! Lower/upper limit in [N/C]
real*4 v_N_C_ ! ? N/C abundance ratio
real*4 e__N_C_ ! Error on [N/C]
character*1 l__N_O_ ! Lower/upper limit in [N/O]
real*4 v_N_O_ ! N/O abundance ratio
real*4 e__N_O_ ! Error on [N/O]
*Note (1): Note on Name ars follows:
* 1 = because the UCLES spectrum only covers the wavelength range ~4350-6800{AA},
* neither H{epsilon}, H{delta}, H{gamma}, He I 4026, nor the C and O lines
* in the 4060-4082{AA} region were used. In order to check the reliability
* of our results, we have determined the atmospheric parameters of HD 172367
* (whose spectral type is similar to HD 53755) considering either only one
* (H{beta}), or four (H{epsilon}, H{delta}, H{gamma}, H{beta}) Balmer lines.
* No significant differences were found between the two sets of results,
* ensuring that our parameter derivation for HD 53755 is secure.
* 2 = stars studied with both DETAIL/SURFACE and CMFGEN.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hot.dat' ! Hotter stars stellar parameters and assumed wind
parameters (table F2, part 2)
integer*4 nr__1
parameter (nr__1=23) ! Number of records
character*165 ar__1 ! Full-size record
character*8 Name_1 ! Star name
integer*4 n_Name_1 ! [1]? Note on Name (1)
integer*4 vsini_1 ! (km/s) Projected rotational velocity
integer*4 e_vsini_1 ! (km/s) Error on vsini
character*1 l_vmac_1 ! Upper limit sign for vmac
integer*4 vmac_1 ! (km/s) Macroturbulence velocity
character*12 mult_1 ! Multiplicity from spectroscopy
character*3 runStat_1 ! Runaway?
integer*4 Teff_1 ! (K) Effective temperature
integer*4 e_Teff_1 ! (K) Error on Teff
real*4 logg_1 ! ([cm/s2]) Surface gravity
real*4 e_logg_1 ! ([cm/s2]) Error on logg
real*4 loggC_1 ! ([cm/s2]) Surface gravity corrected for centrifuge
* forces
real*4 e_loggC_1 ! ([cm/s2]) Error on loggC
real*4 logdM_dt ! ([Msun/yr]) Mass-loss rate
integer*4 vinf ! (km/s) Terminal velocity
real*4 beta ! Wind beta parameter
real*4 Y_1 ! Helium abundance
real*4 e_Y_1 ! Error on y
character*2 l_logE_C_1 ! [<>= ] Lower/upper limit in logE(C)
real*4 logE_C_1 ! Carbon abundance
real*4 e_logE_C_1 ! Error on logE(C)
character*2 l_logE_N_1 ! [<>= ] Lower/upper limit in logE(N)
real*4 logE_N_1 ! Nitrogen abundance
real*4 e_logE_N_1 ! Error on logE(N)
character*2 l_logE_O ! [<>= ] Lower/upper limit in logE(O)
real*4 logE_O_1 ! Oxygen abundance
real*4 e_logE_O_1 ! Error on logE(O)
character*2 l__N_C__1 ! [<>= ] Lower/upper limit in [N/C]
real*4 v_N_C__1 ! N/C abundance ratio
real*4 e__N_C__1 ! Error on [N/C] (2)
character*2 l__N_O__1 ! [<>= ] Lower/upper limit in [N/O]
real*4 v_N_O__1 ! ? N/O abundance ratio
real*4 e__N_O__1 ! ? Error on [N/O] (2)
*Note (1): Note on Name as follows:
* 1 = atmospheric parameters and surface abundances must be considered with
* caution, as they were derived from a low S/N spectrum that only ranges
* from 4075 to 4920{AA} (i.e., with fewer diagnostic lines).
*Note (2): quadratic sum of the CNO abundance uncertainties.
C=============================================================================
C Loading file 'cool.dat' ! Cooler stars stellar parameters (table F2, part 1)
C Format for file interpretation
1 format(
+ A8,1X,I1,1X,I3,1X,I2,1X,A1,1X,I3,1X,A12,1X,A3,1X,I5,1X,I4,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.3,1X,F5.3,1X,A1,1X,F4.2,1X,
+ F4.2,1X,A1,1X,F4.2,1X,F4.2,1X,A1,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'cool.dat')
write(6,*) '....Loading file: cool.dat'
do i__=1,17
read(1,'(A155)')ar__
read(ar__,1)
+ Name,n_Name,vsini,e_vsini,l_vmac,vmac,mult,runStat,Teff,
+ e_Teff,logg,e_logg,loggC,e_loggC,Y,e_Y,l_logE_C,logE_C,
+ e_logE_C,l_logE_N,logE_N,e_logE_N,l_logE_O1,logE_O1,logE_O2,
+ logE_O,e_logE_O,l__N_C_,v_N_C_,e__N_C_,l__N_O_,v_N_O_,e__N_O_
if(ar__(10:10) .EQ. '') n_Name = iNULL__
if(ar__(87:90) .EQ. '') logE_C = rNULL__
if(ar__(111:114) .EQ. '') logE_O1 = rNULL__
if(ar__(133:137) .EQ. '') v_N_C_ = rNULL__
c ..............Just test output...........
write(6,1)
+ Name,n_Name,vsini,e_vsini,l_vmac,vmac,mult,runStat,Teff,
+ e_Teff,logg,e_logg,loggC,e_loggC,Y,e_Y,l_logE_C,logE_C,
+ e_logE_C,l_logE_N,logE_N,e_logE_N,l_logE_O1,logE_O1,logE_O2,
+ logE_O,e_logE_O,l__N_C_,v_N_C_,e__N_C_,l__N_O_,v_N_O_,e__N_O_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hot.dat' ! Hotter stars stellar parameters and assumed wind
* parameters (table F2, part 2)
C Format for file interpretation
2 format(
+ A8,1X,I1,1X,I3,1X,I2,1X,A1,1X,I3,1X,A12,1X,A3,1X,I5,1X,I4,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.1,1X,I4,1X,F4.2,1X,F5.3,1X,
+ F5.3,1X,A2,1X,F4.2,1X,F4.2,1X,A2,1X,F4.2,1X,F4.2,1X,A2,1X,
+ F4.2,1X,F4.2,1X,A2,1X,F5.2,1X,F4.2,1X,A2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'hot.dat')
write(6,*) '....Loading file: hot.dat'
do i__=1,23
read(1,'(A165)')ar__1
read(ar__1,2)
+ Name_1,n_Name_1,vsini_1,e_vsini_1,l_vmac_1,vmac_1,mult_1,
+ runStat_1,Teff_1,e_Teff_1,logg_1,e_logg_1,loggC_1,e_loggC_1,
+ logdM_dt,vinf,beta,Y_1,e_Y_1,l_logE_C_1,logE_C_1,e_logE_C_1,
+ l_logE_N_1,logE_N_1,e_logE_N_1,l_logE_O,logE_O_1,e_logE_O_1,
+ l__N_C__1,v_N_C__1,e__N_C__1,l__N_O__1,v_N_O__1,e__N_O__1
if(ar__1(10:10) .EQ. '') n_Name_1 = iNULL__
if(ar__1(156:160) .EQ. '') v_N_O__1 = rNULL__
if(ar__1(162:165) .EQ. '') e__N_O__1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,n_Name_1,vsini_1,e_vsini_1,l_vmac_1,vmac_1,mult_1,
+ runStat_1,Teff_1,e_Teff_1,logg_1,e_logg_1,loggC_1,e_loggC_1,
+ logdM_dt,vinf,beta,Y_1,e_Y_1,l_logE_C_1,logE_C_1,e_logE_C_1,
+ l_logE_N_1,logE_N_1,e_logE_N_1,l_logE_O,logE_O_1,e_logE_O_1,
+ l__N_C__1,v_N_C__1,e__N_C__1,l__N_O__1,v_N_O__1,e__N_O__1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end