Conversion of standardized ReadMe file for
file /./ftp/cats/VII/155 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 2026-May-19
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/155 Third Reference Cat. of Bright Galaxies (RC3) (de Vaucouleurs+ 1991)
*================================================================================
*Third Reference Catalogue of Bright Galaxies (RC3)
* de Vaucouleurs G., de Vaucouleurs A., Corwin H.G., Buta R.J.,
* Paturel G., Fouque P.
* <Springer-Verlag: New York, (1991)>
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rc3' ! The revised RC3 catalog
integer*4 nr__
parameter (nr__=23011) ! Number of records
character*363 ar__ ! Full-size record
C B2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension B2000
real*8 DEdeg ! (deg) Declination B2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 RAh ! (h) *Right Ascension B2000 (hours)
integer*4 RAm ! (min) *Right Ascension B2000 (minutes)
real*4 RAs ! (s) *Right Ascension B2000 (sec. or min.)
character*1 DE_ ! [-+] Sign of declination
integer*4 DEd ! (deg) *Declination B2000 (degrees)
integer*4 DEm ! (arcmin) *Declination B2000 (minutes)
integer*4 DEs ! (arcsec) *? Declination B2000 (seconds)
integer*4 RA1950h ! (h) Right Ascension 1950 (hours)
integer*4 RA1950m ! (min) Right Ascension 1950 (minutes)
real*4 RA1950s ! (s) *Right Ascension 1950 (sec. or min.)
character*1 DE1950_ ! [-+] Sign of declination 1950
integer*4 DE1950d ! (deg) Declination 1950 (degrees)
integer*4 DE1950m ! (arcmin) Declination 1950 (minutes)
integer*4 DE1950s ! (arcsec) *? Declination 1950 (seconds)
real*4 GLON ! (deg) *Galactic longitude
real*4 GLAT ! (deg) *Galactic latitude
real*4 SGLON ! (deg) Supergalactic longitude on RC2 sys.
real*4 SGLAT ! (deg) Supergalactic latitude on RC2 sys.
character*12 name ! *Name
character*15 altname ! *Alternate name
character*14 desig ! Other designation
character*11 PGC ! *PGC number
character*7 type ! *Type
character*5 typesr ! *Source of revised type
real*4 T ! *? T (Hubble stage)
real*4 e_T ! ? Mean error of T
real*4 lumcl ! ? Luminosity class (file "mtypes.tex")
real*4 e_lumcl ! ? Mean error on luminosity class
integer*4 o_lumcl ! ? Num. of luminosity class estimates
real*4 D25 ! ([0.1arcmin]) *? Log D25
character*1 u_D25 ! [?]? indicates Log D25 very uncer.
real*4 e_D25 ! ([0.1arcmin]) *? Mean error on log D25
real*4 R25 ! *? Log R25 (isophotal diameter ratio)
character*1 u_R25 ! [?]? indicates Log R25 very uncertain
real*4 e_R25 ! *? Mean error on log R25
real*4 Do ! ([0.1arcmin]) *? Log Do
real*4 Ae ! ([0.1arcmin]) *? Log Ae (effective aperture)
real*4 e_Ae ! ([0.1arcmin]) ? Mean error on log Ae
integer*4 PA ! (deg) *? Position angle of the major axis
real*4 BT ! (mag) *? BT (total B magnitude)
character*1 BT_code ! *? BT code
real*4 e_BT ! (mag) ? Mean error on BT
real*4 Bmag ! (mag) *? mB (photographic magnitude)
real*4 e_Bmag ! (mag) *? Mean error on Bmag
real*4 BoT ! (mag) *? BoT
real*4 m_25 ! (mag/arcmin2) *? m'25
real*4 e_m_25 ! (mag/arcmin2) ? Mean error m'25
real*4 m_e ! (mag/arcmin2) *? m'e
real*4 e_m_e ! (mag/arcmin2) ? Mean error m'e
real*4 mFIR ! (mag) *? mFIR
real*4 m21 ! (mag) *? m21
real*4 e_m21 ! (mag) ? Mean error on m21
real*4 B_VT ! (mag) *? (B-V)T (total (B-V))
real*4 e_B_VT ! (mag) ? Mean error on (B-V)T
real*4 B_Ve ! (mag) *? (B-V)e (mean B-V within Ae)
real*4 e_B_Ve ! (mag) ? Mean error on (B-V)e
real*4 B_VoT ! (mag) *? (B-V)oT
real*4 U_BT ! (mag) *? (U-B)T (total (U-B))
real*4 e_U_BT ! (mag) ? Mean error on (U-B)T
real*4 U_Be ! (mag) *? (U-B)e
real*4 e_U_Be ! (mag) ? Mean error on (U-B)e
real*4 U_BoT ! (mag) *? (U-B)To
real*4 HI ! (mag) *? HI (neutral hydrogen index)
real*4 Ai ! (mag) *? Ai (internal B extinction)
real*4 A21 ! (mag) *? A21 (HI self absorption)
real*4 Ag ! (mag) *? Ag (galactic extinction in B)
integer*4 W20 ! (km/s) *? W20
integer*4 e_W20 ! (km/s) ? Mean error on W20
integer*4 W50 ! (km/s) *? W50
integer*4 e_W50 ! (km/s) ? Mean error of W50
integer*4 V21 ! (km/s) *? V21 (heliocentric radial velocity)
integer*4 e_V21 ! (km/s) ? Mean error on V21
integer*4 cz ! (km/s) *? Vopt= cz
integer*4 e_cz ! (km/s) ? Mean error on Vopt
integer*4 VGSR ! (km/s) *? VGSR
integer*4 V3K ! (km/s) *? V3K
*Note on RAh, RAm, RAs, DEd, DEm, DEs, RA1950s and DE1950s:
* The right ascensions and declinations for the equinox 2000.0, are
* precessed from the 1950.0 position using the 1976 IAU constants
* (see, e.g., the Supplement to The Astronomical Almanac for 1984).
* These are given to 0.1 second of time and 1 arcsec when available,
* and to 0.1 minute of time and 1 arcmin otherwise (Section 3.1.a,
* page 11).
* The right ascensions originally given with an accuracy of 0.1min
* have been converted to seconds; such converted RAs and RAs1950
* values have no decimal, and the corresponding DEs and DEs1950
* are blank.
*Note on GLON and GLAT:
* The galactic coordinates given to 0.01, are calculated following
* the IAU 1958 prescription (Blaauw et al. 1960) with the North
* Galactic Pole at alpha = 12h49m , delta = +27 24 (1950), and the
* origin at alpha = 17h 42.4m, delta = -28 55 (1950).
*Note on name:
* Names (e.g., LMC, SMC) or NGC and IC designations.
*Note on altname:
* UGC (Nilson 1973), ESO (Lauberts 1982), MCG (Vorontsov-Velyaminov
* et al. 1962-1974), UGCA (Nilson 1974), and CGCG (Zwicky et al.
* 1961-1968) designations, given in that order of preference. MCG
* designations not listed here are given in UGC and ESO.
*Note on PGC:
* PGC (Paturel et al. 1989a,b) designation. For cross identifications
* of various catalogues with the PGC, see Appendix 10, page 561.
*Note on type:
* Mean revised morphological type in the RC2 system, coded as in RC2
* (Section 3.3.a, page 13; see also file "mtypes.tex").
*Note on typesr:
* "V" indicates a classification by de Vaucouleurs; "4" indicates a
* ring galaxy classified by Buta on 4-meter plates.
*Note on T:
* Mean numerical index of stage along the Hubble sequence in RC2
* system (coded as explained in Section 3.3.c, page 16;
* see also file "mtypes.tex").
*Note on D25:
* Mean decimal logarithm of the apparent major isophotal diameter
* measured at or reduced to surface brightness level muB = 25.0 B/mag2,
* as explained in Section 3.4.a, page 21. The unit of D is 0.1
* arcmin to avoid negative entries.
*Note on e_D25:
* A "?" indicates a large, uncertain error.
*Note on R25:
* Mean decimal logarithm of the ratio of the major isophotal
* diameter, D25, to the minor isophotal diameter, D25, measured at or
* reduced to the surface brightness level muB = 25.0 B/mag2, and its
* mean error as explained in Section 3.4.b, page 26.
*Note on e_R25:
* A "?" indicates a large, uncertain error.
*Note on Ae:
* Decimal logarithm of the apparent diameter (in 0.1 arcmin) of the
* "effective aperture," the circle centered on the nucleus within
* which one-half of the total B-band flux is emitted, derived as
* explained in Section 3.4.c, page 28.
*Note on Do:
* Decimal logarithm of the isophotal major diameter corrected to
* "face-on" (i = 0), and corrected for galactic extinction to Ag = 0,
* but not for redshift, as explained in Section 3.4.d, page 29.
*Note on PA:
* Position angle, measured in degrees from north through east (all
* <180), taken when available from UGC, ESO, and ESGC (and in a few
* cases from H I data) (Section 3.5.a, page 30).
*Note on Ag:
* Galactic extinction in B-band magnitudes, calculated following
* Burstein and Heiles (1978a,b, 1982, 1984) as explained in Section
* 3.5.b, page 30.
*Note on Ai:
* Internal extinction in B-band magnitudes (for correction to
* face-on), calculated from log R and T as explained in Section
* 3.5.c, page 31.
*Note on A21:
* A21 = H I line self-absorption in magnitudes (for correction to
* face-on), calculated from log R and T greater or equal to 1 as
* explained in Section 3.5.d, page 32.
*Note on BT:
* BT = total (asymptotic) magnitude in the B system derived by
* extrapolation from photoelectric aperture- magnitude data, BAT, and
* from surface photometry with photoelectric zero point, BSTas
* explained in Section 3.6.a, page 32.
*Note on BT_code:
* "M" = BT is the weighted mean of BAT and BST,
* "S" = only surface photometry used
* "V" = BT is a V-band magnitude rather than a B-band magnitude
* "v" = the nucleus of the galaxy is variable.
* "*" indicates that deriving BAT would have required an
* extrapolation in excess of 0.75 mag.
*Note on Bmag and e_Bmag:
* Photographic magnitude and its mean error from Ames (1930), Shapley
* and Ames (1932), CGCG, Buta and Corwin (1986), and/or Lauberts and
* Valentijn (1989) reduced to the BT system as explained in Section
* 3.6.b, page 37.
*Note on mFIR:
* Calculated from mFIR = -20.0 - 2.5 logFIR, where FIR is the far
* infrared continuum flux measured at 60 and 100 microns as listed in
* the IRAS Point Source Catalog (1987) <II/125>. For galaxies larger
* than 80 in RC2 and for the Virgo cluster area, resolved by the IRAS
* beam, integrated fluxes are taken from Rice et al.(1988) or Helou et
* al.(1988). See Section 3.6.c, page 43, for details.
*Note on BoT:
* Total "face-on" magnitude corrected for galactic and internal
* extinction, and for redshift as explained in Section 3.6.d, page
* 44.
*Note on B-VT and U-BT:
* Total (asymptotic) color index in the Johnson B-V and U-B system,
* respectively, derived by extrapolation from photoelectric
* color-aperture data, and/or from surface photometry with a
* photoelectric zero point as explained in Section 3.7.a, page 45.
*Note on B-VoT and U-BoT:
* Total B-V and U-B color indices, respectively, corrected for
* galactic and internal extinction, and for redshift, as explained in
* Section 3.7.b, page 47.
*Note on B-Ve and U-Be:
* Mean B-V and U-B color indices, respectively, within the effective
* aperture Ae, derived by interpolation from photoelectric
* color-aperture data as explained in Section 3.7.a, page 45.
*Note on m'e:
* Mean B-band surface brightness in magnitudes per square arcmin
* (B-m/sm) within the effective aperture Ae, and its mean error,
* calculated by the relation m'e= BT + 0.75 + 5 logAe- 5.26. This m'e
* is statistically related to the effective mean surface brightness,
* mu'e (RC2, p. 31; Olson and de Vaucouleurs 1981), with which it
* coincides when log R = 0 (i = 0) (Section 3.8.a, page 49).
*Note on m'25:
* The mean surface brightness in magnitudes per square arcmin
* (B-m/sm) within the muB = 25.0 B-m/ss elliptical isophote of major
* axis log D25 and axis ratio log R25, defined as in RC2 (Equation
* 21) by:
* m'25= BT + delta(m25) + 5 logD25- 2.5 logR25- 5.26,
* where delta(m25) = 2.5 log LT/L25 = B25 - BT and is the magnitude
* increment contributed by the outer regions of a galaxy fainter than
* muB = 25.0B-m/ss and, its mean error. For details, see Section
* 3.8.b, page 50.
*Note on m21:
* 21-cm emission line magnitude defined by m21 = 21.6 - 2.5 log SH ,
* where SH is the measured neutral hydrogen flux density in units of
* 10-24W/m2. For details, see Section 3.9.a, page 51.
*Note on W20 and W50:
* Neutral hydrogen line full width (in km/s) measured at the 20%
* level (I20/Imax) and the 50% level (I50/Imax), respectively, as
* explained in Section 3.9.b, page 51.
*Note on HI:
* Corrected neutral hydrogen index, which is the difference
* mo21-BoT between the corrected (face-on) 21-cm emission line
* magnitude and the similarly corrected magnitude in the BT system.
* Details are given in Section 3.9.c, page 52. Since m21 and BT are
* listed separately in columns 6 and 9, record 1, there is no need to
* print the uncorrected index.
*Note on V21:
* The mean heliocentric radial velocity derived from neutral hydrogen
* observations, as explained in Section 3.10.a, page 52.
*Note on cz:
* The mean heliocentric radial velocity, derived from optical
* observations, as explained in Section 3.10.b, page 53.
*Note on VGSR:
* The weighted mean of the neutral hydrogen and optical velocities,
* corrected to the "Galactic standard of rest", as explained in
* Section 3.10.c, page 54.
*Note on V3K:
* The weighted mean velocity corrected to the reference frame defined
* by the 3K microwave background radiation, as explained in Section
* 3.10.d, page 55.
C=============================================================================
C Loading file 'rc3' ! The revised RC3 catalog
C Format for file interpretation
1 format(
+ I2,I2,F4.1,1X,A1,I2,I2,I2,1X,I2,I2,F4.1,1X,A1,I2,I2,I2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,A12,A15,1X,A14,1X,A11,1X,A7,
+ 1X,A5,1X,F4.1,1X,F3.1,1X,F4.1,1X,F3.1,1X,I1,1X,F4.2,A1,1X,
+ F3.2,1X,F4.2,A1,1X,F3.2,1X,F4.2,1X,F4.2,1X,F3.2,1X,I3,1X,F5.2,
+ A1,1X,F3.2,1X,F5.2,1X,F3.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
+ F3.2,1X,F5.2,1X,F5.2,1X,F2.1,1X,F4.2,1X,F3.2,1X,F4.2,1X,F3.2,
+ 1X,F4.2,1X,F5.2,1X,F3.2,1X,F4.2,1X,F3.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F3.2,1X,F4.2,1X,I3,1X,I2,1X,I3,1X,I2,1X,I5,1X,I2,1X,
+ I5,1X,I3,1X,I5,1X,I5)
C Effective file loading
open(unit=1,status='old',file=
+'rc3')
write(6,*) '....Loading file: rc3'
do i__=1,23011
read(1,'(A363)')ar__
read(ar__,1)
+ RAh,RAm,RAs,DE_,DEd,DEm,DEs,RA1950h,RA1950m,RA1950s,DE1950_,
+ DE1950d,DE1950m,DE1950s,GLON,GLAT,SGLON,SGLAT,name,altname,
+ desig,PGC,type,typesr,T,e_T,lumcl,e_lumcl,o_lumcl,D25,u_D25,
+ e_D25,R25,u_R25,e_R25,Do,Ae,e_Ae,PA,BT,BT_code,e_BT,Bmag,
+ e_Bmag,BoT,m_25,e_m_25,m_e,e_m_e,mFIR,m21,e_m21,B_VT,e_B_VT,
+ B_Ve,e_B_Ve,B_VoT,U_BT,e_U_BT,U_Be,e_U_Be,U_BoT,HI,Ai,A21,Ag,
+ W20,e_W20,W50,e_W50,V21,e_V21,cz,e_cz,VGSR,V3K
if(ar__(15:16) .EQ. '') DEs = iNULL__
if(ar__(32:33) .EQ. '') DE1950s = iNULL__
if(ar__(132:135) .EQ. '') T = rNULL__
if(ar__(137:139) .EQ. '') e_T = rNULL__
if(ar__(141:144) .EQ. '') lumcl = rNULL__
if(ar__(146:148) .EQ. '') e_lumcl = rNULL__
if(ar__(150:150) .EQ. '') o_lumcl = iNULL__
if(ar__(152:155) .EQ. '') D25 = rNULL__
if(ar__(158:160) .EQ. '') e_D25 = rNULL__
if(ar__(162:165) .EQ. '') R25 = rNULL__
if(ar__(168:170) .EQ. '') e_R25 = rNULL__
if(ar__(172:175) .EQ. '') Do = rNULL__
if(ar__(177:180) .EQ. '') Ae = rNULL__
if(ar__(182:184) .EQ. '') e_Ae = rNULL__
if(ar__(186:188) .EQ. '') PA = iNULL__
if(ar__(190:194) .EQ. '') BT = rNULL__
if(ar__(197:199) .EQ. '') e_BT = rNULL__
if(ar__(201:205) .EQ. '') Bmag = rNULL__
if(ar__(207:209) .EQ. '') e_Bmag = rNULL__
if(ar__(211:215) .EQ. '') BoT = rNULL__
if(ar__(217:221) .EQ. '') m_25 = rNULL__
if(ar__(223:226) .EQ. '') e_m_25 = rNULL__
if(ar__(228:232) .EQ. '') m_e = rNULL__
if(ar__(234:236) .EQ. '') e_m_e = rNULL__
if(ar__(238:242) .EQ. '') mFIR = rNULL__
if(ar__(244:248) .EQ. '') m21 = rNULL__
if(ar__(250:251) .EQ. '') e_m21 = rNULL__
if(ar__(253:256) .EQ. '') B_VT = rNULL__
if(ar__(258:260) .EQ. '') e_B_VT = rNULL__
if(ar__(262:265) .EQ. '') B_Ve = rNULL__
if(ar__(267:269) .EQ. '') e_B_Ve = rNULL__
if(ar__(271:274) .EQ. '') B_VoT = rNULL__
if(ar__(276:280) .EQ. '') U_BT = rNULL__
if(ar__(282:284) .EQ. '') e_U_BT = rNULL__
if(ar__(286:289) .EQ. '') U_Be = rNULL__
if(ar__(291:293) .EQ. '') e_U_Be = rNULL__
if(ar__(295:298) .EQ. '') U_BoT = rNULL__
if(ar__(300:304) .EQ. '') HI = rNULL__
if(ar__(306:309) .EQ. '') Ai = rNULL__
if(ar__(311:313) .EQ. '') A21 = rNULL__
if(ar__(315:318) .EQ. '') Ag = rNULL__
if(ar__(320:322) .EQ. '') W20 = iNULL__
if(ar__(324:325) .EQ. '') e_W20 = iNULL__
if(ar__(327:329) .EQ. '') W50 = iNULL__
if(ar__(331:332) .EQ. '') e_W50 = iNULL__
if(ar__(334:338) .EQ. '') V21 = iNULL__
if(ar__(340:341) .EQ. '') e_V21 = iNULL__
if(ar__(343:347) .EQ. '') cz = iNULL__
if(ar__(349:351) .EQ. '') e_cz = iNULL__
if(ar__(353:357) .EQ. '') VGSR = iNULL__
if(ar__(359:363) .EQ. '') V3K = iNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ RAh,RAm,RAs,DE_,DEd,DEm,DEs,RA1950h,RA1950m,RA1950s,DE1950_,
+ DE1950d,DE1950m,DE1950s,GLON,GLAT,SGLON,SGLAT,name,altname,
+ desig,PGC,type,typesr,T,e_T,lumcl,e_lumcl,o_lumcl,D25,u_D25,
+ e_D25,R25,u_R25,e_R25,Do,Ae,e_Ae,PA,BT,BT_code,e_BT,Bmag,
+ e_Bmag,BoT,m_25,e_m_25,m_e,e_m_e,mFIR,m21,e_m21,B_VT,e_B_VT,
+ B_Ve,e_B_Ve,B_VoT,U_BT,e_U_BT,U_Be,e_U_Be,U_BoT,HI,Ai,A21,Ag,
+ W20,e_W20,W50,e_W50,V21,e_V21,cz,e_cz,VGSR,V3K
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end