FORTRAN Generation
(/./ftp/cats/J/A_A/625/A115)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/625/A115 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-Jun-08
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/625/A115            Ruprecht 147 DANCe. I.               (Olivares+, 2019)
*================================================================================
*Ruprecht 147 DANCe.
*I. Members, empirical isochrone, luminosity, and mass distributions.
*    Olivares J., Bouy H., Sarro L.M., Miret-Roig N., Berihuete A., Bertin E.,
*    Barrado D., Huelamo N., Tamura M., Allen L., Beletsky Y., Serre S.,
*    Cuillandre J.-C.
*    <Astron. Astrophys. 625, A115 (2019)>
*    =2019A&A...625A.115O        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table6.dat'	! Gaia DR2 ids and membership probabilities

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

      integer*8     GaiaDR2    (nr__) ! Gaia DR2 source id (source_id)
      integer*4     ID         (nr__) ! [1/16025095] Internal id (ID)
      real*4        Ppin0_6    (nr__) ! Cluster membership probability at pin 0.6
*                                   (probability_pin_0.6)
      real*4        Ppin0_7    (nr__) ! Cluster membership probability at pin 0.7
*                                   (probability_pin_0.7)
      real*4        Ppin0_8    (nr__) ! Cluster membership probability at pin 0.8
*                                   (probability_pin_0.8)
      real*4        Ppin0_9    (nr__) ! Cluster membership probability at pin 0.9
*                                   (probability_pin_0.9)

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

C  Declarations for 'table7.dat'	! DANCe astrometry, photometry and
                                 membership probabilities

      integer*4 nr__1
      parameter (nr__1=6776857)	! Number of records
      character*477 ar__1  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     ID_1       (nr__1) ! [1/16025095] Internal id (ID)
      real*8        RAdeg      (nr__1) ! (deg) Right ascension (J2000) (RAJ2000)
      real*8        e_RAdeg    (nr__1) ! (mas) Standard error of right ascension
*                                      (e_RAJ2000)
      real*8        DEdeg      (nr__1) ! (deg) Declination (J2000) (DEJ2000)
      real*8        e_DEdeg    (nr__1) ! (mas) Standard error of declination (e_DEJ2000)
      real*8        pmRA_      (nr__1) ! (mas/yr) Proper motion in right ascension (pmRA)
*                                      direction (pmRA*cosDE) (pmRA)
      real*4        e_pmRA_    (nr__1) ! (mas/yr) Standard error of pmRA (e_pmRA)
      real*4        e_pmDE     (nr__1) ! (mas/yr) Standard error of proper motion in
*                                      declination direction  (e_pmDec)
      real*8        pmDE       (nr__1) ! (mas/yr) Proper motion in declination direction
*                                      (pmDec)
      real*8        Plx        (nr__1) ! (mas) ? Absolute stellar parallax (parallax)
      real*4        e_Plx      (nr__1) ! (mas) ? Standard error of parallax
*                                     (parallax_error)
      real*8        PlxpmDEcor (nr__1) ! ? Correlation between parallax and proper
*                                      motion in declination
*                                      (parallax_pmdec_corr)
      real*8        RApmDEcor  (nr__1) ! ? Correlation between right ascension and
*                                      proper motion in declination
*                                      (ra_pmdec_corr)
      real*8        pmRApmDEcor(nr__1) ! ? Correlation between proper motion in
*                                      right ascension and proper motion in
*                                      declination (pmra_pmdec_corr)
      real*8        PlxpmRAcor (nr__1) ! ? Correlation between parallax and proper
*                                      motion in right ascension
*                                      (parallax_pmra_corr)
      real*8        DEPlxcor   (nr__1) ! ? Correlation between parallax and proper
*                                      motion in declination (dec_parallax_corr)
      real*8        DEpmDEcor  (nr__1) ! ? Correlation between declination and
*                                      proper motion in declination
*                                      (dec_pmdec_corr)
      real*8        DEpmRAcor  (nr__1) ! ? Correlation between declination and
*                                      proper motion in right ascension
*                                      (dec_pmra_corr)
      real*8        RAPlxcor   (nr__1) ! ? Correlation between right ascension and
*                                      parallax (ra_parallax_corr)
      real*8        RApmRAcor  (nr__1) ! ? Correlation between right ascension and
