Conversion of standardized ReadMe file for
catalog 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-May-20
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/50 CO Radial Velocities Toward Galactic H II Regions (Blitz+ 1982)
*================================================================================
*Catalog of CO Radial Velocities toward Galactic H II Regions
* Blitz L., Fich M., Stark A.A.
* <Astrophys. Journ. Suppl., 49, 183 (1982)>
* =1982ApJS...49..183B
*Catalog of Carbon Monoxide Observations of Southern Hemisphere H II Regions
* Gillespie A.R., Huggins P.J., Sollner T.C.L.G., Phillips T.G.,
* Gardner F.F., Knowles S.H.
* <Astron. and Astrophys., 60, 221 (1977)>
* =1977A&A....60..221G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cosouth.dat' ! CO at Southern Hemisphere H II Regions
integer*4 nr__
parameter (nr__=37) ! Number of records
character*138 ar__ ! Full-size record
C B1950 position composed of: RAh RAm RAs DE- DEd DEm
real*8 RAdeg ! (deg) Right Ascension B1950
real*8 DEdeg ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
character*14 Sname ! Source designation (1)
integer*4 RAh ! (h) Right Ascension (Equatorial coord. B1950)
integer*4 RAm ! (min) Right Ascension Minutes
integer*4 RAs ! (s) Right Ascension Seconds
character*1 DE_ ! Sign, always '-', of Declination
integer*4 DEd ! (deg) Declination Degrees of arc (equinox B1950)
real*4 DEm ! (arcmin) Declination Minutes of arc
real*8 GLON ! (deg) Galactic longitude, l(II) of the observation
real*8 GLAT ! (deg) Galactic latitude, b(II) of the observation
real*4 VLSR ! (km/s) Median CO velocity (Local Standard of Rest)
real*4 EW ! (km/s) Equivalent line width, in km/s (2)
real*4 TAp ! (K) Peak antenna temperature (3)
real*4 TAc ! (K) Continuum antenna temperature at 5GHz (4)
real*8 V_H109 ! (km/s) ? Velocity of the H 109{alpha} line (5)
real*8 V_OH1 ! (km/s) ? Velocity of the OH line (5)
real*8 V_OH2 ! (km/s) ? Velocity of another strong OH line (5)
real*8 V_OH3 ! (km/s) ? Velocity of another strong OH line (5)
real*8 V_H2CO1 ! (km/s) ? Velocity of the H2CO line (5)
real*4 V_H2CO2 ! (km/s) ? Velocity of another strong H2CO line (5)
character*1 u_V_H2CO2 ! [) ] ')' - uncertainty about the V-H2CO2 data
real*4 Dist ! (kpc) ?=0.0 Distance (6)
real*4 e_Dist ! (kpc) ?=0 Error of distance measurement
character*5 Com ! [A-E, ] Remarks, see "notes.dat" (7)
*Note (1): Source designation.
* Names have been given for 24 of the 37 sources appearing in this file.
* NGC = Sulentic and Tifft 1973 <VII/1>,
* RCW = Rogers, Campbell, and Whiteoak (1960MNRAS.121..103R),
* W = Westerhout (1958BAN....14..215W).
*Note (2):
* Defined as (Integration of TA* over dv)/TA*(max) where TA* is the
* corrected antenna temperature, bytes 63-66.
*Note (3):
* Peak antenna temperature, TA*, corrected for atmospheric extinction
* and telescope efficiency.
*Note (4):
* Continuum antenna temperature at 5 GHz, Tc. This value was taken from
* Whiteoak and Gardner (1974A&A....37..389W).
*Note (5):
* Comparison values for the velocities of the H 109{alpha} line and the
* strongest OH and H2CO lines (Wilson, et al. 1970A&A.....6..364W,
* Whiteoak and Gardner 1974A&A....37..389W,
* Caswell and Robinson 1974AuJPh..27..597C).
* These values will be 0.0 if no data are available.
* For V-OH3, a value is available only for the source listed in record 27,
* all other records have 0.0 in this field.
*Note (6): Distance to the star exciting the HII region, in kpc.
* This value is 0.0 if no distance is available.
*Note (7):
* This field contains a list of letters A through E separated by commas.
* The note inferred by each letter may be found in the "notes.dat" file.
C=============================================================================
C Loading file 'cosouth.dat' ! CO at Southern Hemisphere H II Regions
C Format for file interpretation
1 format(
+ A14,I3,I3,I3,1X,A1,I2,1X,F5.2,1X,F7.2,F7.2,2X,F5.1,1X,F4.1,2X,
+ F4.1,3X,F4.1,F7.1,F7.1,F7.1,F7.1,F10.2,1X,F5.1,A1,2X,F5.2,
+ F5.2,3X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'VII/50/cosouth.dat')
write(6,*) '....Loading file: cosouth.dat'
do i__=1,37
read(1,'(A138)')ar__
read(ar__,1)
+ Sname,RAh,RAm,RAs,DE_,DEd,DEm,GLON,GLAT,VLSR,EW,TAp,TAc,
+ V_H109,V_OH1,V_OH2,V_OH3,V_H2CO1,V_H2CO2,u_V_H2CO2,Dist,
+ e_Dist,Com
if(ar__(74:80) .EQ. '') V_H109 = rNULL__
if(ar__(81:87) .EQ. '') V_OH1 = rNULL__
if(ar__(88:94) .EQ. '') V_OH2 = rNULL__
if(ar__(95:101) .EQ. '') V_OH3 = rNULL__
if(ar__(102:111) .EQ. '') V_H2CO1 = rNULL__
if(ar__(113:117) .EQ. '') V_H2CO2 = 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(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ Sname,RAh,RAm,RAs,DE_,DEd,DEm,GLON,GLAT,VLSR,EW,TAp,TAc,
+ V_H109,V_OH1,V_OH2,V_OH3,V_H2CO1,V_H2CO2,u_V_H2CO2,Dist,
+ e_Dist,Com
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end