Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/517/A81 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/517/A81 2nd and 3d parameters of HB of globular clusters (Gratton+, 2010) *================================================================================ *The second and third parameters of the horizontal branch of globular clusters. * Gratton R.G., Carretta E., Bragaglia A., Lucatello S., S'Orazii V. * <Astron. Astrophys. 517, A81 (2010)> * =2010A&A...517A..81G C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Basic parameters for globular clusters considered in this paper integer*4 nr__ parameter (nr__=94) ! Number of records character*98 ar__ ! Full-size record character*7 NGC_IC ! Cluster name character*7 Name ! Other name real*4 v_Fe_H_ ! ([Sun]) Metallicity real*8 v_m_M_V ! (mag) Distance modulus calculated from V magnitude real*8 E_B_V ! (mag) Redenning real*8 VMAG ! (mag) Absolute V magnitude real*8 HBR ! ? Horizontal-branch ratio = (B-R)/(B+V+R) real*8 Age ! (Gyr) ? Estimated cluster age real*8 M_TO ! (solMass) ? Mass at turn-off real*8 M_RGB ! (solMass) ? Mass at the tip of the red-giant branch real*8 IQR1 ! ? Interquartile of the [Na/O] * (anti-)correlations real*8 IQR2 ! ? Interquartile of the [Mg/Al] * (anti-)correlations real*8 T_HB ! ([K]) ? Horizontal-branch effective temperature * (log scale) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Photometric data from the HST-snapshot (Piotto et al., 2002, Cat. J/A+A/391/945) integer*4 nr__1 parameter (nr__1=74) ! Number of records character*108 ar__1 ! Full-size record character*7 NGC_IC_1 ! Cluster name integer*4 n_HB ! Number of Horizon Branch stars integer*4 n_fld ! [0/125] Number of field stars (G1) real*4 F555 ! (mag) Median F555W magnitude real*4 e_F555 ! (mag) rms uncertainty on F555 real*4 v439_555 ! (mag) Median F439W-F555W colour index real*4 e_439_555 ! (mag) rms uncertainty on F439-F555 real*4 F555m ! (mag) Minimum F555W magnitude real*4 e_F555m ! (mag) rms uncertainty on F555m real*4 v439_555m ! (mag) Minimum F439W-F555W colour index real*4 e_439_555m ! (mag) rms uncertainty on (F439-F555)m real*4 v439_555M_1 ! (mag) Maximum F439W-F555W colour index real*4 e_439_555M_1 ! (mag) rms uncertainty on (F439-F555)M integer*4 n_RGB ! Number of RGB stars real*4 v439_555G ! (mag) F439W-F555W colour index for RGB real*4 e_439_555G ! (mag) rms uncertainty on (F439-F555)RGB c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Photometric data from ground-based data (Rosenberg et al. 2000, Cat. J/A+AS/144/5 and J/A+AS/145/451) integer*4 nr__2 parameter (nr__2=49) ! Number of records character*108 ar__2 ! Full-size record character*7 NGC_IC_2 ! Cluster name integer*4 n_HB_1 ! Number of HB stars in the cluster integer*4 n_fld_1 ! [0/47] Number of field stars (G1) real*4 Vmag_1 ! (mag) Median V magnitude real*4 e_Vmag ! (mag) rms uncertainty on real*4 V_I ! (mag) Median V-I colour index real*4 e_V_I ! (mag) rms uncertainty on V-I real*4 Vmagm ! (mag) Minimum V magnitude real*4 e_Vmagm ! (mag) rms uncertainty on real*4 v_V_I_m ! (mag) Minimum V magnitude real*4 e__V_I_m ! (mag) rms uncertainty on (V-I)m real*4 v_V_I_M_1 ! (mag) Maximum V-I colour index real*4 e__V_I_M_1 ! (mag) rms uncertainty on (V-I)M integer*4 n_RGB_1 ! Number of RGB stars real*4 v_V_I_G ! (mag) V-I colour index for RGB real*4 e__V_I_G ! (mag) rms uncertainty on (V-I)RGB c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Intrinsic HB photometric parameters for GCs with HST data integer*4 nr__3 parameter (nr__3=74) ! Number of records character*57 ar__3 ! Full-size record character*7 NGC_IC_3 ! Cluster name real*4 v_B_V_0G ! (mag) Intrinsic B-V colour index value, for RGB real*4 v_B_V_0m ! (mag) Minimim B-V colour index value real*4 v_B_V_0 ! (mag) Median B-V colour index value real*4 v_B_V_0M_1 ! (mag) Maximum B-V colour index value real*4 VMAGM_1 ! (mag) Maximum absolute V magnitude real*4 VMAG_2 ! (mag) Median absolute V magnitude c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Intrinsic HB photometric parameters for GCs with ground-based data integer*4 nr__4 parameter (nr__4=49) ! Number of records character*57 ar__4 ! Full-size record character*7 NGC_IC_4 ! Cluster name real*4 v_V_I_0G ! (mag) Intrinsic V-I colour index value, for RGB real*4 v_V_I_0m ! (mag) Mimimum V-I colour index value real*4 v_V_I_0 ! (mag) Median V-I colour index value real*4 v_V_I_0M_1 ! (mag) Maximum V-I colour index value real*4 VMAGM_2 ! (mag) Maximum absolute V magnitude real*4 VMAG_3 ! (mag) Median absolute V magnitude c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Masses on HB for GCs with HST data integer*4 nr__5 parameter (nr__5=74) ! Number of records character*61 ar__5 ! Full-size record character*7 NGC_IC_5 ! Cluster name real*4 Mmin ! (solMass) Mimimum mass real*4 e_Mmin ! (solMass) ? rms uncertainty on Mmin real*4 Mmed ! (solMass) Median mass real*4 e_Mmed ! (solMass) ? rms uncertainty on Mmed real*4 Mmax ! (solMass) ? Maximum mass real*4 e_Mmax ! (solMass) ? rms uncertainty on Mmax real*4 dMmed ! (solMass) ? Difference between median and minimum masses c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Masses on HB for GCs with ground-based data integer*4 nr__6 parameter (nr__6=49) ! Number of records character*61 ar__6 ! Full-size record character*7 NGC_IC_6 ! Cluster name real*4 Mmin_1 ! (solMass) Mimimum mass real*4 e_Mmin_1 ! (solMass) ? rms uncertainty on Mmin real*4 Mmed_1 ! (solMass) Median mass real*4 e_Mmed_1 ! (solMass) ? rms uncertainty on Mmed real*4 Mmax_1 ! (solMass) ? Maximum mass real*4 e_Mmax_1 ! (solMass) ? rms uncertainty on Mmax real*4 dMmed_1 ! (solMass) ? Difference between median and minimum masses c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Properties of the RR Lyrae variables for each cluster integer*4 nr__7 parameter (nr__7=66) ! Number of records character*73 ar__7 ! Full-size record character*7 NGC_IC_7 ! Cluster name real*4 S_RR ! ? RR Lyrae richness parameter (Harris, 1996, * Cat. VII/195) (1) integer*4 n_RR_H ! ? Number of RR Lyrae stars for HST sample real*4 f_RR_H ! ? Fraction of RR Lyrae stars for HST sample integer*4 n_RR_g ! ? Number of RR Lyrae for ground-based sample real*4 f_RR_g ! ? Fraction of RR Lyrae for ground-based sample real*4 f_RR ! ? Fraction of RR Lyrae stars real*4 v_Pab_ ! (d) ? Mean period integer*4 n_Pab ! ? Number of stars used to determine the period *Note (1): Number of RR Lyrae per unit cluster luminosity c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table10.dat' ! Clusters ages from main sequence fitting and from HB both for HST and ground based data-set integer*4 nr__8 parameter (nr__8=94) ! Number of records character*40 ar__8 ! Full-size record character*7 NGC_IC_8 ! Cluster name real*8 Age_1 ! (Gyr) ? Age from main-sequence fitting real*8 Age_HB_H ! (Gyr) ? Age from horizontal branch fitting, * for HST data-set real*8 Age_HB_g ! (Gyr) ? Age from horizontal branch fitting, * for ground-based data-set real*8 Age_Y ! (Gyr) ? Age value, determined with Ymed c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table11.dat' ! R' parameter as derived from both HST and ground based data integer*4 nr__9 parameter (nr__9=94) ! Number of records character*61 ar__9 ! Full-size record character*7 NGC_IC_9 ! Cluster name real*8 R__HST ! ? R' parameter as derived from HST data (1) real*4 e_R__HST ! ? rms uncertainty on R'(HST) real*8 R__GR ! ? R' parameter as derived from ground-based * data (1) real*4 e_R__GR ! ? rms uncertainty on R'(GR) real*8 v_R__ ! ? Mean R' parameter (1) (2) real*4 e__R__ ! ? rms uncertainty on <R'> real*8 Y__R__ ! ? He abundance for R'=<R'> real*8 e_Y__R__ ! ? rms uncertainty on Y(<R'>) *Note (1): R' = N_HB_/N'_RGB_, where N'_RGB_ is the number of stars on the * RGB brighter than V(HB)+1. *Note (2): Final adopted value <R'> is the weighted average between the two * estimates, when available. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table12.dat' ! Spread in He and mass required explaining the HB integer*4 nr__10 parameter (nr__10=64) ! Number of records character*55 ar__10 ! Full-size record character*7 NGC_IC_10 ! Cluster name real*8 dM_Y ! (solMass) Spread in mass {delta}M(Y) real*8 Ymed ! Median He abundance value real*8 Ymax ! Maximum He abundance value real*8 dY ! Difference between maximum and median Y values real*8 d_B_V ! (mag) Colour spread in {delta}(B-V) real*8 d_V_I ! (mag) Colour spread in {delta}(V-I) C============================================================================= C Loading file 'table1.dat' ! Basic parameters for globular clusters considered * in this paper C Format for file interpretation 1 format( + A7,2X,A7,F6.2,F7.2,F7.2,F7.2,F8.2,F7.2,F8.3,F8.3,F7.3,F8.3, + F9.3) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,94 read(1,'(A98)')ar__ read(ar__,1) + NGC_IC,Name,v_Fe_H_,v_m_M_V,E_B_V,VMAG,HBR,Age,M_TO,M_RGB, + IQR1,IQR2,T_HB if(ar__(44:51) .EQ. '') HBR = rNULL__ if(ar__(52:58) .EQ. '') Age = rNULL__ if(ar__(59:66) .EQ. '') M_TO = rNULL__ if(ar__(67:74) .EQ. '') M_RGB = rNULL__ if(ar__(75:81) .EQ. '') IQR1 = rNULL__ if(ar__(82:89) .EQ. '') IQR2 = rNULL__ if(ar__(90:98) .EQ. '') T_HB = rNULL__ c ..............Just test output........... write(6,1) + NGC_IC,Name,v_Fe_H_,v_m_M_V,E_B_V,VMAG,HBR,Age,M_TO,M_RGB, + IQR1,IQR2,T_HB c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Photometric data from the HST-snapshot * (Piotto et al., 2002, Cat. J/A+A/391/945) C Format for file interpretation 2 format( + A7,1X,I4,1X,I3,3X,F6.3,2X,F5.3,1X,F6.3,2X,F5.3,2X,F6.3,2X, + F5.3,1X,F6.3,2X,F5.3,1X,F6.3,2X,F5.3,1X,I4,1X,F6.3,2X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,74 read(1,'(A108)')ar__1 read(ar__1,2) + NGC_IC_1,n_HB,n_fld,F555,e_F555,v439_555,e_439_555,F555m, + e_F555m,v439_555m,e_439_555m,v439_555M_1,e_439_555M_1,n_RGB, + v439_555G,e_439_555G c ..............Just test output........... write(6,2) + NGC_IC_1,n_HB,n_fld,F555,e_F555,v439_555,e_439_555,F555m, + e_F555m,v439_555m,e_439_555m,v439_555M_1,e_439_555M_1,n_RGB, + v439_555G,e_439_555G c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Photometric data from ground-based data (Rosenberg * et al. 2000, Cat. J/A+AS/144/5 and J/A+AS/145/451) C Format for file interpretation 3 format( + A7,I4,1X,I3,4X,F6.3,2X,F5.3,1X,F6.3,2X,F5.3,2X,F6.3,2X,F5.3, + 1X,F6.3,2X,F5.3,1X,F6.3,2X,F5.3,1X,I4,1X,F6.3,2X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,49 read(1,'(A108)')ar__2 read(ar__2,3) + NGC_IC_2,n_HB_1,n_fld_1,Vmag_1,e_Vmag,V_I,e_V_I,Vmagm,e_Vmagm, + v_V_I_m,e__V_I_m,v_V_I_M_1,e__V_I_M_1,n_RGB_1,v_V_I_G,e__V_I_G c ..............Just test output........... write(6,3) + NGC_IC_2,n_HB_1,n_fld_1,Vmag_1,e_Vmag,V_I,e_V_I,Vmagm,e_Vmagm, + v_V_I_m,e__V_I_m,v_V_I_M_1,e__V_I_M_1,n_RGB_1,v_V_I_G,e__V_I_G c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Intrinsic HB photometric parameters for GCs * with HST data C Format for file interpretation 4 format(A7,3X,F5.3,3X,F6.3,3X,F6.3,2X,F6.3,3X,F5.3,3X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,74 read(1,'(A57)')ar__3 read(ar__3,4) + NGC_IC_3,v_B_V_0G,v_B_V_0m,v_B_V_0,v_B_V_0M_1,VMAGM_1,VMAG_2 c ..............Just test output........... write(6,4) + NGC_IC_3,v_B_V_0G,v_B_V_0m,v_B_V_0,v_B_V_0M_1,VMAGM_1,VMAG_2 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! Intrinsic HB photometric parameters for GCs * with ground-based data C Format for file interpretation 5 format(A7,2X,F6.3,3X,F6.3,3X,F6.3,2X,F6.3,3X,F5.3,3X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,49 read(1,'(A57)')ar__4 read(ar__4,5) + NGC_IC_4,v_V_I_0G,v_V_I_0m,v_V_I_0,v_V_I_0M_1,VMAGM_2,VMAG_3 c ..............Just test output........... write(6,5) + NGC_IC_4,v_V_I_0G,v_V_I_0m,v_V_I_0,v_V_I_0M_1,VMAGM_2,VMAG_3 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table6.dat' ! Masses on HB for GCs with HST data C Format for file interpretation 6 format( + A7,3X,F5.3,2X,F5.3,3X,F5.3,2X,F5.3,3X,F5.3,3X,F5.3,3X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,74 read(1,'(A61)')ar__5 read(ar__5,6) + NGC_IC_5,Mmin,e_Mmin,Mmed,e_Mmed,Mmax,e_Mmax,dMmed if(ar__5(18:22) .EQ. '') e_Mmin = rNULL__ if(ar__5(33:37) .EQ. '') e_Mmed = rNULL__ if(ar__5(41:45) .EQ. '') Mmax = rNULL__ if(ar__5(49:53) .EQ. '') e_Mmax = rNULL__ if(ar__5(57:61) .EQ. '') dMmed = rNULL__ c ..............Just test output........... write(6,6) + NGC_IC_5,Mmin,e_Mmin,Mmed,e_Mmed,Mmax,e_Mmax,dMmed c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! Masses on HB for GCs with ground-based data C Format for file interpretation 7 format( + A7,3X,F5.3,2X,F5.3,3X,F5.3,2X,F5.3,3X,F5.3,3X,F5.3,3X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,49 read(1,'(A61)')ar__6 read(ar__6,7) + NGC_IC_6,Mmin_1,e_Mmin_1,Mmed_1,e_Mmed_1,Mmax_1,e_Mmax_1, + dMmed_1 if(ar__6(18:22) .EQ. '') e_Mmin_1 = rNULL__ if(ar__6(33:37) .EQ. '') e_Mmed_1 = rNULL__ if(ar__6(41:45) .EQ. '') Mmax_1 = rNULL__ if(ar__6(49:53) .EQ. '') e_Mmax_1 = rNULL__ if(ar__6(57:61) .EQ. '') dMmed_1 = rNULL__ c ..............Just test output........... write(6,7) + NGC_IC_6,Mmin_1,e_Mmin_1,Mmed_1,e_Mmed_1,Mmax_1,e_Mmax_1, + dMmed_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! Properties of the RR Lyrae variables for each cluster C Format for file interpretation 8 format( + A7,2X,F5.1,3X,I2,7X,F4.2,6X,I3,8X,F4.2,4X,F4.2,2X,F5.3,3X,I4) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,66 read(1,'(A73)')ar__7 read(ar__7,8) + NGC_IC_7,S_RR,n_RR_H,f_RR_H,n_RR_g,f_RR_g,f_RR,v_Pab_,n_Pab if(ar__7(10:14) .EQ. '') S_RR = rNULL__ if(ar__7(18:19) .EQ. '') n_RR_H = iNULL__ if(ar__7(27:30) .EQ. '') f_RR_H = rNULL__ if(ar__7(37:39) .EQ. '') n_RR_g = iNULL__ if(ar__7(48:51) .EQ. '') f_RR_g = rNULL__ if(ar__7(56:59) .EQ. '') f_RR = rNULL__ if(ar__7(62:66) .EQ. '') v_Pab_ = rNULL__ if(ar__7(70:73) .EQ. '') n_Pab = iNULL__ c ..............Just test output........... write(6,8) + NGC_IC_7,S_RR,n_RR_H,f_RR_H,n_RR_g,f_RR_g,f_RR,v_Pab_,n_Pab c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table10.dat' ! Clusters ages from main sequence fitting and * from HB both for HST and ground based data-set C Format for file interpretation 9 format(A7,2X,F7.2,F8.3,F8.3,F8.3) C Effective file loading open(unit=1,status='old',file= +'table10.dat') write(6,*) '....Loading file: table10.dat' do i__=1,94 read(1,'(A40)')ar__8 read(ar__8,9)NGC_IC_8,Age_1,Age_HB_H,Age_HB_g,Age_Y if(ar__8(10:16) .EQ. '') Age_1 = rNULL__ if(ar__8(17:24) .EQ. '') Age_HB_H = rNULL__ if(ar__8(25:32) .EQ. '') Age_HB_g = rNULL__ if(ar__8(33:40) .EQ. '') Age_Y = rNULL__ c ..............Just test output........... write(6,9)NGC_IC_8,Age_1,Age_HB_H,Age_HB_g,Age_Y c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table11.dat' ! R' parameter as derived from both HST and * ground based data C Format for file interpretation 10 format(A7,F7.2,F6.2,F7.2,F6.2,F7.2,F6.2,F8.3,F7.3) C Effective file loading open(unit=1,status='old',file= +'table11.dat') write(6,*) '....Loading file: table11.dat' do i__=1,94 read(1,'(A61)')ar__9 read(ar__9,10) + NGC_IC_9,R__HST,e_R__HST,R__GR,e_R__GR,v_R__,e__R__,Y__R__, + e_Y__R__ if(ar__9(8:14) .EQ. '') R__HST = rNULL__ if(ar__9(15:20) .EQ. '') e_R__HST = rNULL__ if(ar__9(21:27) .EQ. '') R__GR = rNULL__ if(ar__9(28:33) .EQ. '') e_R__GR = rNULL__ if(ar__9(34:40) .EQ. '') v_R__ = rNULL__ if(ar__9(41:46) .EQ. '') e__R__ = rNULL__ if(ar__9(47:54) .EQ. '') Y__R__ = rNULL__ if(ar__9(55:61) .EQ. '') e_Y__R__ = rNULL__ c ..............Just test output........... write(6,10) + NGC_IC_9,R__HST,e_R__HST,R__GR,e_R__GR,v_R__,e__R__,Y__R__, + e_Y__R__ c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table12.dat' ! Spread in He and mass required explaining the HB C Format for file interpretation 11 format(A7,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3) C Effective file loading open(unit=1,status='old',file= +'table12.dat') write(6,*) '....Loading file: table12.dat' do i__=1,64 read(1,'(A55)')ar__10 read(ar__10,11)NGC_IC_10,dM_Y,Ymed,Ymax,dY,d_B_V,d_V_I c ..............Just test output........... write(6,11)NGC_IC_10,dM_Y,Ymed,Ymax,dY,d_B_V,d_V_I c .......End.of.Just test output........... end do close(1) C============================================================================= stop end