FORTRAN Generation
(/./ftp/cats/J/ApJ/700/209)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/700/209 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 2024-Mar-29
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/ApJ/700/209   Spectroscopy of the B[e] star VES 695   (Miroshnichenko+, 2009)
*================================================================================
*Toward understanding the B[e] phenomenon.
*III. Properties of the optical counterpart of IRAS 00470+6429.
*    Miroshnichenko A.S., Chentsov E.L., Klochkova V.G., Zharikov S.V.,
*    Grankin K.N., Kusakin A.V., Gandet T.L., Klingenberg G., Kildahl S.,
*    Rudy R.J., Lynch D.K., Venturini C.C., Mazuk S., Puetter R.C., Perry R.B.,
*    Carciofi A.C., Bjorkman K.S., Gray R.O., Bernabei S., Polcaro V.F.,
*    Viotti R.F., Norci L.
*   <Astrophys. J., 700, 209-220 (2009)>
*   =2009ApJ...700..209M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Line identification, intensities and radial
                             velocities V_r_ in the optical spectrum of
                             IRAS 00470+6429

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

      character*9   Line        ! Line identification (1)
      real*8        lambda      ! (0.1nm) ? Laboratory wavelength in Angstroms
*--------------------------------------------------------------------------------
      real*4        I030814a    ! ?=- 2003-08-14 peak relative intensity (2)
      real*4        I030814e    ! ?=- 2003-08-14 peak relative intensity (2)
      integer*4     V030814a    ! (km/s) ?=- 2003-08-14 heliocentric radial velocity (2)
      character*1   u_V030814a  ! Uncertainty flag (:) on V030814a
      integer*4     V030814e    ! (km/s) ?=- 2003-08-14 heliocentric radial velocity (2)
      character*1   u_V030814e  ! Uncertainty flag (:) on V030814e
*--------------------------------------------------------------------------------
      real*4        I040309a    ! ?=- 2004-03-09 peak relative intensity (2)
      character*1   u_I040309a  ! Uncertainty flag (:) on I040309a
      real*4        I040309e    ! ?=- 2004-03-09 peak relative intensity (2)
      character*1   u_I040309e  ! Uncertainty flag (:) on I040309e
      integer*4     V040309a    ! (km/s) ?=- 2004-03-09 heliocentric radial velocity (2)
      character*1   u_V040309a  ! Uncertainty flag (:) on V040309a
      integer*4     V040309e    ! (km/s) ?=- 2004-03-09 heliocentric radial velocity (2)
      character*1   u_V040309e  ! Uncertainty flag (:) on V040309e
*--------------------------------------------------------------------------------
      real*4        I041006a    ! ?=- 2004-10-06 peak relative intensity (2)
      real*4        I041006e    ! ?=- 2004-10-06 peak relative intensity (2)
      integer*4     V041006a    ! (km/s) ?=- 2004-10-06 heliocentric radial velocity (2)
      character*1   u_V041006a  ! Uncertainty flag (:) on V041006a
      integer*4     V041006e    ! (km/s) ?=- 2004-10-06 heliocentric radial velocity (2)
      character*1   u_V041006e  ! Uncertainty flag (:) on V041006e
*--------------------------------------------------------------------------------
      real*4        I051111a    ! ?=- 2005-11-11 peak relative intensity (2)
      character*1   u_I051111a  ! Uncertainty flag (:) on I051111a
      real*4        I051111e    ! ?=- 2005-11-11 peak relative intensity (2)
      character*1   u_I051111e  ! Uncertainty flag (:) on I051111e
      integer*4     V051111a    ! (km/s) ?=- 2005-11-11 heliocentric radial velocity (2)
      character*1   u_V051111a  ! Uncertainty flag (:) on V051111a
      integer*4     V051111e    ! (km/s) ?=- 2005-11-11 heliocentric radial velocity (2)
      character*1   u_V051111e  ! Uncertainty flag (:) on V051111e
*--------------------------------------------------------------------------------
      real*4        I051114a    ! ?=- 2005-11-14 peak relative intensity (2)
      character*1   u_I051114a  ! Uncertainty flag (:) on I051114a
      real*4        I051114e    ! ?=- 2005-11-14 peak relative intensity (2)
      character*1   u_I051114e  ! Uncertainty flag (:) on I051114e
      integer*4     V051114a    ! (km/s) ?=- 2005-11-14 heliocentric radial velocity (2)
      character*1   u_V051114a  ! Uncertainty flag (:) on V051114a
      integer*4     V051114e    ! (km/s) ?=- 2005-11-14 heliocentric radial velocity (2)
      character*1   u_V051114e  ! Uncertainty flag (:) on V051114e
