FORTRAN Generation
(/./ftp/cats/J/ApJ/705/14)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/705/14 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-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/ApJ/705/14       Seyfert galaxies in the mid-IR          (Deo+, 2009)
*================================================================================
*The mid-infrared continua of Seyfert galaxies.
*    Deo R.P., Richards G.T., Crenshaw D.M., Kraemer S.B.
*   <Astrophys. J., 705, 14-31 (2009)>
*   =2009ApJ...705...14D
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Observation summary, redshifts, optical Seyfert
                             types, and extraction apertures

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

C  J2000 position composed of: RAdeg DEdeg
      character*19  Name        ! Galaxy name
      real*8        z           ! NED redshift
      character*4   SyT         ! Seyfert type (1)
      real*4        b_a         ! Axial ratio (2)
      integer*4     AOR         ! Spitzer archive AOR number
      character*1   Obs         ! [SM] Observing mode for Spitzer/IRS spectrum:
*                                 (S=staring or M=mapping)
      real*8        RAdeg       ! (deg) ? Extraction rectangle for mapping mode:
*                                 right ascension in decimal degrees (J2000) (3)
      real*8        DEdeg       ! (deg) ? Extraction rectangle for mapping mode:
*                                 declination in decimal degrees (J2000) (3)
      real*4        Ap1         ! (arcsec) ? First extraction aperture
      real*4        Ap2         ! (arcsec) ? Second extraction aperture
      real*8        scale       ! (pc/arcsec) Radial extent of the extraction aperture for
*                                 the galaxy
*Note (1): Seyfert type from Veron-Cetty & Veron (2006, Cat. VII/248; superseded
*          by Cat. VII/258) as follows:
*  S1  = Seyfert 1 optical spectrum;
*  S1h = broad polarized Balmer lines detected;
*  S1i = broad Paschen lines observed in the infrared;
*  S1n = narrow-line Seyfert 1;
*  S1.0, S1.2, S1.5, S1.8, and S1.9 = intermediate Seyfert galaxies (Note: In
*        this paper, we consider S1.0-1.5 to be broadly Seyfert 1s in all
*        analysis.);
*  S2  = Seyfert 2 spectrum;
*  S3  = LINER;
*  S3b = LINER with broad Balmer lines;
*  S3h = LINER with broad polarized Balmer lines detected;
*  H2  = nuclear HII region
*Note (2): Axial ratio, minor-to-major axis ratio of host galaxy obtained from
*     NED, for host galaxies where b/a was not available we assumed it to be 1.
*Note (3): Extraction rectangle for mapping-mode spectra: R.A. and Decl. of
*     center-point.

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

C  Declarations for 'table2.dat'	! Continuum flux densities and apparent optical
                             depth at 9.7um

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

      character*19  Name_1      ! Galaxy name
      real*4        v5_5um      ! (Jy) ? Spitzer/IRS 5.5um continuum flux density
      real*4        v10um       ! (Jy) ? Spitzer/IRS 10um continuum flux density
      real*4        v14_7um     ! (Jy) ? Spitzer/IRS 14.7um continuum flux density
      real*4        v20um       ! (Jy) ? Spitzer/IRS 20um continuum flux density
      real*4        v30um       ! (Jy) ? Spitzer/IRS 30um continuum flux density
      real*8        tau         ! ? Observed optical depth (apparent) at 9.7um
      real*4        e_5_5um     ! (Jy) ? 1{sigma} error in 5.5um
      real*4        e_10um      ! (Jy) ? 1{sigma} error in 10um
      real*4        e_14_7um    ! (Jy) ? 1{sigma} error in 14.7um
      real*4        e_20um      ! (Jy) ? 1{sigma} error in 20um
      real*4        e_30um      ! (Jy) ? 1{sigma} error in 30um

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

C  Declarations for 'table3.dat'	! Starburst-to-AGN flux density ratios

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

      character*19  Name_2      ! Galaxy name
      real*4        R5_5        ! ? Starbust-to-AGN flux density ratio at 5.5um (1)
      real*8        R10         ! ? Starbust-to-AGN flux density ratio at 10um (1)
      real*8        R14_7       ! ? Starbust-to-AGN flux density ratio at 14.7um (1)
      real*8        R20         ! ? Starbust-to-AGN flux density ratio at 20um (1)
      real*8        R30         ! ? Starbust-to-AGN flux density ratio at 30um (1)