*                                      proper motion in right ascension
*                                      (ra_pmra_corr)
      real*8        Gmag       (nr__1) ! (mag) Gaia DR2 G magnitude (G)
      real*8        e_Gmag     (nr__1) ! (mag) Standard error of G mean magnitude (e_G)
      real*8        BPmag      (nr__1) ! (mag) Gaia DR2 BP magnitude (BP)
      real*8        e_BPmag    (nr__1) ! (mag) Standard error of BP mean magnitude (e_BP)
      real*8        RPmag      (nr__1) ! (mag) Gaia DR2 RP magnitude (RP)
      real*8        e_RPmag    (nr__1) ! (mag) Standard error of RP mean magnitude (e_RP)
      real*4        umag       (nr__1) ! (mag) ?=99 Pan-STARRS PS1 u magnitude (u_sdss)
      real*4        e_umag     (nr__1) ! (mag) ?=99 Standard error of Pan-STARRS PS1 u
*                                      magnitude (e_u_sdss)
      real*4        gmag_1     (nr__1) ! (mag) ?=99 Pan-STARRS PS1 g magnitude (g_sdss)
      real*4        e_gmag_1   (nr__1) ! (mag) ?=99 Standard error of Pan-STARRS PS1 g
*                                      magnitude (e_g_sdss)
      real*4        rmag       (nr__1) ! (mag) ?=99 Pan-STARRS PS1 r magnitude (r_sdss)
      real*4        e_rmag     (nr__1) ! (mag) ?=99 Standard error of Pan-STARRS PS1 r
*                                      magnitude (e_r_sdss)
      real*4        imag       (nr__1) ! (mag) ?=99 Pan-STARRS PS1 i magnitude (i_sdss)
      real*4        e_imag     (nr__1) ! (mag) ?=99 Standard error of Pan-STARRS PS1 i
*                                      magnitude (e_i_sdss)
      real*4        zmag       (nr__1) ! (mag) ?=99 Pan-STARRS PS1 z magnitude (z_sdss)
      real*4        e_zmag     (nr__1) ! (mag) ?=99 Standard error of Pan-STARRS PS1 z
*                                      magnitude (e_z_sdss)
      real*4        Ymag       (nr__1) ! (mag) ?=99 Pan-STARRS PS1 Y magnitude (Y)
      real*4        e_Ymag     (nr__1) ! (mag) ?=99 Standard error of Pan-STARRS PS1 Y
*                                      magnitude (e_Y)
      real*4        Jmag       (nr__1) ! (mag) ?=99 2MASS J magnitude (J)
      real*4        e_Jmag     (nr__1) ! (mag) ?=99 Standard error of J magnitude (e_J)
      real*4        Ksmag      (nr__1) ! (mag) ?=99 2MASS Ks magnitude (Ks)
      real*4        e_Hmag     (nr__1) ! (mag) ?=99 Standard error of H magnitude (e_H)
      real*4        Hmag       (nr__1) ! (mag) ?=99 2MASS H magnitude (H)
      real*4        e_Ksmag    (nr__1) ! (mag) ?=99 Standard error of 2MASS Ks 
*                                      magnitude (e_Ks)
      real*4        Ppin0_5    (nr__1) ! Cluster membership probability at pin 0.5
*                                      (probability_pin_0.5)
      real*4        Ppin0_6_1  (nr__1) ! Cluster membership probability at pin 0.6
*                                      (probability_pin_0.6)
      real*4        Ppin0_7_1  (nr__1) ! Cluster membership probability at pin 0.7
*                                      (probability_pin_0.7)
      real*4        Ppin0_8_1  (nr__1) ! Cluster membership probability at pin 0.8
*                                      (probability_pin_0.8)
      real*4        Ppin0_9_1  (nr__1) ! Cluster membership probability at pin 0.9
*                                      (probability_pin_0.9)

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

C  Declarations for 'table8.dat'	! DANCe Central id and membership probabilities

      integer*4 nr__2
      parameter (nr__2=2246660)	! Number of records
      character*122 ar__2  	! Full-size record

      integer*4     ID_2       (nr__2) ! Internal id (ID)
      real*4        Ppin0_5_1  (nr__2) ! Cluster membership probability at pin 0.5
*                                   (probability_pin_0.5)
      real*4        Ppin0_6_2  (nr__2) ! Cluster membership probability at pin 0.6
*                                   (probability_pin_0.6)
      real*4        Ppin0_7_2  (nr__2) ! Cluster membership probability at pin 0.7
*                                   (probability_pin_0.7)
      real*4        Ppin0_8_2  (nr__2) ! Cluster membership probability at pin 0.8
*                                   (probability_pin_0.8)
      real*4        Ppin0_9_2  (nr__2) ! Cluster membership probability at pin 0.9
*                                   (probability_pin_0.9)

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

C  Declarations for 'table9.dat'	! Cluster candidate members: id, astrometry,
                                  photometry, membership probability

      integer*4 nr__3
      parameter (nr__3=259)	! Number of records
      character*674 ar__3  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     IDmember   (nr__3) ! [1/259] ID member (ID_member)
      real*8        RAdeg_1    (nr__3) ! (deg) Right ascension (J2000) (RAJ2000)
      real*4        e_RAdeg_1  (nr__3) ! (mas) ? Standard error of right ascension
