Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/836/168 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 2024-Apr-16
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/ApJ/836/168 Abundances & RVs for stars near (or in) NGC6273 (Johnson+, 2017)
*================================================================================
*A chemical composition survey of the iron-complex globular cluster NGC 6273
*(M19).
* Johnson C.I., Caldwell N., Rich R.M., Mateo M., Bailey III J.I.,
* Clarkson W.I., Olszewski E.W., Walker M.G.
* <Astrophys. J., 836, 168-168 (2017)>
* =2017ApJ...836..168J (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Observing log
integer*4 nr__
parameter (nr__=18) ! Number of records
character*66 ar__ ! Full-size record
character*1 Obs (nr__) ! Observation (P=Photometry; S=Spectroscopy)
character*2 Field (nr__) ! Field (1)
character*13 Tel (nr__) ! Telescope/Instrument
character*9 Setup (nr__) ! Instrument setup
character*14 Date (nr__) ! ("Y/M/D") UT Date of Observation
character*22 Exp (nr__) ! (s) Exposure sequence(s)
*Note (1): Fields with different designations indicate different telescope
* pointings. The "a" and "b" designations for the VLT--FLAMES setups
* correspond to the "HIERARCH ESO OBS NAME" keyword in the original image
* headers. The "a" and "b" fields with the same numbers typically
* observed the same stars, but sometimes with different fibers.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Star identifiers, coordinates, photometry, and
radial velocities for NGC 6273 members
integer*4 nr__1
parameter (nr__1=333) ! Number of records
character*69 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Set (nr__1) ! Sample (G1)
character*16 v2MASS (nr__1) ! 2MASS identifier (HHMMSSss+DDMMSSs; J2000)
character*1 f_2MASS (nr__1) ! Observation notes (G2)
real*8 RAdeg (nr__1) ! (deg) Right Ascension, decimal degrees (J2000) (1)
real*8 DEdeg (nr__1) ! (deg) Declination, decimal degrees (J2000) (1)
real*4 Jmag (nr__1) ! (mag) [10.3/16.2] 2MASS J magnitude (1)
real*4 Ksmag (nr__1) ! (mag) [9.3/16.5] 2MASS Ks magnitude (1)
real*8 HRV (nr__1) ! (km/s) [-351/347] Heliocentric radial velocity
* (range between ~120-165km/s for the members)
real*4 e_HRV (nr__1) ! (km/s) [0.1/6] Uncertainty in HRV
*Note (1): Source positions and J, Ks magnitudes taken from 2MASS.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Star identifiers, coordinates, photometry, and
radial velocities for non-members
integer*4 nr__2
parameter (nr__2=515) ! Number of records
character*69 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Set_1 (nr__2) ! Sample (G1)
character*16 v2MASS_1 (nr__2) ! 2MASS identifier (HHMMSSss+DDMMSSs; J2000)
character*1 f_2MASS_1 (nr__2) ! Observation notes (G2)
real*8 RAdeg_1 (nr__2) ! (deg) Right Ascension, decimal degrees (J2000) (1)
real*8 DEdeg_1 (nr__2) ! (deg) Declination, decimal degrees (J2000) (1)
real*4 Jmag_1 (nr__2) ! (mag) [10.3/16.2] 2MASS J magnitude (1)
real*4 Ksmag_1 (nr__2) ! (mag) [9.3/16.5] 2MASS Ks magnitude (1)
real*8 HRV_1 (nr__2) ! (km/s) [-351/347] Heliocentric radial velocity
* (range between ~120-165km/s for the members)
real*4 e_HRV_1 (nr__2) ! (km/s) [0.1/6] Uncertainty in HRV
*Note (1): Source positions and J, Ks magnitudes taken from 2MASS.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Model atmosphere parameters for NGC 6273 members
integer*4 nr__3
parameter (nr__3=59) ! Number of records
character*39 ar__3 ! Full-size record
character*16 v2MASS_2 (nr__3) ! 2MASS identifier (HHMMSSss+DDMMSSs; J2000)
character*1 n_2MASS (nr__3) ! [abd] Observation notes (G2)
integer*4 Teff (nr__3) ! (K) [4200/4900]? Effective temperature
real*4 log_g (nr__3) ! (cm/s2) [0.2/2.4]? Surface gravity
real*4 v_Fe_H_ (nr__3) ! ([Sun]) [-2/-1]? Metallicity
real*4 xi (nr__3) ! (km/s) [1.3/2.3]? Microturbulent velocity
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Chemical abundances for NGC 6273 members: Na-Cr
integer*4 nr__4
parameter (nr__4=59) ! Number of records
character*82 ar__4 ! Full-size record
character*16 v2MASS_3 (nr__4) ! 2MASS identifier (HHMMSSss+DDMMSSs; J2000)
character*1 f_2MASS_2 (nr__4) ! [a-c] Observation notes (G2)
real*4 v_NaI_Fe_ (nr__4) ! ([-]) [-0.2/0.7]? [NaI/Fe] abundance
real*4 e__NaI_Fe_ (nr__4) ! ([-]) ? Uncertainty for [NaI/Fe]
real*4 v_MgI_Fe_ (nr__4) ! ([-]) [0.2/0.6]? [MgI/Fe] abundance
real*4 e__MgI_Fe_ (nr__4) ! ([-]) ? Uncertainty for [MgI/Fe]
real*4 v_AlI_Fe_ (nr__4) ! ([-]) [0.1/1.3]? [AlI/Fe] abundance
real*4 e__AlI_Fe_ (nr__4) ! ([-]) ? Uncertainty for [AlI/Fe]
real*4 v_SiI_Fe_ (nr__4) ! ([-]) [-0.2/0.6]? [SiI/Fe] abundance
real*4 e__SiI_Fe_ (nr__4) ! ([-]) ? Uncertainty for [SiI/Fe]
real*4 v_CaI_Fe_ (nr__4) ! ([-]) [-0.07/0.5]? [CaI/Fe] abundance
real*4 e__CaI_Fe_ (nr__4) ! ([-]) ? Uncertainty for [CaI/Fe]
real*4 v_CrI_Fe_ (nr__4) ! ([-]) [-0.2/0.3]? [CrI/Fe] abundance
real*4 e__CrI_Fe_ (nr__4) ! ([-]) ? Uncertainty for [CrI/Fe]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Chemical abundances for NGC 6273 members: Fe-Eu
integer*4 nr__5
parameter (nr__5=59) ! Number of records
character*73 ar__5 ! Full-size record
character*16 v2MASS_4 (nr__5) ! 2MASS identifier (HHMMSSss+DDMMSSs; J2000)
character*1 f_2MASS_3 (nr__5) ! [a-c] Observation notes (G2)
real*4 v_FeI_H_ (nr__5) ! ([Sun]) [-2.1/-1]? [FeI/H] abundance
real*4 e__FeI_H_ (nr__5) ! ([Sun]) ? Uncertainty in [FeI/H]
real*4 v_FeII_H_ (nr__5) ! ([Sun]) [-2/-1.1]? [FeII/H] abundance
real*4 e__FeII_H_ (nr__5) ! ([Sun]) ? Uncertainty in [FeII/H]
real*4 v_NiI_Fe_ (nr__5) ! ([Sun]) [-0.4/0.2]? [NiI/Fe] abundance
real*4 e__NiI_Fe_ (nr__5) ! ([Sun]) ? Uncertainty in [NiI/Fe]
real*4 v_LaII_Fe_ (nr__5) ! ([Sun]) [-0.4/1.2]? [LaII/Fe] abundance
real*4 e__LaII_Fe_(nr__5) ! ([Sun]) ? Uncertainty in [LaII/Fe]
real*4 v_EuII_Fe_ (nr__5) ! ([Sun]) [0.1/0.8]? [EuII/Fe] abundance
real*4 e__EuII_Fe_(nr__5) ! ([Sun]) ? Uncertainty in [EuII/Fe]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Calcium triplet metallicity data
integer*4 nr__6
parameter (nr__6=274) ! Number of records
character*51 ar__6 ! Full-size record
integer*4 Set_2 (nr__6) ! Sample (G1)
character*16 v2MASS_5 (nr__6) ! 2MASS identifier (HHMMSSss+DDMMSSs; J2000)
character*1 f_2MASS_4 (nr__6) ! Observation notes (G2)
real*4 W_8542 (nr__6) ! (0.1nm) [1.5/4]? Equivalent width of CaII{lambda}8542
real*4 W_8662 (nr__6) ! (0.1nm) [1.2/3]? Equivalent width of CaII{lambda}8662
real*4 sW (nr__6) ! (0.1nm) [2.9/7]? Summed equivalent width
real*4 rW (nr__6) ! (0.1nm) [2.6/6]? Reduced equivalent width
real*4 v_Fe_H__1 (nr__6) ! ([Sun]) [-2.3/-0.5]? Metallicity values derived
* from Calcium Triplet
real*4 e__Fe_H_ (nr__6) ! ([Sun]) [0.1/0.6]? Uncertainty in [Fe/H]
C=============================================================================
C Loading file 'table1.dat' ! Observing log
C Format for file interpretation
1 format(A1,1X,A2,1X,A13,1X,A9,1X,A14,1X,A22)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,18
read(1,'(A66)')ar__
read(ar__,1)
+ Obs(i__),Field(i__),Tel(i__),Setup(i__),Date(i__),Exp(i__)
c ..............Just test output...........
write(6,1)
+ Obs(i__),Field(i__),Tel(i__),Setup(i__),Date(i__),Exp(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Star identifiers, coordinates, photometry, and
* radial velocities for NGC 6273 members
C Format for file interpretation
2 format(
+ I1,1X,A16,1X,A1,1X,F10.6,1X,F10.6,1X,F6.3,1X,F6.3,1X,F7.2,1X,
+ F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,333
read(1,'(A69)')ar__1
read(ar__1,2)
+ Set(i__),v2MASS(i__),f_2MASS(i__),RAdeg(i__),DEdeg(i__),
+ Jmag(i__),Ksmag(i__),HRV(i__),e_HRV(i__)
c ..............Just test output...........
write(6,2)
+ Set(i__),v2MASS(i__),f_2MASS(i__),RAdeg(i__),DEdeg(i__),
+ Jmag(i__),Ksmag(i__),HRV(i__),e_HRV(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Star identifiers, coordinates, photometry, and
* radial velocities for non-members
C Format for file interpretation
3 format(
+ I1,1X,A16,1X,A1,1X,F10.6,1X,F10.6,1X,F6.3,1X,F6.3,1X,F7.2,1X,
+ F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,515
read(1,'(A69)')ar__2
read(ar__2,3)
+ Set_1(i__),v2MASS_1(i__),f_2MASS_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),Jmag_1(i__),Ksmag_1(i__),HRV_1(i__),e_HRV_1(i__)
c ..............Just test output...........
write(6,3)
+ Set_1(i__),v2MASS_1(i__),f_2MASS_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),Jmag_1(i__),Ksmag_1(i__),HRV_1(i__),e_HRV_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Model atmosphere parameters for NGC 6273 members
C Format for file interpretation
4 format(A16,1X,A1,1X,I4,1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,59
read(1,'(A39)')ar__3
read(ar__3,4)
+ v2MASS_2(i__),n_2MASS(i__),Teff(i__),log_g(i__),v_Fe_H_(i__),
+ xi(i__)
if(ar__3(20:23) .EQ. '') Teff(i__) = iNULL__
if(ar__3(25:28) .EQ. '') log_g(i__) = rNULL__
if(ar__3(30:34) .EQ. '') v_Fe_H_(i__) = rNULL__
if(ar__3(36:39) .EQ. '') xi(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ v2MASS_2(i__),n_2MASS(i__),Teff(i__),log_g(i__),v_Fe_H_(i__),
+ xi(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Chemical abundances for NGC 6273 members: Na-Cr
C Format for file interpretation
5 format(
+ A16,1X,A1,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,59
read(1,'(A82)')ar__4
read(ar__4,5)
+ v2MASS_3(i__),f_2MASS_2(i__),v_NaI_Fe_(i__),e__NaI_Fe_(i__),
+ v_MgI_Fe_(i__),e__MgI_Fe_(i__),v_AlI_Fe_(i__),e__AlI_Fe_(i__),
+ v_SiI_Fe_(i__),e__SiI_Fe_(i__),v_CaI_Fe_(i__),e__CaI_Fe_(i__),
+ v_CrI_Fe_(i__),e__CrI_Fe_(i__)
if(ar__4(20:24) .EQ. '') v_NaI_Fe_(i__) = rNULL__
if(ar__4(26:29) .EQ. '') e__NaI_Fe_(i__) = rNULL__
if(ar__4(31:34) .EQ. '') v_MgI_Fe_(i__) = rNULL__
if(ar__4(36:39) .EQ. '') e__MgI_Fe_(i__) = rNULL__
if(ar__4(41:44) .EQ. '') v_AlI_Fe_(i__) = rNULL__
if(ar__4(46:49) .EQ. '') e__AlI_Fe_(i__) = rNULL__
if(ar__4(51:55) .EQ. '') v_SiI_Fe_(i__) = rNULL__
if(ar__4(57:60) .EQ. '') e__SiI_Fe_(i__) = rNULL__
if(ar__4(62:66) .EQ. '') v_CaI_Fe_(i__) = rNULL__
if(ar__4(68:71) .EQ. '') e__CaI_Fe_(i__) = rNULL__
if(ar__4(73:77) .EQ. '') v_CrI_Fe_(i__) = rNULL__
if(ar__4(79:82) .EQ. '') e__CrI_Fe_(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ v2MASS_3(i__),f_2MASS_2(i__),v_NaI_Fe_(i__),e__NaI_Fe_(i__),
+ v_MgI_Fe_(i__),e__MgI_Fe_(i__),v_AlI_Fe_(i__),e__AlI_Fe_(i__),
+ v_SiI_Fe_(i__),e__SiI_Fe_(i__),v_CaI_Fe_(i__),e__CaI_Fe_(i__),
+ v_CrI_Fe_(i__),e__CrI_Fe_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Chemical abundances for NGC 6273 members: Fe-Eu
C Format for file interpretation
6 format(
+ A16,1X,A1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F6.3,1X,F4.2,1X,
+ F5.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,59
read(1,'(A73)')ar__5
read(ar__5,6)
+ v2MASS_4(i__),f_2MASS_3(i__),v_FeI_H_(i__),e__FeI_H_(i__),
+ v_FeII_H_(i__),e__FeII_H_(i__),v_NiI_Fe_(i__),e__NiI_Fe_(i__),
+ v_LaII_Fe_(i__),e__LaII_Fe_(i__),v_EuII_Fe_(i__),
+ e__EuII_Fe_(i__)
if(ar__5(20:24) .EQ. '') v_FeI_H_(i__) = rNULL__
if(ar__5(26:29) .EQ. '') e__FeI_H_(i__) = rNULL__
if(ar__5(31:35) .EQ. '') v_FeII_H_(i__) = rNULL__
if(ar__5(37:40) .EQ. '') e__FeII_H_(i__) = rNULL__
if(ar__5(42:47) .EQ. '') v_NiI_Fe_(i__) = rNULL__
if(ar__5(49:52) .EQ. '') e__NiI_Fe_(i__) = rNULL__
if(ar__5(54:58) .EQ. '') v_LaII_Fe_(i__) = rNULL__
if(ar__5(60:63) .EQ. '') e__LaII_Fe_(i__) = rNULL__
if(ar__5(65:68) .EQ. '') v_EuII_Fe_(i__) = rNULL__
if(ar__5(70:73) .EQ. '') e__EuII_Fe_(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ v2MASS_4(i__),f_2MASS_3(i__),v_FeI_H_(i__),e__FeI_H_(i__),
+ v_FeII_H_(i__),e__FeII_H_(i__),v_NiI_Fe_(i__),e__NiI_Fe_(i__),
+ v_LaII_Fe_(i__),e__LaII_Fe_(i__),v_EuII_Fe_(i__),
+ e__EuII_Fe_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Calcium triplet metallicity data
C Format for file interpretation
7 format(
+ I1,1X,A16,1X,A1,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,274
read(1,'(A51)')ar__6
read(ar__6,7)
+ Set_2(i__),v2MASS_5(i__),f_2MASS_4(i__),W_8542(i__),
+ W_8662(i__),sW(i__),rW(i__),v_Fe_H__1(i__),e__Fe_H_(i__)
if(ar__6(22:25) .EQ. '') W_8542(i__) = rNULL__
if(ar__6(27:30) .EQ. '') W_8662(i__) = rNULL__
if(ar__6(32:35) .EQ. '') sW(i__) = rNULL__
if(ar__6(37:40) .EQ. '') rW(i__) = rNULL__
if(ar__6(42:46) .EQ. '') v_Fe_H__1(i__) = rNULL__
if(ar__6(48:51) .EQ. '') e__Fe_H_(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ Set_2(i__),v2MASS_5(i__),f_2MASS_4(i__),W_8542(i__),
+ W_8662(i__),sW(i__),rW(i__),v_Fe_H__1(i__),e__Fe_H_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end