Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_AS/98/589 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 2026-Oct-08
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/A+AS/98/589 IRAS Sources behind the Solar circle (Wouterloot+ 1993)
*================================================================================
*IRAS Sources behind the Solar circle. III. Observations of H2O, OH, CH3OH and CO
* Wouterloot J.G.A., Brand J., Fiegle K.
* <Astron. Astrophys. Suppl. Ser. 98, 589 (1993)>
* =1993A&AS...98..589W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1' ! Line parameters towards IRAS sources observed
in H2O
integer*4 nr__
parameter (nr__=1759) ! Number of records
character*197 ar__ ! Full-size record
character*11 IRAS ! IRAS denomination
character*6 Name ! other name
real*4 sigmaH2O ! (Jy) rms per channel of the spectrum taken on a
* given day (H2O result)
character*1 n_sigmaH2O ! resolution of the spectra
* a: 0.64 km/s b: 0.32 km/s c: 0.16 km/s
* d: 0.08 km/s e: 0.04 km/s f: 0.02 km/s
real*4 V_H2O ! (km/s) []? LSR velocity of the strongest maser
* feature (H2O result)
character*1 n_V_H2O ! [*] a '*' indicates value unknown
real*8 F_L_H2O ! (Jy) []? peak flux density of the strongest
* feature (H2O result)
real*4 Vmin_H2O ! (km/s) []? velocity range V_low, where emission is
* detected (H2O result)
character*1 n_Vmin_H2O ! [*] a '*' indicates value unknown
real*4 Vmax_H2O ! (km/s) []? velocity range V_high where emission is
* detected (H2O result)
character*1 n_Vmax_H2O ! [*] a '*' indicates value unknown
real*8 I_F_L_H2O_ ! (Jy.km/s) []? integrated flux density (H2O result)
real*4 L_H2O ! (solLum) []? H2O luminosity
character*8 Obs ! date of observation (as DD/MM/YY)
integer*4 r_Obs ! references for observation:
* (1): this paper; (2): WW; (3): WBH;
* (4):KJB; (5): CWC; (6): PBC
real*4 n_CO ! []? rms value if no CO has been detected
real*4 V_CO ! (km/s) []? radial velocity (CO result)
real*4 d_CO ! (kpc) []? kinematic distance (CO result)
real*4 LFIR ! (solLum) []? total FIR luminosity
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2' ! Positions of H2O masers
integer*4 nr__1
parameter (nr__1=32) ! Number of records
character*62 ar__1 ! Full-size record
C 1950 position composed of: RA1h RA1m RA1s DE1- DE1d DE1m DE1s
real*8 RAdeg ! (deg) Right Ascension 1950
real*8 DEdeg ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*10 IRAS_1 ! IRAS denomination
character*1 n_IRAS ! a 'n' denotes no IRAS
integer*4 RA1h ! (h) []? right ascension 1950 for IRAS position
integer*4 RA1m ! (min) []? right ascension 1950 for IRAS position
real*4 RA1s ! (s) []? right ascension 1950 for IRAS position
character*1 DE1_ ! Declination sign
integer*4 DE1d ! (deg) []? declination 1950 for IRAS position
integer*4 DE1m ! (arcmin) []? declination 1950 for IRAS position
integer*4 DE1s ! (arcsec) []? declination 1950 for IRAS position
integer*4 RA2h ! (h) right ascension 1950 for H2O maser position
integer*4 RA2m ! (min) right ascension 1950 for H2O maser position
real*4 RA2s ! (s) right ascension 1950 for H2O maser position
character*1 DE2_ ! Declination sign
integer*4 DE2d ! (deg) declination 1950 for H2O maser position
integer*4 DE2m ! (arcmin) declination 1950 for H2O maser position
integer*4 DE2s ! (arcsec) declination 1950 for H2O maser position
character*4 Rem_ ! remarks
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3' ! Sources observed in OH
integer*4 nr__2
parameter (nr__2=302) ! Number of records
character*61 ar__2 ! Full-size record
character*10 IRAS_2 ! IRAS denomination
character*21 Name_1 ! other name
character*1 H2O ! [-+] H2O detected (+) or not detected (-)
integer*4 r_H2O ! reference for H2O detection,
* (1): this paper, (2): WW, (3): WBH, (4): Be
character*3 v1665 ! OH result for 1665 MHz
* M: maser emission, A: absorption,
* -: not detected
character*3 v1667 ! OH result for 1667 MHz
* M: maser emission, A: absorption,
* E: thermal emission, -: not detected
* blank: not observed
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4' ! Parameters of OH maser emission
integer*4 nr__3
parameter (nr__3=24) ! Number of records
character*109 ar__3 ! Full-size record
character*14 IRAS_3 ! IRAS name
integer*4 Freq_ ! (MHz) frequency
real*4 FL1 ! (Jy) []? LCP peak flux density
real*4 VLSR1 ! (km/s) []? LCP LSR velocity of the strongest maser
* feature
real*4 Vmin1 ! (km/s) []? LCP velocity range V_low where emission is
* detected
real*4 Vmax1 ! (km/s) []? LCP velocity range V_high where emission
* is detected
real*4 I_FL1 ! (Jy.km/s) []? LCP integrated flux density
real*4 FL2 ! (Jy) []? RCP peak flux density
real*4 VLSR2 ! (km/s) []? RCP LSR velocity of the strongest maser
* feature
real*4 Vmin2 ! (km/s) []? RCP velocity range V_low where emission is
* detected
real*4 Vmax2 ! (km/s) []? RCP velocity range V_high where emission
* is detected
real*4 Fint2 ! (Jy.km/s) []? RCP integrated flux density
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5' ! Positions of OH masers
integer*4 nr__4
parameter (nr__4=8) ! Number of records
character*63 ar__4 ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension 1950
real*8 DEdeg_1 ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*10 IRAS_4 ! IRAS name
integer*4 RAh ! (h) right ascension 1950 for OH maser position
integer*4 RAm ! (min) right ascension 1950 for OH maser position
real*4 RAs ! (s) right ascension 1950 for OH maser position
character*1 DE_ ! declination sign
integer*4 DEd ! (deg) declination 1950 for OH maser position
integer*4 DEm ! (arcmin) declination 1950 for OH maser position
integer*4 DEs ! (arcsec) declination 1950 for OH maser position
character*30 Rem__1 ! remarks
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6' ! Parameters of OH thermal absorption or emission
integer*4 nr__5
parameter (nr__5=69) ! Number of records
character*89 ar__5 ! Full-size record
character*15 Name_2 ! IRAS or BE denomination
integer*4 Freq__1 ! (MHz) frequency
real*4 Tmb ! (K) brightness temperature
real*4 VLSR ! (km/s) LSR velocity
real*4 e_VLSR ! (km/s) uncertainty of LSR velocity
real*4 DeltaV ! (km/s) linewidth
real*4 e_DeltaV ! (km/s) uncertainty of linewidth
real*4 I_Tmb ! (K.km/s) brightness temperature integrated over
* velocity
real*4 e_I_Tmb ! (K.km/s) rms uncertainty on I(Tmb)
character*25 Rem__2 ! remarks
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7' ! CH3OH observations
integer*4 nr__6
parameter (nr__6=22) ! Number of records
character*45 ar__6 ! Full-size record
C Position composed of: GLON GLAT
character*10 IRAS_5 ! IRAS denomination
real*4 GLON ! (deg) galactic longitude
real*4 GLAT ! (deg) galactic latitude
real*4 sigma ! (Jy) rms
real*4 DeltaV_1 ! (km/s) resolution of spectra
character*8 Date ! ("DD/MM/YY") observation date
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8' ! 12CO(2-1) observations
integer*4 nr__7
parameter (nr__7=34) ! Number of records
character*91 ar__7 ! Full-size record
character*10 IRAS_6 ! IRAS denomination
real*4 TR_ ! (K) []? temperature
character*1 n_TR_ ! a '*' denotes that rms = 0.2 when TR* is not
* detected
real*4 e_TR_ ! (K) []? rms uncertainty on TR*
real*4 VLSR_1 ! (km/s) []? LSR velocity
real*4 e_VLSR_1 ! (km/s) []? uncertainty of LSR velocity
real*4 DeltaV_2 ! (km/s) []? linewidth
real*4 e_DeltaV_1 ! (km/s) []? uncertainty of linewidth
character*17 Rem__3 ! remarks (1)
real*4 R ! (kpc) []? distance from galactic centre
real*4 d ! (kpc) []? distance from Sun
integer*4 z ! (pc) []? distance from galactic plane
real*4 LFIR_1 ! (solLum) []? FIR luminosity
character*5 Name_3 ! other name
*Note (1):
* as dist: for these sources we assumed a distance of 0.35 kpc
* nk dist: for these sources we used the near kinematic distances
C=============================================================================
C Loading file 'table1' ! Line parameters towards IRAS sources observed
* in H2O
C Format for file interpretation
1 format(
+ 1X,A11,3X,A6,8X,F3.1,A1,10X,F6.1,A1,6X,F7.1,11X,F6.1,A1,8X,
+ F6.1,A1,6X,F7.1,15X,E7.1,10X,A8,6X,I1,1X,F3.1,6X,F6.1,10X,
+ F4.1,6X,E5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1')
write(6,*) '....Loading file: table1'
do i__=1,1759
read(1,'(A197)')ar__
read(ar__,1)
+ IRAS,Name,sigmaH2O,n_sigmaH2O,V_H2O,n_V_H2O,F_L_H2O,Vmin_H2O,
+ n_Vmin_H2O,Vmax_H2O,n_Vmax_H2O,I_F_L_H2O_,L_H2O,Obs,r_Obs,
+ n_CO,V_CO,d_CO,LFIR
if(ar__(44:49) .EQ. '') V_H2O = rNULL__
if(ar__(57:63) .EQ. '') F_L_H2O = rNULL__
if(ar__(75:80) .EQ. '') Vmin_H2O = rNULL__
if(ar__(90:95) .EQ. '') Vmax_H2O = rNULL__
if(ar__(103:109) .EQ. '') I_F_L_H2O_ = rNULL__
if(ar__(125:131) .EQ. '') L_H2O = rNULL__
if(ar__(158:160) .EQ. '') n_CO = rNULL__
if(ar__(167:172) .EQ. '') V_CO = rNULL__
if(ar__(183:186) .EQ. '') d_CO = rNULL__
if(ar__(193:197) .EQ. '') LFIR = rNULL__
c ..............Just test output...........
write(6,1)
+ IRAS,Name,sigmaH2O,n_sigmaH2O,V_H2O,n_V_H2O,F_L_H2O,Vmin_H2O,
+ n_Vmin_H2O,Vmax_H2O,n_Vmax_H2O,I_F_L_H2O_,L_H2O,Obs,r_Obs,
+ n_CO,V_CO,d_CO,LFIR
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2' ! Positions of H2O masers
C Format for file interpretation
2 format(
+ 1X,A10,A1,2X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,2X,I2,1X,
+ I2,1X,F4.1,2X,A1,I2,1X,I2,1X,I2,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table2')
write(6,*) '....Loading file: table2'
do i__=1,32
read(1,'(A62)')ar__1
read(ar__1,2)
+ IRAS_1,n_IRAS,RA1h,RA1m,RA1s,DE1_,DE1d,DE1m,DE1s,RA2h,RA2m,
+ RA2s,DE2_,DE2d,DE2m,DE2s,Rem_
if(ar__1(15:16) .EQ. '') RA1h = iNULL__
if(ar__1(18:19) .EQ. '') RA1m = iNULL__
if(ar__1(21:24) .EQ. '') RA1s = rNULL__
if(ar__1(27:28) .EQ. '') DE1d = iNULL__
if(ar__1(30:31) .EQ. '') DE1m = iNULL__
if(ar__1(33:34) .EQ. '') DE1s = 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(RA1h .GT. -180) RAdeg=RA1h*15.
if(RA1m .GT. -180) RAdeg=RAdeg+RA1m/4.
if(RA1s .GT. -180) RAdeg=RAdeg+RA1s/240.
if(DE1d .GE. 0) DEdeg=DE1d
if(DE1m .GE. 0) DEdeg=DEdeg+DE1m/60.
if(DE1s .GE. 0) DEdeg=DEdeg+DE1s/3600.
if(DE1_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,2)
+ IRAS_1,n_IRAS,RA1h,RA1m,RA1s,DE1_,DE1d,DE1m,DE1s,RA2h,RA2m,
+ RA2s,DE2_,DE2d,DE2m,DE2s,Rem_
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3' ! Sources observed in OH
C Format for file interpretation
3 format(1X,A10,1X,A21,2X,A1,8X,I1,6X,A3,3X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table3')
write(6,*) '....Loading file: table3'
do i__=1,302
read(1,'(A61)')ar__2
read(ar__2,3)IRAS_2,Name_1,H2O,r_H2O,v1665,v1667
c ..............Just test output...........
write(6,3)IRAS_2,Name_1,H2O,r_H2O,v1665,v1667
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4' ! Parameters of OH maser emission
C Format for file interpretation
4 format(
+ 1X,A14,2X,I4,2X,F4.2,4X,F5.1,4X,F5.1,4X,F5.1,7X,F4.2,2X,F4.2,
+ 4X,F5.1,4X,F5.1,4X,F5.1,7X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4')
write(6,*) '....Loading file: table4'
do i__=1,24
read(1,'(A109)')ar__3
read(ar__3,4)
+ IRAS_3,Freq_,FL1,VLSR1,Vmin1,Vmax1,I_FL1,FL2,VLSR2,Vmin2,
+ Vmax2,Fint2
if(ar__3(24:27) .EQ. '') FL1 = rNULL__
if(ar__3(32:36) .EQ. '') VLSR1 = rNULL__
if(ar__3(41:45) .EQ. '') Vmin1 = rNULL__
if(ar__3(50:54) .EQ. '') Vmax1 = rNULL__
if(ar__3(62:65) .EQ. '') I_FL1 = rNULL__
if(ar__3(68:71) .EQ. '') FL2 = rNULL__
if(ar__3(76:80) .EQ. '') VLSR2 = rNULL__
if(ar__3(85:89) .EQ. '') Vmin2 = rNULL__
if(ar__3(94:98) .EQ. '') Vmax2 = rNULL__
if(ar__3(106:109) .EQ. '') Fint2 = rNULL__
c ..............Just test output...........
write(6,4)
+ IRAS_3,Freq_,FL1,VLSR1,Vmin1,Vmax1,I_FL1,FL2,VLSR2,Vmin2,
+ Vmax2,Fint2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5' ! Positions of OH masers
C Format for file interpretation
5 format(1X,A10,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A30)
C Effective file loading
open(unit=1,status='old',file=
+'table5')
write(6,*) '....Loading file: table5'
do i__=1,8
read(1,'(A63)')ar__4
read(ar__4,5)IRAS_4,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Rem__1
RAdeg_1 = rNULL__
DEdeg_1 = 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 .GT. -180) RAdeg_1=RAh*15.
if(RAm .GT. -180) RAdeg_1=RAdeg_1+RAm/4.
if(RAs .GT. -180) RAdeg_1=RAdeg_1+RAs/240.
if(DEd .GE. 0) DEdeg_1=DEd
if(DEm .GE. 0) DEdeg_1=DEdeg_1+DEm/60.
if(DEs .GE. 0) DEdeg_1=DEdeg_1+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,5)IRAS_4,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Rem__1
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6' ! Parameters of OH thermal absorption or emission
C Format for file interpretation
6 format(
+ 1X,A15,I4,1X,F6.2,1X,F5.1,2X,F3.1,2X,F4.1,2X,F3.1,2X,F6.1,2X,
+ F3.1,2X,A25)
C Effective file loading
open(unit=1,status='old',file=
+'table6')
write(6,*) '....Loading file: table6'
do i__=1,69
read(1,'(A89)')ar__5
read(ar__5,6)
+ Name_2,Freq__1,Tmb,VLSR,e_VLSR,DeltaV,e_DeltaV,I_Tmb,e_I_Tmb,
+ Rem__2
c ..............Just test output...........
write(6,6)
+ Name_2,Freq__1,Tmb,VLSR,e_VLSR,DeltaV,e_DeltaV,I_Tmb,e_I_Tmb,
+ Rem__2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7' ! CH3OH observations
C Format for file interpretation
7 format(1X,A10,1X,F6.2,2X,F5.2,3X,F3.1,1X,F4.2,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'table7')
write(6,*) '....Loading file: table7'
do i__=1,22
read(1,'(A45)')ar__6
read(ar__6,7)IRAS_5,GLON,GLAT,sigma,DeltaV_1,Date
c ..............Just test output...........
write(6,7)IRAS_5,GLON,GLAT,sigma,DeltaV_1,Date
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8' ! 12CO(2-1) observations
C Format for file interpretation
8 format(
+ 1X,A10,1X,F4.1,A1,1X,F2.1,2X,F5.1,2X,F2.1,2X,F4.1,2X,F2.1,2X,
+ A17,1X,F5.2,1X,F5.2,1X,I4,1X,E5.1,3X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'table8')
write(6,*) '....Loading file: table8'
do i__=1,34
read(1,'(A91)')ar__7
read(ar__7,8)
+ IRAS_6,TR_,n_TR_,e_TR_,VLSR_1,e_VLSR_1,DeltaV_2,e_DeltaV_1,
+ Rem__3,R,d,z,LFIR_1,Name_3
if(ar__7(13:16) .EQ. '') TR_ = rNULL__
if(ar__7(19:20) .EQ. '') e_TR_ = rNULL__
if(ar__7(23:27) .EQ. '') VLSR_1 = rNULL__
if(ar__7(30:31) .EQ. '') e_VLSR_1 = rNULL__
if(ar__7(34:37) .EQ. '') DeltaV_2 = rNULL__
if(ar__7(40:41) .EQ. '') e_DeltaV_1 = rNULL__
if(ar__7(62:66) .EQ. '') R = rNULL__
if(ar__7(68:72) .EQ. '') d = rNULL__
if(ar__7(74:77) .EQ. '') z = iNULL__
if(ar__7(79:83) .EQ. '') LFIR_1 = rNULL__
c ..............Just test output...........
write(6,8)
+ IRAS_6,TR_,n_TR_,e_TR_,VLSR_1,e_VLSR_1,DeltaV_2,e_DeltaV_1,
+ Rem__3,R,d,z,LFIR_1,Name_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end