FORTRAN Generation
(/./ftp/cats/J/A_A/406/75)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/406/75 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-Apr-23
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/A+A/406/75        JHKs magnitudes of variables in NGC 5128    (Rejkuba+, 2003)
*================================================================================
*Long period variables in NGC 5128. I. Catalogue.
*     Rejkuba M., Minniti D., Silva D.R.
*    <Astron. Astrophys. 406, 75 (2003)>
*    =2003A&A...406...75R
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Near-IR multi-epoch photometry observing log

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

      character*2   Field      (nr__) ! Field (F1 or F2)
      real*8        JD         (nr__) ! (d) Julian date
      integer*4     Exp        (nr__) ! (min) Exposure time
      real*4        X          (nr__) ! Mean airmass during the observing sequence
      real*4        Seeing     (nr__) ! (arcsec) Seeing
      character*5   FilEp      (nr__) ! Filter epoch
      integer*4     Instr      (nr__) ! [1/2]? Instrument (1)
*Note (1): the codes are as follows:
*    1 = Observed with SOFI at NTT.
*    2 = Observed with LW arm of ISAAC at VLT.

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

C  Declarations for 'table2.dat'	! Data for a sample of variables with light curves
                            plotted in Fig. 5

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

      character*8   v_RMS2003_ (nr__1) ! Identification, FN NNNNN
      real*8        xpos       (nr__1) ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos       (nr__1) ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_1    (nr__1) ! Reference epoch, defined in table1.dat
      integer*4     Per        (nr__1) ! (d) ? Best fitting period
      real*4        A          (nr__1) ! (mag) ? Semi-amplitude
      real*4        chi2       (nr__1) ! ? Reduced chi2 of the sine-curve fit
      real*4        signif     (nr__1) ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag      (nr__1) ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag       (nr__1) ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag   (nr__1) ! (mag) Mean Ks-band magnitude from the sine-curve fit
*Note (1): Magnitude 99.99 denotes that no measurement for that star was obtained

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

C  Declarations for 'table3.dat'	! LPVs in Field 1

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

      character*8   v_RMS2003__1(nr__2) ! Identification, FN NNNNN
      real*8        xpos_1     (nr__2) ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos_1     (nr__2) ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_2    (nr__2) ! Reference epoch, defined in table1.dat
      integer*4     Per_1      (nr__2) ! (d) ? Best fitting period
      real*4        A_1        (nr__2) ! (mag) ? Semi-amplitude
      real*4        chi2_1     (nr__2) ! ? Reduced chi2 of the sine-curve fit
      real*4        signif_1   (nr__2) ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag_1    (nr__2) ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag_1     (nr__2) ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag_1 (nr__2) ! (mag) Mean Ks-band magnitude from the sine-curve fit
*Note (1): Magnitude 99.99 denotes that no measurement for that star was obtained

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

C  Declarations for 'table4.dat'	! LPVs in Field 2

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

      character*8   v_RMS2003__2(nr__3) ! Identification, FN NNNNN
      real*8        xpos_2     (nr__3) ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos_2     (nr__3) ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_3    (nr__3) ! Reference epoch, defined in table1.dat
      integer*4     Per_2      (nr__3) ! (d) ? Best fitting period
      real*4        A_2        (nr__3) ! (mag) ? Semi-amplitude
      real*4        chi2_2     (nr__3) ! ? Reduced chi2 of the sine-curve fit
      real*4        signif_2   (nr__3) ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag_2    (nr__3) ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag_2     (nr__3) ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag_2 (nr__3) ! (mag) Mean Ks-band magnitude from the sine-curve fit
*Note (1): Magnitude 99.99 denotes that no measurement for that star was obtained

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

