FORTRAN Generation
(/./ftp/cats/J/A_A/615/A49)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/615/A49 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 2026-Aug-11
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/615/A49  TGAS stars membership in 128 open clusters (Cantat-Gaudin+, 2018)
*================================================================================
*Characterising open clusters in the solar neigbourhood with the
*Tycho-Gaia astrometric solution.
*    Cantat-Gaudin T., Vallenari A., Sordo R., Pensabene F., Krone-Martins A.,
*    Moitinho A., Jordi C., Casamiquela L., Balaguer-Nunez L., Soubiran C.,
*    Brouillet N.
*    <Astron. Astrophys. 615, A49 (2018)>
*    =2018A&A...615A..49C        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'members.dat'	! Membership probability for TGAS stars in
                                 128 open clusters

      integer*4 nr__
      parameter (nr__=241711)	! Number of records
      character*596 ar__   	! Full-size record

C  J2015.0 position composed of: RAdeg DEdeg
      real*8        RAdeg      (nr__) ! (deg) GaiaDR1 Right ascension (J2015.0) (ra)
      real*4        DEdeg      (nr__) ! (deg) GaiaDR1 Declination (J2015.0) (dec)
      real*8        e_RAdeg    (nr__) ! (mas) Uncertainty on Gaia Right ascension
*                                      (ra_error)
      real*8        e_DEdeg    (nr__) ! (mas) Uncertainty on Gaia Declination (dec_error)
      integer*4     HIP        (nr__) ! ? Hipparcos ID (hip)
      character*12  TYC2       (nr__) ! Tycho-2 ID (tycho2_id)
      integer*8     Source     (nr__) ! GaiaDR1 Source ID (source_id)
      real*4        plx        (nr__) ! (mas) GaiaDR1 parallax (parallax)
      real*8        e_plx      (nr__) ! (mas) GaiaDR1 parallax error (parallax_error)
      real*4        pmRA       (nr__) ! (mas/yr) GaiaDR1 proper motion in RA (pmra)
      real*8        e_pmRA     (nr__) ! (mas/yr) GaiaDR1 pmra error (pmra_error)
      real*4        pmDE       (nr__) ! (mas/yr) GaiaDR1 proper motion in DE (pmdec)
      real*8        e_pmDE     (nr__) ! (mas/yr) GaiaDR1 pmdec error (pmdec_error)
      integer*4     NobsG      (nr__) ! GaiaDR1 number of G observations (Gnobs)
      real*8        Gmag       (nr__) ! (mag) GaiaDR1 G-band magnitude (G)
      integer*4     v2MASSID   (nr__) ! 2MASS ID (tmass_oid)
      real*8        RA2deg     (nr__) ! (deg) 2MASS Right ascension (J2000) (ra_twomass)
      real*4        DE2deg     (nr__) ! (deg) 2MASS Declination (J2000) (dec_twomass)
      character*9   UCAC4      (nr__) ! UCAC4 ID (ucac4_oid)
      real*8        RAUdeg     (nr__) ! (deg) UCAC4 Right ascension (J2000) (ra_ucac)
      real*4        DEUdeg     (nr__) ! (deg) UCAC4 Declination (J2000) (dec_ucac)
      real*4        pmRAU      (nr__) ! (mas/yr) ? UCAC4 pmra (pmra_ucac)
      real*4        e_pmRAU    (nr__) ! (mas/yr) ? UCAC4 pmra_error (pmra_error_ucac)
      real*4        pmDEU      (nr__) ! (mas/yr) ? UCAC4 pmdec (pmdec_ucac)
      real*4        e_pmDEU    (nr__) ! (mas/yr) ? UCAC4 pmdec_error (pmdec_error_ucac)
      real*4        Jmag       (nr__) ! (mag) ? 2MASS J magnitude (1.2um) (Jmag)
      real*4        Hmag       (nr__) ! (mag) ? 2MASS H magnitude (1.6um) (Hmag)
      real*4        Kmag       (nr__) ! (mag) ? 2MASS Ks magnitude (2.2um) (Kmag)
      real*4        Bmag       (nr__) ! (mag) ? B magnitude from APASS (Bmag)
      real*4        Vmag       (nr__) ! (mag) ? V magnitude from APASS (Vmag)
      real*4        gmag_1     (nr__) ! (mag) ? g magnitude from APASS (gmag)
      real*4        rmag       (nr__) ! (mag) ? r magnitude from APASS (rmag)
      real*4        imag       (nr__) ! (mag) ? i magnitude from APASS (imag)
      real*4        RADEcor    (nr__) ! Correlation between GaiaDR1 ra and dec
*                                      (ra_dec_corr)
      real*4        RAplxcor   (nr__) ! Correlation between GaiaDR1 ra and parallax
*                                      (ra_parallax_corr)
      real*4        RApmRAcor  (nr__) ! Correlation between GaiaDR1 ra and pmra
*                                      (ra_pmra_corr)
      real*4        RApmDEcor  (nr__) ! Correlation between GaiaDR1 ra and pmdec
*                                      (ra_pmdec_corr)
      real*4        DEplxcor   (nr__) ! Correlation between GaiaDR1 dec and
*                                      parallax (dec_parallax_corr)
      real*4        DEpmRAcor  (nr__) ! Correlation between GaiaDR1 dec and pmra
*                                      (dec_pmra_corr)
      real*4        DEpmDEcor  (nr__) ! Correlation between GaiaDR1 dec and pmdec
