Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/230/1 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-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/ApJS/230/1 Herschel SPIRE/FTS 194-671um survey of GOALS LIRGs (Lu+, 2017)
*================================================================================
*A Herschel Space Observatory spectral line survey of local luminous infrared
*galaxies from 194 to 671 microns.
* Lu N., Zhao Y., Diaz-Santos T., Xu C.K., Gao Y., Armus L., Isaak K.G.,
* Mazzarella J.M., van der Werf P.P., Appleton P.N., Charmandaris V.,
* Evans A.S., Howell J., Iwasawa K., Leech J., Lord S., Petric A.O.,
* Privon G.C., Sanders D.B., Schulz B., Surace J.A.
* <Astrophys. J. Suppl. Ser., 230, 1-1 (2017)>
* =2017ApJS..230....1L (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! SPIRE/FTS observations
integer*4 nr__
parameter (nr__=127) ! Number of records
character*121 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*23 ID ! Object identifier (1)
character*7 OName ! Other identifier
character*1 n_ID ! Note on object (2)
integer*4 RAh ! (h) Hour of Right Ascension (J2000) (3)
integer*4 RAm ! (min) Minute of Right Ascension (J2000)
real*4 RAs ! (s) Second of Right Ascension (J2000)
character*1 DE_ ! Sign of the Declination (J2000) (3)
integer*4 DEd ! (deg) Degree of Declination (J2000) (3)
integer*4 DEm ! (arcmin) Arcminute of Declination (J2000)
integer*4 DEs ! (arcsec) Arcsecond of Declination (J2000)
real*4 Offset ! (arcsec) [0/5.7] Systematic pointing offset
real*4 logLIR ! ([Lsun]) [10.8/12.6]? Log adopted IR luminosity
character*1 f_logLIR ! [*] *: Target in a multiple system of 2 or
* more galaxies
real*4 CFIR ! [0.3/1.4]? FIR color; 60/100 micron band ratio
real*8 DL ! (Mpc) [15/271]? Luminosity distance
integer*4 HVel ! (km/s) [100/17880] Heliocentric velocity
character*10 ObsId ! Herschel observation identifier
integer*4 Exp ! (s) [1584/16996] Exposure time
character*17 FTS ! Fourier transform spectrometer (FTS) of SPIRE
* program identifier
*Note (1): That is spatially closest to the actual pointing of the
* SPIRE/FTS observation.
*Note (2): Flag as follows:
* a = A member of the subsample of powerful AGNs defined in Section 2.2.
* b = IRAS 05223+1908 consists of three objects. The SPIRE/FTS spectrum is
* likely dominated by the brightest northern object, which is a Galactic
* source at a heliocentric velocity of 100 km/s (see Appendix A for
* more details). While the SPIRE/FTS data of this target are included
* here, this target was not used in our statistical analysis.
* c = NGC 2146 has 3 independent observations, placed at the galaxy center
* and two offset positions, respectively. The total IRAS LIR is listed
* for each observation in the table.
* d = These 3 independent observations targeted the 3 surface brightness
* peaks of the interacting galaxy pair system NGC 3690. In the literature,
* UGC 06471 may also be referred to as NGC 3690 or Arp 299B, and
* UGC 06472 as IC 0694 or NGC 3690A or Arp 299A, The total IRAS LIR of
* the system is listed for each individual observation in the table.
* e = NGC 5010 was included in the SPIRE/FTS sample based on its LIR from
* Sanders+, 2003, Cat. J/AJ/126/1607, where its systemic velocity used is
* incorrect. While the SPIRE/FTS data of this target are still documented
* here, this target was not used in any statistical analysis in this paper.
*Note (3): The J2000 R.A. and decl. of the actual pointing of the SPIRE/FTS
* observation.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Fluxes of the CO, [CI] and [NII] lines
integer*4 nr__1
parameter (nr__1=127) ! Number of records
character*860 ar__1 ! Full-size record
character*23 ID_1 ! Object identifier
character*1 f_ID ! Flag on ID (4)
real*4 FCO4_3 ! (10-17W/m2) [-4.4/25.6] CO (4-3) line flux (5)
real*4 FCO5_4 ! (10-17W/m2) [-2.4/23.2] CO (5-4) line flux (5)
real*4 FCO6_5 ! (10-17W/m2) [-1.1/22.1] CO (6-5) line flux (5)
real*4 FCO7_6 ! (10-17W/m2) [-0.9/20.2] CO (7-6) line flux (5)
real*4 FCO8_7 ! (10-17W/m2) [-1.8/20.2] CO (8-7) line flux (5)
real*4 FCO9_8 ! (10-17W/m2) [-2/17.2] CO (9-8) line flux (5)
real*4 FCO10_9 ! (10-17W/m2) [-1.9/153.2] CO (10-9) line flux (5)
real*4 FCO11_10 ! (10-17W/m2) [-2.1/12.7] CO (11-10) line flux (5)
real*4 FCO12_11 ! (10-17W/m2) [-1.6/11] CO (12-11) line flux (5)
real*4 FCO13_12 ! (10-17W/m2) [-4.1/9.4] CO (13-12) line flux (5)
real*4 FCI609 ! (10-17W/m2) [-3.5/16.4] The [CI] 609um line flux (5)
real*4 FCI370 ! (10-17W/m2) [-0.7/29.6] The [CI] 370um line flux (5)
real*4 FNII ! (10-17W/m2) [-1.2/102.5] The [NII] 205um line flux
real*4 e_FCO4_3 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CO4-3 (6)
real*4 e_FCO5_4 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CO5-4 (6)
real*4 e_FCO6_5 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CO6-5 (6)
real*4 e_FCO7_6 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CO7-6 (6)
real*4 e_FCO8_7 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CO8-7 (6)
real*4 e_FCO9_8 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CO9-8 (6)
real*8 e_FCO10_9 ! (10-17W/m2) [0.06/6.756e+06]?=-99.99 Uncertainty
* in F-CO10-9 (6)
real*4 e_FCO11_10 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CO11-10 (6)
real*4 e_FCO12_11 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CO12-11 (6)
real*4 e_FCO13_12 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CO13-12 (6)
real*4 e_FCI609 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CI609 (6)
real*4 e_FCI370 ! (10-17W/m2) ?=-99.99 Uncertainty in F-CI370 (6)
real*4 e_FNII ! (10-17W/m2) ?=-99.99 Uncertainty in F-NII (6)
real*4 CFCO4_3 ! (GHz) ?=-99.99 CO (4-3) central frequency (7)
real*4 CFCO5_4 ! (GHz) ?=-99.99 CO (5-4) central frequency (7)
real*4 CFCO6_5 ! (GHz) ?=-99.99 CO (6-5) central frequency (7)
real*4 CFCO7_6 ! (GHz) ?=-99.99 CO (7-6) central frequency (7)
real*4 CFCO8_7 ! (GHz) ?=-99.99 CO (8-7) central frequency (7)
real*8 CFCO9_8 ! (GHz) ?=-99.99 CO (9-8) central frequency (7)
real*8 CFCO10_9 ! (GHz) ?=-99.99 CO (10-9) central frequency (7)
real*8 CFCO11_10 ! (GHz) ?=-99.99 CO (11-10) central frequency (7)
real*8 CFCO12_11 ! (GHz) ?=-99.99 CO (12-11) central frequency (7)
real*8 CFCO13_12 ! (GHz) ?=-99.99 CO (13-12) central frequency (7)
real*4 CFCI609 ! (GHz) ?=-99.99 The [CI] 609 micron line
* central frequency (7)
real*4 CFCI370 ! (GHz) ?=-99.99 The [CI] 370 micron line
* central frequency (7)
real*8 CFNII ! (GHz) ?=-99.99 The [NII] 205 micron line
* central frequency (7)
real*4 e_CFCO4_3 ! (GHz) ?=-99.99 Uncertainty in CF-CO4-3
real*4 e_CFCO5_4 ! (GHz) ?=-99.99 Uncertainty in CF-CO5-4
real*4 e_CFCO6_5 ! (GHz) ?=-99.99 Uncertainty in CF-CO6-5
real*4 e_CFCO7_6 ! (GHz) ?=-99.99 Uncertainty in CF-CO7-6
real*4 e_CFCO8_7 ! (GHz) ?=-99.99 Uncertainty in CF-CO8-7
real*4 e_CFCO9_8 ! (GHz) ?=-99.99 Uncertainty in CF-CO9-8
real*4 e_CFCO10_9 ! (GHz) ?=-99.99 Uncertainty in CF-CO10-9
real*4 e_CFCO11_10 ! (GHz) ?=-99.99 Uncertainty in CF-CO11-10
real*4 e_CFCO12_11 ! (GHz) ?=-99.99 Uncertainty in CF-CO12-11
real*4 e_CFCO13_12 ! (GHz) ?=-99.99 Uncertainty in CF-CO13-12
real*4 e_CFCI609 ! (GHz) ?=-99.99 Uncertainty in CF-CI609
real*4 e_CFCI370 ! (GHz) ?=-99.99 Uncertainty in CF-CI370
real*4 e_CFNII ! (GHz) ?=-99.99 Uncertainty in CF-NII
real*4 DCO4_3 ! (GHz) ?=-99.99 CO (4-3) line difference (8)
real*4 DCO5_4 ! (GHz) ?=-99.99 CO (5-4) line difference (8)
real*4 DCO6_5 ! (GHz) ?=-99.99 CO (6-5) line difference (8)
real*4 DCO7_6 ! (GHz) ?=-99.99 CO (7-6) line difference (8)
real*4 DCO8_7 ! (GHz) ?=-99.99 CO (8-7) line difference (8)
real*4 DCO9_8 ! (GHz) ?=-99.99 CO (9-8) line difference (8)
real*4 DCO10_9 ! (GHz) ?=-99.99 CO (10-9) line difference (8)
real*4 DCO11_10 ! (GHz) ?=-99.99 CO (11-10) line difference (8)
real*4 DCO12_11 ! (GHz) ?=-99.99 CO (12-11) line difference (8)
real*4 DCO13_12 ! (GHz) ?=-99.99 CO (13-12) line difference (8)
real*4 DCI609 ! (GHz) ?=-99.99 The [CI] 609um
* line difference (8)
real*4 DCI370 ! (GHz) ?=-99.99 The [CI] 370um
* line difference (8)
real*4 DNII ! (GHz) ?=-99.99 The [NII] 205um
* line difference (8)
real*4 PCO4_3 ! (Jy) ?=-99.99 Peak CO (4-3) line flux density
real*4 PCO5_4 ! (Jy) ?=-99.99 Peak CO (5-4) line flux density
real*4 PCO6_5 ! (Jy) ?=-99.99 Peak CO (6-5) line flux density
real*4 PCO7_6 ! (Jy) ?=-99.99 Peak CO (7-6) line flux density
real*4 PCO8_7 ! (Jy) ?=-99.99 Peak CO (8-7) line flux density
real*4 PCO9_8 ! (Jy) ?=-99.99 Peak CO (9-8) line flux density
real*4 PCO10_9 ! (Jy) ?=-99.99 Peak CO (10-9) line flux density
real*4 PCO11_10 ! (Jy) ?=-99.99 Peak CO (11-10) line flux density
real*4 PCO12_11 ! (Jy) ?=-99.99 Peak CO (12-11) line flux density
real*4 PCO13_12 ! (Jy) ?=-99.99 Peak CO (13-12) line flux density
real*4 PCI609 ! (Jy) ?=-99.99 Peak [CI] 609um line flux density
real*4 PCI370 ! (Jy) ?=-99.99 Peak [CI] 370um line flux density
real*4 PNII ! (Jy) ?=-99.99 Peak [NII] 205um line
* flux density
real*4 S_NCO4_3 ! ?=-99.99 Signal-to-Noise for
* CO (4-3) line (9)
real*4 S_NCO5_4 ! ?=-99.99 Signal-to-Noise for
* CO (5-4) line (9)
real*4 S_NCO6_5 ! ?=-99.99 Signal-to-Noise for
* CO (6-5) line (9)
real*4 S_NCO7_6 ! ?=-99.99 Signal-to-Noise for
* CO (7-6) line (9)
real*4 S_NCO8_7 ! ?=-99.99 Signal-to-Noise for
* CO (8-7) line (9)
real*4 S_NCO9_8 ! ?=-99.99 Signal-to-Noise for
* CO (9-8) line (9)
real*4 S_NCO10_9 ! ?=-99.99 Signal-to-Noise for
* CO (10-9) line (9)
real*4 S_NCO11_10 ! ?=-99.99 Signal-to-Noise for
* CO (11-10) line (9)
real*4 S_NCO12_11 ! ?=-99.99 Signal-to-Noise for
* CO (12-11) line (9)
real*4 S_NCO13_12 ! ?=-99.99 Signal-to-Noise for
* CO (13-12) line (9)
real*4 S_NCI609 ! ?=-99.99 Signal-to-Noise for
* [CI] 609um line (9)
real*4 S_NCI370 ! ?=-99.99 Signal-to-Noise for
* [CI] 370um line (9)
real*4 S_NNII ! ?=-99.99 Signal-to-Noise for
* [NII] 205um line (9)
real*4 q_CO4_3 ! ?=-99.99 Quality flag for
* CO (4-3) line (1=robust detection) (G1)
real*4 q_CO5_4 ! ?=-99.99 Quality flag for
* CO (5-4) line (1=robust detection) (G1)
real*4 q_CO6_5 ! ?=-99.99 Quality flag for
* CO (6-5) line (1=robust detection) (G1)
real*4 q_CO7_6 ! ?=-99.99 Quality flag for
* CO (7-6) line (1=robust detection) (G1)
real*4 q_CO8_7 ! ?=-99.99 Quality flag for
* CO (8-7) line (1=robust detection) (G1)
real*4 q_CO9_8 ! ?=-99.99 Quality flag for
* CO (9-8) line (1=robust detection) (G1)
real*4 q_CO10_9 ! ?=-99.99 Quality flag for
* CO (10-9) line (1=robust detection) (G1)
real*4 q_CO11_10 ! ?=-99.99 Quality flag for
* CO (11-10) line (1=robust detection) (G1)
real*4 q_CO12_11 ! ?=-99.99 Quality flag for
* CO (12-11) line (1=robust detection) (G1)
real*4 q_CO13_12 ! ?=-99.99 Quality flag for
* CO (13-12) line (1=robust detection) (G1)
real*4 q_CI609 ! ?=-99.99 Quality flag for
* [CI] 609um line (1=robust detection) (G1)
real*4 q_CI370 ! ?=-99.99 Quality flag for
* [CI] 370um line (1=robust detection) (G1)
real*4 q_NII ! ?=-99.99 Quality flag for
* [NII] 205um line (1=robust detection) (G1)
real*4 WCO4_3 ! (km/s) [-99.99] CO (4-3) line FWHM (10)
real*4 WCO5_4 ! (km/s) [-99.99] CO (5-4) line FWHM (10)
real*4 WCO6_5 ! (km/s) [-99.99] CO (6-5) line FWHM (10)
real*4 WCO7_6 ! (km/s) [-99.99] CO (7-6) line FWHM (10)
real*4 WCO8_7 ! (km/s) [-99.99] CO (8-7) line FWHM (10)
real*4 WCO9_8 ! (km/s) ?=-99.99 CO (9-8) line FWHM (10)
real*4 WCO10_9 ! (km/s) ?=-99.99 CO (10-9) line FWHM (10)
real*4 WCO11_10 ! (km/s) ?=-99.99 CO (11-10) line FWHM (10)
real*4 WCO12_11 ! (km/s) ?=-99.99 CO (12-11) line FWHM (10)
real*4 WCO13_12 ! (km/s) ?=-99.99 CO (13-12) line FWHM (10)
real*4 WCI609 ! (km/s) [-99.99] The [CI] 609um line FWHM (10)
real*4 WCI370 ! (km/s) [-99.99] The [CI] 370um line FWHM (10)
real*4 WNII ! (km/s) ?=-99.99 The [NII] 205um line FWHM (10)
real*4 F70_35 ! ? The 70 micron fractional flux
* with 35" beam (11)
real*4 F70_30 ! ? The 70 micron fractional flux
* with 30" beam (12)
real*4 F70_17 ! ? The 70 micron fractional flux
* with 17" beam (13)
*Note (4): Flag as follows:
* a = The spectrum of IRAS 05223+1908 shows strong CO lines that has an
* inferred heliocentric velocity of 100 km/s. This velocity was used
* for other line detections.
* b = These f_70micron_({theta}) values are all with respect to the
* total 70 micron flux of the galaxy NGC 2146.
* c = These f_70micron_({theta}) values are all with respect to the
* integrated 70 micron flux of the whole NGC 3690 system.
*Note (5): A positive number indicates that the line is detected, whilst a
* negative number represents a non-detection with its absolute value
* being the 3{sigma} upper limit, where {sigma} was set using a sinc
* line profile together with the local r.m.s. noise around the expected
* line frequency. For CO (4-3), a flux value of zero means that the
* line is redshifted out of the low frequency end of the SLW coverage.
* This is the case for a total of 17 targets with V_h_{<=}9558km/s.
*Note (6): This is the uncertainty from the line-fitting procedure.
*Note (7): The observed central frequency of the line in the
* Local Standard of Rest.
*Note (8): The difference between the observed line central frequency and the
* expected line frequency based on the heliocentric velocity of
* the target in Table 1.
*Note (9): The S/N ratio of the peak line flux density to the local
* r.m.s. noise {sigma}^t^_local_.
*Note (10): The FWHM of the Gaussian component when a sinc-Gaussian profile
* was used for the line fitting.
*Note (11): Appropriate for the SPIRE/FTS beam sizes of the CO (5-4),
* CO (7-6) or [CI] 370 micron lines.
*Note (12): Appropriate for the SPIRE/FTS beam size of the CO (6-5) line.
* Local Standard of Rest.
*Note (13): Appropriate for the SPIRE/FTS beam size of the [NII] line and the
* higher-J CO lines covered by the SSW spectral segment,
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Additional detected lines
integer*4 nr__2
parameter (nr__2=513) ! Number of records
character*105 ar__2 ! Full-size record
character*23 ID_2 ! Object identifier
character*32 Line ! Line description
real*4 Flux ! (10-17W/m2) [-11.5/14.3] Line flux (1)
character*1 f_Flux ! Flag on Flux (2)
real*4 e_Flux ! (10-17W/m2) [0.04/1.2] Uncertainty in Flux
real*8 nu ! (GHz) [515.9/1491.9] Observed line frequency
real*4 e_nu ! (GHz) [0.02/0.8] Uncertainty in nu
real*4 nudiff ! (GHz) [-2.2/2.1] Observed minus expected
* frequencies (3)
real*4 fnu ! (Jy) [-9.7/12.1] Peak line flux density (4)
real*4 S_N ! [2.9/34.5] Signal-to-Noise for
* peak line flux density
integer*4 Q ! [1/4] Quality flag (G1)
*Note (1): A negative value here indicates an absorption. Note that the flux
* derivation assumed a point-source case. See Section 4.2 for a
* prescription of flux aperture correction in case the target is
* moderately extended with respect to the SPIRE/FTS beam.
*Note (2): The flux values of the OH^+^ (1_22_-0_11_) and H_2_O(2_02_-1_11_)
* lines were:
* L = Measured in the SLW detector array.
* S = Measured in the SSW detector array.
*Note (3): This equals the observed frequency minus the expected line frequency
* based on the redshift of the target.
*Note (4): A negative value here also indicates an absorption.
C=============================================================================
C Loading file 'table1.dat' ! SPIRE/FTS observations
C Format for file interpretation
1 format(
+ A23,1X,A7,1X,A1,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,
+ F3.1,1X,F5.2,1X,A1,1X,F5.3,1X,F7.3,1X,I5,1X,A10,1X,I5,1X,A17)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,127
read(1,'(A121)')ar__
read(ar__,1)
+ ID,OName,n_ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Offset,logLIR,
+ f_logLIR,CFIR,DL,HVel,ObsId,Exp,FTS
if(ar__(60:64) .EQ. '') logLIR = rNULL__
if(ar__(68:72) .EQ. '') CFIR = rNULL__
if(ar__(74:80) .EQ. '') DL = 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(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ ID,OName,n_ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Offset,logLIR,
+ f_logLIR,CFIR,DL,HVel,ObsId,Exp,FTS
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Fluxes of the CO, [CI] and [NII] lines
C Format for file interpretation
2 format(
+ A23,3X,A1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F10.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F6.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,
+ F6.2,1X,F6.2,1X,F7.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,127
read(1,'(A860)')ar__1
read(ar__1,2)
+ ID_1,f_ID,FCO4_3,FCO5_4,FCO6_5,FCO7_6,FCO8_7,FCO9_8,FCO10_9,
+ FCO11_10,FCO12_11,FCO13_12,FCI609,FCI370,FNII,e_FCO4_3,
+ e_FCO5_4,e_FCO6_5,e_FCO7_6,e_FCO8_7,e_FCO9_8,e_FCO10_9,
+ e_FCO11_10,e_FCO12_11,e_FCO13_12,e_FCI609,e_FCI370,e_FNII,
+ CFCO4_3,CFCO5_4,CFCO6_5,CFCO7_6,CFCO8_7,CFCO9_8,CFCO10_9,
+ CFCO11_10,CFCO12_11,CFCO13_12,CFCI609,CFCI370,CFNII,e_CFCO4_3,
+ e_CFCO5_4,e_CFCO6_5,e_CFCO7_6,e_CFCO8_7,e_CFCO9_8,e_CFCO10_9,
+ e_CFCO11_10,e_CFCO12_11,e_CFCO13_12,e_CFCI609,e_CFCI370,
+ e_CFNII,DCO4_3,DCO5_4,DCO6_5,DCO7_6,DCO8_7,DCO9_8,DCO10_9,
+ DCO11_10,DCO12_11,DCO13_12,DCI609,DCI370,DNII,PCO4_3,PCO5_4,
+ PCO6_5,PCO7_6,PCO8_7,PCO9_8,PCO10_9,PCO11_10,PCO12_11,
+ PCO13_12,PCI609,PCI370,PNII,S_NCO4_3,S_NCO5_4,S_NCO6_5,
+ S_NCO7_6,S_NCO8_7,S_NCO9_8,S_NCO10_9,S_NCO11_10,S_NCO12_11,
+ S_NCO13_12,S_NCI609,S_NCI370,S_NNII,q_CO4_3,q_CO5_4,q_CO6_5,
+ q_CO7_6,q_CO8_7,q_CO9_8,q_CO10_9,q_CO11_10,q_CO12_11,
+ q_CO13_12,q_CI609,q_CI370,q_NII,WCO4_3,WCO5_4,WCO6_5,WCO7_6,
+ WCO8_7,WCO9_8,WCO10_9,WCO11_10,WCO12_11,WCO13_12,WCI609,
+ WCI370,WNII,F70_35,F70_30,F70_17
if(ar__1(847:850) .EQ. '') F70_35 = rNULL__
if(ar__1(852:855) .EQ. '') F70_30 = rNULL__
if(ar__1(857:860) .EQ. '') F70_17 = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1,f_ID,FCO4_3,FCO5_4,FCO6_5,FCO7_6,FCO8_7,FCO9_8,FCO10_9,
+ FCO11_10,FCO12_11,FCO13_12,FCI609,FCI370,FNII,e_FCO4_3,
+ e_FCO5_4,e_FCO6_5,e_FCO7_6,e_FCO8_7,e_FCO9_8,e_FCO10_9,
+ e_FCO11_10,e_FCO12_11,e_FCO13_12,e_FCI609,e_FCI370,e_FNII,
+ CFCO4_3,CFCO5_4,CFCO6_5,CFCO7_6,CFCO8_7,CFCO9_8,CFCO10_9,
+ CFCO11_10,CFCO12_11,CFCO13_12,CFCI609,CFCI370,CFNII,e_CFCO4_3,
+ e_CFCO5_4,e_CFCO6_5,e_CFCO7_6,e_CFCO8_7,e_CFCO9_8,e_CFCO10_9,
+ e_CFCO11_10,e_CFCO12_11,e_CFCO13_12,e_CFCI609,e_CFCI370,
+ e_CFNII,DCO4_3,DCO5_4,DCO6_5,DCO7_6,DCO8_7,DCO9_8,DCO10_9,
+ DCO11_10,DCO12_11,DCO13_12,DCI609,DCI370,DNII,PCO4_3,PCO5_4,
+ PCO6_5,PCO7_6,PCO8_7,PCO9_8,PCO10_9,PCO11_10,PCO12_11,
+ PCO13_12,PCI609,PCI370,PNII,S_NCO4_3,S_NCO5_4,S_NCO6_5,
+ S_NCO7_6,S_NCO8_7,S_NCO9_8,S_NCO10_9,S_NCO11_10,S_NCO12_11,
+ S_NCO13_12,S_NCI609,S_NCI370,S_NNII,q_CO4_3,q_CO5_4,q_CO6_5,
+ q_CO7_6,q_CO8_7,q_CO9_8,q_CO10_9,q_CO11_10,q_CO12_11,
+ q_CO13_12,q_CI609,q_CI370,q_NII,WCO4_3,WCO5_4,WCO6_5,WCO7_6,
+ WCO8_7,WCO9_8,WCO10_9,WCO11_10,WCO12_11,WCO13_12,WCI609,
+ WCI370,WNII,F70_35,F70_30,F70_17
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Additional detected lines
C Format for file interpretation
3 format(
+ A23,3X,A32,1X,F6.2,1X,A1,1X,F4.2,1X,F7.2,1X,F4.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,513
read(1,'(A105)')ar__2
read(ar__2,3)
+ ID_2,Line,Flux,f_Flux,e_Flux,nu,e_nu,nudiff,fnu,S_N,Q
c ..............Just test output...........
write(6,3)
+ ID_2,Line,Flux,f_Flux,e_Flux,nu,e_nu,nudiff,fnu,S_N,Q
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end