Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/620/A128 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-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 ! Name of the Herbig Ae/Be star
integer*4 RAh ! (h) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 RAm ! (min) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
real*4 RAs ! (s) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
character*1 DE_ ! Declination sign from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEd ! (deg) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEm ! (arcmin) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEs ! (arcsec) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
real*8 plx ! (mas) Parallax from Gaia DR2
real*8 e_plx ! (mas) Parallax uncertainty from Gaia DR2
real*4 Dist ! (pc) Distance derived using an exponentially
* decreasing density prior
real*4 E_Dist ! (pc) Upper error on the Distance, to the 95th
* percentile confidence interval on the
* Distance estimate
real*4 e_Dist_1 ! (pc) Lower error on the Distance, to the 5th
* percentile confidence interval on the
* Distance estimate
integer*4 Teff ! (K) Effective temperature from the literature (1)
integer*4 E_Teff ! (K) Upper error of Teff
integer*4 e_Teff_1 ! (K) Lower error of Teff
real*4 LogL ! ([Lsun]) Decimal logarithm of the luminosity
real*4 E_LogL ! ([Lsun]) Upper error of LogL
real*4 e_LogL_1 ! ([Lsun]) Lower error of LogL
real*4 Av ! (mag) Total extinction (1)
real*4 E_Av ! (mag) Upper error of Av
real*4 e_Av_1 ! (mag) Lower error of Av
real*4 V ! (mag) Johnson V band magnitude (1)
character*3 Bin ! [Yes ] Binarity status: "Yes" if the source
* has been classified as binary,
* nothing otherwise
integer*4 r_Bin ! ?=- References for Binary column
* (in ref.dat file)
real*4 E_NIR ! ?=- Infrared excess in the 1.24-3.4um
* interval. Derived following Equation 2 of
* the paper
real*4 E_E_NIR ! ?=- Upper error of Near_IR_excess
real*4 e_E_NIR_1 ! ?=- Lower error of Near_IR_excess
real*4 E_MIR ! ?=- Infrared excess in the 3.4-22um
* interval. Derived following Equation 2 of
* the paper
real*4 E_E_MIR ! ?=- Upper error of Mid_IR_excess
real*4 e_E_MIR_1 ! ?=- Lower error of Mid_IR_excess
real*8 EWHa ! (0.1nm) ?=- Halpha equivalent width, from the
* literature. Negative numbers indicate
* emission
real*4 e_EWHa ! (0.1nm) ?=- Uncertainty of EWHa, from the literature
integer*4 r_EWHa ! ?=- References for EWHa and e_EWHa columns
* (in refs.dat file)
character*1 Hashape ! Line profile of the Halpha line, from
* the literature (2)
integer*4 r_Hashape ! ?=- References for Hashape
* (in refs.dat file)
real*4 Vi ! ?=- Variability indicator. Derived following
* Equation 4 of the paper (3)
character*3 UXOR ! [Yes ] UXOR status: "Yes" if the source has
* been classified as UXOR type,
* nothing otherwise
real*4 Mass ! (Msun) ?=- Mass of the source
real*4 E_Mass ! (Msun) ?=- Upper error of the Mass
real*4 e_Mass_1 ! (Msun) ?=- Lower error of the Mass
real*4 Age ! (Myr) ?=- Age of the source
real*4 E_Age ! (Myr) ?=- Upper error of the Age
real*4 e_Age_1 ! (Myr) ?=- Lower error of the Age
real*4 E_J ! ?=- Infrared excess at the J band (1.24um) (4)
real*4 E_H ! ?=- Infrared excess at the H band (1.66um) (4)
real*8 E_Ks ! ?=- Infrared excess at the
* Ks band (2.16um) (4)
real*8 E_W1 ! ?=- Infrared excess at the W1 band (3.4um) (4)
real*8 E_W2 ! ?=- Infrared excess at the W2 band (4.6um) (4)
real*8 E_W3 ! ?=- Infrared excess at the W3 band (12um) (4)
real*8 E_W4 ! ?=- Infrared excess at the W4 band (22um) (4)
character*3 The ! [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 ! Name of the Herbig Ae/Be star
integer*4 RAh_1 ! (h) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 RAm_1 ! (min) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
real*4 RAs_1 ! (s) Right ascension from Gaia DR2
* (ICRS at epoch 2015.5)
character*1 DE__1 ! Declination sign from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEd_1 ! (deg) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEm_1 ! (arcmin) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
integer*4 DEs_1 ! (arcsec) Declination from Gaia DR2
* (ICRS at epoch 2015.5)
real*8 plx_1 ! (mas) Parallax from Gaia DR2
real*8 e_plx_1 ! (mas) Parallax uncertainty from Gaia DR2
real*4 Dist_1 ! (pc) Distance derived using an exponentially
* decreasing density prior
real*4 E_Dist_2 ! (pc) Upper error on the Distance, to the 95th
* percentile confidence interval on the
* Distance estimate
real*4 e_Dist_3 ! (pc) Lower error on the Distance, to the 5th
* percentile confidence interval on the
* Distance estimate
integer*4 Teff_1 ! (K) Effective temperature from the literature (1)
integer*4 E_Teff_2 ! (K) Upper error of Teff
integer*4 e_Teff_3 ! (K) Lower error of Teff
real*4 LogL_1 ! ([Lsun]) Decimal logarithm of the luminosity
real*4 E_LogL_2 ! ([Lsun]) Upper error of LogL
real*4 e_LogL_3 ! ([Lsun]) Lower error of LogL
real*4 Av_1 ! (mag) Total extinction (1)
real*4 E_Av_2 ! (mag) Upper error of Av
real*4 e_Av_3 ! (mag) Lower error of Av
real*4 V_1 ! (mag) Johnson V band magnitude (1)
character*3 Bin_1 ! [Yes ] Binarity status: "Yes" if the source
* has been classified as binary,
* nothing otherwise
integer*4 r_Bin_1 ! ?=- References for Binary column
* (in ref.dat file)
real*4 E_NIR_1 ! ?=- Infrared excess in the 1.24-3.4um
* interval. Derived following Equation 2 of
* the paper
real*4 E_E_NIR_1 ! ?=- Upper error of Near_IR_excess
real*4 e_E_NIR_1 ! ?=- Lower error of Near_IR_excess
real*4 E_MIR_1 ! ?=- Infrared excess in the 3.4-22um
* interval. Derived following Equation 2 of
* the paper
real*4 E_E_MIR_1 ! ?=- Upper error of Mid_IR_excess
real*4 e_E_MIR_1 ! ?=- Lower error of Mid_IR_excess
real*8 EWHa_1 ! (0.1nm) ?=- Halpha equivalent width, from the
* literature. Negative numbers indicate
* emission
real*4 e_EWHa_1 ! (0.1nm) ?=- Uncertainty of EWHa, from the literature
integer*4 r_EWHa_1 ! ?=- References for EWHa and e_EWHa columns
* (in refs.dat file)
character*1 Hashape_1 ! Line profile of the Halpha line, from
* the literature (2)
integer*4 r_Hashape_1 ! ?=- References for Hashape
* (in refs.dat file)
real*4 Vi_1 ! ?=- Variability indicator. Derived following
* Equation 4 of the paper (3)
character*3 UXOR_1 ! [Yes ] UXOR status: "Yes" if the source has
* been classified as UXOR type,
* nothing otherwise
real*4 Mass_1 ! (Msun) ?=- Mass of the source
real*4 E_Mass_2 ! (Msun) ?=- Upper error of the Mass
real*4 e_Mass_3 ! (Msun) ?=- Lower error of the Mass
real*4 Age_1 ! (Myr) ?=- Age of the source
real*4 E_Age_2 ! (Myr) ?=- Upper error of the Age
real*4 e_Age_3 ! (Myr) ?=- Lower error of the Age
real*4 E_J_1 ! ?=- Infrared excess at the J band (1.24um) (4)
real*4 E_H_1 ! ?=- Infrared excess at the H band (1.66um) (4)
real*8 E_Ks_1 ! ?=- Infrared excess at the
* Ks band (2.16um) (4)
real*8 E_W1_1 ! ?=- Infrared excess at the W1 band (3.4um) (4)
real*8 E_W2_1 ! ?=- Infrared excess at the W2 band (4.6um) (4)
real*8 E_W3_1 ! ?=- Infrared excess at the W3 band (12um) (4)
real*8 E_W4_1 ! ?=- Infrared excess at the W4 band (22um) (4)
character*3 The_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 ! Reference number
character*19 BibCode ! BibCode
character*33 Aut ! Author's name
character*75 Com ! 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,RAh,RAm,RAs,DE_,DEd,DEm,DEs,plx,e_plx,Dist,E_Dist,
+ e_Dist_1,Teff,E_Teff,e_Teff_1,LogL,E_LogL,e_LogL_1,Av,E_Av,
+ e_Av_1,V,Bin,r_Bin,E_NIR,E_E_NIR,e_E_NIR_1,E_MIR,E_E_MIR,
+ e_E_MIR_1,EWHa,e_EWHa,r_EWHa,Hashape,r_Hashape,Vi,UXOR,Mass,
+ E_Mass,e_Mass_1,Age,E_Age,e_Age_1,E_J,E_H,E_Ks,E_W1,E_W2,E_W3,
+ E_W4,The
if (r_Bin .EQ. 45) r_Bin = iNULL__
if (idig(ar__(140:148)).EQ.0) E_NIR = rNULL__
if (idig(ar__(150:158)).EQ.0) E_E_NIR = rNULL__
if (idig(ar__(160:168)).EQ.0) e_E_NIR_1 = rNULL__
if (idig(ar__(170:178)).EQ.0) E_MIR = rNULL__
if (idig(ar__(180:188)).EQ.0) E_E_MIR = rNULL__
if (idig(ar__(190:198)).EQ.0) e_E_MIR_1 = rNULL__
if (idig(ar__(200:207)).EQ.0) EWHa = rNULL__
if (idig(ar__(209:214)).EQ.0) e_EWHa = rNULL__
if (r_EWHa .EQ. 45) r_EWHa = iNULL__
if (r_Hashape .EQ. 45) r_Hashape = iNULL__
if (idig(ar__(224:228)).EQ.0) Vi = rNULL__
if (idig(ar__(234:239)).EQ.0) Mass = rNULL__
if (idig(ar__(241:246)).EQ.0) E_Mass = rNULL__
if (idig(ar__(248:253)).EQ.0) e_Mass_1 = rNULL__
if (idig(ar__(255:263)).EQ.0) Age = rNULL__
if (idig(ar__(265:273)).EQ.0) E_Age = rNULL__
if (idig(ar__(275:283)).EQ.0) e_Age_1 = rNULL__
if (idig(ar__(285:289)).EQ.0) E_J = rNULL__
if (idig(ar__(291:296)).EQ.0) E_H = rNULL__
if (idig(ar__(298:304)).EQ.0) E_Ks = rNULL__
if (idig(ar__(306:312)).EQ.0) E_W1 = rNULL__
if (idig(ar__(314:321)).EQ.0) E_W2 = rNULL__
if (idig(ar__(323:331)).EQ.0) E_W3 = rNULL__
if (idig(ar__(333:342)).EQ.0) E_W4 = rNULL__
c ..............Just test output...........
write(6,1)
+ Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,plx,e_plx,Dist,E_Dist,
+ e_Dist_1,Teff,E_Teff,e_Teff_1,LogL,E_LogL,e_LogL_1,Av,E_Av,
+ e_Av_1,V,Bin,r_Bin,E_NIR,E_E_NIR,e_E_NIR_1,E_MIR,E_E_MIR,
+ e_E_MIR_1,EWHa,e_EWHa,r_EWHa,Hashape,r_Hashape,Vi,UXOR,Mass,
+ E_Mass,e_Mass_1,Age,E_Age,e_Age_1,E_J,E_H,E_Ks,E_W1,E_W2,E_W3,
+ E_W4,The
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,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,plx_1,
+ e_plx_1,Dist_1,E_Dist_2,e_Dist_3,Teff_1,E_Teff_2,e_Teff_3,
+ LogL_1,E_LogL_2,e_LogL_3,Av_1,E_Av_2,e_Av_3,V_1,Bin_1,r_Bin_1,
+ E_NIR_1,E_E_NIR_1,e_E_NIR_1,E_MIR_1,E_E_MIR_1,e_E_MIR_1,
+ EWHa_1,e_EWHa_1,r_EWHa_1,Hashape_1,r_Hashape_1,Vi_1,UXOR_1,
+ Mass_1,E_Mass_2,e_Mass_3,Age_1,E_Age_2,e_Age_3,E_J_1,E_H_1,
+ E_Ks_1,E_W1_1,E_W2_1,E_W3_1,E_W4_1,The_1
if (r_Bin_1 .EQ. 45) r_Bin_1 = iNULL__
if (idig(ar__1(140:148)).EQ.0) E_NIR_1 = rNULL__
if (idig(ar__1(150:158)).EQ.0) E_E_NIR_1 = rNULL__
if (idig(ar__1(160:168)).EQ.0) e_E_NIR_1 = rNULL__
if (idig(ar__1(170:178)).EQ.0) E_MIR_1 = rNULL__
if (idig(ar__1(180:188)).EQ.0) E_E_MIR_1 = rNULL__
if (idig(ar__1(190:198)).EQ.0) e_E_MIR_1 = rNULL__
if (idig(ar__1(200:207)).EQ.0) EWHa_1 = rNULL__
if (idig(ar__1(209:214)).EQ.0) e_EWHa_1 = rNULL__
if (r_EWHa_1 .EQ. 45) r_EWHa_1 = iNULL__
if (r_Hashape_1 .EQ. 45) r_Hashape_1 = iNULL__
if (idig(ar__1(224:228)).EQ.0) Vi_1 = rNULL__
if (idig(ar__1(234:239)).EQ.0) Mass_1 = rNULL__
if (idig(ar__1(241:246)).EQ.0) E_Mass_2 = rNULL__
if (idig(ar__1(248:253)).EQ.0) e_Mass_3 = rNULL__
if (idig(ar__1(255:263)).EQ.0) Age_1 = rNULL__
if (idig(ar__1(265:273)).EQ.0) E_Age_2 = rNULL__
if (idig(ar__1(275:283)).EQ.0) e_Age_3 = rNULL__
if (idig(ar__1(285:289)).EQ.0) E_J_1 = rNULL__
if (idig(ar__1(291:296)).EQ.0) E_H_1 = rNULL__
if (idig(ar__1(298:304)).EQ.0) E_Ks_1 = rNULL__
if (idig(ar__1(306:312)).EQ.0) E_W1_1 = rNULL__
if (idig(ar__1(314:321)).EQ.0) E_W2_1 = rNULL__
if (idig(ar__1(323:331)).EQ.0) E_W3_1 = rNULL__
if (idig(ar__1(333:342)).EQ.0) E_W4_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,plx_1,
+ e_plx_1,Dist_1,E_Dist_2,e_Dist_3,Teff_1,E_Teff_2,e_Teff_3,
+ LogL_1,E_LogL_2,e_LogL_3,Av_1,E_Av_2,e_Av_3,V_1,Bin_1,r_Bin_1,
+ E_NIR_1,E_E_NIR_1,e_E_NIR_1,E_MIR_1,E_E_MIR_1,e_E_MIR_1,
+ EWHa_1,e_EWHa_1,r_EWHa_1,Hashape_1,r_Hashape_1,Vi_1,UXOR_1,
+ Mass_1,E_Mass_2,e_Mass_3,Age_1,E_Age_2,e_Age_3,E_J_1,E_H_1,
+ E_Ks_1,E_W1_1,E_W2_1,E_W3_1,E_W4_1,The_1
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,BibCode,Aut,Com
c ..............Just test output...........
write(6,3)Ref,BibCode,Aut,Com
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