Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/612/A79 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-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 ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*4 Name ! B30 source ID (LBNN)
character*1 n_Name ! [a-i] Note on Name (1)
integer*4 RAh ! (h) Maximum emission of the core
* right ascension (J2000) (2)
integer*4 RAm ! (min) Maximum emission of the core
* right ascension (J2000) (2)
real*4 RAs ! (s) Maximum emission of the core
* right ascension (J2000) (2)
character*1 DE_ ! Maximum emission of the core
* declination sign (J2000) (2)
integer*4 DEd ! (deg) Maximum emission of the core
* declination (J2000) (2)
integer*4 DEm ! (arcmin) Maximum emission of the core
* declination (J2000) (2)
real*4 DEs ! (arcsec) Maximum emission of the core
* declination (J2000) (2)
integer*4 S870 ! (mJy) APEX/LABOCA flux density at 870um computed
* above 3{sigma} emission (3)
integer*4 Ipeak ! (mJy/beam) APEX/LABOCA peak intensity at 870um
* (maximum of the emission) (4)
real*4 AV ! (mag) ?=- Values derived from the extinction map,
* with a resolution of 1.5 arcmin, three times
* larger the LABOCA beam-size
real*4 Mall ! (Msun) Envelope masses (5)
real*4 Mpeak ! (Msun) Envelope masses (5)
character*1 l_Lint ! Limit flag on Lint
real*4 Lint ! (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 ! [j] Note on Lint (1)
character*1 l_F70 ! Limit flag on F70
real*4 F70 ! (Jy) ?=- MIPS flux at 70um
character*1 n_F70 ! [j] Note on F70 (1)
integer*4 AngSize ! (arcsec) ? Deconvolved size after fitting a Gaussian to
* the cores, the averaged of the major and
* minor axis is represented
character*2 n_AngSize ! [PL ] "PL" stands for point-like source
integer*4 rms ! (mJy) rms at 870um
character*3 Gr ! Groups based on the detection at 70 and 24um
* and presence of counterparts, as discussed
* in Sect. 2.7 (6)
character*8 Type ! Tentative classification (7)
character*1 BD_T ! [Y ] Presence of BD candidates, Y indicates
* whether there is a BD candidate within
* APEX/LABOCA beam
character*6 BD_5 ! 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 ! B30 source ID (LBNN A)
real*8 RAdeg_1 ! (deg) Source right ascension (J2000)
real*8 DEdeg_1 ! (deg) Source declination (J2000)
real*4 Jmag ! (mag) ?=- J band photometry from CAHA/OMEGA2000
real*4 e_Jmag ! (mag) ?=- Uncertainty in J band photometry
* from CAHA/OMEGA2000
real*4 Hmag ! (mag) ?=- H band photometry from CAHA/OMEGA2000
real*4 e_Hmag ! (mag) ?=- Uncertainty in H band photometry
* from CAHA/OMEGA2000
real*4 Kmag ! (mag) ?=- K band photometry from CAHA/OMEGA2000
real*4 e_Kmag ! (mag) ?=- Uncertainty in K band photometry
* from CAHA/OMEGA2000
real*4 I1mag ! (mag) ?=- IRAC band 1 photometry
real*4 e_I1mag ! (mag) ?=- Uncertainty in IRAC band 1 photometry
real*4 I2mag ! (mag) ?=- IRAC band 2 photometry
real*4 e_I2mag ! (mag) ?=- Uncertainty in IRAC band 2 photometry
real*4 I3mag ! (mag) ?=- IRAC band 3 photometry
real*4 e_I3mag ! (mag) ?=- Uncertainty in IRAC band 3 photometry
real*4 I4mag ! (mag) ?=- IRAC band 4 photometry
real*4 e_I4mag ! (mag) ?=- Uncertainty in IRAC band 4 photometry
real*4 M1mag ! (mag) ?=- MIPS band 1 photometry
real*4 e_M1mag ! (mag) ?=- Uncertainty in MIPS band 1 photometry
character*1 Note ! [*] 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 ! B30 source ID (LBNN A)
real*8 RAdeg_2 ! (deg) Source right ascension (J2000)
real*8 DEdeg_2 ! (deg) Source declination (J2000)
real*4 Dist ! (arcsec) Distance to optical counterpart
character*19 WISE ! Source ID in the WISE catalogue
real*4 W1mag ! (mag) WISE band 1 photometry
character*1 n_W1mag ! [)] ) for values that have not been used
real*4 e_W1mag ! (mag) Uncertainty in WISE band 1 photometry
character*1 n_e_W1mag ! [)] ) for values that have not been used
character*1 l_W2mag ! Limit flag on W2mag
real*4 W2mag ! (mag) WISE band 2 photometry
character*1 n_W2mag ! [)] ) for values that have not been used
real*4 e_W2mag ! (mag) ?=- Uncertainty in WISE band 2 photometry
character*1 n_e_W2mag ! [)] ) for values that have not been used
character*1 l_W3mag ! Limit flag on W3mag
real*4 W3mag ! (mag) WISE band 3 photometry
character*1 n_W3mag ! [)] ) for values that have not been used
real*4 e_W3mag ! (mag) ?=- Uncertainty in WISE band 3 photometry
character*1 n_e_W3mag ! [)] ) for values that have not been used
character*1 l_W4mag ! Limit flag on W4mag
real*4 W4mag ! (mag) WISE band 4 photometry
character*1 n_W4mag ! [)] ) for values that have not been used
real*4 e_W4mag ! (mag) ?=- Uncertainty in WISE band 4 photometry
character*1 n_e_W4mag ! [)] ) for values that have not been used
character*4 Cont ! Contamination and confusion flag,
* one per band (1)
character*1 Ext ! [0/1] Extended source flag (2)
character*4 Qflag ! 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 ! B30 source ID (LBNN A)
real*8 RAdeg_3 ! (deg) Source right ascension (J2000)
real*8 DEdeg_3 ! (deg) Source declination (J2000)
real*4 Vmag ! (mag) ?=- V magnitude from Dolan & Mathieu
* (2001, Cat. J/AJ/121/2124)
real*4 Rmag ! (mag) ?=- R magnitude from Dolan & Mathieu
* (2001, Cat. J/AJ/121/2124)
real*4 Imag ! (mag) ?=- I magnitude from Dolan & Mathieu
* (2001, Cat. J/AJ/121/2124)
real*4 rmag_1 ! (mag) ?=- r photometry from INT/WFC, this work
real*4 e_rmag ! (mag) ?=- Uncertainty in r photometry from INT/WFC,
* this work
real*4 imag_1 ! (mag) ?=- i photometry from INT/WFC, this work
real*4 e_imag ! (mag) ?=- Uncertainty in r photometry from INT/WFC,
* this work
real*4 J2mag ! (mag) ?=- J photometry from 2MASS
real*4 e_J2mag ! (mag) ?=- Uncertainty in J photometry from 2MASS
real*4 H2mag ! (mag) ?=- H photometry from 2MASS
real*4 e_H2mag ! (mag) ?=- Uncertainty in H photometry from 2MASS
real*4 K2mag ! (mag) ?=- K photometry from 2MASS
real*4 e_K2mag ! (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 ! B30 source ID (LBNN A)
real*8 RAdeg_4 ! (deg) Source right ascension (J2000)
real*8 DEdeg_4 ! (deg) Source declination (J2000)
real*4 S09 ! (Jy) ?=- AKARI flux density at 9 microns
real*4 e_S09 ! (Jy) ?=- Uncertainty in AKARI flux density at 9um
real*4 S18 ! (Jy) ?=- AKARI flux density at 18 microns
real*8 e_S18 ! (Jy) ?=- Uncertainty in AKARI flux density at 18um
real*4 S65 ! (Jy) ?=- AKARI flux density at 65 microns
real*4 e_S65 ! (Jy) ?=- Uncertainty in AKARI flux density at 65um
real*4 S90 ! (Jy) ?=- AKARI flux density at 90 microns
real*4 e_S90 ! (Jy) ?=- Uncertainty in AKARI flux density at 90um
real*4 S140 ! (Jy) ?=- AKARI flux density at 140 microns
real*4 e_S140 ! (Jy) ?=- Uncertainty in AKARI flux density at 140um
real*4 S160 ! (Jy) ?=- AKARI flux density at 160 microns
real*4 e_S160 ! (Jy) ?=- Uncertainty in AKARI flux density at 165um
character*14 AKARI ! AKARI name
character*1 Note_1 ! [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 ! B30 source ID (LBNN A)
real*8 RAdeg_5 ! (deg) Source right ascension (J2000)
real*8 DEdeg_5 ! (deg) Source declination (J2000)
integer*4 Teff ! (K) ?=- Source effective temperature
* assuming Av=0.322
real*4 Lbol ! (Lsun) ?=- Source bolometric luminosity
* assuming Av=0.322
real*4 e_Lbol ! (Lsun) ?=- Uncertainty in source bolometric
* luminosity assuming Av=0.322
real*4 Av_1 ! (mag) ?=- Visual extinction
integer*4 Teffvar ! (K) ?=- Source effective temperature with
* variable Av
real*4 Lbolvar ! (Lsun) ?=- Source bolometric luminosity with
* variable Av
real*4 e_Lbolvar ! (Lsun) ?=- Uncertainty in source bolometric
* luminosity with variable Av
character*16 Memb1 ! [YN?- ] Membership criteria based on the
* location on Herzprung-Russell, optical-IR
* Color-Color and Color-Magnitude Diagrams (1)
character*2 MembFinal ! Final membership assessment
integer*4 NpFit ! ?=- Number of photometric points
* on the SED: Fit
integer*4 Nptot ! ?=- Number of photometric points
* on the SED: total
character*4 Class ! Evolutionary class based on a Spitzer/IRAC
* (I1-I2, I3-I4) CCD
character*3 IRACSlope ! IRAC Slope
character*11 IRExcess ! [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 ! Name (LBNN A)
integer*4 n_B30ID ! [5/8]? Note on B30ID (1)
real*8 RAdeg_6 ! (deg) Right ascension (J2000)
real*8 DEdeg_6 ! (deg) Declination (J2000)
integer*4 Tefff ! (K) ?=- Effective temperature, AV=0.322
real*4 Lbolf ! (Lsun) ?=- Bolometric luminosity, AV=0.322
character*10 Group ! 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 ! (K) ?=- Effective temperature, AV=variable
real*4 Lbolv ! (Lsun) ?=- Bolometric luminosity, AV=variable
character*3 MemF ! Membership status (3)
integer*4 Tbol ! (K) ?=- Bolometric temperature
character*2 n_Tbol ! [)* ] Note on Tbol (4)
integer*4 Tbol2 ! (K) ?=- Bolometric temperature, including the flux
* at 870um from APEX/LABOCA and, if present,
* MIPS M2
character*2 n_Tbol2 ! [)* ] Note on Tbol2 (4)
real*4 Lbol_1 ! (Lsun) ?=- Bolometric luminosity
character*2 n_Lbol ! [)* ]Note on Lbol (4)
character*2 l_Lbol2 ! [<= ] Limit flag on Lbol2
real*4 Lbol2 ! (Lsun) ?=- Bolometric luminosity, including the flux
* at 870um from APEX/LABOCA and, if present,
* MIPS M2
character*2 n_Lbol2 ! [)* ] Note on Lbol2 (4)
character*2 l_Lint_1 ! [<= ] Limit flag on Lint
real*4 Lint_1 ! (Lsun) ?=- Internal luminosity, based on MIPS flux at
* 70um, after Dunham et al.
* (2008, Cat. J/ApJS/179/249)
character*6 Class_1 ! SED classification based on the Tbol
* versus Lbol diagram
character*12 Type_1 ! Type
integer*4 Npoint ! 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 ! (deg) Right Ascension of center (J2000)
real*8 DEdeg_7 ! (deg) Declination of center (J2000)
real*4 Scale ! (arcsec/pix) ? Scale of the image
integer*4 Nx ! Number of pixels along X-axis
integer*4 Ny ! Number of pixels along Y-axis
character*24 Obs_Date ! Observation date (YYYY-MM-DDThh:mm:ss.ssss)
integer*4 size ! (Kibyte) Size of FITS file
character*31 FileName ! Name of FITS file, in subdirectory fits
character*32 Title ! 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,n_Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,S870,Ipeak,AV,Mall,
+ Mpeak,l_Lint,Lint,n_Lint,l_F70,F70,n_F70,AngSize,n_AngSize,
+ rms,Gr,Type,BD_T,BD_5
if (idig(ar__(39:43)).EQ.0) AV = rNULL__
if (idig(ar__(58:62)).EQ.0) Lint = rNULL__
if (idig(ar__(65:68)).EQ.0) F70 = rNULL__
if(ar__(70:71) .EQ. '') AngSize = iNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ Name,n_Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,S870,Ipeak,AV,Mall,
+ Mpeak,l_Lint,Lint,n_Lint,l_F70,F70,n_F70,AngSize,n_AngSize,
+ rms,Gr,Type,BD_T,BD_5
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
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,RAdeg_1,DEdeg_1,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,e_Kmag,
+ I1mag,e_I1mag,I2mag,e_I2mag,I3mag,e_I3mag,I4mag,e_I4mag,M1mag,
+ e_M1mag,Note
if (idig(ar__1(35:40)).EQ.0) Jmag = rNULL__
if (idig(ar__1(42:46)).EQ.0) e_Jmag = rNULL__
if (idig(ar__1(48:53)).EQ.0) Hmag = rNULL__
if (idig(ar__1(55:59)).EQ.0) e_Hmag = rNULL__
if (idig(ar__1(61:66)).EQ.0) Kmag = rNULL__
if (idig(ar__1(68:72)).EQ.0) e_Kmag = rNULL__
if (idig(ar__1(75:80)).EQ.0) I1mag = rNULL__
if (idig(ar__1(82:86)).EQ.0) e_I1mag = rNULL__
if (idig(ar__1(89:94)).EQ.0) I2mag = rNULL__
if (idig(ar__1(96:100)).EQ.0) e_I2mag = rNULL__
if (idig(ar__1(103:108)).EQ.0) I3mag = rNULL__
if (idig(ar__1(110:114)).EQ.0) e_I3mag = rNULL__
if (idig(ar__1(117:122)).EQ.0) I4mag = rNULL__
if (idig(ar__1(124:128)).EQ.0) e_I4mag = rNULL__
if (idig(ar__1(131:135)).EQ.0) M1mag = rNULL__
if (idig(ar__1(141:145)).EQ.0) e_M1mag = rNULL__
c ..............Just test output...........
write(6,2)
+ B30ID,RAdeg_1,DEdeg_1,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,e_Kmag,
+ I1mag,e_I1mag,I2mag,e_I2mag,I3mag,e_I3mag,I4mag,e_I4mag,M1mag,
+ e_M1mag,Note
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,RAdeg_2,DEdeg_2,Dist,WISE,W1mag,n_W1mag,e_W1mag,
+ n_e_W1mag,l_W2mag,W2mag,n_W2mag,e_W2mag,n_e_W2mag,l_W3mag,
+ W3mag,n_W3mag,e_W3mag,n_e_W3mag,l_W4mag,W4mag,n_W4mag,e_W4mag,
+ n_e_W4mag,Cont,Ext,Qflag
if (idig(ar__2(88:92)).EQ.0) e_W2mag = rNULL__
if (idig(ar__2(106:110)).EQ.0) e_W3mag = rNULL__
if (idig(ar__2(121:125)).EQ.0) e_W4mag = rNULL__
c ..............Just test output...........
write(6,3)
+ B30ID_1,RAdeg_2,DEdeg_2,Dist,WISE,W1mag,n_W1mag,e_W1mag,
+ n_e_W1mag,l_W2mag,W2mag,n_W2mag,e_W2mag,n_e_W2mag,l_W3mag,
+ W3mag,n_W3mag,e_W3mag,n_e_W3mag,l_W4mag,W4mag,n_W4mag,e_W4mag,
+ n_e_W4mag,Cont,Ext,Qflag
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,RAdeg_3,DEdeg_3,Vmag,Rmag,Imag,rmag_1,e_rmag,imag_1,
+ e_imag,J2mag,e_J2mag,H2mag,e_H2mag,K2mag,e_K2mag
if (idig(ar__3(34:39)).EQ.0) Vmag = rNULL__
if (idig(ar__3(41:46)).EQ.0) Rmag = rNULL__
if (idig(ar__3(48:53)).EQ.0) Imag = rNULL__
if (idig(ar__3(55:60)).EQ.0) rmag_1 = rNULL__
if (idig(ar__3(62:66)).EQ.0) e_rmag = rNULL__
if (idig(ar__3(68:73)).EQ.0) imag_1 = rNULL__
if (idig(ar__3(75:79)).EQ.0) e_imag = rNULL__
if (idig(ar__3(81:86)).EQ.0) J2mag = rNULL__
if (idig(ar__3(88:92)).EQ.0) e_J2mag = rNULL__
if (idig(ar__3(94:99)).EQ.0) H2mag = rNULL__
if (idig(ar__3(101:105)).EQ.0) e_H2mag = rNULL__
if (idig(ar__3(107:112)).EQ.0) K2mag = rNULL__
if (idig(ar__3(114:118)).EQ.0) e_K2mag = rNULL__
c ..............Just test output...........
write(6,4)
+ B30ID_2,RAdeg_3,DEdeg_3,Vmag,Rmag,Imag,rmag_1,e_rmag,imag_1,
+ e_imag,J2mag,e_J2mag,H2mag,e_H2mag,K2mag,e_K2mag
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,RAdeg_4,DEdeg_4,S09,e_S09,S18,e_S18,S65,e_S65,S90,
+ e_S90,S140,e_S140,S160,e_S160,AKARI,Note_1
if (idig(ar__4(32:36)).EQ.0) S09 = rNULL__
if (idig(ar__4(39:44)).EQ.0) e_S09 = rNULL__
if (idig(ar__4(46:50)).EQ.0) S18 = rNULL__
if (idig(ar__4(53:59)).EQ.0) e_S18 = rNULL__
if (idig(ar__4(61:65)).EQ.0) S65 = rNULL__
if (idig(ar__4(68:72)).EQ.0) e_S65 = rNULL__
if (idig(ar__4(74:79)).EQ.0) S90 = rNULL__
if (idig(ar__4(82:87)).EQ.0) e_S90 = rNULL__
if (idig(ar__4(89:94)).EQ.0) S140 = rNULL__
if (idig(ar__4(96:100)).EQ.0) e_S140 = rNULL__
if (idig(ar__4(103:107)).EQ.0) S160 = rNULL__
if (idig(ar__4(110:114)).EQ.0) e_S160 = rNULL__
c ..............Just test output...........
write(6,5)
+ B30ID_3,RAdeg_4,DEdeg_4,S09,e_S09,S18,e_S18,S65,e_S65,S90,
+ e_S90,S140,e_S140,S160,e_S160,AKARI,Note_1
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,RAdeg_5,DEdeg_5,Teff,Lbol,e_Lbol,Av_1,Teffvar,Lbolvar,
+ e_Lbolvar,Memb1,MembFinal,NpFit,Nptot,Class,IRACSlope,IRExcess
if (Teff .EQ. 45) Teff = iNULL__
if (idig(ar__5(41:46)).EQ.0) Lbol = rNULL__
if (idig(ar__5(48:53)).EQ.0) e_Lbol = rNULL__
if (idig(ar__5(55:59)).EQ.0) Av_1 = rNULL__
if (Teffvar .EQ. 45) Teffvar = iNULL__
if (idig(ar__5(67:73)).EQ.0) Lbolvar = rNULL__
if (idig(ar__5(75:80)).EQ.0) e_Lbolvar = rNULL__
if (NpFit .EQ. 45) NpFit = iNULL__
if (Nptot .EQ. 45) Nptot = iNULL__
c ..............Just test output...........
write(6,6)
+ B30ID_4,RAdeg_5,DEdeg_5,Teff,Lbol,e_Lbol,Av_1,Teffvar,Lbolvar,
+ e_Lbolvar,Memb1,MembFinal,NpFit,Nptot,Class,IRACSlope,IRExcess
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,n_B30ID,RAdeg_6,DEdeg_6,Tefff,Lbolf,Group,Teffv,Lbolv,
+ MemF,Tbol,n_Tbol,Tbol2,n_Tbol2,Lbol_1,n_Lbol,l_Lbol2,Lbol2,
+ n_Lbol2,l_Lint_1,Lint_1,Class_1,Type_1,Npoint
if(ar__6(7:7) .EQ. '') n_B30ID = iNULL__
if (Tefff .EQ. 45) Tefff = iNULL__
if (idig(ar__6(37:42)).EQ.0) Lbolf = rNULL__
if (Teffv .EQ. 45) Teffv = iNULL__
if (idig(ar__6(62:67)).EQ.0) Lbolv = rNULL__
if (Tbol .EQ. 45) Tbol = iNULL__
if (Tbol2 .EQ. 45) Tbol2 = iNULL__
if (idig(ar__6(87:91)).EQ.0) Lbol_1 = rNULL__
if (idig(ar__6(96:101)).EQ.0) Lbol2 = rNULL__
if (idig(ar__6(106:110)).EQ.0) Lint_1 = rNULL__
c ..............Just test output...........
write(6,7)
+ B30ID_5,n_B30ID,RAdeg_6,DEdeg_6,Tefff,Lbolf,Group,Teffv,Lbolv,
+ MemF,Tbol,n_Tbol,Tbol2,n_Tbol2,Lbol_1,n_Lbol,l_Lbol2,Lbol2,
+ n_Lbol2,l_Lint_1,Lint_1,Class_1,Type_1,Npoint
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,DEdeg_7,Scale,Nx,Ny,Obs_Date,size,FileName,Title
if(ar__7(20:24) .EQ. '') Scale = rNULL__
c ..............Just test output...........
write(6,8)
+ RAdeg_7,DEdeg_7,Scale,Nx,Ny,Obs_Date,size,FileName,Title
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