FORTRAN Generation
(/./ftp/cats/J/AJ/133/791)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/133/791 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 2024-Mar-28
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/AJ/133/791 Spitzer photometry of normal & interacting galaxies (Smith+, 2007)
*================================================================================
*The Spitzer spirals, bridges, and tails interacting galaxy survey:
*interaction-induced star formation in the mid-infrared.
*    Smith B.J., Struck C., Hancock M., Appleton P.N., Charmandaris V.,
*    Reach W.T.
*   <Astron. J., 133, 791-817 (2007)>
*   =2007AJ....133..791S
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Interacting Galaxy Sample

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

      character*8   Arp        (nr__) ! Arp designation of the system
      character*17  OName      (nr__) ! Other Name
      real*4        Dist       (nr__) ! (Mpc) Distance of the system
      real*4        Sep        (nr__) ! (kpc) Separation
      real*4        logL_FIR   (nr__) ! ([solLum]) ?=- FIR luminosity (1)
      character*51  Notes      (nr__) ! Notes on morphology
      character*27  SpType     (nr__) ! Nuclear Spectral Class
      character*3   Flags      (nr__) ! Individual notes (2)
*Note (1): Includes the IRAS flux from both galaxies in the pair.
*Note (2): Flags as follows:
*      b = Nuclear spectral type from Dahari (1985ApJS...57..643D).
*      c = Nuclear spectral type from Keel et al. (1985AJ.....90..708K).
*      d = The bright galaxy in the Arp (1966, Cat. <VII/74>) Atlas
*          photograph of Arp 136 is NGC 5820. This galaxy has a nearby
*          similar-mass companion, NGC 5821, at the same redshift, outside
*          of the Arp Atlas field of view. The second small galaxy in the
*          Arp Atlas image is a background galaxy.

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

C  Declarations for 'table2.dat'	! Normal Galaxy Sample

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

      character*1   Type       (nr__1) ! [SEI] Type of galaxy (G1)
      character*11  Name       (nr__1) ! Galaxy name
      character*5   MType      (nr__1) ! Morphological Type
      real*4        Dist_1     (nr__1) ! (Mpc) Distance
      character*1   l_logL_IR  (nr__1) ! Limit flag on logL(IR)
      real*4        logL_IR    (nr__1) ! ([solLum]) Infrared luminosity

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

C  Declarations for 'table5.dat'	! Spitzer Flux Densities for Interacting Galaxy Sample

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

      character*8   Arp_1      (nr__2) ! The Arp designation
      character*9   Comp       (nr__2) ! Component
      character*23  OName_1    (nr__2) ! Other name
      real*8        F3_6um     (nr__2) ! (mJy) IRAC 3.6 micron flux density
      real*4        e_F3_6um   (nr__2) ! (mJy) Uncertainty in F3.6um
      character*2   l_F4_5um   (nr__2) ! [<= ] Limit flag on F4.5um
      real*8        F4_5um     (nr__2) ! (mJy) ? IRAC 4.5 micron flux density
      real*4        e_F4_5um   (nr__2) ! (mJy) ? Uncertainty in F4.5um
      character*2   l_F5_8um   (nr__2) ! [<= ] Limit flag on F5.8um
      real*4        F5_8um     (nr__2) ! (mJy) IRAC 5.8 micron flux density
      real*4        e_F5_8um   (nr__2) ! (mJy) ? Uncertainty in F5.8 um
      character*2   l_F8_0um   (nr__2) ! [<= ] Limit flag on F8.0um
      real*8        F8_0um     (nr__2) ! (mJy) ? IRAC 8.0 micron flux density
      real*4        e_F8_0um   (nr__2) ! (mJy) ? Uncertainty in F8.0um
      character*2   l_F24um    (nr__2) ! [<= ] Limit flag on F24um
      real*8        F24um      (nr__2) ! (mJy) ? IRAC 24 micron flux density
      real*4        e_F24um    (nr__2) ! (mJy) ? Uncertainty in F24 um

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

C  Declarations for 'table6.dat'	! Spitzer Flux Densities for Normal Galaxy Sample

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

      character*1   Type_1     (nr__3) ! Galaxy type (G1)
      character*12  Name_1     (nr__3) ! Source name
      real*4        F3_6um_1   (nr__3) ! (mJy) IRAC 3.6 micron flux density
      real*4        e_F3_6um_1 (nr__3) ! (mJy) Uncertainty in F3.6um
      real*4        F4_5um_1   (nr__3) ! (mJy) IRAC 4.5 micron flux density
      real*4        e_F4_5um_1 (nr__3) ! (mJy) Uncertainty in F4.5um
      character*2   l_F5_8um_1 (nr__3) ! [<= ] Limit flag on F5.8um
      real*4        F5_8um_1   (nr__3) ! (mJy) IRAC 5.8 micron flux density
      real*4        e_F5_8um_1 (nr__3) ! (mJy) ? Uncertainty in F5.8um
      character*2   l_F8_0um_1 (nr__3) ! [<= ] Limit flag on F8.0um
      real*8        F8_0um_1   (nr__3) ! (mJy) IRAC 8.0 micron flux density
      real*4        e_F8_0um_1 (nr__3) ! (mJy) ? Uncertainty in F8.0um
      character*2   l_F24um_1  (nr__3) ! [<= ] Limit flag on F24um
      real*8        F24um_1    (nr__3) ! (mJy) IRAC 24 micron flux density
      real*4        e_F24um_1  (nr__3) ! (mJy) ? Uncertainty in F24um

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