*--------------------------------------------------------------------------------
      real*4        I051218a    ! ?=- 2005-12-18 peak relative intensity (2)
      real*4        I051218e    ! ?=- 2005-12-18 peak relative intensity (2)
      integer*4     V051218a    ! (km/s) ?=- 2005-12-18 heliocentric radial velocity (2)
      character*1   u_V051218a  ! Uncertainty flag (:) on V051218a
      integer*4     V051218e    ! (km/s) ?=- 2005-12-18 heliocentric radial velocity (2)
      character*1   u_V051218e  ! Uncertainty flag (:) on V051218e
*--------------------------------------------------------------------------------
      real*4        I061227a    ! ?=- 2006-12-27 peak relative intensity (2)
      character*1   u_I061227a  ! Uncertainty flag (:) on I061227a
      real*4        I061227e    ! ?=- 2006-12-27 peak relative intensity (2)
      character*1   u_I061227e  ! Uncertainty flag (:) on I061227e
      integer*4     V061227a    ! (km/s) ?=- 2006-12-27 heliocentric radial velocity (2)
      character*1   u_V061227a  ! Uncertainty flag (:) on V061227a
      integer*4     V061227e    ! (km/s) ?=- 2006-12-27 heliocentric radial velocity (2)
      character*1   u_V061227e  ! Uncertainty flag (:) on V061227e
*--------------------------------------------------------------------------------
      real*4        I081214a    ! ?=- 2008-12-14 peak relative intensity (2)
      character*1   u_I081214a  ! Uncertainty flag (:) on I081214a
      real*4        I081214e    ! ?=- 2008-12-14 peak relative intensity (2)
      character*1   u_I081214e  ! Uncertainty flag (:) on I081214e
      integer*4     V081214a    ! (km/s) ?=- 2008-12-14 heliocentric radial velocity (2)
      character*1   u_V081214a  ! Uncertainty flag (:) on V081214a
      integer*4     V081214e    ! (km/s) ?=- 2008-12-14 heliocentric radial velocity (2)
      character*1   u_V081214e  ! Uncertainty flag (:) on V081214e
*Note (1): Multiplet number is given in brackets.
*Note (2): Columns contain residual intensities of the absorption / emission
*     peak in continuum units (r) and corresponding heliocentric radial
*     velocities for the specified date. Dashes means that the line
*     component is either not recognized due to noise or its parameters can
*     hardly be measured. If parameters of a line were hard to measure in a
*     particular spectrum, no information about this line was entered.

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

C  Loading file 'table2.dat'	! Line identification, intensities and radial
*                             velocities V_r_ in the optical spectrum of
*                             IRAS 00470+6429

