Conversion of standardized ReadMe file for
file /./ftp/cats/V/117A/ 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-Mar-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. V/117A Geneva-Copenhagen Survey of Solar neighbourhood (Holmberg+, 2007) *================================================================================ *The Geneva-Copenhagen survey of the Solar neighbourhood. *I. Ages, metallicities, and kinematic properties of ~14000 F and G dwarfs. * Nordstrom B., Mayor M., Andersen J., Holmberg J., Pont F., Jorgensen B.R., * Olsen E.H., Udry S., Mowlavi N. * <Astron. Astrophys. 419, 989 (2004)> * =2004A&A...418..989N *II. New uvby calibrations and rediscussion of stellar ages, the G dwarf problem, *age-metallicity diagram, and heating mechanisms of the disk. * Holmberg J., Nordstroem B., Andersen J. * <Astron. Astrophys. 475, 519 (2007)> * =2007A&A...475..519H C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'newcat.dat' ! Main catalogue recalibrated data (2007) integer*4 nr__ parameter (nr__=16682) ! Number of records character*96 ar__ ! Full-size record C Position composed of: RAh RAm RAs DE- DEd DEm DEs integer*4 HIP (nr__) ! ? Hipparcos number if available character*11 Name (nr__) ! Object name (HD,BD,CD or CP) (1) character*4 Comp (nr__) ! Components included in photometry (2) integer*4 RAh (nr__) ! (h) Right ascension (ICRF) integer*4 RAm (nr__) ! (min) Right ascension (ICRF) real*4 RAs (nr__) ! (s) Right ascension (ICRF) character*1 DE_ (nr__) ! Declination sign (ICRF) integer*4 DEd (nr__) ! (deg) Declination (ICRF) integer*4 DEm (nr__) ! (arcmin) Declination (ICRF) integer*4 DEs (nr__) ! (arcsec) Declination (ICRF) real*4 logTe (nr__) ! ([K]) ? Effective temperature real*4 v_Fe_H_ (nr__) ! ([Sun]) ? Metallicity integer*4 Dist (nr__) ! (pc) ? Distance real*4 VMAG (nr__) ! (mag) ? Absolute magnitude real*4 Age (nr__) ! (Gyr) ? Age, in Gigayears real*4 clAge (nr__) ! (Gyr) ? Lower confidence limit on Age real*4 chAge (nr__) ! (Gyr) ? Upper confidence limit on Age integer*4 UVel (nr__) ! (km/s) ? Heliocentric space velocity component U (3) integer*4 VVel (nr__) ! (km/s) ? Heliocentric space velocity component V (3) integer*4 WVel (nr__) ! (km/s) ? Heliocentric space velocity component W (3) *Note (1): Object name given in order of preference. A second number * following a slash indicates that a double star with separate HD * numbers has been observed together. *Note (2): If the star is a member of a multiple system the component(s) * included in the photometry are identified here. *Note (3): Heliocentric space velocity components given in a right handed * coordinate system with U towards the galactic centre. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! First version of the catalog (2004) integer*4 nr__1 parameter (nr__1=16682) ! Number of records character*287 ar__1 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__1) ! (deg) Right Ascension J2000.0 real*8 DEdeg (nr__1) ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) integer*4 HIP_1 (nr__1) ! ? Hipparcos number if available character*11 Name_1 (nr__1) ! Object name (HD,BD,CD or CP) (1) character*4 Comp_1 (nr__1) ! Components included in photometry (2) character*1 fb (nr__1) ! Flag for binaries of all types (3) character*1 fs (nr__1) ! Flag for cool dwarfs south of -26deg integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000.0) (4) integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000.0) (4) real*4 RAs_1 (nr__1) ! (s) Right ascension (J2000.0) (4) character*1 DE__1 (nr__1) ! Declination sign (J2000.0) (4) integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000.0) (4) integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000.0) (4) integer*4 DEs_1 (nr__1) ! (arcsec) Declination (J2000.0) (4) integer*4 GLON (nr__1) ! (deg) Galactic longitude integer*4 GLAT (nr__1) ! (deg) Galactic latitude real*4 Vmag_1 (nr__1) ! (mag) Johnson V magnitude real*4 b_y (nr__1) ! (mag) Stromgren b-y colour real*4 Hbeta (nr__1) ! (mag) ? Stromgren Hbeta index real*4 E_b_y (nr__1) ! (mag) ? Colour excess (5) real*4 logTe_1 (nr__1) ! ([K]) ? Effective temperature (6) real*4 v_Fe_H__1 (nr__1) ! ([Sun]) ? Metallicity (7) integer*4 Dist_1 (nr__1) ! (pc) ? Distance (8) real*4 VMAG_2 (nr__1) ! (mag) ? Absolute magnitude real*4 dVMag (nr__1) ! (mag) ? Magnitude difference to ZAMS character*1 fr (nr__1) ! Source for the distance (8) character*1 fg (nr__1) ! Flag for suspected giants (9) real*4 Age_1 (nr__1) ! (Gyr) ? Age, in Gigayears real*4 clAge_1 (nr__1) ! (Gyr) ? Lower confidence limit on Age real*4 chAge_1 (nr__1) ! (Gyr) ? Upper confidence limit on Age real*4 mass (nr__1) ! (solMass) ? Mass, in solar masses real*4 clmass (nr__1) ! (solMass) ? Lower confidence limit on mass real*4 chmass (nr__1) ! (solMass) ? Upper confidence limit on mass real*4 RVel (nr__1) ! (km/s) ? Mean radial velocity (10) real*4 meRVel (nr__1) ! (km/s) ? Mean error of radial velocity real*4 e_RVel (nr__1) ! (km/s) ? Standard deviation of RVel measurements integer*4 o_RVel (nr__1) ! ? Number of RVel measurements integer*4 dT (nr__1) ! (d) ? Time-span of RVel measurements real*4 P_chi2 (nr__1) ! ? Probability of constant RVel character*1 fd (nr__1) ! Flag for spectroscopic binaries (10) character*1 fv (nr__1) ! Source of RVel (11) integer*4 vsini (nr__1) ! (km/s) ? Rotational velocity integer*4 pmRA (nr__1) ! (mas/a) ? Proper motion in right ascension (4) integer*4 pmDE (nr__1) ! (mas/a) ? Proper motion in declination (4) integer*4 e_pm (nr__1) ! (mas/a) ? Standard error of total proper motion (4) real*4 plx (nr__1) ! (mas) ? Hipparcos parallax real*4 e_plx (nr__1) ! (mas) ? Standard error of parallax integer*4 UVel_1 (nr__1) ! (km/s) ? Heliocentric space velocity component U (12) integer*4 VVel_1 (nr__1) ! (km/s) ? Heliocentric space velocity component V (12) integer*4 WVel_1 (nr__1) ! (km/s) ? Heliocentric space velocity component W (12) real*4 Rgal (nr__1) ! (kpc) ? Galactic radial position (13) real*4 zgal (nr__1) ! (kpc) ? Galactic vertical position (13) real*8 Rmin (nr__1) ! (kpc) ? Perigalactic distance real*8 Rmax (nr__1) ! (kpc) ? Apogalactic distance real*4 ecc (nr__1) ! ? Eccentricity of galactic orbit real*8 zmax (nr__1) ! (kpc) ? Maximum distance from galactic plane character*4 Notes (nr__1) ! General note (14) *Note (1): Object name given in order of preference. * A second number following a slash indicates that a double star with * separate HD numbers has been observed together. *Note (2): If the star is a member of a multiple system the component(s) * included in the photometry are identified here. *Note (3): This flag identifies confirmed and suspected binaries. * The information can come from one or several sources such as * photometry, radial velocity or astrometry. *Note (4): From the Tycho-2 catalogue. *Note (5): From the calibration of Olsen (1988A&A...189..173O) * if E(b-y)>=0.02 and Dist>40pc; otherwise the star is assumed * to be unreddened. *Note (6): From the calibration of Alonso et al. (1996A&A...313..873A) *Note (7): From the calibrations by Schuster & Nissen (1989A&A...221...65S), * Edvardsson et al. (1993A&A...275..101E), or the one defined in Sect * 4.3 of the present paper. *Note (8): Source for the distance: * H = Hipparcos parallax, * F,G = F or G-star photometric distance. * If a Hipparcos distance with relative error below 13% is available, * this is used, otherwise a photometric distance is used. *Note (9): Flag for suspected giants. * Indicates a disagreement between the photometric distance * determination and the Hipparcos parallax at the 3 sigma level, * suggesting that the star is a giant not detected from the photometry. *Note (10): Mean radial velocity. For double lined binaries, the * computed systemic velocity is given if so indicated by the fd flag. *Note (11): Source of the radial velocity: * C: Coravel, presented in the present paper * A: CfA from Nordstrom et al. (1997, Cat. <J/A+AS/126/21>) * L: literature from Barbier-Brossat & Figon (2000A&AS..142..217B, * Cat. <III/213>) *Note (12): Heliocentric space velocity components given in a right handed * coordinate system with U towards the galactic centre. *Note (13): With the solar position 8.000kpc from the galactic centre, * and 7pc above the plane. *Note (14): Global notes: * a = Double star with Delta_mag_<5mag * b = Double star with Hipparcos measurement of Delta_mag_<5mag * c = Variable star * d = Simbad note * e = Olsen (1983A&AS...54...55O) * f = Olsen (1979A&AS...37..367O) * g = Olsen (1980A&AS...39..205O) * h = Olsen (1993A&AS..102...89O) * i = Olsen (1994A&AS..104..429O) * j = Olsen (1994A&AS..106..257O) * k = Abt et al. (1979PASP...91..176A) * l = Abt (1984ApJ...285..247A) * m = Abt (1986ApJ...309..260A) * n = Gray & Garrison (1989ApJS...69..301G) * o = Gray (1989AJ.....89.1049G) * p = Henry et al. (1996AJ....111..439H) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Mass ratios for doubled-lined binaries (2004) integer*4 nr__2 parameter (nr__2=510) ! Number of records character*23 ar__2 ! Full-size record character*11 Name_2 (nr__2) ! HD number real*4 mratio (nr__2) ! Mass ratio (1) real*4 e_mratio (nr__2) ! ? Estimated uncertainty of the mass ratio (1) *Note (1): For spectroscopic binaries with multiple observations containing * double correlation peaks, the centre-of-mass velocity and the mass * ratio of the binary may be computed by the method of Wilson * (1941ApJ....93...29W). If only two observations are available, * no uncertainty can be estimated. C============================================================================= C Loading file 'newcat.dat' ! Main catalogue recalibrated data (2007) C Format for file interpretation 1 format( + I6,1X,A11,1X,A4,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + F5.3,1X,F5.2,1X,I3,1X,F5.2,1X,F4.1,1X,F4.1,1X,F4.1,1X,I4,1X, + I4,1X,I4) C Effective file loading open(unit=1,status='old',file= +'newcat.dat') write(6,*) '....Loading file: newcat.dat' do i__=1,16682 read(1,'(A96)')ar__ read(ar__,1) + HIP(i__),Name(i__),Comp(i__),RAh(i__),RAm(i__),RAs(i__), + DE_(i__),DEd(i__),DEm(i__),DEs(i__),logTe(i__),v_Fe_H_(i__), + Dist(i__),VMAG(i__),Age(i__),clAge(i__),chAge(i__),UVel(i__), + VVel(i__),WVel(i__) if(ar__(1:6) .EQ. '') HIP(i__) = iNULL__ if(ar__(46:50) .EQ. '') logTe(i__) = rNULL__ if(ar__(52:56) .EQ. '') v_Fe_H_(i__) = rNULL__ if(ar__(58:60) .EQ. '') Dist(i__) = iNULL__ if(ar__(62:66) .EQ. '') VMAG(i__) = rNULL__ if(ar__(68:71) .EQ. '') Age(i__) = rNULL__ if(ar__(73:76) .EQ. '') clAge(i__) = rNULL__ if(ar__(78:81) .EQ. '') chAge(i__) = rNULL__ if(ar__(83:86) .EQ. '') UVel(i__) = iNULL__ if(ar__(88:91) .EQ. '') VVel(i__) = iNULL__ if(ar__(93:96) .EQ. '') WVel(i__) = iNULL__ c ..............Just test output........... write(6,1) + HIP(i__),Name(i__),Comp(i__),RAh(i__),RAm(i__),RAs(i__), + DE_(i__),DEd(i__),DEm(i__),DEs(i__),logTe(i__),v_Fe_H_(i__), + Dist(i__),VMAG(i__),Age(i__),clAge(i__),chAge(i__),UVel(i__), + VVel(i__),WVel(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table1.dat' ! First version of the catalog (2004) C Format for file interpretation 2 format( + I6,1X,A11,1X,A4,1X,A1,1X,A1,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X, + I2,1X,I2,1X,I3,1X,I3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3, + 1X,F5.2,1X,I3,1X,F5.2,1X,F5.2,1X,A1,1X,A1,1X,F4.1,1X,F4.1,1X, + F4.1,1X,F4.2,1X,F4.2,1X,F4.2,13X,F6.1,1X,F4.1,1X,F5.1,1X,I3, + 1X,I4,1X,F5.3,1X,A1,1X,A1,1X,I3,1X,I5,1X,I5,1X,I2,1X,F5.1,1X, + F4.1,1X,I4,1X,I4,1X,I4,1X,F6.3,1X,F6.3,1X,F7.2,1X,F7.2,1X, + F5.2,1X,F7.2,1X,A4) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,16682 read(1,'(A287)')ar__1 read(ar__1,2) + HIP_1(i__),Name_1(i__),Comp_1(i__),fb(i__),fs(i__),RAh_1(i__), + RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__), + DEs_1(i__),GLON(i__),GLAT(i__),Vmag_1(i__),b_y(i__), + Hbeta(i__),E_b_y(i__),logTe_1(i__),v_Fe_H__1(i__),Dist_1(i__), + VMAG_2(i__),dVMag(i__),fr(i__),fg(i__),Age_1(i__), + clAge_1(i__),chAge_1(i__),mass(i__),clmass(i__),chmass(i__), + RVel(i__),meRVel(i__),e_RVel(i__),o_RVel(i__),dT(i__), + P_chi2(i__),fd(i__),fv(i__),vsini(i__),pmRA(i__),pmDE(i__), + e_pm(i__),plx(i__),e_plx(i__),UVel_1(i__),VVel_1(i__), + WVel_1(i__),Rgal(i__),zgal(i__),Rmin(i__),Rmax(i__),ecc(i__), + zmax(i__),Notes(i__) if(ar__1(1:6) .EQ. '') HIP_1(i__) = iNULL__ if(ar__1(71:75) .EQ. '') Hbeta(i__) = rNULL__ if(ar__1(77:82) .EQ. '') E_b_y(i__) = rNULL__ if(ar__1(84:88) .EQ. '') logTe_1(i__) = rNULL__ if(ar__1(90:94) .EQ. '') v_Fe_H__1(i__) = rNULL__ if(ar__1(96:98) .EQ. '') Dist_1(i__) = iNULL__ if(ar__1(100:104) .EQ. '') VMAG_2(i__) = rNULL__ if(ar__1(106:110) .EQ. '') dVMag(i__) = rNULL__ if(ar__1(116:119) .EQ. '') Age_1(i__) = rNULL__ if(ar__1(121:124) .EQ. '') clAge_1(i__) = rNULL__ if(ar__1(126:129) .EQ. '') chAge_1(i__) = rNULL__ if(ar__1(131:134) .EQ. '') mass(i__) = rNULL__ if(ar__1(136:139) .EQ. '') clmass(i__) = rNULL__ if(ar__1(141:144) .EQ. '') chmass(i__) = rNULL__ if(ar__1(158:163) .EQ. '') RVel(i__) = rNULL__ if(ar__1(165:168) .EQ. '') meRVel(i__) = rNULL__ if(ar__1(170:174) .EQ. '') e_RVel(i__) = rNULL__ if(ar__1(176:178) .EQ. '') o_RVel(i__) = iNULL__ if(ar__1(180:183) .EQ. '') dT(i__) = iNULL__ if(ar__1(185:189) .EQ. '') P_chi2(i__) = rNULL__ if(ar__1(195:197) .EQ. '') vsini(i__) = iNULL__ if(ar__1(199:203) .EQ. '') pmRA(i__) = iNULL__ if(ar__1(205:209) .EQ. '') pmDE(i__) = iNULL__ if(ar__1(211:212) .EQ. '') e_pm(i__) = iNULL__ if(ar__1(214:218) .EQ. '') plx(i__) = rNULL__ if(ar__1(220:223) .EQ. '') e_plx(i__) = rNULL__ if(ar__1(225:228) .EQ. '') UVel_1(i__) = iNULL__ if(ar__1(230:233) .EQ. '') VVel_1(i__) = iNULL__ if(ar__1(235:238) .EQ. '') WVel_1(i__) = iNULL__ if(ar__1(240:245) .EQ. '') Rgal(i__) = rNULL__ if(ar__1(247:252) .EQ. '') zgal(i__) = rNULL__ if(ar__1(254:260) .EQ. '') Rmin(i__) = rNULL__ if(ar__1(262:268) .EQ. '') Rmax(i__) = rNULL__ if(ar__1(270:274) .EQ. '') ecc(i__) = rNULL__ if(ar__1(276:282) .EQ. '') zmax(i__) = rNULL__ RAdeg(i__) = rNULL__ DEdeg(i__) = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh_1(i__) .GT. -180) RAdeg(i__)=RAh_1(i__)*15. if(RAm_1(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm_1(i__)/4. if(RAs_1(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs_1(i__)/240. if(DEd_1(i__) .GE. 0) DEdeg(i__)=DEd_1(i__) if(DEm_1(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm_1(i__)/60. if(DEs_1(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs_1(i__)/3600. if(DE__1(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,2) + HIP_1(i__),Name_1(i__),Comp_1(i__),fb(i__),fs(i__),RAh_1(i__), + RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__), + DEs_1(i__),GLON(i__),GLAT(i__),Vmag_1(i__),b_y(i__), + Hbeta(i__),E_b_y(i__),logTe_1(i__),v_Fe_H__1(i__),Dist_1(i__), + VMAG_2(i__),dVMag(i__),fr(i__),fg(i__),Age_1(i__), + clAge_1(i__),chAge_1(i__),mass(i__),clmass(i__),chmass(i__), + RVel(i__),meRVel(i__),e_RVel(i__),o_RVel(i__),dT(i__), + P_chi2(i__),fd(i__),fv(i__),vsini(i__),pmRA(i__),pmDE(i__), + e_pm(i__),plx(i__),e_plx(i__),UVel_1(i__),VVel_1(i__), + WVel_1(i__),Rgal(i__),zgal(i__),Rmin(i__),Rmax(i__),ecc(i__), + zmax(i__),Notes(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Mass ratios for doubled-lined binaries (2004) C Format for file interpretation 3 format(A11,1X,F5.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,510 read(1,'(A23)')ar__2 read(ar__2,3)Name_2(i__),mratio(i__),e_mratio(i__) if(ar__2(19:23) .EQ. '') e_mratio(i__) = rNULL__ c ..............Just test output........... write(6,3)Name_2(i__),mratio(i__),e_mratio(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end