Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/530/A90 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 2024-Apr-20
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/530/A90 Optically bright post-AGB population of LMC (van Aarle+, 2011)
*================================================================================
*The optically bright post-AGB population of the LMC.
* van Aarle E., Van Winckel H., Lloyd Evans T., Ueta T., Wood P.R.,
* Ginsburg A.G.
* <Astron. Astrophys. 530, A90 (2011)>
* =2011A&A...530A..90V
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat_pagb.dat' ! Post-AGB candidates
integer*4 nr__
parameter (nr__=1407) ! Number of records
character*606 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*19 Name (nr__) ! Name (IRAC, JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
character*13 Type (nr__) ! SED type (cat_pagb and cat_spec only)
real*8 Umagz (nr__) ! (mag) ? Apparent U magnitude from Zaritsky et al.
* (2004AJ....128.1606Z )
real*8 e_Umagz (nr__) ! (mag) ? Mean error on Umagz
real*8 Umagm (nr__) ! (mag) ? Apparent U magnitude from Massey
* (2002, Cat II/236)
real*8 e_Umagm (nr__) ! (mag) ? Mean error on Umagm
real*8 Bmagz (nr__) ! (mag) ? Apparent B magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Bmagz (nr__) ! (mag) ? Mean error on Bmagz
real*8 Bmagm (nr__) ! (mag) ? Apparent B magnitude from Massey
* (2002, Cat II/236)
real*8 e_Bmagm (nr__) ! (mag) ? Mean error on Bmagm
real*8 Bmagg (nr__) ! (mag) ? Apparent B magnitude from GSC 2.3.2 (I/305)
real*8 e_Bmagg (nr__) ! (mag) ? Mean error on Bmagg
real*8 Vmagz (nr__) ! (mag) ? Apparent V magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Vmagz (nr__) ! (mag) ? Mean error on Vmagz
real*8 Vmagm (nr__) ! (mag) ? Apparent V magnitude from Massey
* (2002, Cat II/236)
real*8 e_Vmagm (nr__) ! (mag) ? Mean error on Vmagm
real*8 Vmagg (nr__) ! (mag) ? Apparent V magnitude from GSC 2.3.2 (I/305)
real*8 e_Vmagg (nr__) ! (mag) ? Mean error on Vmagg
real*8 Rmagm (nr__) ! (mag) ? Apparent R magnitude from Massey
* (2002, Cat II/236)
real*8 e_Rmagm (nr__) ! (mag) ? Mean error on Rmagm
real*8 Rfmagg (nr__) ! (mag) ? Apparent Rf magnitude from GSC 2.3.2 (I/305)
real*8 e_Rfmagg (nr__) ! (mag) ? Mean error on Rfmagg
real*8 Imagz (nr__) ! (mag) ? Apparent I magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Imagz (nr__) ! (mag) ? Mean error on Imagz
real*8 Jmag (nr__) ! (mag) ? 2MASS J apparent magnitude (II/246)
real*8 e_Jmag (nr__) ! (mag) ? Mean error on Jmag
real*8 Hmag (nr__) ! (mag) ? 2MASS H apparent magnitude (II/246)
real*8 e_Hmag (nr__) ! (mag) ? Mean error on Hmag
real*8 Kmag (nr__) ! (mag) ? 2MASS Ks apparent magnitude (II/246)
real*8 e_Kmag (nr__) ! (mag) ? Mean error on Kmag
real*8 F3_6 (nr__) ! (mJy) ? Spitzer/IRAC 3.6 micron flux
real*8 e_F3_6 (nr__) ! (mJy) ? Mean error on F3.6
real*8 F4_5 (nr__) ! (mJy) ? Spitzer/IRAC 4.5 micron flux
real*8 e_F4_5 (nr__) ! (mJy) ? Mean error on F4.5
real*8 F5_8 (nr__) ! (mJy) ? Spitzer/IRAC 5.8 micron flux
real*8 e_F5_8 (nr__) ! (mJy) ? Mean error on F5.8
real*8 F8_0 (nr__) ! (mJy) ? Spitzer/IRAC 8.0 micron flux
real*8 e_F8_0 (nr__) ! (mJy) ? Mean error on F8.0
real*8 F24 (nr__) ! (mJy) ? Spitzer/MIPS 24 micron flux
real*8 e_F24 (nr__) ! (mJy) ? Mean error on F24
real*8 F70 (nr__) ! (mJy) ? Spitzer/MIPS 70 micron flux
real*8 e_F70 (nr__) ! (mJy) ? Mean error on F70
real*8 F160 (nr__) ! (mJy) ? Spitzer/MIPS 160 micron flux
real*8 e_F160 (nr__) ! (mJy) ? Mean error on F160
integer*4 LumBB (nr__) ! (solLum) ? Black-body-based luminosity
* (cat_pag and cat_spec) or
* Estimated luminosity (cat_sg and cat_yso)
integer*4 LumSED (nr__) ! (solLum) ? SED-based luminosity (1)
integer*4 e_LumSED (nr__) ! (solLum) ? Mean error on LumSED (1)
integer*4 TeffSED (nr__) ! (K) ? SED-based effective temperature (1)
integer*4 e_TeffSED (nr__) ! (K) ? Mean error on TeffSED (1)
real*8 E_B_V (nr__) ! (mag) ? SED-based E(B-V) (1)
real*8 e_E_B_V (nr__) ! (mag) ? Mean error on E(B-V) (1)
character*7 Xcat (nr__) ! ? Cross-correlation with catalogues
* of other types of objects (1) (3)
character*30 SpType (nr__) ! ? MK spectral type (1) (4)
character*30 Var (nr__) ! ? Variability remarks (1) (5)
character*1 n_Name (nr__) ! Reason why this object was discarded
* from our post-AGB sample (2)
*Note (1): For cat_pagb and cat_spec only.
*Note (2): Reason codes as follows (for cat_spec only):
* A = Not in the current version of the SAGE catalogue
* B = Does not obey the colour criterion
* C = Does not survive the black-body luminosity-cut (Sect. 2.3)
* D = Is listed as a cool carbon star by Kontizas et al.
* (2001, Cat. J/A+A/369/932) or an AGB star with spectral type by
* Trams et al. (1999A&A...346..843T)
* E = Is not a post-AGB star according to the spectral type
* F = Does not survive the SED luminosity-cut (Sect. 4.1)
*Note (3): Empirically confirmed cross-matches of the object in the catalogues of
* Sect. 5.1.3 with
* A = the RVT stars listed in the recent OGLE-III Catalogue of Variable
* Stars (Soszynski et al., 2008, Cat J/AcA/58/293)
* B = the RCrB stars in Soszynski et al. (2009, Cat J/AcA/59/335)
* C = the MSX post-AGB stars of Wood & Cohen (2001ASSL..265...71W)
* D = the YSOs from Seale et al. (2009, Cat. J/AJ/699/150)
* E = the AGB stars from Trams et al. (1999A&A...346..843T)
* F = the PNe in Reid & Parker (2006, Cat. J/MNRAS/373/521)
* G = the galaxies in the line of sight of the LMC from NED
*Note (4): All spectral types indicated with (A) were found in the catalogue of
* massive stars of Bonanos et al. (2009AJ....138.1003B), the others were
* determined based on our low-resolution, optical spectra.
*Note (5): For cat_spec only. The variability remarks are based on both the
* luminosity curves as discussed in Sect. 5.2 and cross-correlation with
* the catalogue of Fraser et al. (2008, Cat. J/AJ/136/1242). The latter
* is indicated with 'Fr_' followed by the correct sequence in subscript.
* Sequences 9, 0 and -99 indicate subsequently that the star is
* identified with the One-Year Artifact which is caused by the annual
* observing schedule of the MACHO project, that the star is outside the
* boundaries of any period-luminosity classification, or that Fraser et
* al. (2008, Cat. J/AJ/136/1242) were unable to classify this object.
* The other abbreviations used to characterise the variability of the
* different objects are Population II Cepheids (Ceph.) which include the
* subclass of RV Tauri stars (RVT), semi-regular variables (SRV) with
* the subclass of sequence-D variables (Seq. D) that were recognized
* by us, objects with lightcurves in Fig. 18 that show long-term, slow
* variations (SV), RCrB stars (RCrB), objects in regions of very strong
* nebular emission (em), Ae and Be stars (AB) and objects that do not
* show strong evidence of variability (NV). For all Population II
* Cepheids, the periods we computed between deep minima or two cycles of
* the light curve if the minima are of equal depth, are listed in days.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat_spec.dat' ! Discarded objects with a spectral type
integer*4 nr__1
parameter (nr__1=89) ! Number of records
character*606 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*19 Name_1 (nr__1) ! Name (IRAC, JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg_1 (nr__1) ! (deg) Right ascension (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Declination (J2000)
character*13 Type_1 (nr__1) ! SED type (cat_pagb and cat_spec only)
real*8 Umagz_1 (nr__1) ! (mag) ? Apparent U magnitude from Zaritsky et al.
* (2004AJ....128.1606Z )
real*8 e_Umagz_1 (nr__1) ! (mag) ? Mean error on Umagz
real*8 Umagm_1 (nr__1) ! (mag) ? Apparent U magnitude from Massey
* (2002, Cat II/236)
real*8 e_Umagm_1 (nr__1) ! (mag) ? Mean error on Umagm
real*8 Bmagz_1 (nr__1) ! (mag) ? Apparent B magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Bmagz_1 (nr__1) ! (mag) ? Mean error on Bmagz
real*8 Bmagm_1 (nr__1) ! (mag) ? Apparent B magnitude from Massey
* (2002, Cat II/236)
real*8 e_Bmagm_1 (nr__1) ! (mag) ? Mean error on Bmagm
real*8 Bmagg_1 (nr__1) ! (mag) ? Apparent B magnitude from GSC 2.3.2 (I/305)
real*8 e_Bmagg_1 (nr__1) ! (mag) ? Mean error on Bmagg
real*8 Vmagz_1 (nr__1) ! (mag) ? Apparent V magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Vmagz_1 (nr__1) ! (mag) ? Mean error on Vmagz
real*8 Vmagm_1 (nr__1) ! (mag) ? Apparent V magnitude from Massey
* (2002, Cat II/236)
real*8 e_Vmagm_1 (nr__1) ! (mag) ? Mean error on Vmagm
real*8 Vmagg_1 (nr__1) ! (mag) ? Apparent V magnitude from GSC 2.3.2 (I/305)
real*8 e_Vmagg_1 (nr__1) ! (mag) ? Mean error on Vmagg
real*8 Rmagm_1 (nr__1) ! (mag) ? Apparent R magnitude from Massey
* (2002, Cat II/236)
real*8 e_Rmagm_1 (nr__1) ! (mag) ? Mean error on Rmagm
real*8 Rfmagg_1 (nr__1) ! (mag) ? Apparent Rf magnitude from GSC 2.3.2 (I/305)
real*8 e_Rfmagg_1 (nr__1) ! (mag) ? Mean error on Rfmagg
real*8 Imagz_1 (nr__1) ! (mag) ? Apparent I magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Imagz_1 (nr__1) ! (mag) ? Mean error on Imagz
real*8 Jmag_1 (nr__1) ! (mag) ? 2MASS J apparent magnitude (II/246)
real*8 e_Jmag_1 (nr__1) ! (mag) ? Mean error on Jmag
real*8 Hmag_1 (nr__1) ! (mag) ? 2MASS H apparent magnitude (II/246)
real*8 e_Hmag_1 (nr__1) ! (mag) ? Mean error on Hmag
real*8 Kmag_1 (nr__1) ! (mag) ? 2MASS Ks apparent magnitude (II/246)
real*8 e_Kmag_1 (nr__1) ! (mag) ? Mean error on Kmag
real*8 F3_6_1 (nr__1) ! (mJy) ? Spitzer/IRAC 3.6 micron flux
real*8 e_F3_6_1 (nr__1) ! (mJy) ? Mean error on F3.6
real*8 F4_5_1 (nr__1) ! (mJy) ? Spitzer/IRAC 4.5 micron flux
real*8 e_F4_5_1 (nr__1) ! (mJy) ? Mean error on F4.5
real*8 F5_8_1 (nr__1) ! (mJy) ? Spitzer/IRAC 5.8 micron flux
real*8 e_F5_8_1 (nr__1) ! (mJy) ? Mean error on F5.8
real*8 F8_0_1 (nr__1) ! (mJy) ? Spitzer/IRAC 8.0 micron flux
real*8 e_F8_0_1 (nr__1) ! (mJy) ? Mean error on F8.0
real*8 F24_1 (nr__1) ! (mJy) ? Spitzer/MIPS 24 micron flux
real*8 e_F24_1 (nr__1) ! (mJy) ? Mean error on F24
real*8 F70_1 (nr__1) ! (mJy) ? Spitzer/MIPS 70 micron flux
real*8 e_F70_1 (nr__1) ! (mJy) ? Mean error on F70
real*8 F160_1 (nr__1) ! (mJy) ? Spitzer/MIPS 160 micron flux
real*8 e_F160_1 (nr__1) ! (mJy) ? Mean error on F160
integer*4 LumBB_1 (nr__1) ! (solLum) ? Black-body-based luminosity
* (cat_pag and cat_spec) or
* Estimated luminosity (cat_sg and cat_yso)
integer*4 LumSED_1 (nr__1) ! (solLum) ? SED-based luminosity (1)
integer*4 e_LumSED_1 (nr__1) ! (solLum) ? Mean error on LumSED (1)
integer*4 TeffSED_1 (nr__1) ! (K) ? SED-based effective temperature (1)
integer*4 e_TeffSED_1(nr__1) ! (K) ? Mean error on TeffSED (1)
real*8 E_B_V_1 (nr__1) ! (mag) ? SED-based E(B-V) (1)
real*8 e_E_B_V_1 (nr__1) ! (mag) ? Mean error on E(B-V) (1)
character*7 Xcat_1 (nr__1) ! ? Cross-correlation with catalogues
* of other types of objects (1) (3)
character*30 SpType_1 (nr__1) ! ? MK spectral type (1) (4)
character*30 Var_1 (nr__1) ! ? Variability remarks (1) (5)
character*1 n_Name_1 (nr__1) ! Reason why this object was discarded
* from our post-AGB sample (2)
*Note (1): For cat_pagb and cat_spec only.
*Note (2): Reason codes as follows (for cat_spec only):
* A = Not in the current version of the SAGE catalogue
* B = Does not obey the colour criterion
* C = Does not survive the black-body luminosity-cut (Sect. 2.3)
* D = Is listed as a cool carbon star by Kontizas et al.
* (2001, Cat. J/A+A/369/932) or an AGB star with spectral type by
* Trams et al. (1999A&A...346..843T)
* E = Is not a post-AGB star according to the spectral type
* F = Does not survive the SED luminosity-cut (Sect. 4.1)
*Note (3): Empirically confirmed cross-matches of the object in the catalogues of
* Sect. 5.1.3 with
* A = the RVT stars listed in the recent OGLE-III Catalogue of Variable
* Stars (Soszynski et al., 2008, Cat J/AcA/58/293)
* B = the RCrB stars in Soszynski et al. (2009, Cat J/AcA/59/335)
* C = the MSX post-AGB stars of Wood & Cohen (2001ASSL..265...71W)
* D = the YSOs from Seale et al. (2009, Cat. J/AJ/699/150)
* E = the AGB stars from Trams et al. (1999A&A...346..843T)
* F = the PNe in Reid & Parker (2006, Cat. J/MNRAS/373/521)
* G = the galaxies in the line of sight of the LMC from NED
*Note (4): All spectral types indicated with (A) were found in the catalogue of
* massive stars of Bonanos et al. (2009AJ....138.1003B), the others were
* determined based on our low-resolution, optical spectra.
*Note (5): For cat_spec only. The variability remarks are based on both the
* luminosity curves as discussed in Sect. 5.2 and cross-correlation with
* the catalogue of Fraser et al. (2008, Cat. J/AJ/136/1242). The latter
* is indicated with 'Fr_' followed by the correct sequence in subscript.
* Sequences 9, 0 and -99 indicate subsequently that the star is
* identified with the One-Year Artifact which is caused by the annual
* observing schedule of the MACHO project, that the star is outside the
* boundaries of any period-luminosity classification, or that Fraser et
* al. (2008, Cat. J/AJ/136/1242) were unable to classify this object.
* The other abbreviations used to characterise the variability of the
* different objects are Population II Cepheids (Ceph.) which include the
* subclass of RV Tauri stars (RVT), semi-regular variables (SRV) with
* the subclass of sequence-D variables (Seq. D) that were recognized
* by us, objects with lightcurves in Fig. 18 that show long-term, slow
* variations (SV), RCrB stars (RCrB), objects in regions of very strong
* nebular emission (em), Ae and Be stars (AB) and objects that do not
* show strong evidence of variability (NV). For all Population II
* Cepheids, the periods we computed between deep minima or two cycles of
* the light curve if the minima are of equal depth, are listed in days.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat_yso.dat' ! Discarded possible YSOs
integer*4 nr__2
parameter (nr__2=6823) ! Number of records
character*606 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*19 Name_2 (nr__2) ! Name (IRAC, JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg_2 (nr__2) ! (deg) Right ascension (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Declination (J2000)
character*13 Type_2 (nr__2) ! SED type (cat_pagb and cat_spec only)
real*8 Umagz_2 (nr__2) ! (mag) ? Apparent U magnitude from Zaritsky et al.
* (2004AJ....128.1606Z )
real*8 e_Umagz_2 (nr__2) ! (mag) ? Mean error on Umagz
real*8 Umagm_2 (nr__2) ! (mag) ? Apparent U magnitude from Massey
* (2002, Cat II/236)
real*8 e_Umagm_2 (nr__2) ! (mag) ? Mean error on Umagm
real*8 Bmagz_2 (nr__2) ! (mag) ? Apparent B magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Bmagz_2 (nr__2) ! (mag) ? Mean error on Bmagz
real*8 Bmagm_2 (nr__2) ! (mag) ? Apparent B magnitude from Massey
* (2002, Cat II/236)
real*8 e_Bmagm_2 (nr__2) ! (mag) ? Mean error on Bmagm
real*8 Bmagg_2 (nr__2) ! (mag) ? Apparent B magnitude from GSC 2.3.2 (I/305)
real*8 e_Bmagg_2 (nr__2) ! (mag) ? Mean error on Bmagg
real*8 Vmagz_2 (nr__2) ! (mag) ? Apparent V magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Vmagz_2 (nr__2) ! (mag) ? Mean error on Vmagz
real*8 Vmagm_2 (nr__2) ! (mag) ? Apparent V magnitude from Massey
* (2002, Cat II/236)
real*8 e_Vmagm_2 (nr__2) ! (mag) ? Mean error on Vmagm
real*8 Vmagg_2 (nr__2) ! (mag) ? Apparent V magnitude from GSC 2.3.2 (I/305)
real*8 e_Vmagg_2 (nr__2) ! (mag) ? Mean error on Vmagg
real*8 Rmagm_2 (nr__2) ! (mag) ? Apparent R magnitude from Massey
* (2002, Cat II/236)
real*8 e_Rmagm_2 (nr__2) ! (mag) ? Mean error on Rmagm
real*8 Rfmagg_2 (nr__2) ! (mag) ? Apparent Rf magnitude from GSC 2.3.2 (I/305)
real*8 e_Rfmagg_2 (nr__2) ! (mag) ? Mean error on Rfmagg
real*8 Imagz_2 (nr__2) ! (mag) ? Apparent I magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Imagz_2 (nr__2) ! (mag) ? Mean error on Imagz
real*8 Jmag_2 (nr__2) ! (mag) ? 2MASS J apparent magnitude (II/246)
real*8 e_Jmag_2 (nr__2) ! (mag) ? Mean error on Jmag
real*8 Hmag_2 (nr__2) ! (mag) ? 2MASS H apparent magnitude (II/246)
real*8 e_Hmag_2 (nr__2) ! (mag) ? Mean error on Hmag
real*8 Kmag_2 (nr__2) ! (mag) ? 2MASS Ks apparent magnitude (II/246)
real*8 e_Kmag_2 (nr__2) ! (mag) ? Mean error on Kmag
real*8 F3_6_2 (nr__2) ! (mJy) ? Spitzer/IRAC 3.6 micron flux
real*8 e_F3_6_2 (nr__2) ! (mJy) ? Mean error on F3.6
real*8 F4_5_2 (nr__2) ! (mJy) ? Spitzer/IRAC 4.5 micron flux
real*8 e_F4_5_2 (nr__2) ! (mJy) ? Mean error on F4.5
real*8 F5_8_2 (nr__2) ! (mJy) ? Spitzer/IRAC 5.8 micron flux
real*8 e_F5_8_2 (nr__2) ! (mJy) ? Mean error on F5.8
real*8 F8_0_2 (nr__2) ! (mJy) ? Spitzer/IRAC 8.0 micron flux
real*8 e_F8_0_2 (nr__2) ! (mJy) ? Mean error on F8.0
real*8 F24_2 (nr__2) ! (mJy) ? Spitzer/MIPS 24 micron flux
real*8 e_F24_2 (nr__2) ! (mJy) ? Mean error on F24
real*8 F70_2 (nr__2) ! (mJy) ? Spitzer/MIPS 70 micron flux
real*8 e_F70_2 (nr__2) ! (mJy) ? Mean error on F70
real*8 F160_2 (nr__2) ! (mJy) ? Spitzer/MIPS 160 micron flux
real*8 e_F160_2 (nr__2) ! (mJy) ? Mean error on F160
integer*4 LumBB_2 (nr__2) ! (solLum) ? Black-body-based luminosity
* (cat_pag and cat_spec) or
* Estimated luminosity (cat_sg and cat_yso)
integer*4 LumSED_2 (nr__2) ! (solLum) ? SED-based luminosity (1)
integer*4 e_LumSED_2 (nr__2) ! (solLum) ? Mean error on LumSED (1)
integer*4 TeffSED_2 (nr__2) ! (K) ? SED-based effective temperature (1)
integer*4 e_TeffSED_2(nr__2) ! (K) ? Mean error on TeffSED (1)
real*8 E_B_V_1 (nr__2) ! (mag) ? SED-based E(B-V) (1)
real*8 e_E_B_V_1 (nr__2) ! (mag) ? Mean error on E(B-V) (1)
character*7 Xcat_2 (nr__2) ! ? Cross-correlation with catalogues
* of other types of objects (1) (3)
character*30 SpType_2 (nr__2) ! ? MK spectral type (1) (4)
character*30 Var_2 (nr__2) ! ? Variability remarks (1) (5)
character*1 n_Name_2 (nr__2) ! Reason why this object was discarded
* from our post-AGB sample (2)
*Note (1): For cat_pagb and cat_spec only.
*Note (2): Reason codes as follows (for cat_spec only):
* A = Not in the current version of the SAGE catalogue
* B = Does not obey the colour criterion
* C = Does not survive the black-body luminosity-cut (Sect. 2.3)
* D = Is listed as a cool carbon star by Kontizas et al.
* (2001, Cat. J/A+A/369/932) or an AGB star with spectral type by
* Trams et al. (1999A&A...346..843T)
* E = Is not a post-AGB star according to the spectral type
* F = Does not survive the SED luminosity-cut (Sect. 4.1)
*Note (3): Empirically confirmed cross-matches of the object in the catalogues of
* Sect. 5.1.3 with
* A = the RVT stars listed in the recent OGLE-III Catalogue of Variable
* Stars (Soszynski et al., 2008, Cat J/AcA/58/293)
* B = the RCrB stars in Soszynski et al. (2009, Cat J/AcA/59/335)
* C = the MSX post-AGB stars of Wood & Cohen (2001ASSL..265...71W)
* D = the YSOs from Seale et al. (2009, Cat. J/AJ/699/150)
* E = the AGB stars from Trams et al. (1999A&A...346..843T)
* F = the PNe in Reid & Parker (2006, Cat. J/MNRAS/373/521)
* G = the galaxies in the line of sight of the LMC from NED
*Note (4): All spectral types indicated with (A) were found in the catalogue of
* massive stars of Bonanos et al. (2009AJ....138.1003B), the others were
* determined based on our low-resolution, optical spectra.
*Note (5): For cat_spec only. The variability remarks are based on both the
* luminosity curves as discussed in Sect. 5.2 and cross-correlation with
* the catalogue of Fraser et al. (2008, Cat. J/AJ/136/1242). The latter
* is indicated with 'Fr_' followed by the correct sequence in subscript.
* Sequences 9, 0 and -99 indicate subsequently that the star is
* identified with the One-Year Artifact which is caused by the annual
* observing schedule of the MACHO project, that the star is outside the
* boundaries of any period-luminosity classification, or that Fraser et
* al. (2008, Cat. J/AJ/136/1242) were unable to classify this object.
* The other abbreviations used to characterise the variability of the
* different objects are Population II Cepheids (Ceph.) which include the
* subclass of RV Tauri stars (RVT), semi-regular variables (SRV) with
* the subclass of sequence-D variables (Seq. D) that were recognized
* by us, objects with lightcurves in Fig. 18 that show long-term, slow
* variations (SV), RCrB stars (RCrB), objects in regions of very strong
* nebular emission (em), Ae and Be stars (AB) and objects that do not
* show strong evidence of variability (NV). For all Population II
* Cepheids, the periods we computed between deep minima or two cycles of
* the light curve if the minima are of equal depth, are listed in days.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat_sg.dat' ! Discarded possible supergiants
integer*4 nr__3
parameter (nr__3=286) ! Number of records
character*606 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*19 Name_3 (nr__3) ! Name (IRAC, JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg_3 (nr__3) ! (deg) Right ascension (J2000)
real*8 DEdeg_3 (nr__3) ! (deg) Declination (J2000)
character*13 Type_3 (nr__3) ! SED type (cat_pagb and cat_spec only)
real*8 Umagz_3 (nr__3) ! (mag) ? Apparent U magnitude from Zaritsky et al.
* (2004AJ....128.1606Z )
real*8 e_Umagz_3 (nr__3) ! (mag) ? Mean error on Umagz
real*8 Umagm_3 (nr__3) ! (mag) ? Apparent U magnitude from Massey
* (2002, Cat II/236)
real*8 e_Umagm_3 (nr__3) ! (mag) ? Mean error on Umagm
real*8 Bmagz_3 (nr__3) ! (mag) ? Apparent B magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Bmagz_3 (nr__3) ! (mag) ? Mean error on Bmagz
real*8 Bmagm_3 (nr__3) ! (mag) ? Apparent B magnitude from Massey
* (2002, Cat II/236)
real*8 e_Bmagm_3 (nr__3) ! (mag) ? Mean error on Bmagm
real*8 Bmagg_3 (nr__3) ! (mag) ? Apparent B magnitude from GSC 2.3.2 (I/305)
real*8 e_Bmagg_3 (nr__3) ! (mag) ? Mean error on Bmagg
real*8 Vmagz_3 (nr__3) ! (mag) ? Apparent V magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Vmagz_3 (nr__3) ! (mag) ? Mean error on Vmagz
real*8 Vmagm_3 (nr__3) ! (mag) ? Apparent V magnitude from Massey
* (2002, Cat II/236)
real*8 e_Vmagm_3 (nr__3) ! (mag) ? Mean error on Vmagm
real*8 Vmagg_3 (nr__3) ! (mag) ? Apparent V magnitude from GSC 2.3.2 (I/305)
real*8 e_Vmagg_3 (nr__3) ! (mag) ? Mean error on Vmagg
real*8 Rmagm_3 (nr__3) ! (mag) ? Apparent R magnitude from Massey
* (2002, Cat II/236)
real*8 e_Rmagm_3 (nr__3) ! (mag) ? Mean error on Rmagm
real*8 Rfmagg_3 (nr__3) ! (mag) ? Apparent Rf magnitude from GSC 2.3.2 (I/305)
real*8 e_Rfmagg_3 (nr__3) ! (mag) ? Mean error on Rfmagg
real*8 Imagz_3 (nr__3) ! (mag) ? Apparent I magnitude from Zaritsky et al.
* (2004AJ....128.1606Z)
real*8 e_Imagz_3 (nr__3) ! (mag) ? Mean error on Imagz
real*8 Jmag_3 (nr__3) ! (mag) ? 2MASS J apparent magnitude (II/246)
real*8 e_Jmag_3 (nr__3) ! (mag) ? Mean error on Jmag
real*8 Hmag_3 (nr__3) ! (mag) ? 2MASS H apparent magnitude (II/246)
real*8 e_Hmag_3 (nr__3) ! (mag) ? Mean error on Hmag
real*8 Kmag_3 (nr__3) ! (mag) ? 2MASS Ks apparent magnitude (II/246)
real*8 e_Kmag_3 (nr__3) ! (mag) ? Mean error on Kmag
real*8 F3_6_3 (nr__3) ! (mJy) ? Spitzer/IRAC 3.6 micron flux
real*8 e_F3_6_3 (nr__3) ! (mJy) ? Mean error on F3.6
real*8 F4_5_3 (nr__3) ! (mJy) ? Spitzer/IRAC 4.5 micron flux
real*8 e_F4_5_3 (nr__3) ! (mJy) ? Mean error on F4.5
real*8 F5_8_3 (nr__3) ! (mJy) ? Spitzer/IRAC 5.8 micron flux
real*8 e_F5_8_3 (nr__3) ! (mJy) ? Mean error on F5.8
real*8 F8_0_3 (nr__3) ! (mJy) ? Spitzer/IRAC 8.0 micron flux
real*8 e_F8_0_3 (nr__3) ! (mJy) ? Mean error on F8.0
real*8 F24_3 (nr__3) ! (mJy) ? Spitzer/MIPS 24 micron flux
real*8 e_F24_3 (nr__3) ! (mJy) ? Mean error on F24
real*8 F70_3 (nr__3) ! (mJy) ? Spitzer/MIPS 70 micron flux
real*8 e_F70_3 (nr__3) ! (mJy) ? Mean error on F70
real*8 F160_3 (nr__3) ! (mJy) ? Spitzer/MIPS 160 micron flux
real*8 e_F160_3 (nr__3) ! (mJy) ? Mean error on F160
integer*4 LumBB_3 (nr__3) ! (solLum) ? Black-body-based luminosity
* (cat_pag and cat_spec) or
* Estimated luminosity (cat_sg and cat_yso)
integer*4 LumSED_3 (nr__3) ! (solLum) ? SED-based luminosity (1)
integer*4 e_LumSED_3 (nr__3) ! (solLum) ? Mean error on LumSED (1)
integer*4 TeffSED_3 (nr__3) ! (K) ? SED-based effective temperature (1)
integer*4 e_TeffSED_3(nr__3) ! (K) ? Mean error on TeffSED (1)
real*8 E_B_V_1 (nr__3) ! (mag) ? SED-based E(B-V) (1)
real*8 e_E_B_V_1 (nr__3) ! (mag) ? Mean error on E(B-V) (1)
character*7 Xcat_3 (nr__3) ! ? Cross-correlation with catalogues
* of other types of objects (1) (3)
character*30 SpType_3 (nr__3) ! ? MK spectral type (1) (4)
character*30 Var_3 (nr__3) ! ? Variability remarks (1) (5)
character*1 n_Name_3 (nr__3) ! Reason why this object was discarded
* from our post-AGB sample (2)
*Note (1): For cat_pagb and cat_spec only.
*Note (2): Reason codes as follows (for cat_spec only):
* A = Not in the current version of the SAGE catalogue
* B = Does not obey the colour criterion
* C = Does not survive the black-body luminosity-cut (Sect. 2.3)
* D = Is listed as a cool carbon star by Kontizas et al.
* (2001, Cat. J/A+A/369/932) or an AGB star with spectral type by
* Trams et al. (1999A&A...346..843T)
* E = Is not a post-AGB star according to the spectral type
* F = Does not survive the SED luminosity-cut (Sect. 4.1)
*Note (3): Empirically confirmed cross-matches of the object in the catalogues of
* Sect. 5.1.3 with
* A = the RVT stars listed in the recent OGLE-III Catalogue of Variable
* Stars (Soszynski et al., 2008, Cat J/AcA/58/293)
* B = the RCrB stars in Soszynski et al. (2009, Cat J/AcA/59/335)
* C = the MSX post-AGB stars of Wood & Cohen (2001ASSL..265...71W)
* D = the YSOs from Seale et al. (2009, Cat. J/AJ/699/150)
* E = the AGB stars from Trams et al. (1999A&A...346..843T)
* F = the PNe in Reid & Parker (2006, Cat. J/MNRAS/373/521)
* G = the galaxies in the line of sight of the LMC from NED
*Note (4): All spectral types indicated with (A) were found in the catalogue of
* massive stars of Bonanos et al. (2009AJ....138.1003B), the others were
* determined based on our low-resolution, optical spectra.
*Note (5): For cat_spec only. The variability remarks are based on both the
* luminosity curves as discussed in Sect. 5.2 and cross-correlation with
* the catalogue of Fraser et al. (2008, Cat. J/AJ/136/1242). The latter
* is indicated with 'Fr_' followed by the correct sequence in subscript.
* Sequences 9, 0 and -99 indicate subsequently that the star is
* identified with the One-Year Artifact which is caused by the annual
* observing schedule of the MACHO project, that the star is outside the
* boundaries of any period-luminosity classification, or that Fraser et
* al. (2008, Cat. J/AJ/136/1242) were unable to classify this object.
* The other abbreviations used to characterise the variability of the
* different objects are Population II Cepheids (Ceph.) which include the
* subclass of RV Tauri stars (RVT), semi-regular variables (SRV) with
* the subclass of sequence-D variables (Seq. D) that were recognized
* by us, objects with lightcurves in Fig. 18 that show long-term, slow
* variations (SV), RCrB stars (RCrB), objects in regions of very strong
* nebular emission (em), Ae and Be stars (AB) and objects that do not
* show strong evidence of variability (NV). For all Population II
* Cepheids, the periods we computed between deep minima or two cycles of
* the light curve if the minima are of equal depth, are listed in days.
C=============================================================================
C Loading file 'cat_pagb.dat' ! Post-AGB candidates
C Format for file interpretation
1 format(
+ A19,2X,F10.6,2X,F10.6,2X,A13,1X,F9.3,1X,F9.3,1X,F9.2,1X,F9.2,
+ 1X,F9.3,1X,F9.3,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.3,1X,
+ F9.3,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,
+ 1X,F9.2,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,
+ F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,
+ 1X,F9.3,1X,F9.3,1X,F9.4,1X,F9.4,1X,F9.2,1X,F9.3,1X,F9.1,1X,
+ F9.2,1X,I7,1X,I7,1X,I7,1X,I7,1X,I7,1X,F7.2,1X,F7.2,1X,A7,1X,
+ A30,1X,A30,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'cat_pagb.dat')
write(6,*) '....Loading file: cat_pagb.dat'
do i__=1,1407
read(1,'(A606)')ar__
read(ar__,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),Type(i__),Umagz(i__),
+ e_Umagz(i__),Umagm(i__),e_Umagm(i__),Bmagz(i__),e_Bmagz(i__),
+ Bmagm(i__),e_Bmagm(i__),Bmagg(i__),e_Bmagg(i__),Vmagz(i__),
+ e_Vmagz(i__),Vmagm(i__),e_Vmagm(i__),Vmagg(i__),e_Vmagg(i__),
+ Rmagm(i__),e_Rmagm(i__),Rfmagg(i__),e_Rfmagg(i__),Imagz(i__),
+ e_Imagz(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),
+ Kmag(i__),e_Kmag(i__),F3_6(i__),e_F3_6(i__),F4_5(i__),
+ e_F4_5(i__),F5_8(i__),e_F5_8(i__),F8_0(i__),e_F8_0(i__),
+ F24(i__),e_F24(i__),F70(i__),e_F70(i__),F160(i__),e_F160(i__),
+ LumBB(i__),LumSED(i__),e_LumSED(i__),TeffSED(i__),
+ e_TeffSED(i__),E_B_V(i__),e_E_B_V(i__),Xcat(i__),SpType(i__),
+ Var(i__),n_Name(i__)
if(ar__(60:68) .EQ. '') Umagz(i__) = rNULL__
if(ar__(70:78) .EQ. '') e_Umagz(i__) = rNULL__
if(ar__(80:88) .EQ. '') Umagm(i__) = rNULL__
if(ar__(90:98) .EQ. '') e_Umagm(i__) = rNULL__
if(ar__(100:108) .EQ. '') Bmagz(i__) = rNULL__
if(ar__(110:118) .EQ. '') e_Bmagz(i__) = rNULL__
if(ar__(120:128) .EQ. '') Bmagm(i__) = rNULL__
if(ar__(130:138) .EQ. '') e_Bmagm(i__) = rNULL__
if(ar__(140:148) .EQ. '') Bmagg(i__) = rNULL__
if(ar__(150:158) .EQ. '') e_Bmagg(i__) = rNULL__
if(ar__(160:168) .EQ. '') Vmagz(i__) = rNULL__
if(ar__(170:178) .EQ. '') e_Vmagz(i__) = rNULL__
if(ar__(180:188) .EQ. '') Vmagm(i__) = rNULL__
if(ar__(190:198) .EQ. '') e_Vmagm(i__) = rNULL__
if(ar__(200:208) .EQ. '') Vmagg(i__) = rNULL__
if(ar__(210:218) .EQ. '') e_Vmagg(i__) = rNULL__
if(ar__(220:228) .EQ. '') Rmagm(i__) = rNULL__
if(ar__(230:238) .EQ. '') e_Rmagm(i__) = rNULL__
if(ar__(240:248) .EQ. '') Rfmagg(i__) = rNULL__
if(ar__(250:258) .EQ. '') e_Rfmagg(i__) = rNULL__
if(ar__(260:268) .EQ. '') Imagz(i__) = rNULL__
if(ar__(270:278) .EQ. '') e_Imagz(i__) = rNULL__
if(ar__(280:288) .EQ. '') Jmag(i__) = rNULL__
if(ar__(290:298) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__(300:308) .EQ. '') Hmag(i__) = rNULL__
if(ar__(310:318) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__(320:328) .EQ. '') Kmag(i__) = rNULL__
if(ar__(330:338) .EQ. '') e_Kmag(i__) = rNULL__
if(ar__(340:348) .EQ. '') F3_6(i__) = rNULL__
if(ar__(350:358) .EQ. '') e_F3_6(i__) = rNULL__
if(ar__(360:368) .EQ. '') F4_5(i__) = rNULL__
if(ar__(370:378) .EQ. '') e_F4_5(i__) = rNULL__
if(ar__(380:388) .EQ. '') F5_8(i__) = rNULL__
if(ar__(390:398) .EQ. '') e_F5_8(i__) = rNULL__
if(ar__(400:408) .EQ. '') F8_0(i__) = rNULL__
if(ar__(410:418) .EQ. '') e_F8_0(i__) = rNULL__
if(ar__(420:428) .EQ. '') F24(i__) = rNULL__
if(ar__(430:438) .EQ. '') e_F24(i__) = rNULL__
if(ar__(440:448) .EQ. '') F70(i__) = rNULL__
if(ar__(450:458) .EQ. '') e_F70(i__) = rNULL__
if(ar__(460:468) .EQ. '') F160(i__) = rNULL__
if(ar__(470:478) .EQ. '') e_F160(i__) = rNULL__
if(ar__(480:486) .EQ. '') LumBB(i__) = iNULL__
if(ar__(488:494) .EQ. '') LumSED(i__) = iNULL__
if(ar__(496:502) .EQ. '') e_LumSED(i__) = iNULL__
if(ar__(504:510) .EQ. '') TeffSED(i__) = iNULL__
if(ar__(512:518) .EQ. '') e_TeffSED(i__) = iNULL__
if(ar__(520:526) .EQ. '') E_B_V(i__) = rNULL__
if(ar__(528:534) .EQ. '') e_E_B_V(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),Type(i__),Umagz(i__),
+ e_Umagz(i__),Umagm(i__),e_Umagm(i__),Bmagz(i__),e_Bmagz(i__),
+ Bmagm(i__),e_Bmagm(i__),Bmagg(i__),e_Bmagg(i__),Vmagz(i__),
+ e_Vmagz(i__),Vmagm(i__),e_Vmagm(i__),Vmagg(i__),e_Vmagg(i__),
+ Rmagm(i__),e_Rmagm(i__),Rfmagg(i__),e_Rfmagg(i__),Imagz(i__),
+ e_Imagz(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),
+ Kmag(i__),e_Kmag(i__),F3_6(i__),e_F3_6(i__),F4_5(i__),
+ e_F4_5(i__),F5_8(i__),e_F5_8(i__),F8_0(i__),e_F8_0(i__),
+ F24(i__),e_F24(i__),F70(i__),e_F70(i__),F160(i__),e_F160(i__),
+ LumBB(i__),LumSED(i__),e_LumSED(i__),TeffSED(i__),
+ e_TeffSED(i__),E_B_V(i__),e_E_B_V(i__),Xcat(i__),SpType(i__),
+ Var(i__),n_Name(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat_spec.dat' ! Discarded objects with a spectral type
C Format for file interpretation
2 format(
+ A19,2X,F10.6,2X,F10.6,2X,A13,1X,F9.3,1X,F9.3,1X,F9.2,1X,F9.2,
+ 1X,F9.3,1X,F9.3,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.3,1X,
+ F9.3,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,
+ 1X,F9.2,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,
+ F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,
+ 1X,F9.3,1X,F9.3,1X,F9.4,1X,F9.4,1X,F9.2,1X,F9.3,1X,F9.1,1X,
+ F9.2,1X,I7,1X,I7,1X,I7,1X,I7,1X,I7,1X,F7.2,1X,F7.2,1X,A7,1X,
+ A30,1X,A30,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'cat_spec.dat')
write(6,*) '....Loading file: cat_spec.dat'
do i__=1,89
read(1,'(A606)')ar__1
read(ar__1,2)
+ Name_1(i__),RAdeg_1(i__),DEdeg_1(i__),Type_1(i__),
+ Umagz_1(i__),e_Umagz_1(i__),Umagm_1(i__),e_Umagm_1(i__),
+ Bmagz_1(i__),e_Bmagz_1(i__),Bmagm_1(i__),e_Bmagm_1(i__),
+ Bmagg_1(i__),e_Bmagg_1(i__),Vmagz_1(i__),e_Vmagz_1(i__),
+ Vmagm_1(i__),e_Vmagm_1(i__),Vmagg_1(i__),e_Vmagg_1(i__),
+ Rmagm_1(i__),e_Rmagm_1(i__),Rfmagg_1(i__),e_Rfmagg_1(i__),
+ Imagz_1(i__),e_Imagz_1(i__),Jmag_1(i__),e_Jmag_1(i__),
+ Hmag_1(i__),e_Hmag_1(i__),Kmag_1(i__),e_Kmag_1(i__),
+ F3_6_1(i__),e_F3_6_1(i__),F4_5_1(i__),e_F4_5_1(i__),
+ F5_8_1(i__),e_F5_8_1(i__),F8_0_1(i__),e_F8_0_1(i__),
+ F24_1(i__),e_F24_1(i__),F70_1(i__),e_F70_1(i__),F160_1(i__),
+ e_F160_1(i__),LumBB_1(i__),LumSED_1(i__),e_LumSED_1(i__),
+ TeffSED_1(i__),e_TeffSED_1(i__),E_B_V_1(i__),e_E_B_V_1(i__),
+ Xcat_1(i__),SpType_1(i__),Var_1(i__),n_Name_1(i__)
if(ar__1(60:68) .EQ. '') Umagz_1(i__) = rNULL__
if(ar__1(70:78) .EQ. '') e_Umagz_1(i__) = rNULL__
if(ar__1(80:88) .EQ. '') Umagm_1(i__) = rNULL__
if(ar__1(90:98) .EQ. '') e_Umagm_1(i__) = rNULL__
if(ar__1(100:108) .EQ. '') Bmagz_1(i__) = rNULL__
if(ar__1(110:118) .EQ. '') e_Bmagz_1(i__) = rNULL__
if(ar__1(120:128) .EQ. '') Bmagm_1(i__) = rNULL__
if(ar__1(130:138) .EQ. '') e_Bmagm_1(i__) = rNULL__
if(ar__1(140:148) .EQ. '') Bmagg_1(i__) = rNULL__
if(ar__1(150:158) .EQ. '') e_Bmagg_1(i__) = rNULL__
if(ar__1(160:168) .EQ. '') Vmagz_1(i__) = rNULL__
if(ar__1(170:178) .EQ. '') e_Vmagz_1(i__) = rNULL__
if(ar__1(180:188) .EQ. '') Vmagm_1(i__) = rNULL__
if(ar__1(190:198) .EQ. '') e_Vmagm_1(i__) = rNULL__
if(ar__1(200:208) .EQ. '') Vmagg_1(i__) = rNULL__
if(ar__1(210:218) .EQ. '') e_Vmagg_1(i__) = rNULL__
if(ar__1(220:228) .EQ. '') Rmagm_1(i__) = rNULL__
if(ar__1(230:238) .EQ. '') e_Rmagm_1(i__) = rNULL__
if(ar__1(240:248) .EQ. '') Rfmagg_1(i__) = rNULL__
if(ar__1(250:258) .EQ. '') e_Rfmagg_1(i__) = rNULL__
if(ar__1(260:268) .EQ. '') Imagz_1(i__) = rNULL__
if(ar__1(270:278) .EQ. '') e_Imagz_1(i__) = rNULL__
if(ar__1(280:288) .EQ. '') Jmag_1(i__) = rNULL__
if(ar__1(290:298) .EQ. '') e_Jmag_1(i__) = rNULL__
if(ar__1(300:308) .EQ. '') Hmag_1(i__) = rNULL__
if(ar__1(310:318) .EQ. '') e_Hmag_1(i__) = rNULL__
if(ar__1(320:328) .EQ. '') Kmag_1(i__) = rNULL__
if(ar__1(330:338) .EQ. '') e_Kmag_1(i__) = rNULL__
if(ar__1(340:348) .EQ. '') F3_6_1(i__) = rNULL__
if(ar__1(350:358) .EQ. '') e_F3_6_1(i__) = rNULL__
if(ar__1(360:368) .EQ. '') F4_5_1(i__) = rNULL__
if(ar__1(370:378) .EQ. '') e_F4_5_1(i__) = rNULL__
if(ar__1(380:388) .EQ. '') F5_8_1(i__) = rNULL__
if(ar__1(390:398) .EQ. '') e_F5_8_1(i__) = rNULL__
if(ar__1(400:408) .EQ. '') F8_0_1(i__) = rNULL__
if(ar__1(410:418) .EQ. '') e_F8_0_1(i__) = rNULL__
if(ar__1(420:428) .EQ. '') F24_1(i__) = rNULL__
if(ar__1(430:438) .EQ. '') e_F24_1(i__) = rNULL__
if(ar__1(440:448) .EQ. '') F70_1(i__) = rNULL__
if(ar__1(450:458) .EQ. '') e_F70_1(i__) = rNULL__
if(ar__1(460:468) .EQ. '') F160_1(i__) = rNULL__
if(ar__1(470:478) .EQ. '') e_F160_1(i__) = rNULL__
if(ar__1(480:486) .EQ. '') LumBB_1(i__) = iNULL__
if(ar__1(488:494) .EQ. '') LumSED_1(i__) = iNULL__
if(ar__1(496:502) .EQ. '') e_LumSED_1(i__) = iNULL__
if(ar__1(504:510) .EQ. '') TeffSED_1(i__) = iNULL__
if(ar__1(512:518) .EQ. '') e_TeffSED_1(i__) = iNULL__
if(ar__1(520:526) .EQ. '') E_B_V_1(i__) = rNULL__
if(ar__1(528:534) .EQ. '') e_E_B_V_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),RAdeg_1(i__),DEdeg_1(i__),Type_1(i__),
+ Umagz_1(i__),e_Umagz_1(i__),Umagm_1(i__),e_Umagm_1(i__),
+ Bmagz_1(i__),e_Bmagz_1(i__),Bmagm_1(i__),e_Bmagm_1(i__),
+ Bmagg_1(i__),e_Bmagg_1(i__),Vmagz_1(i__),e_Vmagz_1(i__),
+ Vmagm_1(i__),e_Vmagm_1(i__),Vmagg_1(i__),e_Vmagg_1(i__),
+ Rmagm_1(i__),e_Rmagm_1(i__),Rfmagg_1(i__),e_Rfmagg_1(i__),
+ Imagz_1(i__),e_Imagz_1(i__),Jmag_1(i__),e_Jmag_1(i__),
+ Hmag_1(i__),e_Hmag_1(i__),Kmag_1(i__),e_Kmag_1(i__),
+ F3_6_1(i__),e_F3_6_1(i__),F4_5_1(i__),e_F4_5_1(i__),
+ F5_8_1(i__),e_F5_8_1(i__),F8_0_1(i__),e_F8_0_1(i__),
+ F24_1(i__),e_F24_1(i__),F70_1(i__),e_F70_1(i__),F160_1(i__),
+ e_F160_1(i__),LumBB_1(i__),LumSED_1(i__),e_LumSED_1(i__),
+ TeffSED_1(i__),e_TeffSED_1(i__),E_B_V_1(i__),e_E_B_V_1(i__),
+ Xcat_1(i__),SpType_1(i__),Var_1(i__),n_Name_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat_yso.dat' ! Discarded possible YSOs
C Format for file interpretation
3 format(
+ A19,2X,F10.6,2X,F10.6,2X,A13,1X,F9.3,1X,F9.3,1X,F9.2,1X,F9.2,
+ 1X,F9.3,1X,F9.3,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.3,1X,
+ F9.3,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,
+ 1X,F9.2,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,
+ F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,
+ 1X,F9.3,1X,F9.3,1X,F9.4,1X,F9.4,1X,F9.2,1X,F9.3,1X,F9.1,1X,
+ F9.2,1X,I7,1X,I7,1X,I7,1X,I7,1X,I7,1X,F7.2,1X,F7.2,1X,A7,1X,
+ A30,1X,A30,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'cat_yso.dat')
write(6,*) '....Loading file: cat_yso.dat'
do i__=1,6823
read(1,'(A606)')ar__2
read(ar__2,3)
+ Name_2(i__),RAdeg_2(i__),DEdeg_2(i__),Type_2(i__),
+ Umagz_2(i__),e_Umagz_2(i__),Umagm_2(i__),e_Umagm_2(i__),
+ Bmagz_2(i__),e_Bmagz_2(i__),Bmagm_2(i__),e_Bmagm_2(i__),
+ Bmagg_2(i__),e_Bmagg_2(i__),Vmagz_2(i__),e_Vmagz_2(i__),
+ Vmagm_2(i__),e_Vmagm_2(i__),Vmagg_2(i__),e_Vmagg_2(i__),
+ Rmagm_2(i__),e_Rmagm_2(i__),Rfmagg_2(i__),e_Rfmagg_2(i__),
+ Imagz_2(i__),e_Imagz_2(i__),Jmag_2(i__),e_Jmag_2(i__),
+ Hmag_2(i__),e_Hmag_2(i__),Kmag_2(i__),e_Kmag_2(i__),
+ F3_6_2(i__),e_F3_6_2(i__),F4_5_2(i__),e_F4_5_2(i__),
+ F5_8_2(i__),e_F5_8_2(i__),F8_0_2(i__),e_F8_0_2(i__),
+ F24_2(i__),e_F24_2(i__),F70_2(i__),e_F70_2(i__),F160_2(i__),
+ e_F160_2(i__),LumBB_2(i__),LumSED_2(i__),e_LumSED_2(i__),
+ TeffSED_2(i__),e_TeffSED_2(i__),E_B_V_1(i__),e_E_B_V_1(i__),
+ Xcat_2(i__),SpType_2(i__),Var_2(i__),n_Name_2(i__)
if(ar__2(60:68) .EQ. '') Umagz_2(i__) = rNULL__
if(ar__2(70:78) .EQ. '') e_Umagz_2(i__) = rNULL__
if(ar__2(80:88) .EQ. '') Umagm_2(i__) = rNULL__
if(ar__2(90:98) .EQ. '') e_Umagm_2(i__) = rNULL__
if(ar__2(100:108) .EQ. '') Bmagz_2(i__) = rNULL__
if(ar__2(110:118) .EQ. '') e_Bmagz_2(i__) = rNULL__
if(ar__2(120:128) .EQ. '') Bmagm_2(i__) = rNULL__
if(ar__2(130:138) .EQ. '') e_Bmagm_2(i__) = rNULL__
if(ar__2(140:148) .EQ. '') Bmagg_2(i__) = rNULL__
if(ar__2(150:158) .EQ. '') e_Bmagg_2(i__) = rNULL__
if(ar__2(160:168) .EQ. '') Vmagz_2(i__) = rNULL__
if(ar__2(170:178) .EQ. '') e_Vmagz_2(i__) = rNULL__
if(ar__2(180:188) .EQ. '') Vmagm_2(i__) = rNULL__
if(ar__2(190:198) .EQ. '') e_Vmagm_2(i__) = rNULL__
if(ar__2(200:208) .EQ. '') Vmagg_2(i__) = rNULL__
if(ar__2(210:218) .EQ. '') e_Vmagg_2(i__) = rNULL__
if(ar__2(220:228) .EQ. '') Rmagm_2(i__) = rNULL__
if(ar__2(230:238) .EQ. '') e_Rmagm_2(i__) = rNULL__
if(ar__2(240:248) .EQ. '') Rfmagg_2(i__) = rNULL__
if(ar__2(250:258) .EQ. '') e_Rfmagg_2(i__) = rNULL__
if(ar__2(260:268) .EQ. '') Imagz_2(i__) = rNULL__
if(ar__2(270:278) .EQ. '') e_Imagz_2(i__) = rNULL__
if(ar__2(280:288) .EQ. '') Jmag_2(i__) = rNULL__
if(ar__2(290:298) .EQ. '') e_Jmag_2(i__) = rNULL__
if(ar__2(300:308) .EQ. '') Hmag_2(i__) = rNULL__
if(ar__2(310:318) .EQ. '') e_Hmag_2(i__) = rNULL__
if(ar__2(320:328) .EQ. '') Kmag_2(i__) = rNULL__
if(ar__2(330:338) .EQ. '') e_Kmag_2(i__) = rNULL__
if(ar__2(340:348) .EQ. '') F3_6_2(i__) = rNULL__
if(ar__2(350:358) .EQ. '') e_F3_6_2(i__) = rNULL__
if(ar__2(360:368) .EQ. '') F4_5_2(i__) = rNULL__
if(ar__2(370:378) .EQ. '') e_F4_5_2(i__) = rNULL__
if(ar__2(380:388) .EQ. '') F5_8_2(i__) = rNULL__
if(ar__2(390:398) .EQ. '') e_F5_8_2(i__) = rNULL__
if(ar__2(400:408) .EQ. '') F8_0_2(i__) = rNULL__
if(ar__2(410:418) .EQ. '') e_F8_0_2(i__) = rNULL__
if(ar__2(420:428) .EQ. '') F24_2(i__) = rNULL__
if(ar__2(430:438) .EQ. '') e_F24_2(i__) = rNULL__
if(ar__2(440:448) .EQ. '') F70_2(i__) = rNULL__
if(ar__2(450:458) .EQ. '') e_F70_2(i__) = rNULL__
if(ar__2(460:468) .EQ. '') F160_2(i__) = rNULL__
if(ar__2(470:478) .EQ. '') e_F160_2(i__) = rNULL__
if(ar__2(480:486) .EQ. '') LumBB_2(i__) = iNULL__
if(ar__2(488:494) .EQ. '') LumSED_2(i__) = iNULL__
if(ar__2(496:502) .EQ. '') e_LumSED_2(i__) = iNULL__
if(ar__2(504:510) .EQ. '') TeffSED_2(i__) = iNULL__
if(ar__2(512:518) .EQ. '') e_TeffSED_2(i__) = iNULL__
if(ar__2(520:526) .EQ. '') E_B_V_1(i__) = rNULL__
if(ar__2(528:534) .EQ. '') e_E_B_V_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),RAdeg_2(i__),DEdeg_2(i__),Type_2(i__),
+ Umagz_2(i__),e_Umagz_2(i__),Umagm_2(i__),e_Umagm_2(i__),
+ Bmagz_2(i__),e_Bmagz_2(i__),Bmagm_2(i__),e_Bmagm_2(i__),
+ Bmagg_2(i__),e_Bmagg_2(i__),Vmagz_2(i__),e_Vmagz_2(i__),
+ Vmagm_2(i__),e_Vmagm_2(i__),Vmagg_2(i__),e_Vmagg_2(i__),
+ Rmagm_2(i__),e_Rmagm_2(i__),Rfmagg_2(i__),e_Rfmagg_2(i__),
+ Imagz_2(i__),e_Imagz_2(i__),Jmag_2(i__),e_Jmag_2(i__),
+ Hmag_2(i__),e_Hmag_2(i__),Kmag_2(i__),e_Kmag_2(i__),
+ F3_6_2(i__),e_F3_6_2(i__),F4_5_2(i__),e_F4_5_2(i__),
+ F5_8_2(i__),e_F5_8_2(i__),F8_0_2(i__),e_F8_0_2(i__),
+ F24_2(i__),e_F24_2(i__),F70_2(i__),e_F70_2(i__),F160_2(i__),
+ e_F160_2(i__),LumBB_2(i__),LumSED_2(i__),e_LumSED_2(i__),
+ TeffSED_2(i__),e_TeffSED_2(i__),E_B_V_1(i__),e_E_B_V_1(i__),
+ Xcat_2(i__),SpType_2(i__),Var_2(i__),n_Name_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat_sg.dat' ! Discarded possible supergiants
C Format for file interpretation
4 format(
+ A19,2X,F10.6,2X,F10.6,2X,A13,1X,F9.3,1X,F9.3,1X,F9.2,1X,F9.2,
+ 1X,F9.3,1X,F9.3,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.3,1X,
+ F9.3,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,
+ 1X,F9.2,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,
+ F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,
+ 1X,F9.3,1X,F9.3,1X,F9.4,1X,F9.4,1X,F9.2,1X,F9.3,1X,F9.1,1X,
+ F9.2,1X,I7,1X,I7,1X,I7,1X,I7,1X,I7,1X,F7.2,1X,F7.2,1X,A7,1X,
+ A30,1X,A30,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'cat_sg.dat')
write(6,*) '....Loading file: cat_sg.dat'
do i__=1,286
read(1,'(A606)')ar__3
read(ar__3,4)
+ Name_3(i__),RAdeg_3(i__),DEdeg_3(i__),Type_3(i__),
+ Umagz_3(i__),e_Umagz_3(i__),Umagm_3(i__),e_Umagm_3(i__),
+ Bmagz_3(i__),e_Bmagz_3(i__),Bmagm_3(i__),e_Bmagm_3(i__),
+ Bmagg_3(i__),e_Bmagg_3(i__),Vmagz_3(i__),e_Vmagz_3(i__),
+ Vmagm_3(i__),e_Vmagm_3(i__),Vmagg_3(i__),e_Vmagg_3(i__),
+ Rmagm_3(i__),e_Rmagm_3(i__),Rfmagg_3(i__),e_Rfmagg_3(i__),
+ Imagz_3(i__),e_Imagz_3(i__),Jmag_3(i__),e_Jmag_3(i__),
+ Hmag_3(i__),e_Hmag_3(i__),Kmag_3(i__),e_Kmag_3(i__),
+ F3_6_3(i__),e_F3_6_3(i__),F4_5_3(i__),e_F4_5_3(i__),
+ F5_8_3(i__),e_F5_8_3(i__),F8_0_3(i__),e_F8_0_3(i__),
+ F24_3(i__),e_F24_3(i__),F70_3(i__),e_F70_3(i__),F160_3(i__),
+ e_F160_3(i__),LumBB_3(i__),LumSED_3(i__),e_LumSED_3(i__),
+ TeffSED_3(i__),e_TeffSED_3(i__),E_B_V_1(i__),e_E_B_V_1(i__),
+ Xcat_3(i__),SpType_3(i__),Var_3(i__),n_Name_3(i__)
if(ar__3(60:68) .EQ. '') Umagz_3(i__) = rNULL__
if(ar__3(70:78) .EQ. '') e_Umagz_3(i__) = rNULL__
if(ar__3(80:88) .EQ. '') Umagm_3(i__) = rNULL__
if(ar__3(90:98) .EQ. '') e_Umagm_3(i__) = rNULL__
if(ar__3(100:108) .EQ. '') Bmagz_3(i__) = rNULL__
if(ar__3(110:118) .EQ. '') e_Bmagz_3(i__) = rNULL__
if(ar__3(120:128) .EQ. '') Bmagm_3(i__) = rNULL__
if(ar__3(130:138) .EQ. '') e_Bmagm_3(i__) = rNULL__
if(ar__3(140:148) .EQ. '') Bmagg_3(i__) = rNULL__
if(ar__3(150:158) .EQ. '') e_Bmagg_3(i__) = rNULL__
if(ar__3(160:168) .EQ. '') Vmagz_3(i__) = rNULL__
if(ar__3(170:178) .EQ. '') e_Vmagz_3(i__) = rNULL__
if(ar__3(180:188) .EQ. '') Vmagm_3(i__) = rNULL__
if(ar__3(190:198) .EQ. '') e_Vmagm_3(i__) = rNULL__
if(ar__3(200:208) .EQ. '') Vmagg_3(i__) = rNULL__
if(ar__3(210:218) .EQ. '') e_Vmagg_3(i__) = rNULL__
if(ar__3(220:228) .EQ. '') Rmagm_3(i__) = rNULL__
if(ar__3(230:238) .EQ. '') e_Rmagm_3(i__) = rNULL__
if(ar__3(240:248) .EQ. '') Rfmagg_3(i__) = rNULL__
if(ar__3(250:258) .EQ. '') e_Rfmagg_3(i__) = rNULL__
if(ar__3(260:268) .EQ. '') Imagz_3(i__) = rNULL__
if(ar__3(270:278) .EQ. '') e_Imagz_3(i__) = rNULL__
if(ar__3(280:288) .EQ. '') Jmag_3(i__) = rNULL__
if(ar__3(290:298) .EQ. '') e_Jmag_3(i__) = rNULL__
if(ar__3(300:308) .EQ. '') Hmag_3(i__) = rNULL__
if(ar__3(310:318) .EQ. '') e_Hmag_3(i__) = rNULL__
if(ar__3(320:328) .EQ. '') Kmag_3(i__) = rNULL__
if(ar__3(330:338) .EQ. '') e_Kmag_3(i__) = rNULL__
if(ar__3(340:348) .EQ. '') F3_6_3(i__) = rNULL__
if(ar__3(350:358) .EQ. '') e_F3_6_3(i__) = rNULL__
if(ar__3(360:368) .EQ. '') F4_5_3(i__) = rNULL__
if(ar__3(370:378) .EQ. '') e_F4_5_3(i__) = rNULL__
if(ar__3(380:388) .EQ. '') F5_8_3(i__) = rNULL__
if(ar__3(390:398) .EQ. '') e_F5_8_3(i__) = rNULL__
if(ar__3(400:408) .EQ. '') F8_0_3(i__) = rNULL__
if(ar__3(410:418) .EQ. '') e_F8_0_3(i__) = rNULL__
if(ar__3(420:428) .EQ. '') F24_3(i__) = rNULL__
if(ar__3(430:438) .EQ. '') e_F24_3(i__) = rNULL__
if(ar__3(440:448) .EQ. '') F70_3(i__) = rNULL__
if(ar__3(450:458) .EQ. '') e_F70_3(i__) = rNULL__
if(ar__3(460:468) .EQ. '') F160_3(i__) = rNULL__
if(ar__3(470:478) .EQ. '') e_F160_3(i__) = rNULL__
if(ar__3(480:486) .EQ. '') LumBB_3(i__) = iNULL__
if(ar__3(488:494) .EQ. '') LumSED_3(i__) = iNULL__
if(ar__3(496:502) .EQ. '') e_LumSED_3(i__) = iNULL__
if(ar__3(504:510) .EQ. '') TeffSED_3(i__) = iNULL__
if(ar__3(512:518) .EQ. '') e_TeffSED_3(i__) = iNULL__
if(ar__3(520:526) .EQ. '') E_B_V_1(i__) = rNULL__
if(ar__3(528:534) .EQ. '') e_E_B_V_1(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3(i__),RAdeg_3(i__),DEdeg_3(i__),Type_3(i__),
+ Umagz_3(i__),e_Umagz_3(i__),Umagm_3(i__),e_Umagm_3(i__),
+ Bmagz_3(i__),e_Bmagz_3(i__),Bmagm_3(i__),e_Bmagm_3(i__),
+ Bmagg_3(i__),e_Bmagg_3(i__),Vmagz_3(i__),e_Vmagz_3(i__),
+ Vmagm_3(i__),e_Vmagm_3(i__),Vmagg_3(i__),e_Vmagg_3(i__),
+ Rmagm_3(i__),e_Rmagm_3(i__),Rfmagg_3(i__),e_Rfmagg_3(i__),
+ Imagz_3(i__),e_Imagz_3(i__),Jmag_3(i__),e_Jmag_3(i__),
+ Hmag_3(i__),e_Hmag_3(i__),Kmag_3(i__),e_Kmag_3(i__),
+ F3_6_3(i__),e_F3_6_3(i__),F4_5_3(i__),e_F4_5_3(i__),
+ F5_8_3(i__),e_F5_8_3(i__),F8_0_3(i__),e_F8_0_3(i__),
+ F24_3(i__),e_F24_3(i__),F70_3(i__),e_F70_3(i__),F160_3(i__),
+ e_F160_3(i__),LumBB_3(i__),LumSED_3(i__),e_LumSED_3(i__),
+ TeffSED_3(i__),e_TeffSED_3(i__),E_B_V_1(i__),e_E_B_V_1(i__),
+ Xcat_3(i__),SpType_3(i__),Var_3(i__),n_Name_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end