FORTRAN Generation
(/./ftp/cats/J/AJ/161/193)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/161/193 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-Aug-16
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/AJ/161/193          BVRI light curve of RR Lyrae V* AX UMa         (Li+, 2021)
*================================================================================
*AX UMa, an ab-type RR Lyrae star with a rapidly decreasing pulsation period, and
*its binarity.
*    Li L.-J., Qian S.-B., Zhu L.-Y., He J.-J., Fang X.-H.
*   <Astron. J., 161, 193-193 (2021)>
*   =2021AJ....161..193L    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Times of light maximum and corresponding errors
                                  of AX UMa determined from photometric surveys
                                  and our observations

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

      real*8        HJD        (nr__) ! (d) [27291/58512] Times of light maximum, HJD-2400000
      real*4        e_HJD      (nr__) ! (d) [0.0002/0.01]? Uncertainty in HJD
      character*3   Method     (nr__) ! Method used (1)
      integer*4     Epoch      (nr__) ! [-33808/24550] Epoch (cycle)
      real*8        O_C        (nr__) ! (d) [-6.66/0.08] O-C, Equation 2
      character*7   Ref        (nr__) ! Reference or source (2)
*Note (1): Methods as follows:
*    vis = visual (13 occurrences)
*    pg = photographic (134 occurrences)
*    ccd = Charge Coupled Device (24 occurrences)
*Note (2): References and sources as follows:
*    Rom64 = Romano, 1964MmSAI..35..101R
*    I2862 = Romano, 1986IBVS.2862....1R
*    Bon85 = Boninsegna, 1985GEOCR...8.....B
*    Sch95 = Schmidt+, 1995AJ....109.1239S
*    Wils06 = Wils+, 2006, J/MNRAS/368/1757
*    I5854 = Le Borgne+, 2008IBVS.5854....1L
*    VS10 = Vandenbroere J., Salmon G., 2010, GEOS Circular RR45;
*    I5959 = Hubscher & Monninger, 2011, J/other/IBVS/5959
*    I6149 = Hubscher & Lehmann, 2015IBVS.6149....1H
*    GaiaDr2 = Gaia Collaboration+, 2018, I/345
*    Fer15 = data from the GEOS RR Lyr Database observed by V. Ferrand, 2015.
*    AAVSO = American Association of Variable Star Observers
*    CRTS = Drake+, 2009, J/ApJ/696/870
*    DASCH = Grindlay+, 2009ASPC..410..101G
*    NSVS = Wozniak+, 2004AJ....127.2436W
*    QES = Alsubai+, 2013AcA....63..465A

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

C  Declarations for 'fig9.dat'	! Color curves of AX UMa obtained by Sino-Thai
                                  70cm telescope

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

      real*4        Phase      (nr__1) ! [-0.25/1.25] Phase
      real*4        Color      (nr__1) ! [-0.01/0.63] Color Index
      character*2   Ref_1      (nr__1) ! Color Code; B-V, R-I, V-I or V-R

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

C  Loading file 'table1.dat'	! Times of light maximum and corresponding errors
*                                  of AX UMa determined from photometric surveys
*                                  and our observations

C  Format for file interpretation

    1 format(F10.4,1X,F6.4,1X,A3,1X,I6,1X,F7.4,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,171
        read(1,'(A44)')ar__
        read(ar__,1)
     +  HJD(i__),e_HJD(i__),Method(i__),Epoch(i__),O_C(i__),Ref(i__)
        if(ar__(12:17) .EQ. '') e_HJD(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  HJD(i__),e_HJD(i__),Method(i__),Epoch(i__),O_C(i__),Ref(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig9.dat'	! Color curves of AX UMa obtained by Sino-Thai
*                                  70cm telescope

C  Format for file interpretation

    2 format(F6.3,1X,F6.3,1X,A2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig9.dat')
      write(6,*) '....Loading file: fig9.dat'
      do i__=1,1574
        read(1,'(A16)')ar__1
        read(ar__1,2)Phase(i__),Color(i__),Ref_1(i__)
c    ..............Just test output...........
        write(6,2)Phase(i__),Color(i__),Ref_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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