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

Conversion of standardized ReadMe file for file /./ftp/cats/VI/137 into FORTRAN code for reading data files line by line.

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          ! source identifier (sourceExtendedId)
      real*8        RA2kdeg     ! (deg) Right ascension ICRS, Epoch=J2000 (G1)
      real*8        DE2kdeg     ! (deg) Declination ICRS, Epoch=J2000 (G1)
      real*4        VMAG        ! (mag) V-band absolute magnitude (meanAbsoluteV)
      real*4        Gmag        ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag       ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag       ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag       ! (mag) Gaia Grvs magnitude (860+/-28nm) (magGRvs)
      real*8        RAdeg       ! (deg) Right ascension (ICRS, Epoch=J2010) (alpha)
      real*8        DEdeg       ! (deg) Declination (ICRS, Epoch=J2010) (delta)
      real*8        r           ! (pc) Barycentric distance (distance)
      real*8        pmRA        ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*8        pmDE        ! (mas/yr) Proper motion in Declination (muDelta)
      real*8        RV          ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        V_I         ! (mag) Intrinsic color V-I (colorVminusI)
      real*4        Av          ! (mag) Absorption (Av)
      real*4        Age         ! (Gyr) Estimated age (age)
      real*4        v_a_Fe_     ! ([Sun]) {alpha}-elements enhancement (alphaFe)
      real*4        bAlb        ! [0/1] Bond albedo (for exoplanet) (bondAlbedo)
      real*4        e           ! Eccentricity (for binaries) (eccentricity)
      real*4        v_Fe_H_     ! ([Sun]) Metallicity [Fe/H] (feH)
      integer*4     fI          ! [0,1] 1 (true) if interacting with a companion
*                                 (flagInteracting)
      real*4        gAlb        ! Geometric albedo (for exoplanet) (geomAlbedo)
      integer*4     fm          ! [0,1] 1 (true) if photocenter moves
*                                 (hasPhotocenterMotion)
      integer*4     Host        ! Host galaxy: 1=Milky Way, 2=SMC, 3=LMC (host)
      real*4        i           ! (deg) Orbit Inclination (inclination)
      real*4        logg        ! ([cm/s2]) Gravity (log) (logg)
      real*4        Omega       ! (deg) Longitude of ascending node {Omega}
*                                 (longitudeAscendingNode)
      real*8        Mass        ! (Sun) Stellar mass (mass)
      real*4        Mbol        ! (mag) Bolometric absolute magnitude (mbol)
      integer*4     nc          ! [1-5] Number of components in the system
      integer*4     nt          ! [1-5] Total number of objects in the system
      real*4        P           ! (d) Orbital period (orbitPeriod)
      real*4        omega_1     ! (deg) Periastron argument {omega}(periastronArgument)
      real*4        T0          ! (d) Periastron date (MJD) (periastronDate)
      real*4        Phase       ! [0/1] Phase if an exoplanet (phase)
      integer*4     Pop         ! [1/3]?=99 Stellar population the object
*                                        belongs to (population)
      real*8        BeEnv       ! Envelope characteristic for Be star (Renv/R*)
*                                 (rEnvRStar)
      real*4        Radius      ! (Rsun) Stellar radius (radius)
      real*4        a           ! (au) Semi-major axis (in binary system)
