Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/498/949 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-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/498/949 Velocities of 1253 F-K dwarfs in open clusters (Mermilliod+, 2009) *================================================================================ *Radial- and rotational velocities for 1253 F-K dwarfs in 13 nearby open *clusters. * Mermilliod J.-C., Mayor M., Udry S. * <Astron. Astrophys. 498, 949 (2009)> * =2009A&A...498..949M C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Orbital elements integer*4 nr__ parameter (nr__=61) ! Number of records character*182 ar__ ! Full-size record character*9 Cluster (nr__) ! Cluster designation integer*4 WEBDA (nr__) ! Star designation (WEBDA) real*8 P0 (nr__) ! (d) Orbital period real*8 T0 (nr__) ! (d) T at periastron - 240000. real*8 gamma (nr__) ! (km/s) Systemic radial velocity real*4 e (nr__) ! Eccentricity real*4 Omega (nr__) ! (deg) ? Angle omega real*4 K1 (nr__) ! (km/s) Amplitude (A) real*4 K2 (nr__) ! (km/s) ? Amplitude (B) real*8 f_m (nr__) ! (solMass) Mass function real*8 asini (nr__) ! (Gm) semi major axis real*4 O_C (nr__) ! (km/s) sigma (observed - calculated) integer*4 Nobs (nr__) ! Number of observations real*8 e_P0 (nr__) ! (d) rms uncertainty on P0 real*8 e_T0 (nr__) ! (d) rms uncertainty on T0 real*4 e_gamma (nr__) ! (km/s) rms uncertainty on gamma real*4 e_e (nr__) ! ? rms uncertainty on eccentricity real*4 e_Omega (nr__) ! (deg) ? rms uncertainty on omega real*4 e_K1 (nr__) ! (km/s) rms uncertainty on K1 real*4 e_K2 (nr__) ! (km/s) ? rms uncertainty on K2 real*8 e_f_m (nr__) ! (solMass) ? rms uncertainty on f(m) real*8 e_asini (nr__) ! (Gm) rms uncertainty on asini real*4 e_O_C (nr__) ! (km/s) ? rms uncertainty on O-C c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table11.dat' ! Mean radial and rotational velocities integer*4 nr__1 parameter (nr__1=1259) ! Number of records character*123 ar__1 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__1) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__1) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*9 Cluster_1 (nr__1) ! Cluster designation integer*4 WEBDA_1 (nr__1) ! Star designation (WEBDA) character*11 OName (nr__1) ! Alternate designation integer*4 RAh (nr__1) ! (h) Right ascension (J2000) integer*4 RAm (nr__1) ! (min) Right ascension (J2000) real*4 RAs (nr__1) ! (s) Right ascension (J2000) character*1 DE_ (nr__1) ! Declination sign (J2000) integer*4 DEd (nr__1) ! (deg) Declination (J2000) integer*4 DEm (nr__1) ! (arcmin) Declination (J2000) real*4 DEs (nr__1) ! (arcsec) Declination (J2000) real*8 RV (nr__1) ! (km/s) Mean radial velocity real*4 e_RV (nr__1) ! (km/s) Error on the radial velocity real*4 sigma (nr__1) ! (km/s) RMS on velocity real*4 E_I (nr__1) ! ? External/Internal error ratio integer*4 N (nr__1) ! Number of measurements real*4 Vsini (nr__1) ! (km/s) ? Rotational velocity real*4 e_Vsini (nr__1) ! (km/s) ? Error on Vsini integer*4 DeltaT (nr__1) ! (d) ? Time interval spent by the observations real*4 Pchi2 (nr__1) ! ? Probability of variability P(chi2) character*8 Notes (nr__1) ! Membership, or duplicity c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table10.dat' ! Individual radial velocities integer*4 nr__2 parameter (nr__2=8149) ! Number of records character*80 ar__2 ! Full-size record character*9 Cluster_2 (nr__2) ! Cluster designation integer*4 WEBDA_2 (nr__2) ! Star designation (WEBDA) real*8 HJD (nr__2) ! Heliocentric Julian date character*1 Comp (nr__2) ! [ABCVS] Component in SB2 real*8 RV_1 (nr__2) ! (km/s) Radial velocity real*4 e_RV_1 (nr__2) ! (km/s) Error on the radial velocity character*10 Obs_Date (nr__2) ! ("DD/MM/YYYY") Date of observation integer*4 Integ (nr__2) ! (s) ? Integration length integer*4 Counts (nr__2) ! (ct) ? Number of counts character*3 Obs (nr__2) ! Observatory C============================================================================= C Loading file 'table8.dat' ! Orbital elements C Format for file interpretation 1 format( + A9,3X,I4,2X,F11.6,2X,F9.3,1X,F7.2,2X,F5.3,2X,F6.2,2X,F5.2,2X, + F5.2,2X,F8.6,2X,F8.4,3X,F4.2,3X,I2,1X,F10.6,1X,F7.3,1X,F4.2, + 1X,F5.3,1X,F5.2,1X,F4.2,1X,F4.2,2X,F8.6,1X,F8.5,2X,F6.4) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,61 read(1,'(A182)')ar__ read(ar__,1) + Cluster(i__),WEBDA(i__),P0(i__),T0(i__),gamma(i__),e(i__), + Omega(i__),K1(i__),K2(i__),f_m(i__),asini(i__),O_C(i__), + Nobs(i__),e_P0(i__),e_T0(i__),e_gamma(i__),e_e(i__), + e_Omega(i__),e_K1(i__),e_K2(i__),e_f_m(i__),e_asini(i__), + e_O_C(i__) if(ar__(58:63) .EQ. '') Omega(i__) = rNULL__ if(ar__(73:77) .EQ. '') K2(i__) = rNULL__ if(ar__(135:139) .EQ. '') e_e(i__) = rNULL__ if(ar__(141:145) .EQ. '') e_Omega(i__) = rNULL__ if(ar__(152:155) .EQ. '') e_K2(i__) = rNULL__ if(ar__(158:165) .EQ. '') e_f_m(i__) = rNULL__ if(ar__(177:182) .EQ. '') e_O_C(i__) = rNULL__ c ..............Just test output........... write(6,1) + Cluster(i__),WEBDA(i__),P0(i__),T0(i__),gamma(i__),e(i__), + Omega(i__),K1(i__),K2(i__),f_m(i__),asini(i__),O_C(i__), + Nobs(i__),e_P0(i__),e_T0(i__),e_gamma(i__),e_e(i__), + e_Omega(i__),e_K1(i__),e_K2(i__),e_f_m(i__),e_asini(i__), + e_O_C(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table11.dat' ! Mean radial and rotational velocities C Format for file interpretation 2 format( + A9,3X,I4,2X,A11,1X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,1X, + F7.2,2X,F5.2,2X,F5.2,2X,F5.2,3X,I2,2X,F4.1,2X,F4.1,2X,I4,2X, + F5.3,2X,A8) C Effective file loading open(unit=1,status='old',file= +'table11.dat') write(6,*) '....Loading file: table11.dat' do i__=1,1259 read(1,'(A123)')ar__1 read(ar__1,2) + Cluster_1(i__),WEBDA_1(i__),OName(i__),RAh(i__),RAm(i__), + RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),RV(i__), + e_RV(i__),sigma(i__),E_I(i__),N(i__),Vsini(i__),e_Vsini(i__), + DeltaT(i__),Pchi2(i__),Notes(i__) if(ar__1(79:83) .EQ. '') E_I(i__) = rNULL__ if(ar__1(91:94) .EQ. '') Vsini(i__) = rNULL__ if(ar__1(97:100) .EQ. '') e_Vsini(i__) = rNULL__ if(ar__1(103:106) .EQ. '') DeltaT(i__) = iNULL__ if(ar__1(109:113) .EQ. '') Pchi2(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(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15. if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4. if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240. if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__) if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60. if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,2) + Cluster_1(i__),WEBDA_1(i__),OName(i__),RAh(i__),RAm(i__), + RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),RV(i__), + e_RV(i__),sigma(i__),E_I(i__),N(i__),Vsini(i__),e_Vsini(i__), + DeltaT(i__),Pchi2(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 'table10.dat' ! Individual radial velocities C Format for file interpretation 3 format( + A9,3X,I4,2X,F11.3,2X,A1,2X,F7.2,2X,F4.2,3X,A10,2X,I4,3X,I6,2X, + A3) C Effective file loading open(unit=1,status='old',file= +'table10.dat') write(6,*) '....Loading file: table10.dat' do i__=1,8149 read(1,'(A80)')ar__2 read(ar__2,3) + Cluster_2(i__),WEBDA_2(i__),HJD(i__),Comp(i__),RV_1(i__), + e_RV_1(i__),Obs_Date(i__),Integ(i__),Counts(i__),Obs(i__) if(ar__2(63:66) .EQ. '') Integ(i__) = iNULL__ if(ar__2(70:75) .EQ. '') Counts(i__) = iNULL__ c ..............Just test output........... write(6,3) + Cluster_2(i__),WEBDA_2(i__),HJD(i__),Comp(i__),RV_1(i__), + e_RV_1(i__),Obs_Date(i__),Integ(i__),Counts(i__),Obs(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end