Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/251/26 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Oct-14
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/251/26 CO outflow cand. toward W3/4/5 complex. II. (Li+, 2020)
*================================================================================
*CO outflow candidates toward the W3/4/5 complex.
*II. Feedback from candidate outflows.
* Li Y., Xu Ye, Sun Y., Yang Ji
* <Astrophys. J. Suppl. Ser., 251, 26 (2020)>
* =2020ApJS..251...26L
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Correlation among clouds and outflow candidates
integer*4 nr__
parameter (nr__=268) ! Number of records
character*301 ar__ ! Full-size record
C Position composed of: GLON GLAT
character*3 Reg (nr__) ! sub-region (1)
character*6 Seq (nr__) ! Cloud Index
real*8 GLON (nr__) ! (deg) [129.9/140]? Degree of Galactic Longitude
real*8 GLAT (nr__) ! (deg) [-4.5/4.9]? Degree of Galactic Latitude
real*4 Vlow (nr__) ! (km/s) [-66/0]? Lower velocity of integral interval
real*4 Vhigh (nr__) ! (km/s) [-58/8]? Upper velocity of integral interval
real*4 Mflow (nr__) ! (Msun) [0.03/186.3] Mass of outflow candidate (2)
real*8 Pflow (nr__) ! (Msun.km/s) [0.05/1470.5] Momentum of outflow candidate (2)
real*4 Eflow (nr__) ! (10-7J) [9.2e+41/1.4e+47] Kinetic energy of outflow
* candidate in units of erg (2)
real*4 Lflow (nr__) ! (10-7W) [2.9e+29/4.5e+34] Luminosity of outflow
* candidate in units of erg/s (2)
character*229 COI (nr__) ! Candidate outflow index(es)
*Note (1): Sub-regions as follows:
* Per = Perseus arm (cloud indexes 1-136 and "Median")
* Loc = Local arm (cloud indexes 137-260 and "Median")
* Int = Interarm 1 (cloud indexes 261-265 and "Median")
*Note (2): The total outflow mass, momentum, kinetic energy and luminosity
* corresponding to each cloud. They are from (Li+ 2019, J/ApJS/242/19),
* but here we newly apply the corrections for inclination angle,
* blending effect and optical depth. All values have been multiplied by
* 1.36 to take into consideration the mean molecular weight
* (Brunt 2010A&A...513A..67B) instead of the weight of only molecular H2
* given in (Li+ 2019, J/ApJS/242/19). The Median values are presented in
* the end of these three sub-regions with index of "Median".
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Global physical properties of the clouds
integer*4 nr__1
parameter (nr__1=277) ! Number of records
character*93 ar__1 ! Full-size record
character*3 Reg_1 (nr__1) ! sub-region (G1)
character*6 Seq_1 (nr__1) ! Cloud Index
real*8 Mcloud (nr__1) ! (Msun) [0.1/44942.5] Mass of cloud
real*4 Rcloud (nr__1) ! (pc) [0.04/17.12] Radius of cloud
real*4 DVcloud (nr__1) ! (km/s) [0.5/4] Line FWHM of cloud
real*4 e_DVcloud (nr__1) ! (km/s) [0.02/0.2] Uncertainty in DVcloud
real*4 Tcloud (nr__1) ! (K) [5.3/18.3] Excitation temperature of
* cloud
real*4 Taucloud (nr__1) ! [0.1/0.7] Optical Depth of cloud
real*4 Eturb (nr__1) ! (10-7J) [6.6e+41/3.7e+48] Turbulent energy of
* cloud in units of erg
real*8 Pturb (nr__1) ! (Msun.km/s) [0.1/129467.5] Turbulent momentum of
* cloud
real*4 tff (nr__1) ! (Myr) [0.5/5.5] Free-fall timescale of cloud
real*4 tdiss (nr__1) ! (Myr) [0.2/7.1] Turbulent dissipation
* timescale of cloud
real*4 Lturb (nr__1) ! (10-7W) [1.7e+28/1.4e+35] Turbulent dissipation
* rate of cloud in units of erg/s
real*4 vesc (nr__1) ! (km/s) [0.1/4.8] Escape velocity of cloud
real*4 Mesc (nr__1) ! (Msun) [0.2/309.4] Escape mass of cloud
real*4 Egrav (nr__1) ! (10-7J) [1.3e+40/1e+49] Gravitational binding
* energy of cloud in units of erg
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! The effect of feedback of outflow candidates on
their parent cloud
integer*4 nr__2
parameter (nr__2=277) ! Number of records
character*55 ar__2 ! Full-size record
character*3 Reg_2 (nr__2) ! sub-region (G1)
character*6 Seq_2 (nr__2) ! Cloud Index
real*4 Ef_Et (nr__2) ! [0.0008/47] Ratio of outflow to turbulent
* energies
real*4 Pf_Pt (nr__2) ! [0.00039/12] Ratio of outflow to turbulent
* momenta
real*4 Lf_Lt (nr__2) ! [0.003/350] Ratio of outflow to turbulent
* luminosities
real*4 Ef_Eg (nr__2) ! [0.0017/2400] Ratio of outflow to gravitational
* binding energies
real*4 Me_Mc (nr__2) ! [0.0005/84] Ratio of escape to total cloud
* masses
real*4 Me_Mf (nr__2) ! [1.6/34] Ratio of escape to outflow masses
C=============================================================================
C Loading file 'table1.dat' ! Correlation among clouds and outflow candidates
C Format for file interpretation
1 format(
+ A3,1X,A6,1X,F7.3,1X,F7.3,1X,F5.1,1X,F5.1,1X,F6.2,1X,F7.2,1X,
+ E8.2,1X,E8.2,1X,A229)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,268
read(1,'(A301)')ar__
read(ar__,1)
+ Reg(i__),Seq(i__),GLON(i__),GLAT(i__),Vlow(i__),Vhigh(i__),
+ Mflow(i__),Pflow(i__),Eflow(i__),Lflow(i__),COI(i__)
if(ar__(12:18) .EQ. '') GLON(i__) = rNULL__
if(ar__(20:26) .EQ. '') GLAT(i__) = rNULL__
if(ar__(28:32) .EQ. '') Vlow(i__) = rNULL__
if(ar__(34:38) .EQ. '') Vhigh(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Reg(i__),Seq(i__),GLON(i__),GLAT(i__),Vlow(i__),Vhigh(i__),
+ Mflow(i__),Pflow(i__),Eflow(i__),Lflow(i__),COI(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Global physical properties of the clouds
C Format for file interpretation
2 format(
+ A3,1X,A6,1X,F7.1,1X,F5.2,1X,F3.1,1X,F4.2,1X,F4.1,1X,F3.1,1X,
+ E7.1,1X,F8.1,1X,F3.1,1X,F3.1,1X,E7.1,1X,F3.1,1X,F5.1,1X,E7.1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,277
read(1,'(A93)')ar__1
read(ar__1,2)
+ Reg_1(i__),Seq_1(i__),Mcloud(i__),Rcloud(i__),DVcloud(i__),
+ e_DVcloud(i__),Tcloud(i__),Taucloud(i__),Eturb(i__),
+ Pturb(i__),tff(i__),tdiss(i__),Lturb(i__),vesc(i__),Mesc(i__),
+ Egrav(i__)
c ..............Just test output...........
write(6,2)
+ Reg_1(i__),Seq_1(i__),Mcloud(i__),Rcloud(i__),DVcloud(i__),
+ e_DVcloud(i__),Tcloud(i__),Taucloud(i__),Eturb(i__),
+ Pturb(i__),tff(i__),tdiss(i__),Lturb(i__),vesc(i__),Mesc(i__),
+ Egrav(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! The effect of feedback of outflow candidates on
* their parent cloud
C Format for file interpretation
3 format(A3,1X,A6,1X,E7.1,1X,E7.1,1X,E7.1,1X,E7.1,1X,E7.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,277
read(1,'(A55)')ar__2
read(ar__2,3)
+ Reg_2(i__),Seq_2(i__),Ef_Et(i__),Pf_Pt(i__),Lf_Lt(i__),
+ Ef_Eg(i__),Me_Mc(i__),Me_Mf(i__)
c ..............Just test output...........
write(6,3)
+ Reg_2(i__),Seq_2(i__),Ef_Et(i__),Pf_Pt(i__),Lf_Lt(i__),
+ Ef_Eg(i__),Me_Mc(i__),Me_Mf(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end