Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/267/515 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-Sep-07
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+A/267/515 Very cold C-rich circumstellar envelopes (Omont+ 1993)
*================================================================================
*Characterization and proportion of very cold C-rich circumstellar envelopes
* Omont A., Loup C., Forveille T., te Lintel Hekkert P.,
* Habing H., Sivagnanam P.
* <Astron. Astrophys. 267, 515, (1993)>
* =1993A&A...267..515O
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Supergiants
integer*4 nr__
parameter (nr__=19) ! Number of records
character*58 ar__ ! Full-size record
character*3 No ! Identification number in tablea1 and tablea2
character*10 Name ! Usual name
character*1 Sp ! Spectral type
character*5 Region ! Region in the IRAS color-color diagram (G1)
character*3 LRS ! LRS spectral class
integer*4 FIRAS ! (solLum) IRAS total flux expressed in solar
* luminosity at 1kpc
character*1 n ! [+?] A + indicates the post-AGB nature of
* the star
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Sources with CO(2-1) or HCN(1-0) emission
integer*4 nr__1
parameter (nr__1=57) ! Number of records
character*120 ar__1 ! 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*3 No_1 ! Identification number in tablea1 and tablea2
character*10 Name_1 ! Usual name
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
character*1 Chem ! [COS] Chemical type (C, O or S)
character*1 n_Chem ! [*#] (1)
character*1 u_Chem ! [?] uncertainty flag on chemical type
character*5 Region_1 ! Region in the IRAS color-color diagram (G1)
character*9 Trans ! Observed transitions
integer*4 Vlsr ! (km/s) ? Star Local Standard of Rest velocity
real*4 Ve ! (km/s) ? Expansion velocity as the half width at
* zero power
real*4 TMB ! (K) ? Main-beam brightness temperature measured
* at the peak of the line
real*4 e_TMB ! (K) ? rms noise of T_MB
real*4 I ! (K.km/s) ? Integrated area
character*1 u_I ! [?] uncertainty flag on I
character*1 l_CO2 ! limit flag on ICO/IHCN
real*4 CO2 ! ? Intensity ratio CO(2-1)/HCN(1-0)
real*4 S60_CO ! ? Ratio 60{mu}m to CO(2-1) emissions
character*24 Rem ! Remarks and references (3)
*Note (1):
* * indicates that the determination of the chemical type comes from this
* paper,
* # indicates that the chemical type was confirmed in this paper
*Note (3):
* Obervation date, mapping, interstellar contamination,
* numbers in brackets indicate references of previous detections
* (1) Zuckerman & Dyck 1986ApJ...304..394Z
* (2) Zuckerman et al. 1986ApJ...304..401Z
* (3) Zuckerman & Dyck 1986ApJ...311..345Z
* (4) Arquilla et al. 1986MNRAS.220..125A
* (5) Likkel et al. 1987A&A...173L..11L
* (6) Lucas et al. 1988A&A...194..230L
* (7) Lindqvist et al 1988A&A...205L..15L
* (8) Knapp et al. 1989ApJ...336..822K
* (9) Zuckerman & Dyck 1989A&A...209..119Z
* (10) Woodsworth et al. 1990A&A...228..503W
* (11) Likkel et al. 1991A&A...246..153L
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Non-detected sources
integer*4 nr__2
parameter (nr__2=27) ! Number of records
character*120 ar__2 ! 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*3 No_2 ! Identification number in tablea1 and tablea2
character*10 Name_2 ! Usual name
integer*4 RAh_1 ! (h) ? Right ascension 1950
integer*4 RAm_1 ! (min) ? Right ascension 1950
real*4 RAs_1 ! (s) ? Right ascension 1950
character*1 DE__1 ! ? Declination sign
integer*4 DEd_1 ! (deg) ? Declination 1950
integer*4 DEm_1 ! (arcmin) ? Declination 1950
integer*4 DEs_1 ! (arcsec) ? Declination 1950
character*1 Chem_1 ! [COS] Chemical type (C, O or S)
character*1 n_Chem_1 ! [*#] (1)
character*1 u_Chem_1 ! [?] uncertainty flag on chemical type
character*5 Region_2 ! Region in the IRAS color-color diagram (G1)
character*9 Trans_1 ! Observed transitions
integer*4 Vlsr_1 ! (km/s) ? Star Local Standard of Rest velocity
real*4 Ve_1 ! (km/s) ? Expansion velocity as the half width at
* zero power
real*4 TMB_1 ! (K) ? Main-beam brightness temperature measured
* at the peak of the line
real*4 e_TMB_1 ! (K) ? rms noise of T_MB
real*4 I_1 ! (K.km/s) ? Integrated area
character*1 u_I_1 ! [?] uncertainty flag on I
character*1 l_CO2_1 ! limit flag on ICO/IHCN
real*4 CO2_1 ! ? Intensity ratio CO(2-1)/HCN(1-0)
real*4 S60_CO_1 ! ? Ratio 60{mu}m to CO(2-1) emissions
character*24 Rem_1 ! Remarks and references (3)
*Note (1):
* * indicates that the determination of the chemical type comes from this
* paper,
* # indicates that the chemical type was confirmed in this paper
*Note (3):
* Obervation date, mapping, interstellar contamination,
* numbers in brackets indicate references of previous detections
* (1) Zuckerman & Dyck 1986ApJ...304..394Z
* (2) Zuckerman et al. 1986ApJ...304..401Z
* (3) Zuckerman & Dyck 1986ApJ...311..345Z
* (4) Arquilla et al. 1986MNRAS.220..125A
* (5) Likkel et al. 1987A&A...173L..11L
* (6) Lucas et al. 1988A&A...194..230L
* (7) Lindqvist et al 1988A&A...205L..15L
* (8) Knapp et al. 1989ApJ...336..822K
* (9) Zuckerman & Dyck 1989A&A...209..119Z
* (10) Woodsworth et al. 1990A&A...228..503W
* (11) Likkel et al. 1991A&A...246..153L
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! CO/HCN line ratios
integer*4 nr__3
parameter (nr__3=28) ! Number of records
character*77 ar__3 ! Full-size record
character*3 No_3 ! Identification number in tablea1 & tablea2
character*10 Name_3 ! Usual name
character*5 Region_3 ! Region in the IRAS color-color diagram (G1)
character*6 LRS_1 ! LRS spectral class
character*1 Chem_2 ! [COS] Chemical type
character*1 u_Chem_2 ! [?] uncertainty flag on chemical type
character*1 l_CO1 ! [~] limit flag on CO(1-0)/HCN
real*4 CO1 ! ? CO(1-0)/HCN line ratio
real*4 CO2_2 ! ? CO(2-1)/HCN line ratio
character*5 Ref ! References (2)
character*15 Com ! Comments
*Note (2):
* (0) present paper (Table 2)
* (13) Nguyen-Q-Rieu et al. 1987A&A...180..117N
* (14) Sopka et al. 1989A&A...210...78S
* (15) Knapp & Morris 1985ApJ...292..640K
* (16) Nercessian et al. 1989A&A...210..225N
* (17) Bujarrabal et al. 1988A&A...204..242B
* (18) Morris et al. 1987ApJ...321..888M
* (19) Bachiller et al. 1988A&A...196L...5B
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Estimates of mass-loss rates
integer*4 nr__4
parameter (nr__4=24) ! Number of records
character*80 ar__4 ! Full-size record
character*3 No_4 ! Identification number in tablea1 & tablea2
character*10 Name_4 ! Usual name
character*5 Region_4 ! Region in the IRAS color-color diagram (G1)
character*6 LRS_2 ! LRS spectral class
character*1 Chem_3 ! [COS] Chemical type
character*1 u_Chem_3 ! [?] uncertainty flag on chemical type
real*4 BC ! Bolometric correction F_tot_/F(IRAS)
character*1 r_BC ! references for Ftot (2)
real*4 D ! (kpc) Distance
character*1 n_D ! note for D (3)
real*4 dM1_dt ! (10-5Msun/yr) ? mass-loss rate from CO(1-0) (4)
character*3 r_dM1_dt ! references for J=1-0 CO observations
real*4 dM2_dt ! (10-5Msun/yr) ? mass-loss rate from CO(2-1) (4)
character*1 f_dM2_dt ! [pq] problematic dMCO(2-1)/dt (3)
real*4 dMIR_dt ! ? mass-loss rates deduced from infrared
* observations following Jura & Kleinmann
* (1989ApJ...341..359J)
*Note (2): References as follows:
* a = van der Veen & Breukers (1989A&A...213..133V)
* from actual 1-60{mu}m photometry
* b = applying to S12 the bolometric corrections proposed by
* van der Veen & Breukers (1989A&A...213..133V)
* c = assuming the following values for Ftot/Firas:
* 2 for region IIIa1, 1.2 for region IIIb2,
* 1.0 for region V without visible star,
* 2.5 for region V where a visible star is reported
* d = Kwok et al. 1989.
*Note (3):
* 'p' inappropriate if the source is a supergiant
* 'q' extended CO(2-1) source, dM2/dt should be multiplied by ~2
*Note (4):
* mass-loss rate are deduced from Eqs (3) and (4) assuming the abundance
* of CO relative to H_2_ f(CO): f(CO)=10^-3^ for C-rich sources and
* f(CO)=5.10^-4^ for other sources
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Master list of the main sample
integer*4 nr__5
parameter (nr__5=228) ! Number of records
character*113 ar__5 ! Full-size record
character*1 n_No ! [*] A "*" preceding the number indicates that
* the source was observed in the millimeter
* lines in the present study
character*3 No_5 ! Classification number
character*10 IRAS ! IRAS-PSC (Cat. II/125) name
real*4 C21 ! 25/12{mu}m index, C21=log(12F_25_/25F_12_)
character*1 u_C21 ! [)] uncertainty flag on C21
real*4 C32 ! 60/25{mu}m index, C32=log(25F_60_/12F_25_)
character*1 u_C32 ! [)] uncertainty flag on C32
real*4 FIRAS_1 ! (solLum) ? IRAS total flux expressed in solar
* luminosity at 1kpc
integer*4 Var_ ! ? variability index from the IRAS-PSC
* catalogue (Cat. II/125)
character*7 Region_5 ! Region of the IRAS color-color diagram (G1)
character*1 OH ! [-+.] OH maser detection (2)
character*4 r_OH ! References for OH maser observations (3)
character*6 LRS_3 ! LRS spectral class (4)
character*1 u_LRS ! [?] uncertainty flag on LRS
character*1 Chem_4 ! [COS] Proposed chemical type (C, O or S)
character*1 u_Chem_4 ! [?] uncertainty flag on chemical type
character*1 CO ! CO detections (2)
character*1 HCN ! [-+.*] HCN detections (2)
character*1 H2O ! [-+.] H2O detections (+) or non-detections (-)
* from the catalogue of Comoretto et al.
* (1990A&AS...84..179C)
character*29 Rem_2 ! Remarks (5)
*Note (2): code has the following meaning:
* + = detected with circumstellar characteristics;
* - = undetected or interstellar characteristics;
* . = not searched
* * = new detection of the present paper
*Note (3): the symbols refer to:
* t = catalogue of OH/IR stars by te Lintel Hekkert et al.
* (1989A&AS...78..399T)
* D,B,P = observations at Dwingeloo, Bonn-Effelsberg, and Parkes,
* respectively, in te Lintel Hekkert et al. (1991A&AS...90..327T)
* n,N,N'= observations at Nancay:
* n = Sivagnanam et al. (1990A&A...233..112S)
* N = Le Squeren et al. (1992A&A...254..133L)
* N' = te Lintel Hekkert et al. (1991, in preparation)
* L = Likkel (1989ApJ...344..350L)
* a = observations at Arecibo
* W = catalogue of observations of stellar masers by Benson et al.
* (1990ApJS...74..911B)
*Note (4): LRS spectral class, from the IRAS-LRS atlas (1986, Cat. III/197)
* when available; otherwise, from the IRAS-LRS database, without
* determination of the subclass n (Appendix).
* * Most sources classified 5n would have appeared as 50 or 05 in the
* LRS atlas.
* * Spectra quoted as 4no are 4n spectra with 18 micrometer silicate
* emission (Fig.2, rows II and III).
* * "kvh" indicates the presence of the 21 micrometer emission
* (Kwok et al., 1989ApJ...345L..51K, Fig.5, row II).
* * "int" denotes intermediate spectra with a maximum in the middle
* of the LRS range (Fig.5, row I), probably characteristics of the
* very early stages of the detachment of PPN envelopes (Loup, 1991, thesis).
*Note (5): The notes include:
* * Most usual name(s):
* T = IRC from the Two Micron Sky Survey (Neugebauer and Leighton, 1969,
* see catalog II/2)
* GL = RAFGL catalogue (Price and Murdock, 1983; see catalog II/94)
* SAO = SAO catalogue (Cat. I/131)
* V = General Catalogue of Variables Stars (Kholopov 1985-1987, Cat. II/139)
* DO = Deaborn catalogue (Cat. II/68)
* * visible spectral types mostly from the IRAS PSC catalogue (Cat. II/125)
* * zeta44 autoclass from the auto-classification of LRS spectra by
* Cheeseman et al. (1989, NASA Ref. Publ. 1217) for objects of
* autoclass zeta44
* * Sources marked with a * present a special interest and deserve further
* peculiar studies. They include (objects xa belong to Table A2):
* a = double-shell envelopes of region VIb, Nos 29, 57, 74, 89, 124, 130,
* 141, 16a
* b = 21{mu}m sources (kvh spectra), Nos 33, 216, 4a, 17a, possibly 13a
* c = "int" spectra, Nos 85, 136, 159, 179 (Fig.5), 10a
* d = 4n OH emitters without 18 micrometer feature (Le Squeren et al.,
* 1992A&A...254..133L), Nos 67, 142;
* e = No. 8a with a very strong silicate absorption and no OH emission
* f = OH+HCN emitter, No.208
* g = The peculiar OH spectrum 17376-3021 (No. 68) (te Lintel Hekkert et
* al. 991A&AS...90..327T)
* h = Sources associated with visible stars: Be (Nos. 24b, 26, 30, 102,
* 109, 195) or others (Nos. 36, 38, 5a, 11a, 15a)
* i = Peculiar objects of region V:17441-2411 (GL 2385, No.72), which
* displays extended CO emission; 18095+2704 (No.95) with a very
* strong silicate absorption and a very strong 25 micrometer excess
* (Volg & Kwok 1988ApJ...331..435V)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Additional sources
integer*4 nr__6
parameter (nr__6=18) ! Number of records
character*113 ar__6 ! Full-size record
character*1 n_No_1 ! [*] A "*" preceding the number indicates that
* the source was observed in the millimeter
* lines in the present study
character*3 No_6 ! Classification number
character*10 IRAS_1 ! IRAS-PSC (Cat. II/125) name
real*4 C21_1 ! 25/12{mu}m index, C21=log(12F_25_/25F_12_)
character*1 u_C21_1 ! [)] uncertainty flag on C21
real*4 C32_1 ! 60/25{mu}m index, C32=log(25F_60_/12F_25_)
character*1 u_C32_1 ! [)] uncertainty flag on C32
real*4 FIRAS_2 ! (solLum) ? IRAS total flux expressed in solar
* luminosity at 1kpc
integer*4 Var__1 ! ? variability index from the IRAS-PSC
* catalogue (Cat. II/125)
character*7 Region_6 ! Region of the IRAS color-color diagram (G1)
character*1 OH_1 ! [-+.] OH maser detection (2)
character*4 r_OH_1 ! References for OH maser observations (3)
character*6 LRS_4 ! LRS spectral class (4)
character*1 u_LRS_1 ! [?] uncertainty flag on LRS
character*1 Chem_5 ! [COS] Proposed chemical type (C, O or S)
character*1 u_Chem_5 ! [?] uncertainty flag on chemical type
character*1 CO_1 ! CO detections (2)
character*1 HCN_1 ! [-+.*] HCN detections (2)
character*1 H2O_1 ! [-+.] H2O detections (+) or non-detections (-)
* from the catalogue of Comoretto et al.
* (1990A&AS...84..179C)
character*29 Rem_3 ! Remarks (5)
*Note (2): code has the following meaning:
* + = detected with circumstellar characteristics;
* - = undetected or interstellar characteristics;
* . = not searched
* * = new detection of the present paper
*Note (3): the symbols refer to:
* t = catalogue of OH/IR stars by te Lintel Hekkert et al.
* (1989A&AS...78..399T)
* D,B,P = observations at Dwingeloo, Bonn-Effelsberg, and Parkes,
* respectively, in te Lintel Hekkert et al. (1991A&AS...90..327T)
* n,N,N'= observations at Nancay:
* n = Sivagnanam et al. (1990A&A...233..112S)
* N = Le Squeren et al. (1992A&A...254..133L)
* N' = te Lintel Hekkert et al. (1991, in preparation)
* L = Likkel (1989ApJ...344..350L)
* a = observations at Arecibo
* W = catalogue of observations of stellar masers by Benson et al.
* (1990ApJS...74..911B)
*Note (4): LRS spectral class, from the IRAS-LRS atlas (1986, Cat. III/197)
* when available; otherwise, from the IRAS-LRS database, without
* determination of the subclass n (Appendix).
* * Most sources classified 5n would have appeared as 50 or 05 in the
* LRS atlas.
* * Spectra quoted as 4no are 4n spectra with 18 micrometer silicate
* emission (Fig.2, rows II and III).
* * "kvh" indicates the presence of the 21 micrometer emission
* (Kwok et al., 1989ApJ...345L..51K, Fig.5, row II).
* * "int" denotes intermediate spectra with a maximum in the middle
* of the LRS range (Fig.5, row I), probably characteristics of the
* very early stages of the detachment of PPN envelopes (Loup, 1991, thesis).
*Note (5): The notes include:
* * Most usual name(s):
* T = IRC from the Two Micron Sky Survey (Neugebauer and Leighton, 1969,
* see catalog II/2)
* GL = RAFGL catalogue (Price and Murdock, 1983; see catalog II/94)
* SAO = SAO catalogue (Cat. I/131)
* V = General Catalogue of Variables Stars (Kholopov 1985-1987, Cat. II/139)
* DO = Deaborn catalogue (Cat. II/68)
* * visible spectral types mostly from the IRAS PSC catalogue (Cat. II/125)
* * zeta44 autoclass from the auto-classification of LRS spectra by
* Cheeseman et al. (1989, NASA Ref. Publ. 1217) for objects of
* autoclass zeta44
* * Sources marked with a * present a special interest and deserve further
* peculiar studies. They include (objects xa belong to Table A2):
* a = double-shell envelopes of region VIb, Nos 29, 57, 74, 89, 124, 130,
* 141, 16a
* b = 21{mu}m sources (kvh spectra), Nos 33, 216, 4a, 17a, possibly 13a
* c = "int" spectra, Nos 85, 136, 159, 179 (Fig.5), 10a
* d = 4n OH emitters without 18 micrometer feature (Le Squeren et al.,
* 1992A&A...254..133L), Nos 67, 142;
* e = No. 8a with a very strong silicate absorption and no OH emission
* f = OH+HCN emitter, No.208
* g = The peculiar OH spectrum 17376-3021 (No. 68) (te Lintel Hekkert et
* al. 991A&AS...90..327T)
* h = Sources associated with visible stars: Be (Nos. 24b, 26, 30, 102,
* 109, 195) or others (Nos. 36, 38, 5a, 11a, 15a)
* i = Peculiar objects of region V:17441-2411 (GL 2385, No.72), which
* displays extended CO emission; 18095+2704 (No.95) with a very
* strong silicate absorption and a very strong 25 micrometer excess
* (Volg & Kwok 1988ApJ...331..435V)
C=============================================================================
C Loading file 'table1.dat' ! Supergiants
C Format for file interpretation
1 format(A3,3X,A10,3X,A1,5X,A5,3X,A3,4X,I5,12X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,19
read(1,'(A58)')ar__
read(ar__,1)No,Name,Sp,Region,LRS,FIRAS,n
c ..............Just test output...........
write(6,1)No,Name,Sp,Region,LRS,FIRAS,n
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Sources with CO(2-1) or HCN(1-0) emission
C Format for file interpretation
2 format(
+ A3,1X,A10,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A1,A1,
+ A1,1X,A5,1X,A9,1X,I3,1X,F4.1,1X,F5.3,1X,F5.3,1X,F5.1,A1,A1,
+ F4.1,1X,F5.2,2X,A24)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,57
read(1,'(A120)')ar__1
read(ar__1,2)
+ No_1,Name_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Chem,n_Chem,u_Chem,
+ Region_1,Trans,Vlsr,Ve,TMB,e_TMB,I,u_I,l_CO2,CO2,S60_CO,Rem
if(ar__1(16:17) .EQ. '') RAh = iNULL__
if(ar__1(19:20) .EQ. '') RAm = iNULL__
if(ar__1(22:25) .EQ. '') RAs = rNULL__
if(ar__1(28:29) .EQ. '') DEd = iNULL__
if(ar__1(31:32) .EQ. '') DEm = iNULL__
if(ar__1(34:35) .EQ. '') DEs = iNULL__
if(ar__1(57:59) .EQ. '') Vlsr = iNULL__
if(ar__1(61:64) .EQ. '') Ve = rNULL__
if(ar__1(66:70) .EQ. '') TMB = rNULL__
if(ar__1(72:76) .EQ. '') e_TMB = rNULL__
if(ar__1(78:82) .EQ. '') I = rNULL__
if(ar__1(85:88) .EQ. '') CO2 = rNULL__
if(ar__1(90:94) .EQ. '') S60_CO = 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,2)
+ No_1,Name_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Chem,n_Chem,u_Chem,
+ Region_1,Trans,Vlsr,Ve,TMB,e_TMB,I,u_I,l_CO2,CO2,S60_CO,Rem
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Non-detected sources
C Format for file interpretation
3 format(
+ A3,1X,A10,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A1,A1,
+ A1,1X,A5,1X,A9,1X,I3,1X,F4.1,1X,F5.3,1X,F5.3,1X,F5.1,A1,A1,
+ F4.1,1X,F5.2,2X,A24)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,27
read(1,'(A120)')ar__2
read(ar__2,3)
+ No_2,Name_2,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Chem_1,
+ n_Chem_1,u_Chem_1,Region_2,Trans_1,Vlsr_1,Ve_1,TMB_1,e_TMB_1,
+ I_1,u_I_1,l_CO2_1,CO2_1,S60_CO_1,Rem_1
if(ar__2(16:17) .EQ. '') RAh_1 = iNULL__
if(ar__2(19:20) .EQ. '') RAm_1 = iNULL__
if(ar__2(22:25) .EQ. '') RAs_1 = rNULL__
if(ar__2(28:29) .EQ. '') DEd_1 = iNULL__
if(ar__2(31:32) .EQ. '') DEm_1 = iNULL__
if(ar__2(34:35) .EQ. '') DEs_1 = iNULL__
if(ar__2(57:59) .EQ. '') Vlsr_1 = iNULL__
if(ar__2(61:64) .EQ. '') Ve_1 = rNULL__
if(ar__2(66:70) .EQ. '') TMB_1 = rNULL__
if(ar__2(72:76) .EQ. '') e_TMB_1 = rNULL__
if(ar__2(78:82) .EQ. '') I_1 = rNULL__
if(ar__2(85:88) .EQ. '') CO2_1 = rNULL__
if(ar__2(90:94) .EQ. '') S60_CO_1 = rNULL__
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_1 .GT. -180) RAdeg_1=RAh_1*15.
if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
if(DEd_1 .GE. 0) DEdeg_1=DEd_1
if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,3)
+ No_2,Name_2,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Chem_1,
+ n_Chem_1,u_Chem_1,Region_2,Trans_1,Vlsr_1,Ve_1,TMB_1,e_TMB_1,
+ I_1,u_I_1,l_CO2_1,CO2_1,S60_CO_1,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 'table4.dat' ! CO/HCN line ratios
C Format for file interpretation
4 format(
+ A3,1X,A10,1X,A5,2X,A6,2X,A1,A1,3X,A1,F4.1,4X,F4.1,5X,A5,4X,A15)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,28
read(1,'(A77)')ar__3
read(ar__3,4)
+ No_3,Name_3,Region_3,LRS_1,Chem_2,u_Chem_2,l_CO1,CO1,CO2_2,
+ Ref,Com
if(ar__3(37:40) .EQ. '') CO1 = rNULL__
if(ar__3(45:48) .EQ. '') CO2_2 = rNULL__
c ..............Just test output...........
write(6,4)
+ No_3,Name_3,Region_3,LRS_1,Chem_2,u_Chem_2,l_CO1,CO1,CO2_2,
+ Ref,Com
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Estimates of mass-loss rates
C Format for file interpretation
5 format(
+ A3,1X,A10,1X,A5,2X,A6,1X,A1,A1,3X,F4.2,2X,A1,3X,F4.2,A1,3X,
+ F4.2,3X,A3,3X,F6.4,A1,3X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,24
read(1,'(A80)')ar__4
read(ar__4,5)
+ No_4,Name_4,Region_4,LRS_2,Chem_3,u_Chem_3,BC,r_BC,D,n_D,
+ dM1_dt,r_dM1_dt,dM2_dt,f_dM2_dt,dMIR_dt
if(ar__4(53:56) .EQ. '') dM1_dt = rNULL__
if(ar__4(66:71) .EQ. '') dM2_dt = rNULL__
if(ar__4(76:80) .EQ. '') dMIR_dt = rNULL__
c ..............Just test output...........
write(6,5)
+ No_4,Name_4,Region_4,LRS_2,Chem_3,u_Chem_3,BC,r_BC,D,n_D,
+ dM1_dt,r_dM1_dt,dM2_dt,f_dM2_dt,dMIR_dt
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea1.dat' ! Master list of the main sample
C Format for file interpretation
6 format(
+ A1,A3,2X,A10,2X,F5.2,A1,1X,F5.2,A1,1X,E5.2,3X,I2,3X,A7,1X,A1,
+ 3X,A4,2X,A6,A1,1X,A1,A1,3X,A1,A1,3X,A1,2X,A29)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,228
read(1,'(A113)')ar__5
read(ar__5,6)
+ n_No,No_5,IRAS,C21,u_C21,C32,u_C32,FIRAS_1,Var_,Region_5,OH,
+ r_OH,LRS_3,u_LRS,Chem_4,u_Chem_4,CO,HCN,H2O,Rem_2
if(ar__5(33:37) .EQ. '') FIRAS_1 = rNULL__
if(ar__5(41:42) .EQ. '') Var_ = iNULL__
c ..............Just test output...........
write(6,6)
+ n_No,No_5,IRAS,C21,u_C21,C32,u_C32,FIRAS_1,Var_,Region_5,OH,
+ r_OH,LRS_3,u_LRS,Chem_4,u_Chem_4,CO,HCN,H2O,Rem_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Additional sources
C Format for file interpretation
7 format(
+ A1,A3,2X,A10,2X,F5.2,A1,1X,F5.2,A1,1X,E5.2,3X,I2,3X,A7,1X,A1,
+ 3X,A4,2X,A6,A1,1X,A1,A1,3X,A1,A1,3X,A1,2X,A29)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,18
read(1,'(A113)')ar__6
read(ar__6,7)
+ n_No_1,No_6,IRAS_1,C21_1,u_C21_1,C32_1,u_C32_1,FIRAS_2,Var__1,
+ Region_6,OH_1,r_OH_1,LRS_4,u_LRS_1,Chem_5,u_Chem_5,CO_1,HCN_1,
+ H2O_1,Rem_3
if(ar__6(33:37) .EQ. '') FIRAS_2 = rNULL__
if(ar__6(41:42) .EQ. '') Var__1 = iNULL__
c ..............Just test output...........
write(6,7)
+ n_No_1,No_6,IRAS_1,C21_1,u_C21_1,C32_1,u_C32_1,FIRAS_2,Var__1,
+ Region_6,OH_1,r_OH_1,LRS_4,u_LRS_1,Chem_5,u_Chem_5,CO_1,HCN_1,
+ H2O_1,Rem_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end