Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/612/A79 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/A+A/612/A79 Barnard 30 dark cloud IR and submm data (Barrado+, 2018)
*================================================================================
*Early phases in the stellar and substellar formation and evolution.
*Infrared and submillimeter data in the Barnard 30 dark cloud.
* Barrado D., de Gregorio Monsalvo I., Huelamo N., Morales-Calderon M.,
* Bayo A., Palau A., Ruiz M.T., Riviere-Marichalar P., Bouy H., Morata O.,
* Stauffer J.R., Eiroa C., Noriega-Crespo A.
* <Astron. Astrophys. 612, A79 (2018)>
* =2018A&A...612A..79B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! APEX/LABOCA sources in Barnard 30
integer*4 nr__
parameter (nr__=34) ! Number of records
character*98 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*4 Name (nr__) ! B30 source ID (LBNN)
character*1 n_Name (nr__) ! [a-i] Note on Name (1)
integer*4 RAh (nr__) ! (h) Maximum emission of the core
* right ascension (J2000) (2)
integer*4 RAm (nr__) ! (min) Maximum emission of the core
* right ascension (J2000) (2)
real*4 RAs (nr__) ! (s) Maximum emission of the core
* right ascension (J2000) (2)
character*1 DE_ (nr__) ! Maximum emission of the core
* declination sign (J2000) (2)
integer*4 DEd (nr__) ! (deg) Maximum emission of the core
* declination (J2000) (2)
integer*4 DEm (nr__) ! (arcmin) Maximum emission of the core
* declination (J2000) (2)
real*4 DEs (nr__) ! (arcsec) Maximum emission of the core
* declination (J2000) (2)
integer*4 S870 (nr__) ! (mJy) APEX/LABOCA flux density at 870um computed
* above 3{sigma} emission (3)
integer*4 Ipeak (nr__) ! (mJy/beam) APEX/LABOCA peak intensity at 870um
* (maximum of the emission) (4)
real*4 AV (nr__) ! (mag) ?=- Values derived from the extinction map,
* with a resolution of 1.5 arcmin, three times
* larger the LABOCA beam-size
real*4 Mall (nr__) ! (Msun) Envelope masses (5)
real*4 Mpeak (nr__) ! (Msun) Envelope masses (5)
character*1 l_Lint (nr__) ! Limit flag on Lint
real*4 Lint (nr__) ! (Lsun) ?=- Internal luminosity of the core, based on
* MIPS flux at 70um, after Dunham et al.
* (2008, Cat. J/ApJS/179/249)
character*1 n_Lint (nr__) ! [j] Note on Lint (1)
character*1 l_F70 (nr__) ! Limit flag on F70
real*4 F70 (nr__) ! (Jy) ?=- MIPS flux at 70um
character*1 n_F70 (nr__) ! [j] Note on F70 (1)
integer*4 AngSize (nr__) ! (arcsec) ? Deconvolved size after fitting a Gaussian to
* the cores, the averaged of the major and
* minor axis is represented
character*2 n_AngSize (nr__) ! [PL ] "PL" stands for point-like source
integer*4 rms (nr__) ! (mJy) rms at 870um
character*3 Gr (nr__) ! Groups based on the detection at 70 and 24um
* and presence of counterparts, as discussed
* in Sect. 2.7 (6)
character*8 Type (nr__) ! Tentative classification (7)
character*1 BD_T (nr__) ! [Y ] Presence of BD candidates, Y indicates
* whether there is a BD candidate within
* APEX/LABOCA beam
character*6 BD_5 (nr__) ! Indicates whether there is a optical/IR
* counterpart within 5 arcsec from
* APEX/LABOCA peak (8)
*Note (1): Notes as follows:
* a = B30-LB03 and B30-LB04 very close to each other, see Sect. A.1.
* b = IRAS 05286+1203, RA=82.8671, Dec=+12.0899, [12]=0.779+/-0.0779Jy,
* [25]=1.91+/-0.2101Jy, [60]=9.35+/-1.7765Jy, [100]=70.4e+/-14.784Jy.
* LB30-LB19d has a MIPS 2 flux at 70um of 0.200 Jy, which translate to
* L(int)=0.094L_{sun}_ (i.e., a VeLLO). Note that within the LABOCA beam
* there is a APEX/SABOCA source detected at 350um
* (B30-SB08, Huelamo et al. 2017A&A...597A..17H).
* c = B30-LB20 contains a APEX/SABOCA at 350um
* (B30-SB09, Huelamo et al. 2017A&A...597A..17H).
* d = There are two APEX/SABOCA sources within the B30-LB22 beam
* (B30-SB03 and B30-SB04, Huelamo et al. 2017A&A...597A..17H).
* e = There are two APEX/SABOCA sources within the B30-LB23 beam
* (B30-SB05 and B30-SB06, Huelamo et al. 2017A&A...597A..17H).
* f = B30-LB27 contains a APEX/SABOCA at 350um
* (B30-SB12, Huelamo et al. 2017A&A...597A..17H).
* g = B30-LB28 contains a APEX/SABOCA at 350um
* (B30-SB16, Huelamo et al. 2017A&A...597A..17H).
* h = IRAS 05293+1207, RA=83.0450, Dec=+12.1629, [12]<0.250Jy, [25]< 0.271Jy,
* [60]=1.15+/-0.1725Jy, [100]< 26.6Jy.
* i = B30-LB33 and B30-LB34 very close to each other, see Sect. A.1.
* j = The MIPS M2 flux probably corresponds to the counterpart identified
* as LB25c.
*Note (2): The coordinates have an error of approximately 1 arcsec
* (pointing rms, but each APEX/LABOCA cell has a 4 arcsec size).
*Note (3): The absolute uncertainty in the flux scale is 8%.
*Note (4): For point-like sources, the value is equal to S870.
*Note (5): derived using the 870um flux densities and assuming a dust temperature
* of 15K (Sanchez-Monge et al. 2013MNRAS.432.3288S), and a dust (and gas) mass
* opacity coefficient of 0.0175cm^2^/g (obtained by interpolating the tabulated
* values of Ossenkopf & Henning (1994, Cat. J/A+A/291/943), see Sect. 3.1).
* The uncertainty in the masses due to the dust temperature and opacity law is
* estimated to be a factor of 4.
*Note (6): Groups as follows:
* A1 = detection at 70 and 24um
* A2 = detection at 70 and upper limit at 24um
* A2* = detection at 70 and upper limit at 24um, extended emission at 24um
* B1 = upper limit at 70um and detection at 24um
* B2 = upper limit both at 70 and 24um
* B2* = upper limit both at 70 and 24um, extended emission at 24um
* C1 = no data at 70um and detection at 24um
* C2 = no data at 70um and upper limit at 24um
* C2* = no data at 70um and upper limit at 24um, extended emission at 24um
*Note (7): Tentative classification as follows:
* YSO = young stellar object, with a optical and/or near-IR sources within
* 5 arcsec of APEX/LABOCA central coordinate. Our tentative
* interpretation is that most of them are proto-stars or proto-BDs.
* Excess = optical and/or near-IR sources with excesses within the APEX/LABOCA
* beam but farther than 5 arcsec. The submm source cannot be assigned
* unambiguously to any counterpart so its nature remains unknown.
* Starless = possible starless core, since there is neither counterparts closer
* than 5 arcsec nor an optical/IR source farther away with excess.
*Note (8): the letter for the identification and the BD for the substellar
* candidates.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Photometric data for all stellar counterparts
integer*4 nr__1
parameter (nr__1=224) ! Number of records
character*151 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*6 B30ID (nr__1) ! B30 source ID (LBNN A)
real*8 RAdeg_1 (nr__1) ! (deg) Source right ascension (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Source declination (J2000)
real*4 Jmag (nr__1) ! (mag) ?=- J band photometry from CAHA/OMEGA2000
real*4 e_Jmag (nr__1) ! (mag) ?=- Uncertainty in J band photometry
* from CAHA/OMEGA2000
real*4 Hmag (nr__1) ! (mag) ?=- H band photometry from CAHA/OMEGA2000
real*4 e_Hmag (nr__1) ! (mag) ?=- Uncertainty in H band photometry
* from CAHA/OMEGA2000
real*4 Kmag (nr__1) ! (mag) ?=- K band photometry from CAHA/OMEGA2000
real*4 e_Kmag (nr__1) ! (mag) ?=- Uncertainty in K band photometry
* from CAHA/OMEGA2000
real*4 I1mag (nr__1) ! (mag) ?=- IRAC band 1 photometry
real*4 e_I1mag (nr__1) ! (mag) ?=- Uncertainty in IRAC band 1 photometry
real*4 I2mag (nr__1) ! (mag) ?=- IRAC band 2 photometry
real*4 e_I2mag (nr__1) ! (mag) ?=- Uncertainty in IRAC band 2 photometry
real*4 I3mag (nr__1) ! (mag) ?=- IRAC band 3 photometry
real*4 e_I3mag (nr__1) ! (mag) ?=- Uncertainty in IRAC band 3 photometry
real*4 I4mag (nr__1) ! (mag) ?=- IRAC band 4 photometry
real*4 e_I4mag (nr__1) ! (mag) ?=- Uncertainty in IRAC band 4 photometry
real*4 M1mag (nr__1) ! (mag) ?=- MIPS band 1 photometry
real*4 e_M1mag (nr__1) ! (mag) ?=- Uncertainty in MIPS band 1 photometry
character*1 Note (nr__1) ! [*] Note for individual entries (1)
*Note (1): Notes as follows:
* * = MIPS photometry might be affected by the inhomogeneous extended emission
* by the nebulosity
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Photometric data, including WISE, release 2013
integer*4 nr__2
parameter (nr__2=67) ! Number of records
character*139 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*6 B30ID_1 (nr__2) ! B30 source ID (LBNN A)
real*8 RAdeg_2 (nr__2) ! (deg) Source right ascension (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Source declination (J2000)
real*4 Dist (nr__2) ! (arcsec) Distance to optical counterpart
character*19 WISE (nr__2) ! Source ID in the WISE catalogue
real*4 W1mag (nr__2) ! (mag) WISE band 1 photometry
character*1 n_W1mag (nr__2) ! [)] ) for values that have not been used
real*4 e_W1mag (nr__2) ! (mag) Uncertainty in WISE band 1 photometry
character*1 n_e_W1mag (nr__2) ! [)] ) for values that have not been used
character*1 l_W2mag (nr__2) ! Limit flag on W2mag
real*4 W2mag (nr__2) ! (mag) WISE band 2 photometry
character*1 n_W2mag (nr__2) ! [)] ) for values that have not been used
real*4 e_W2mag (nr__2) ! (mag) ?=- Uncertainty in WISE band 2 photometry
character*1 n_e_W2mag (nr__2) ! [)] ) for values that have not been used
character*1 l_W3mag (nr__2) ! Limit flag on W3mag
real*4 W3mag (nr__2) ! (mag) WISE band 3 photometry
character*1 n_W3mag (nr__2) ! [)] ) for values that have not been used
real*4 e_W3mag (nr__2) ! (mag) ?=- Uncertainty in WISE band 3 photometry
character*1 n_e_W3mag (nr__2) ! [)] ) for values that have not been used
character*1 l_W4mag (nr__2) ! Limit flag on W4mag
real*4 W4mag (nr__2) ! (mag) WISE band 4 photometry
character*1 n_W4mag (nr__2) ! [)] ) for values that have not been used
real*4 e_W4mag (nr__2) ! (mag) ?=- Uncertainty in WISE band 4 photometry
character*1 n_e_W4mag (nr__2) ! [)] ) for values that have not been used
character*4 Cont (nr__2) ! Contamination and confusion flag,
* one per band (1)
character*1 Ext (nr__2) ! [0/1] Extended source flag (2)
character*4 Qflag (nr__2) ! Photometric quality flag (3)
*Note (1): Contamination and confusion flag as follows:
* "0" for OK, any other value indicate a problem.
* See source catalog for details (Cutri et al., 2013, Cat. II/328).
*Note (2): Extended source flag as follows:
* "0" for OK, any other value indicate a problem.
* See source catalog for details (Cutri et al., 2013, Cat. II/328).
*Note (3): Photometric quality flag as follows:
* "A", "B", "C" for SNR >10, in the range 3-10 and <3
* "U" stand for upper limit
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Optical and near-IR photometry from 2MASS
integer*4 nr__3
parameter (nr__3=68) ! Number of records
character*118 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*6 B30ID_2 (nr__3) ! B30 source ID (LBNN A)
real*8 RAdeg_3 (nr__3) ! (deg) Source right ascension (J2000)
real*8 DEdeg_3 (nr__3) ! (deg) Source declination (J2000)
real*4 Vmag (nr__3) ! (mag) ?=- V magnitude from Dolan & Mathieu
* (2001, Cat. J/AJ/121/2124)
real*4 Rmag (nr__3) ! (mag) ?=- R magnitude from Dolan & Mathieu
* (2001, Cat. J/AJ/121/2124)
real*4 Imag (nr__3) ! (mag) ?=- I magnitude from Dolan & Mathieu
* (2001, Cat. J/AJ/121/2124)
real*4 rmag_1 (nr__3) ! (mag) ?=- r photometry from INT/WFC, this work
real*4 e_rmag (nr__3) ! (mag) ?=- Uncertainty in r photometry from INT/WFC,
* this work
real*4 imag_1 (nr__3) ! (mag) ?=- i photometry from INT/WFC, this work
real*4 e_imag (nr__3) ! (mag) ?=- Uncertainty in r photometry from INT/WFC,
* this work
real*4 J2mag (nr__3) ! (mag) ?=- J photometry from 2MASS
real*4 e_J2mag (nr__3) ! (mag) ?=- Uncertainty in J photometry from 2MASS
real*4 H2mag (nr__3) ! (mag) ?=- H photometry from 2MASS
real*4 e_H2mag (nr__3) ! (mag) ?=- Uncertainty in H photometry from 2MASS
real*4 K2mag (nr__3) ! (mag) ?=- K photometry from 2MASS
real*4 e_K2mag (nr__3) ! (mag) ?=- Uncertainty in K photometry from 2MASS
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! AKARI photometry for APEX/LABOCA sources
integer*4 nr__4
parameter (nr__4=6) ! Number of records
character*132 ar__4 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*6 B30ID_3 (nr__4) ! B30 source ID (LBNN A)
real*8 RAdeg_4 (nr__4) ! (deg) Source right ascension (J2000)
real*8 DEdeg_4 (nr__4) ! (deg) Source declination (J2000)
real*4 S09 (nr__4) ! (Jy) ?=- AKARI flux density at 9 microns
real*4 e_S09 (nr__4) ! (Jy) ?=- Uncertainty in AKARI flux density at 9um
real*4 S18 (nr__4) ! (Jy) ?=- AKARI flux density at 18 microns
real*8 e_S18 (nr__4) ! (Jy) ?=- Uncertainty in AKARI flux density at 18um
real*4 S65 (nr__4) ! (Jy) ?=- AKARI flux density at 65 microns
real*4 e_S65 (nr__4) ! (Jy) ?=- Uncertainty in AKARI flux density at 65um
real*4 S90 (nr__4) ! (Jy) ?=- AKARI flux density at 90 microns
real*4 e_S90 (nr__4) ! (Jy) ?=- Uncertainty in AKARI flux density at 90um
real*4 S140 (nr__4) ! (Jy) ?=- AKARI flux density at 140 microns
real*4 e_S140 (nr__4) ! (Jy) ?=- Uncertainty in AKARI flux density at 140um
real*4 S160 (nr__4) ! (Jy) ?=- AKARI flux density at 160 microns
real*4 e_S160 (nr__4) ! (Jy) ?=- Uncertainty in AKARI flux density at 165um
character*14 AKARI (nr__4) ! AKARI name
character*1 Note_1 (nr__4) ! [ab- ] Note for individual entries (1)
*Note (1): Notes as follows:
* a: IRAS05286+1203, RA=82.8671, DEC=+12.0899,
* [12]=7.79e-01+/-0.0779 Jy, [25]=1.91e+00+/-0.2101 Jy,
* [60]=9.35e+00+/-1.7765 Jy, [100]=7.04e+01+/-14.784 Jy
* b: IRAS05293+1207, RA=83.0450, DEC=+12.1629,
* [12]<2.50e-01 Jy, [25] <2.71e-01 Jy,
* [60]=1.15e+00+/-0.1725 Jy, [100] <2.66e+01 Jy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Properties and memberships criteria for all
counterparts within the APEX/LABOCA beam
integer*4 nr__5
parameter (nr__5=224) ! Number of records
character*128 ar__5 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*6 B30ID_4 (nr__5) ! B30 source ID (LBNN A)
real*8 RAdeg_5 (nr__5) ! (deg) Source right ascension (J2000)
real*8 DEdeg_5 (nr__5) ! (deg) Source declination (J2000)
integer*4 Teff (nr__5) ! (K) ?=- Source effective temperature
* assuming Av=0.322
real*4 Lbol (nr__5) ! (Lsun) ?=- Source bolometric luminosity
* assuming Av=0.322
real*4 e_Lbol (nr__5) ! (Lsun) ?=- Uncertainty in source bolometric
* luminosity assuming Av=0.322
real*4 Av_1 (nr__5) ! (mag) ?=- Visual extinction
integer*4 Teffvar (nr__5) ! (K) ?=- Source effective temperature with
* variable Av
real*4 Lbolvar (nr__5) ! (Lsun) ?=- Source bolometric luminosity with
* variable Av
real*4 e_Lbolvar (nr__5) ! (Lsun) ?=- Uncertainty in source bolometric
* luminosity with variable Av
character*16 Memb1 (nr__5) ! [YN?- ] Membership criteria based on the
* location on Herzprung-Russell, optical-IR
* Color-Color and Color-Magnitude Diagrams (1)
character*2 MembFinal (nr__5) ! Final membership assessment
integer*4 NpFit (nr__5) ! ?=- Number of photometric points
* on the SED: Fit
integer*4 Nptot (nr__5) ! ?=- Number of photometric points
* on the SED: total
character*4 Class (nr__5) ! Evolutionary class based on a Spitzer/IRAC
* (I1-I2, I3-I4) CCD
character*3 IRACSlope (nr__5) ! IRAC Slope
character*11 IRExcess (nr__5) ! [YNTG?- ] IR excess based on optical,
* near-IR and Spitzer Color-Color and
* Color-Magnitude Diagrams (2)
*Note (1): Membership criteria A?B?C?D?E?F?G?H? for
* HRD(AV=0.322), HRD(AV variable), (J, J-I1), (J, J-I2, (H, H-I1), (H, H-I2),
* (K, K-I1), (K, K-I2).
*Note (2): IR excess A?B?C?D?E?F? for
* (I1, I1-I4), (H-K, K-M1), (I1-I4, I4-M1), (I2-I3, I3-I4), (I1-I3, I2-I4),
* (I2, I2-I4).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Temperatures and luminosities for Barnard 30
objects with significant infrared excesses
integer*4 nr__6
parameter (nr__6=23) ! Number of records
character*133 ar__6 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*6 B30ID_5 (nr__6) ! Name (LBNN A)
integer*4 n_B30ID (nr__6) ! [5/8]? Note on B30ID (1)
real*8 RAdeg_6 (nr__6) ! (deg) Right ascension (J2000)
real*8 DEdeg_6 (nr__6) ! (deg) Declination (J2000)
integer*4 Tefff (nr__6) ! (K) ?=- Effective temperature, AV=0.322
real*4 Lbolf (nr__6) ! (Lsun) ?=- Bolometric luminosity, AV=0.322
character*10 Group (nr__6) ! Grouping based on the detection at 70 and 24um
* and presence of counterparts, as discussed
* in Sect. 2.7 (Gr and Type) (2)
integer*4 Teffv (nr__6) ! (K) ?=- Effective temperature, AV=variable
real*4 Lbolv (nr__6) ! (Lsun) ?=- Bolometric luminosity, AV=variable
character*3 MemF (nr__6) ! Membership status (3)
integer*4 Tbol (nr__6) ! (K) ?=- Bolometric temperature
character*2 n_Tbol (nr__6) ! [)* ] Note on Tbol (4)
integer*4 Tbol2 (nr__6) ! (K) ?=- Bolometric temperature, including the flux
* at 870um from APEX/LABOCA and, if present,
* MIPS M2
character*2 n_Tbol2 (nr__6) ! [)* ] Note on Tbol2 (4)
real*4 Lbol_1 (nr__6) ! (Lsun) ?=- Bolometric luminosity
character*2 n_Lbol (nr__6) ! [)* ]Note on Lbol (4)
character*2 l_Lbol2 (nr__6) ! [<= ] Limit flag on Lbol2
real*4 Lbol2 (nr__6) ! (Lsun) ?=- Bolometric luminosity, including the flux
* at 870um from APEX/LABOCA and, if present,
* MIPS M2
character*2 n_Lbol2 (nr__6) ! [)* ] Note on Lbol2 (4)
character*2 l_Lint_1 (nr__6) ! [<= ] Limit flag on Lint
real*4 Lint_1 (nr__6) ! (Lsun) ?=- Internal luminosity, based on MIPS flux at
* 70um, after Dunham et al.
* (2008, Cat. J/ApJS/179/249)
character*6 Class_1 (nr__6) ! SED classification based on the Tbol
* versus Lbol diagram
character*12 Type_1 (nr__6) ! Type
integer*4 Npoint (nr__6) ! Number of points
*Note (1): Notes as follows:
* 5 = There is a submm source nearby, namely SB08, detected with APEX/SABOCA
* at 350um (Huelamo et al. 2017A&A...597A..17H)
* 6 = There is a submm source nearby, namely SB05, detected with APEX/SABOCA
* at 350um (Huelamo et al. 2017A&A...597A..17H)
* 7 = There is a submm source nearby, namely SB04, detected with APEX/SABOCA
* at 350um (Huelamo et al. 2017A&A...597A..17H)
* 8 = There is a submm source nearby, namely SB03, detected with APEX/SABOCA
* at 350um (Huelamo et al. 2017A&A...597A..17H)
*Note (2): Codes as follows:
* A1 = Detection at 70 and 24um
* A2 = detection at 70 and upper limit at 24um
* A2* = detection at 70 and upper limit at 24um, extended emission at 24um
* B1 = upper limit at 70um and detection at 24um
* C1 = no data at 70um and detection at 24um
* YSO = Young Stellar Object, with a optical and/or nearIR sources within
* 5 arcsec of APEX/LABOCA central coordinate. Our tentative
* interpretation is that most of them are proto-stars or proto-BDs.
* excess = Optical and/or nearIR sources with excesses within the APEX/LABOCA
* beam but farther than 5 arcsec. The submm source cannot be assigned
* unambiguously to any counterpart so its nature remains unknown.
*
* Objects without detected emission at 24 and/or 70um are not listed here
* (B2 and C2).
*Note (3): + for membership status updated from Table 7 to Y? based on all the
* available data and our analysis.
*Note (4): )*: Possible detection at MIPS 24um, but if real, it might be affected
* by the inhomogeneous extended emission by the nebulosity.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'list.dat' ! List of fits images
integer*4 nr__7
parameter (nr__7=12) ! Number of records
character*131 ar__7 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 RAdeg_7 (nr__7) ! (deg) Right Ascension of center (J2000)
real*8 DEdeg_7 (nr__7) ! (deg) Declination of center (J2000)
real*4 Scale (nr__7) ! (arcsec/pix) ? Scale of the image
integer*4 Nx (nr__7) ! Number of pixels along X-axis
integer*4 Ny (nr__7) ! Number of pixels along Y-axis
character*24 Obs_Date (nr__7) ! Observation date (YYYY-MM-DDThh:mm:ss.ssss)
integer*4 size (nr__7) ! (Kibyte) Size of FITS file
character*31 FileName (nr__7) ! Name of FITS file, in subdirectory fits
character*32 Title (nr__7) ! Title of the FITS file
C=============================================================================
C Loading file 'table2.dat' ! APEX/LABOCA sources in Barnard 30
C Format for file interpretation
1 format(
+ A4,A1,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,I3,1X,I3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,A1,F5.3,A1,A1,F4.2,A1,I2,A2,1X,I2,
+ 1X,A3,1X,A8,1X,A1,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,34
read(1,'(A98)')ar__
read(ar__,1)
+ Name(i__),n_Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),S870(i__),Ipeak(i__),AV(i__),
+ Mall(i__),Mpeak(i__),l_Lint(i__),Lint(i__),n_Lint(i__),
+ l_F70(i__),F70(i__),n_F70(i__),AngSize(i__),n_AngSize(i__),
+ rms(i__),Gr(i__),Type(i__),BD_T(i__),BD_5(i__)
if (idig(ar__(39:43)).EQ.0) AV(i__) = rNULL__
if (idig(ar__(58:62)).EQ.0) Lint(i__) = rNULL__
if (idig(ar__(65:68)).EQ.0) F70(i__) = rNULL__
if(ar__(70:71) .EQ. '') AngSize(i__) = iNULL__
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,1)
+ Name(i__),n_Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),S870(i__),Ipeak(i__),AV(i__),
+ Mall(i__),Mpeak(i__),l_Lint(i__),Lint(i__),n_Lint(i__),
+ l_F70(i__),F70(i__),n_F70(i__),AngSize(i__),n_AngSize(i__),
+ rms(i__),Gr(i__),Type(i__),BD_T(i__),BD_5(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' ! Photometric data for all stellar counterparts
C Format for file interpretation
2 format(
+ A6,3X,F11.8,1X,F10.7,3X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,
+ 1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,
+ F6.3,1X,F5.3,2X,F5.3,5X,F5.3,5X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,224
read(1,'(A151)')ar__1
read(ar__1,2)
+ B30ID(i__),RAdeg_1(i__),DEdeg_1(i__),Jmag(i__),e_Jmag(i__),
+ Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),I1mag(i__),
+ e_I1mag(i__),I2mag(i__),e_I2mag(i__),I3mag(i__),e_I3mag(i__),
+ I4mag(i__),e_I4mag(i__),M1mag(i__),e_M1mag(i__),Note(i__)
if (idig(ar__1(35:40)).EQ.0) Jmag(i__) = rNULL__
if (idig(ar__1(42:46)).EQ.0) e_Jmag(i__) = rNULL__
if (idig(ar__1(48:53)).EQ.0) Hmag(i__) = rNULL__
if (idig(ar__1(55:59)).EQ.0) e_Hmag(i__) = rNULL__
if (idig(ar__1(61:66)).EQ.0) Kmag(i__) = rNULL__
if (idig(ar__1(68:72)).EQ.0) e_Kmag(i__) = rNULL__
if (idig(ar__1(75:80)).EQ.0) I1mag(i__) = rNULL__
if (idig(ar__1(82:86)).EQ.0) e_I1mag(i__) = rNULL__
if (idig(ar__1(89:94)).EQ.0) I2mag(i__) = rNULL__
if (idig(ar__1(96:100)).EQ.0) e_I2mag(i__) = rNULL__
if (idig(ar__1(103:108)).EQ.0) I3mag(i__) = rNULL__
if (idig(ar__1(110:114)).EQ.0) e_I3mag(i__) = rNULL__
if (idig(ar__1(117:122)).EQ.0) I4mag(i__) = rNULL__
if (idig(ar__1(124:128)).EQ.0) e_I4mag(i__) = rNULL__
if (idig(ar__1(131:135)).EQ.0) M1mag(i__) = rNULL__
if (idig(ar__1(141:145)).EQ.0) e_M1mag(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ B30ID(i__),RAdeg_1(i__),DEdeg_1(i__),Jmag(i__),e_Jmag(i__),
+ Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),I1mag(i__),
+ e_I1mag(i__),I2mag(i__),e_I2mag(i__),I3mag(i__),e_I3mag(i__),
+ I4mag(i__),e_I4mag(i__),M1mag(i__),e_M1mag(i__),Note(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Photometric data, including WISE, release 2013
C Format for file interpretation
3 format(
+ A6,3X,F10.7,2X,F10.7,2X,F5.3,1X,A19,1X,F6.3,A1,2X,F5.3,A1,1X,
+ A1,F6.3,A1,4X,F5.3,A1,1X,A1,F6.3,A1,3X,F5.3,A1,1X,A1,F5.3,A1,
+ 1X,F5.3,A1,2X,A4,1X,A1,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,67
read(1,'(A139)')ar__2
read(ar__2,3)
+ B30ID_1(i__),RAdeg_2(i__),DEdeg_2(i__),Dist(i__),WISE(i__),
+ W1mag(i__),n_W1mag(i__),e_W1mag(i__),n_e_W1mag(i__),
+ l_W2mag(i__),W2mag(i__),n_W2mag(i__),e_W2mag(i__),
+ n_e_W2mag(i__),l_W3mag(i__),W3mag(i__),n_W3mag(i__),
+ e_W3mag(i__),n_e_W3mag(i__),l_W4mag(i__),W4mag(i__),
+ n_W4mag(i__),e_W4mag(i__),n_e_W4mag(i__),Cont(i__),Ext(i__),
+ Qflag(i__)
if (idig(ar__2(88:92)).EQ.0) e_W2mag(i__) = rNULL__
if (idig(ar__2(106:110)).EQ.0) e_W3mag(i__) = rNULL__
if (idig(ar__2(121:125)).EQ.0) e_W4mag(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ B30ID_1(i__),RAdeg_2(i__),DEdeg_2(i__),Dist(i__),WISE(i__),
+ W1mag(i__),n_W1mag(i__),e_W1mag(i__),n_e_W1mag(i__),
+ l_W2mag(i__),W2mag(i__),n_W2mag(i__),e_W2mag(i__),
+ n_e_W2mag(i__),l_W3mag(i__),W3mag(i__),n_W3mag(i__),
+ e_W3mag(i__),n_e_W3mag(i__),l_W4mag(i__),W4mag(i__),
+ n_W4mag(i__),e_W4mag(i__),n_e_W4mag(i__),Cont(i__),Ext(i__),
+ Qflag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Optical and near-IR photometry from 2MASS
C Format for file interpretation
4 format(
+ A6,3X,F10.7,2X,F10.7,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.3,
+ 1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
+ F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,68
read(1,'(A118)')ar__3
read(ar__3,4)
+ B30ID_2(i__),RAdeg_3(i__),DEdeg_3(i__),Vmag(i__),Rmag(i__),
+ Imag(i__),rmag_1(i__),e_rmag(i__),imag_1(i__),e_imag(i__),
+ J2mag(i__),e_J2mag(i__),H2mag(i__),e_H2mag(i__),K2mag(i__),
+ e_K2mag(i__)
if (idig(ar__3(34:39)).EQ.0) Vmag(i__) = rNULL__
if (idig(ar__3(41:46)).EQ.0) Rmag(i__) = rNULL__
if (idig(ar__3(48:53)).EQ.0) Imag(i__) = rNULL__
if (idig(ar__3(55:60)).EQ.0) rmag_1(i__) = rNULL__
if (idig(ar__3(62:66)).EQ.0) e_rmag(i__) = rNULL__
if (idig(ar__3(68:73)).EQ.0) imag_1(i__) = rNULL__
if (idig(ar__3(75:79)).EQ.0) e_imag(i__) = rNULL__
if (idig(ar__3(81:86)).EQ.0) J2mag(i__) = rNULL__
if (idig(ar__3(88:92)).EQ.0) e_J2mag(i__) = rNULL__
if (idig(ar__3(94:99)).EQ.0) H2mag(i__) = rNULL__
if (idig(ar__3(101:105)).EQ.0) e_H2mag(i__) = rNULL__
if (idig(ar__3(107:112)).EQ.0) K2mag(i__) = rNULL__
if (idig(ar__3(114:118)).EQ.0) e_K2mag(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ B30ID_2(i__),RAdeg_3(i__),DEdeg_3(i__),Vmag(i__),Rmag(i__),
+ Imag(i__),rmag_1(i__),e_rmag(i__),imag_1(i__),e_imag(i__),
+ J2mag(i__),e_J2mag(i__),H2mag(i__),e_H2mag(i__),K2mag(i__),
+ e_K2mag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! AKARI photometry for APEX/LABOCA sources
C Format for file interpretation
5 format(
+ A6,1X,F10.7,2X,F10.7,2X,F5.3,2X,F6.4,1X,F5.3,2X,F7.5,1X,F5.3,
+ 2X,F5.3,1X,F6.3,2X,F6.4,1X,F6.3,1X,F5.3,2X,F5.2,2X,F5.2,2X,
+ A14,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,6
read(1,'(A132)')ar__4
read(ar__4,5)
+ B30ID_3(i__),RAdeg_4(i__),DEdeg_4(i__),S09(i__),e_S09(i__),
+ S18(i__),e_S18(i__),S65(i__),e_S65(i__),S90(i__),e_S90(i__),
+ S140(i__),e_S140(i__),S160(i__),e_S160(i__),AKARI(i__),
+ Note_1(i__)
if (idig(ar__4(32:36)).EQ.0) S09(i__) = rNULL__
if (idig(ar__4(39:44)).EQ.0) e_S09(i__) = rNULL__
if (idig(ar__4(46:50)).EQ.0) S18(i__) = rNULL__
if (idig(ar__4(53:59)).EQ.0) e_S18(i__) = rNULL__
if (idig(ar__4(61:65)).EQ.0) S65(i__) = rNULL__
if (idig(ar__4(68:72)).EQ.0) e_S65(i__) = rNULL__
if (idig(ar__4(74:79)).EQ.0) S90(i__) = rNULL__
if (idig(ar__4(82:87)).EQ.0) e_S90(i__) = rNULL__
if (idig(ar__4(89:94)).EQ.0) S140(i__) = rNULL__
if (idig(ar__4(96:100)).EQ.0) e_S140(i__) = rNULL__
if (idig(ar__4(103:107)).EQ.0) S160(i__) = rNULL__
if (idig(ar__4(110:114)).EQ.0) e_S160(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ B30ID_3(i__),RAdeg_4(i__),DEdeg_4(i__),S09(i__),e_S09(i__),
+ S18(i__),e_S18(i__),S65(i__),e_S65(i__),S90(i__),e_S90(i__),
+ S140(i__),e_S140(i__),S160(i__),e_S160(i__),AKARI(i__),
+ Note_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Properties and memberships criteria for all
* counterparts within the APEX/LABOCA beam
C Format for file interpretation
6 format(
+ A6,5X,F11.8,1X,F10.7,2X,I4,1X,E6.3,1X,E6.2,1X,F5.3,1X,I5,1X,
+ E7.4,1X,E6.3,1X,A16,1X,A2,1X,I2,1X,I2,2X,A4,1X,A3,1X,A11)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,224
read(1,'(A128)')ar__5
read(ar__5,6)
+ B30ID_4(i__),RAdeg_5(i__),DEdeg_5(i__),Teff(i__),Lbol(i__),
+ e_Lbol(i__),Av_1(i__),Teffvar(i__),Lbolvar(i__),
+ e_Lbolvar(i__),Memb1(i__),MembFinal(i__),NpFit(i__),
+ Nptot(i__),Class(i__),IRACSlope(i__),IRExcess(i__)
if (Teff(i__) .EQ. 45) Teff(i__) = iNULL__
if (idig(ar__5(41:46)).EQ.0) Lbol(i__) = rNULL__
if (idig(ar__5(48:53)).EQ.0) e_Lbol(i__) = rNULL__
if (idig(ar__5(55:59)).EQ.0) Av_1(i__) = rNULL__
if (Teffvar(i__) .EQ. 45) Teffvar(i__) = iNULL__
if (idig(ar__5(67:73)).EQ.0) Lbolvar(i__) = rNULL__
if (idig(ar__5(75:80)).EQ.0) e_Lbolvar(i__) = rNULL__
if (NpFit(i__) .EQ. 45) NpFit(i__) = iNULL__
if (Nptot(i__) .EQ. 45) Nptot(i__) = iNULL__
c ..............Just test output...........
write(6,6)
+ B30ID_4(i__),RAdeg_5(i__),DEdeg_5(i__),Teff(i__),Lbol(i__),
+ e_Lbol(i__),Av_1(i__),Teffvar(i__),Lbolvar(i__),
+ e_Lbolvar(i__),Memb1(i__),MembFinal(i__),NpFit(i__),
+ Nptot(i__),Class(i__),IRACSlope(i__),IRExcess(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! Temperatures and luminosities for Barnard 30
* objects with significant infrared excesses
C Format for file interpretation
7 format(
+ A6,I1,2X,F10.7,1X,F10.7,1X,I4,1X,F6.4,1X,A10,2X,I5,1X,F6.4,1X,
+ A3,1X,I4,A2,1X,I4,A2,1X,F5.3,A2,A2,F6.3,A2,A2,F5.3,1X,A6,1X,
+ A12,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,23
read(1,'(A133)')ar__6
read(ar__6,7)
+ B30ID_5(i__),n_B30ID(i__),RAdeg_6(i__),DEdeg_6(i__),
+ Tefff(i__),Lbolf(i__),Group(i__),Teffv(i__),Lbolv(i__),
+ MemF(i__),Tbol(i__),n_Tbol(i__),Tbol2(i__),n_Tbol2(i__),
+ Lbol_1(i__),n_Lbol(i__),l_Lbol2(i__),Lbol2(i__),n_Lbol2(i__),
+ l_Lint_1(i__),Lint_1(i__),Class_1(i__),Type_1(i__),Npoint(i__)
if(ar__6(7:7) .EQ. '') n_B30ID(i__) = iNULL__
if (Tefff(i__) .EQ. 45) Tefff(i__) = iNULL__
if (idig(ar__6(37:42)).EQ.0) Lbolf(i__) = rNULL__
if (Teffv(i__) .EQ. 45) Teffv(i__) = iNULL__
if (idig(ar__6(62:67)).EQ.0) Lbolv(i__) = rNULL__
if (Tbol(i__) .EQ. 45) Tbol(i__) = iNULL__
if (Tbol2(i__) .EQ. 45) Tbol2(i__) = iNULL__
if (idig(ar__6(87:91)).EQ.0) Lbol_1(i__) = rNULL__
if (idig(ar__6(96:101)).EQ.0) Lbol2(i__) = rNULL__
if (idig(ar__6(106:110)).EQ.0) Lint_1(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ B30ID_5(i__),n_B30ID(i__),RAdeg_6(i__),DEdeg_6(i__),
+ Tefff(i__),Lbolf(i__),Group(i__),Teffv(i__),Lbolv(i__),
+ MemF(i__),Tbol(i__),n_Tbol(i__),Tbol2(i__),n_Tbol2(i__),
+ Lbol_1(i__),n_Lbol(i__),l_Lbol2(i__),Lbol2(i__),n_Lbol2(i__),
+ l_Lint_1(i__),Lint_1(i__),Class_1(i__),Type_1(i__),Npoint(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'list.dat' ! List of fits images
C Format for file interpretation
8 format(F9.5,F9.5,1X,F5.3,1X,I4,1X,I4,1X,A24,1X,I6,1X,A31,1X,A32)
C Effective file loading
open(unit=1,status='old',file=
+'list.dat')
write(6,*) '....Loading file: list.dat'
do i__=1,12
read(1,'(A131)')ar__7
read(ar__7,8)
+ RAdeg_7(i__),DEdeg_7(i__),Scale(i__),Nx(i__),Ny(i__),
+ Obs_Date(i__),size(i__),FileName(i__),Title(i__)
if(ar__7(20:24) .EQ. '') Scale(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ RAdeg_7(i__),DEdeg_7(i__),Scale(i__),Nx(i__),Ny(i__),
+ Obs_Date(i__),size(i__),FileName(i__),Title(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end