C  Loading file 'table1.dat'	! Interacting Galaxy Sample

C  Format for file interpretation

    1 format(A8,1X,A17,F5.1,1X,F5.1,1X,F5.2,1X,A51,1X,A27,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,35
        read(1,'(A126)')ar__
        read(ar__,1)
     +  Arp(i__),OName(i__),Dist(i__),Sep(i__),logL_FIR(i__),
     +  Notes(i__),SpType(i__),Flags(i__)
        if (idig(ar__(39:43)).EQ.0) logL_FIR(i__) =  rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Arp(i__),OName(i__),Dist(i__),Sep(i__),logL_FIR(i__),
     +  Notes(i__),SpType(i__),Flags(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Normal Galaxy Sample

C  Format for file interpretation

    2 format(A1,1X,A11,3X,A5,3X,F5.2,2X,A1,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,42
        read(1,'(A37)')ar__1
        read(ar__1,2)
     +  Type(i__),Name(i__),MType(i__),Dist_1(i__),l_logL_IR(i__),
     +  logL_IR(i__)
c    ..............Just test output...........
        write(6,2)
     +  Type(i__),Name(i__),MType(i__),Dist_1(i__),l_logL_IR(i__),
     +  logL_IR(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Spitzer Flux Densities for Interacting Galaxy Sample

C  Format for file interpretation

    3 format(
     +  A8,1X,A9,1X,A23,1X,F7.2,1X,F4.2,1X,A2,F7.2,1X,F4.2,1X,A2,F6.1,
     +  1X,F5.1,1X,A2,F8.2,1X,F6.2,1X,A2,F7.1,1X,F5.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,116
        read(1,'(A119)')ar__2
        read(ar__2,3)
     +  Arp_1(i__),Comp(i__),OName_1(i__),F3_6um(i__),e_F3_6um(i__),
     +  l_F4_5um(i__),F4_5um(i__),e_F4_5um(i__),l_F5_8um(i__),
     +  F5_8um(i__),e_F5_8um(i__),l_F8_0um(i__),F8_0um(i__),
     +  e_F8_0um(i__),l_F24um(i__),F24um(i__),e_F24um(i__)
        if(ar__2(59:65) .EQ. '') F4_5um(i__) = rNULL__
        if(ar__2(67:70) .EQ. '') e_F4_5um(i__) = rNULL__
        if(ar__2(81:85) .EQ. '') e_F5_8um(i__) = rNULL__
        if(ar__2(89:96) .EQ. '') F8_0um(i__) = rNULL__
        if(ar__2(98:103) .EQ. '') e_F8_0um(i__) = rNULL__
        if(ar__2(107:113) .EQ. '') F24um(i__) = rNULL__
        if(ar__2(115:119) .EQ. '') e_F24um(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Arp_1(i__),Comp(i__),OName_1(i__),F3_6um(i__),e_F3_6um(i__),
     +  l_F4_5um(i__),F4_5um(i__),e_F4_5um(i__),l_F5_8um(i__),
     +  F5_8um(i__),e_F5_8um(i__),l_F8_0um(i__),F8_0um(i__),
     +  e_F8_0um(i__),l_F24um(i__),F24um(i__),e_F24um(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Spitzer Flux Densities for Normal Galaxy Sample

C  Format for file interpretation

    4 format(
     +  A1,1X,A12,1X,F6.1,1X,F5.1,1X,F6.1,1X,F5.1,1X,A2,F6.1,1X,F5.1,
     +  1X,A2,F7.1,1X,F5.1,1X,A2,F7.1,1X,F5.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,40
        read(1,'(A87)')ar__3
        read(ar__3,4)
     +  Type_1(i__),Name_1(i__),F3_6um_1(i__),e_F3_6um_1(i__),
     +  F4_5um_1(i__),e_F4_5um_1(i__),l_F5_8um_1(i__),F5_8um_1(i__),
     +  e_F5_8um_1(i__),l_F8_0um_1(i__),F8_0um_1(i__),e_F8_0um_1(i__),
     +  l_F24um_1(i__),F24um_1(i__),e_F24um_1(i__)
        if(ar__3(51:55) .EQ. '') e_F5_8um_1(i__) = rNULL__
        if(ar__3(67:71) .EQ. '') e_F8_0um_1(i__) = rNULL__
        if(ar__3(83:87) .EQ. '') e_F24um_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Type_1(i__),Name_1(i__),F3_6um_1(i__),e_F3_6um_1(i__),
     +  F4_5um_1(i__),e_F4_5um_1(i__),l_F5_8um_1(i__),F5_8um_1(i__),
     +  e_F5_8um_1(i__),l_F8_0um_1(i__),F8_0um_1(i__),e_F8_0um_1(i__),
     +  l_F24um_1(i__),F24um_1(i__),e_F24um_1(i__)
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