Conversion of standardized ReadMe file for
file /./ftp/cats/VII/213 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-Apr-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__ parameter (rNULL__=--2147483648.) ! NULL real number parameter (iNULL__=(-2147483647-1)) ! NULL int number integer idig ! testing NULL number C============================================================================= Cat. VII/213 Hickson's Compact groups of Galaxies (Hickson+ 1982-1994) *================================================================================ *Systematic Properties of Compact Groups of Galaxies * Hickson P. * <Astrophys. J. 255, 382 + Erratum 259, 930 (1982)> * =1982ApJ...255..382H *A photometric catalog of compact groups of galaxies. * Hickson P., Kindl E., Auman J.R. * <Astrophys. J. Suppl. Ser. 70, 687 (1989)> * =1989ApJS...70..687H *The luminosity function of compact groups of galaxies * Mendes de Oliveira C., Hickson P. * <Astrophys. J. 380, 30 (1991)> * =1991ApJ...380...30M *Dynamical properties of compact groups of galaxies. * Hickson P., Mendes de Oliveira C., Huchra J.P., Palumbo G.G.C. * <Astrophys. J. 399, 353 (1992)> * =1992ApJ...399..353H *Morphology of galaxies in compact groups. * Mendes de Oliveira C., Hickson P. * <Astrophys. J., 427, 684 (1994)> * =1994ApJ...427..684M C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'groups.dat' ! Catalog of groups (1982 paper, and errata) integer*4 nr__ parameter (nr__=100) ! Number of records character*72 ar__ ! Full-size record C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__) ! (deg) Right Ascension 1950 real*8 DEdeg (nr__) ! (deg) Declination 1950 C ---------------------------------- ! (position vector(s) in degrees) integer*4 HCG (nr__) ! [1/100]+ Number of the group in this catalog integer*4 RAh (nr__) ! (h) Right Ascension 1950 hours integer*4 RAm (nr__) ! (min) Right Ascension minutes integer*4 RAs (nr__) ! (s) Right Ascension seconds character*1 DE_ (nr__) ! Declination sign integer*4 DEd (nr__) ! (deg) Declination degrees (1950) integer*4 DEm (nr__) ! (arcmin) Declination arcminutes integer*4 DEs (nr__) ! (arcsec) [0,60] Declination arcseconds character*1 Dyna (nr__) ! [*] '*' when group has dynamical parameters character*2 Type (nr__) ! *Group type integer*4 MCount (nr__) ! *Member count real*4 AngSize (nr__) ! (arcmin) *Angular size real*4 TotMag (nr__) ! (mag) *Total magnitude of group real*4 Emag (nr__) ! (mag) *Brightest member magnitude real*4 z (nr__) ! *? Redshift character*6 Aname1 (nr__) ! *Alternate name 1 character*6 Aname2 (nr__) ! *Alternate name 2 character*6 Aname3 (nr__) ! *Alternate name 3 character*6 Aname4 (nr__) ! *Alternate name 4 *Note on Type: * Two-character code (a letter followed by a number) * classifying the group according to its two brightest members: * S The brightest galaxy is a spiral galaxy. * E The brightest galaxy is not a spiral galaxy. * 1 m(b) - m(a) >= 1.0 * 2 0.5 <= m(b) - m(a) < 1.0 * 3 m(b) - m(a) < 0.5 * where m(a) and m(b) are the estimated red magnitudes of the brightest * and second brightest members, respectively, of the cluster. *Note on MCount: * Number of galaxies in the group. Only those galaxies within three * magnitudes of the brightest galaxy are counted as members. *Note on AngSize: * Angular diameter, in arcminutes, of the smallest circle containing the * geometric centers of all group members. *Note on TotMag: * Total estimated red magnitude of those galaxies counted as group members * (see member count, above.) *Note on Emag: * Estimated red magnitude of the brightest galaxy in the group, in units * of magnitude. *Note on z: * Corrected redshift of the brightest galaxy. In all but ten cases this * field is blank; prefer the group redshift supplied in the file * "dynamics.dat". *Note on Aname1, Aname2, Aname3, Aname4: * Other designations of the group, if any (at most, four). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dynamics.dat' ! Dynamical properties (1992 paper, table3) integer*4 nr__1 parameter (nr__1=92) ! Number of records character*94 ar__1 ! Full-size record integer*4 HCG_1 (nr__1) ! [1/100]+= Group number real*4 z_1 (nr__1) ! Group Redshift derived from the median * galaxy heliocentric velocity integer*4 n (nr__1) ! Number of galaxies with accordant velocities real*4 log_sV (nr__1) ! ([km/s]) Radial velocity dispersion real*4 log_V (nr__1) ! ([km/s]) ? Estimated instrinsic 3D velocity dispersion real*4 log_R (nr__1) ! ([kpc]) Median projected separation real*4 log_Ho_tc (nr__1) ! Crossing time (expressed in 1/Ho units) real*4 log_Mv (nr__1) ! ([g]) ? Virial mass of group in log. scale real*4 log_Mp (nr__1) ! ([g]) ? Projected mass real*4 log_Ma (nr__1) ! ([g]) ? Average mass real*4 log_Mm (nr__1) ! ([g]) ? Median mass real*4 log_M (nr__1) ! ([g]) ? Adopted mass real*4 log_L (nr__1) ! ([W]) Total blue luminosity * of the accordant galaxies real*4 log_M_L (nr__1) ! ([Sun]) ? Mass-to-light ratio real*4 log_rho (nr__1) ! ([g/cm3]) ? Mass density c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'galaxies.dat' ! Photometric and Velocity data for galaxies (combined tables from 1989 and 1992 papers) integer*4 nr__2 parameter (nr__2=463) ! Number of records character*130 ar__2 ! Full-size record C B1950 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_1 (nr__2) ! (deg) Right Ascension B1950 real*8 DEdeg_1 (nr__2) ! (deg) Declination B1950 C ---------------------------------- ! (position vector(s) in degrees) integer*4 HCG_2 (nr__2) ! [1/100]+= Hickson Group number (groups.dat) character*1 m_HCG (nr__2) ! [a-i]! Galaxy in HCG group integer*4 RAh_1 (nr__2) ! (h) Right Ascension B1950 (hours) integer*4 RAm_1 (nr__2) ! (min) Right Ascension B1950 (minutes) real*4 RAs_1 (nr__2) ! (s) Right Ascension B1950 (seconds) character*1 DE__1 (nr__2) ! Declination B1950 (sign) integer*4 DEd_1 (nr__2) ! (deg) Declination B1950 (degrees) integer*4 DEm_1 (nr__2) ! (arcmin) Declination B1950 (minutes) real*4 DEs_1 (nr__2) ! (arcsec) Declination B1950 (seconds) real*4 a (nr__2) ! (arcsec) Semi-major axis at 25.0mag/arcsec^2^ isophote real*4 b (nr__2) ! (arcsec) Semi-minor axis at 25.0mag/arcsec^2^ isophote character*4 MType (nr__2) ! Hubble Morphological type integer*4 T (nr__2) ! *Numeric morphological type character*1 morph (nr__2) ! [m] when more morphology in file "morpho.dat" integer*4 RVhel (nr__2) ! (km/s) ? Heliocentric velocity integer*4 e_RVhel (nr__2) ! (km/s) ? Estimated rms in RVhel integer*4 q_RVhel (nr__2) ! *[0/2]? Confidence code on RVhel real*4 Bmag (nr__2) ! (mag) B magnitude within 24.5mag/arcsec^2^ isophote integer*4 q_Bmag (nr__2) ! *[0,4] Confidence level for Bmag real*4 DBmag (nr__2) ! (arcsec) *Diameter of B=24.5mag/arcsec^2^ isophote real*4 Rmag (nr__2) ! (mag) R magnitude within 24.0mag/arcsec^2^ isophote integer*4 q_Rmag (nr__2) ! *[0,4] Confidence level for Rmag real*4 DRmag (nr__2) ! *Diameter of R=24.0mag/arcsec^2^ isophote real*4 B_R (nr__2) ! (mag) Colour within the 24.5mag/arcsec^2^ isophote real*4 BTmag (nr__2) ! (mag) B_T_ asymptotic magnitude real*4 BTc (nr__2) ! (mag) B_T_ asymptotic magnitude corrected for * internal and external extinction real*4 e_BTmag (nr__2) ! (mag) Mean error on BTmag real*4 MB_5log_h (nr__2) ! (mag) ? Absolute magnitude of individual galaxies, * h=Ho/100, Ho=Hubble constant character*18 Names (nr__2) ! Other names *Note on T: the morphological type codes are: * --------------------------- * Hubble Code RC2 * --------------------------- * E0-E7 0 -4, -5 * S0, SB0 1 -2 * S0a, SB0a 2 0 * Sa, SBa 3 1 * Sab, SBab 4 2 * Sb, SBb 5 3 * Sbc, SBbc 6 4 * Sc, SBc 7 5 * Scd, SBcd 8 6 * Sd, SBd 9 7 * Sdm, SBdm 10 8 * Sm, SDm 11 9 * Im, IBm 12 10 * cI 13 11 * --------------------------- *Note on q_RVhel: the confidence code is * 0 = highest confidence: spectral features clearly visible * 1 = lower confidence: spectral features not clearly visible * 2 = lowest confidence: very noisy spectra with no obvious features. * Only 3 galaxies have q_RVhel=2. *Note on q_Bmag, q_Rmag: confidence levels are: * 0 = good * 1 = uncertain due to overlapping isophotes * 2 = uncertain due to contamination * 3 = uncertain due to poor calibration * 4 = uncertain due to clouds. * (value 5 for HCG 64c unexplained) *Note on DBmag, DRmag: * the diameter is defined as sqrt(Area/{pi}) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'morpho.dat' ! Morphology (1994 paper, table2) integer*4 nr__3 parameter (nr__3=210) ! Number of records character*92 ar__3 ! Full-size record integer*4 HCG_3 (nr__3) ! [1/100]+= Hickson Group number (groups.dat) character*1 m_HCG_1 (nr__3) ! [a-i]! Galaxy in HCG group integer*4 RVhel_1 (nr__3) ! (km/s) ? Heliocentric velocity character*4 MType_1 (nr__3) ! Morphological type character*69 Comments (nr__3) ! *Comments *Note on Comments: Abbreviations used in comments: * [string 'ansac' in 58a is unexplained] * asym. = asymmetric * emission = galaxy with emission lines in its nuclear spectrum * IR = infrared emitter (Hickson et al. 1989ApJ...341..679H, * Sulentic & de Mello Rabaca 1993ApJ...410..520S) * morpho = morphology * n.c.i. = nonconcentric isophotes * r.c. = rotation curve * radio = radio-loud galaxies (Menon & Hickson 1985; Menon 1993) * tail-like= galaxy has tail-like features * vel. = velocity * WR = galaxy exhibits Wolf-Rayet spectral features C============================================================================= C Loading file 'groups.dat' ! Catalog of groups (1982 paper, and errata) C Format for file interpretation 1 format( + I3,1X,I2,I2,I2,A1,I2,I2,I2,A1,1X,A2,I2,F5.1,F5.1,F5.1,F6.3,1X, + A6,1X,A6,1X,A6,1X,A6) C Effective file loading open(unit=1,status='old',file= +'groups.dat') write(6,*) '....Loading file: groups.dat' do i__=1,100 read(1,'(A72)')ar__ read(ar__,1) + HCG(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),Dyna(i__),Type(i__),MCount(i__), + AngSize(i__),TotMag(i__),Emag(i__),z(i__),Aname1(i__), + Aname2(i__),Aname3(i__),Aname4(i__) if(ar__(39:44) .EQ. '') z(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) + HCG(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),Dyna(i__),Type(i__),MCount(i__), + AngSize(i__),TotMag(i__),Emag(i__),z(i__),Aname1(i__), + Aname2(i__),Aname3(i__),Aname4(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 'dynamics.dat' ! Dynamical properties (1992 paper, table3) C Format for file interpretation 2 format( + I3,2X,F6.4,2X,I1,2X,F4.2,1X,F5.2,1X,F5.2,2X,F5.2,2X,F5.2,2X, + F5.2,2X,F5.2,2X,F5.2,2X,F5.2,2X,F5.2,1X,F5.2,1X,F6.2) C Effective file loading open(unit=1,status='old',file= +'dynamics.dat') write(6,*) '....Loading file: dynamics.dat' do i__=1,92 read(1,'(A94)')ar__1 read(ar__1,2) + HCG_1(i__),z_1(i__),n(i__),log_sV(i__),log_V(i__),log_R(i__), + log_Ho_tc(i__),log_Mv(i__),log_Mp(i__),log_Ma(i__), + log_Mm(i__),log_M(i__),log_L(i__),log_M_L(i__),log_rho(i__) if(ar__1(22:26) .EQ. '') log_V(i__) = rNULL__ if(ar__1(42:46) .EQ. '') log_Mv(i__) = rNULL__ if(ar__1(49:53) .EQ. '') log_Mp(i__) = rNULL__ if(ar__1(56:60) .EQ. '') log_Ma(i__) = rNULL__ if(ar__1(63:67) .EQ. '') log_Mm(i__) = rNULL__ if(ar__1(70:74) .EQ. '') log_M(i__) = rNULL__ if(ar__1(83:87) .EQ. '') log_M_L(i__) = rNULL__ if(ar__1(89:94) .EQ. '') log_rho(i__) = rNULL__ c ..............Just test output........... write(6,2) + HCG_1(i__),z_1(i__),n(i__),log_sV(i__),log_V(i__),log_R(i__), + log_Ho_tc(i__),log_Mv(i__),log_Mp(i__),log_Ma(i__), + log_Mm(i__),log_M(i__),log_L(i__),log_M_L(i__),log_rho(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'galaxies.dat' ! Photometric and Velocity data for galaxies * (combined tables from 1989 and 1992 papers) C Format for file interpretation 3 format( + I3,A1,1X,I2,I2,F5.2,A1,I2,I2,F4.1,1X,F5.1,1X,F4.1,2X,A4,1X,I2, + A1,I5,1X,I3,1X,I1,1X,F5.2,1X,I1,1X,F5.1,1X,F5.2,1X,I1,1X,F5.1, + 1X,F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.1,1X,A18) C Effective file loading open(unit=1,status='old',file= +'galaxies.dat') write(6,*) '....Loading file: galaxies.dat' do i__=1,463 read(1,'(A130)')ar__2 read(ar__2,3) + HCG_2(i__),m_HCG(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__), + DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),a(i__),b(i__), + MType(i__),T(i__),morph(i__),RVhel(i__),e_RVhel(i__), + q_RVhel(i__),Bmag(i__),q_Bmag(i__),DBmag(i__),Rmag(i__), + q_Rmag(i__),DRmag(i__),B_R(i__),BTmag(i__),BTc(i__), + e_BTmag(i__),MB_5log_h(i__),Names(i__) if(ar__2(45:49) .EQ. '') RVhel(i__) = iNULL__ if(ar__2(51:53) .EQ. '') e_RVhel(i__) = iNULL__ if(ar__2(55:55) .EQ. '') q_RVhel(i__) = iNULL__ if(ar__2(107:111) .EQ. '') MB_5log_h(i__) = rNULL__ RAdeg_1(i__) = rNULL__ DEdeg_1(i__) = rNULL__ c Derive coordinates RAdeg_1 and DEdeg_1 from input data c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown) if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15. if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4. if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240. if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__) if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60. if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600. if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__) c ..............Just test output........... write(6,3) + HCG_2(i__),m_HCG(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__), + DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),a(i__),b(i__), + MType(i__),T(i__),morph(i__),RVhel(i__),e_RVhel(i__), + q_RVhel(i__),Bmag(i__),q_Bmag(i__),DBmag(i__),Rmag(i__), + q_Rmag(i__),DRmag(i__),B_R(i__),BTmag(i__),BTc(i__), + e_BTmag(i__),MB_5log_h(i__),Names(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'morpho.dat' ! Morphology (1994 paper, table2) C Format for file interpretation 4 format(I3,A1,3X,I5,3X,A4,4X,A69) C Effective file loading open(unit=1,status='old',file= +'morpho.dat') write(6,*) '....Loading file: morpho.dat' do i__=1,210 read(1,'(A92)')ar__3 read(ar__3,4) + HCG_3(i__),m_HCG_1(i__),RVhel_1(i__),MType_1(i__), + Comments(i__) if(ar__3(8:12) .EQ. '') RVhel_1(i__) = iNULL__ c ..............Just test output........... write(6,4) + HCG_3(i__),m_HCG_1(i__),RVhel_1(i__),MType_1(i__), + Comments(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end