*                                      (e_RAJ2000)
      real*8        DEdeg_1    (nr__3) ! (deg) Declination (J2000) (DEJ2000)
      real*8        e_DEdeg_1  (nr__3) ! (mas) ? Standard error of declination (e_DEJ2000)
      integer*4     mGDR2      (nr__3) ! [0/1] Gaia DR2 member ? (member_GDR2)
      integer*4     mDANCe     (nr__3) ! [0/1] DANCe member ? (member_DANCe)
      integer*4     mDANCeC    (nr__3) ! [0/1] DANCe Central member ?
*                                      (member_DANCe_Central)
      integer*4     mNew       (nr__3) ! [0/1] New member ? (new_member)
      integer*4     Fregion    (nr__3) ! [0/1] In full region ? (Full_region)
      integer*4     Mregion    (nr__3) ! [0/1] In medium region ? (Middle_region)
      integer*4     Cregion    (nr__3) ! [0/1] In core region ? (Core_region)
      integer*4     WD         (nr__3) ! [0/1] White dwarf ? (WD)
      integer*4     EPIC       (nr__3) ! ? EPIC identification number (EPIC)
      integer*4     CWW        (nr__3) ! ? Curtis et al. (2013, Cat. J/AJ/145/134)
*                                      identification number (CWW)
      character*22  Com        (nr__3) ! Comments (Comments)
      real*4        pDANCeC    (nr__3) ! ? DANCe Central probability
*                                      (probability_DANCe_Central)
      integer*4     DANCe      (nr__3) ! ? DANCe identification number (ID_DANCe)
      real*4        pDANCe     (nr__3) ! ? DANCe probability (probability_DANCe)
      real*4        umag_1     (nr__3) ! (mag) ? Pan-STARR PS1 u magnitude (u_sdss)
      real*4        e_umag_1   (nr__3) ! (mag) ? Standard error of u magnitude (e_u_sdss)
      real*4        gmag_2     (nr__3) ! (mag) ? Pan-STARR PS1 g magnitude (g_sdss)
      real*4        e_gmag_2   (nr__3) ! (mag) ? Standard error of g magnitude (e_g_sdss)
      real*4        rmag_1     (nr__3) ! (mag) ? Pan-STARR PS1 r magnitude (r_sdss)
      real*4        e_rmag_1   (nr__3) ! (mag) ? Standard error of r magnitude (e_r_sdss)
      real*4        imag_1     (nr__3) ! (mag) ? Pan-STARR PS1 i magnitude (i_sdss)
      real*4        e_imag_1   (nr__3) ! (mag) ? Standard error of i magnitude (e_i_sdss)
      real*4        zmag_1     (nr__3) ! (mag) ? Pan-STARR PS1 z magnitude (z_sdss)
      real*4        e_zmag_1   (nr__3) ! (mag) ? Standard error of z magnitude (e_z_sdss)
      real*4        Ymag_1     (nr__3) ! (mag) ? Pan-STARR PS1 Y magnitude (Y)
      real*4        e_Ymag_1   (nr__3) ! (mag) ? Standard error of Y magnitude (e_Y)
      real*4        Jmag_1     (nr__3) ! (mag) ? 2MASS J magnitude (J)
      real*4        e_Jmag_1   (nr__3) ! (mag) ? Standard error of J magnitude (e_J)
      real*4        Hmag_1     (nr__3) ! (mag) ? 2MASS H magnitude (H)
      real*4        e_Hmag_1   (nr__3) ! (mag) ? Standard error of H magnitude (e_H)
      real*4        Ksmag_1    (nr__3) ! (mag) ? 2MASS Ks magnitude (Ks)
      real*4        e_Ksmag_1  (nr__3) ! (mag) ? Standard error of Ks magnitude (e_Ks)
      integer*8     GaiaDR2_1  (nr__3) ! ? Gaia DR2 source ID (source_id)
      real*8        RAGdeg     (nr__3) ! (deg) ? Gaia DR2 right ascension (ICRS)
*                                      at Ep=2015.5 (ra)
      real*4        e_RAGdeg   (nr__3) ! (mas) ? Gaia DR2 Standard error of
*                                      right ascension (e_RA*cosDE) (ra_error)
      real*8        DEGdeg     (nr__3) ! (deg) ? Gaia DR2 declination (ICRS)
*                                      at Ep=2015.5 (dec)
      real*4        e_DEGdeg   (nr__3) ! (mas) ? Gaia DR2 Standard error of declination
*                                      (dec_error)
      real*8        pmRA__1    (nr__3) ! (mas/yr) Gaia DR2 Proper motion in right ascension
