FORTRAN Generation
(/./ftp/cats/VI/137)

Conversion of standardized ReadMe file for file /./ftp/cats/VI/137 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 2018-May-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. VI/137            GaiaSimu Universe Model Snapshot               (Robin+, 2012)
*================================================================================
*GaiaSimu Universe Model Snapshot.
*   Robin A.C., Luri X., Reyle C., Isasi, Grux E., Blanco S., Arenou F.,
*   Babusiaux C., Belcheva M., Drimmel R., Jordi C., Krone-Martins A.,
*   Masana E., Mauduit J.C., Mignard F., Mowlawi N., Rocca-Volmerange B.,
*   Sartoretti P., Slezak E., Sozzetti A.
*  <Astron. Astrophys. 543, A100 (2012)> 
*  =2012A&A...543A.100R
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'gum_mw.sam'	! Gaia Universe Model Snapshot (GUMS): Milky Way stars
                            (among 2,143,475,885 stars)

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

C  Position composed of: RAdeg DEdeg (Epoch=J2010)
      character*19  ID         (nr__) ! source identifier (sourceExtendedId)
      real*8        RA2kdeg    (nr__) ! (deg) Right ascension ICRS, Epoch=J2000 (G1)
      real*8        DE2kdeg    (nr__) ! (deg) Declination ICRS, Epoch=J2000 (G1)
      real*4        VMAG       (nr__) ! (mag) V-band absolute magnitude (meanAbsoluteV)
      real*4        Gmag       (nr__) ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag      (nr__) ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag      (nr__) ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag      (nr__) ! (mag) Gaia Grvs magnitude (860+/-28nm) (magGRvs)
      real*8        RAdeg      (nr__) ! (deg) Right ascension (ICRS, Epoch=J2010) (alpha)
      real*8        DEdeg      (nr__) ! (deg) Declination (ICRS, Epoch=J2010) (delta)
      real*8        r          (nr__) ! (pc) Barycentric distance (distance)
      real*8        pmRA       (nr__) ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*8        pmDE       (nr__) ! (mas/yr) Proper motion in Declination (muDelta)
      real*8        RV         (nr__) ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        V_I        (nr__) ! (mag) Intrinsic color V-I (colorVminusI)
      real*4        Av         (nr__) ! (mag) Absorption (Av)
      real*4        Age        (nr__) ! (Gyr) Estimated age (age)
      real*4        v_a_Fe_    (nr__) ! ([Sun]) {alpha}-elements enhancement (alphaFe)
      real*4        bAlb       (nr__) ! [0/1] Bond albedo (for exoplanet) (bondAlbedo)
      real*4        e          (nr__) ! Eccentricity (for binaries) (eccentricity)
      real*4        v_Fe_H_    (nr__) ! ([Sun]) Metallicity [Fe/H] (feH)
      integer*4     fI         (nr__) ! [0,1] 1 (true) if interacting with a companion
*                                 (flagInteracting)
      real*4        gAlb       (nr__) ! Geometric albedo (for exoplanet) (geomAlbedo)
      integer*4     fm         (nr__) ! [0,1] 1 (true) if photocenter moves
*                                 (hasPhotocenterMotion)
      integer*4     Host       (nr__) ! Host galaxy: 1=Milky Way, 2=SMC, 3=LMC (host)
      real*4        i          (nr__) ! (deg) Orbit Inclination (inclination)
      real*4        logg       (nr__) ! ([cm/s2]) Gravity (log) (logg)
      real*4        Omega      (nr__) ! (deg) Longitude of ascending node {Omega}
*                                 (longitudeAscendingNode)
      real*8        Mass       (nr__) ! (Sun) Stellar mass (mass)
      real*4        Mbol       (nr__) ! (mag) Bolometric absolute magnitude (mbol)
      integer*4     nc         (nr__) ! [1-5] Number of components in the system
      integer*4     nt         (nr__) ! [1-5] Total number of objects in the system
      real*4        P          (nr__) ! (d) Orbital period (orbitPeriod)
      real*4        omega_1    (nr__) ! (deg) Periastron argument {omega}(periastronArgument)
      real*4        T0         (nr__) ! (d) Periastron date (MJD) (periastronDate)
      real*4        Phase      (nr__) ! [0/1] Phase if an exoplanet (phase)
      integer*4     Pop        (nr__) ! [1/3]?=99 Stellar population the object
*                                        belongs to (population)
      real*8        BeEnv      (nr__) ! Envelope characteristic for Be star (Renv/R*)
*                                 (rEnvRStar)
      real*4        Radius     (nr__) ! (Rsun) Stellar radius (radius)
      real*4        a          (nr__) ! (au) Semi-major axis (in binary system)
