Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/117/2039 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/117/2039 Sc galaxies I photometry, 21cm data (Haynes+, 1999)
*================================================================================
*The I-band Tully-Fisher relation for Sc galaxies: 21 centimeter H I line data.
* Haynes M.P., Giovanelli R., Chamaraux P., Da Costa L.N., Freudling W.,
* Salzer J.J., Wegner G.
* <Astron. J. 117, 2039 (1999)>
* =1999AJ....117.2039H (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Properties of galaxies derived from spectra in
digital archive
integer*4 nr__
parameter (nr__=1201) ! Number of records
character*105 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 UGC_AGC (nr__) ! UGC (Cat. <VII/26>) number, where applicable,
* or else in our private database, referred to
* as the Arecibo General Catalog (AGC) (1)
character*15 Name (nr__) ! Other galaxy designation (2)
integer*4 RAh (nr__) ! (h) Right ascension (1950)
integer*4 RAm (nr__) ! (min) Right ascension (1950)
real*4 RAs (nr__) ! (s) Right ascension (1950)
character*1 DE_ (nr__) ! Declination sign (1950)
integer*4 DEd (nr__) ! (deg) Declination (1950)
integer*4 DEm (nr__) ! (arcmin) Declination (1950)
integer*4 DEs (nr__) ! (arcsec) Declination (1950)
real*4 a (nr__) ! (arcmin) Blue major diameter
real*4 b (nr__) ! (arcmin) Blue minor diameter
integer*4 TType (nr__) ! Morphological type code following the RC3
* (Cat. <VII/155>) system
real*4 S21cm (nr__) ! (Jy.km/s) Observed integrated 21 cm H I line flux
real*4 S21cmc (nr__) ! (Jy.km/s) Corrected integrated 21 cm H I line flux
real*4 rms (nr__) ! (mJy) rms noise per channel of the spectrum
real*4 S_N (nr__) ! Emission profile signal-to-noise ratio taken
* as the ratio of the peak flux to the rms noise
integer*4 Vhel (nr__) ! (km/s) Heliocentric velocity of the HI line signal
integer*4 W21cm (nr__) ! (km/s) Full velocity width of the HI line
real*4 Wc (nr__) ! (km/s) Corrected velocity width
real*4 e_Wc (nr__) ! (km/s) Estimated error in Wc
character*4 Code (nr__) ! Codes indicating the correlator mode,
* smoothing, width quality, and data source (3)
character*1 Note (nr__) ! [*] Note detailed in notes.dat file
*Note (1): 14 to 12901, UGC number; 20471 to 630257, AGC number
*Note (2): MCG number were corrected from a bad format in the original version.
* For AGC 530053 (MCG -33-40-14), the coordinates correspond to
* MGC-03-34-014, which is not AGC 530053
* For UGC 12355 (MCG +25-80-63), the coordinates do not correspond
* to MCG+02-58-063
*Note (3): The first code refers to the origin and configuration of the spectrum
* from which parameters are derived:
* (a) Arecibo spectrum, 20 MHz/512 channels
* (b) Arecibo spectrum, 10 MHz/252 channels
* (c) Arecibo spectrum, 10 MHz/512 channels
* (d) Arecibo spectrum, 10 MHz/504 channels
* (e) Arecibo spectrum, 5 MHz/512 channels
* (f) Effelsberg spectrum, 25 MHz/512 channels
* (g) Effelsberg spectrum, 12.5 MHz/512 channels
* (h) Green Bank 91 m spectrum, 10 MHz/192 channels
* (i) Green Bank 91 m spectrum, 5 MHz/192 channels
* (j) Green Bank 43 m spectrum, 20 MHz/512 channels
* (k) Green Bank 43 m spectrum, 10 MHz/512 channels
* (l) Green Bank 43 m spectrum, 5 MHz/512 channels
* (m) Nangay spectrum, 6.4 MHz/256 channels
* The second code is the smoothing code:
* (H) single Hanning only
* (B) Hanning plus three-channel boxcar
* The third is a qualitative assessment of the quality of the profile for
* TF applications:
* (G) good
* (F) fair
* (S) single peak
* (P) poor quality for TF applications
* (M) marginal detection
* (C) confused
* Velocity widths measured from profiles classified as good detections
* should be useful for TF applications; those measured on profiles
* denoted "F" should be used with caution. Because the width measuring
* algorithm is designed for application to two-horned profiles, the
* widths measured on single peaked profiles should also be used with
* caution, as they may underestimate the true rotation width.
* Widths in the last three categories cannot be used for TF purposes. A
* designation of marginal detection is given to those cases in which
* the signal has not been verified through adequate reobservation.
* Confusion is identified only when contamination from another galaxy
* in the beam is certain.
* Another code follows where applicable that indicates inclusion in
* other compilations:
* "c" for the cluster compilation of Giovanelli et al. (1997,
* Cat. <J/AJ/113/22>)
* "a" for the study of asymmetry in H I profiles by Haynes et al.
* (1998AJ....115...62H)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'notes.dat' ! Individual notes of table1.dat
integer*4 nr__1
parameter (nr__1=174) ! Number of records
character*80 ar__1 ! Full-size record
character*1 n_UGC_AGC (nr__1) ! [AU] A: AGC; U: UGC
integer*4 UGC_AGC_1 (nr__1) ! UGC or AGC number
character*72 Note_1 (nr__1) ! Text of the note
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Properties of galaxies derived from data in
literature
integer*4 nr__2
parameter (nr__2=72) ! Number of records
character*64 ar__2 ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__2) ! (deg) Right Ascension 1950
real*8 DEdeg_1 (nr__2) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 UGC_AGC_2 (nr__2) ! UGC (Cat. <VII/26>) number, where applicable,
* or else in our private database, referred to
* as the Arecibo General Catalog (AGC) (1)
character*8 Name_1 (nr__2) ! Other name
integer*4 RAh_1 (nr__2) ! (h) Right ascension (1950)
integer*4 RAm_1 (nr__2) ! (min) Right ascension (1950)
real*4 RAs_1 (nr__2) ! (s) Right ascension (1950)
character*1 DE__1 (nr__2) ! Declination sign (1950)
integer*4 DEd_1 (nr__2) ! (deg) Declination (1950)
integer*4 DEm_1 (nr__2) ! (arcmin) Declination (1950)
integer*4 DEs_1 (nr__2) ! (arcsec) Declination (1950)
real*4 a_1 (nr__2) ! (arcmin) Blue major diameter
real*4 b_1 (nr__2) ! (arcmin) Blue minor diameter
integer*4 TType_1 (nr__2) ! RC3 (Cat. <VII/155>) morphological type code
integer*4 Vhel_1 (nr__2) ! (km/s) Heliocentric velocity
real*4 Wc_1 (nr__2) ! (km/s) Corrected velocity width
integer*4 Ref (nr__2) ! References (2)
*Note (1): 14 to 12901, UGC number; 20471 to 630257, AGC number
*Note (2): the references are:
* 1: Theureau et al., 1998, Cat. <J/A+AS/130/333>
* 2: Fisher & Tully, 1981ApJS...47..139F
* 3: These observations but digital spectrum lost
* 4: Richter & Huchtmeier, 1989; General Catalog of H I Observations
* of Galaxies
* 5: Haynes & Giovanelli, 1984AJ.....89..758H
* 6: Eder et al., 1991AJ....102..572E
* 7: Mould et al., 1995ApJS...96....1M
* 8: Mathewson & Ford, 1996, Cat. <J/ApJS/107/97>
* 9: Fouque et al., 1990A&AS...86..473F
* 10: Magri, 1994, Cat. <J/AJ/108/896>
* 11: Gavazzi et al., 1996A&A...312..397G
* 12: Gavazzi, 1989ApJ...346...59G
* 13: Giovanelli & Haynes, 1985ApJ...292..404G
C=============================================================================
C Loading file 'table1.dat' ! Properties of galaxies derived from spectra in
* digital archive
C Format for file interpretation
1 format(
+ I6,1X,A15,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F3.1,1X,
+ F3.1,1X,I1,1X,F5.2,1X,F6.2,1X,F5.2,1X,F5.1,1X,I5,1X,I3,1X,
+ F5.1,1X,F4.1,1X,A4,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,1201
read(1,'(A105)')ar__
read(ar__,1)
+ UGC_AGC(i__),Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),a(i__),b(i__),TType(i__),
+ S21cm(i__),S21cmc(i__),rms(i__),S_N(i__),Vhel(i__),W21cm(i__),
+ Wc(i__),e_Wc(i__),Code(i__),Note(i__)
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)
+ UGC_AGC(i__),Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),a(i__),b(i__),TType(i__),
+ S21cm(i__),S21cmc(i__),rms(i__),S_N(i__),Vhel(i__),W21cm(i__),
+ Wc(i__),e_Wc(i__),Code(i__),Note(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 'notes.dat' ! Individual notes of table1.dat
C Format for file interpretation
2 format(A1,I6,1X,A72)
C Effective file loading
open(unit=1,status='old',file=
+'notes.dat')
write(6,*) '....Loading file: notes.dat'
do i__=1,174
read(1,'(A80)')ar__1
read(ar__1,2)n_UGC_AGC(i__),UGC_AGC_1(i__),Note_1(i__)
c ..............Just test output...........
write(6,2)n_UGC_AGC(i__),UGC_AGC_1(i__),Note_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Properties of galaxies derived from data in
* literature
C Format for file interpretation
3 format(
+ I6,2X,A8,1X,I2,1X,I2,1X,F4.1,2X,A1,I2,1X,I2,1X,I2,1X,F3.1,1X,
+ F3.1,1X,I2,1X,I5,1X,F5.1,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,72
read(1,'(A64)')ar__2
read(ar__2,3)
+ UGC_AGC_2(i__),Name_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),a_1(i__),b_1(i__),
+ TType_1(i__),Vhel_1(i__),Wc_1(i__),Ref(i__)
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)
+ UGC_AGC_2(i__),Name_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),a_1(i__),b_1(i__),
+ TType_1(i__),Vhel_1(i__),Wc_1(i__),Ref(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end