Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/620/A128 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-Sep-09
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/620/A128 Gaia DR2 study of Herbig Ae/Be stars (Vioque+, 2018)
*================================================================================
*Gaia DR2 study of Herbig Ae/Be stars.
* Vioque M., Oudmaijer R.D., Baines D., Mendigutia I., Perez-Martinez R.
* <Astron. Astrophys. 620, A128 (2018)>
* =2018A&A...620A.128V (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hqsample.dat' ! Main parameters of each Herbig Ae/Be star belonging
to the high quality sample of 218 sources (table 1)
integer*4 nr__
parameter (nr__=218) ! Number of records
character*346 ar__ ! Full-size record
C Position composed of: RAh RAm RAs DE- DEd DEm DEs (Epoch=2015.5)
character*13 Name (nr__) ! Name of the Herbig Ae/Be star
integer*4 RAh (nr__) ! (h) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 RAm (nr__) ! (min) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
real*4 RAs (nr__) ! (s) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
character*1 DE_ (nr__) ! Declination sign from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEd (nr__) ! (deg) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEm (nr__) ! (arcmin) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEs (nr__) ! (arcsec) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
real*8 plx (nr__) ! (mas) Parallax from Gaia DR2
real*8 e_plx (nr__) ! (mas) Parallax uncertainty from Gaia DR2
real*4 Dist (nr__) ! (pc) Distance derived using an exponentially
* decreasing density prior
real*4 E_Dist (nr__) ! (pc) Upper error on the Distance, to the 95th
* percentile confidence interval on the
* Distance estimate
real*4 e_Dist_1 (nr__) ! (pc) Lower error on the Distance, to the 5th
* percentile confidence interval on the
* Distance estimate
integer*4 Teff (nr__) ! (K) Effective temperature from the literature (1)
integer*4 E_Teff (nr__) ! (K) Upper error of Teff
integer*4 e_Teff_1 (nr__) ! (K) Lower error of Teff
real*4 LogL (nr__) ! ([Lsun]) Decimal logarithm of the luminosity
real*4 E_LogL (nr__) ! ([Lsun]) Upper error of LogL
real*4 e_LogL_1 (nr__) ! ([Lsun]) Lower error of LogL
real*4 Av (nr__) ! (mag) Total extinction (1)
real*4 E_Av (nr__) ! (mag) Upper error of Av
real*4 e_Av_1 (nr__) ! (mag) Lower error of Av
real*4 V (nr__) ! (mag) Johnson V band magnitude (1)
character*3 Bin (nr__) ! [Yes ] Binarity status: "Yes" if the source
* has been classified as binary,
* nothing otherwise
integer*4 r_Bin (nr__) ! ?=- References for Binary column
* (in ref.dat file)
real*4 E_NIR (nr__) ! ?=- Infrared excess in the 1.24-3.4um
* interval. Derived following Equation 2 of
* the paper
real*4 E_E_NIR (nr__) ! ?=- Upper error of Near_IR_excess
real*4 e_E_NIR_1 (nr__) ! ?=- Lower error of Near_IR_excess
real*4 E_MIR (nr__) ! ?=- Infrared excess in the 3.4-22um
* interval. Derived following Equation 2 of
* the paper
real*4 E_E_MIR (nr__) ! ?=- Upper error of Mid_IR_excess
real*4 e_E_MIR_1 (nr__) ! ?=- Lower error of Mid_IR_excess
real*8 EWHa (nr__) ! (0.1nm) ?=- Halpha equivalent width, from the
* literature. Negative numbers indicate
* emission
real*4 e_EWHa (nr__) ! (0.1nm) ?=- Uncertainty of EWHa, from the literature
integer*4 r_EWHa (nr__) ! ?=- References for EWHa and e_EWHa columns
* (in refs.dat file)
character*1 Hashape (nr__) ! Line profile of the Halpha line, from
* the literature (2)
integer*4 r_Hashape (nr__) ! ?=- References for Hashape
* (in refs.dat file)
real*4 Vi (nr__) ! ?=- Variability indicator. Derived following
* Equation 4 of the paper (3)
character*3 UXOR (nr__) ! [Yes ] UXOR status: "Yes" if the source has
* been classified as UXOR type,
* nothing otherwise
real*4 Mass (nr__) ! (Msun) ?=- Mass of the source
real*4 E_Mass (nr__) ! (Msun) ?=- Upper error of the Mass
real*4 e_Mass_1 (nr__) ! (Msun) ?=- Lower error of the Mass
real*4 Age (nr__) ! (Myr) ?=- Age of the source
real*4 E_Age (nr__) ! (Myr) ?=- Upper error of the Age
real*4 e_Age_1 (nr__) ! (Myr) ?=- Lower error of the Age
real*4 E_J (nr__) ! ?=- Infrared excess at the J band (1.24um) (4)
real*4 E_H (nr__) ! ?=- Infrared excess at the H band (1.66um) (4)
real*8 E_Ks (nr__) ! ?=- Infrared excess at the
* Ks band (2.16um) (4)
real*8 E_W1 (nr__) ! ?=- Infrared excess at the W1 band (3.4um) (4)
real*8 E_W2 (nr__) ! ?=- Infrared excess at the W2 band (4.6um) (4)
real*8 E_W3 (nr__) ! ?=- Infrared excess at the W3 band (12um) (4)
real*8 E_W4 (nr__) ! ?=- Infrared excess at the W4 band (22um) (4)
character*3 The (nr__) ! [Yes ] Yes if the source appears in Table 1
* of The et al. (1994A&AS..104..315T),
* nothing otherwise
*Note (1): Atmospheric parameters Teff, AV and V taken from the following
* sources in order of choice:
* Fairlamb et al., 2015MNRAS.453..976F; Montesinos et al., 2009A&A...495..901M;
* Hernandez et al., 2004AJ....127.1682H; Mendigutia et al., 2012A&A...543A..59M;
* Carmona et al., 2010A&A...517A..67C; Chen et al., 2016NewA...44....1C;
* Alecian et al., 2013MNRAS.429.1001A; Sartori et al., 2010AJ....139...27S;
* Manoj et al., 2006ApJ...653..657M; Hernandez et al., 2005AJ....129..856H;
* Vieira et al., 2003AJ....126.2971V; APASS Data Release 9 (cat. II/336) and
* the SIMBAD database.
* If not available they were derived as described in Sect. 2.2.
* See Sect. 3 for derivation of L, Mass and Age.
*Note (2): Classification as follows:
* s = single-peaked
* d = double-peaked
* P = P-Cygni profile, both regular or inverse
*Note (3): in low-quality sample, variability indicator values (Vi) for these
* objects could not be derived as they are not astrometrically well behaved.
* Similarly, many of these sources fall outside the Pre-Main Sequence tracks
* and isochrones in the HR diagram and no masses or ages could be derived for
* them. We decided to present the masses and ages in the cases they were
* computable but these values have to be taken with caution.
*Note (4): Infrared excess at individual bandpasses is defined as Fobserved/FCK,
* being Fobserved the dereddened observed monochromatic flux and FCK the
* expected flux according to the Castelli-Kurucz model
* (Castelli & Kurucz, 2004, arXiv:astro-ph/0405087).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'lqsample.dat' ! Main parameters of each Herbig Ae/Be star belonging
to the low quality sample of 34 sources (table 2)
integer*4 nr__1
parameter (nr__1=34) ! Number of records
character*346 ar__1 ! Full-size record
C Position composed of: RAh RAm RAs DE- DEd DEm DEs (Epoch=2015.5)
character*13 Name_1 (nr__1) ! Name of the Herbig Ae/Be star
integer*4 RAh_1 (nr__1) ! (h) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 RAm_1 (nr__1) ! (min) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
real*4 RAs_1 (nr__1) ! (s) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
character*1 DE__1 (nr__1) ! Declination sign from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEd_1 (nr__1) ! (deg) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEs_1 (nr__1) ! (arcsec) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
real*8 plx_1 (nr__1) ! (mas) Parallax from Gaia DR2
real*8 e_plx_1 (nr__1) ! (mas) Parallax uncertainty from Gaia DR2
real*4 Dist_1 (nr__1) ! (pc) Distance derived using an exponentially
* decreasing density prior
real*4 E_Dist_2 (nr__1) ! (pc) Upper error on the Distance, to the 95th
* percentile confidence interval on the
* Distance estimate
real*4 e_Dist_3 (nr__1) ! (pc) Lower error on the Distance, to the 5th
* percentile confidence interval on the
* Distance estimate
integer*4 Teff_1 (nr__1) ! (K) Effective temperature from the literature (1)
integer*4 E_Teff_2 (nr__1) ! (K) Upper error of Teff
integer*4 e_Teff_3 (nr__1) ! (K) Lower error of Teff
real*4 LogL_1 (nr__1) ! ([Lsun]) Decimal logarithm of the luminosity
real*4 E_LogL_2 (nr__1) ! ([Lsun]) Upper error of LogL
real*4 e_LogL_3 (nr__1) ! ([Lsun]) Lower error of LogL
real*4 Av_1 (nr__1) ! (mag) Total extinction (1)
real*4 E_Av_2 (nr__1) ! (mag) Upper error of Av
real*4 e_Av_3 (nr__1) ! (mag) Lower error of Av
real*4 V_1 (nr__1) ! (mag) Johnson V band magnitude (1)
character*3 Bin_1 (nr__1) ! [Yes ] Binarity status: "Yes" if the source
* has been classified as binary,
* nothing otherwise
integer*4 r_Bin_1 (nr__1) ! ?=- References for Binary column
* (in ref.dat file)
real*4 E_NIR_1 (nr__1) ! ?=- Infrared excess in the 1.24-3.4um
* interval. Derived following Equation 2 of
* the paper
real*4 E_E_NIR_1 (nr__1) ! ?=- Upper error of Near_IR_excess
real*4 e_E_NIR_1 (nr__1) ! ?=- Lower error of Near_IR_excess
real*4 E_MIR_1 (nr__1) ! ?=- Infrared excess in the 3.4-22um
* interval. Derived following Equation 2 of
* the paper
real*4 E_E_MIR_1 (nr__1) ! ?=- Upper error of Mid_IR_excess
real*4 e_E_MIR_1 (nr__1) ! ?=- Lower error of Mid_IR_excess
real*8 EWHa_1 (nr__1) ! (0.1nm) ?=- Halpha equivalent width, from the
* literature. Negative numbers indicate
* emission
real*4 e_EWHa_1 (nr__1) ! (0.1nm) ?=- Uncertainty of EWHa, from the literature
integer*4 r_EWHa_1 (nr__1) ! ?=- References for EWHa and e_EWHa columns
* (in refs.dat file)
character*1 Hashape_1 (nr__1) ! Line profile of the Halpha line, from
* the literature (2)
integer*4 r_Hashape_1(nr__1) ! ?=- References for Hashape
* (in refs.dat file)
real*4 Vi_1 (nr__1) ! ?=- Variability indicator. Derived following
* Equation 4 of the paper (3)
character*3 UXOR_1 (nr__1) ! [Yes ] UXOR status: "Yes" if the source has
* been classified as UXOR type,
* nothing otherwise
real*4 Mass_1 (nr__1) ! (Msun) ?=- Mass of the source
real*4 E_Mass_2 (nr__1) ! (Msun) ?=- Upper error of the Mass
real*4 e_Mass_3 (nr__1) ! (Msun) ?=- Lower error of the Mass
real*4 Age_1 (nr__1) ! (Myr) ?=- Age of the source
real*4 E_Age_2 (nr__1) ! (Myr) ?=- Upper error of the Age
real*4 e_Age_3 (nr__1) ! (Myr) ?=- Lower error of the Age
real*4 E_J_1 (nr__1) ! ?=- Infrared excess at the J band (1.24um) (4)
real*4 E_H_1 (nr__1) ! ?=- Infrared excess at the H band (1.66um) (4)
real*8 E_Ks_1 (nr__1) ! ?=- Infrared excess at the
* Ks band (2.16um) (4)
real*8 E_W1_1 (nr__1) ! ?=- Infrared excess at the W1 band (3.4um) (4)
real*8 E_W2_1 (nr__1) ! ?=- Infrared excess at the W2 band (4.6um) (4)
real*8 E_W3_1 (nr__1) ! ?=- Infrared excess at the W3 band (12um) (4)
real*8 E_W4_1 (nr__1) ! ?=- Infrared excess at the W4 band (22um) (4)
character*3 The_1 (nr__1) ! [Yes ] Yes if the source appears in Table 1
* of The et al. (1994A&AS..104..315T),
* nothing otherwise
*Note (1): Atmospheric parameters Teff, AV and V taken from the following
* sources in order of choice:
* Fairlamb et al., 2015MNRAS.453..976F; Montesinos et al., 2009A&A...495..901M;
* Hernandez et al., 2004AJ....127.1682H; Mendigutia et al., 2012A&A...543A..59M;
* Carmona et al., 2010A&A...517A..67C; Chen et al., 2016NewA...44....1C;
* Alecian et al., 2013MNRAS.429.1001A; Sartori et al., 2010AJ....139...27S;
* Manoj et al., 2006ApJ...653..657M; Hernandez et al., 2005AJ....129..856H;
* Vieira et al., 2003AJ....126.2971V; APASS Data Release 9 (cat. II/336) and
* the SIMBAD database.
* If not available they were derived as described in Sect. 2.2.
* See Sect. 3 for derivation of L, Mass and Age.
*Note (2): Classification as follows:
* s = single-peaked
* d = double-peaked
* P = P-Cygni profile, both regular or inverse
*Note (3): in low-quality sample, variability indicator values (Vi) for these
* objects could not be derived as they are not astrometrically well behaved.
* Similarly, many of these sources fall outside the Pre-Main Sequence tracks
* and isochrones in the HR diagram and no masses or ages could be derived for
* them. We decided to present the masses and ages in the cases they were
* computable but these values have to be taken with caution.
*Note (4): Infrared excess at individual bandpasses is defined as Fobserved/FCK,
* being Fobserved the dereddened observed monochromatic flux and FCK the
* expected flux according to the Castelli-Kurucz model
* (Castelli & Kurucz, 2004, arXiv:astro-ph/0405087).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__2
parameter (nr__2=69) ! Number of records
character*132 ar__2 ! Full-size record
integer*4 Ref (nr__2) ! Reference number
character*19 BibCode (nr__2) ! BibCode
character*33 Aut (nr__2) ! Author's name
character*75 Com (nr__2) ! Comments
C=============================================================================
C Loading file 'hqsample.dat' ! Main parameters of each Herbig Ae/Be star belonging
* to the high quality sample of 218 sources (table 1)
C Format for file interpretation
1 format(
+ A13,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F9.6,1X,F8.6,
+ 1X,F6.1,1X,F5.1,1X,F5.1,1X,I5,1X,I4,1X,I4,1X,F5.2,1X,F4.2,1X,
+ F4.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.2,1X,A3,1X,I2,1X,E9.4,1X,
+ E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,F8.3,1X,F6.3,1X,I2,1X,
+ A1,1X,I2,1X,F5.2,1X,A3,1X,F6.3,1X,F6.3,1X,F6.3,1X,E9.4,1X,
+ E9.4,1X,E9.4,1X,F5.2,1X,F6.2,1X,F7.2,1X,F7.2,1X,F8.2,1X,F9.2,
+ 1X,F10.2,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'hqsample.dat')
write(6,*) '....Loading file: hqsample.dat'
do i__=1,218
read(1,'(A346)')ar__
read(ar__,1)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),plx(i__),e_plx(i__),Dist(i__),E_Dist(i__),
+ e_Dist_1(i__),Teff(i__),E_Teff(i__),e_Teff_1(i__),LogL(i__),
+ E_LogL(i__),e_LogL_1(i__),Av(i__),E_Av(i__),e_Av_1(i__),
+ V(i__),Bin(i__),r_Bin(i__),E_NIR(i__),E_E_NIR(i__),
+ e_E_NIR_1(i__),E_MIR(i__),E_E_MIR(i__),e_E_MIR_1(i__),
+ EWHa(i__),e_EWHa(i__),r_EWHa(i__),Hashape(i__),r_Hashape(i__),
+ Vi(i__),UXOR(i__),Mass(i__),E_Mass(i__),e_Mass_1(i__),
+ Age(i__),E_Age(i__),e_Age_1(i__),E_J(i__),E_H(i__),E_Ks(i__),
+ E_W1(i__),E_W2(i__),E_W3(i__),E_W4(i__),The(i__)
if (r_Bin(i__) .EQ. 45) r_Bin(i__) = iNULL__
if (idig(ar__(140:148)).EQ.0) E_NIR(i__) = rNULL__
if (idig(ar__(150:158)).EQ.0) E_E_NIR(i__) = rNULL__
if (idig(ar__(160:168)).EQ.0) e_E_NIR_1(i__) = rNULL__
if (idig(ar__(170:178)).EQ.0) E_MIR(i__) = rNULL__
if (idig(ar__(180:188)).EQ.0) E_E_MIR(i__) = rNULL__
if (idig(ar__(190:198)).EQ.0) e_E_MIR_1(i__) = rNULL__
if (idig(ar__(200:207)).EQ.0) EWHa(i__) = rNULL__
if (idig(ar__(209:214)).EQ.0) e_EWHa(i__) = rNULL__
if (r_EWHa(i__) .EQ. 45) r_EWHa(i__) = iNULL__
if (r_Hashape(i__) .EQ. 45) r_Hashape(i__) = iNULL__
if (idig(ar__(224:228)).EQ.0) Vi(i__) = rNULL__
if (idig(ar__(234:239)).EQ.0) Mass(i__) = rNULL__
if (idig(ar__(241:246)).EQ.0) E_Mass(i__) = rNULL__
if (idig(ar__(248:253)).EQ.0) e_Mass_1(i__) = rNULL__
if (idig(ar__(255:263)).EQ.0) Age(i__) = rNULL__
if (idig(ar__(265:273)).EQ.0) E_Age(i__) = rNULL__
if (idig(ar__(275:283)).EQ.0) e_Age_1(i__) = rNULL__
if (idig(ar__(285:289)).EQ.0) E_J(i__) = rNULL__
if (idig(ar__(291:296)).EQ.0) E_H(i__) = rNULL__
if (idig(ar__(298:304)).EQ.0) E_Ks(i__) = rNULL__
if (idig(ar__(306:312)).EQ.0) E_W1(i__) = rNULL__
if (idig(ar__(314:321)).EQ.0) E_W2(i__) = rNULL__
if (idig(ar__(323:331)).EQ.0) E_W3(i__) = rNULL__
if (idig(ar__(333:342)).EQ.0) E_W4(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),plx(i__),e_plx(i__),Dist(i__),E_Dist(i__),
+ e_Dist_1(i__),Teff(i__),E_Teff(i__),e_Teff_1(i__),LogL(i__),
+ E_LogL(i__),e_LogL_1(i__),Av(i__),E_Av(i__),e_Av_1(i__),
+ V(i__),Bin(i__),r_Bin(i__),E_NIR(i__),E_E_NIR(i__),
+ e_E_NIR_1(i__),E_MIR(i__),E_E_MIR(i__),e_E_MIR_1(i__),
+ EWHa(i__),e_EWHa(i__),r_EWHa(i__),Hashape(i__),r_Hashape(i__),
+ Vi(i__),UXOR(i__),Mass(i__),E_Mass(i__),e_Mass_1(i__),
+ Age(i__),E_Age(i__),e_Age_1(i__),E_J(i__),E_H(i__),E_Ks(i__),
+ E_W1(i__),E_W2(i__),E_W3(i__),E_W4(i__),The(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'lqsample.dat' ! Main parameters of each Herbig Ae/Be star belonging
* to the low quality sample of 34 sources (table 2)
C Format for file interpretation
2 format(
+ A13,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F9.6,1X,F8.6,
+ 1X,F6.1,1X,F5.1,1X,F5.1,1X,I5,1X,I4,1X,I4,1X,F5.2,1X,F4.2,1X,
+ F4.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.2,1X,A3,1X,I2,1X,E9.4,1X,
+ E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,F8.3,1X,F6.3,1X,I2,1X,
+ A1,1X,I2,1X,F5.2,1X,A3,1X,F6.3,1X,F6.3,1X,F6.3,1X,E9.4,1X,
+ E9.4,1X,E9.4,1X,F5.2,1X,F6.2,1X,F7.2,1X,F7.2,1X,F8.2,1X,F9.2,
+ 1X,F10.2,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'lqsample.dat')
write(6,*) '....Loading file: lqsample.dat'
do i__=1,34
read(1,'(A346)')ar__1
read(ar__1,2)
+ Name_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),plx_1(i__),e_plx_1(i__),
+ Dist_1(i__),E_Dist_2(i__),e_Dist_3(i__),Teff_1(i__),
+ E_Teff_2(i__),e_Teff_3(i__),LogL_1(i__),E_LogL_2(i__),
+ e_LogL_3(i__),Av_1(i__),E_Av_2(i__),e_Av_3(i__),V_1(i__),
+ Bin_1(i__),r_Bin_1(i__),E_NIR_1(i__),E_E_NIR_1(i__),
+ e_E_NIR_1(i__),E_MIR_1(i__),E_E_MIR_1(i__),e_E_MIR_1(i__),
+ EWHa_1(i__),e_EWHa_1(i__),r_EWHa_1(i__),Hashape_1(i__),
+ r_Hashape_1(i__),Vi_1(i__),UXOR_1(i__),Mass_1(i__),
+ E_Mass_2(i__),e_Mass_3(i__),Age_1(i__),E_Age_2(i__),
+ e_Age_3(i__),E_J_1(i__),E_H_1(i__),E_Ks_1(i__),E_W1_1(i__),
+ E_W2_1(i__),E_W3_1(i__),E_W4_1(i__),The_1(i__)
if (r_Bin_1(i__) .EQ. 45) r_Bin_1(i__) = iNULL__
if (idig(ar__1(140:148)).EQ.0) E_NIR_1(i__) = rNULL__
if (idig(ar__1(150:158)).EQ.0) E_E_NIR_1(i__) = rNULL__
if (idig(ar__1(160:168)).EQ.0) e_E_NIR_1(i__) = rNULL__
if (idig(ar__1(170:178)).EQ.0) E_MIR_1(i__) = rNULL__
if (idig(ar__1(180:188)).EQ.0) E_E_MIR_1(i__) = rNULL__
if (idig(ar__1(190:198)).EQ.0) e_E_MIR_1(i__) = rNULL__
if (idig(ar__1(200:207)).EQ.0) EWHa_1(i__) = rNULL__
if (idig(ar__1(209:214)).EQ.0) e_EWHa_1(i__) = rNULL__
if (r_EWHa_1(i__) .EQ. 45) r_EWHa_1(i__) = iNULL__
if (r_Hashape_1(i__) .EQ. 45) r_Hashape_1(i__) = iNULL__
if (idig(ar__1(224:228)).EQ.0) Vi_1(i__) = rNULL__
if (idig(ar__1(234:239)).EQ.0) Mass_1(i__) = rNULL__
if (idig(ar__1(241:246)).EQ.0) E_Mass_2(i__) = rNULL__
if (idig(ar__1(248:253)).EQ.0) e_Mass_3(i__) = rNULL__
if (idig(ar__1(255:263)).EQ.0) Age_1(i__) = rNULL__
if (idig(ar__1(265:273)).EQ.0) E_Age_2(i__) = rNULL__
if (idig(ar__1(275:283)).EQ.0) e_Age_3(i__) = rNULL__
if (idig(ar__1(285:289)).EQ.0) E_J_1(i__) = rNULL__
if (idig(ar__1(291:296)).EQ.0) E_H_1(i__) = rNULL__
if (idig(ar__1(298:304)).EQ.0) E_Ks_1(i__) = rNULL__
if (idig(ar__1(306:312)).EQ.0) E_W1_1(i__) = rNULL__
if (idig(ar__1(314:321)).EQ.0) E_W2_1(i__) = rNULL__
if (idig(ar__1(323:331)).EQ.0) E_W3_1(i__) = rNULL__
if (idig(ar__1(333:342)).EQ.0) E_W4_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),plx_1(i__),e_plx_1(i__),
+ Dist_1(i__),E_Dist_2(i__),e_Dist_3(i__),Teff_1(i__),
+ E_Teff_2(i__),e_Teff_3(i__),LogL_1(i__),E_LogL_2(i__),
+ e_LogL_3(i__),Av_1(i__),E_Av_2(i__),e_Av_3(i__),V_1(i__),
+ Bin_1(i__),r_Bin_1(i__),E_NIR_1(i__),E_E_NIR_1(i__),
+ e_E_NIR_1(i__),E_MIR_1(i__),E_E_MIR_1(i__),e_E_MIR_1(i__),
+ EWHa_1(i__),e_EWHa_1(i__),r_EWHa_1(i__),Hashape_1(i__),
+ r_Hashape_1(i__),Vi_1(i__),UXOR_1(i__),Mass_1(i__),
+ E_Mass_2(i__),e_Mass_3(i__),Age_1(i__),E_Age_2(i__),
+ e_Age_3(i__),E_J_1(i__),E_H_1(i__),E_Ks_1(i__),E_W1_1(i__),
+ E_W2_1(i__),E_W3_1(i__),E_W4_1(i__),The_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
3 format(I2,1X,A19,1X,A33,1X,A75)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,69
read(1,'(A132)')ar__2
read(ar__2,3)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,3)Ref(i__),BibCode(i__),Aut(i__),Com(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