*                                 (semimajorAxis)
      integer*4     Teff        ! (K) Effective temperature (teff)
      real*4        Vamp        ! (mag) Variability amplitude (variabilityAmplitude)
      real*8        Vper        ! (d) Variability period (variabilityPeriod)
      real*4        Vphase      ! [0/1] Variability phase (variabilityPhase)
      character*4   Vtype       ! Variability type (variabilityType) (2)
      real*4        vsini       ! (km/s) Rotational velocity
      integer*8     recno       ! [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        ! source identifier (sourceExtendedId)
      real*8        RA2kdeg_1   ! (deg) Right ascension ICRS, Epoch=J2000 (G1)
      real*8        DE2kdeg_1   ! (deg) Declination ICRS, Epoch=J2000 (G1)
      real*4        VMAG_1      ! (mag) V-band absolute magnitude (meanAbsoluteV)
      real*4        Gmag_1      ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag_1     ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag_1     ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag_1     ! (mag) Gaia Grvs magnitude (860+/-28nm) (magGRvs)
      real*8        RAdeg_1     ! (deg) Right ascension (ICRS, Epoch=J2010) (alpha)
      real*8        DEdeg_1     ! (deg) Declination (ICRS, Epoch=J2010) (delta)
      real*8        r_1         ! (pc) Barycentric distance (distance)
      real*8        pmRA_1      ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*8        pmDE_1      ! (mas/yr) Proper motion in Declination (muDelta)
      real*8        RV_1        ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        V_I_1       ! (mag) Intrinsic color V-I (colorVminusI)
      real*4        Av_1        ! (mag) Absorption (Av)
      real*4        Age_1       ! (Gyr) Estimated age (age)
      real*4        v_a_Fe__1   ! ([Sun]) {alpha}-elements enhancement (alphaFe)
      real*4        bAlb_1      ! [0/1] Bond albedo (for exoplanet) (bondAlbedo)
      real*4        e_1         ! Eccentricity (for binaries) (eccentricity)
      real*4        v_Fe_H__1   ! ([Sun]) Metallicity [Fe/H] (feH)
      integer*4     fI_1        ! [0,1] 1 (true) if interacting with a companion
*                                 (flagInteracting)
      real*4        gAlb_1      ! Geometric albedo (for exoplanet) (geomAlbedo)
      integer*4     fm_1        ! [0,1] 1 (true) if photocenter moves
*                                 (hasPhotocenterMotion)
      integer*4     Host_1      ! Host galaxy: 1=Milky Way, 2=SMC, 3=LMC (host)
      real*4        i_1         ! (deg) Orbit Inclination (inclination)
      real*4        logg_1      ! ([cm/s2]) Gravity (log) (logg)
      real*4        Omega_2     ! (deg) Longitude of ascending node {Omega}
*                                 (longitudeAscendingNode)
      real*8        Mass_1      ! (Sun) Stellar mass (mass)
      real*4        Mbol_1      ! (mag) Bolometric absolute magnitude (mbol)
      integer*4     nc_1        ! [1-5] Number of components in the system
      integer*4     nt_1        ! [1-5] Total number of objects in the system
      real*4        P_1         ! (d) Orbital period (orbitPeriod)
      real*4        omega_3     ! (deg) Periastron argument {omega}(periastronArgument)
      real*4        T0_1        ! (d) Periastron date (MJD) (periastronDate)
      real*4        Phase_1     ! [0/1] Phase if an exoplanet (phase)
      integer*4     Pop_1       ! [1/3]?=99 Stellar population the object
*                                        belongs to (population)
      real*8        BeEnv_1     ! Envelope characteristic for Be star (Renv/R*)
*                                 (rEnvRStar)
      real*4        Radius_1    ! (Rsun) Stellar radius (radius)
      real*4        a_1         ! (au) Semi-major axis (in binary system)
*                                 (semimajorAxis)
      integer*4     Teff_1      ! (K) Effective temperature (teff)
      real*4        Vamp_1      ! (mag) Variability amplitude (variabilityAmplitude)
      real*8        Vper_1      ! (d) Variability period (variabilityPeriod)
      real*4        Vphase_1    ! [0/1] Variability phase (variabilityPhase)
      character*4   Vtype_1     ! Variability type (variabilityType) (2)
      real*4        vsini_1     ! (km/s) Rotational velocity
      integer*8     recno_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        ! source identifier (sourceExtendedId)
      real*8        RA2kdeg_2   ! (deg) Right ascension ICRS, Epoch=J2000 (G1)
      real*8        DE2kdeg_2   ! (deg) Declination ICRS, Epoch=J2000 (G1)
      real*4        VMAG_2      ! (mag) V-band absolute magnitude (meanAbsoluteV)
      real*4        Gmag_2      ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag_2     ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag_2     ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag_2     ! (mag) Gaia Grvs magnitude (860+/-28nm) (magGRvs)
      real*8        RAdeg_2     ! (deg) Right ascension (ICRS, Epoch=J2010) (alpha)
      real*8        DEdeg_2     ! (deg) Declination (ICRS, Epoch=J2010) (delta)
      real*8        r_2         ! (pc) Barycentric distance (distance)
      real*8        pmRA_2      ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*8        pmDE_2      ! (mas/yr) Proper motion in Declination (muDelta)
      real*8        RV_2        ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        V_I_2       ! (mag) Intrinsic color V-I (colorVminusI)
      real*4        Av_2        ! (mag) Absorption (Av)
      real*4        Age_2       ! (Gyr) Estimated age (age)
      real*4        v_a_Fe__2   ! ([Sun]) {alpha}-elements enhancement (alphaFe)
      real*4        bAlb_2      ! [0/1] Bond albedo (for exoplanet) (bondAlbedo)
      real*4        e_2         ! Eccentricity (for binaries) (eccentricity)
      real*4        v_Fe_H__2   ! ([Sun]) Metallicity [Fe/H] (feH)
      integer*4     fI_2        ! [0,1] 1 (true) if interacting with a companion
*                                 (flagInteracting)
      real*4        gAlb_2      ! Geometric albedo (for exoplanet) (geomAlbedo)
      integer*4     fm_2        ! [0,1] 1 (true) if photocenter moves
*                                 (hasPhotocenterMotion)
      integer*4     Host_2      ! Host galaxy: 1=Milky Way, 2=SMC, 3=LMC (host)
      real*4        i_2         ! (deg) Orbit Inclination (inclination)
      real*4        logg_2      ! ([cm/s2]) Gravity (log) (logg)
      real*4        Omega_4     ! (deg) Longitude of ascending node {Omega}
*                                 (longitudeAscendingNode)
      real*8        Mass_2      ! (Sun) Stellar mass (mass)
      real*4        Mbol_2      ! (mag) Bolometric absolute magnitude (mbol)
      integer*4     nc_2        ! [1-5] Number of components in the system
      integer*4     nt_2        ! [1-5] Total number of objects in the system
      real*4        P_2         ! (d) Orbital period (orbitPeriod)
      real*4        omega_5     ! (deg) Periastron argument {omega}(periastronArgument)
      real*4        T0_2        ! (d) Periastron date (MJD) (periastronDate)
      real*4        Phase_2     ! [0/1] Phase if an exoplanet (phase)
      integer*4     Pop_2       ! [1/3]?=99 Stellar population the object
*                                        belongs to (population)
      real*8        BeEnv_2     ! Envelope characteristic for Be star (Renv/R*)
*                                 (rEnvRStar)
      real*4        Radius_2    ! (Rsun) Stellar radius (radius)
      real*4        a_2         ! (au) Semi-major axis (in binary system)
*                                 (semimajorAxis)
      integer*4     Teff_2      ! (K) Effective temperature (teff)
      real*4        Vamp_2      ! (mag) Variability amplitude (variabilityAmplitude)
      real*8        Vper_2      ! (d) Variability period (variabilityPeriod)
      real*4        Vphase_2    ! [0/1] Variability phase (variabilityPhase)
      character*4   Vtype_2     ! Variability type (variabilityType) (2)
      real*4        vsini_2     ! (km/s) Rotational velocity
      integer*8     recno_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        ! Source identifier (sourceExtendedId)
      real*8        RA2kdeg_3   ! (deg) Right ascension ICRS, Epoch=J2000 (G1)
      real*8        DE2kdeg_3   ! (deg) Declination ICRS, Epoch=J2000 (G1)
      real*8        Vmag_3      ! (mag) V-band intrinsic (extinction free) apparent
*                                  mean magnitude (meanAbsoluteV)
      real*4        Gmag_3      ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag_3     ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag_3     ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag_3     ! (mag) Gaia Grvs magnitude (860+/-28nm) (magGRvs)nm)
      real*8        RAdeg_3     ! (deg) Right ascension (ICRS, Epoch=J2010) (alpha)
      real*8        DEdeg_3     ! (deg) Declination (ICRS, Epoch=J2010) (delta)
      real*4        r_3         ! (pc) Barycentric distance (distance)
      real*4        pmRA_3      ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*4        pmDE_3      ! (mas/yr) Proper motion in Declination (muDelta)
      real*8        RV_3        ! (km/s) Barycentric radial velocity (radialVelocity)
      integer*4     SNID        ! [-1] Attached supernova ID (not filled)
      real*4        V_I_3       ! (mag) Intrinsic color V-I (colorVminusI)
      integer*4     Htype       ! [0-8] Hubble type (hubbleType) (3)
      real*4        Av_3        ! (mag) Absorption in the V-band
      real*8        z           ! Redshift (redShift)
      real*4        B_T         ! [0,1] Bulge to total flux ratio (4)
      real*4        Brad        ! (arcsec) Bulge equivalent radius (4)
      real*4        Bb_a        ! [0,1] Bulge aspect ratio b/a (4)
      real*8        Bpa         ! (deg) [-180/180] Bulge position angle (4)
      real*8        Dsize       ! (arcsec) Disk scale length (4)
      real*4        Db_a        ! [0,1] Disk aspect ratio b/a (4)
      real*8        Dpa         ! (deg) Disk position angle (4)
      integer*4     recno_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   ! (deg) Right ascension (ICRS), at Epoch J2000 (G1)
      real*8        DE2kdeg_4   ! (deg) Declination (ICRS), at Epoch J2000 (G1)
      real*8        RAdeg_4     ! (deg) Right ascension (ICRS), at Epoch J2010 (alpha)
      real*8        DEdeg_4     ! (deg) Right ascension (ICRS), at Epoch J2010 (delta)
      integer*4     r_4         ! (Mpc) Barycentric distance (distance)
      real*8        pmRA_4      ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*8        pmDE_4      ! (mas/yr) Proper motion in Declination (muDelta)
      integer*4     RV_4        ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        Gmag_4      ! (mag) Gaia G-band magnitude (350-1050nm) (magG)
      real*4        GBmag_4     ! (mag) Gaia Gbp band magnitude (350-770nm) (magGBp)
      real*4        GRmag_4     ! (mag) Gaia Grp band magnitude (650-1050nm) (magGRp)
      real*4        Gsmag_4     ! (mag) Gaia/Grvs band magnitude (850-880nm) (magGRvs)
      real*4        Av_4        ! (mag) Galactic absorption in the V-band
      real*4        z_1         ! Redshift (redShift)
      real*4        slope       ! Slope of the continuum
      integer*4     W           ! Total equivalent width of the emission lines
      character*18  ID_4        ! Source identifier (sourceExtendedId)
      character*19  objID       ! 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   ! (deg) Right ascension (ICRS), at Epoch J2000 (G1)
      real*8        DE2kdeg_5   ! (deg) Declination (ICRS), at Epoch J2000 (G1)
      real*8        RAdeg_5     ! (deg) Right ascension (ICRS), at Epoch J2010 (alpha)
      real*8        DEdeg_5     ! (deg) Right ascension (ICRS), at Epoch J2010 (delta)
      real*8        r_5         ! (Mpc) Barycentric distance distance)
      real*4        pmRA_5      ! (mas/yr) Proper motion in RA (*cos(DE)) (muAlpha)
      real*4        pmDE_5      ! (mas/yr) Proper motion in Declination (muDelta)
      integer*4     RV_5        ! (km/s) Barycentric radial velocity (radialVelocity)
      real*4        VMAG_4      ! (mag) V-band absolute magnitude (meanAbsoluteV)
      real*8        z_2         ! Redshift (redShift)
      character*4   SNt         ! Type of the supernova (gmdbtype)
      character*18  ID_5        ! Source identifier (sourceExtendedId)
      character*19  objID_1     ! 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,RA2kdeg,DE2kdeg,VMAG,Gmag,GBmag,GRmag,Gsmag,RAdeg,DEdeg,r,
     +  pmRA,pmDE,RV,V_I,Av,Age,v_a_Fe_,bAlb,e,v_Fe_H_,fI,gAlb,fm,
     +  Host,i,logg,Omega,Mass,Mbol,nc,nt,P,omega_1,T0,Phase,Pop,
     +  BeEnv,Radius,a,Teff,Vamp,Vper,Vphase,Vtype,vsini,recno