*Note (1): A missing ratio implies either the spectrum was incomplete at
*     that wavelength or the AGN contribution was negative after starburst
*     subtraction. A very large ratio indicates very weak AGN contribution,
*     these typically occur at or beyond 14.7um in a few galaxies. See
*     Figure 10 for a visual representation of these ratios, and Table 1 for
*     aperture size in parsecs. Even in nearby AGN, Spitzer spectra sample
*     regions of size ~1kpc. See Section 5 for discussion of the starburst
*     subtraction process and results.

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

C  Loading file 'table1.dat'	! Observation summary, redshifts, optical Seyfert
*                             types, and extraction apertures

C  Format for file interpretation

    1 format(
     +  A19,1X,F8.6,1X,A4,1X,F4.2,1X,I8,1X,A1,1X,F10.6,1X,F10.6,1X,
     +  F6.3,1X,F6.3,1X,F8.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,109
        read(1,'(A94)')ar__
        read(ar__,1)Name,z,SyT,b_a,AOR,Obs,RAdeg,DEdeg,Ap1,Ap2,scale
        if(ar__(51:60) .EQ. '') RAdeg = rNULL__
        if(ar__(62:71) .EQ. '') DEdeg = rNULL__
        if(ar__(73:78) .EQ. '') Ap1 = rNULL__
        if(ar__(80:85) .EQ. '') Ap2 = rNULL__
c    ..............Just test output...........
        write(6,1)Name,z,SyT,b_a,AOR,Obs,RAdeg,DEdeg,Ap1,Ap2,scale
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Continuum flux densities and apparent optical
*                             depth at 9.7um

C  Format for file interpretation

    2 format(
     +  A19,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,F7.4,1X,E8.3,
     +  1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,108
        read(1,'(A117)')ar__1
        read(ar__1,2)
     +  Name_1,v5_5um,v10um,v14_7um,v20um,v30um,tau,e_5_5um,e_10um,
     +  e_14_7um,e_20um,e_30um
        if(ar__1(21:28) .EQ. '') v5_5um = rNULL__
        if(ar__1(30:37) .EQ. '') v10um = rNULL__
        if(ar__1(39:46) .EQ. '') v14_7um = rNULL__
        if(ar__1(48:55) .EQ. '') v20um = rNULL__
        if(ar__1(57:64) .EQ. '') v30um = rNULL__
        if(ar__1(66:72) .EQ. '') tau = rNULL__
        if(ar__1(74:81) .EQ. '') e_5_5um = rNULL__
        if(ar__1(83:90) .EQ. '') e_10um = rNULL__
        if(ar__1(92:99) .EQ. '') e_14_7um = rNULL__
        if(ar__1(101:108) .EQ. '') e_20um = rNULL__
        if(ar__1(110:117) .EQ. '') e_30um = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Name_1,v5_5um,v10um,v14_7um,v20um,v30um,tau,e_5_5um,e_10um,
     +  e_14_7um,e_20um,e_30um
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Starburst-to-AGN flux density ratios

C  Format for file interpretation

    3 format(A19,1X,F6.4,1X,F7.4,1X,F9.4,1X,F8.4,1X,F7.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,93
        read(1,'(A61)')ar__2
        read(ar__2,3)Name_2,R5_5,R10,R14_7,R20,R30
        if(ar__2(21:26) .EQ. '') R5_5 = rNULL__
        if(ar__2(28:34) .EQ. '') R10 = rNULL__
        if(ar__2(36:44) .EQ. '') R14_7 = rNULL__
        if(ar__2(46:53) .EQ. '') R20 = rNULL__
        if(ar__2(55:61) .EQ. '') R30 = rNULL__
c    ..............Just test output...........
        write(6,3)Name_2,R5_5,R10,R14_7,R20,R30
c    .......End.of.Just test output...........
      end do
      close(1)

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