*                                      direction (pmRA*cosDE) (pmra)
      real*8        e_pmRA__1  (nr__3) ! (mas/yr) Gaia DR2 standard error of proper motion in
*                                      right ascension direction (pmra_error)
      real*8        pmDE_1     (nr__3) ! (mas/yr) Gaia DR2 Proper motion in declination
*                                      direction (pmdec)
      real*8        e_pmDE_1   (nr__3) ! (mas/yr) Gaia DR2 standard error of proper motion in
*                                      declination (pmdec_error)
      real*4        Plx_1      (nr__3) ! (mas) ? Gaia DR2 absolute stellar parallax
*                                      (parallax)
      real*4        e_Plx_1    (nr__3) ! (mas) ? Gaia DR2 standard error of parallax
*                                      (parallax_error)
      real*8        PlxpmRAcor_1(nr__3) ! ? Gaia DR2 Correlation between parallax and
*                                      proper motion in right ascension
*                                      (parallax_pmra_corr)
      real*8        PlxpmDEcor_1(nr__3) ! ? Gaia DR2 Correlation between parallax and
*                                      proper motion in right ascension
*                                      (parallax_pmdec_corr)
      real*8        pmRApmDEcor_1(nr__3) ! ? Gaia DR2 Correlation between parallax and
*                                      proper motion in right ascension
*                                      (pmra_pmdec_corr)
      integer*4     Dup        (nr__3) ! [0/1]? Gaia DR2 Source with duplicate
*                                      sources (duplicated_source)
      real*8        Gmag_3     (nr__3) ! (mag) ? Gaia DR2 G-band mean magnitude (Vega)
*                                      (phot_g_mean_mag)
      real*4        e_Gmag_3   (nr__3) ! (mag) ? Gaia DR2 Standard error of G-band mean
*                                      magnitude (phot_g_mean_mag_error)
      real*8        BPmag_1    (nr__3) ! (mag) ? Gaia DR2 Integrated BP mean magnitude
*                                      (Vega) (phot_bp_mean_mag)
      real*4        e_BPmag_1  (nr__3) ! (mag) ? Gaia DR2 Standard error of BP-band mean
*                                      magnitude (phot_bp_mean_mag_error)
      real*8        RPmag_1    (nr__3) ! (mag) ? Gaia DR2 Integrated RP mean magnitude
*                                      (Vega) (phot_rp_mean_mag)
      real*4        e_RPmag_1  (nr__3) ! (mag) ? Gaia DR2 Standard error of RP-band mean
*                                      magnitude (phot_rp_mean_mag_error)
      real*4        RV         (nr__3) ! (km/s) ? Gaia DR2 radial velocity
*                                      (radial_velocity)
      real*4        e_RV       (nr__3) ! (km/s) ? Gaia DR2 Radial velocity error
*                                      (radial_velocity_error)
      real*8        epsi       (nr__3) ! ? Gaia DR2 Excess noise of the source
*                                      (astrometric_excess_noise)
      real*4        RFG        (nr__3) ! ? Gaia DR2 G-band mean flux divided by its
*                                      error  (phot_g_mean_flux_over_error)
      real*4        RFBP       (nr__3) ! ? Gaia DR2 Integrated BP mean flux divided
*                                      by its error
*                                      (phot_bp_mean_flux_over_error)
      real*4        RFRP       (nr__3) ! ? Gaia DR2 Integrated RP Integrated RP mean
*                                      flux divided by its error
*                                      (phot_rp_mean_flux_over_error)
      character*13  Var        (nr__3) ! [NOT_AVAILABLE] Gaia DR2 Photometric
*                                      variability flag (phot_variable_flag)
      real*4        AG         (nr__3) ! ? Gaia DR2 Estimate of extinction in the
*                                      G band from Apsis-Priam (a_g_val)
      real*8        RA2000     (nr__3) ! (deg) ? Gaia DR2 right ascension (ICRS) at
*                                      Ep=2000.0 (ra_epoch2000)
      real*8        e_RA2000   (nr__3) ! (mas) ? Gaia DR2 Standard error of right
*                                      ascension (e_RA*cosDE)
*                                      (ra_epoch2000_error)
      real*8        DE2000     (nr__3) ! (deg) ? Gaia DR2 declination (ICRS)
*                                      at Ep=2000.0 (dec_epoch2000)
      real*8        e_DE2000   (nr__3) ! (mas) ? Gaia DR2 Standard error of declination
*                                      (dec_epoch2000_error)
      integer*4     ID_3       (nr__3) ! ? [1/16025095] Internal id (ID_GDR2)
      real*4        pGDR2      (nr__3) ! ? Gaia DR2 probability (probability_GDR2)

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