c    ..............Just test output...........
        write(6,1)
     +  ID,RA2kdeg,DE2kdeg,VMAG,Gmag,GBmag,GRmag,Gsmag,RAdeg,DEdeg,r,
     +  pmRA,pmDE,RV,V_I,Av,Age,v_a_Fe_,bAlb,e,v_Fe_H_,fI,gAlb,fm,
     +  Host,i,logg,Omega,Mass,Mbol,nc,nt,P,omega_1,T0,Phase,Pop,
     +  BeEnv,Radius,a,Teff,Vamp,Vper,Vphase,Vtype,vsini,recno
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,RA2kdeg_1,DE2kdeg_1,VMAG_1,Gmag_1,GBmag_1,GRmag_1,
     +  Gsmag_1,RAdeg_1,DEdeg_1,r_1,pmRA_1,pmDE_1,RV_1,V_I_1,Av_1,
     +  Age_1,v_a_Fe__1,bAlb_1,e_1,v_Fe_H__1,fI_1,gAlb_1,fm_1,Host_1,
     +  i_1,logg_1,Omega_2,Mass_1,Mbol_1,nc_1,nt_1,P_1,omega_3,T0_1,
     +  Phase_1,Pop_1,BeEnv_1,Radius_1,a_1,Teff_1,Vamp_1,Vper_1,
     +  Vphase_1,Vtype_1,vsini_1,recno_1
