Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/607/A135 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 2026-Sep-12
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/607/A135 NGC 104, 6121 & 6809 AGB and RGB stars Na abundance (Wang+, 2017)
*================================================================================
*Sodium abundances of AGB and RGB stars in Galactic globular clusters
*II. Analysis and results of NGC 104, NGC 6121, and NGC 6809.
* Wang Y., Primas F., Charbonnel C., Van der Swaelmen M., Bono G.,
* Chantereau W., Zhao G.
* <Astron. Astrophys. 607, A135 (2017)>
* =2017A&A...607A.135W (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Basic information of our sample stars
integer*4 nr__
parameter (nr__=282) ! Number of records
character*176 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*4 NGC ! Globular cluster identification
character*9 StarID ! Star identification (1)
character*3 EvolPh ! Evolutionary phase
character*12 Inst ! Instrument used to obtain the spectrum
integer*4 RAh ! (h) Right Ascension J2000 (hours)
integer*4 RAm ! (min) Right Ascension J2000 (minutes)
real*4 RAs ! (s) Right Ascension J2000 (seconds)
character*1 DE_ ! Declination J2000 (sign)
integer*4 DEd ! (deg) Declination J2000 (degrees)
integer*4 DEm ! (arcmin) Declination J2000 (minutes)
real*4 DEs ! (arcsec) Declination J2000 (seconds)
real*4 Bmag ! (mag) Apparent B magnitude (2)
real*4 e_Bmag ! (mag) Error of B magnitude (2)
real*4 Vmag ! (mag) Apparent V magnitude (2)
real*4 e_Vmag ! (mag) Error of V magnitude (2)
real*4 Imag ! (mag) Apparent I magnitude (2)
real*4 e_Imag ! (mag) Error of I magnitude (2)
real*4 Jmag ! (mag) ? J magnitude (3)
real*4 e_Jmag ! (mag) ? Error of J magnitude (3)
real*4 Hmag ! (mag) ? H magnitude (3)
real*4 e_Hmag ! (mag) ? Error of H magnitude (3)
real*4 Kmag ! (mag) ? K magnitude (3)
real*4 e_Kmag ! (mag) ? Error of K magnitude (3)
real*4 RV ! (km/s) Barycentric radial velocity
*Note (1): The 'StarID' reports the original ID from the photometric catalogue,
* to which we added the suffix AGB/RGB to ease our own data handling.
*Note (2): The optical photometry (B, V and I) were obtained from the
* Johnson-Morgan photometric database which is part of the project
* described in Stetson (2000PASP..112..925S, 2005PASP..117.1325S).
*Note (3): The infrared J, H and K magnitudes were obtained from the 2MASS
* catalogue (Skrutskie et al., 2006, Cat. VII/233) by coordinates cross-matching.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Stellar parameters of our sample stars
integer*4 nr__1
parameter (nr__1=282) ! Number of records
character*101 ar__1 ! Full-size record
character*4 NGC_1 ! Globular cluster identification
character*9 StarID_1 ! Star identification
character*3 EvolPh_1 ! Evolutionary phase
integer*4 Teff ! (K) Effective temperature
real*4 s_Teff ! (K) Scatter of the temperature (1)
real*4 logg ! Surface gravity
real*4 Vturb ! (km/s) Microturbulence velocity
real*4 v_FeI_H_LTE ! ([-]) LTE [FeI/H] ratio
real*4 s_FeIline ! ([-]) Abundance dispersion of FeI lines
real*4 v_FeII_H_ ! ([-]) Line ratio [FeII/H]
real*4 s_FeIIline ! ([-]) Abundance dispersion of FeII lines
real*4 v_FeI_H_NLTE ! ([-]) NLTE [FeI/H] ratio
*Note (1): The s_Teff ({sigma}_Teff_) is the scatter of the temperatures derived
* from the five colors we considered.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Na abundances of our sample stars
integer*4 nr__2
parameter (nr__2=282) ! Number of records
character*50 ar__2 ! Full-size record
character*4 NGC_2 ! Globular cluster identification
character*9 StarID_2 ! Star identification
character*3 EvolPh_2 ! Evolutionary phase
real*4 v_Na_H_LTE ! ([-]) ? LTE [Na/H] (1)
real*4 v_Na_H_NLTE ! ([-]) ? NLTE [Na/H] (1)
real*4 v_Na_FeI_NLTE ! ([-]) ? NLTE [Na/FeI] (1)
*Note (1): The abundance value is left blank for the stars whose Na line/lines
* are of poor quality or saturate, and no reliable abundance could be derived.
C=============================================================================
C Loading file 'table3.dat' ! Basic information of our sample stars
C Format for file interpretation
1 format(
+ 1X,A4,3X,A9,4X,A3,2X,A12,4X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,
+ 1X,F5.2,4X,F6.3,2X,F6.4,3X,F6.3,2X,F6.4,3X,F6.3,2X,F6.4,3X,
+ F6.3,2X,F5.3,3X,F6.3,2X,F5.3,3X,F6.3,2X,F5.3,3X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,282
read(1,'(A176)')ar__
read(ar__,1)
+ NGC,StarID,EvolPh,Inst,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Bmag,
+ e_Bmag,Vmag,e_Vmag,Imag,e_Imag,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,
+ e_Kmag,RV
if(ar__(123:128) .EQ. '') Jmag = rNULL__
if(ar__(131:135) .EQ. '') e_Jmag = rNULL__
if(ar__(139:144) .EQ. '') Hmag = rNULL__
if(ar__(147:151) .EQ. '') e_Hmag = rNULL__
if(ar__(155:160) .EQ. '') Kmag = rNULL__
if(ar__(163:167) .EQ. '') e_Kmag = rNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ NGC,StarID,EvolPh,Inst,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Bmag,
+ e_Bmag,Vmag,e_Vmag,Imag,e_Imag,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,
+ e_Kmag,RV
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Stellar parameters of our sample stars
C Format for file interpretation
2 format(
+ 1X,A4,3X,A9,3X,A3,4X,I4,2X,F5.1,5X,F4.2,5X,F4.2,4X,F5.2,5X,
+ F4.2,4X,F5.2,5X,F4.2,4X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,282
read(1,'(A101)')ar__1
read(ar__1,2)
+ NGC_1,StarID_1,EvolPh_1,Teff,s_Teff,logg,Vturb,v_FeI_H_LTE,
+ s_FeIline,v_FeII_H_,s_FeIIline,v_FeI_H_NLTE
c ..............Just test output...........
write(6,2)
+ NGC_1,StarID_1,EvolPh_1,Teff,s_Teff,logg,Vturb,v_FeI_H_LTE,
+ s_FeIline,v_FeII_H_,s_FeIIline,v_FeI_H_NLTE
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Na abundances of our sample stars
C Format for file interpretation
3 format(1X,A4,3X,A9,3X,A3,4X,F5.2,4X,F5.2,4X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,282
read(1,'(A50)')ar__2
read(ar__2,3)
+ NGC_2,StarID_2,EvolPh_2,v_Na_H_LTE,v_Na_H_NLTE,v_Na_FeI_NLTE
if(ar__2(28:32) .EQ. '') v_Na_H_LTE = rNULL__
if(ar__2(37:41) .EQ. '') v_Na_H_NLTE = rNULL__
if(ar__2(46:50) .EQ. '') v_Na_FeI_NLTE = rNULL__
c ..............Just test output...........
write(6,3)
+ NGC_2,StarID_2,EvolPh_2,v_Na_H_LTE,v_Na_H_NLTE,v_Na_FeI_NLTE
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end