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 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 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       ! Field (F1 or F2)
      real*8        JD          ! (d) Julian date
      integer*4     Exp         ! (min) Exposure time
      real*4        X           ! Mean airmass during the observing sequence
      real*4        Seeing      ! (arcsec) Seeing
      character*5   FilEp       ! Filter epoch
      integer*4     Instr       ! [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_  ! Identification, FN NNNNN
      real*8        xpos        ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos        ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_1     ! Reference epoch, defined in table1.dat
      integer*4     Per         ! (d) ? Best fitting period
      real*4        A           ! (mag) ? Semi-amplitude
      real*4        chi2        ! ? Reduced chi2 of the sine-curve fit
      real*4        signif      ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag       ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag        ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag    ! (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 ! Identification, FN NNNNN
      real*8        xpos_1      ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos_1      ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_2     ! Reference epoch, defined in table1.dat
      integer*4     Per_1       ! (d) ? Best fitting period
      real*4        A_1         ! (mag) ? Semi-amplitude
      real*4        chi2_1      ! ? Reduced chi2 of the sine-curve fit
      real*4        signif_1    ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag_1     ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag_1      ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag_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 '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 ! Identification, FN NNNNN
      real*8        xpos_2      ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos_2      ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_3     ! Reference epoch, defined in table1.dat
      integer*4     Per_2       ! (d) ? Best fitting period
      real*4        A_2         ! (mag) ? Semi-amplitude
      real*4        chi2_2      ! ? Reduced chi2 of the sine-curve fit
      real*4        signif_2    ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag_2     ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag_2      ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag_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 '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 ! Identification, FN NNNNN
      real*8        xpos_3      ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos_3      ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_4     ! Reference epoch, defined in table1.dat
      integer*4     Per_3       ! (d) ? Best fitting period
      real*4        A_3         ! (mag) ? Semi-amplitude
      real*4        chi2_3      ! ? Reduced chi2 of the sine-curve fit
      real*4        signif_3    ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag_3     ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag_3      ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag_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 '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 ! Identification, FN NNNNN
      real*8        xpos_4      ! (pix) x position with respect to the
*                                 reference epoch (G1)
      real*8        ypos_4      ! (pix) y position with respect to the
*                                 reference epoch (G1)
      character*5   FilEp_5     ! Reference epoch, defined in table1.dat
      integer*4     Per_4       ! (d) ? Best fitting period
      real*4        A_4         ! (mag) ? Semi-amplitude
      real*4        chi2_4      ! ? Reduced chi2 of the sine-curve fit
      real*4        signif_4    ! ? Significance of the period from
*                                   Fourier analysis
      real*4        Jsmag_4     ! (mag) ?=99.99 Single epoch Js-band magnitude (1)
      real*4        Hmag_4      ! (mag) ?=99.99 Single epoch H-band magnitude (1)
      real*4        v_Ks_mag_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 '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 ! Identification, FN NNNN
      real*8        xpos_5      ! (pix) x position (G1) (1)
      real*8        ypos_5      ! (pix) y position (G1) (1)
      character*5   FilEp_6     ! Filter epoch
      real*4        Kmag        ! (mag) ?=99.99 Filter epoch K-band magnitude (2)
      real*4        e_Kmag      ! (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 ! Identification, FN NNNN
      real*8        xpos_6      ! (pix) x position (G1) (1)
      real*8        ypos_6      ! (pix) y position (G1) (1)
      character*5   FilEp_7     ! Filter epoch
      real*4        Kmag_1      ! (mag) ?=99.99 Filter epoch K-band magnitude (2)
      real*4        e_Kmag_1    ! (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,JD,Exp,X,Seeing,FilEp,Instr
        if(ar__(37:37) .EQ. '') Instr = iNULL__
c    ..............Just test output...........
        write(6,1)Field,JD,Exp,X,Seeing,FilEp,Instr
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_,xpos,ypos,FilEp_1,Per,A,chi2,signif,Jsmag,Hmag,
     +  v_Ks_mag
        if(ar__1(32:35) .EQ. '') Per = iNULL__
        if(ar__1(37:41) .EQ. '') A = rNULL__
        if(ar__1(43:48) .EQ. '') chi2 = rNULL__
        if(ar__1(50:54) .EQ. '') signif = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  v_RMS2003_,xpos,ypos,FilEp_1,Per,A,chi2,signif,Jsmag,Hmag,
     +  v_Ks_mag
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,xpos_1,ypos_1,FilEp_2,Per_1,A_1,chi2_1,signif_1,
     +  Jsmag_1,Hmag_1,v_Ks_mag_1
        if(ar__2(32:35) .EQ. '') Per_1 = iNULL__
        if(ar__2(37:41) .EQ. '') A_1 = rNULL__
        if(ar__2(43:48) .EQ. '') chi2_1 = rNULL__
        if(ar__2(50:54) .EQ. '') signif_1 = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  v_RMS2003__1,xpos_1,ypos_1,FilEp_2,Per_1,A_1,chi2_1,signif_1,
     +  Jsmag_1,Hmag_1,v_Ks_mag_1
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,xpos_2,ypos_2,FilEp_3,Per_2,A_2,chi2_2,signif_2,
     +  Jsmag_2,Hmag_2,v_Ks_mag_2
        if(ar__3(32:35) .EQ. '') Per_2 = iNULL__
        if(ar__3(37:41) .EQ. '') A_2 = rNULL__
        if(ar__3(43:48) .EQ. '') chi2_2 = rNULL__
        if(ar__3(50:54) .EQ. '') signif_2 = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  v_RMS2003__2,xpos_2,ypos_2,FilEp_3,Per_2,A_2,chi2_2,signif_2,
     +  Jsmag_2,Hmag_2,v_Ks_mag_2
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,xpos_3,ypos_3,FilEp_4,Per_3,A_3,chi2_3,signif_3,
     +  Jsmag_3,Hmag_3,v_Ks_mag_3
        if(ar__4(32:35) .EQ. '') Per_3 = iNULL__
        if(ar__4(37:41) .EQ. '') A_3 = rNULL__
        if(ar__4(43:48) .EQ. '') chi2_3 = rNULL__
        if(ar__4(50:54) .EQ. '') signif_3 = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  v_RMS2003__3,xpos_3,ypos_3,FilEp_4,Per_3,A_3,chi2_3,signif_3,
     +  Jsmag_3,Hmag_3,v_Ks_mag_3
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,xpos_4,ypos_4,FilEp_5,Per_4,A_4,chi2_4,signif_4,
     +  Jsmag_4,Hmag_4,v_Ks_mag_4
        if(ar__5(32:35) .EQ. '') Per_4 = iNULL__
        if(ar__5(37:41) .EQ. '') A_4 = rNULL__
        if(ar__5(43:48) .EQ. '') chi2_4 = rNULL__
        if(ar__5(50:54) .EQ. '') signif_4 = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  v_RMS2003__4,xpos_4,ypos_4,FilEp_5,Per_4,A_4,chi2_4,signif_4,
     +  Jsmag_4,Hmag_4,v_Ks_mag_4
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,xpos_5,ypos_5,FilEp_6,Kmag,e_Kmag
c    ..............Just test output...........
        write(6,7)v_RMS2003__5,xpos_5,ypos_5,FilEp_6,Kmag,e_Kmag
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,xpos_6,ypos_6,FilEp_7,Kmag_1,e_Kmag_1
c    ..............Just test output...........
        write(6,8)
     +  v_RMS2003__6,xpos_6,ypos_6,FilEp_7,Kmag_1,e_Kmag_1
c    .......End.of.Just test output...........
      end do
      close(1)

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