C  Declarations for 'table10.dat'	! Cluster empirical sequence in Gaia photometry

      integer*4 nr__4
      parameter (nr__4=237)	! Number of records
      character*23 ar__4  	! Full-size record

      real*8        Gmag_4     (nr__4) ! (mag) Gaia DR2 G-band magnitude (G)
      real*8        BPmag_2    (nr__4) ! (mag) Gaia DR2 BP-band magnitude (BP)
      real*8        RPmag_2    (nr__4) ! (mag) Gaia DR2 RP-band magnitude (RP)

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

C  Declarations for 'table11.dat'	! Cluster empirical sequence in DANCe photometry

      integer*4 nr__5
      parameter (nr__5=137)	! Number of records
      character*63 ar__5  	! Full-size record

      real*8        gmag_5     (nr__5) ! (mag) Pan-STARR PS1 g magnitude (g)
      real*8        rmag_2     (nr__5) ! (mag) Pan-STARR PS1 r magnitude (r)
      real*8        imag_2     (nr__5) ! (mag) Pan-STARR PS1 i magnitude (i)
      real*8        zmag_2     (nr__5) ! (mag) Pan-STARR PS1 z magnitude (z)
      real*8        Ymag_2     (nr__5) ! (mag) Pan-STARR PS1 Y magnitude (Y)
      real*8        Jmag_2     (nr__5) ! (mag) 2MASS J magnitude (J)
      real*8        Hmag_2     (nr__5) ! (mag) 2MASS H magnitude (H)
      real*8        Ksmag_2    (nr__5) ! (mag) 2MASS Ks magnitude (Ks)

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

C  Declarations for 'table12.dat'	! Cluster mass distribution

      integer*4 nr__6
      parameter (nr__6=100)	! Number of records
      character*29 ar__6  	! Full-size record

      real*8        logMass    (nr__6) ! ([Msun]) Mass (log_mass)
      real*8        Dens       (nr__6) ! Density (dens)
      real*8        e_Dens     (nr__6) ! rms uncertainty on Dens (e_dens)

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

C  Loading file 'table6.dat'	! Gaia DR2 ids and membership probabilities