C  Format for file interpretation

    1 format(
     +  A9,1X,F8.2,2X,F4.2,2X,F5.2,3X,I4,A1,2X,I4,A1,1X,F5.2,A1,1X,
     +  F5.2,A1,2X,I4,A1,2X,I4,A1,1X,F5.2,2X,F5.2,3X,I4,A1,2X,I4,A1,
     +  1X,F5.2,A1,1X,F5.2,A1,2X,I4,A1,2X,I4,A1,1X,F5.2,A1,1X,F5.2,A1,
     +  2X,I4,A1,2X,I4,A1,1X,F5.2,2X,F5.2,3X,I4,A1,2X,I4,A1,1X,F5.2,
     +  A1,1X,F5.2,A1,2X,I4,A1,2X,I4,A1,1X,F5.2,A1,1X,F5.2,A1,2X,I4,
     +  A1,2X,I4,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,161
        read(1,'(A242)')ar__
        read(ar__,1)
     +  Line,lambda,I030814a,I030814e,V030814a,u_V030814a,V030814e,
     +  u_V030814e,I040309a,u_I040309a,I040309e,u_I040309e,V040309a,
     +  u_V040309a,V040309e,u_V040309e,I041006a,I041006e,V041006a,
     +  u_V041006a,V041006e,u_V041006e,I051111a,u_I051111a,I051111e,
     +  u_I051111e,V051111a,u_V051111a,V051111e,u_V051111e,I051114a,
     +  u_I051114a,I051114e,u_I051114e,V051114a,u_V051114a,V051114e,
     +  u_V051114e,I051218a,I051218e,V051218a,u_V051218a,V051218e,
     +  u_V051218e,I061227a,u_I061227a,I061227e,u_I061227e,V061227a,
     +  u_V061227a,V061227e,u_V061227e,I081214a,u_I081214a,I081214e,
     +  u_I081214e,V081214a,u_V081214a,V081214e,u_V081214e
        if(ar__(11:18) .EQ. '') lambda = rNULL__
        if (idig(ar__(21:24)).EQ.0) I030814a =  rNULL__
        if (idig(ar__(27:31)).EQ.0) I030814e =  rNULL__
        if (V030814a .EQ. 45) V030814a =  iNULL__
        if (V030814e .EQ. 45) V030814e =  iNULL__
        if (idig(ar__(48:52)).EQ.0) I040309a =  rNULL__
        if (idig(ar__(55:59)).EQ.0) I040309e =  rNULL__
        if (V040309a .EQ. 45) V040309a =  iNULL__
        if (V040309e .EQ. 45) V040309e =  iNULL__
        if (idig(ar__(76:80)).EQ.0) I041006a =  rNULL__
        if (idig(ar__(83:87)).EQ.0) I041006e =  rNULL__
        if (V041006a .EQ. 45) V041006a =  iNULL__
        if (V041006e .EQ. 45) V041006e =  iNULL__
        if (idig(ar__(104:108)).EQ.0) I051111a =  rNULL__
        if (idig(ar__(111:115)).EQ.0) I051111e =  rNULL__
        if (V051111a .EQ. 45) V051111a =  iNULL__
        if (V051111e .EQ. 45) V051111e =  iNULL__
        if (idig(ar__(132:136)).EQ.0) I051114a =  rNULL__
        if (idig(ar__(139:143)).EQ.0) I051114e =  rNULL__
        if (V051114a .EQ. 45) V051114a =  iNULL__
        if (V051114e .EQ. 45) V051114e =  iNULL__
        if (idig(ar__(160:164)).EQ.0) I051218a =  rNULL__
        if (idig(ar__(167:171)).EQ.0) I051218e =  rNULL__
        if (V051218a .EQ. 45) V051218a =  iNULL__
        if (V051218e .EQ. 45) V051218e =  iNULL__
        if (idig(ar__(188:192)).EQ.0) I061227a =  rNULL__
        if (idig(ar__(195:199)).EQ.0) I061227e =  rNULL__
        if (V061227a .EQ. 45) V061227a =  iNULL__
        if (V061227e .EQ. 45) V061227e =  iNULL__
        if (idig(ar__(216:220)).EQ.0) I081214a =  rNULL__
        if (idig(ar__(223:227)).EQ.0) I081214e =  rNULL__
        if (V081214a .EQ. 45) V081214a =  iNULL__
        if (V081214e .EQ. 45) V081214e =  iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Line,lambda,I030814a,I030814e,V030814a,u_V030814a,V030814e,
     +  u_V030814e,I040309a,u_I040309a,I040309e,u_I040309e,V040309a,
     +  u_V040309a,V040309e,u_V040309e,I041006a,I041006e,V041006a,
     +  u_V041006a,V041006e,u_V041006e,I051111a,u_I051111a,I051111e,
     +  u_I051111e,V051111a,u_V051111a,V051111e,u_V051111e,I051114a,
     +  u_I051114a,I051114e,u_I051114e,V051114a,u_V051114a,V051114e,
     +  u_V051114e,I051218a,I051218e,V051218a,u_V051218a,V051218e,
     +  u_V051218e,I061227a,u_I061227a,I061227e,u_I061227e,V061227a,
     +  u_V061227a,V061227e,u_V061227e,I081214a,u_I081214a,I081214e,
     +  u_I081214e,V081214a,u_V081214a,V081214e,u_V081214e
c    .......End.of.Just test output...........
      end do
      close(1)

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

C Locate position of first digit in string; or return 0
      integer function idig(c)
      character*(*) c
      character*1 c1
      integer lc,i
      lc=len(c)
      idig=0
      do i=1,lc
         if(c(i:i).ne.' ') go to 1
      end do
    1 if(i.gt.lc) return
      c1=c(i:i)
      if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
      if(i.gt.lc) return
      c1=c(i:i)
      if(c1.ge.'0'.and.c1.le.'9') idig=i
      return
      end