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