Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/93/549 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-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. J/ApJS/93/549 Arecibo 1612 MHz survey supplement (Lewis 1994)
*================================================================================
*Supplement to the Arecibo 1612 MHz survey of color-selected IRAS sources
* Lewis B.M.
* <Astrophys. J. Suppl. Ser. 93, 549 (1994)>
* =1994ApJS...93..549L
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Profile parameters for 1612 MHz detections
integer*4 nr__
parameter (nr__=56) ! Number of records
character*133 ar__ ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension 1950
real*8 DEdeg ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*10 IRAS ! IRAS name
character*1 n_IRAS ! [*] *: part of color-selected Arecibo sample
character*1 Note ! [s] s: single to double peaked profiles
real*4 GLON ! (deg) Galactic longitude
real*4 GLAT ! (deg) Galactic latitude
integer*4 RAh ! (h) Right ascension (1950)
integer*4 RAm ! (min) Right ascension (1950)
real*4 RAs ! (s) Right ascension (1950)
character*1 DE_ ! Declination sign
integer*4 DEd ! (deg) Declination (1950)
integer*4 DEm ! (arcmin) Declination (1950)
integer*4 DEs ! (arcsec) Declination (1950)
real*8 S1612 ! (mJy.km/s) Total integrated 1612MHz flux
real*8 S1612_l ! (mJy.km/s) Integrated 1612MHz flux in the low-velocity
* half of the spectrum
real*8 S1612_h ! (mJy.km/s) ? Integrated 1612MHz flux in the high-velocity
* half of the spectrum
real*4 DV100 ! (km/s) ? Velocity width at 100% of the peak velocity
* {Delta}V_100_
character*1 u_DV100 ! [)] Uncertainty flag on DV100
real*4 DV50 ! (km/s) ? Velocity width at 50% of the peak velocity
* {Delta}V_50_
character*1 u_DV50 ! [)] Uncertainty flag on DV50
real*4 Vl ! (km/s) LSR velocity of the low-velocity peak
real*4 Vh ! (km/s) ? LSR velocity of the high-velocity peak
real*4 v_V_ ! (km/s) Systemic LSR velocity
real*4 v_25_12 ! (mag) (25-12){mu}m infrared colour (1)
real*4 v_60_25 ! (mag) (60-25){mu}m infrared colour (1)
integer*4 LRS ! ? IRAS low-resolution spectral-type
* (see Cat. III/197, Note 1)
character*10 Name ! Other name
*Note (1): (25-12)um = log{S(25*12*0.89/S(12)*25*1.09}
* (60-25)um = log{S(60)*25*0.82/S(25)*60*0.89}
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! IRAS sources without 1612 MHz emission
integer*4 nr__1
parameter (nr__1=210) ! Number of records
character*45 ar__1 ! Full-size record
C Position composed of: GLON GLAT
character*1 n_IRAS_1 ! [*] * : part of flux-limited Arecibo sample
character*10 IRAS_1 ! IRAS name
real*4 GLON_1 ! (deg) Galactic longitude
real*4 GLAT_1 ! (deg) Galactic latitude
real*8 v_12_25 ! (mag) (12-25){mu}m infrared color
real*8 v_60_25_1 ! (mag) (60-25){mu}m infrared color
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Main line parameters of 1612 MHz detections
integer*4 nr__2
parameter (nr__2=34) ! Number of records
character*69 ar__2 ! Full-size record
character*10 IRAS_2 ! IRAS name (1)
real*8 Sp1667 ! (mJy) 1667 MHz first peak flux
character*1 n_Sp1667 ! [s] Note on Sp1667 (2)
real*4 Vp1667 ! (km/s) ? 1667 MHz first peak velocity
real*4 Sp1667b ! (mJy) ? 1667 MHz second peak flux
real*4 Vp1667b ! (km/s) ? 1667 MHz second peak velocity
real*8 Sp1665 ! (mJy) ? 1665 MHz first peak flux
character*1 n_Sp1665 ! [sn)] Note on Sp1665 (2)
real*4 Vp1665 ! (km/s) ? 1665 MHz first peak velocity
character*1 u_Vp1665 ! [)] Uncertainty flag on Vp1665
real*4 Sp1665b ! (mJy) ? 1665 MHz second peak flux
character*1 n_Sp1665b ! [s] Note on Sp1665b (2)
real*4 Vp1665b ! (km/s) ? 1665 MHz second peak velocity
*Note (1): IRAS 19114+0002, 19142+1034, 19246+1637, and 20406+2953
* previously detected, the remaining 1612 MHz detections are
* mainline non-detections.
*Note (2): Flag as follows:
* s = the Sp value is {sigma} value
* n = not analyzed
* ) = uncertainty flag
C=============================================================================
C Loading file 'table1.dat' ! Profile parameters for 1612 MHz detections
C Format for file interpretation
1 format(
+ A10,A1,A1,1X,F5.1,1X,F5.1,2X,I2,I2,F4.1,2X,A1,I2,I2,I2,1X,
+ F8.1,1X,F7.1,1X,F8.1,2X,F4.1,A1,1X,F4.1,A1,1X,F5.1,2X,F5.1,1X,
+ F5.1,2X,F6.3,1X,F6.3,3X,I2,2X,A10)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,56
read(1,'(A133)')ar__
read(ar__,1)
+ IRAS,n_IRAS,Note,GLON,GLAT,RAh,RAm,RAs,DE_,DEd,DEm,DEs,S1612,
+ S1612_l,S1612_h,DV100,u_DV100,DV50,u_DV50,Vl,Vh,v_V_,v_25_12,
+ v_60_25,LRS,Name
if(ar__(62:69) .EQ. '') S1612_h = rNULL__
if(ar__(72:75) .EQ. '') DV100 = rNULL__
if(ar__(78:81) .EQ. '') DV50 = rNULL__
if(ar__(91:95) .EQ. '') Vh = rNULL__
if(ar__(120:121) .EQ. '') LRS = 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)
+ IRAS,n_IRAS,Note,GLON,GLAT,RAh,RAm,RAs,DE_,DEd,DEm,DEs,S1612,
+ S1612_l,S1612_h,DV100,u_DV100,DV50,u_DV50,Vl,Vh,v_V_,v_25_12,
+ v_60_25,LRS,Name
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! IRAS sources without 1612 MHz emission
C Format for file interpretation
2 format(A1,A10,3X,F5.1,1X,F5.1,3X,F7.4,3X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,210
read(1,'(A45)')ar__1
read(ar__1,2)n_IRAS_1,IRAS_1,GLON_1,GLAT_1,v_12_25,v_60_25_1
c ..............Just test output...........
write(6,2)n_IRAS_1,IRAS_1,GLON_1,GLAT_1,v_12_25,v_60_25_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Main line parameters of 1612 MHz detections
C Format for file interpretation
3 format(
+ A10,1X,F7.1,A1,1X,F6.2,F6.1,2X,F5.1,1X,F8.2,A1,F6.2,A1,F6.1,
+ A1,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,34
read(1,'(A69)')ar__2
read(ar__2,3)
+ IRAS_2,Sp1667,n_Sp1667,Vp1667,Sp1667b,Vp1667b,Sp1665,n_Sp1665,
+ Vp1665,u_Vp1665,Sp1665b,n_Sp1665b,Vp1665b
if(ar__2(21:26) .EQ. '') Vp1667 = rNULL__
if(ar__2(27:32) .EQ. '') Sp1667b = rNULL__
if(ar__2(35:39) .EQ. '') Vp1667b = rNULL__
if(ar__2(41:48) .EQ. '') Sp1665 = rNULL__
if(ar__2(50:55) .EQ. '') Vp1665 = rNULL__
if(ar__2(57:62) .EQ. '') Sp1665b = rNULL__
if(ar__2(65:69) .EQ. '') Vp1665b = rNULL__
c ..............Just test output...........
write(6,3)
+ IRAS_2,Sp1667,n_Sp1667,Vp1667,Sp1667b,Vp1667b,Sp1665,n_Sp1665,
+ Vp1665,u_Vp1665,Sp1665b,n_Sp1665b,Vp1665b
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end