FORTRAN Generation
()

Conversion of standardized ReadMe file for catalog 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.5, on 2013-Jun-18
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__
ter  (rNULL__=-1.e37)     	! NULL real number
      parameter  (iNULL__=-2147483647)	! NULL int  number

C=============================================================================
Cat. J/other/Sci/292.2290    Chandra compact binaries in 47 Tuc     (Grindlay+, 2001)
*================================================================================
*High-resolution X-ray imaging of a globular cluster core:
*compact binaries in 47 Tuc
*    Grindlay J.E., Heinke C., Edmonds P.D., Murray S.S.
*   <Science, 292, 2290 (2001)>
*   =2001Sci...292.2290G
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! 47 Tuc X-Ray Source Parameters (MS 61135)

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

      real*8        RA_         ! (deg) Right Ascension J2000
      real*8        Dec         ! (deg)     Declination J2000
      character*4   v_GHE2001_  ! Chandra source number (WNNN)
      integer*4     RAh         ! (h) Right ascension (J2000)
      integer*4     RAm         ! (min) Right ascension (J2000)
      real*4        RAs         ! (s) Right ascension (J2000)
      real*4        e_RAs       ! (s) ? rms uncertainty on RA (1)
      character*1   DE_         ! Declination sign (J2000)
      integer*4     DEd         ! (deg) Declination (J2000)
      integer*4     DEm         ! (arcmin) Declination (J2000)
      real*4        DEs         ! (arcsec) DEclination (J2000)
      real*4        e_DEs       ! (arcsec) ? rms uncertainty on DE (1)
      real*8        Medcts      ! (ct) Counts in 0.5-4.5 keV band (2)
      real*4        Sfcts       ! (ct) Counts in 0.5-1.5 keV band (2)
      real*4        Hdcts       ! (ct) Counts in 1.5-6 keV band (2)
      real*4        Sf_Hd       ! Softness ratio, defined as Sfcts/Hdcts
      real*4        LX          ! (10-7W) Intrinsic luminosities in ROSAT (0.5-2.5keV)
*                                    band (3)
      character*10  ID          ! Identification (4)
*Note (1): Errors are formal 1{sigma} centroiding errors from Chandra
*    X-ray detection routine wavdetect. The final five millisecond pulsars
*    were not formally detected by wavdetect, so the radio positions of
*    Freire et al., 2001MNRAS.326..901F are given without errors.
*Note (2): Values of 0.5 counts are estimated upper limits.
*Note (3): Intrinsic luminosities in ROSAT (0.5-2.5 keV) band, with
*    multiplicative exposure corrections (=1.4, 2.0, 1.7, 2.6, 1.2, 2.5 and
*    1.9 for sources W63(MSPJ), W65, W88, W89, W90, MSPC and MSPM,
*    respectively), based on 1keV thermal bremsstrahlung
*    spectrum,N_H_=2.4x10^20^cm^-2^, d=5.0kpc
*Note (4): X number from ROSAT (Verbunt & Hasinger, 1998A&A...336..895V)
*    observations (e.g. X5), certain optical or radio counterparts (e.g.
*    V1, AKO9, MSP-D), and class identifications:
*      qLMXB=quiescent LMXB,
*        M-S=main sequence binary,
*         CV=cataclysmic variable.
*    Detailed discussion of individual X-ray source spectra and variability
*    is deferred to Heinke et al. (2003ApJ...588..452H) and probable
*    optical counterparts are discussed by Edmonds et al.
*    (2002ApJ...579..741E)

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

C  Loading file 'table1.dat'	! 47 Tuc X-Ray Source Parameters (MS 61135)

C  Format for file interpretation

    1 format(
     +  A4,1X,I1,1X,I2,1X,F6.3,1X,F5.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F4.2,
     +  1X,F7.1,1X,F6.1,1X,F6.1,1X,F5.2,1X,E6.2,1X,A10)

C  Effective file loading

      open(unit=1,file='table1.dat', status='old')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,108
        read(1,'(A86)')ar__
        read(ar__,1)
     +  v_GHE2001_,RAh,RAm,RAs,e_RAs,DE_,DEd,DEm,DEs,e_DEs,Medcts,
     +  Sfcts,Hdcts,Sf_Hd,LX,ID
        if(ar__(18:22) .EQ. '') e_RAs = rNULL__
        if(ar__(37:40) .EQ. '') e_DEs = rNULL__
c       Derive coordinates RA_ and Dec from input data
c       (RA_ and Dec are set to rNULL_ when unknown)
        RA_ = RAh
        if(RAh .GE. 0) RA_=RAh*15.
        if(RAm .GE. 0) RA_=RA_+RAm/4.
        if(RAs .GE. 0) RA_=RA_+RAs/240.
        Dec = DEd
        if(DEm .GE. 0) Dec=Dec+DEm/60.
        if(DEs .GE. 0) Dec=Dec+DEs/3600.
        if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c    ..............Just test output...........
        write(6,1)
     +  v_GHE2001_,RAh,RAm,RAs,e_RAs,DE_,DEd,DEm,DEs,e_DEs,Medcts,
     +  Sfcts,Hdcts,Sf_Hd,LX,ID
        write(6,'(6H Pos: 2F8.4)') RA_,Dec
c    .......End.of.Just test output...........
      end do
      close(1)

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