C  Declarations for 'table5.dat'	! Additional variable stars in Field 1 for which
                            periods could not be determined

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

      character*8   v_RMS2003__3(nr__4) ! Identification, FN NNNNN
      real*8        xpos_3     (nr__4) ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos_3     (nr__4) ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_4    (nr__4) ! Reference epoch, defined in table1.dat
      integer*4     Per_3      (nr__4) ! (d) ? Best fitting period
      real*4        A_3        (nr__4) ! (mag) ? Semi-amplitude
      real*4        chi2_3     (nr__4) ! ? Reduced chi2 of the sine-curve fit
      real*4        signif_3   (nr__4) ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag_3    (nr__4) ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag_3     (nr__4) ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag_3 (nr__4) ! (mag) Mean Ks-band magnitude from the sine-curve fit
*Note (1): Magnitude 99.99 denotes that no measurement for that star was obtained

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

C  Declarations for 'table6.dat'	! Additional variable stars in Field 2 for which
                            periods could not be determined

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

      character*8   v_RMS2003__4(nr__5) ! Identification, FN NNNNN
      real*8        xpos_4     (nr__5) ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos_4     (nr__5) ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_5    (nr__5) ! Reference epoch, defined in table1.dat
      integer*4     Per_4      (nr__5) ! (d) ? Best fitting period
      real*4        A_4        (nr__5) ! (mag) ? Semi-amplitude
      real*4        chi2_4     (nr__5) ! ? Reduced chi2 of the sine-curve fit
      real*4        signif_4   (nr__5) ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag_4    (nr__5) ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag_4     (nr__5) ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag_4 (nr__5) ! (mag) Mean Ks-band magnitude from the sine-curve fit
*Note (1): Magnitude 99.99 denotes that no measurement for that star was obtained

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

C  Declarations for 'table7.dat'	! Individual measurements of the K-band magnitudes
                            for the 3 red variable stars detected in Field 1

      integer*4 nr__6
      parameter (nr__6=12020)	! Number of records
      character*41 ar__6  	! Full-size record

      character*8   v_RMS2003__5(nr__6) ! Identification, FN NNNN
      real*8        xpos_5     (nr__6) ! (pix) x position (G1) (1)
      real*8        ypos_5     (nr__6) ! (pix) y position (G1) (1)
      character*5   FilEp_6    (nr__6) ! Filter epoch
      real*4        Kmag       (nr__6) ! (mag) ?=99.99 Filter epoch K-band magnitude (2)
      real*4        e_Kmag     (nr__6) ! (mag) ?=9.99 rms uncertainty on Kmag  (2)
*Note (1): Positions with respect to the reference epoch,
*    1Ks07 for table7.dat (Field 1), 2Ks10 for table8.dat (Field 2)
*Note (2): Individual K-band magnitudes and the errors of the measurements
*    as given by ALLFRAME to which 0.04 mag calibration error have been
*    added in quadrature.
*    Magnitude 99.99 and the error 9.99 denote that no measurement for
*    that star was obtained.

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

C  Declarations for 'table8.dat'	! Individual measurements of the K-band magnitudes
                            for the 3 red variable stars detected in Field 2

      integer*4 nr__7
      parameter (nr__7=21672)	! Number of records
      character*41 ar__7  	! Full-size record

      character*8   v_RMS2003__6(nr__7) ! Identification, FN NNNN
      real*8        xpos_6     (nr__7) ! (pix) x position (G1) (1)
      real*8        ypos_6     (nr__7) ! (pix) y position (G1) (1)
      character*5   FilEp_7    (nr__7) ! Filter epoch
      real*4        Kmag_1     (nr__7) ! (mag) ?=99.99 Filter epoch K-band magnitude (2)
      real*4        e_Kmag_1   (nr__7) ! (mag) ?=9.99 rms uncertainty on Kmag  (2)
*Note (1): Positions with respect to the reference epoch,
*    1Ks07 for table7.dat (Field 1), 2Ks10 for table8.dat (Field 2)
*Note (2): Individual K-band magnitudes and the errors of the measurements
*    as given by ALLFRAME to which 0.04 mag calibration error have been
*    added in quadrature.
*    Magnitude 99.99 and the error 9.99 denote that no measurement for
*    that star was obtained.

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

C  Loading file 'table1.dat'	! Near-IR multi-epoch photometry observing log