c    ..............Just test output...........
        write(6,2)
     +  ID_1,RA2kdeg_1,DE2kdeg_1,VMAG_1,Gmag_1,GBmag_1,GRmag_1,
     +  Gsmag_1,RAdeg_1,DEdeg_1,r_1,pmRA_1,pmDE_1,RV_1,V_I_1,Av_1,
     +  Age_1,v_a_Fe__1,bAlb_1,e_1,v_Fe_H__1,fI_1,gAlb_1,fm_1,Host_1,
     +  i_1,logg_1,Omega_2,Mass_1,Mbol_1,nc_1,nt_1,P_1,omega_3,T0_1,
     +  Phase_1,Pop_1,BeEnv_1,Radius_1,a_1,Teff_1,Vamp_1,Vper_1,
     +  Vphase_1,Vtype_1,vsini_1,recno_1
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,RA2kdeg_2,DE2kdeg_2,VMAG_2,Gmag_2,GBmag_2,GRmag_2,
     +  Gsmag_2,RAdeg_2,DEdeg_2,r_2,pmRA_2,pmDE_2,RV_2,V_I_2,Av_2,
     +  Age_2,v_a_Fe__2,bAlb_2,e_2,v_Fe_H__2,fI_2,gAlb_2,fm_2,Host_2,
     +  i_2,logg_2,Omega_4,Mass_2,Mbol_2,nc_2,nt_2,P_2,omega_5,T0_2,
     +  Phase_2,Pop_2,BeEnv_2,Radius_2,a_2,Teff_2,Vamp_2,Vper_2,
     +  Vphase_2,Vtype_2,vsini_2,recno_2
