FORTRAN Generation
(/./ftp/cats/J/ApJ/459/278)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/459/278 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.5, on 2013-Jun-19
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__=-1.e37)     	! NULL real number
      parameter  (iNULL__=-2147483647)	! NULL int  number

C=============================================================================
Cat. J/ApJ/459/278            Photometry of HD 3831 (Matthews+ 1996)
*================================================================================
*Opening a new window on Ap star atmospheres: a T-Tauri relation for HR 3831
*from its limb-darkened pulsation amplitudes
*     Matthews J.M., Wehlau W.H., Rice J., Walker G.A.H.
*    <Astrophys. J. 459, 278 (1996)>
*    =1996ApJ...459..278M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'hr3831b.dat'	! Differential photometry, b bandpass

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

      real*8        HJD_2448000(nr__) ! (d) Heliocentric Julian Date - 2448000
      real*8        dMag       (nr__) ! (mag) Differential magnitude normalized to the mean

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

C  Declarations for 'hr3831v.dat'	! Differential photometry, v bandpass

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

      real*8        HJD_2448000_1(nr__1) ! (d) Heliocentric Julian Date - 2448000
      real*8        dMag_1     (nr__1) ! (mag) Differential magnitude normalized to the mean

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

C  Declarations for 'hr3831y.dat'	! Differential photometry, y bandpass

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

      real*8        HJD_2448000_2(nr__2) ! (d) Heliocentric Julian Date - 2448000
      real*8        dMag_2     (nr__2) ! (mag) Differential magnitude normalized to the mean

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

C  Declarations for 'hr3831r.dat'	! Differential photometry, R bandpass

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

      real*8        HJD_2448000_3(nr__3) ! (d) Heliocentric Julian Date - 2448000
      real*8        dMag_3     (nr__3) ! (mag) Differential magnitude normalized to the mean

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

C  Declarations for 'hr3831i.dat'	! Differential photometry, I bandpass

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

      real*8        HJD_2448000_4(nr__4) ! (d) Heliocentric Julian Date - 2448000
      real*8        dMag_4     (nr__4) ! (mag) Differential magnitude normalized to the mean

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

C  Declarations for 'hr3831j.dat'	! Differential photometry, J bandpass

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

      real*8        HJD_2448000_5(nr__5) ! (d) Heliocentric Julian Date - 2448000
      real*8        dMag_5     (nr__5) ! (mag) Differential magnitude normalized to the mean

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

C  Declarations for 'hr3831h.dat'	! Differential photometry, H bandpass

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

      real*8        HJD_2448000_6(nr__6) ! (d) Heliocentric Julian Date - 2448000
      real*8        dMag_6     (nr__6) ! (mag) Differential magnitude normalized to the mean

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

C  Declarations for 'hr3831k.dat'	! Differential photometry, K bandpass

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

      real*8        HJD_2448000_7(nr__7) ! (d) Heliocentric Julian Date - 2448000
      real*8        dMag_7     (nr__7) ! (mag) Differential magnitude normalized to the mean

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

C  Loading file 'hr3831b.dat'	! Differential photometry, b bandpass

C  Format for file interpretation

    1 format(F9.5,1X,F7.4)

C  Effective file loading

      open(unit=1,file='hr3831b.dat', status='old')
      write(6,*) '....Loading file: hr3831b.dat'
      do i__=1,151
        read(1,'(A17)')ar__
        read(ar__,1)HJD_2448000(i__),dMag(i__)
c    ..............Just test output...........
        write(6,1)HJD_2448000(i__),dMag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hr3831v.dat'	! Differential photometry, v bandpass

C  Format for file interpretation

    2 format(F9.5,1X,F7.4)

C  Effective file loading

      open(unit=1,file='hr3831v.dat', status='old')
      write(6,*) '....Loading file: hr3831v.dat'
      do i__=1,152
        read(1,'(A17)')ar__1
        read(ar__1,2)HJD_2448000_1(i__),dMag_1(i__)
c    ..............Just test output...........
        write(6,2)HJD_2448000_1(i__),dMag_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hr3831y.dat'	! Differential photometry, y bandpass

C  Format for file interpretation

    3 format(F9.5,1X,F7.4)

C  Effective file loading

      open(unit=1,file='hr3831y.dat', status='old')
      write(6,*) '....Loading file: hr3831y.dat'
      do i__=1,152
        read(1,'(A17)')ar__2
        read(ar__2,3)HJD_2448000_2(i__),dMag_2(i__)
c    ..............Just test output...........
        write(6,3)HJD_2448000_2(i__),dMag_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hr3831r.dat'	! Differential photometry, R bandpass

C  Format for file interpretation

    4 format(F9.5,1X,F7.4)

C  Effective file loading

      open(unit=1,file='hr3831r.dat', status='old')
      write(6,*) '....Loading file: hr3831r.dat'
      do i__=1,144
        read(1,'(A17)')ar__3
        read(ar__3,4)HJD_2448000_3(i__),dMag_3(i__)
c    ..............Just test output...........
        write(6,4)HJD_2448000_3(i__),dMag_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hr3831i.dat'	! Differential photometry, I bandpass

C  Format for file interpretation

    5 format(F9.5,1X,F7.4)

C  Effective file loading

      open(unit=1,file='hr3831i.dat', status='old')
      write(6,*) '....Loading file: hr3831i.dat'
      do i__=1,141
        read(1,'(A17)')ar__4
        read(ar__4,5)HJD_2448000_4(i__),dMag_4(i__)
c    ..............Just test output...........
        write(6,5)HJD_2448000_4(i__),dMag_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hr3831j.dat'	! Differential photometry, J bandpass

C  Format for file interpretation

    6 format(F9.5,1X,F7.4)

C  Effective file loading

      open(unit=1,file='hr3831j.dat', status='old')
      write(6,*) '....Loading file: hr3831j.dat'
      do i__=1,91
        read(1,'(A17)')ar__5
        read(ar__5,6)HJD_2448000_5(i__),dMag_5(i__)
c    ..............Just test output...........
        write(6,6)HJD_2448000_5(i__),dMag_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hr3831h.dat'	! Differential photometry, H bandpass

C  Format for file interpretation

    7 format(F9.5,1X,F7.4)

C  Effective file loading

      open(unit=1,file='hr3831h.dat', status='old')
      write(6,*) '....Loading file: hr3831h.dat'
      do i__=1,92
        read(1,'(A17)')ar__6
        read(ar__6,7)HJD_2448000_6(i__),dMag_6(i__)
c    ..............Just test output...........
        write(6,7)HJD_2448000_6(i__),dMag_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hr3831k.dat'	! Differential photometry, K bandpass

C  Format for file interpretation

    8 format(F9.5,1X,F7.4)

C  Effective file loading

      open(unit=1,file='hr3831k.dat', status='old')
      write(6,*) '....Loading file: hr3831k.dat'
      do i__=1,92
        read(1,'(A17)')ar__7
        read(ar__7,8)HJD_2448000_7(i__),dMag_7(i__)
c    ..............Just test output...........
        write(6,8)HJD_2448000_7(i__),dMag_7(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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