*                                 (semimajorAxis)
      integer*4     Teff       (nr__) ! (K) Effective temperature (teff)
      real*4        Vamp       (nr__) ! (mag) Variability amplitude (variabilityAmplitude)
      real*8        Vper       (nr__) ! (d) Variability period (variabilityPeriod)
      real*4        Vphase     (nr__) ! [0/1] Variability phase (variabilityPhase)
      character*4   Vtype      (nr__) ! Variability type (variabilityType) (2)
      real*4        vsini      (nr__) ! (km/s) Rotational velocity
      integer*8     recno      (nr__) ! [1/2143475885] Original record number
*Note (2): variability types are:
*   Ceph = Cepheids (pulsating supergiants)
*   delt = {delta} Sct-type (fainter and shorter periods that Cepheids)
*   RRLy = RR Lyr-type, sub-classified in RRab (asymmetric) or RRc (symmetric)
*   gamm = {gamma} Dor-type (non-radial pulsations)
*   RoAp = rapidly oscillating Ap stars
*   ZZce = pulsating white dwarfs
*   Mira = Mira Ceti type (Cool red supergiants)
*   semi = semiregular, usually red giants or supergiants
*   ACV  = {alpha}^2^ CVn-type (strong magnetic field)
*   Flar = flare stars (magnetic reconnection in dM stars)
*   Nova = classical nova (white dwarf accreting matter from its companion)
*   DN   = dwarf novae (white dwarf accreting matter with instable disk)
*   EB   = eclipsing binary
*   micr = microlensing event

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

C  Declarations for 'gum_lmc.sam'	! Gaia Universe Model Snapshot (GUMS): LMC stars
                            (among 7,559,826 stars)

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

C  Position composed of: RAdeg DEdeg (Epoch=J2010)
      character*19  ID_1       (nr__1) ! source identifier (sourceExtendedId)
      real*8        RA2kdeg_1  (nr__1) ! (deg) Right ascension ICRS, Epoch=J2000 (G1)
      real*8        DE2kdeg_1  (nr__1) ! (deg) Declination ICRS, Epoch=J2000 (G1)
      real*4        VMAG_1     (nr__1) ! (mag) V-band absolute magnitude (meanAbsoluteV)
      real*4        Gmag_1     (nr__1) ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag_1    (nr__1) ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag_1    (nr__1) ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag_1    (nr__1) ! (mag) Gaia Grvs magnitude (860+/-28nm) (magGRvs)
      real*8        RAdeg_1    (nr__1) ! (deg) Right ascension (ICRS, Epoch=J2010) (alpha)
      real*8        DEdeg_1    (nr__1) ! (deg) Declination (ICRS, Epoch=J2010) (delta)
      real*8        r_1        (nr__1) ! (pc) Barycentric distance (distance)
      real*8        pmRA_1     (nr__1) ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*8        pmDE_1     (nr__1) ! (mas/yr) Proper motion in Declination (muDelta)
      real*8        RV_1       (nr__1) ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        V_I_1      (nr__1) ! (mag) Intrinsic color V-I (colorVminusI)
      real*4        Av_1       (nr__1) ! (mag) Absorption (Av)
      real*4        Age_1      (nr__1) ! (Gyr) Estimated age (age)
      real*4        v_a_Fe__1  (nr__1) ! ([Sun]) {alpha}-elements enhancement (alphaFe)
      real*4        bAlb_1     (nr__1) ! [0/1] Bond albedo (for exoplanet) (bondAlbedo)
      real*4        e_1        (nr__1) ! Eccentricity (for binaries) (eccentricity)
      real*4        v_Fe_H__1  (nr__1) ! ([Sun]) Metallicity [Fe/H] (feH)
      integer*4     fI_1       (nr__1) ! [0,1] 1 (true) if interacting with a companion
*                                 (flagInteracting)
      real*4        gAlb_1     (nr__1) ! Geometric albedo (for exoplanet) (geomAlbedo)
      integer*4     fm_1       (nr__1) ! [0,1] 1 (true) if photocenter moves
*                                 (hasPhotocenterMotion)
      integer*4     Host_1     (nr__1) ! Host galaxy: 1=Milky Way, 2=SMC, 3=LMC (host)
      real*4        i_1        (nr__1) ! (deg) Orbit Inclination (inclination)
      real*4        logg_1     (nr__1) ! ([cm/s2]) Gravity (log) (logg)
      real*4        Omega_2    (nr__1) ! (deg) Longitude of ascending node {Omega}
*                                 (longitudeAscendingNode)
      real*8        Mass_1     (nr__1) ! (Sun) Stellar mass (mass)
      real*4        Mbol_1     (nr__1) ! (mag) Bolometric absolute magnitude (mbol)
      integer*4     nc_1       (nr__1) ! [1-5] Number of components in the system
      integer*4     nt_1       (nr__1) ! [1-5] Total number of objects in the system
      real*4        P_1        (nr__1) ! (d) Orbital period (orbitPeriod)
      real*4        omega_3    (nr__1) ! (deg) Periastron argument {omega}(periastronArgument)
      real*4        T0_1       (nr__1) ! (d) Periastron date (MJD) (periastronDate)
      real*4        Phase_1    (nr__1) ! [0/1] Phase if an exoplanet (phase)
      integer*4     Pop_1      (nr__1) ! [1/3]?=99 Stellar population the object