*                                      (dec_pmdec_corr)
      real*4        plxpmRAcor (nr__) ! Correlation between GaiaDR1 parallax and
*                                      pmra (parallax_pmra_corr)
      real*4        plxpmDEcor (nr__) ! Correlation between GaiaDR1 parallax and
*                                      pmdec (parallax_pmdec_corr)
      real*4        pmRApmDEcor(nr__) ! Correlation between GaiaDR1 pmra and pmdec
*                                      (pmra_pmdec_corr)
      character*13  Cluster    (nr__) ! Open cluster investigated (ocname) (1)
      character*1   PMused     (nr__) ! [TU] Proper motion catalogue used to
*                                      determine membership (PMused) (2)
      real*4        Pfinal     (nr__) ! Probability that the star belongs to this
*                                      open cluster (finalProba)
*Note (1): a star might be listed multiple times with different ocnames and
*     membership probabilities if it is present in the investigated fields of
*     neighbouring clusters.
*Note (2): T stands for TGAS and U stands for UCAC4

C=============================================================================

C  Loading file 'members.dat'	! Membership probability for TGAS stars in
*                                 128 open clusters

C  Format for file interpretation

    1 format(
     +  F19.15,1X,E21.15,1X,F21.18,1X,F20.17,1X,I6,1X,A12,1X,I19,1X,
     +  E22.18,1X,F19.17,1X,E22.18,1X,F21.18,1X,E22.18,1X,F21.18,1X,
     +  I4,1X,F19.16,1X,I9,1X,F10.6,1X,E10.6,1X,A9,1X,F13.9,1X,E13.9,
     +  1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,E13.9,1X,E14.9,1X,
     +  E13.9,1X,E13.9,1X,E13.9,1X,E13.9,1X,E14.9,1X,E14.9,1X,E13.9,
     +  1X,E13.9,1X,A13,1X,A1,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'members.dat')
      write(6,*) '....Loading file: members.dat'
      do i__=1,241711
        read(1,'(A596)')ar__
        read(ar__,1)
     +  RAdeg(i__),DEdeg(i__),e_RAdeg(i__),e_DEdeg(i__),HIP(i__),
     +  TYC2(i__),Source(i__),plx(i__),e_plx(i__),pmRA(i__),
     +  e_pmRA(i__),pmDE(i__),e_pmDE(i__),NobsG(i__),Gmag(i__),
     +  v2MASSID(i__),RA2deg(i__),DE2deg(i__),UCAC4(i__),RAUdeg(i__),
     +  DEUdeg(i__),pmRAU(i__),e_pmRAU(i__),pmDEU(i__),e_pmDEU(i__),
     +  Jmag(i__),Hmag(i__),Kmag(i__),Bmag(i__),Vmag(i__),gmag_1(i__),
     +  rmag(i__),imag(i__),RADEcor(i__),RAplxcor(i__),RApmRAcor(i__),
     +  RApmDEcor(i__),DEplxcor(i__),DEpmRAcor(i__),DEpmDEcor(i__),
     +  plxpmRAcor(i__),plxpmDEcor(i__),pmRApmDEcor(i__),Cluster(i__),
     +  PMused(i__),Pfinal(i__)
        if(ar__(86:91) .EQ. '') HIP(i__) = iNULL__
        if(ar__(354:359) .EQ. '') pmRAU(i__) = rNULL__
        if(ar__(361:364) .EQ. '') e_pmRAU(i__) = rNULL__
        if(ar__(366:371) .EQ. '') pmDEU(i__) = rNULL__
        if(ar__(373:376) .EQ. '') e_pmDEU(i__) = rNULL__
        if(ar__(378:383) .EQ. '') Jmag(i__) = rNULL__
        if(ar__(385:390) .EQ. '') Hmag(i__) = rNULL__
        if(ar__(392:397) .EQ. '') Kmag(i__) = rNULL__
        if(ar__(399:404) .EQ. '') Bmag(i__) = rNULL__
        if(ar__(406:411) .EQ. '') Vmag(i__) = rNULL__
        if(ar__(413:418) .EQ. '') gmag_1(i__) = rNULL__
        if(ar__(420:425) .EQ. '') rmag(i__) = rNULL__
        if(ar__(427:432) .EQ. '') imag(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  RAdeg(i__),DEdeg(i__),e_RAdeg(i__),e_DEdeg(i__),HIP(i__),
     +  TYC2(i__),Source(i__),plx(i__),e_plx(i__),pmRA(i__),
     +  e_pmRA(i__),pmDE(i__),e_pmDE(i__),NobsG(i__),Gmag(i__),
     +  v2MASSID(i__),RA2deg(i__),DE2deg(i__),UCAC4(i__),RAUdeg(i__),
     +  DEUdeg(i__),pmRAU(i__),e_pmRAU(i__),pmDEU(i__),e_pmDEU(i__),
     +  Jmag(i__),Hmag(i__),Kmag(i__),Bmag(i__),Vmag(i__),gmag_1(i__),
     +  rmag(i__),imag(i__),RADEcor(i__),RAplxcor(i__),RApmRAcor(i__),
     +  RApmDEcor(i__),DEplxcor(i__),DEpmRAcor(i__),DEpmDEcor(i__),
     +  plxpmRAcor(i__),plxpmDEcor(i__),pmRApmDEcor(i__),Cluster(i__),
     +  PMused(i__),Pfinal(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end