Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/105/1271 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-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. J/AJ/105/1271 Pisces-Perseus supercluster. VI (Giovanelli+, 1993)
*================================================================================
*A survey of the Pisces-Perseus supercluster. VI. The declination zone
*+15.5 deg to +21.5 deg
* Giovanelli R., Haynes M.P.
* <Astron. J. 105, 1271 (1993)>
* =1993AJ....105.1271G (SIMBAD/NED Reference)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'notes1.dat' ! Notes to table1
integer*4 nr__
parameter (nr__=249) ! Number of records
character*79 ar__ ! Full-size record
character*8 Galaxy (nr__) ! Galaxy name, either CGCG or U followed by UGC
* number, repeated if note spans on lines.
character*71 Text (nr__) ! Text of Note
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Parameters of observed objects
integer*4 nr__1
parameter (nr__1=765) ! Number of records
character*129 ar__1 ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__1) ! (deg) Right Ascension 1950
real*8 DEdeg (nr__1) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*8 CGCG (nr__1) ! Identification in the Catalog of Galaxies and
* Clusters of Galaxies (Zwicky et al. 1961-68:
* (CGCG)), that is the field number followed by
* the entry number within that field
integer*4 UGC (nr__1) ! [1/12921]? Uppsala General Catalog (UGC)
* entry number (Nilson 1973)
character*1 m_UGC (nr__1) ! UGC component designation
character*1 n_NGC_IC (nr__1) ! New General Catalogue (NGC) or Index Catalog
* (IC) indicator, 'N' or 'I', respectively
character*5 NGC_IC (nr__1) ! NGC/IC number
integer*4 RAh (nr__1) ! (h) Right ascension 1950
integer*4 RAm (nr__1) ! (min) Right ascension 1950
real*4 RAs (nr__1) ! (s) Right ascension 1950
character*1 DE_ (nr__1) ! Declination sign
integer*4 DEd (nr__1) ! (deg) Declination 1950
integer*4 DEm (nr__1) ! (arcmin) Declination 1950
integer*4 DEs (nr__1) ! (arcsec) Declination 1950
character*3 T (nr__1) ! Morphological type code index (1)
real*4 a (nr__1) ! (arcmin) Major diameter from the UGC where available
* or else measured by eye from the PSS prints
real*4 b (nr__1) ! (arcmin) Minor diameter
integer*4 i (nr__1) ! (deg) [0/180[? Inclination derived from the axial
* ratio following Haynes & Giovanelli (HG84).
real*4 m_z (nr__1) ! (mag) Magnitude given in the CGCG or the UGC.
real*4 m_c (nr__1) ! (mag) Corrected magnitude, after applying the
* corrections for galactic and internal
* extinction, redshift and systematic effects
* as discussed by HG84.
integer*4 RV (nr__1) ! (km/s) 21 cm line heliocentric velocity, measured
* as the midpoint of the emission profile at a
* level of 50% of the mean signal intensity.
integer*4 V_0 (nr__1) ! (km/s) ? Adopted Local Group velocity (5)
integer*4 W1 (nr__1) ! (km/s) ? Observed 21 cm profile width, measured at
* a level of 50% of the mean signal intensity.
integer*4 W2 (nr__1) ! (km/s) ? Observed 21 cm profile width, measured at
* a level of 20% of the peak signal intensity.
integer*4 Wc (nr__1) ! (km/s) ? 21 cm profile width corrected for viewing
* inclination and for redshift broadening
real*4 Fobs (nr__1) ! (Jy.km/s) ? Observed 21 cm line flux integral
real*4 Fc (nr__1) ! (Jy.km/s) ? 21 cm line flux integral, after
* correction for random pointing errors,
* beam dilution, and HI internal absorption
* following HG84.
real*4 Sigma (nr__1) ! (mJy) ? rms noise per channel measured in signal
* and interference-free portions of the spectrum
real*4 snr (nr__1) ! ? Signal-to-noise parameter of HI emission
* spectrum, measured as the ratio of peak
* signal flux to rms noise, after smoothing.
real*4 log_L (nr__1) ! ([solLum]) ? Logarithm of optical luminosity,
* log10(h^2 L), derived from m(c) and scaled
* by the dimensionless Hubble parameter,
* h = H0/100 km/s/Mpc.
real*4 log_M_H_ (nr__1) ! ([solMass]) ? Logarithm of HI mass , log10(h^2 M(H)).
* derived from Fc and V(0).
character*1 Tel (nr__1) ! [017]? Telescope identification code for HI
* observations (2)
integer*4 Q (nr__1) ! [0/5]? Quality code for HI detection (3)
integer*4 Ref (nr__1) ! [1/5]? Reference code for HI data. (4)
character*1 Note (nr__1) ! Note flag (see file notes1.dat)
*Note (1): Morphological type code index:
* 0 - Elliptical
* 1 - SO
* 3 - Sa
* 5 - Sb
* 7 - Sc
* 9 - Irregular
* 10 - Peculiar
* 12 - Generic spiral (S...)
* 13 - Compact, unresolved morphology
* 14 - Compact spiral or irregular, with intermediate types allowed
* A "B" indicates the presence of a bar.
*Note (2): Telescope identification code:
* 0 - Arecibo 305 m with dual circular feed
* 1 - 300-foot telescope at Green Bank
* 7 - Arecibo 305 m with flat (linear polarization) feed
*Note (3): Quality code:
* 0 - Object considered to be undetected at 21 cm. 21 cm line
* observations were performed, so we list the rms noise per
* channel of the observations. But no signal was detected; hence
* no flux is given.
* 1 - High quality detection
* 2 - Marginal to poor detection
* 3 - Flux and velocity ok, but profile unsuitable for use in
* T-F relation
* 4 - HI seen in absorption
* 5 - Confused with neighbor
*Note (4): The HI data reported here come from one of the following sources:
* 1 - This survey
* 2 - Isolated galaxy survey (HG84)
* 3 - Survey of large angular diameter galaxies (Hewitt et al. 1983)
* 4 - Galaxy observed by the authors and I. Karachentsev as part of
* the thin galaxy survey based on the catalog of Karachentsev
* (1989); visual magnitude estimate on PSS.
* 5 - Not explained in the article
*Note (5): Adopted Local Group velocity, assuming a correction
* for the Sun's motion of 300 sin l cos b, where l and b are the
* galactic longitude and latitude. The adopted heliocentric velocity is
* the 21 cm velocity if the object is detected When the object was
* considered to be undetected at 21 cm (Q=0), the velocity corresponds
* to the frequency of the central channel of the spectrometer
C=============================================================================
C Loading file 'notes1.dat' ! Notes to table1
C Format for file interpretation
1 format(A8,A71)
C Effective file loading
open(unit=1,status='old',file=
+'notes1.dat')
write(6,*) '....Loading file: notes1.dat'
do i__=1,249
read(1,'(A79)')ar__
read(ar__,1)Galaxy(i__),Text(i__)
c ..............Just test output...........
write(6,1)Galaxy(i__),Text(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table1.dat' ! Parameters of observed objects
C Format for file interpretation
2 format(
+ A8,I5,A1,A1,A5,I2,I2,F4.1,1X,A1,I2,I2,I2,1X,A3,F5.2,F5.2,I3,
+ 1X,F4.1,1X,F5.2,I6,I6,I4,I4,I4,F6.2,F6.2,F5.2,1X,F5.2,F5.2,
+ F6.2,1X,A1,1X,I1,1X,I1,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,765
read(1,'(A129)')ar__1
read(ar__1,2)
+ CGCG(i__),UGC(i__),m_UGC(i__),n_NGC_IC(i__),NGC_IC(i__),
+ RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),T(i__),a(i__),b(i__),i(i__),m_z(i__),m_c(i__),
+ RV(i__),V_0(i__),W1(i__),W2(i__),Wc(i__),Fobs(i__),Fc(i__),
+ Sigma(i__),snr(i__),log_L(i__),log_M_H_(i__),Tel(i__),Q(i__),
+ Ref(i__),Note(i__)
if(ar__1(9:13) .EQ. '') UGC(i__) = iNULL__
if(ar__1(51:53) .EQ. '') i(i__) = iNULL__
if(ar__1(71:76) .EQ. '') V_0(i__) = iNULL__
if(ar__1(77:80) .EQ. '') W1(i__) = iNULL__
if(ar__1(81:84) .EQ. '') W2(i__) = iNULL__
if(ar__1(85:88) .EQ. '') Wc(i__) = iNULL__
if(ar__1(89:94) .EQ. '') Fobs(i__) = rNULL__
if(ar__1(95:100) .EQ. '') Fc(i__) = rNULL__
if(ar__1(101:105) .EQ. '') Sigma(i__) = rNULL__
if(ar__1(107:111) .EQ. '') snr(i__) = rNULL__
if(ar__1(112:116) .EQ. '') log_L(i__) = rNULL__
if(ar__1(117:122) .EQ. '') log_M_H_(i__) = rNULL__
if(ar__1(126:126) .EQ. '') Q(i__) = iNULL__
if(ar__1(128:128) .EQ. '') Ref(i__) = iNULL__
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,2)
+ CGCG(i__),UGC(i__),m_UGC(i__),n_NGC_IC(i__),NGC_IC(i__),
+ RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),T(i__),a(i__),b(i__),i(i__),m_z(i__),m_c(i__),
+ RV(i__),V_0(i__),W1(i__),W2(i__),Wc(i__),Fobs(i__),Fc(i__),
+ Sigma(i__),snr(i__),log_L(i__),log_M_H_(i__),Tel(i__),Q(i__),
+ Ref(i__),Note(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end