*                                        belongs to (population)
      real*8        BeEnv_1    (nr__1) ! Envelope characteristic for Be star (Renv/R*)
*                                 (rEnvRStar)
      real*4        Radius_1   (nr__1) ! (Rsun) Stellar radius (radius)
      real*4        a_1        (nr__1) ! (au) Semi-major axis (in binary system)
*                                 (semimajorAxis)
      integer*4     Teff_1     (nr__1) ! (K) Effective temperature (teff)
      real*4        Vamp_1     (nr__1) ! (mag) Variability amplitude (variabilityAmplitude)
      real*8        Vper_1     (nr__1) ! (d) Variability period (variabilityPeriod)
      real*4        Vphase_1   (nr__1) ! [0/1] Variability phase (variabilityPhase)
      character*4   Vtype_1    (nr__1) ! Variability type (variabilityType) (2)
      real*4        vsini_1    (nr__1) ! (km/s) Rotational velocity
      integer*8     recno_1    (nr__1) ! [1/2143475885] Original record number
*Note (2): variability types are:
*   Ceph = Cepheids (pulsating supergiants)
*   delt = {delta} Sct-type (fainter and shorter periods that Cepheids)
*   RRLy = RR Lyr-type, sub-classified in RRab (asymmetric) or RRc (symmetric)
*   gamm = {gamma} Dor-type (non-radial pulsations)
*   RoAp = rapidly oscillating Ap stars
*   ZZce = pulsating white dwarfs
*   Mira = Mira Ceti type (Cool red supergiants)
*   semi = semiregular, usually red giants or supergiants
*   ACV  = {alpha}^2^ CVn-type (strong magnetic field)
*   Flar = flare stars (magnetic reconnection in dM stars)
*   Nova = classical nova (white dwarf accreting matter from its companion)
*   DN   = dwarf novae (white dwarf accreting matter with instable disk)
*   EB   = eclipsing binary
*   micr = microlensing event

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

C  Declarations for 'gum_smc.sam'	! Gaia Universe Model Snapshot (GUMS): SMC stars
                            (among 1,250,384 stars)

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

C  Position composed of: RAdeg DEdeg (Epoch=J2010)
      character*19  ID_2       (nr__2) ! source identifier (sourceExtendedId)
      real*8        RA2kdeg_2  (nr__2) ! (deg) Right ascension ICRS, Epoch=J2000 (G1)
      real*8        DE2kdeg_2  (nr__2) ! (deg) Declination ICRS, Epoch=J2000 (G1)
      real*4        VMAG_2     (nr__2) ! (mag) V-band absolute magnitude (meanAbsoluteV)
      real*4        Gmag_2     (nr__2) ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag_2    (nr__2) ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag_2    (nr__2) ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag_2    (nr__2) ! (mag) Gaia Grvs magnitude (860+/-28nm) (magGRvs)
      real*8        RAdeg_2    (nr__2) ! (deg) Right ascension (ICRS, Epoch=J2010) (alpha)
      real*8        DEdeg_2    (nr__2) ! (deg) Declination (ICRS, Epoch=J2010) (delta)
      real*8        r_2        (nr__2) ! (pc) Barycentric distance (distance)
      real*8        pmRA_2     (nr__2) ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*8        pmDE_2     (nr__2) ! (mas/yr) Proper motion in Declination (muDelta)
      real*8        RV_2       (nr__2) ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        V_I_2      (nr__2) ! (mag) Intrinsic color V-I (colorVminusI)
      real*4        Av_2       (nr__2) ! (mag) Absorption (Av)
      real*4        Age_2      (nr__2) ! (Gyr) Estimated age (age)
      real*4        v_a_Fe__2  (nr__2) ! ([Sun]) {alpha}-elements enhancement (alphaFe)
      real*4        bAlb_2     (nr__2) ! [0/1] Bond albedo (for exoplanet) (bondAlbedo)
      real*4        e_2        (nr__2) ! Eccentricity (for binaries) (eccentricity)
      real*4        v_Fe_H__2  (nr__2) ! ([Sun]) Metallicity [Fe/H] (feH)
      integer*4     fI_2       (nr__2) ! [0,1] 1 (true) if interacting with a companion
*                                 (flagInteracting)
      real*4        gAlb_2     (nr__2) ! Geometric albedo (for exoplanet) (geomAlbedo)
      integer*4     fm_2       (nr__2) ! [0,1] 1 (true) if photocenter moves
