CDS.headStuff2 FORTRAN Generation
()

Conversion of standardized ReadMe file for catalog 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-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. VII/173             Catalogue of Seyfert Galaxies    (Lipovetsky+, 1988)
*================================================================================
*A Catalogue of Seyfert Galaxies
*     Lipovetsky V.A., Neizvestny S.I., Neizvestnaya O.M.
*    <Comm. SAO (USSR) 55 (1988)>
*    =1988SoSAO..55....5L
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1'	! General list of objects

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

C  1950.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg      (nr__) ! (deg) Right Ascension 1950.0
      real*8        DEdeg      (nr__) ! (deg)     Declination 1950.0
C  ---------------------------------- ! (position vector(s) in degrees)

      character*8   IAU        (nr__) ! IAU name
      character*1   m_IAU      (nr__) ! Multiplicity index on IAU name
      integer*4     RAh        (nr__) ! (h) Right ascension 1950.0
      integer*4     RAm        (nr__) ! (min) Right ascension 1950.0
      real*4        RAs        (nr__) ! (s) Right ascension 1950.0
      character*1   DE_        (nr__) ! Declination sign
      integer*4     DEd        (nr__) ! (deg) Declination 1950.0
      integer*4     DEm        (nr__) ! (arcmin) Declination 1950.0
      integer*4     DEs        (nr__) ! (arcsec) Declination 1950.0
      real*8        z0         (nr__) ! []? Galactocentric redshift
      integer*4     r_z0       (nr__) ! []? Redshift reference
      real*4        Type       (nr__) ! Seyfert type
      character*1   n_Type     (nr__) ! [:?*] Note on Seyfert type (1)
      integer*4     r_Type     (nr__) ! Seyfert type reference
      real*4        Bmag       (nr__) ! (mag) []? B magnitude
      integer*4     r_Bmag     (nr__) ! []? B magnitude reference
      real*4        E_B_V      (nr__) ! (mag) []? Galactic reddening
      real*4        BmagAbs    (nr__) ! (mag) []? Absolute B magnitude
*Note (1):     :  = probable Sy type
*              ?  = uncertain entries of Sy galaxies
*              *  = nucleus of bright galaxy with Sy feature

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

C  Declarations for 'table2'	! Morphological and photoelectric data

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

      character*8   IAU_1      (nr__1) ! IAU name
      character*1   m_IAU_1    (nr__1) ! Multiplicity index on IAU name
      character*10  Mtype      (nr__1) ! Morphological type
      character*7   Memb       (nr__1) ! Membership (1)
      real*4        MajDiam    (nr__1) ! (arcmin) []? Isophotal major diameter
      character*1   u_MajDiam  (nr__1) ! Uncertainty flag on MajDiam
      real*4        b_a        (nr__1) ! []? Axis ratio
      character*1   u_b_a      (nr__1) ! Uncertainty flag on a/b
      integer*4     Ref        (nr__1) ! []? Data reference or finding chart reference
      real*4        Aper       (nr__1) ! (arcmin) []? Aperture diameter of photoelectric
*                                       measurements
      real*4        Vmag       (nr__1) ! (mag) []? V magnitude
      real*4        U_B        (nr__1) ! (mag) []? U-B colour index
      real*4        B_V        (nr__1) ! (mag) []? B-V colour index
      real*4        V_R        (nr__1) ! (mag) []? V-R colour index
      character*1   u_V_R      (nr__1) ! Uncertainty flag on V-R
      integer*4     r_Vmag     (nr__1) ! []? Reference for UBV photometry
*Note (1):  CL    - a member of a cluster of galaxies;
*           GR    - a member of a group of galaxies;
*           PA    - a member of a pair of galaxies;
*          GR/CL  - it is reported in literature about SyG belonging to
*                    a cluster, which is studied not sufficiently well and there
*                    is no concrete information about its richness.
*           Possible membership is noted ":"

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

C  Declarations for 'table3'	! Different spectral range fluxes
                                  (X-ray, optics, IR, radio)

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

      character*8   IAU_2      (nr__2) ! IAU name
      character*1   m_IAU_2    (nr__2) ! Multiplicity index on IAU name
      character*1   l_Xhard    (nr__2) ! [<] Limit flag on Xhard
      real*4        Xhard      (nr__2) ! (uJy) []? Hard X-ray monochromatic flux