C  Format for file interpretation

    1 format(A2,1X,F9.1,2X,I2,2X,F4.2,2X,F4.2,2X,A5,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,48
        read(1,'(A39)')ar__
        read(ar__,1)
     +  Field(i__),JD(i__),Exp(i__),X(i__),Seeing(i__),FilEp(i__),
     +  Instr(i__)
        if(ar__(37:37) .EQ. '') Instr(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Field(i__),JD(i__),Exp(i__),X(i__),Seeing(i__),FilEp(i__),
     +  Instr(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Data for a sample of variables with light curves
*                            plotted in Fig. 5

C  Format for file interpretation

    2 format(
     +  A8,1X,F7.3,1X,F7.3,1X,A5,1X,I4,1X,F5.3,1X,F6.2,1X,F5.3,1X,
     +  F5.2,1X,F5.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,40
        read(1,'(A72)')ar__1
        read(ar__1,2)
     +  v_RMS2003_(i__),xpos(i__),ypos(i__),FilEp_1(i__),Per(i__),
     +  A(i__),chi2(i__),signif(i__),Jsmag(i__),Hmag(i__),
     +  v_Ks_mag(i__)
        if(ar__1(32:35) .EQ. '') Per(i__) = iNULL__
        if(ar__1(37:41) .EQ. '') A(i__) = rNULL__
        if(ar__1(43:48) .EQ. '') chi2(i__) = rNULL__
        if(ar__1(50:54) .EQ. '') signif(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  v_RMS2003_(i__),xpos(i__),ypos(i__),FilEp_1(i__),Per(i__),
     +  A(i__),chi2(i__),signif(i__),Jsmag(i__),Hmag(i__),
     +  v_Ks_mag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! LPVs in Field 1

C  Format for file interpretation

    3 format(
     +  A8,1X,F7.3,1X,F7.3,1X,A5,1X,I4,1X,F5.3,1X,F6.2,1X,F5.3,1X,
     +  F5.2,1X,F5.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,437
        read(1,'(A72)')ar__2
        read(ar__2,3)
     +  v_RMS2003__1(i__),xpos_1(i__),ypos_1(i__),FilEp_2(i__),
     +  Per_1(i__),A_1(i__),chi2_1(i__),signif_1(i__),Jsmag_1(i__),
     +  Hmag_1(i__),v_Ks_mag_1(i__)
        if(ar__2(32:35) .EQ. '') Per_1(i__) = iNULL__
        if(ar__2(37:41) .EQ. '') A_1(i__) = rNULL__
        if(ar__2(43:48) .EQ. '') chi2_1(i__) = rNULL__
        if(ar__2(50:54) .EQ. '') signif_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  v_RMS2003__1(i__),xpos_1(i__),ypos_1(i__),FilEp_2(i__),
     +  Per_1(i__),A_1(i__),chi2_1(i__),signif_1(i__),Jsmag_1(i__),
     +  Hmag_1(i__),v_Ks_mag_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! LPVs in Field 2

C  Format for file interpretation

    4 format(
     +  A8,1X,F7.3,1X,F7.3,1X,A5,1X,I4,1X,F5.3,1X,F6.2,1X,F5.3,1X,
     +  F5.2,1X,F5.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,709
        read(1,'(A72)')ar__3
        read(ar__3,4)
     +  v_RMS2003__2(i__),xpos_2(i__),ypos_2(i__),FilEp_3(i__),
     +  Per_2(i__),A_2(i__),chi2_2(i__),signif_2(i__),Jsmag_2(i__),
     +  Hmag_2(i__),v_Ks_mag_2(i__)
        if(ar__3(32:35) .EQ. '') Per_2(i__) = iNULL__
        if(ar__3(37:41) .EQ. '') A_2(i__) = rNULL__
        if(ar__3(43:48) .EQ. '') chi2_2(i__) = rNULL__
        if(ar__3(50:54) .EQ. '') signif_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  v_RMS2003__2(i__),xpos_2(i__),ypos_2(i__),FilEp_3(i__),
     +  Per_2(i__),A_2(i__),chi2_2(i__),signif_2(i__),Jsmag_2(i__),
     +  Hmag_2(i__),v_Ks_mag_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Additional variable stars in Field 1 for which
*                            periods could not be determined

C  Format for file interpretation

    5 format(
     +  A8,1X,F7.3,1X,F7.3,1X,A5,1X,I4,1X,F5.3,1X,F6.2,1X,F5.3,1X,
     +  F5.2,1X,F5.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,163
        read(1,'(A72)')ar__4
        read(ar__4,5)
     +  v_RMS2003__3(i__),xpos_3(i__),ypos_3(i__),FilEp_4(i__),
     +  Per_3(i__),A_3(i__),chi2_3(i__),signif_3(i__),Jsmag_3(i__),
     +  Hmag_3(i__),v_Ks_mag_3(i__)
        if(ar__4(32:35) .EQ. '') Per_3(i__) = iNULL__
        if(ar__4(37:41) .EQ. '') A_3(i__) = rNULL__
        if(ar__4(43:48) .EQ. '') chi2_3(i__) = rNULL__
        if(ar__4(50:54) .EQ. '') signif_3(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  v_RMS2003__3(i__),xpos_3(i__),ypos_3(i__),FilEp_4(i__),
     +  Per_3(i__),A_3(i__),chi2_3(i__),signif_3(i__),Jsmag_3(i__),
     +  Hmag_3(i__),v_Ks_mag_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Additional variable stars in Field 2 for which
*                            periods could not be determined

C  Format for file interpretation

    6 format(
     +  A8,1X,F7.3,1X,F7.3,1X,A5,1X,I4,1X,F5.3,1X,F6.2,1X,F5.3,1X,
     +  F5.2,1X,F5.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,194
        read(1,'(A72)')ar__5
        read(ar__5,6)
     +  v_RMS2003__4(i__),xpos_4(i__),ypos_4(i__),FilEp_5(i__),
     +  Per_4(i__),A_4(i__),chi2_4(i__),signif_4(i__),Jsmag_4(i__),
     +  Hmag_4(i__),v_Ks_mag_4(i__)
        if(ar__5(32:35) .EQ. '') Per_4(i__) = iNULL__
        if(ar__5(37:41) .EQ. '') A_4(i__) = rNULL__
        if(ar__5(43:48) .EQ. '') chi2_4(i__) = rNULL__
        if(ar__5(50:54) .EQ. '') signif_4(i__) = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  v_RMS2003__4(i__),xpos_4(i__),ypos_4(i__),FilEp_5(i__),
     +  Per_4(i__),A_4(i__),chi2_4(i__),signif_4(i__),Jsmag_4(i__),
     +  Hmag_4(i__),v_Ks_mag_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Individual measurements of the K-band magnitudes
*                            for the 3 red variable stars detected in Field 1

C  Format for file interpretation

    7 format(A8,1X,F7.3,1X,F7.3,1X,A5,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,12020
        read(1,'(A41)')ar__6
        read(ar__6,7)
     +  v_RMS2003__5(i__),xpos_5(i__),ypos_5(i__),FilEp_6(i__),
     +  Kmag(i__),e_Kmag(i__)
c    ..............Just test output...........
        write(6,7)
     +  v_RMS2003__5(i__),xpos_5(i__),ypos_5(i__),FilEp_6(i__),
     +  Kmag(i__),e_Kmag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Individual measurements of the K-band magnitudes
*                            for the 3 red variable stars detected in Field 2

C  Format for file interpretation

    8 format(A8,1X,F7.3,1X,F7.3,1X,A5,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,21672
        read(1,'(A41)')ar__7
        read(ar__7,8)
     +  v_RMS2003__6(i__),xpos_6(i__),ypos_6(i__),FilEp_7(i__),
     +  Kmag_1(i__),e_Kmag_1(i__)
c    ..............Just test output...........
        write(6,8)
     +  v_RMS2003__6(i__),xpos_6(i__),ypos_6(i__),FilEp_7(i__),
     +  Kmag_1(i__),e_Kmag_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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