*                                 (hasPhotocenterMotion)
      integer*4     Host_2     (nr__2) ! Host galaxy: 1=Milky Way, 2=SMC, 3=LMC (host)
      real*4        i_2        (nr__2) ! (deg) Orbit Inclination (inclination)
      real*4        logg_2     (nr__2) ! ([cm/s2]) Gravity (log) (logg)
      real*4        Omega_4    (nr__2) ! (deg) Longitude of ascending node {Omega}
*                                 (longitudeAscendingNode)
      real*8        Mass_2     (nr__2) ! (Sun) Stellar mass (mass)
      real*4        Mbol_2     (nr__2) ! (mag) Bolometric absolute magnitude (mbol)
      integer*4     nc_2       (nr__2) ! [1-5] Number of components in the system
      integer*4     nt_2       (nr__2) ! [1-5] Total number of objects in the system
      real*4        P_2        (nr__2) ! (d) Orbital period (orbitPeriod)
      real*4        omega_5    (nr__2) ! (deg) Periastron argument {omega}(periastronArgument)
      real*4        T0_2       (nr__2) ! (d) Periastron date (MJD) (periastronDate)
      real*4        Phase_2    (nr__2) ! [0/1] Phase if an exoplanet (phase)
      integer*4     Pop_2      (nr__2) ! [1/3]?=99 Stellar population the object
*                                        belongs to (population)
      real*8        BeEnv_2    (nr__2) ! Envelope characteristic for Be star (Renv/R*)
*                                 (rEnvRStar)
      real*4        Radius_2   (nr__2) ! (Rsun) Stellar radius (radius)
      real*4        a_2        (nr__2) ! (au) Semi-major axis (in binary system)
*                                 (semimajorAxis)
      integer*4     Teff_2     (nr__2) ! (K) Effective temperature (teff)
      real*4        Vamp_2     (nr__2) ! (mag) Variability amplitude (variabilityAmplitude)
      real*8        Vper_2     (nr__2) ! (d) Variability period (variabilityPeriod)
      real*4        Vphase_2   (nr__2) ! [0/1] Variability phase (variabilityPhase)
      character*4   Vtype_2    (nr__2) ! Variability type (variabilityType) (2)
      real*4        vsini_2    (nr__2) ! (km/s) Rotational velocity
      integer*8     recno_2    (nr__2) ! [1/2143475885] Original record number
*Note (2): variability types are:
*   Ceph = Cepheids (pulsating supergiants)
*   delt = {delta} Sct-type (fainter and shorter periods that Cepheids)
*   RRLy = RR Lyr-type, sub-classified in RRab (asymmetric) or RRc (symmetric)
*   gamm = {gamma} Dor-type (non-radial pulsations)
*   RoAp = rapidly oscillating Ap stars
*   ZZce = pulsating white dwarfs
*   Mira = Mira Ceti type (Cool red supergiants)
*   semi = semiregular, usually red giants or supergiants
*   ACV  = {alpha}^2^ CVn-type (strong magnetic field)
*   Flar = flare stars (magnetic reconnection in dM stars)
*   Nova = classical nova (white dwarf accreting matter from its companion)
*   DN   = dwarf novae (white dwarf accreting matter with instable disk)
*   EB   = eclipsing binary
*   micr = microlensing event

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

C  Declarations for 'gum_gal.sam'	! Gaia Universe Model Snapshot (GUMS) (galaxies)
                            (among 37,831,197 sources)

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

C  Position composed of: RAdeg DEdeg (Epoch=J2010)
      character*19  ID_3       (nr__3) ! Source identifier (sourceExtendedId)
      real*8        RA2kdeg_3  (nr__3) ! (deg) Right ascension ICRS, Epoch=J2000 (G1)
      real*8        DE2kdeg_3  (nr__3) ! (deg) Declination ICRS, Epoch=J2000 (G1)
      real*8        Vmag_3     (nr__3) ! (mag) V-band intrinsic (extinction free) apparent
*                                  mean magnitude (meanAbsoluteV)
      real*4        Gmag_3     (nr__3) ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag_3    (nr__3) ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag_3    (nr__3) ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag_3    (nr__3) ! (mag) Gaia Grvs magnitude (860+/-28nm) (magGRvs)nm)
      real*8        RAdeg_3    (nr__3) ! (deg) Right ascension (ICRS, Epoch=J2010) (alpha)
      real*8        DEdeg_3    (nr__3) ! (deg) Declination (ICRS, Epoch=J2010) (delta)
      real*4        r_3        (nr__3) ! (pc) Barycentric distance (distance)
      real*4        pmRA_3     (nr__3) ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*4        pmDE_3     (nr__3) ! (mas/yr) Proper motion in Declination (muDelta)
      real*8        RV_3       (nr__3) ! (km/s) Barycentric radial velocity (radialVelocity)
      integer*4     SNID       (nr__3) ! [-1] Attached supernova ID (not filled)
      real*4        V_I_3      (nr__3) ! (mag) Intrinsic color V-I (colorVminusI)
      integer*4     Htype      (nr__3) ! [0-8] Hubble type (hubbleType) (3)
      real*4        Av_3       (nr__3) ! (mag) Absorption in the V-band
      real*8        z          (nr__3) ! Redshift (redShift)
      real*4        B_T        (nr__3) ! [0,1] Bulge to total flux ratio (4)
      real*4        Brad       (nr__3) ! (arcsec) Bulge equivalent radius (4)
      real*4        Bb_a       (nr__3) ! [0,1] Bulge aspect ratio b/a (4)
      real*8        Bpa        (nr__3) ! (deg) [-180/180] Bulge position angle (4)
      real*8        Dsize      (nr__3) ! (arcsec) Disk scale length (4)
      real*4        Db_a       (nr__3) ! [0,1] Disk aspect ratio b/a (4)
      real*8        Dpa        (nr__3) ! (deg) Disk position angle (4)
      integer*4     recno_3    (nr__3) ! [1/37831197] Original record number
