Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/267/515 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-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 (nr__) ! Identification number in tablea1 and tablea2
character*10 Name (nr__) ! Usual name
character*1 Sp (nr__) ! Spectral type
character*5 Region (nr__) ! Region in the IRAS color-color diagram (G1)
character*3 LRS (nr__) ! LRS spectral class
integer*4 FIRAS (nr__) ! (solLum) IRAS total flux expressed in solar
* luminosity at 1kpc
character*1 n (nr__) ! [+?] 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 (nr__1) ! (deg) Right Ascension 1950
real*8 DEdeg (nr__1) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*3 No_1 (nr__1) ! Identification number in tablea1 and tablea2
character*10 Name_1 (nr__1) ! Usual name
integer*4 RAh (nr__1) ! (h) ? Right ascension 1950
integer*4 RAm (nr__1) ! (min) ? Right ascension 1950
real*4 RAs (nr__1) ! (s) ? Right ascension 1950
character*1 DE_ (nr__1) ! ? Declination sign
integer*4 DEd (nr__1) ! (deg) ? Declination 1950
integer*4 DEm (nr__1) ! (arcmin) ? Declination 1950
integer*4 DEs (nr__1) ! (arcsec) ? Declination 1950
character*1 Chem (nr__1) ! [COS] Chemical type (C, O or S)
character*1 n_Chem (nr__1) ! [*#] (1)
character*1 u_Chem (nr__1) ! [?] uncertainty flag on chemical type
character*5 Region_1 (nr__1) ! Region in the IRAS color-color diagram (G1)
character*9 Trans (nr__1) ! Observed transitions
integer*4 Vlsr (nr__1) ! (km/s) ? Star Local Standard of Rest velocity
real*4 Ve (nr__1) ! (km/s) ? Expansion velocity as the half width at
* zero power
real*4 TMB (nr__1) ! (K) ? Main-beam brightness temperature measured
* at the peak of the line
real*4 e_TMB (nr__1) ! (K) ? rms noise of T_MB
real*4 I (nr__1) ! (K.km/s) ? Integrated area
character*1 u_I (nr__1) ! [?] uncertainty flag on I
character*1 l_CO2 (nr__1) ! limit flag on ICO/IHCN
real*4 CO2 (nr__1) ! ? Intensity ratio CO(2-1)/HCN(1-0)
real*4 S60_CO (nr__1) ! ? Ratio 60{mu}m to CO(2-1) emissions
character*24 Rem (nr__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 '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 (nr__2) ! (deg) Right Ascension 1950
real*8 DEdeg_1 (nr__2) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*3 No_2 (nr__2) ! Identification number in tablea1 and tablea2
character*10 Name_2 (nr__2) ! Usual name
integer*4 RAh_1 (nr__2) ! (h) ? Right ascension 1950
integer*4 RAm_1 (nr__2) ! (min) ? Right ascension 1950
real*4 RAs_1 (nr__2) ! (s) ? Right ascension 1950
character*1 DE__1 (nr__2) ! ? Declination sign
integer*4 DEd_1 (nr__2) ! (deg) ? Declination 1950
integer*4 DEm_1 (nr__2) ! (arcmin) ? Declination 1950
integer*4 DEs_1 (nr__2) ! (arcsec) ? Declination 1950
character*1 Chem_1 (nr__2) ! [COS] Chemical type (C, O or S)
character*1 n_Chem_1 (nr__2) ! [*#] (1)
character*1 u_Chem_1 (nr__2) ! [?] uncertainty flag on chemical type
character*5 Region_2 (nr__2) ! Region in the IRAS color-color diagram (G1)
character*9 Trans_1 (nr__2) ! Observed transitions
integer*4 Vlsr_1 (nr__2) ! (km/s) ? Star Local Standard of Rest velocity
real*4 Ve_1 (nr__2) ! (km/s) ? Expansion velocity as the half width at
* zero power
real*4 TMB_1 (nr__2) ! (K) ? Main-beam brightness temperature measured
* at the peak of the line
real*4 e_TMB_1 (nr__2) ! (K) ? rms noise of T_MB
real*4 I_1 (nr__2) ! (K.km/s) ? Integrated area
character*1 u_I_1 (nr__2) ! [?] uncertainty flag on I
character*1 l_CO2_1 (nr__2) ! limit flag on ICO/IHCN
real*4 CO2_1 (nr__2) ! ? Intensity ratio CO(2-1)/HCN(1-0)
real*4 S60_CO_1 (nr__2) ! ? Ratio 60{mu}m to CO(2-1) emissions
character*24 Rem_1 (nr__2) ! 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 (nr__3) ! Identification number in tablea1 & tablea2
character*10 Name_3 (nr__3) ! Usual name
character*5 Region_3 (nr__3) ! Region in the IRAS color-color diagram (G1)
character*6 LRS_1 (nr__3) ! LRS spectral class
character*1 Chem_2 (nr__3) ! [COS] Chemical type
character*1 u_Chem_2 (nr__3) ! [?] uncertainty flag on chemical type
character*1 l_CO1 (nr__3) ! [~] limit flag on CO(1-0)/HCN
real*4 CO1 (nr__3) ! ? CO(1-0)/HCN line ratio
real*4 CO2_2 (nr__3) ! ? CO(2-1)/HCN line ratio
character*5 Ref (nr__3) ! References (2)
character*15 Com (nr__3) ! 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 (nr__4) ! Identification number in tablea1 & tablea2
character*10 Name_4 (nr__4) ! Usual name
character*5 Region_4 (nr__4) ! Region in the IRAS color-color diagram (G1)
character*6 LRS_2 (nr__4) ! LRS spectral class
character*1 Chem_3 (nr__4) ! [COS] Chemical type
character*1 u_Chem_3 (nr__4) ! [?] uncertainty flag on chemical type
real*4 BC (nr__4) ! Bolometric correction F_tot_/F(IRAS)
character*1 r_BC (nr__4) ! references for Ftot (2)
real*4 D (nr__4) ! (kpc) Distance
character*1 n_D (nr__4) ! note for D (3)
real*4 dM1_dt (nr__4) ! (10-5Msun/yr) ? mass-loss rate from CO(1-0) (4)
character*3 r_dM1_dt (nr__4) ! references for J=1-0 CO observations
real*4 dM2_dt (nr__4) ! (10-5Msun/yr) ? mass-loss rate from CO(2-1) (4)
character*1 f_dM2_dt (nr__4) ! [pq] problematic dMCO(2-1)/dt (3)
real*4 dMIR_dt (nr__4) ! ? 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 (nr__5) ! [*] A "*" preceding the number indicates that
* the source was observed in the millimeter
* lines in the present study
character*3 No_5 (nr__5) ! Classification number
character*10 IRAS (nr__5) ! IRAS-PSC (Cat. II/125) name
real*4 C21 (nr__5) ! 25/12{mu}m index, C21=log(12F_25_/25F_12_)
character*1 u_C21 (nr__5) ! [)] uncertainty flag on C21
real*4 C32 (nr__5) ! 60/25{mu}m index, C32=log(25F_60_/12F_25_)
character*1 u_C32 (nr__5) ! [)] uncertainty flag on C32
real*4 FIRAS_1 (nr__5) ! (solLum) ? IRAS total flux expressed in solar
* luminosity at 1kpc
integer*4 Var_ (nr__5) ! ? variability index from the IRAS-PSC
* catalogue (Cat. II/125)
character*7 Region_5 (nr__5) ! Region of the IRAS color-color diagram (G1)
character*1 OH (nr__5) ! [-+.] OH maser detection (2)
character*4 r_OH (nr__5) ! References for OH maser observations (3)
character*6 LRS_3 (nr__5) ! LRS spectral class (4)
character*1 u_LRS (nr__5) ! [?] uncertainty flag on LRS
character*1 Chem_4 (nr__5) ! [COS] Proposed chemical type (C, O or S)
character*1 u_Chem_4 (nr__5) ! [?] uncertainty flag on chemical type
character*1 CO (nr__5) ! CO detections (2)
character*1 HCN (nr__5) ! [-+.*] HCN detections (2)
character*1 H2O (nr__5) ! [-+.] H2O detections (+) or non-detections (-)
* from the catalogue of Comoretto et al.
* (1990A&AS...84..179C)
character*29 Rem_2 (nr__5) ! 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 (nr__6) ! [*] A "*" preceding the number indicates that
* the source was observed in the millimeter
* lines in the present study
character*3 No_6 (nr__6) ! Classification number
character*10 IRAS_1 (nr__6) ! IRAS-PSC (Cat. II/125) name
real*4 C21_1 (nr__6) ! 25/12{mu}m index, C21=log(12F_25_/25F_12_)
character*1 u_C21_1 (nr__6) ! [)] uncertainty flag on C21
real*4 C32_1 (nr__6) ! 60/25{mu}m index, C32=log(25F_60_/12F_25_)
character*1 u_C32_1 (nr__6) ! [)] uncertainty flag on C32
real*4 FIRAS_2 (nr__6) ! (solLum) ? IRAS total flux expressed in solar
* luminosity at 1kpc
integer*4 Var__1 (nr__6) ! ? variability index from the IRAS-PSC
* catalogue (Cat. II/125)
character*7 Region_6 (nr__6) ! Region of the IRAS color-color diagram (G1)
character*1 OH_1 (nr__6) ! [-+.] OH maser detection (2)
character*4 r_OH_1 (nr__6) ! References for OH maser observations (3)
character*6 LRS_4 (nr__6) ! LRS spectral class (4)
character*1 u_LRS_1 (nr__6) ! [?] uncertainty flag on LRS
character*1 Chem_5 (nr__6) ! [COS] Proposed chemical type (C, O or S)
character*1 u_Chem_5 (nr__6) ! [?] uncertainty flag on chemical type
character*1 CO_1 (nr__6) ! CO detections (2)
character*1 HCN_1 (nr__6) ! [-+.*] HCN detections (2)
character*1 H2O_1 (nr__6) ! [-+.] H2O detections (+) or non-detections (-)
* from the catalogue of Comoretto et al.
* (1990A&AS...84..179C)
character*29 Rem_3 (nr__6) ! 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(i__),Name(i__),Sp(i__),Region(i__),LRS(i__),FIRAS(i__),
+ n(i__)
c ..............Just test output...........
write(6,1)
+ No(i__),Name(i__),Sp(i__),Region(i__),LRS(i__),FIRAS(i__),
+ n(i__)
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(i__),Name_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Chem(i__),n_Chem(i__),u_Chem(i__),
+ Region_1(i__),Trans(i__),Vlsr(i__),Ve(i__),TMB(i__),
+ e_TMB(i__),I(i__),u_I(i__),l_CO2(i__),CO2(i__),S60_CO(i__),
+ Rem(i__)
if(ar__1(16:17) .EQ. '') RAh(i__) = iNULL__
if(ar__1(19:20) .EQ. '') RAm(i__) = iNULL__
if(ar__1(22:25) .EQ. '') RAs(i__) = rNULL__
if(ar__1(28:29) .EQ. '') DEd(i__) = iNULL__
if(ar__1(31:32) .EQ. '') DEm(i__) = iNULL__
if(ar__1(34:35) .EQ. '') DEs(i__) = iNULL__
if(ar__1(57:59) .EQ. '') Vlsr(i__) = iNULL__
if(ar__1(61:64) .EQ. '') Ve(i__) = rNULL__
if(ar__1(66:70) .EQ. '') TMB(i__) = rNULL__
if(ar__1(72:76) .EQ. '') e_TMB(i__) = rNULL__
if(ar__1(78:82) .EQ. '') I(i__) = rNULL__
if(ar__1(85:88) .EQ. '') CO2(i__) = rNULL__
if(ar__1(90:94) .EQ. '') S60_CO(i__) = rNULL__
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(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,2)
+ No_1(i__),Name_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Chem(i__),n_Chem(i__),u_Chem(i__),
+ Region_1(i__),Trans(i__),Vlsr(i__),Ve(i__),TMB(i__),
+ e_TMB(i__),I(i__),u_I(i__),l_CO2(i__),CO2(i__),S60_CO(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.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(i__),Name_2(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Chem_1(i__),
+ n_Chem_1(i__),u_Chem_1(i__),Region_2(i__),Trans_1(i__),
+ Vlsr_1(i__),Ve_1(i__),TMB_1(i__),e_TMB_1(i__),I_1(i__),
+ u_I_1(i__),l_CO2_1(i__),CO2_1(i__),S60_CO_1(i__),Rem_1(i__)
if(ar__2(16:17) .EQ. '') RAh_1(i__) = iNULL__
if(ar__2(19:20) .EQ. '') RAm_1(i__) = iNULL__
if(ar__2(22:25) .EQ. '') RAs_1(i__) = rNULL__
if(ar__2(28:29) .EQ. '') DEd_1(i__) = iNULL__
if(ar__2(31:32) .EQ. '') DEm_1(i__) = iNULL__
if(ar__2(34:35) .EQ. '') DEs_1(i__) = iNULL__
if(ar__2(57:59) .EQ. '') Vlsr_1(i__) = iNULL__
if(ar__2(61:64) .EQ. '') Ve_1(i__) = rNULL__
if(ar__2(66:70) .EQ. '') TMB_1(i__) = rNULL__
if(ar__2(72:76) .EQ. '') e_TMB_1(i__) = rNULL__
if(ar__2(78:82) .EQ. '') I_1(i__) = rNULL__
if(ar__2(85:88) .EQ. '') CO2_1(i__) = rNULL__
if(ar__2(90:94) .EQ. '') S60_CO_1(i__) = rNULL__
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_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,3)
+ No_2(i__),Name_2(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Chem_1(i__),
+ n_Chem_1(i__),u_Chem_1(i__),Region_2(i__),Trans_1(i__),
+ Vlsr_1(i__),Ve_1(i__),TMB_1(i__),e_TMB_1(i__),I_1(i__),
+ u_I_1(i__),l_CO2_1(i__),CO2_1(i__),S60_CO_1(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 '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(i__),Name_3(i__),Region_3(i__),LRS_1(i__),Chem_2(i__),
+ u_Chem_2(i__),l_CO1(i__),CO1(i__),CO2_2(i__),Ref(i__),Com(i__)
if(ar__3(37:40) .EQ. '') CO1(i__) = rNULL__
if(ar__3(45:48) .EQ. '') CO2_2(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ No_3(i__),Name_3(i__),Region_3(i__),LRS_1(i__),Chem_2(i__),
+ u_Chem_2(i__),l_CO1(i__),CO1(i__),CO2_2(i__),Ref(i__),Com(i__)
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(i__),Name_4(i__),Region_4(i__),LRS_2(i__),Chem_3(i__),
+ u_Chem_3(i__),BC(i__),r_BC(i__),D(i__),n_D(i__),dM1_dt(i__),
+ r_dM1_dt(i__),dM2_dt(i__),f_dM2_dt(i__),dMIR_dt(i__)
if(ar__4(53:56) .EQ. '') dM1_dt(i__) = rNULL__
if(ar__4(66:71) .EQ. '') dM2_dt(i__) = rNULL__
if(ar__4(76:80) .EQ. '') dMIR_dt(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ No_4(i__),Name_4(i__),Region_4(i__),LRS_2(i__),Chem_3(i__),
+ u_Chem_3(i__),BC(i__),r_BC(i__),D(i__),n_D(i__),dM1_dt(i__),
+ r_dM1_dt(i__),dM2_dt(i__),f_dM2_dt(i__),dMIR_dt(i__)
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(i__),No_5(i__),IRAS(i__),C21(i__),u_C21(i__),C32(i__),
+ u_C32(i__),FIRAS_1(i__),Var_(i__),Region_5(i__),OH(i__),
+ r_OH(i__),LRS_3(i__),u_LRS(i__),Chem_4(i__),u_Chem_4(i__),
+ CO(i__),HCN(i__),H2O(i__),Rem_2(i__)
if(ar__5(33:37) .EQ. '') FIRAS_1(i__) = rNULL__
if(ar__5(41:42) .EQ. '') Var_(i__) = iNULL__
c ..............Just test output...........
write(6,6)
+ n_No(i__),No_5(i__),IRAS(i__),C21(i__),u_C21(i__),C32(i__),
+ u_C32(i__),FIRAS_1(i__),Var_(i__),Region_5(i__),OH(i__),
+ r_OH(i__),LRS_3(i__),u_LRS(i__),Chem_4(i__),u_Chem_4(i__),
+ CO(i__),HCN(i__),H2O(i__),Rem_2(i__)
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(i__),No_6(i__),IRAS_1(i__),C21_1(i__),u_C21_1(i__),
+ C32_1(i__),u_C32_1(i__),FIRAS_2(i__),Var__1(i__),
+ Region_6(i__),OH_1(i__),r_OH_1(i__),LRS_4(i__),u_LRS_1(i__),
+ Chem_5(i__),u_Chem_5(i__),CO_1(i__),HCN_1(i__),H2O_1(i__),
+ Rem_3(i__)
if(ar__6(33:37) .EQ. '') FIRAS_2(i__) = rNULL__
if(ar__6(41:42) .EQ. '') Var__1(i__) = iNULL__
c ..............Just test output...........
write(6,7)
+ n_No_1(i__),No_6(i__),IRAS_1(i__),C21_1(i__),u_C21_1(i__),
+ C32_1(i__),u_C32_1(i__),FIRAS_2(i__),Var__1(i__),
+ Region_6(i__),OH_1(i__),r_OH_1(i__),LRS_4(i__),u_LRS_1(i__),
+ Chem_5(i__),u_Chem_5(i__),CO_1(i__),HCN_1(i__),H2O_1(i__),
+ Rem_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end