*                                       (2-10keV or 2-6keV region)
      character*1   n_Xhard    (nr__2) ! [:V] ':' for uncertainty, 'V' for variable
      integer*4     r_Xhard    (nr__2) ! []? Reference for Xhard
      character*1   l_Xsoft    (nr__2) ! [<] Limit flag on Xsoft
      real*4        Xsoft      (nr__2) ! (uJy) []? Soft X-ray monochromatic flux (0.4-4.5keV
*                                       or 0.3-3.5keV) in 10-29erg/cm2/s/Hz
      integer*4     r_Xsoft    (nr__2) ! []? Reference for Xsoft
      character*1   l_Hbeta    (nr__2) ! [<] Limit flag on Hbeta
      real*8        Hbeta      (nr__2) ! (10-16W/m2) []? Hbeta-line flux in 10-13erg/cm2/s
      character*1   l_O_H      (nr__2) ! [<>] Limit flag on O/Hbeta
      real*4        O_H        (nr__2) ! []? Nebular(5007A[OIII])/Hbeta ratio
      character*1   n_O_H      (nr__2) ! [:V] ':' for uncertainty, 'V' for variable
      integer*4     r_HFlux    (nr__2) ! []? References for H-line flux and O/H ratio
      real*4        JFlux      (nr__2) ! (mJy) []? Flux in Johnson J band (1)
      real*8        HFlux      (nr__2) ! (mJy) []? Flux in Johnson H band (1)
      real*8        KFlux      (nr__2) ! (mJy) []? Flux in Johnson K band (1)
      integer*4     r_JFlux    (nr__2) ! []? References for J,H,K- fluxes
      character*1   l_LFlux    (nr__2) ! [<] Limit flag on LF
      real*8        LFlux      (nr__2) ! (mJy) []? Johnson L-band or 3.5{mu}m-band flux
      integer*4     r_LFlux    (nr__2) ! []? L-flux reference
      character*1   l_NFlux    (nr__2) ! [<] Limit flag on NF
      real*8        NFlux      (nr__2) ! (mJy) []? Johnson N-band (10{mu}m-band) flux
      integer*4     r_NFlux    (nr__2) ! []? N-flux reference
      character*1   l_F6cm     (nr__2) ! [<] Limit flag on F6cm
      real*8        F6cm       (nr__2) ! (Jy) []? 6-cm radio flux (or 7.6 or 8.2cm)
      character*1   n_F6cm     (nr__2) ! [*#] * indicates that F6cm contains 7.6cm flux
*                                       # indicates that F6cm contains 8.2cm flux
      integer*4     r_F6cm     (nr__2) ! []? 6-cm flux references
      character*1   l_F11cm    (nr__2) ! [<] Limit flag on F11cm
      real*8        F11cm      (nr__2) ! (Jy) []? 11-cm radio flux (or 21, 31 or 73cm)
      character*1   n_F11cm    (nr__2) ! [#+*] # means that F11cm contains 21cm flux
*                                        + means that F11cm contains 31cm flux
*                                        * means that F11cm contains 73cm flux
      integer*4     r_F11cm    (nr__2) ! []? 11-cm flux references
*Note (1):  When the data published included measurements of the object
*            obtained with several apertures, we give those obtained with the
*            maximum one.

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

C  Declarations for 'table3a'	! Far-infrared fluxes

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

      character*8   IAU_3      (nr__3) ! IAU name
      character*1   m_IAU_3    (nr__3) ! Multiplicity index on IAU name
      character*1   l_F12      (nr__3) ! Limit flag on F12
      real*4        F12        (nr__3) ! (Jy) []? Far IR flux at 12 {mu}m
      character*1   l_F25      (nr__3) ! Limit flag on F25
      real*4        F25        (nr__3) ! (Jy) Far IR flux at 25 {mu}m
      character*1   l_F60      (nr__3) ! Limit flag on F60
      real*4        F60        (nr__3) ! (Jy) Far IR flux at 60 {mu}m
      character*1   l_F100     (nr__3) ! Limit flag on F100
      real*8        F100       (nr__3) ! (Jy) Far IR flux at 100 {mu}m

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

C  Declarations for 'table4'	! List of common names

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

      character*8   IAU_4      (nr__4) ! IAU name
      character*1   m_IAU_4    (nr__4) ! Multiplicity index on IAU name
      character*101 Name       (nr__4) ! Other names

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

C  Declarations for 'table5'	! Finding list of common names

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

      character*18  Name_1     (nr__5) ! Name
      character*8   IAU_5      (nr__5) ! IAU name
      character*1   u_IAU      (nr__5) ! [?] uncertainty flag on IAU name

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

C  Loading file 'table1'	! General list of objects

C  Format for file interpretation

    1 format(
     +  A8,A1,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,F7.4,1X,I3,1X,
     +  F3.1,A1,I3,1X,F5.2,1X,I3,1X,F4.2,1X,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1')
      write(6,*) '....Loading file: table1'
      do i__=1,959
        read(1,'(A70)')ar__
        read(ar__,1)
     +  IAU(i__),m_IAU(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),z0(i__),r_z0(i__),Type(i__),
     +  n_Type(i__),r_Type(i__),Bmag(i__),r_Bmag(i__),E_B_V(i__),
     +  BmagAbs(i__)
        if(ar__(30:36) .EQ. '') z0(i__) = rNULL__
        if(ar__(38:40) .EQ. '') r_z0(i__) = iNULL__
        if(ar__(50:54) .EQ. '') Bmag(i__) = rNULL__
        if(ar__(56:58) .EQ. '') r_Bmag(i__) = iNULL__
        if(ar__(60:63) .EQ. '') E_B_V(i__) = rNULL__
        if(ar__(65:70) .EQ. '') BmagAbs(i__) = rNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  IAU(i__),m_IAU(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),z0(i__),r_z0(i__),Type(i__),
     +  n_Type(i__),r_Type(i__),Bmag(i__),r_Bmag(i__),E_B_V(i__),
     +  BmagAbs(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2'	! Morphological and photoelectric data

C  Format for file interpretation

    2 format(
     +  A8,A1,1X,A10,2X,A7,1X,F5.2,A1,F4.2,A1,1X,I3,1X,F5.2,1X,F5.2,
     +  1X,F5.2,F5.2,1X,F4.2,A1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2')
      write(6,*) '....Loading file: table2'
      do i__=1,704
        read(1,'(A78)')ar__1
        read(ar__1,2)
     +  IAU_1(i__),m_IAU_1(i__),Mtype(i__),Memb(i__),MajDiam(i__),
     +  u_MajDiam(i__),b_a(i__),u_b_a(i__),Ref(i__),Aper(i__),
     +  Vmag(i__),U_B(i__),B_V(i__),V_R(i__),u_V_R(i__),r_Vmag(i__)
        if(ar__1(31:35) .EQ. '') MajDiam(i__) = rNULL__
        if(ar__1(37:40) .EQ. '') b_a(i__) = rNULL__
        if(ar__1(43:45) .EQ. '') Ref(i__) = iNULL__
        if(ar__1(47:51) .EQ. '') Aper(i__) = rNULL__
        if(ar__1(53:57) .EQ. '') Vmag(i__) = rNULL__
        if(ar__1(59:63) .EQ. '') U_B(i__) = rNULL__
        if(ar__1(64:68) .EQ. '') B_V(i__) = rNULL__
        if(ar__1(70:73) .EQ. '') V_R(i__) = rNULL__
        if(ar__1(76:78) .EQ. '') r_Vmag(i__) = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  IAU_1(i__),m_IAU_1(i__),Mtype(i__),Memb(i__),MajDiam(i__),
     +  u_MajDiam(i__),b_a(i__),u_b_a(i__),Ref(i__),Aper(i__),
     +  Vmag(i__),U_B(i__),B_V(i__),V_R(i__),u_V_R(i__),r_Vmag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3'	! Different spectral range fluxes
*                                  (X-ray, optics, IR, radio)

C  Format for file interpretation

    3 format(
     +  A8,A1,1X,A1,F6.3,A1,I3,A1,F6.3,1X,I3,A1,F7.3,A1,F5.2,A1,I3,1X,
     +  F6.2,F7.2,F7.2,1X,I3,A1,F7.2,1X,I3,1X,A1,F7.1,1X,I3,1X,A1,
     +  F7.3,A1,I3,A1,F7.3,A1,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3')
      write(6,*) '....Loading file: table3'
      do i__=1,542
        read(1,'(A125)')ar__2
        read(ar__2,3)
     +  IAU_2(i__),m_IAU_2(i__),l_Xhard(i__),Xhard(i__),n_Xhard(i__),
     +  r_Xhard(i__),l_Xsoft(i__),Xsoft(i__),r_Xsoft(i__),
     +  l_Hbeta(i__),Hbeta(i__),l_O_H(i__),O_H(i__),n_O_H(i__),
     +  r_HFlux(i__),JFlux(i__),HFlux(i__),KFlux(i__),r_JFlux(i__),
     +  l_LFlux(i__),LFlux(i__),r_LFlux(i__),l_NFlux(i__),NFlux(i__),
     +  r_NFlux(i__),l_F6cm(i__),F6cm(i__),n_F6cm(i__),r_F6cm(i__),
     +  l_F11cm(i__),F11cm(i__),n_F11cm(i__),r_F11cm(i__)
        if(ar__2(12:17) .EQ. '') Xhard(i__) = rNULL__
        if(ar__2(19:21) .EQ. '') r_Xhard(i__) = iNULL__
        if(ar__2(23:28) .EQ. '') Xsoft(i__) = rNULL__
        if(ar__2(30:32) .EQ. '') r_Xsoft(i__) = iNULL__
        if(ar__2(34:40) .EQ. '') Hbeta(i__) = rNULL__
        if(ar__2(42:46) .EQ. '') O_H(i__) = rNULL__
        if(ar__2(48:50) .EQ. '') r_HFlux(i__) = iNULL__
        if(ar__2(52:57) .EQ. '') JFlux(i__) = rNULL__
        if(ar__2(58:64) .EQ. '') HFlux(i__) = rNULL__
        if(ar__2(65:71) .EQ. '') KFlux(i__) = rNULL__
        if(ar__2(73:75) .EQ. '') r_JFlux(i__) = iNULL__
        if(ar__2(77:83) .EQ. '') LFlux(i__) = rNULL__
        if(ar__2(85:87) .EQ. '') r_LFlux(i__) = iNULL__
        if(ar__2(90:96) .EQ. '') NFlux(i__) = rNULL__
        if(ar__2(98:100) .EQ. '') r_NFlux(i__) = iNULL__
        if(ar__2(103:109) .EQ. '') F6cm(i__) = rNULL__
        if(ar__2(111:113) .EQ. '') r_F6cm(i__) = iNULL__
        if(ar__2(115:121) .EQ. '') F11cm(i__) = rNULL__
        if(ar__2(123:125) .EQ. '') r_F11cm(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  IAU_2(i__),m_IAU_2(i__),l_Xhard(i__),Xhard(i__),n_Xhard(i__),
     +  r_Xhard(i__),l_Xsoft(i__),Xsoft(i__),r_Xsoft(i__),
     +  l_Hbeta(i__),Hbeta(i__),l_O_H(i__),O_H(i__),n_O_H(i__),
     +  r_HFlux(i__),JFlux(i__),HFlux(i__),KFlux(i__),r_JFlux(i__),
     +  l_LFlux(i__),LFlux(i__),r_LFlux(i__),l_NFlux(i__),NFlux(i__),
     +  r_NFlux(i__),l_F6cm(i__),F6cm(i__),n_F6cm(i__),r_F6cm(i__),
     +  l_F11cm(i__),F11cm(i__),n_F11cm(i__),r_F11cm(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3a'	! Far-infrared fluxes

C  Format for file interpretation

    4 format(A8,A1,2X,A1,F5.2,3X,A1,F6.2,A1,F6.2,3X,A1,F7.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3a')
      write(6,*) '....Loading file: table3a'
      do i__=1,433
        read(1,'(A45)')ar__3
        read(ar__3,4)
     +  IAU_3(i__),m_IAU_3(i__),l_F12(i__),F12(i__),l_F25(i__),
     +  F25(i__),l_F60(i__),F60(i__),l_F100(i__),F100(i__)
        if(ar__3(13:17) .EQ. '') F12(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  IAU_3(i__),m_IAU_3(i__),l_F12(i__),F12(i__),l_F25(i__),
     +  F25(i__),l_F60(i__),F60(i__),l_F100(i__),F100(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4'	! List of common names

C  Format for file interpretation

    5 format(1X,A8,A1,1X,A101)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4')
      write(6,*) '....Loading file: table4'
      do i__=1,955
        read(1,'(A112)')ar__4
        read(ar__4,5)IAU_4(i__),m_IAU_4(i__),Name(i__)
c    ..............Just test output...........
        write(6,5)IAU_4(i__),m_IAU_4(i__),Name(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5'	! Finding list of common names

C  Format for file interpretation

    6 format(A18,12X,A8,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5')
      write(6,*) '....Loading file: table5'
      do i__=1,2549
        read(1,'(A39)')ar__5
        read(ar__5,6)Name_1(i__),IAU_5(i__),u_IAU(i__)
c    ..............Just test output...........
        write(6,6)Name_1(i__),IAU_5(i__),u_IAU(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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