*Note (3): the 8 classes are:
*   0..7 = E2, E-S0, Sa, Sb, Sbc, Sc, Sd, Im
*      8 = QSPG: Quenched Star Forming Galaxies (2009A&A...504.1071T)
*Note (4): these parameters represent the SkyMaker shape parameters.
*     Note that the position angles are measured with respect to the RA axis.

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

C  Declarations for 'gum_qso.dat'	! Gaia Universe Model Snapshot (GUMS) (quasars)

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

C  Position composed of: RAdeg DEdeg (Epoch=J2010)
      real*8        RA2kdeg_4  (nr__4) ! (deg) Right ascension (ICRS), at Epoch J2000 (G1)
      real*8        DE2kdeg_4  (nr__4) ! (deg) Declination (ICRS), at Epoch J2000 (G1)
      real*8        RAdeg_4    (nr__4) ! (deg) Right ascension (ICRS), at Epoch J2010 (alpha)
      real*8        DEdeg_4    (nr__4) ! (deg) Right ascension (ICRS), at Epoch J2010 (delta)
      integer*4     r_4        (nr__4) ! (Mpc) Barycentric distance (distance)
      real*8        pmRA_4     (nr__4) ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*8        pmDE_4     (nr__4) ! (mas/yr) Proper motion in Declination (muDelta)
      integer*4     RV_4       (nr__4) ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        Gmag_4     (nr__4) ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag_4    (nr__4) ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag_4    (nr__4) ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag_4    (nr__4) ! (mag) Gaia/Grvs band magnitude (850-880nm) (magGRvs)
      real*4        Av_4       (nr__4) ! (mag) Galactic absorption in the V-band
      real*4        z_1        (nr__4) ! Redshift (redShift)
      real*4        slope      (nr__4) ! Slope of the continuum
      integer*4     W          (nr__4) ! Total equivalent width of the emission lines
      character*18  ID_4       (nr__4) ! Source identifier (sourceExtendedId)
      character*19  objID      (nr__4) ! GUMS source identifier (sourceId)

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

C  Declarations for 'gum_sn.dat'	! Gaia Universe Model Snapshot (GUMS (supernovae)

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

C  Position composed of: RAdeg DEdeg (Epoch=J2010)
      real*8        RA2kdeg_5  (nr__5) ! (deg) Right ascension (ICRS), at Epoch J2000 (G1)
      real*8        DE2kdeg_5  (nr__5) ! (deg) Declination (ICRS), at Epoch J2000 (G1)
      real*8        RAdeg_5    (nr__5) ! (deg) Right ascension (ICRS), at Epoch J2010 (alpha)
      real*8        DEdeg_5    (nr__5) ! (deg) Right ascension (ICRS), at Epoch J2010 (delta)
      real*8        r_5        (nr__5) ! (Mpc) Barycentric distance distance)
      real*4        pmRA_5     (nr__5) ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*4        pmDE_5     (nr__5) ! (mas/yr) Proper motion in Declination (muDelta)
      integer*4     RV_5       (nr__5) ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        VMAG_4     (nr__5) ! (mag) V-band absolute magnitude (meanAbsoluteV)
      real*8        z_2        (nr__5) ! Redshift (redShift)
      character*4   SNt        (nr__5) ! Type of the supernova (gmdbtype)
      character*18  ID_5       (nr__5) ! Source identifier (sourceExtendedId)
      character*19  objID_1    (nr__5) ! GUMS extended source identifier (sourceId)

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

C  Loading file 'gum_mw.sam'	! Gaia Universe Model Snapshot (GUMS): Milky Way stars
*                            (among 2,143,475,885 stars)

