FORTRAN Generation
(/./ftp/cats/J/ApJ/959/102)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/959/102 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-May-14
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/ApJ/959/102  Physical properties of cool SGs in LMC  (Dorn-Wallenstein+, 2023)
*================================================================================
*Physical properties of 5000 cool Large Magellanic Cloud supergiants with Gaia XP
*spectra: a detailed portrait of the upper H-R diagram hints at missing supernova
*progenitors.
*    Dorn-Wallenstein T.Z., Neugent K.F., Levesque E.M.
*   <Astrophys. J., 959, 102 (2023)>
*   =2023ApJ...959..102D
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'tableb1.dat'	! Temperatures and luminosities for 6146 AFGKM
                             supergiants

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

      integer*8     Gaia       (nr__) ! Gaia DR3 source ID
      real*4        logTeff    (nr__) ! ([K]) [3.3/4.4] Log effective temperature
      real*4        logL       (nr__) ! ([Lsun]) [3/5.72] Log luminosity
      character*1   Ref        (nr__) ! [NY] Source of known stellar properties (1)
      integer*4     Qual       (nr__) ! [0/2]? Quality of predicted values (2)
*Note (1): Reference code as follows:
*    N = stars from N12 (Neugent+ 2012, J/ApJ/749/177; 641 occurrences)
*    Y = stars from Y21 (Yang+ 2021, J/A+A/646/A141), for which logTeff and logL
*          are derived by our regressor (5505 occurrences)
*Note (2): Code as follows:
*    0 = stars from N12 (Neugent+ 2012, J/ApJ/749/177) or
*         Y21 (Yang+ 2021, J/A+A/646/A141) that fall within the temperature
*         and luminosity boundaries of the training set (1272 occurrences)
*    1 = low-L stars (4449 occurrences)
*    2 = out-of-boundary (OOB) stars (425 occurrences)
*
*           See Section 3.

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

C  Loading file 'tableb1.dat'	! Temperatures and luminosities for 6146 AFGKM
*                             supergiants

C  Format for file interpretation

    1 format(I19,1X,F5.3,1X,F5.3,1X,A1,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tableb1.dat')
      write(6,*) '....Loading file: tableb1.dat'
      do i__=1,6146
        read(1,'(A35)')ar__
        read(ar__,1)
     +  Gaia(i__),logTeff(i__),logL(i__),Ref(i__),Qual(i__)
        if(ar__(35:35) .EQ. '') Qual(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Gaia(i__),logTeff(i__),logL(i__),Ref(i__),Qual(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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