Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/710/150 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.5, on 2013-May-26
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/ApJ/710/150 Molecular lines in EGOs (Chen+, 2010)
*================================================================================
*A search for infall evidence in EGOs.
*I. The northern sample.
* Chen X., Shen Z.-Q., Li J.-J., Xu Y., He J.-H.
* <Astrophys. J., 710, 150-169 (2010)>
* =2010ApJ...710..150C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Sample parameters
integer*4 nr__
parameter (nr__=88) ! Number of records
character*42 ar__ ! Full-size record
real*8 RA_ (nr__) ! (deg) Right Ascension J2000
real*8 Dec (nr__) ! (deg) Declination J2000
character*11 Name (nr__) ! Source name (GLL.ll+B.bb)
integer*4 RAh (nr__) ! (h) Hour of right ascension (J2000) (1)
integer*4 RAm (nr__) ! (min) Minute of right ascension (J2000)
real*4 RAs (nr__) ! (s) Second of right ascension (J2000)
character*1 DE_ (nr__) ! Sign of declination (J2000) (1)
integer*4 DEd (nr__) ! (deg) Degree of declination (J2000) (1)
integer*4 DEm (nr__) ! (arcmin) Arcminute of declination (J2000)
integer*4 DEs (nr__) ! (arcsec) Arcsecond of declination (J2000)
character*1 IRDC (nr__) ! [Y/N] Association with IR dark clouds (2)
character*1 CH3OH (nr__) ! [Y/N]? Associations with 6.7GHz class II
* methanol masers (3)
character*1 UCHII (nr__) ! [Y/N]? Associations with ultra-compact
* HII regions (3)
character*1 v1_1mm (nr__) ! [Y/N]? Associations with 1.1mm continuum
* sources (3)
integer*4 C08 (nr__) ! [1-5] Sources in Cyganowski et al. 2008 (4)
*Note (1): The targeted positions in the observations are from the EGO positions
* presented by the EGO catalog of Cyganowski et al. (2008,
* Cat. J/AJ/136/2391; <EGO GLLL.ll+B.bb(a)> in Simbad).
*Note (2): Association with IR dark clouds : "Y", Yes; "N", No; presented by
* Cyganowski et al. (2008, Cat. J/AJ/136/2391).
*Note (3): Associations with 6.7GHz class II methanol masers, UC HII regions and
* 1.1mm continuum sources within 30": "Y", Yes; "N", No; identified from
* the 6.7GHz maser catalogs (Szymczak et al. 2007A&A...468..617S;
* Cyganowski et al. 2008, Cat. J/AJ/136/2391, 2009, Cat. J/ApJ/702/1615;
* Caswell 2009PASA...26..454C; Xu et al. 2009A&A...507.1117X) and UC HII
* catalogs (Wood & Churchwell 1989ApJS...69..831W; Becker et al. 1994,
* Cat. J/ApJS/91/347; Kurtz et al. 1994, Cat. J/ApJS/91/659; Walsh et
* al. 1998, Cat. J/MNRAS/301/640; Forster & Caswell 2000ApJ...530..371F;
* Wu et al. 2007ApJ...669L..37W), and 1.1mm continuum BOLOCAM GPS
* archive (Rosolowsky et al. 2010, Cat. J/ApJS/188/123).
*Note (4): Remarks: "1"-"5" represent that the sources are selected from Tables
* 1-5 of Cyganowski et al. (2008, Cat. J/AJ/136/2391), respectively.
* "1", "2", and "5" are classified as "likely" outflow candidates, while
* 3 and 4 are classified as "possible" outflow candidates.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Observed HCO^+^(1-0) line properties of all 88
sample sources
integer*4 nr__1
parameter (nr__1=119) ! Number of records
character*68 ar__1 ! Full-size record
character*11 Name_1 (nr__1) ! Source name (GLL.ll+B.bb)
character*1 f_Name (nr__1) ! [bcd] Peculiar profile (1)
real*4 TdV (nr__1) ! (K.km/s) ? Integrated intensity ({int}T_MB_dV)
real*4 e_TdV (nr__1) ! (K.km/s) ? TdV uncertainty
real*4 Vlsr (nr__1) ! (km/s) ? Velocity at peak V_LSR_
real*4 e_Vlsr (nr__1) ! (km/s) ? Vlsr uncertainty
real*4 DelV (nr__1) ! (km/s) ? Line width {Delta}V
real*4 e_DelV (nr__1) ! (km/s) ? DelV uncertainty
real*4 Tmb (nr__1) ! (K) ? Peak main beam temperature T_MB_ (2)
real*4 Tmb2 (nr__1) ! (K) ? Real observed T_MB_ (3)
real*4 Vlsrpk (nr__1) ! (km/s) ? Observed V_LSR_ at HCO^+^ spectral peaks
* for double-peaked and skewed profiles (3)
real*4 rms (nr__1) ! (K) ? 1{sigma} noise in the observed HCO^+^
* spectrum at T_MB_ scale
*Note (1): Flag as follows:
* b = Sources with non-asymmetric profiles (21 in total).
* c = Sources with complex HCO^+^ spectral profiles. They were not used in
* analysis for the line asymmetry (see Section 3.2).
* d = No HCO^+^ emission detected (16 sources in total).
*Note (2): The peak main beam temperature T_MB_ estimated from Gaussian fits to
* HCO^+^ lines after masking the absorption dips and subtracting
* Gaussian fits to any line wings (see Table 5).
*Note (3): The peak position is used to analyze line asymmetry (see Section 3.2).
* For non-asymmetric profile sources, the real observed T_MB_ and V_LSR_
* adopt the Gaussian fit values in Columns "Tmb" and "Vlsr". Their
* uncertainties are assumed to be the 1{sigma} noise of HCO^+^ spectrum
* listed in Column "rms" and one-half the channel width, respectively.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! *Gaussian fits to HCO^+^ line wings
integer*4 nr__2
parameter (nr__2=31) ! Number of records
character*53 ar__2 ! Full-size record
character*11 Name_2 (nr__2) ! Source name (GLL.ll+B.bb)
character*1 f_Name_1 (nr__2) ! [c] Peculiar profile (1)
real*4 TdV_1 (nr__2) ! (K.km/s) Integrated intensity ({int}T_MB_dV)
real*4 e_TdV_1 (nr__2) ! (K.km/s) TdV uncertainty
real*4 Vlsr_1 (nr__2) ! (km/s) Velocity at peak V_LSR_
real*4 e_Vlsr_1 (nr__2) ! (km/s) Vlsr uncertainty
real*4 DelV_1 (nr__2) ! (km/s) Line width ({Delta}V)
real*4 e_DelV_1 (nr__2) ! (km/s) DelV uncertainty
real*4 Tmb_1 (nr__2) ! (K) Peak main beam temperature (T_MB_)
* determined from single broad Gaussian fit to
* HCO^+^ line wing
character*1 line (nr__2) ! [DBR] Double/Blue/Red line wing shape (2)
*Note (1):
* c = As there are strong emissions at the blue- and redshifted wings of this
* source (see Figure 1), two separate components (rather than a single
* broad component) were used to fit the strong blue- and redshifted wings.
*Note (2): Line wing shape classified as:
* D = the line wing emission is from both blue- and redshifted (double) wings
* B = the line wing emission is mainly from blueshifted wing
* R = the line wing emission is mainly from redshifted wing
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Observed CO(1-0) line properties of 72 sources
integer*4 nr__3
parameter (nr__3=77) ! Number of records
character*146 ar__3 ! Full-size record
character*11 Name_3 (nr__3) ! Source name (GLL.ll+B.bb)
character*1 f_Name_2 (nr__3) ! [ab] Flag on Name (1)
real*4 TdV12 (nr__3) ! (K.km/s) ^12^CO integrated intensity {int}T_MB_dV
real*4 e_TdV12 (nr__3) ! (K.km/s) TdV12 uncertainty
real*4 Vlsr12 (nr__3) ! (km/s) ^12^CO velocity at peak V_LSR_
real*4 e_Vlsr12 (nr__3) ! (km/s) Vlsr12 uncertainty
real*4 DelV12 (nr__3) ! (km/s) ^12^CO line width {Delta}V
real*4 e_DelV12 (nr__3) ! (km/s) DelV12 uncertainty
real*4 Tmb12 (nr__3) ! (K) Peak main beam temperature T_MB_ for ^12^CO
* line determined from a single-Gaussian fit
real*4 rms12 (nr__3) ! (K) ? 1{sigma} noise in the observed ^12^CO
* spectrum at T_MB_ scale
real*4 TdV13 (nr__3) ! (K.km/s) ^13^CO integrated intensity {int}T_MB_dV
real*4 e_TdV13 (nr__3) ! (K.km/s) TdV13 uncertainty
real*4 Vlsr13 (nr__3) ! (km/s) ^13^CO velocity at peak V_LSR_
real*4 e_Vlsr13 (nr__3) ! (km/s) Vlsr13 uncertainty
real*4 DelV13 (nr__3) ! (km/s) ^13^CO line width {Delta}V
real*4 e_DelV13 (nr__3) ! (km/s) DelV13 uncertainty
real*4 Tmb13 (nr__3) ! (K) Peak main beam temperature T_MB_ for ^13^CO
* line determined from a single-Gaussian fit
real*4 rms13 (nr__3) ! (K) ? 1{sigma} noise in the observed ^13^CO
* spectrum at T_MB_ scale
real*4 TdV18 (nr__3) ! (K.km/s) C^18^O integrated intensity {int}T_MB_dV
real*4 e_TdV18 (nr__3) ! (K.km/s) TdV18 uncertainty
real*4 Vlsr18 (nr__3) ! (km/s) C^18^O velocity at peak V_LSR_
real*4 e_Vlsr18 (nr__3) ! (km/s) Vlsr18 uncertainty
real*4 DelV18 (nr__3) ! (km/s) C^18^O line width {Delta}V
real*4 e_DelV18 (nr__3) ! (km/s) DelV18 uncertainty
real*4 Tmb18 (nr__3) ! (K) Peak main beam temperature T_MB_ for C^18^O
* line determined from a single-Gaussian fit
real*4 rms18 (nr__3) ! (K) ? 1{sigma} noise in the observed C^18^O
* spectrum at T_MB_ scale
*Note (1): Flag as follows:
* a = Sources with complex C^18^O spectral profiles. They were not used in
* analysis for the line asymmetry (see Section 3.2).
* b = ^12^CO was used to calculate the optical depth of C^18^O with typical
* LTE method for these sources (see Section 3.2).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Derived line asymmetric parameters and
source physical properties for 69 sources
integer*4 nr__4
parameter (nr__4=69) ! Number of records
character*73 ar__4 ! Full-size record
character*11 Name_4 (nr__4) ! Source name (GLL.ll+B.bb)
real*4 delV_2 (nr__4) ! HCO^+^ line asymmetry parameter {delta}V (2)
real*4 e_delV_2 (nr__4) ! delV uncertainty (1)
real*4 TB_TR (nr__4) ! ? HCO^+^ line asymmetry parameters determined
* by the T_MB_(B)/T_MB_(R) method (Section 3.2)
real*4 e_TB_TR (nr__4) ! ? TB/TR uncertainty (1)
character*1 P (nr__4) ! [NBR] Non-asymmetric, Blue or Red profile (2)
real*4 tau (nr__4) ! Optical depth of C^18^O line
real*4 Dist (nr__4) ! (kpc) Kinematic distance to source (3)
real*4 f (nr__4) ! Beam filling factor
real*4 Size (nr__4) ! (pc) Linear size
real*4 NC18O (nr__4) ! (10+16cm-2) Column density of C^18^O
real*4 NH2 (nr__4) ! (10+23cm-2) Column density of H_2_
real*4 nH2_1 (nr__4) ! (10+5cm-3) Volume density of H_2_
integer*4 Mass (nr__4) ! (Msun) Core mass
*Note (1): Uncertainties calculated from the rms uncertainties in temperatures
* and one-half the channel widths.
*Note (2): {delta}V is the ratio (V_thick_-V_thin_)/{Delta}V_thin_ where
* V_thick_ and V_thin_ are the peak velocities of the optically thick
* and optically thin lines, and {Delta}V_thin_ the line width of the
* optically thin line. We adopt {delta}V of +/-0.25 as the threshold
* to define the line asymmetry. For double-peaked profile, the velocity
* at the brightest peak is adopted as the value of V_thick_ to calculate
* {delta}V, and we consider an asymmetry to be significant if the
* difference (or sum) between T_MB_(B)/T_MB_(R) and its uncertainty is
* still larger (or less) than 1 for blue (or red) double-peaked profile.
* The line profile is then classified as Blue, Red or Non-asymmetric.
*Note (3): The near kinematic distance was adopted for source with near/far
* distance ambiguity. The distances to G49.07-0.33, G49.27-0.32,
* G49.27-0.34 and G59.79+0.63 are assumed to be 5kpc.
C=============================================================================
C Loading file 'table1.dat' ! Sample parameters
C Format for file interpretation
1 format(
+ A11,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A1,1X,A1,1X,
+ A1,1X,A1,1X,I1)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,88
read(1,'(A42)')ar__
read(ar__,1)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),IRDC(i__),CH3OH(i__),UCHII(i__),v1_1mm(i__),
+ C08(i__)
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_(i__) = RAh(i__)
if(RAh(i__) .GE. 0) RA_(i__)=RAh(i__)*15.
if(RAm(i__) .GE. 0) RA_(i__)=RA_(i__)+RAm(i__)/4.
if(RAs(i__) .GE. 0) RA_(i__)=RA_(i__)+RAs(i__)/240.
Dec(i__) = DEd(i__)
if(DEm(i__) .GE. 0) Dec(i__)=Dec(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) Dec(i__)=Dec(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),IRDC(i__),CH3OH(i__),UCHII(i__),v1_1mm(i__),
+ C08(i__)
write(6,'(6H Pos: 2F8.4)') RA_(i__),Dec(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Observed HCO^+^(1-0) line properties of all 88
* sample sources
C Format for file interpretation
2 format(
+ A11,A1,1X,F5.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F6.2,1X,F4.2)
C Effective file loading
open(unit=1,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,119
read(1,'(A68)')ar__1
read(ar__1,2)
+ Name_1(i__),f_Name(i__),TdV(i__),e_TdV(i__),Vlsr(i__),
+ e_Vlsr(i__),DelV(i__),e_DelV(i__),Tmb(i__),Tmb2(i__),
+ Vlsrpk(i__),rms(i__)
if(ar__1(14:18) .EQ. '') TdV(i__) = rNULL__
if(ar__1(20:23) .EQ. '') e_TdV(i__) = rNULL__
if(ar__1(25:30) .EQ. '') Vlsr(i__) = rNULL__
if(ar__1(32:35) .EQ. '') e_Vlsr(i__) = rNULL__
if(ar__1(37:41) .EQ. '') DelV(i__) = rNULL__
if(ar__1(43:46) .EQ. '') e_DelV(i__) = rNULL__
if(ar__1(48:51) .EQ. '') Tmb(i__) = rNULL__
if(ar__1(53:56) .EQ. '') Tmb2(i__) = rNULL__
if(ar__1(58:63) .EQ. '') Vlsrpk(i__) = rNULL__
if(ar__1(65:68) .EQ. '') rms(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),f_Name(i__),TdV(i__),e_TdV(i__),Vlsr(i__),
+ e_Vlsr(i__),DelV(i__),e_DelV(i__),Tmb(i__),Tmb2(i__),
+ Vlsrpk(i__),rms(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! *Gaussian fits to HCO^+^ line wings
C Format for file interpretation
3 format(
+ A11,A1,1X,F5.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,
+ F4.2,1X,A1)
C Effective file loading
open(unit=1,file='table4.dat', status='old')
write(6,*) '....Loading file: table4.dat'
do i__=1,31
read(1,'(A53)')ar__2
read(ar__2,3)
+ Name_2(i__),f_Name_1(i__),TdV_1(i__),e_TdV_1(i__),Vlsr_1(i__),
+ e_Vlsr_1(i__),DelV_1(i__),e_DelV_1(i__),Tmb_1(i__),line(i__)
c ..............Just test output...........
write(6,3)
+ Name_2(i__),f_Name_1(i__),TdV_1(i__),e_TdV_1(i__),Vlsr_1(i__),
+ e_Vlsr_1(i__),DelV_1(i__),e_DelV_1(i__),Tmb_1(i__),line(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Observed CO(1-0) line properties of 72 sources
C Format for file interpretation
4 format(
+ A11,A1,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,file='table6.dat', status='old')
write(6,*) '....Loading file: table6.dat'
do i__=1,77
read(1,'(A146)')ar__3
read(ar__3,4)
+ Name_3(i__),f_Name_2(i__),TdV12(i__),e_TdV12(i__),Vlsr12(i__),
+ e_Vlsr12(i__),DelV12(i__),e_DelV12(i__),Tmb12(i__),rms12(i__),
+ TdV13(i__),e_TdV13(i__),Vlsr13(i__),e_Vlsr13(i__),DelV13(i__),
+ e_DelV13(i__),Tmb13(i__),rms13(i__),TdV18(i__),e_TdV18(i__),
+ Vlsr18(i__),e_Vlsr18(i__),DelV18(i__),e_DelV18(i__),
+ Tmb18(i__),rms18(i__)
if(ar__3(55:58) .EQ. '') rms12(i__) = rNULL__
if(ar__3(99:102) .EQ. '') rms13(i__) = rNULL__
if(ar__3(143:146) .EQ. '') rms18(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3(i__),f_Name_2(i__),TdV12(i__),e_TdV12(i__),Vlsr12(i__),
+ e_Vlsr12(i__),DelV12(i__),e_DelV12(i__),Tmb12(i__),rms12(i__),
+ TdV13(i__),e_TdV13(i__),Vlsr13(i__),e_Vlsr13(i__),DelV13(i__),
+ e_DelV13(i__),Tmb13(i__),rms13(i__),TdV18(i__),e_TdV18(i__),
+ Vlsr18(i__),e_Vlsr18(i__),DelV18(i__),e_DelV18(i__),
+ Tmb18(i__),rms18(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Derived line asymmetric parameters and
* source physical properties for 69 sources
C Format for file interpretation
5 format(
+ A11,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A1,1X,F4.2,1X,F4.1,1X,
+ F4.2,1X,F4.2,1X,F4.1,1X,F3.1,1X,F3.1,1X,I5)
C Effective file loading
open(unit=1,file='table7.dat', status='old')
write(6,*) '....Loading file: table7.dat'
do i__=1,69
read(1,'(A73)')ar__4
read(ar__4,5)
+ Name_4(i__),delV_2(i__),e_delV_2(i__),TB_TR(i__),e_TB_TR(i__),
+ P(i__),tau(i__),Dist(i__),f(i__),Size(i__),NC18O(i__),
+ NH2(i__),nH2_1(i__),Mass(i__)
if(ar__4(24:27) .EQ. '') TB_TR(i__) = rNULL__
if(ar__4(29:32) .EQ. '') e_TB_TR(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Name_4(i__),delV_2(i__),e_delV_2(i__),TB_TR(i__),e_TB_TR(i__),
+ P(i__),tau(i__),Dist(i__),f(i__),Size(i__),NC18O(i__),
+ NH2(i__),nH2_1(i__),Mass(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end