C  Format for file interpretation

    1 format(
     +  A19,1X,F14.10,1X,F14.10,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F14.10,1X,F14.10,1X,F7.1,1X,F9.3,1X,F9.3,1X,F7.1,1X,
     +  F6.3,1X,F5.3,1X,F6.3,1X,F5.2,1X,F4.2,1X,F5.3,1X,F5.2,1X,I1,1X,
     +  F4.2,1X,I1,1X,I1,1X,F4.1,1X,F6.3,1X,F5.1,1X,F9.4,1X,F6.3,1X,
     +  I1,1X,I1,1X,E11.6,1X,F5.1,1X,E11.5,1X,F3.1,1X,I2,1X,F7.4,1X,
     +  E11.4,1X,E11.4,1X,I6,1X,F6.3,1X,F9.5,1X,F5.3,1X,A4,1X,F5.1,1X,
     +  I10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gum_mw.sam')
      write(6,*) '....Loading file: gum_mw.sam'
      do i__=1,1000
        read(1,'(A364)')ar__
        read(ar__,1)
     +  ID(i__),RA2kdeg(i__),DE2kdeg(i__),VMAG(i__),Gmag(i__),
     +  GBmag(i__),GRmag(i__),Gsmag(i__),RAdeg(i__),DEdeg(i__),r(i__),
     +  pmRA(i__),pmDE(i__),RV(i__),V_I(i__),Av(i__),Age(i__),
     +  v_a_Fe_(i__),bAlb(i__),e(i__),v_Fe_H_(i__),fI(i__),gAlb(i__),
     +  fm(i__),Host(i__),i(i__),logg(i__),Omega(i__),Mass(i__),
     +  Mbol(i__),nc(i__),nt(i__),P(i__),omega_1(i__),T0(i__),
     +  Phase(i__),Pop(i__),BeEnv(i__),Radius(i__),a(i__),Teff(i__),
     +  Vamp(i__),Vper(i__),Vphase(i__),Vtype(i__),vsini(i__),
     +  recno(i__)
c    ..............Just test output...........
        write(6,1)
     +  ID(i__),RA2kdeg(i__),DE2kdeg(i__),VMAG(i__),Gmag(i__),
     +  GBmag(i__),GRmag(i__),Gsmag(i__),RAdeg(i__),DEdeg(i__),r(i__),
     +  pmRA(i__),pmDE(i__),RV(i__),V_I(i__),Av(i__),Age(i__),
     +  v_a_Fe_(i__),bAlb(i__),e(i__),v_Fe_H_(i__),fI(i__),gAlb(i__),
     +  fm(i__),Host(i__),i(i__),logg(i__),Omega(i__),Mass(i__),
     +  Mbol(i__),nc(i__),nt(i__),P(i__),omega_1(i__),T0(i__),
     +  Phase(i__),Pop(i__),BeEnv(i__),Radius(i__),a(i__),Teff(i__),
     +  Vamp(i__),Vper(i__),Vphase(i__),Vtype(i__),vsini(i__),
     +  recno(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'gum_lmc.sam'	! Gaia Universe Model Snapshot (GUMS): LMC stars
*                            (among 7,559,826 stars)

C  Format for file interpretation

    2 format(
     +  A19,1X,F14.10,1X,F14.10,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F14.10,1X,F14.10,1X,F7.1,1X,F9.3,1X,F9.3,1X,F7.1,1X,
     +  F6.3,1X,F5.3,1X,F6.3,1X,F5.2,1X,F4.2,1X,F5.3,1X,F5.2,1X,I1,1X,
     +  F4.2,1X,I1,1X,I1,1X,F4.1,1X,F6.3,1X,F5.1,1X,F9.4,1X,F6.3,1X,
     +  I1,1X,I1,1X,E11.6,1X,F5.1,1X,E11.5,1X,F3.1,1X,I2,1X,F7.4,1X,
     +  E11.4,1X,E11.4,1X,I6,1X,F6.3,1X,F9.5,1X,F5.3,1X,A4,1X,F5.1,1X,
     +  I10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gum_lmc.sam')
      write(6,*) '....Loading file: gum_lmc.sam'
      do i__=1,1000
        read(1,'(A364)')ar__1
        read(ar__1,2)
     +  ID_1(i__),RA2kdeg_1(i__),DE2kdeg_1(i__),VMAG_1(i__),
     +  Gmag_1(i__),GBmag_1(i__),GRmag_1(i__),Gsmag_1(i__),
     +  RAdeg_1(i__),DEdeg_1(i__),r_1(i__),pmRA_1(i__),pmDE_1(i__),
     +  RV_1(i__),V_I_1(i__),Av_1(i__),Age_1(i__),v_a_Fe__1(i__),
     +  bAlb_1(i__),e_1(i__),v_Fe_H__1(i__),fI_1(i__),gAlb_1(i__),
     +  fm_1(i__),Host_1(i__),i_1(i__),logg_1(i__),Omega_2(i__),
     +  Mass_1(i__),Mbol_1(i__),nc_1(i__),nt_1(i__),P_1(i__),
     +  omega_3(i__),T0_1(i__),Phase_1(i__),Pop_1(i__),BeEnv_1(i__),
     +  Radius_1(i__),a_1(i__),Teff_1(i__),Vamp_1(i__),Vper_1(i__),
     +  Vphase_1(i__),Vtype_1(i__),vsini_1(i__),recno_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  ID_1(i__),RA2kdeg_1(i__),DE2kdeg_1(i__),VMAG_1(i__),
     +  Gmag_1(i__),GBmag_1(i__),GRmag_1(i__),Gsmag_1(i__),
     +  RAdeg_1(i__),DEdeg_1(i__),r_1(i__),pmRA_1(i__),pmDE_1(i__),
     +  RV_1(i__),V_I_1(i__),Av_1(i__),Age_1(i__),v_a_Fe__1(i__),
     +  bAlb_1(i__),e_1(i__),v_Fe_H__1(i__),fI_1(i__),gAlb_1(i__),
     +  fm_1(i__),Host_1(i__),i_1(i__),logg_1(i__),Omega_2(i__),
     +  Mass_1(i__),Mbol_1(i__),nc_1(i__),nt_1(i__),P_1(i__),
     +  omega_3(i__),T0_1(i__),Phase_1(i__),Pop_1(i__),BeEnv_1(i__),
     +  Radius_1(i__),a_1(i__),Teff_1(i__),Vamp_1(i__),Vper_1(i__),
     +  Vphase_1(i__),Vtype_1(i__),vsini_1(i__),recno_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'gum_smc.sam'	! Gaia Universe Model Snapshot (GUMS): SMC stars
*                            (among 1,250,384 stars)

C  Format for file interpretation

    3 format(
     +  A19,1X,F14.10,1X,F14.10,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F14.10,1X,F14.10,1X,F7.1,1X,F9.3,1X,F9.3,1X,F7.1,1X,
     +  F6.3,1X,F5.3,1X,F6.3,1X,F5.2,1X,F4.2,1X,F5.3,1X,F5.2,1X,I1,1X,
     +  F4.2,1X,I1,1X,I1,1X,F4.1,1X,F6.3,1X,F5.1,1X,F9.4,1X,F6.3,1X,
     +  I1,1X,I1,1X,E11.6,1X,F5.1,1X,E11.5,1X,F3.1,1X,I2,1X,F7.4,1X,
     +  E11.4,1X,E11.4,1X,I6,1X,F6.3,1X,F9.5,1X,F5.3,1X,A4,1X,F5.1,1X,
     +  I10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gum_smc.sam')
      write(6,*) '....Loading file: gum_smc.sam'
      do i__=1,1000
        read(1,'(A364)')ar__2
        read(ar__2,3)
     +  ID_2(i__),RA2kdeg_2(i__),DE2kdeg_2(i__),VMAG_2(i__),
     +  Gmag_2(i__),GBmag_2(i__),GRmag_2(i__),Gsmag_2(i__),
     +  RAdeg_2(i__),DEdeg_2(i__),r_2(i__),pmRA_2(i__),pmDE_2(i__),
     +  RV_2(i__),V_I_2(i__),Av_2(i__),Age_2(i__),v_a_Fe__2(i__),
     +  bAlb_2(i__),e_2(i__),v_Fe_H__2(i__),fI_2(i__),gAlb_2(i__),
     +  fm_2(i__),Host_2(i__),i_2(i__),logg_2(i__),Omega_4(i__),
     +  Mass_2(i__),Mbol_2(i__),nc_2(i__),nt_2(i__),P_2(i__),
     +  omega_5(i__),T0_2(i__),Phase_2(i__),Pop_2(i__),BeEnv_2(i__),
     +  Radius_2(i__),a_2(i__),Teff_2(i__),Vamp_2(i__),Vper_2(i__),
     +  Vphase_2(i__),Vtype_2(i__),vsini_2(i__),recno_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  ID_2(i__),RA2kdeg_2(i__),DE2kdeg_2(i__),VMAG_2(i__),
     +  Gmag_2(i__),GBmag_2(i__),GRmag_2(i__),Gsmag_2(i__),
     +  RAdeg_2(i__),DEdeg_2(i__),r_2(i__),pmRA_2(i__),pmDE_2(i__),
     +  RV_2(i__),V_I_2(i__),Av_2(i__),Age_2(i__),v_a_Fe__2(i__),
     +  bAlb_2(i__),e_2(i__),v_Fe_H__2(i__),fI_2(i__),gAlb_2(i__),
     +  fm_2(i__),Host_2(i__),i_2(i__),logg_2(i__),Omega_4(i__),
     +  Mass_2(i__),Mbol_2(i__),nc_2(i__),nt_2(i__),P_2(i__),
     +  omega_5(i__),T0_2(i__),Phase_2(i__),Pop_2(i__),BeEnv_2(i__),
     +  Radius_2(i__),a_2(i__),Teff_2(i__),Vamp_2(i__),Vper_2(i__),
     +  Vphase_2(i__),Vtype_2(i__),vsini_2(i__),recno_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'gum_gal.sam'	! Gaia Universe Model Snapshot (GUMS) (galaxies)
*                            (among 37,831,197 sources)

C  Format for file interpretation

    4 format(
     +  A19,1X,F14.10,1X,F14.10,1X,F7.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F14.10,1X,F14.10,1X,E12.4,1X,F6.3,1X,F6.3,1X,F8.1,2X,
     +  I2,1X,F5.3,1X,I1,1X,F6.3,1X,F8.6,1X,F6.4,1X,F6.4,1X,F6.4,1X,
     +  F7.2,1X,F7.4,1X,F6.4,1X,F7.2,1X,I8)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gum_gal.sam')
      write(6,*) '....Loading file: gum_gal.sam'
      do i__=1,1000
        read(1,'(A240)')ar__3
        read(ar__3,4)
     +  ID_3(i__),RA2kdeg_3(i__),DE2kdeg_3(i__),Vmag_3(i__),
     +  Gmag_3(i__),GBmag_3(i__),GRmag_3(i__),Gsmag_3(i__),
     +  RAdeg_3(i__),DEdeg_3(i__),r_3(i__),pmRA_3(i__),pmDE_3(i__),
     +  RV_3(i__),SNID(i__),V_I_3(i__),Htype(i__),Av_3(i__),z(i__),
     +  B_T(i__),Brad(i__),Bb_a(i__),Bpa(i__),Dsize(i__),Db_a(i__),
     +  Dpa(i__),recno_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  ID_3(i__),RA2kdeg_3(i__),DE2kdeg_3(i__),Vmag_3(i__),
     +  Gmag_3(i__),GBmag_3(i__),GRmag_3(i__),Gsmag_3(i__),
     +  RAdeg_3(i__),DEdeg_3(i__),r_3(i__),pmRA_3(i__),pmDE_3(i__),
     +  RV_3(i__),SNID(i__),V_I_3(i__),Htype(i__),Av_3(i__),z(i__),
     +  B_T(i__),Brad(i__),Bb_a(i__),Bpa(i__),Dsize(i__),Db_a(i__),
     +  Dpa(i__),recno_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'gum_qso.dat'	! Gaia Universe Model Snapshot (GUMS) (quasars)

C  Format for file interpretation

    5 format(
     +  F12.8,1X,F12.8,1X,F12.8,1X,F12.8,2X,I5,1X,F9.5,1X,F9.5,2X,I6,
     +  1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,2X,F5.3,1X,F6.3,2X,I5,
     +  1X,A18,6X,A19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gum_qso.dat')
      write(6,*) '....Loading file: gum_qso.dat'
      do i__=1,979315
        read(1,'(A186)')ar__4
        read(ar__4,5)
     +  RA2kdeg_4(i__),DE2kdeg_4(i__),RAdeg_4(i__),DEdeg_4(i__),
     +  r_4(i__),pmRA_4(i__),pmDE_4(i__),RV_4(i__),Gmag_4(i__),
     +  GBmag_4(i__),GRmag_4(i__),Gsmag_4(i__),Av_4(i__),z_1(i__),
     +  slope(i__),W(i__),ID_4(i__),objID(i__)
c    ..............Just test output...........
        write(6,5)
     +  RA2kdeg_4(i__),DE2kdeg_4(i__),RAdeg_4(i__),DEdeg_4(i__),
     +  r_4(i__),pmRA_4(i__),pmDE_4(i__),RV_4(i__),Gmag_4(i__),
     +  GBmag_4(i__),GRmag_4(i__),Gsmag_4(i__),Av_4(i__),z_1(i__),
     +  slope(i__),W(i__),ID_4(i__),objID(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'gum_sn.dat'	! Gaia Universe Model Snapshot (GUMS (supernovae)

C  Format for file interpretation

    6 format(
     +  F14.10,1X,F14.10,1X,F14.10,1X,F14.10,1X,F7.2,1X,E9.3,1X,E9.3,
     +  2X,I6,1X,F6.2,1X,F8.5,1X,A4,1X,A18,6X,A19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gum_sn.dat')
      write(6,*) '....Loading file: gum_sn.dat'
      do i__=1,49814
        read(1,'(A160)')ar__5
        read(ar__5,6)
     +  RA2kdeg_5(i__),DE2kdeg_5(i__),RAdeg_5(i__),DEdeg_5(i__),
     +  r_5(i__),pmRA_5(i__),pmDE_5(i__),RV_5(i__),VMAG_4(i__),
     +  z_2(i__),SNt(i__),ID_5(i__),objID_1(i__)
c    ..............Just test output...........
        write(6,6)
     +  RA2kdeg_5(i__),DE2kdeg_5(i__),RAdeg_5(i__),DEdeg_5(i__),
     +  r_5(i__),pmRA_5(i__),pmDE_5(i__),RV_5(i__),VMAG_4(i__),
     +  z_2(i__),SNt(i__),ID_5(i__),objID_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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