c    ..............Just test output...........
        write(6,3)
     +  ID_2,RA2kdeg_2,DE2kdeg_2,VMAG_2,Gmag_2,GBmag_2,GRmag_2,
     +  Gsmag_2,RAdeg_2,DEdeg_2,r_2,pmRA_2,pmDE_2,RV_2,V_I_2,Av_2,
     +  Age_2,v_a_Fe__2,bAlb_2,e_2,v_Fe_H__2,fI_2,gAlb_2,fm_2,Host_2,
     +  i_2,logg_2,Omega_4,Mass_2,Mbol_2,nc_2,nt_2,P_2,omega_5,T0_2,
     +  Phase_2,Pop_2,BeEnv_2,Radius_2,a_2,Teff_2,Vamp_2,Vper_2,
     +  Vphase_2,Vtype_2,vsini_2,recno_2
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,RA2kdeg_3,DE2kdeg_3,Vmag_3,Gmag_3,GBmag_3,GRmag_3,
     +  Gsmag_3,RAdeg_3,DEdeg_3,r_3,pmRA_3,pmDE_3,RV_3,SNID,V_I_3,
     +  Htype,Av_3,z,B_T,Brad,Bb_a,Bpa,Dsize,Db_a,Dpa,recno_3
c    ..............Just test output...........
        write(6,4)
     +  ID_3,RA2kdeg_3,DE2kdeg_3,Vmag_3,Gmag_3,GBmag_3,GRmag_3,
     +  Gsmag_3,RAdeg_3,DEdeg_3,r_3,pmRA_3,pmDE_3,RV_3,SNID,V_I_3,
     +  Htype,Av_3,z,B_T,Brad,Bb_a,Bpa,Dsize,Db_a,Dpa,recno_3
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,DE2kdeg_4,RAdeg_4,DEdeg_4,r_4,pmRA_4,pmDE_4,RV_4,
     +  Gmag_4,GBmag_4,GRmag_4,Gsmag_4,Av_4,z_1,slope,W,ID_4,objID
c    ..............Just test output...........
        write(6,5)
     +  RA2kdeg_4,DE2kdeg_4,RAdeg_4,DEdeg_4,r_4,pmRA_4,pmDE_4,RV_4,
     +  Gmag_4,GBmag_4,GRmag_4,Gsmag_4,Av_4,z_1,slope,W,ID_4,objID
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,DE2kdeg_5,RAdeg_5,DEdeg_5,r_5,pmRA_5,pmDE_5,RV_5,
     +  VMAG_4,z_2,SNt,ID_5,objID_1
c    ..............Just test output...........
        write(6,6)
     +  RA2kdeg_5,DE2kdeg_5,RAdeg_5,DEdeg_5,r_5,pmRA_5,pmDE_5,RV_5,
     +  VMAG_4,z_2,SNt,ID_5,objID_1
c    .......End.of.Just test output...........
      end do
      close(1)

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