FORTRAN Generation
(/./ftp/cats/J/MNRAS/385/823)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/385/823 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-Oct-08
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/MNRAS/385/823     UBVRI light curves of OJ 287             (Zheng+, 2008)
*================================================================================
*Long-term optical spectra variability of BL Lacertae object OJ 287.
*    Zheng Y.G., Zhang X., Bi X.W., Hao J.M., Zhang H.J.
*   <Mon. Not. R. Astron. Soc., 385, 823-829 (2008)>
*   =2008MNRAS.385..823Z
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'tables1.dat'	! Compiled raw data of OJ 287 (0851+202),
                              magnitude values

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

      integer*4     JD         (nr__) ! (d) Julian date
      real*4        Umag       (nr__) ! (mag) U magnitude
      real*4        e_Umag     (nr__) ! (mag) ? rms uncertainty on Umag
      real*4        Bmag       (nr__) ! (mag) B magnitude
      real*4        e_Bmag     (nr__) ! (mag) ? rms uncertainty on Bmag
      real*4        Vmag       (nr__) ! (mag) V magnitude
      real*4        e_Vmag     (nr__) ! (mag) ? rms uncertainty on Vmag
      real*4        Rmag       (nr__) ! (mag) ? R magnitude
      real*4        e_Rmag     (nr__) ! (mag) ? rms uncertainty on Rmag
      real*4        Imag       (nr__) ! (mag) ? I magnitude
      real*4        e_Imag     (nr__) ! (mag) ? rms uncertainty on Imag
      integer*4     Ref        (nr__) ! Reference, in refs.dat file

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

C  Declarations for 'flux.dat'	! Compiled raw data of OJ 287 (0851+202),
                              flux density values

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

      integer*4     JD_1       (nr__1) ! (d) Julian date
      real*4        UFlux      (nr__1) ! (mJy) Flux density in U band
      real*4        e_UFlux    (nr__1) ! (mJy) ? rms uncertainty on UFlux
      real*4        BFlux      (nr__1) ! (mJy) Flux density in B band
      real*4        e_BFlux    (nr__1) ! (mJy) ? rms uncertainty on BFlux
      real*4        VFlux      (nr__1) ! (mJy) Flux density in V band
      real*4        e_VFlux    (nr__1) ! (mJy) ? rms uncertainty on VFlux
      real*4        RFlux      (nr__1) ! (mJy) ? Flux density in R band
      real*4        e_RFlux    (nr__1) ! (mJy) ? rms uncertainty on RFlux
      real*4        IFlux      (nr__1) ! (mJy) ? Flux density in I band
      real*4        e_IFlux    (nr__1) ! (mJy) ? rms uncertainty on IFlux
      integer*4     Ref_1      (nr__1) ! Reference, in refs.dat file

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

C  Declarations for 'refs.dat'	! References

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

      integer*4     Ref_2      (nr__2) ! Reference number
      character*19  BibCode    (nr__2) ! BibCode
      character*19  Aut        (nr__2) ! Author's name
      character*20  Com        (nr__2) ! Comments

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

C  Loading file 'tables1.dat'	! Compiled raw data of OJ 287 (0851+202),
*                              magnitude values

C  Format for file interpretation

    1 format(
     +  I7,1X,F5.2,1X,F4.2,2X,F5.2,1X,F4.2,2X,F5.2,1X,F4.2,2X,F5.2,1X,
     +  F4.2,2X,F5.2,1X,F4.2,3X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tables1.dat')
      write(6,*) '....Loading file: tables1.dat'
      do i__=1,256
        read(1,'(A71)')ar__
        read(ar__,1)
     +  JD(i__),Umag(i__),e_Umag(i__),Bmag(i__),e_Bmag(i__),Vmag(i__),
     +  e_Vmag(i__),Rmag(i__),e_Rmag(i__),Imag(i__),e_Imag(i__),
     +  Ref(i__)
        if(ar__(15:18) .EQ. '') e_Umag(i__) = rNULL__
        if(ar__(27:30) .EQ. '') e_Bmag(i__) = rNULL__
        if(ar__(39:42) .EQ. '') e_Vmag(i__) = rNULL__
        if(ar__(45:49) .EQ. '') Rmag(i__) = rNULL__
        if(ar__(51:54) .EQ. '') e_Rmag(i__) = rNULL__
        if(ar__(57:61) .EQ. '') Imag(i__) = rNULL__
        if(ar__(63:66) .EQ. '') e_Imag(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  JD(i__),Umag(i__),e_Umag(i__),Bmag(i__),e_Bmag(i__),Vmag(i__),
     +  e_Vmag(i__),Rmag(i__),e_Rmag(i__),Imag(i__),e_Imag(i__),
     +  Ref(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'flux.dat'	! Compiled raw data of OJ 287 (0851+202),
*                              flux density values

C  Format for file interpretation

    2 format(
     +  I7,1X,F5.2,1X,F4.2,2X,F5.2,1X,F4.2,2X,F5.2,1X,F4.2,2X,F5.2,1X,
     +  F4.2,2X,F5.2,1X,F4.2,3X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'flux.dat')
      write(6,*) '....Loading file: flux.dat'
      do i__=1,21
        read(1,'(A71)')ar__1
        read(ar__1,2)
     +  JD_1(i__),UFlux(i__),e_UFlux(i__),BFlux(i__),e_BFlux(i__),
     +  VFlux(i__),e_VFlux(i__),RFlux(i__),e_RFlux(i__),IFlux(i__),
     +  e_IFlux(i__),Ref_1(i__)
        if(ar__1(15:18) .EQ. '') e_UFlux(i__) = rNULL__
        if(ar__1(27:30) .EQ. '') e_BFlux(i__) = rNULL__
        if(ar__1(39:42) .EQ. '') e_VFlux(i__) = rNULL__
        if(ar__1(45:49) .EQ. '') RFlux(i__) = rNULL__
        if(ar__1(51:54) .EQ. '') e_RFlux(i__) = rNULL__
        if(ar__1(57:61) .EQ. '') IFlux(i__) = rNULL__
        if(ar__1(63:66) .EQ. '') e_IFlux(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  JD_1(i__),UFlux(i__),e_UFlux(i__),BFlux(i__),e_BFlux(i__),
     +  VFlux(i__),e_VFlux(i__),RFlux(i__),e_RFlux(i__),IFlux(i__),
     +  e_IFlux(i__),Ref_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    3 format(I2,1X,A19,1X,A19,6X,A20)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,28
        read(1,'(A68)')ar__2
        read(ar__2,3)Ref_2(i__),BibCode(i__),Aut(i__),Com(i__)
c    ..............Just test output...........
        write(6,3)Ref_2(i__),BibCode(i__),Aut(i__),Com(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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