Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/138/338 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 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/138/338 Arecibo Legacy Fast ALFA Survey. IX. (Stierwalt+, 2009)
*================================================================================
*The Arecibo Legacy Fast ALFA Survey.
*IX. The Leo region H I catalog, group membership, and the H I mass function for
*the Leo I group.
* Stierwalt S., Haynes M.P., Giovanelli R., Kent B.R., Martin A.M.,
* Saintonge A., Karachentsev I.D., Karachentseva V.E.
* <Astron. J., 138, 338-361 (2009)>
* =2009AJ....138..338S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! HI Candidate Detections
integer*4 nr__
parameter (nr__=549) ! Number of records
character*101 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Seq ! [1/549] Source identification (5-NNN)
* <ALFALFA 5-NNN> in Simbad
integer*4 AGC ! Arecibo General Catalog number
integer*4 RAh ! (h) HI position Right Ascension (J2000) (1)
integer*4 RAm ! (min) HI position Right Ascension (J2000) (1)
real*4 RAs ! (s) HI position Right Ascension (J2000) (1)
character*1 DE_ ! HI position Declination sign (J2000) (1)
integer*4 DEd ! (deg) HI position Declination (J2000) (1)
integer*4 DEm ! (arcmin) HI position Declination (J2000) (1)
integer*4 DEs ! (arcsec) HI position Declination (J2000) (1)
integer*4 RAOh ! (h) ? Optical position Right Ascension (J2000) (2)
integer*4 RAOm ! (min) ? Optical position Right Ascension (J2000) (2)
real*4 RAOs ! (s) ? Optical position Right Ascension (J2000) (2)
character*1 DEO_ ! Optical position Declination sign (J2000) (2)
integer*4 DEOd ! (deg) ? Optical position Declination (J2000) (2)
integer*4 DEOm ! (arcmin) ? Optical position Declination (J2000) (2)
integer*4 DEOs ! (arcsec) ? Optical position Declination (J2000) (2)
integer*4 cz ! (km/s) Heliocentric velocity of HI source (3)
integer*4 W50 ! (km/s) Line profile at 50% level velocity width (4)
integer*4 e_W50 ! (km/s) Error in W50 (4)
real*4 Fc ! (Jy.km/s) Integrated line flux (5)
real*4 S_N ! Signal-to-noise (6)
real*4 rms ! (mJy) Spatially integrated spectral profile noise
* figure (7)
real*4 Dist ! (Mpc) ? Adopted distance (8)
real*4 logM ! ([Msun]) ? Log of the HI mass (9)
integer*4 Code ! [1/9] Object code (10)
character*1 Note ! [*] * indicates a note in notes.dat file
*Note (1): Centroid position after correction for systematic telescope pointing
* errors (see Kent et al. (2008AJ....136..713K) for a description of how
* pointing errors vary with declination for the Arecibo telescope). The
* accuracy of HI positions depends on source strength.
*Note (2): Centroid position of the optical galaxy found to provide the most
* reasonable optical counterpart to the HI detection. Assignments of optical
* identifications are made via the Skyview website and are based on spatial
* proximity, morphology, color, and redshift. Accuracy of centroids is
* estimated to be <=25". For cases with lacking or ambiguous optical
* counterparts, comments are provided as alerted by an asterisk in the
* "Code" column.
*Note (3): Measured as the midpoint between the channels at which the flux
* density drops to 50%. The error can be estimated as half the error on
* the width.
*Note (4): Corrections for broadening but not turbulent motions, disk
* inclination, or cosmological effects are applied. The estimated error on the
* velocity width is estimated by the sum in quadrature of two components: a
* statistical error, principally dependent on the S/N ratio of the feature
* measured, and a systematic error associated with the observer's subjective
* guess at the quality of the chosen spectral extent of the feature. In the
* majority of cases, the statistical error is significantly larger than the
* systematic error; thus the latter is ignored.
*Note (5): This is measured on the integrated spectrum, obtained by spatially
* integrating the source image over a solid angle of at least 7'x7' and
* dividing by the sum of the survey beam values over the same set of image
* pixels (see Shostak & Allen 1980A&A....81..167S).
*Note (6): S/N of the detection, as estimated by
* S/N=(1000Fc/W50){omega}_smo_^1/2^/{sigma}_rms_,
* where Fc is the integrated flux density, as listed in Column 7, the ratio
* of 1000Fc/W50 is the mean flux across the feature in mJy, {omega}_smo_,
* the smoothing width expressed as the number of spectral resolution bins
* of 10km/s bridging half of the signal width, is either W50/(2x10) for
* W50<400km/s or 400/(2x10)=20 for W50>=400km/s, and {sigma}_rms_ is the
* rms noise figure across the spectrum measured in mJy at 10km/s resolution,
* as tabulated in Column "rms".
*Note (7): The noise figure is the r.m.s. as measured over the signal- and
* rfi-free portions of the spectrum, after Hanning smoothing to a spectral
* resolution of 10km/s.
*Note (8): For objects with cz_cmb_>6000, the distance is simply cz_cmb_/H0,
* where cz_cmb_ is the recessional velocity measured in the Cosmic Microwave
* Background reference frame and H0 is the Hubble constant, for which we use
* a value of 70km/s/Mpc. For objects of lower cz_cmb_, we use the
* multiattractor, peculiar velocity model for the local Universe presented in
* Masters (2005, PhD thesis, Cornell Univ). Objects which are thought to be
* parts of clusters or groups (for group membership assignments Springob et al.
* (2007, Cat. J/ApJS/172/599)) are assigned the czcmb of the cluster or group.
* A detailed analysis of group and membership of Leo objects is presented in
* Section 4.
*Note (9): Obtained using M_HI_=2.356x10^5^D^2^_Mpc_F_c_.
*Note (10): object code, defined as follows:
* 1 = refers to sources of S/N and general qualities that make it a reliable
* detection: an approximate S/N threshold of 6.5, a good match between the
* two independent polarizations, and a spatial extent consistent with the
* characteristics of the telescope beam. Thus, some candidate detections
* with S/N 6.5 have been excluded on grounds of polarization mismatch,
* spectral vicinity to RFI features or peculiar spatial properties.
* Likewise, some features of S/N<6.5 are included as reliable detections
* if the source's optical characteristics clearly resemble typical
* galaxies found at the redshift of the HI feature. We estimate that
* detection candidates with S/N<6.5 in Table will be confirmed in
* follow-up observations in better than 95% of cases (Saintonge,
* 2007AJ....133.2087S).
* 2 = refers to sources of low S/N (<6.5), which would ordinarily not be
* considered reliable detections by the criteria set for code 1. However,
* those HI candidate sources are matched with optical counterparts with
* known optical redshifts which, within their respective errors, coincide
* with those measured in the HI line. We refer to these sources
* as "priors."
* 9 = refers to objects assumed to be high-velocity clouds (HVCs) based on
* their low heliocentric velocities (<200km/s) and their lack of an
* optical counterpart; no estimate for their distance is made.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'notes.dat' ! Individual notes
integer*4 nr__1
parameter (nr__1=138) ! Number of records
character*80 ar__1 ! Full-size record
integer*4 Seq_1 ! [1/549] Source identification (5-NNN)
character*74 Note_1 ! Text of the note
C=============================================================================
C Loading file 'table1.dat' ! HI Candidate Detections
C Format for file interpretation
1 format(
+ 2X,I3,1X,I6,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,I2,1X,
+ I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,I5,1X,I3,1X,I3,1X,F5.2,1X,
+ F5.1,1X,F4.2,1X,F5.1,1X,F5.2,1X,I1,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,549
read(1,'(A101)')ar__
read(ar__,1)
+ Seq,AGC,RAh,RAm,RAs,DE_,DEd,DEm,DEs,RAOh,RAOm,RAOs,DEO_,DEOd,
+ DEOm,DEOs,cz,W50,e_W50,Fc,S_N,rms,Dist,logM,Code,Note
if(ar__(35:36) .EQ. '') RAOh = iNULL__
if(ar__(38:39) .EQ. '') RAOm = iNULL__
if(ar__(41:44) .EQ. '') RAOs = rNULL__
if(ar__(47:48) .EQ. '') DEOd = iNULL__
if(ar__(50:51) .EQ. '') DEOm = iNULL__
if(ar__(53:54) .EQ. '') DEOs = iNULL__
if(ar__(87:91) .EQ. '') Dist = rNULL__
if(ar__(93:97) .EQ. '') logM = rNULL__
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)
+ Seq,AGC,RAh,RAm,RAs,DE_,DEd,DEm,DEs,RAOh,RAOm,RAOs,DEO_,DEOd,
+ DEOm,DEOs,cz,W50,e_W50,Fc,S_N,rms,Dist,logM,Code,Note
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'notes.dat' ! Individual notes
C Format for file interpretation
2 format(2X,I3,1X,A74)
C Effective file loading
open(unit=1,status='old',file=
+'notes.dat')
write(6,*) '....Loading file: notes.dat'
do i__=1,138
read(1,'(A80)')ar__1
read(ar__1,2)Seq_1,Note_1
c ..............Just test output...........
write(6,2)Seq_1,Note_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end