Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/181/321 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-Oct-12
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/181/321 Properties of Spitzer c2d dark clouds (Evans+, 2009)
*================================================================================
*The Spitzer c2d legacy results:
*star-formation rates and efficiencies: evolution and lifetimes.
* Evans N.J., Dunham M.M., Jorgensen J.K., Enoch M.L., Merin B.,
* Van Dishoeck E.F., Alcala J.M., Myers P.C., Stapelfeldt K.R., Huard T.L.,
* Allen L.E., Harvey P.M., Van Kempen T., Blake G.A., Koerner D.W.,
* Mundy L.G., Padgett D.L., Sargent A.I.
* <Astrophys. J. Suppl. Ser., 181, 321-350 (2009)>
* =2009ApJS..181..321E
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! 2MASS and Spitzer flux densities of YSOs in the
c2d clouds
integer*4 nr__
parameter (nr__=1024) ! Number of records
character*173 ar__ ! Full-size record
integer*4 v_EDJ2009_ ! Running source index number
real*8 F1_25 ! (mJy) ? 2MASS J band flux density
real*8 e_F1_25 ! (mJy) ? Error in F1.25
real*8 F1_65 ! (mJy) ? 2MASS H band flux density
real*8 e_F1_65 ! (mJy) ? Error in F1.65
real*8 F2_17 ! (mJy) ? 2MASS K_S_ band flux density
real*8 e_F2_17 ! (mJy) ? Error in F2.17
real*8 F3_6 ! (mJy) ? IRAC 3.6 micron flux density
real*8 e_F3_6 ! (mJy) ? Error in F3.6
real*8 F4_5 ! (mJy) ? IRAC 4.5 micron flux density
real*8 e_F4_5 ! (mJy) ? Error in F4.5
real*8 F5_8 ! (mJy) ? IRAC 5.8 micron flux density
real*8 e_F5_8 ! (mJy) ? Error in F5.8
real*8 F8_0 ! (mJy) ? IRAC 8.0 micron flux density
real*8 e_F8_0 ! (mJy) ? Error in F8.0
real*8 F24 ! (mJy) ? MIPS 24 micron flux density
real*8 e_F24 ! (mJy) ? Error in F24
integer*4 F70 ! (mJy) ? MIPS 70 micron flux density
real*8 e_F70 ! (mJy) ? Error in F70
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! 0.36-0.96um flux densities of YSOs in the c2d clouds
integer*4 nr__1
parameter (nr__1=1024) ! Number of records
character*136 ar__1 ! Full-size record
integer*4 v_EDJ2009__1 ! Running source index number
real*4 F0_36 ! (mJy) ? The U band flux density
real*8 e_F0_36 ! (mJy) ? Error in F0.36
real*8 F0_44 ! (mJy) ? The B band flux density
real*8 e_F0_44 ! (mJy) ? Error in F0.44
real*4 F0_55 ! (mJy) ? The V band flux density
real*4 e_F0_55 ! (mJy) ? Error in F0.55
real*8 F0_64 ! (mJy) ? The R band flux density
real*8 e_F0_64 ! (mJy) ? Error in F0.64
real*8 F0_665 ! (mJy) ? The H{alpha}12 flux density
real*8 e_F0_665 ! (mJy) ? Error in F0.665
real*8 F0_79 ! (mJy) ? The I band flux density
real*8 e_F0_79 ! (mJy) ? Error in F0.79
real*8 F0_915 ! (mJy) ? The m914 band flux density
real*8 e_F0_915 ! (mJy) ? Error in F0.915
real*8 F0_96 ! (mJy) ? The z band flux density
real*8 e_F0_96 ! (mJy) ? Error in F0.96
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table10.dat' ! 3.4-100um flux densities of YSOs in the c2d clouds
integer*4 nr__2
parameter (nr__2=1024) ! Number of records
character*104 ar__2 ! Full-size record
integer*4 v_EDJ2009__2 ! Running source index number
integer*4 F3_4 ! (mJy) ? The L band flux density
real*4 e_F3_4 ! (mJy) ? Error in F3.4
integer*4 F5_0 ! (mJy) ? The M band flux density
integer*4 e_F5_0 ! (mJy) ? Error in F5.0
real*4 F6_7 ! (mJy) ? ISO 6.7 micron flux density
real*4 e_F6_7 ! (mJy) ? Error in F6.7
integer*4 F12 ! (mJy) ? IRAS 12 micron flux density
real*4 e_F12 ! (mJy) ? Error in F12
real*4 F14_3 ! (mJy) ? ISO 14.3 micron flux density
real*4 e_F14_3 ! (mJy) ? Error in F14.3
integer*4 F25 ! (mJy) ? IRAS 25 micron flux density
real*4 e_F25 ! (mJy) ? Error in F25
integer*4 F60 ! (mJy) ? IRAS 60 micron flux density
integer*4 e_F60 ! (mJy) ? Error in F60
integer*4 F100 ! (mJy) ? IRAS 100 micron flux density
integer*4 e_F100 ! (mJy) ? Error in F100
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! 160-1300um flux densities of YSOs in the c2d clouds
integer*4 nr__3
parameter (nr__3=1024) ! Number of records
character*79 ar__3 ! Full-size record
integer*4 v_EDJ2009__3 ! Running source index number
integer*4 F160 ! (mJy) ? MIPS 160 micron flux density
integer*4 e_F160 ! (mJy) ? Error in F160
integer*4 F350 ! (mJy) ? The 350 micron flux density
integer*4 e_F350 ! (mJy) ? Error in F350
integer*4 F450 ! (mJy) ? The 450 micron flux density
integer*4 e_F450 ! (mJy) ? Error in F450
integer*4 F850 ! (mJy) ? The 850 micron flux density
integer*4 e_F850 ! (mJy) ? Error in F850
integer*4 F1100 ! (mJy) ? The 1100 micron flux density
real*4 e_F1100 ! (mJy) ? Error in F1100
integer*4 F1200 ! (mJy) ? The 1200 micron flux density
integer*4 e_F1200 ! (mJy) ? Error in F1200
integer*4 F1300 ! (mJy) ? The 1300 micron flux density
real*4 e_F1300 ! (mJy) ? Error in F1300
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table12.dat' ! Properties of YSOs in the c2d clouds
integer*4 nr__4
parameter (nr__4=1024) ! Number of records
character*112 ar__4 ! Full-size record
integer*4 v_EDJ2009__4 ! Running source index number
character*7 Cloud ! Cloud name containing the source
character*19 SSTc2d ! Spitzer source name (JHHMMSS.ss+DDMMSS.s)
character*19 Class ! c2d classification (1)
real*4 AV ! (mag) The V band extinction (2)
real*4 alpha ! ? Observed 2-24{mu}m spectral index
character*2 l_Tbol ! [<>= ] Limit flag on Tbol
integer*4 Tbol ! (K) ? Observed bolometric temperature
character*2 l_Lbol ! [<>= ] Limit flag on Lbol
real*8 Lbol ! (solLum) ? Observed bolometric luminosity
real*4 alphac ! ? Corrected 2-24{mu}m spectral index
character*2 l_Tbolc ! [<>= ] Limit flag on Tbolc
integer*4 Tbolc ! (K) ? Corrected bolometric temperature
character*2 l_Lbolc ! [<>= ] Limit flag on Lbolc
real*8 Lbolc ! (solLum) ? Corrected bolometric luminosity
character*1 Env ! Envelope evidence? (3)
character*4 Qual ! Quality of calculated Lbol and Tbol (4)
*Note (1): Only sources with detections in at least three bands could be
* classified at all. Sources are labelled as follows:
* Galc = Galaxy candidates
* YSOc = Candidates for YSOs. Detections in all four IRAC bands and MIPS-1;
* see selection criteria in section 3.1.
* red = Flux densities at 24um that are at least 3 times the flux density from
* the nearest available IRAC band.
* PAH_em = Have colors indicative of a peak in the 8um band.
* star+dust = The SED is consistent with that of a stellar photosphere for
* wavelengths shorter than a particular band but an excess in that band
* of at least 3{sigma}. This band was appended to make a full name such
* as "YSOc_star+dust(MP1)" where MP1 implies the first MIPS band,
* at 24um or "YSOc_star+dust(IR4)" where IR4 implies the fourth
* IRAC band, at 8.0{mu}m.
* rising = Lack photometry at enough wavelengths to be classified as YSOc, but
* were later added as YSOs based on other information. The available
* flux densities rise significantly to longer wavelengths (15 cases)
* red1 = Lack photometry at enough wavelengths to be classified as YSOc, but
* were later added as YSOs based on other information. Source is
* detected only in IRAC-4 or MIPS-1 (3 cases).
*Note (2): Value of A_V_ used for de-reddening.
*Note (3): Indicates whether or not the source is associated with an
* envelope as traced by extended millimeter emission (based primarily
* on the sample of Enoch et al., 2008ApJ...684.1240E).
*Note (4): Quality flag as follows:
* G = Existing photometry provides sufficient spectral coverage for
* reliable calculations.
* L_I = the source is associated with an envelope and thus likely an
* embedded source, but lacks spectral coverage between 24{mu}m
* and 1mm. The calculated values of Lbol and Tbol are thus
* considered lower and upper limits, respectively.
* L_II = the source is not associated with an envelope and has an observed
* {alpha}{<=}-0.3, thus is not likely an embedded source, but no
* extinction correction could be derived based on available data,
* usually because of a lack of sufficient coverage in the near-IR.
* The calculated values of Lbol and Tbol are thus considered
* lower limits.
* O = Other. These three sources were added to the list of YSOs based
* on the searches for embedded protostars presented by Jorgensen
* et al. (2007ApJ...656..293J; for Perseus) and Jorgensen et al.
* (2008ApJ...683..822J; for Ophiuchus), but they are not associated
* with envelopes according to Enoch et al. (2008ApJ...684.1240).
* They are only detected at one wavelength (24{mu}m), thus
* calculations of {alpha}, Tbol, and Lbol are not possible. They are
* counted as Class I objects when classifying by {alpha} and {alpha}c;
* they are not counted at all when classifying by Tbol and Tbolc.
C=============================================================================
C Loading file 'table8.dat' ! 2MASS and Spitzer flux densities of YSOs in the
* c2d clouds
C Format for file interpretation
1 format(
+ I4,1X,F9.3,1X,F11.5,1X,F9.3,1X,F9.4,1X,F8.2,1X,F8.3,1X,F9.4,
+ 1X,F9.5,1X,F8.3,1X,F9.5,1X,F9.3,1X,F8.3,1X,F8.2,1X,F8.3,1X,
+ F7.2,1X,F8.3,1X,I6,1X,F8.1)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,1024
read(1,'(A173)')ar__
read(ar__,1)
+ v_EDJ2009_,F1_25,e_F1_25,F1_65,e_F1_65,F2_17,e_F2_17,F3_6,
+ e_F3_6,F4_5,e_F4_5,F5_8,e_F5_8,F8_0,e_F8_0,F24,e_F24,F70,e_F70
if(ar__(6:14) .EQ. '') F1_25 = rNULL__
if(ar__(16:26) .EQ. '') e_F1_25 = rNULL__
if(ar__(28:36) .EQ. '') F1_65 = rNULL__
if(ar__(38:46) .EQ. '') e_F1_65 = rNULL__
if(ar__(48:55) .EQ. '') F2_17 = rNULL__
if(ar__(57:64) .EQ. '') e_F2_17 = rNULL__
if(ar__(66:74) .EQ. '') F3_6 = rNULL__
if(ar__(76:84) .EQ. '') e_F3_6 = rNULL__
if(ar__(86:93) .EQ. '') F4_5 = rNULL__
if(ar__(95:103) .EQ. '') e_F4_5 = rNULL__
if(ar__(105:113) .EQ. '') F5_8 = rNULL__
if(ar__(115:122) .EQ. '') e_F5_8 = rNULL__
if(ar__(124:131) .EQ. '') F8_0 = rNULL__
if(ar__(133:140) .EQ. '') e_F8_0 = rNULL__
if(ar__(142:148) .EQ. '') F24 = rNULL__
if(ar__(150:157) .EQ. '') e_F24 = rNULL__
if(ar__(159:164) .EQ. '') F70 = iNULL__
if(ar__(166:173) .EQ. '') e_F70 = rNULL__
c ..............Just test output...........
write(6,1)
+ v_EDJ2009_,F1_25,e_F1_25,F1_65,e_F1_65,F2_17,e_F2_17,F3_6,
+ e_F3_6,F4_5,e_F4_5,F5_8,e_F5_8,F8_0,e_F8_0,F24,e_F24,F70,e_F70
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! 0.36-0.96um flux densities of YSOs in the c2d clouds
C Format for file interpretation
2 format(
+ I4,1X,F4.2,1X,F7.5,1X,F8.5,1X,F7.5,1X,F5.2,1X,F6.4,1X,F9.5,1X,
+ F8.5,1X,F8.5,1X,F8.5,1X,F7.3,1X,F8.5,1X,F7.3,1X,F8.5,1X,F7.3,
+ 1X,F9.5)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,1024
read(1,'(A136)')ar__1
read(ar__1,2)
+ v_EDJ2009__1,F0_36,e_F0_36,F0_44,e_F0_44,F0_55,e_F0_55,F0_64,
+ e_F0_64,F0_665,e_F0_665,F0_79,e_F0_79,F0_915,e_F0_915,F0_96,
+ e_F0_96
if(ar__1(6:9) .EQ. '') F0_36 = rNULL__
if(ar__1(11:17) .EQ. '') e_F0_36 = rNULL__
if(ar__1(19:26) .EQ. '') F0_44 = rNULL__
if(ar__1(28:34) .EQ. '') e_F0_44 = rNULL__
if(ar__1(36:40) .EQ. '') F0_55 = rNULL__
if(ar__1(42:47) .EQ. '') e_F0_55 = rNULL__
if(ar__1(49:57) .EQ. '') F0_64 = rNULL__
if(ar__1(59:66) .EQ. '') e_F0_64 = rNULL__
if(ar__1(68:75) .EQ. '') F0_665 = rNULL__
if(ar__1(77:84) .EQ. '') e_F0_665 = rNULL__
if(ar__1(86:92) .EQ. '') F0_79 = rNULL__
if(ar__1(94:101) .EQ. '') e_F0_79 = rNULL__
if(ar__1(103:109) .EQ. '') F0_915 = rNULL__
if(ar__1(111:118) .EQ. '') e_F0_915 = rNULL__
if(ar__1(120:126) .EQ. '') F0_96 = rNULL__
if(ar__1(128:136) .EQ. '') e_F0_96 = rNULL__
c ..............Just test output...........
write(6,2)
+ v_EDJ2009__1,F0_36,e_F0_36,F0_44,e_F0_44,F0_55,e_F0_55,F0_64,
+ e_F0_64,F0_665,e_F0_665,F0_79,e_F0_79,F0_915,e_F0_915,F0_96,
+ e_F0_96
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table10.dat' ! 3.4-100um flux densities of YSOs in the c2d clouds
C Format for file interpretation
3 format(
+ I4,1X,I5,1X,F5.1,1X,I5,1X,I3,1X,F5.1,1X,F5.2,1X,I5,1X,F6.1,1X,
+ F6.1,1X,F4.1,1X,I6,1X,F6.1,1X,I6,1X,I5,1X,I7,1X,I5)
C Effective file loading
open(unit=1,status='old',file=
+'table10.dat')
write(6,*) '....Loading file: table10.dat'
do i__=1,1024
read(1,'(A104)')ar__2
read(ar__2,3)
+ v_EDJ2009__2,F3_4,e_F3_4,F5_0,e_F5_0,F6_7,e_F6_7,F12,e_F12,
+ F14_3,e_F14_3,F25,e_F25,F60,e_F60,F100,e_F100
if(ar__2(6:10) .EQ. '') F3_4 = iNULL__
if(ar__2(12:16) .EQ. '') e_F3_4 = rNULL__
if(ar__2(18:22) .EQ. '') F5_0 = iNULL__
if(ar__2(24:26) .EQ. '') e_F5_0 = iNULL__
if(ar__2(28:32) .EQ. '') F6_7 = rNULL__
if(ar__2(34:38) .EQ. '') e_F6_7 = rNULL__
if(ar__2(40:44) .EQ. '') F12 = iNULL__
if(ar__2(46:51) .EQ. '') e_F12 = rNULL__
if(ar__2(53:58) .EQ. '') F14_3 = rNULL__
if(ar__2(60:63) .EQ. '') e_F14_3 = rNULL__
if(ar__2(65:70) .EQ. '') F25 = iNULL__
if(ar__2(72:77) .EQ. '') e_F25 = rNULL__
if(ar__2(79:84) .EQ. '') F60 = iNULL__
if(ar__2(86:90) .EQ. '') e_F60 = iNULL__
if(ar__2(92:98) .EQ. '') F100 = iNULL__
if(ar__2(100:104) .EQ. '') e_F100 = iNULL__
c ..............Just test output...........
write(6,3)
+ v_EDJ2009__2,F3_4,e_F3_4,F5_0,e_F5_0,F6_7,e_F6_7,F12,e_F12,
+ F14_3,e_F14_3,F25,e_F25,F60,e_F60,F100,e_F100
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! 160-1300um flux densities of YSOs in the c2d clouds
C Format for file interpretation
4 format(
+ I4,1X,I5,1X,I5,1X,I5,1X,I4,1X,I5,1X,I5,1X,I4,1X,I4,1X,I4,1X,
+ F6.1,1X,I4,1X,I3,1X,I3,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,1024
read(1,'(A79)')ar__3
read(ar__3,4)
+ v_EDJ2009__3,F160,e_F160,F350,e_F350,F450,e_F450,F850,e_F850,
+ F1100,e_F1100,F1200,e_F1200,F1300,e_F1300
if(ar__3(6:10) .EQ. '') F160 = iNULL__
if(ar__3(12:16) .EQ. '') e_F160 = iNULL__
if(ar__3(18:22) .EQ. '') F350 = iNULL__
if(ar__3(24:27) .EQ. '') e_F350 = iNULL__
if(ar__3(29:33) .EQ. '') F450 = iNULL__
if(ar__3(35:39) .EQ. '') e_F450 = iNULL__
if(ar__3(41:44) .EQ. '') F850 = iNULL__
if(ar__3(46:49) .EQ. '') e_F850 = iNULL__
if(ar__3(51:54) .EQ. '') F1100 = iNULL__
if(ar__3(56:61) .EQ. '') e_F1100 = rNULL__
if(ar__3(63:66) .EQ. '') F1200 = iNULL__
if(ar__3(68:70) .EQ. '') e_F1200 = iNULL__
if(ar__3(72:74) .EQ. '') F1300 = iNULL__
if(ar__3(76:79) .EQ. '') e_F1300 = rNULL__
c ..............Just test output...........
write(6,4)
+ v_EDJ2009__3,F160,e_F160,F350,e_F350,F450,e_F450,F850,e_F850,
+ F1100,e_F1100,F1200,e_F1200,F1300,e_F1300
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table12.dat' ! Properties of YSOs in the c2d clouds
C Format for file interpretation
5 format(
+ I4,1X,A7,1X,A19,1X,A19,1X,F4.1,1X,F5.2,1X,A2,I4,1X,A2,F8.5,1X,
+ F5.2,1X,A2,I4,1X,A2,F8.4,1X,A1,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table12.dat')
write(6,*) '....Loading file: table12.dat'
do i__=1,1024
read(1,'(A112)')ar__4
read(ar__4,5)
+ v_EDJ2009__4,Cloud,SSTc2d,Class,AV,alpha,l_Tbol,Tbol,l_Lbol,
+ Lbol,alphac,l_Tbolc,Tbolc,l_Lbolc,Lbolc,Env,Qual
if(ar__4(59:63) .EQ. '') alpha = rNULL__
if(ar__4(67:70) .EQ. '') Tbol = iNULL__
if(ar__4(74:81) .EQ. '') Lbol = rNULL__
if(ar__4(83:87) .EQ. '') alphac = rNULL__
if(ar__4(91:94) .EQ. '') Tbolc = iNULL__
if(ar__4(98:105) .EQ. '') Lbolc = rNULL__
c ..............Just test output...........
write(6,5)
+ v_EDJ2009__4,Cloud,SSTc2d,Class,AV,alpha,l_Tbol,Tbol,l_Lbol,
+ Lbol,alphac,l_Tbolc,Tbolc,l_Lbolc,Lbolc,Env,Qual
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end