Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_AS/98/589 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 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 (nr__) ! IRAS denomination
character*6 Name (nr__) ! other name
real*4 sigmaH2O (nr__) ! (Jy) rms per channel of the spectrum taken on a
* given day (H2O result)
character*1 n_sigmaH2O (nr__) ! 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 (nr__) ! (km/s) []? LSR velocity of the strongest maser
* feature (H2O result)
character*1 n_V_H2O (nr__) ! [*] a '*' indicates value unknown
real*8 F_L_H2O (nr__) ! (Jy) []? peak flux density of the strongest
* feature (H2O result)
real*4 Vmin_H2O (nr__) ! (km/s) []? velocity range V_low, where emission is
* detected (H2O result)
character*1 n_Vmin_H2O (nr__) ! [*] a '*' indicates value unknown
real*4 Vmax_H2O (nr__) ! (km/s) []? velocity range V_high where emission is
* detected (H2O result)
character*1 n_Vmax_H2O (nr__) ! [*] a '*' indicates value unknown
real*8 I_F_L_H2O_ (nr__) ! (Jy.km/s) []? integrated flux density (H2O result)
real*4 L_H2O (nr__) ! (solLum) []? H2O luminosity
character*8 Obs (nr__) ! date of observation (as DD/MM/YY)
integer*4 r_Obs (nr__) ! references for observation:
* (1): this paper; (2): WW; (3): WBH;
* (4):KJB; (5): CWC; (6): PBC
real*4 n_CO (nr__) ! []? rms value if no CO has been detected
real*4 V_CO (nr__) ! (km/s) []? radial velocity (CO result)
real*4 d_CO (nr__) ! (kpc) []? kinematic distance (CO result)
real*4 LFIR (nr__) ! (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 (nr__1) ! (deg) Right Ascension 1950
real*8 DEdeg (nr__1) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*10 IRAS_1 (nr__1) ! IRAS denomination
character*1 n_IRAS (nr__1) ! a 'n' denotes no IRAS
integer*4 RA1h (nr__1) ! (h) []? right ascension 1950 for IRAS position
integer*4 RA1m (nr__1) ! (min) []? right ascension 1950 for IRAS position
real*4 RA1s (nr__1) ! (s) []? right ascension 1950 for IRAS position
character*1 DE1_ (nr__1) ! Declination sign
integer*4 DE1d (nr__1) ! (deg) []? declination 1950 for IRAS position
integer*4 DE1m (nr__1) ! (arcmin) []? declination 1950 for IRAS position
integer*4 DE1s (nr__1) ! (arcsec) []? declination 1950 for IRAS position
integer*4 RA2h (nr__1) ! (h) right ascension 1950 for H2O maser position
integer*4 RA2m (nr__1) ! (min) right ascension 1950 for H2O maser position
real*4 RA2s (nr__1) ! (s) right ascension 1950 for H2O maser position
character*1 DE2_ (nr__1) ! Declination sign
integer*4 DE2d (nr__1) ! (deg) declination 1950 for H2O maser position
integer*4 DE2m (nr__1) ! (arcmin) declination 1950 for H2O maser position
integer*4 DE2s (nr__1) ! (arcsec) declination 1950 for H2O maser position
character*4 Rem_ (nr__1) ! 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 (nr__2) ! IRAS denomination
character*21 Name_1 (nr__2) ! other name
character*1 H2O (nr__2) ! [-+] H2O detected (+) or not detected (-)
integer*4 r_H2O (nr__2) ! reference for H2O detection,
* (1): this paper, (2): WW, (3): WBH, (4): Be
character*3 v1665 (nr__2) ! OH result for 1665 MHz
* M: maser emission, A: absorption,
* -: not detected
character*3 v1667 (nr__2) ! 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 (nr__3) ! IRAS name
integer*4 Freq_ (nr__3) ! (MHz) frequency
real*4 FL1 (nr__3) ! (Jy) []? LCP peak flux density
real*4 VLSR1 (nr__3) ! (km/s) []? LCP LSR velocity of the strongest maser
* feature
real*4 Vmin1 (nr__3) ! (km/s) []? LCP velocity range V_low where emission is
* detected
real*4 Vmax1 (nr__3) ! (km/s) []? LCP velocity range V_high where emission
* is detected
real*4 I_FL1 (nr__3) ! (Jy.km/s) []? LCP integrated flux density
real*4 FL2 (nr__3) ! (Jy) []? RCP peak flux density
real*4 VLSR2 (nr__3) ! (km/s) []? RCP LSR velocity of the strongest maser
* feature
real*4 Vmin2 (nr__3) ! (km/s) []? RCP velocity range V_low where emission is
* detected
real*4 Vmax2 (nr__3) ! (km/s) []? RCP velocity range V_high where emission
* is detected
real*4 Fint2 (nr__3) ! (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 (nr__4) ! (deg) Right Ascension 1950
real*8 DEdeg_1 (nr__4) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*10 IRAS_4 (nr__4) ! IRAS name
integer*4 RAh (nr__4) ! (h) right ascension 1950 for OH maser position
integer*4 RAm (nr__4) ! (min) right ascension 1950 for OH maser position
real*4 RAs (nr__4) ! (s) right ascension 1950 for OH maser position
character*1 DE_ (nr__4) ! declination sign
integer*4 DEd (nr__4) ! (deg) declination 1950 for OH maser position
integer*4 DEm (nr__4) ! (arcmin) declination 1950 for OH maser position
integer*4 DEs (nr__4) ! (arcsec) declination 1950 for OH maser position
character*30 Rem__1 (nr__4) ! 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 (nr__5) ! IRAS or BE denomination
integer*4 Freq__1 (nr__5) ! (MHz) frequency
real*4 Tmb (nr__5) ! (K) brightness temperature
real*4 VLSR (nr__5) ! (km/s) LSR velocity
real*4 e_VLSR (nr__5) ! (km/s) uncertainty of LSR velocity
real*4 DeltaV (nr__5) ! (km/s) linewidth
real*4 e_DeltaV (nr__5) ! (km/s) uncertainty of linewidth
real*4 I_Tmb (nr__5) ! (K.km/s) brightness temperature integrated over
* velocity
real*4 e_I_Tmb (nr__5) ! (K.km/s) rms uncertainty on I(Tmb)
character*25 Rem__2 (nr__5) ! 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 (nr__6) ! IRAS denomination
real*4 GLON (nr__6) ! (deg) galactic longitude
real*4 GLAT (nr__6) ! (deg) galactic latitude
real*4 sigma (nr__6) ! (Jy) rms
real*4 DeltaV_1 (nr__6) ! (km/s) resolution of spectra
character*8 Date (nr__6) ! ("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 (nr__7) ! IRAS denomination
real*4 TR_ (nr__7) ! (K) []? temperature
character*1 n_TR_ (nr__7) ! a '*' denotes that rms = 0.2 when TR* is not
* detected
real*4 e_TR_ (nr__7) ! (K) []? rms uncertainty on TR*
real*4 VLSR_1 (nr__7) ! (km/s) []? LSR velocity
real*4 e_VLSR_1 (nr__7) ! (km/s) []? uncertainty of LSR velocity
real*4 DeltaV_2 (nr__7) ! (km/s) []? linewidth
real*4 e_DeltaV_1 (nr__7) ! (km/s) []? uncertainty of linewidth
character*17 Rem__3 (nr__7) ! remarks (1)
real*4 R (nr__7) ! (kpc) []? distance from galactic centre
real*4 d (nr__7) ! (kpc) []? distance from Sun
integer*4 z (nr__7) ! (pc) []? distance from galactic plane
real*4 LFIR_1 (nr__7) ! (solLum) []? FIR luminosity
character*5 Name_3 (nr__7) ! 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(i__),Name(i__),sigmaH2O(i__),n_sigmaH2O(i__),V_H2O(i__),
+ n_V_H2O(i__),F_L_H2O(i__),Vmin_H2O(i__),n_Vmin_H2O(i__),
+ Vmax_H2O(i__),n_Vmax_H2O(i__),I_F_L_H2O_(i__),L_H2O(i__),
+ Obs(i__),r_Obs(i__),n_CO(i__),V_CO(i__),d_CO(i__),LFIR(i__)
if(ar__(44:49) .EQ. '') V_H2O(i__) = rNULL__
if(ar__(57:63) .EQ. '') F_L_H2O(i__) = rNULL__
if(ar__(75:80) .EQ. '') Vmin_H2O(i__) = rNULL__
if(ar__(90:95) .EQ. '') Vmax_H2O(i__) = rNULL__
if(ar__(103:109) .EQ. '') I_F_L_H2O_(i__) = rNULL__
if(ar__(125:131) .EQ. '') L_H2O(i__) = rNULL__
if(ar__(158:160) .EQ. '') n_CO(i__) = rNULL__
if(ar__(167:172) .EQ. '') V_CO(i__) = rNULL__
if(ar__(183:186) .EQ. '') d_CO(i__) = rNULL__
if(ar__(193:197) .EQ. '') LFIR(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ IRAS(i__),Name(i__),sigmaH2O(i__),n_sigmaH2O(i__),V_H2O(i__),
+ n_V_H2O(i__),F_L_H2O(i__),Vmin_H2O(i__),n_Vmin_H2O(i__),
+ Vmax_H2O(i__),n_Vmax_H2O(i__),I_F_L_H2O_(i__),L_H2O(i__),
+ Obs(i__),r_Obs(i__),n_CO(i__),V_CO(i__),d_CO(i__),LFIR(i__)
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(i__),n_IRAS(i__),RA1h(i__),RA1m(i__),RA1s(i__),
+ DE1_(i__),DE1d(i__),DE1m(i__),DE1s(i__),RA2h(i__),RA2m(i__),
+ RA2s(i__),DE2_(i__),DE2d(i__),DE2m(i__),DE2s(i__),Rem_(i__)
if(ar__1(15:16) .EQ. '') RA1h(i__) = iNULL__
if(ar__1(18:19) .EQ. '') RA1m(i__) = iNULL__
if(ar__1(21:24) .EQ. '') RA1s(i__) = rNULL__
if(ar__1(27:28) .EQ. '') DE1d(i__) = iNULL__
if(ar__1(30:31) .EQ. '') DE1m(i__) = iNULL__
if(ar__1(33:34) .EQ. '') DE1s(i__) = iNULL__
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(RA1h(i__) .GT. -180) RAdeg(i__)=RA1h(i__)*15.
if(RA1m(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RA1m(i__)/4.
if(RA1s(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RA1s(i__)/240.
if(DE1d(i__) .GE. 0) DEdeg(i__)=DE1d(i__)
if(DE1m(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DE1m(i__)/60.
if(DE1s(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DE1s(i__)/3600.
if(DE1_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,2)
+ IRAS_1(i__),n_IRAS(i__),RA1h(i__),RA1m(i__),RA1s(i__),
+ DE1_(i__),DE1d(i__),DE1m(i__),DE1s(i__),RA2h(i__),RA2m(i__),
+ RA2s(i__),DE2_(i__),DE2d(i__),DE2m(i__),DE2s(i__),Rem_(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 '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(i__),Name_1(i__),H2O(i__),r_H2O(i__),v1665(i__),
+ v1667(i__)
c ..............Just test output...........
write(6,3)
+ IRAS_2(i__),Name_1(i__),H2O(i__),r_H2O(i__),v1665(i__),
+ v1667(i__)
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(i__),Freq_(i__),FL1(i__),VLSR1(i__),Vmin1(i__),
+ Vmax1(i__),I_FL1(i__),FL2(i__),VLSR2(i__),Vmin2(i__),
+ Vmax2(i__),Fint2(i__)
if(ar__3(24:27) .EQ. '') FL1(i__) = rNULL__
if(ar__3(32:36) .EQ. '') VLSR1(i__) = rNULL__
if(ar__3(41:45) .EQ. '') Vmin1(i__) = rNULL__
if(ar__3(50:54) .EQ. '') Vmax1(i__) = rNULL__
if(ar__3(62:65) .EQ. '') I_FL1(i__) = rNULL__
if(ar__3(68:71) .EQ. '') FL2(i__) = rNULL__
if(ar__3(76:80) .EQ. '') VLSR2(i__) = rNULL__
if(ar__3(85:89) .EQ. '') Vmin2(i__) = rNULL__
if(ar__3(94:98) .EQ. '') Vmax2(i__) = rNULL__
if(ar__3(106:109) .EQ. '') Fint2(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ IRAS_3(i__),Freq_(i__),FL1(i__),VLSR1(i__),Vmin1(i__),
+ Vmax1(i__),I_FL1(i__),FL2(i__),VLSR2(i__),Vmin2(i__),
+ Vmax2(i__),Fint2(i__)
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(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),Rem__1(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(i__) .GT. -180) RAdeg_1(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg_1(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,5)
+ IRAS_4(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),Rem__1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
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(i__),Freq__1(i__),Tmb(i__),VLSR(i__),e_VLSR(i__),
+ DeltaV(i__),e_DeltaV(i__),I_Tmb(i__),e_I_Tmb(i__),Rem__2(i__)
c ..............Just test output...........
write(6,6)
+ Name_2(i__),Freq__1(i__),Tmb(i__),VLSR(i__),e_VLSR(i__),
+ DeltaV(i__),e_DeltaV(i__),I_Tmb(i__),e_I_Tmb(i__),Rem__2(i__)
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(i__),GLON(i__),GLAT(i__),sigma(i__),DeltaV_1(i__),
+ Date(i__)
c ..............Just test output...........
write(6,7)
+ IRAS_5(i__),GLON(i__),GLAT(i__),sigma(i__),DeltaV_1(i__),
+ Date(i__)
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(i__),TR_(i__),n_TR_(i__),e_TR_(i__),VLSR_1(i__),
+ e_VLSR_1(i__),DeltaV_2(i__),e_DeltaV_1(i__),Rem__3(i__),
+ R(i__),d(i__),z(i__),LFIR_1(i__),Name_3(i__)
if(ar__7(13:16) .EQ. '') TR_(i__) = rNULL__
if(ar__7(19:20) .EQ. '') e_TR_(i__) = rNULL__
if(ar__7(23:27) .EQ. '') VLSR_1(i__) = rNULL__
if(ar__7(30:31) .EQ. '') e_VLSR_1(i__) = rNULL__
if(ar__7(34:37) .EQ. '') DeltaV_2(i__) = rNULL__
if(ar__7(40:41) .EQ. '') e_DeltaV_1(i__) = rNULL__
if(ar__7(62:66) .EQ. '') R(i__) = rNULL__
if(ar__7(68:72) .EQ. '') d(i__) = rNULL__
if(ar__7(74:77) .EQ. '') z(i__) = iNULL__
if(ar__7(79:83) .EQ. '') LFIR_1(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ IRAS_6(i__),TR_(i__),n_TR_(i__),e_TR_(i__),VLSR_1(i__),
+ e_VLSR_1(i__),DeltaV_2(i__),e_DeltaV_1(i__),Rem__3(i__),
+ R(i__),d(i__),z(i__),LFIR_1(i__),Name_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end