C  Format for file interpretation

    1 format(I19,1X,I8,1X,E22.15,1X,E22.15,1X,E22.15,1X,E22.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,12081778
        read(1,'(A120)')ar__
        read(ar__,1)
     +  GaiaDR2(i__),ID(i__),Ppin0_6(i__),Ppin0_7(i__),Ppin0_8(i__),
     +  Ppin0_9(i__)
c    ..............Just test output...........
        write(6,1)
     +  GaiaDR2(i__),ID(i__),Ppin0_6(i__),Ppin0_7(i__),Ppin0_8(i__),
     +  Ppin0_9(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! DANCe astrometry, photometry and
*                                 membership probabilities

C  Format for file interpretation

    2 format(
     +  I8,1X,F15.11,1X,F8.4,1X,F15.11,1X,F8.4,1X,F10.3,1X,F6.3,1X,
     +  F6.3,1X,F10.3,1X,F8.4,1X,F5.3,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F9.6,1X,
     +  F7.4,1X,F10.6,1X,F7.4,1X,F9.6,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,F6.3,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,F6.3,1X,
     +  E22.15,1X,E22.15,1X,E22.15,1X,E22.15,1X,E22.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,6776857
        read(1,'(A477)')ar__1
        read(ar__1,2)
     +  ID_1(i__),RAdeg(i__),e_RAdeg(i__),DEdeg(i__),e_DEdeg(i__),
     +  pmRA_(i__),e_pmRA_(i__),e_pmDE(i__),pmDE(i__),Plx(i__),
     +  e_Plx(i__),PlxpmDEcor(i__),RApmDEcor(i__),pmRApmDEcor(i__),
     +  PlxpmRAcor(i__),DEPlxcor(i__),DEpmDEcor(i__),DEpmRAcor(i__),
     +  RAPlxcor(i__),RApmRAcor(i__),Gmag(i__),e_Gmag(i__),BPmag(i__),
     +  e_BPmag(i__),RPmag(i__),e_RPmag(i__),umag(i__),e_umag(i__),
     +  gmag_1(i__),e_gmag_1(i__),rmag(i__),e_rmag(i__),imag(i__),
     +  e_imag(i__),zmag(i__),e_zmag(i__),Ymag(i__),e_Ymag(i__),
     +  Jmag(i__),e_Jmag(i__),Ksmag(i__),e_Hmag(i__),Hmag(i__),
     +  e_Ksmag(i__),Ppin0_5(i__),Ppin0_6_1(i__),Ppin0_7_1(i__),
     +  Ppin0_8_1(i__),Ppin0_9_1(i__)
        if(ar__1(96:103) .EQ. '') Plx(i__) = rNULL__
        if(ar__1(105:109) .EQ. '') e_Plx(i__) = rNULL__
        if(ar__1(111:117) .EQ. '') PlxpmDEcor(i__) = rNULL__
        if(ar__1(119:125) .EQ. '') RApmDEcor(i__) = rNULL__
        if(ar__1(127:133) .EQ. '') pmRApmDEcor(i__) = rNULL__
        if(ar__1(135:141) .EQ. '') PlxpmRAcor(i__) = rNULL__
        if(ar__1(143:149) .EQ. '') DEPlxcor(i__) = rNULL__
        if(ar__1(151:157) .EQ. '') DEpmDEcor(i__) = rNULL__
        if(ar__1(159:165) .EQ. '') DEpmRAcor(i__) = rNULL__
        if(ar__1(167:173) .EQ. '') RAPlxcor(i__) = rNULL__
        if(ar__1(175:181) .EQ. '') RApmRAcor(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  ID_1(i__),RAdeg(i__),e_RAdeg(i__),DEdeg(i__),e_DEdeg(i__),
     +  pmRA_(i__),e_pmRA_(i__),e_pmDE(i__),pmDE(i__),Plx(i__),
     +  e_Plx(i__),PlxpmDEcor(i__),RApmDEcor(i__),pmRApmDEcor(i__),
     +  PlxpmRAcor(i__),DEPlxcor(i__),DEpmDEcor(i__),DEpmRAcor(i__),
     +  RAPlxcor(i__),RApmRAcor(i__),Gmag(i__),e_Gmag(i__),BPmag(i__),
     +  e_BPmag(i__),RPmag(i__),e_RPmag(i__),umag(i__),e_umag(i__),
     +  gmag_1(i__),e_gmag_1(i__),rmag(i__),e_rmag(i__),imag(i__),
     +  e_imag(i__),zmag(i__),e_zmag(i__),Ymag(i__),e_Ymag(i__),
     +  Jmag(i__),e_Jmag(i__),Ksmag(i__),e_Hmag(i__),Hmag(i__),
     +  e_Ksmag(i__),Ppin0_5(i__),Ppin0_6_1(i__),Ppin0_7_1(i__),
     +  Ppin0_8_1(i__),Ppin0_9_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! DANCe Central id and membership probabilities

C  Format for file interpretation

    3 format(I7,1X,E22.15,1X,E22.15,1X,E22.15,1X,E22.15,1X,E22.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,2246660
        read(1,'(A122)')ar__2
        read(ar__2,3)
     +  ID_2(i__),Ppin0_5_1(i__),Ppin0_6_2(i__),Ppin0_7_2(i__),
     +  Ppin0_8_2(i__),Ppin0_9_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  ID_2(i__),Ppin0_5_1(i__),Ppin0_6_2(i__),Ppin0_7_2(i__),
     +  Ppin0_8_2(i__),Ppin0_9_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! Cluster candidate members: id, astrometry,
*                                  photometry, membership probability

C  Format for file interpretation

    4 format(
     +  I3,1X,F15.11,1X,F6.3,1X,F15.11,1X,F8.4,1X,I1,1X,I1,1X,I1,1X,
     +  I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I9,1X,I3,1X,A22,1X,E22.15,1X,I8,
     +  1X,E22.16,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,F6.3,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,F6.3,1X,I19,1X,F15.11,1X,F6.4,1X,
     +  F15.11,1X,F6.4,1X,F10.3,1X,F7.5,1X,F10.3,1X,F8.6,1X,F6.4,1X,
     +  F6.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,I1,1X,F7.4,1X,E12.6,1X,F7.4,
     +  1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.2,1X,F5.2,1X,F8.6,1X,E22.15,1X,
     +  E22.15,1X,E22.15,2X,A13,1X,F6.4,1X,F15.11,1X,F7.4,1X,F15.11,
     +  1X,F7.4,1X,I8,1X,E22.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,259
        read(1,'(A674)')ar__3
        read(ar__3,4)
     +  IDmember(i__),RAdeg_1(i__),e_RAdeg_1(i__),DEdeg_1(i__),
     +  e_DEdeg_1(i__),mGDR2(i__),mDANCe(i__),mDANCeC(i__),mNew(i__),
     +  Fregion(i__),Mregion(i__),Cregion(i__),WD(i__),EPIC(i__),
     +  CWW(i__),Com(i__),pDANCeC(i__),DANCe(i__),pDANCe(i__),
     +  umag_1(i__),e_umag_1(i__),gmag_2(i__),e_gmag_2(i__),
     +  rmag_1(i__),e_rmag_1(i__),imag_1(i__),e_imag_1(i__),
     +  zmag_1(i__),e_zmag_1(i__),Ymag_1(i__),e_Ymag_1(i__),
     +  Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),
     +  Ksmag_1(i__),e_Ksmag_1(i__),GaiaDR2_1(i__),RAGdeg(i__),
     +  e_RAGdeg(i__),DEGdeg(i__),e_DEGdeg(i__),pmRA__1(i__),
     +  e_pmRA__1(i__),pmDE_1(i__),e_pmDE_1(i__),Plx_1(i__),
     +  e_Plx_1(i__),PlxpmRAcor_1(i__),PlxpmDEcor_1(i__),
     +  pmRApmDEcor_1(i__),Dup(i__),Gmag_3(i__),e_Gmag_3(i__),
     +  BPmag_1(i__),e_BPmag_1(i__),RPmag_1(i__),e_RPmag_1(i__),
     +  RV(i__),e_RV(i__),epsi(i__),RFG(i__),RFBP(i__),RFRP(i__),
     +  Var(i__),AG(i__),RA2000(i__),e_RA2000(i__),DE2000(i__),
     +  e_DE2000(i__),ID_3(i__),pGDR2(i__)
        if(ar__3(21:26) .EQ. '') e_RAdeg_1(i__) = rNULL__
        if(ar__3(44:51) .EQ. '') e_DEdeg_1(i__) = rNULL__
        if(ar__3(69:77) .EQ. '') EPIC(i__) = iNULL__
        if(ar__3(79:81) .EQ. '') CWW(i__) = iNULL__
        if(ar__3(106:127) .EQ. '') pDANCeC(i__) = rNULL__
        if(ar__3(129:136) .EQ. '') DANCe(i__) = iNULL__
        if(ar__3(138:159) .EQ. '') pDANCe(i__) = rNULL__
        if(ar__3(161:166) .EQ. '') umag_1(i__) = rNULL__
        if(ar__3(168:173) .EQ. '') e_umag_1(i__) = rNULL__
        if(ar__3(175:180) .EQ. '') gmag_2(i__) = rNULL__
        if(ar__3(182:187) .EQ. '') e_gmag_2(i__) = rNULL__
        if(ar__3(189:194) .EQ. '') rmag_1(i__) = rNULL__
        if(ar__3(196:201) .EQ. '') e_rmag_1(i__) = rNULL__
        if(ar__3(203:208) .EQ. '') imag_1(i__) = rNULL__
        if(ar__3(210:215) .EQ. '') e_imag_1(i__) = rNULL__
        if(ar__3(217:222) .EQ. '') zmag_1(i__) = rNULL__
        if(ar__3(224:229) .EQ. '') e_zmag_1(i__) = rNULL__
        if(ar__3(231:236) .EQ. '') Ymag_1(i__) = rNULL__
        if(ar__3(238:243) .EQ. '') e_Ymag_1(i__) = rNULL__
        if(ar__3(245:250) .EQ. '') Jmag_1(i__) = rNULL__
        if(ar__3(252:257) .EQ. '') e_Jmag_1(i__) = rNULL__
        if(ar__3(259:264) .EQ. '') Hmag_1(i__) = rNULL__
        if(ar__3(266:271) .EQ. '') e_Hmag_1(i__) = rNULL__
        if(ar__3(273:278) .EQ. '') Ksmag_1(i__) = rNULL__
        if(ar__3(280:285) .EQ. '') e_Ksmag_1(i__) = rNULL__
        if(ar__3(287:305) .EQ. '') GaiaDR2_1(i__) = iNULL__
        if(ar__3(307:321) .EQ. '') RAGdeg(i__) = rNULL__
        if(ar__3(323:328) .EQ. '') e_RAGdeg(i__) = rNULL__
        if(ar__3(330:344) .EQ. '') DEGdeg(i__) = rNULL__
        if(ar__3(346:351) .EQ. '') e_DEGdeg(i__) = rNULL__
        if(ar__3(392:397) .EQ. '') Plx_1(i__) = rNULL__
        if(ar__3(399:404) .EQ. '') e_Plx_1(i__) = rNULL__
        if(ar__3(406:412) .EQ. '') PlxpmRAcor_1(i__) = rNULL__
        if(ar__3(414:420) .EQ. '') PlxpmDEcor_1(i__) = rNULL__
        if(ar__3(422:428) .EQ. '') pmRApmDEcor_1(i__) = rNULL__
        if(ar__3(430:430) .EQ. '') Dup(i__) = iNULL__
        if(ar__3(432:438) .EQ. '') Gmag_3(i__) = rNULL__
        if(ar__3(440:451) .EQ. '') e_Gmag_3(i__) = rNULL__
        if(ar__3(453:459) .EQ. '') BPmag_1(i__) = rNULL__
        if(ar__3(461:466) .EQ. '') e_BPmag_1(i__) = rNULL__
        if(ar__3(468:474) .EQ. '') RPmag_1(i__) = rNULL__
        if(ar__3(476:481) .EQ. '') e_RPmag_1(i__) = rNULL__
        if(ar__3(483:488) .EQ. '') RV(i__) = rNULL__
        if(ar__3(490:494) .EQ. '') e_RV(i__) = rNULL__
        if(ar__3(496:503) .EQ. '') epsi(i__) = rNULL__
        if(ar__3(505:526) .EQ. '') RFG(i__) = rNULL__
        if(ar__3(528:549) .EQ. '') RFBP(i__) = rNULL__
        if(ar__3(551:572) .EQ. '') RFRP(i__) = rNULL__
        if(ar__3(589:594) .EQ. '') AG(i__) = rNULL__
        if(ar__3(596:610) .EQ. '') RA2000(i__) = rNULL__
        if(ar__3(612:618) .EQ. '') e_RA2000(i__) = rNULL__
        if(ar__3(620:634) .EQ. '') DE2000(i__) = rNULL__
        if(ar__3(636:642) .EQ. '') e_DE2000(i__) = rNULL__
        if(ar__3(644:651) .EQ. '') ID_3(i__) = iNULL__
        if(ar__3(653:674) .EQ. '') pGDR2(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  IDmember(i__),RAdeg_1(i__),e_RAdeg_1(i__),DEdeg_1(i__),
     +  e_DEdeg_1(i__),mGDR2(i__),mDANCe(i__),mDANCeC(i__),mNew(i__),
     +  Fregion(i__),Mregion(i__),Cregion(i__),WD(i__),EPIC(i__),
     +  CWW(i__),Com(i__),pDANCeC(i__),DANCe(i__),pDANCe(i__),
     +  umag_1(i__),e_umag_1(i__),gmag_2(i__),e_gmag_2(i__),
     +  rmag_1(i__),e_rmag_1(i__),imag_1(i__),e_imag_1(i__),
     +  zmag_1(i__),e_zmag_1(i__),Ymag_1(i__),e_Ymag_1(i__),
     +  Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),
     +  Ksmag_1(i__),e_Ksmag_1(i__),GaiaDR2_1(i__),RAGdeg(i__),
     +  e_RAGdeg(i__),DEGdeg(i__),e_DEGdeg(i__),pmRA__1(i__),
     +  e_pmRA__1(i__),pmDE_1(i__),e_pmDE_1(i__),Plx_1(i__),
     +  e_Plx_1(i__),PlxpmRAcor_1(i__),PlxpmDEcor_1(i__),
     +  pmRApmDEcor_1(i__),Dup(i__),Gmag_3(i__),e_Gmag_3(i__),
     +  BPmag_1(i__),e_BPmag_1(i__),RPmag_1(i__),e_RPmag_1(i__),
     +  RV(i__),e_RV(i__),epsi(i__),RFG(i__),RFBP(i__),RFRP(i__),
     +  Var(i__),AG(i__),RA2000(i__),e_RA2000(i__),DE2000(i__),
     +  e_DE2000(i__),ID_3(i__),pGDR2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! Cluster empirical sequence in Gaia photometry

C  Format for file interpretation

    5 format(F7.4,1X,F7.4,1X,F7.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,237
        read(1,'(A23)')ar__4
        read(ar__4,5)Gmag_4(i__),BPmag_2(i__),RPmag_2(i__)
c    ..............Just test output...........
        write(6,5)Gmag_4(i__),BPmag_2(i__),RPmag_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table11.dat'	! Cluster empirical sequence in DANCe photometry

C  Format for file interpretation

    6 format(
     +  F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table11.dat')
      write(6,*) '....Loading file: table11.dat'
      do i__=1,137
        read(1,'(A63)')ar__5
        read(ar__5,6)
     +  gmag_5(i__),rmag_2(i__),imag_2(i__),zmag_2(i__),Ymag_2(i__),
     +  Jmag_2(i__),Hmag_2(i__),Ksmag_2(i__)
c    ..............Just test output...........
        write(6,6)
     +  gmag_5(i__),rmag_2(i__),imag_2(i__),zmag_2(i__),Ymag_2(i__),
     +  Jmag_2(i__),Hmag_2(i__),Ksmag_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table12.dat'	! Cluster mass distribution

C  Format for file interpretation

    7 format(F9.6,1X,F9.6,1X,F9.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table12.dat')
      write(6,*) '....Loading file: table12.dat'
      do i__=1,100
        read(1,'(A29)')ar__6
        read(ar__6,7)logMass(i__),Dens(i__),e_Dens(i__)
c    ..............Just test output...........
        write(6,7)logMass(i__),Dens(i__),e_Dens(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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