Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/479/541 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/479/541 VLT-FLAMES survey of massive stars (Hunter+, 2008)
*================================================================================
*The VLT-FLAMES survey of massive stars: atmospheric parameters and rotational
*velocity distributions for B-type stars in the Magellanic Clouds.
* Hunter I., Lennon D.J., Dufton P.L., Trundle C., Simon-Diaz S.,
* Smartt S.J., Ryans R.S.I., Evans C.J.
* <Astron. Astrophys., 479, 541-555 (2008)>
* =2008A&A...479..541H
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Radial distances of the objects in the NGC 346
FLAMES sample
integer*4 nr__
parameter (nr__=116) ! Number of records
character*17 ar__ ! Full-size record
character*12 Name ! Star designation, NGC 346-NNN (G1)
real*4 Dist ! (arcmin) Radial distance (1)
*Note (1): Radial distances of the objects in the NGC 346 FLAMES sample with
* respect to the centre of star formation taken as the position of
* object 435 from the Massey et al. (1989, Cat. J/AJ/98/1305)
* catalogue of OB stars in NGC 346; RA=00:59:04.49, DE=-72:10:24.7.
*
* The centre of NGC 346 adopted in Paper II (Evans et al., Cat.
* J/A+A/456/623) for calculation of radial distances was obtained from
* the SIMBAD database (operated at the CDS, Strasbourg, France) and is
* actually the centre of the ionised shell of gas, which is offset from
* the cluster centre.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Atmospheric parameters and vsini values for the
NGC 346 sample
integer*4 nr__1
parameter (nr__1=100) ! Number of records
character*74 ar__1 ! Full-size record
character*12 Name_1 ! Star name (G1)
character*1 m_Name ! [B] Multiplicity index on Name (1)
character*16 SpType ! MK spectral type
integer*4 Teff ! (K) ? Effective temperature
character*6 n_Teff ! Note on Teff when no measurement (2)
real*4 logg ! ([cm/s2]) ?=- Surface gravity
character*1 u_logg ! Uncertainty flag on logg (3)
character*1 l_vsini ! Limit flag on vsini
integer*4 vsini ! (km/s) Rotational velocity
character*1 u_vsini ! Uncertainty flag on vsini (3)
character*3 Met ! Method used (4)
real*4 Lum ! (solLum) ? Luminosity
integer*4 Mass ! (solMass) ? Mass
character*7 Com ! Comments
*Note (1): Where it was possible to measure the projected rotational
* velocity of the secondary object in a double lined spectroscopic
* binary system, a second entry is given with "B" appended to the
* object identifier.
*Note (2): We do not give atmospheric parameters for those objects listed as
* Be stars (Be), double lined spectroscopic binaries (SB2) or classified
* with an uncertain spectral type (SP?).
*Note (3): Values marked with colons indicate parameters where fitting has
* been difficult and the uncertainties may be larger than normal.
*Note (4): Method as follows:
* M06 = values directly obtained from Mokiem et al.
* (2006A&A...456.1131M, 2007A&A...465.1003M)
* T06 = values directly obtained from Trundle et al.
* (2007, Cat. <J/A+A/471/625>)
* H06 = values directly obtained from Hunter et al.
* (2007, Cat. <J/A+A/466/277>)
* He = indicates the methodology outlined in Sect. 3.2 has been used
* to derived the atmospheric parameters.
* A = the temperature has been assumed based on spectral type
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Atmospheric parameters and vsini values for the
NGC 330 sample
integer*4 nr__2
parameter (nr__2=104) ! Number of records
character*74 ar__2 ! Full-size record
character*12 Name_2 ! Star name (G1)
character*1 m_Name_1 ! [B] Multiplicity index on Name (1)
character*16 SpType_1 ! MK spectral type
integer*4 Teff_1 ! (K) ? Effective temperature
character*6 n_Teff_1 ! Note on Teff when no measurement (2)
real*4 logg_1 ! ([cm/s2]) ?=- Surface gravity
character*1 u_logg_1 ! Uncertainty flag on logg (3)
character*1 l_vsini_1 ! Limit flag on vsini
integer*4 vsini_1 ! (km/s) Rotational velocity
character*1 u_vsini_1 ! Uncertainty flag on vsini (3)
character*3 Met_1 ! Method used (4)
real*4 Lum_1 ! (solLum) ? Luminosity
integer*4 Mass_1 ! (solMass) ? Mass
character*7 Com_1 ! Comments
*Note (1): Where it was possible to measure the projected rotational
* velocity of the secondary object in a double lined spectroscopic
* binary system, a second entry is given with "B" appended to the
* object identifier.
*Note (2): We do not give atmospheric parameters for those objects listed as
* Be stars (Be), double lined spectroscopic binaries (SB2) or classified
* with an uncertain spectral type (SP?).
*Note (3): Values marked with colons indicate parameters where fitting has
* been difficult and the uncertainties may be larger than normal.
*Note (4): Method as follows:
* M06 = values directly obtained from Mokiem et al.
* (2006A&A...456.1131M, 2007A&A...465.1003M)
* T06 = values directly obtained from Trundle et al.
* (2007, Cat. <J/A+A/471/625>)
* H06 = values directly obtained from Hunter et al.
* (2007, Cat. <J/A+A/466/277>)
* He = indicates the methodology outlined in Sect. 3.2 has been used
* to derived the atmospheric parameters.
* A = the temperature has been assumed based on spectral type
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Atmospheric parameters and vsini values for the
N 11 (IC 2116) sample
integer*4 nr__3
parameter (nr__3=103) ! Number of records
character*74 ar__3 ! Full-size record
character*12 Name_3 ! Star name (G1)
character*1 m_Name_2 ! [B] Multiplicity index on Name (1)
character*16 SpType_2 ! MK spectral type
integer*4 Teff_2 ! (K) ? Effective temperature
character*6 n_Teff_2 ! Note on Teff when no measurement (2)
real*4 logg_2 ! ([cm/s2]) ?=- Surface gravity
character*1 u_logg_2 ! Uncertainty flag on logg (3)
character*1 l_vsini_2 ! Limit flag on vsini
integer*4 vsini_2 ! (km/s) Rotational velocity
character*1 u_vsini_2 ! Uncertainty flag on vsini (3)
character*3 Met_2 ! Method used (4)
real*4 Lum_2 ! (solLum) ? Luminosity
integer*4 Mass_2 ! (solMass) ? Mass
character*7 Com_2 ! Comments
*Note (1): Where it was possible to measure the projected rotational
* velocity of the secondary object in a double lined spectroscopic
* binary system, a second entry is given with "B" appended to the
* object identifier.
*Note (2): We do not give atmospheric parameters for those objects listed as
* Be stars (Be), double lined spectroscopic binaries (SB2) or classified
* with an uncertain spectral type (SP?).
*Note (3): Values marked with colons indicate parameters where fitting has
* been difficult and the uncertainties may be larger than normal.
*Note (4): Method as follows:
* M06 = values directly obtained from Mokiem et al.
* (2006A&A...456.1131M, 2007A&A...465.1003M)
* T06 = values directly obtained from Trundle et al.
* (2007, Cat. <J/A+A/471/625>)
* H06 = values directly obtained from Hunter et al.
* (2007, Cat. <J/A+A/466/277>)
* He = indicates the methodology outlined in Sect. 3.2 has been used
* to derived the atmospheric parameters.
* A = the temperature has been assumed based on spectral type
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Atmospheric parameters and vsini values for the
NGC 2004 sample
integer*4 nr__4
parameter (nr__4=101) ! Number of records
character*74 ar__4 ! Full-size record
character*12 Name_4 ! Star name (G1)
character*1 m_Name_3 ! [B] Multiplicity index on Name (1)
character*16 SpType_3 ! MK spectral type
integer*4 Teff_3 ! (K) ? Effective temperature
character*6 n_Teff_3 ! Note on Teff when no measurement (2)
real*4 logg_3 ! ([cm/s2]) ?=- Surface gravity
character*1 u_logg_3 ! Uncertainty flag on logg (3)
character*1 l_vsini_3 ! Limit flag on vsini
integer*4 vsini_3 ! (km/s) Rotational velocity
character*1 u_vsini_3 ! Uncertainty flag on vsini (3)
character*3 Met_3 ! Method used (4)
real*4 Lum_3 ! (solLum) ? Luminosity
integer*4 Mass_3 ! (solMass) ? Mass
character*7 Com_3 ! Comments
*Note (1): Where it was possible to measure the projected rotational
* velocity of the secondary object in a double lined spectroscopic
* binary system, a second entry is given with "B" appended to the
* object identifier.
*Note (2): We do not give atmospheric parameters for those objects listed as
* Be stars (Be), double lined spectroscopic binaries (SB2) or classified
* with an uncertain spectral type (SP?).
*Note (3): Values marked with colons indicate parameters where fitting has
* been difficult and the uncertainties may be larger than normal.
*Note (4): Method as follows:
* M06 = values directly obtained from Mokiem et al.
* (2006A&A...456.1131M, 2007A&A...465.1003M)
* T06 = values directly obtained from Trundle et al.
* (2007, Cat. <J/A+A/471/625>)
* H06 = values directly obtained from Hunter et al.
* (2007, Cat. <J/A+A/466/277>)
* He = indicates the methodology outlined in Sect. 3.2 has been used
* to derived the atmospheric parameters.
* A = the temperature has been assumed based on spectral type
C=============================================================================
C Loading file 'table1.dat' ! Radial distances of the objects in the NGC 346
* FLAMES sample
C Format for file interpretation
1 format(A12,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,116
read(1,'(A17)')ar__
read(ar__,1)Name,Dist
c ..............Just test output...........
write(6,1)Name,Dist
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Atmospheric parameters and vsini values for the
* NGC 346 sample
C Format for file interpretation
2 format(
+ A12,A1,1X,A16,1X,I5,A6,1X,F4.2,A1,1X,A1,I3,A1,1X,A3,1X,F4.2,
+ 1X,I2,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,100
read(1,'(A74)')ar__1
read(ar__1,2)
+ Name_1,m_Name,SpType,Teff,n_Teff,logg,u_logg,l_vsini,vsini,
+ u_vsini,Met,Lum,Mass,Com
if(ar__1(32:36) .EQ. '') Teff = iNULL__
if (idig(ar__1(44:47)).EQ.0) logg = rNULL__
if(ar__1(60:63) .EQ. '') Lum = rNULL__
if(ar__1(65:66) .EQ. '') Mass = iNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,m_Name,SpType,Teff,n_Teff,logg,u_logg,l_vsini,vsini,
+ u_vsini,Met,Lum,Mass,Com
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Atmospheric parameters and vsini values for the
* NGC 330 sample
C Format for file interpretation
3 format(
+ A12,A1,1X,A16,1X,I5,A6,1X,F4.2,A1,1X,A1,I3,A1,1X,A3,1X,F4.2,
+ 1X,I2,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,104
read(1,'(A74)')ar__2
read(ar__2,3)
+ Name_2,m_Name_1,SpType_1,Teff_1,n_Teff_1,logg_1,u_logg_1,
+ l_vsini_1,vsini_1,u_vsini_1,Met_1,Lum_1,Mass_1,Com_1
if(ar__2(32:36) .EQ. '') Teff_1 = iNULL__
if (idig(ar__2(44:47)).EQ.0) logg_1 = rNULL__
if(ar__2(60:63) .EQ. '') Lum_1 = rNULL__
if(ar__2(65:66) .EQ. '') Mass_1 = iNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,m_Name_1,SpType_1,Teff_1,n_Teff_1,logg_1,u_logg_1,
+ l_vsini_1,vsini_1,u_vsini_1,Met_1,Lum_1,Mass_1,Com_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Atmospheric parameters and vsini values for the
* N 11 (IC 2116) sample
C Format for file interpretation
4 format(
+ A12,A1,1X,A16,1X,I5,A6,1X,F4.2,A1,1X,A1,I3,A1,1X,A3,1X,F4.2,
+ 1X,I2,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,103
read(1,'(A74)')ar__3
read(ar__3,4)
+ Name_3,m_Name_2,SpType_2,Teff_2,n_Teff_2,logg_2,u_logg_2,
+ l_vsini_2,vsini_2,u_vsini_2,Met_2,Lum_2,Mass_2,Com_2
if(ar__3(32:36) .EQ. '') Teff_2 = iNULL__
if (idig(ar__3(44:47)).EQ.0) logg_2 = rNULL__
if(ar__3(60:63) .EQ. '') Lum_2 = rNULL__
if(ar__3(65:66) .EQ. '') Mass_2 = iNULL__
c ..............Just test output...........
write(6,4)
+ Name_3,m_Name_2,SpType_2,Teff_2,n_Teff_2,logg_2,u_logg_2,
+ l_vsini_2,vsini_2,u_vsini_2,Met_2,Lum_2,Mass_2,Com_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Atmospheric parameters and vsini values for the
* NGC 2004 sample
C Format for file interpretation
5 format(
+ A12,A1,1X,A16,1X,I5,A6,1X,F4.2,A1,1X,A1,I3,A1,1X,A3,1X,F4.2,
+ 1X,I2,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,101
read(1,'(A74)')ar__4
read(ar__4,5)
+ Name_4,m_Name_3,SpType_3,Teff_3,n_Teff_3,logg_3,u_logg_3,
+ l_vsini_3,vsini_3,u_vsini_3,Met_3,Lum_3,Mass_3,Com_3
if(ar__4(32:36) .EQ. '') Teff_3 = iNULL__
if (idig(ar__4(44:47)).EQ.0) logg_3 = rNULL__
if(ar__4(60:63) .EQ. '') Lum_3 = rNULL__
if(ar__4(65:66) .EQ. '') Mass_3 = iNULL__
c ..............Just test output...........
write(6,5)
+ Name_4,m_Name_3,SpType_3,Teff_3,n_Teff_3,logg_3,u_logg_3,
+ l_vsini_3,vsini_3,u_vsini_3,Met_3,Lum_3,Mass_3,Com_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end