FORTRAN Generation
(/./ftp/cats/J/ApJS/241/6)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/241/6 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-Sep-10
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/ApJS/241/6     Searching for super-fast rotators using PS1     (Chang+, 2019)
*================================================================================
*Searching for super-fast rotators using the Pan-STARRS 1.
*    Chang C.-K., Lin H.-W., Ip W.-H., Chen W.-P., Yeh T.-S., Chambers K.C.,
*    Magnier E.A., Huber M.E., Flewelling H.A., Waters C.Z., Wainscoat R.J.,
*    Schultz A.S.B.
*   <Astrophys. J. Suppl. Ser., 241, 6 (2019)>
*   =2019ApJS..241....6C
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! The 876 reliable rotation periods

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

      character*8   ID         (nr__) ! Compact object identifier (asteroid number or
*                                 MPC packed format)
      character*1   f_ID       (nr__) ! [ab] Flag on ID (1)
      integer*4     Planet     (nr__) ! [1489/399139]? Asteroid number
      character*15  Name       (nr__) ! Name or preliminary designation
      real*4        a          (nr__) ! (AU) [1.8/4] Semi-major axis
      real*4        e          (nr__) ! [0.01/0.5] Eccentricity
      real*4        i          (nr__) ! (deg) [0.3/32] Inclination
      real*4        Delta      (nr__) ! (AU) [4.9/8] Mean heliocentric distance
      real*4        r          (nr__) ! (AU) [4.9/8] Mean geodesic distance
      real*4        alpha      (nr__) ! (deg) [4.9/8] Mean phase angle
      real*4        Diam       (nr__) ! (km) [0.2/28.5] Diameter
      real*4        HMag       (nr__) ! (mag) [11.7/20.5] Absolute magnitude
      real*4        e_HMag     (nr__) ! (mag) [0.03/0.5] Uncertainty in HMag
      real*4        mag        (nr__) ! (mag) [12/22] Apparent magnitude in PS1 w_P1_ filter
      real*4        e_mag      (nr__) ! (mag) [0/0.2] Uncertainty in mag
      real*4        Per        (nr__) ! (h) [0.5/126.4] Derived rotation period
      real*4        e_Per      (nr__) ! (h) [0/9.3] Uncertainty in P
      real*4        Deltam     (nr__) ! (mag) [0.09/1.5] Light-curve amplitude
      integer*4     U          (nr__) ! [2/3] Rotation period quality code
*                                 (3=highly reliable; 2=some ambiguity)
      real*4        a_         (nr__) ! [-5.6/3.1]? SDSS color a^*^ (2)
      real*4        e_a_       (nr__) ! [0/0.6]? Uncertainty in a*
      real*4        i_z        (nr__) ! [-2.4/2]? SDSS (i-z) color index
      real*4        e_i_z      (nr__) ! [0/1.4]? Uncertainty in i-z
      character*1   Type       (nr__) ! Spectral type ("C"=385 occurrences,
*                                 "S"=248 occurrences or "V"=128 occurrences)
*Note (1): Flag as follows:
*    a = Rotation period measurements available in the LCDB.
*    b = Long-period objects with partial coverage on rotational phase.
*Note (2): Where a*=0.89*(g-r)+0.45*(r-i)-0.57.

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

C  Loading file 'table3.dat'	! The 876 reliable rotation periods

C  Format for file interpretation

    1 format(
     +  A8,1X,A1,1X,I6,1X,A15,1X,F4.2,1X,F4.2,1X,F4.1,1X,F3.1,1X,F3.1,
     +  1X,F3.1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F6.2,1X,
     +  F4.2,1X,F4.2,1X,I1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,876
        read(1,'(A130)')ar__
        read(ar__,1)
     +  ID(i__),f_ID(i__),Planet(i__),Name(i__),a(i__),e(i__),i(i__),
     +  Delta(i__),r(i__),alpha(i__),Diam(i__),HMag(i__),e_HMag(i__),
     +  mag(i__),e_mag(i__),Per(i__),e_Per(i__),Deltam(i__),U(i__),
     +  a_(i__),e_a_(i__),i_z(i__),e_i_z(i__),Type(i__)
        if(ar__(12:17) .EQ. '') Planet(i__) = iNULL__
        if(ar__(108:112) .EQ. '') a_(i__) = rNULL__
        if(ar__(114:117) .EQ. '') e_a_(i__) = rNULL__
        if(ar__(119:123) .EQ. '') i_z(i__) = rNULL__
        if(ar__(125:128) .EQ. '') e_i_z(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  ID(i__),f_ID(i__),Planet(i__),Name(i__),a(i__),e(i__),i(i__),
     +  Delta(i__),r(i__),alpha(i__),Diam(i__),HMag(i__),e_HMag(i__),
     +  mag(i__),e_mag(i__),Per(i__),e_Per(i__),Deltam(i__),U(i__),
     +  a_(i__),e_a_(i__),i_z(i__),e_i_z(i__),Type(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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