Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/522/A79 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-19
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/522/A79 Faint, high-Galactic-latitude red clump stars (Valentini+, 2010)
*================================================================================
*A spectroscopic survey of faint, high-galactic-latitude red clump stars.
*I. The high resolution sample.
* Valentini M., Munari U.
* <Astron. Astrophys., 522, A79-79 (2010)>
* =2010A&A...522A..79V
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'arcs_a.dat' ! Parameters of 207 faint red clump stars
integer*4 nr__
parameter (nr__=273) ! Number of records
character*406 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 HD (nr__) ! ? HD number
integer*4 HIP (nr__) ! ? HIP (Cat. I/239) number
integer*4 TYC1 (nr__) ! TYCHO-2 1st identifier (Cat. I/239)
integer*4 TYC2 (nr__) ! TYCHO-2 2nd identifier (Cat. I/239)
integer*4 TYC3 (nr__) ! TYCHO-2 3rd identifier (Cat. I/239)
character*11 SpType (nr__) ! MK spectral type from Michigan catalog (1)
real*8 RAdeg (nr__) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
real*8 Glon (nr__) ! (deg) Galactic longitude
real*8 Glat (nr__) ! (deg) Galactic latitude
real*4 Hpmag (nr__) ! (mag) ? Hipparcos Hp magnitude
real*4 e_Hpmag (nr__) ! (mag) ? Error on Hpmag
real*4 V_I (nr__) ! (mag) ? Cousins V-I from Hipparcos catalog
real*4 e_V_I (nr__) ! (mag) ? Error on V-I
real*8 plx1 (nr__) ! (mas) ? Hipparcos (ESA 1997, I/239) parallax
real*4 e_plx1 (nr__) ! (mas) ? Error on plx1
integer*4 dmax (nr__) ! (pc) ? Maximum distance from Hipparcos
* (ESA 1997, I/239) parallax
integer*4 dmin (nr__) ! (pc) ? Minimum distance from Hipparcos
* (ESA 1997, I/239) parallax
real*8 plx2 (nr__) ! (mas) ? Revised Hipparcos (van Leeuwen 2007, I/311)
* parallax
real*4 e_plx2 (nr__) ! (mas) ? Error on plx2
integer*4 dmax2 (nr__) ! (pc) ? Distance from revised Hipparcos parallax
integer*4 dmin2 (nr__) ! (pc) ? Distance from revised Hipparcos parallax
character*1 n_dmin2 (nr__) ! [AED]
real*4 BTmag (nr__) ! (mag) Tycho BT magnitude
real*4 e_BTmag (nr__) ! (mag) Error on Tycho BT magnitude
real*4 VTmag (nr__) ! (mag) Tycho VT magnitude
real*4 e_VTmag (nr__) ! (mag) Error on Tycho VT magnitude
real*4 pmRA (nr__) ! (mas/yr) ? Tycho-2 (I/259) RA proper motion
real*4 e_pmRA (nr__) ! (mas/yr) ? Error on Tycho-2 RA proper motion
real*4 pmDE (nr__) ! (mas/yr) ? Tycho-2 (I/259) DE proper motion
real*4 e_pmDE (nr__) ! (mas/yr) ? Error on Tycho-2 DE proper motion
real*4 Jmag (nr__) ! (mag) 2MASS J magnitude
real*4 e_Jmag (nr__) ! (mag) Error on 2MASS J magnitude
real*4 Hmag (nr__) ! (mag) 2MASS H magnitude
real*4 e_Hmag (nr__) ! (mag) Error on 2MASS H magnitude
real*4 Kmag (nr__) ! (mag) 2MASS K magnitude
real*4 e_Kmag (nr__) ! (mag) Error on 2MASS K magnitude
character*3 q2M (nr__) ! 2MASS quality index
real*4 Imag (nr__) ! (mag) ? DENIS I magnitude
real*4 e_Imag (nr__) ! (mag) ? Error on IDENIS
integer*4 q_Imag (nr__) ! ? DENIS quality index
real*8 HJD (nr__) ! Heliocentric Julian date of Observation
real*4 HRV (nr__) ! (km/s) Heliocentric radial velocity
real*4 e_HRV (nr__) ! (km/s) Error on heliocentric radial velocity
real*4 TCor (nr__) ! (km/s) Radial velocity of telluric absorption
* spectrum
integer*4 Teff (nr__) ! (K) Effective temperature
integer*4 e_Teff (nr__) ! (K) Error on effective temperature
real*4 logg (nr__) ! ([cm/s2]) Surface gravity
real*4 e_logg (nr__) ! ([cm/s2]) Error on surface gravity
real*4 v_M_H_ (nr__) ! ([Sun]) Metallicity
real*4 e__M_H_ (nr__) ! ([Sun]) Error on metallicity
real*4 Vrot (nr__) ! (km/s) ? Rotational velocity, Vrot*sini
real*4 e_Vrot (nr__) ! (km/s) ? Error on Vrot*sini
integer*4 d (nr__) ! (pc) ? Spectro-photometric distance
integer*4 e_d (nr__) ! (pc) ? Error on spectro-photometric distance
real*8 E_B_V (nr__) ! (mag) Reddening from Arenou et al.
* (1992A&A...258..104A)
integer*4 X (nr__) ! (pc) ? Galactic X coordinate
integer*4 Y (nr__) ! (pc) ? Galactic Y coordinate
integer*4 Z (nr__) ! (pc) ? Galactic Z coordinate
integer*4 Uvel (nr__) ! (km/s) ? U velocity
integer*4 Vvel (nr__) ! (km/s) ? V velocity
integer*4 Wvel (nr__) ! (km/s) ? W velocity
character*7 Com (nr__) ! Comments (2)
*Note (1): Houk et al. 1975-1999, see Catalogs
* III/31, III/51, III/80, III/133, III/214
*Note (2): Comments as follows:
* B = binary candidate
* lcII = luminosity class II
* Praes = in Praesepe cluster
* Coma = in Coma cluster
* Ta08 = in Takeda et al. (2008PASJ...60..781T)
* HM07 = in Hekker & Melendez (2007, Cat. J/A+A/475/1003)
* So05 = in Soubiran et al. (2005, Cat. J/A+A/43/8139)
* RaveDR2 = from Rave DR2
* Lit = from literature
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'arcs_b.dat' ! Parameters of 70 brighter red clump stars
integer*4 nr__1
parameter (nr__1=103) ! Number of records
character*406 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 HD_1 (nr__1) ! ? HD number
integer*4 HIP_1 (nr__1) ! ? HIP (Cat. I/239) number
integer*4 TYC1_1 (nr__1) ! TYCHO-2 1st identifier (Cat. I/239)
integer*4 TYC2_1 (nr__1) ! TYCHO-2 2nd identifier (Cat. I/239)
integer*4 TYC3_1 (nr__1) ! TYCHO-2 3rd identifier (Cat. I/239)
character*11 SpType_1 (nr__1) ! MK spectral type from Michigan catalog (1)
real*8 RAdeg_1 (nr__1) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Declination in decimal degrees (J2000)
real*8 Glon_1 (nr__1) ! (deg) Galactic longitude
real*8 Glat_1 (nr__1) ! (deg) Galactic latitude
real*4 Hpmag_1 (nr__1) ! (mag) ? Hipparcos Hp magnitude
real*4 e_Hpmag_1 (nr__1) ! (mag) ? Error on Hpmag
real*4 V_I_1 (nr__1) ! (mag) ? Cousins V-I from Hipparcos catalog
real*4 e_V_I_1 (nr__1) ! (mag) ? Error on V-I
real*8 plx1_1 (nr__1) ! (mas) ? Hipparcos (ESA 1997, I/239) parallax
real*4 e_plx1_1 (nr__1) ! (mas) ? Error on plx1
integer*4 dmax_1 (nr__1) ! (pc) ? Maximum distance from Hipparcos
* (ESA 1997, I/239) parallax
integer*4 dmin_1 (nr__1) ! (pc) ? Minimum distance from Hipparcos
* (ESA 1997, I/239) parallax
real*8 plx2_1 (nr__1) ! (mas) ? Revised Hipparcos (van Leeuwen 2007, I/311)
* parallax
real*4 e_plx2_1 (nr__1) ! (mas) ? Error on plx2
integer*4 dmax2_1 (nr__1) ! (pc) ? Distance from revised Hipparcos parallax
integer*4 dmin2_1 (nr__1) ! (pc) ? Distance from revised Hipparcos parallax
character*1 n_dmin2_1 (nr__1) ! [AED]
real*4 BTmag_1 (nr__1) ! (mag) Tycho BT magnitude
real*4 e_BTmag_1 (nr__1) ! (mag) Error on Tycho BT magnitude
real*4 VTmag_1 (nr__1) ! (mag) Tycho VT magnitude
real*4 e_VTmag_1 (nr__1) ! (mag) Error on Tycho VT magnitude
real*4 pmRA_1 (nr__1) ! (mas/yr) ? Tycho-2 (I/259) RA proper motion
real*4 e_pmRA_1 (nr__1) ! (mas/yr) ? Error on Tycho-2 RA proper motion
real*4 pmDE_1 (nr__1) ! (mas/yr) ? Tycho-2 (I/259) DE proper motion
real*4 e_pmDE_1 (nr__1) ! (mas/yr) ? Error on Tycho-2 DE proper motion
real*4 Jmag_1 (nr__1) ! (mag) 2MASS J magnitude
real*4 e_Jmag_1 (nr__1) ! (mag) Error on 2MASS J magnitude
real*4 Hmag_1 (nr__1) ! (mag) 2MASS H magnitude
real*4 e_Hmag_1 (nr__1) ! (mag) Error on 2MASS H magnitude
real*4 Kmag_1 (nr__1) ! (mag) 2MASS K magnitude
real*4 e_Kmag_1 (nr__1) ! (mag) Error on 2MASS K magnitude
character*3 q2M_1 (nr__1) ! 2MASS quality index
real*4 Imag_1 (nr__1) ! (mag) ? DENIS I magnitude
real*4 e_Imag_1 (nr__1) ! (mag) ? Error on IDENIS
integer*4 q_Imag_1 (nr__1) ! ? DENIS quality index
real*8 HJD_1 (nr__1) ! Heliocentric Julian date of Observation
real*4 HRV_1 (nr__1) ! (km/s) Heliocentric radial velocity
real*4 e_HRV_1 (nr__1) ! (km/s) Error on heliocentric radial velocity
real*4 TCor_1 (nr__1) ! (km/s) Radial velocity of telluric absorption
* spectrum
integer*4 Teff_1 (nr__1) ! (K) Effective temperature
integer*4 e_Teff_1 (nr__1) ! (K) Error on effective temperature
real*4 logg_1 (nr__1) ! ([cm/s2]) Surface gravity
real*4 e_logg_1 (nr__1) ! ([cm/s2]) Error on surface gravity
real*4 v_M_H__1 (nr__1) ! ([Sun]) Metallicity
real*4 e__M_H__1 (nr__1) ! ([Sun]) Error on metallicity
real*4 Vrot_1 (nr__1) ! (km/s) ? Rotational velocity, Vrot*sini
real*4 e_Vrot_1 (nr__1) ! (km/s) ? Error on Vrot*sini
integer*4 d_1 (nr__1) ! (pc) ? Spectro-photometric distance
integer*4 e_d_1 (nr__1) ! (pc) ? Error on spectro-photometric distance
real*8 E_B_V_1 (nr__1) ! (mag) Reddening from Arenou et al.
* (1992A&A...258..104A)
integer*4 X_1 (nr__1) ! (pc) ? Galactic X coordinate
integer*4 Y_1 (nr__1) ! (pc) ? Galactic Y coordinate
integer*4 Z_1 (nr__1) ! (pc) ? Galactic Z coordinate
integer*4 Uvel_1 (nr__1) ! (km/s) ? U velocity
integer*4 Vvel_1 (nr__1) ! (km/s) ? V velocity
integer*4 Wvel_1 (nr__1) ! (km/s) ? W velocity
character*7 Com_1 (nr__1) ! Comments (2)
*Note (1): Houk et al. 1975-1999, see Catalogs
* III/31, III/51, III/80, III/133, III/214
*Note (2): Comments as follows:
* B = binary candidate
* lcII = luminosity class II
* Praes = in Praesepe cluster
* Coma = in Coma cluster
* Ta08 = in Takeda et al. (2008PASJ...60..781T)
* HM07 = in Hekker & Melendez (2007, Cat. J/A+A/475/1003)
* So05 = in Soubiran et al. (2005, Cat. J/A+A/43/8139)
* RaveDR2 = from Rave DR2
* Lit = from literature
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'arcs_c.dat' ! Parameters of 55 test stars
integer*4 nr__2
parameter (nr__2=55) ! Number of records
character*406 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 HD_2 (nr__2) ! ? HD number
integer*4 HIP_2 (nr__2) ! ? HIP (Cat. I/239) number
integer*4 TYC1_2 (nr__2) ! TYCHO-2 1st identifier (Cat. I/239)
integer*4 TYC2_2 (nr__2) ! TYCHO-2 2nd identifier (Cat. I/239)
integer*4 TYC3_2 (nr__2) ! TYCHO-2 3rd identifier (Cat. I/239)
character*11 SpType_2 (nr__2) ! MK spectral type from Michigan catalog (1)
real*8 RAdeg_2 (nr__2) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Declination in decimal degrees (J2000)
real*8 Glon_2 (nr__2) ! (deg) Galactic longitude
real*8 Glat_2 (nr__2) ! (deg) Galactic latitude
real*4 Hpmag_2 (nr__2) ! (mag) ? Hipparcos Hp magnitude
real*4 e_Hpmag_2 (nr__2) ! (mag) ? Error on Hpmag
real*4 V_I_2 (nr__2) ! (mag) ? Cousins V-I from Hipparcos catalog
real*4 e_V_I_2 (nr__2) ! (mag) ? Error on V-I
real*8 plx1_2 (nr__2) ! (mas) ? Hipparcos (ESA 1997, I/239) parallax
real*4 e_plx1_2 (nr__2) ! (mas) ? Error on plx1
integer*4 dmax_2 (nr__2) ! (pc) ? Maximum distance from Hipparcos
* (ESA 1997, I/239) parallax
integer*4 dmin_2 (nr__2) ! (pc) ? Minimum distance from Hipparcos
* (ESA 1997, I/239) parallax
real*8 plx2_2 (nr__2) ! (mas) ? Revised Hipparcos (van Leeuwen 2007, I/311)
* parallax
real*4 e_plx2_2 (nr__2) ! (mas) ? Error on plx2
integer*4 dmax2_2 (nr__2) ! (pc) ? Distance from revised Hipparcos parallax
integer*4 dmin2_2 (nr__2) ! (pc) ? Distance from revised Hipparcos parallax
character*1 n_dmin2_2 (nr__2) ! [AED]
real*4 BTmag_2 (nr__2) ! (mag) Tycho BT magnitude
real*4 e_BTmag_2 (nr__2) ! (mag) Error on Tycho BT magnitude
real*4 VTmag_2 (nr__2) ! (mag) Tycho VT magnitude
real*4 e_VTmag_2 (nr__2) ! (mag) Error on Tycho VT magnitude
real*4 pmRA_2 (nr__2) ! (mas/yr) ? Tycho-2 (I/259) RA proper motion
real*4 e_pmRA_2 (nr__2) ! (mas/yr) ? Error on Tycho-2 RA proper motion
real*4 pmDE_2 (nr__2) ! (mas/yr) ? Tycho-2 (I/259) DE proper motion
real*4 e_pmDE_2 (nr__2) ! (mas/yr) ? Error on Tycho-2 DE proper motion
real*4 Jmag_2 (nr__2) ! (mag) 2MASS J magnitude
real*4 e_Jmag_2 (nr__2) ! (mag) Error on 2MASS J magnitude
real*4 Hmag_2 (nr__2) ! (mag) 2MASS H magnitude
real*4 e_Hmag_2 (nr__2) ! (mag) Error on 2MASS H magnitude
real*4 Kmag_2 (nr__2) ! (mag) 2MASS K magnitude
real*4 e_Kmag_2 (nr__2) ! (mag) Error on 2MASS K magnitude
character*3 q2M_2 (nr__2) ! 2MASS quality index
real*4 Imag_2 (nr__2) ! (mag) ? DENIS I magnitude
real*4 e_Imag_2 (nr__2) ! (mag) ? Error on IDENIS
integer*4 q_Imag_2 (nr__2) ! ? DENIS quality index
real*8 HJD_2 (nr__2) ! Heliocentric Julian date of Observation
real*4 HRV_2 (nr__2) ! (km/s) Heliocentric radial velocity
real*4 e_HRV_2 (nr__2) ! (km/s) Error on heliocentric radial velocity
real*4 TCor_2 (nr__2) ! (km/s) Radial velocity of telluric absorption
* spectrum
integer*4 Teff_2 (nr__2) ! (K) Effective temperature
integer*4 e_Teff_2 (nr__2) ! (K) Error on effective temperature
real*4 logg_2 (nr__2) ! ([cm/s2]) Surface gravity
real*4 e_logg_2 (nr__2) ! ([cm/s2]) Error on surface gravity
real*4 v_M_H__2 (nr__2) ! ([Sun]) Metallicity
real*4 e__M_H__2 (nr__2) ! ([Sun]) Error on metallicity
real*4 Vrot_2 (nr__2) ! (km/s) ? Rotational velocity, Vrot*sini
real*4 e_Vrot_2 (nr__2) ! (km/s) ? Error on Vrot*sini
integer*4 d_2 (nr__2) ! (pc) ? Spectro-photometric distance
integer*4 e_d_2 (nr__2) ! (pc) ? Error on spectro-photometric distance
real*8 E_B_V_1 (nr__2) ! (mag) Reddening from Arenou et al.
* (1992A&A...258..104A)
integer*4 X_2 (nr__2) ! (pc) ? Galactic X coordinate
integer*4 Y_2 (nr__2) ! (pc) ? Galactic Y coordinate
integer*4 Z_2 (nr__2) ! (pc) ? Galactic Z coordinate
integer*4 Uvel_2 (nr__2) ! (km/s) ? U velocity
integer*4 Vvel_2 (nr__2) ! (km/s) ? V velocity
integer*4 Wvel_2 (nr__2) ! (km/s) ? W velocity
character*7 Com_2 (nr__2) ! Comments (2)
*Note (1): Houk et al. 1975-1999, see Catalogs
* III/31, III/51, III/80, III/133, III/214
*Note (2): Comments as follows:
* B = binary candidate
* lcII = luminosity class II
* Praes = in Praesepe cluster
* Coma = in Coma cluster
* Ta08 = in Takeda et al. (2008PASJ...60..781T)
* HM07 = in Hekker & Melendez (2007, Cat. J/A+A/475/1003)
* So05 = in Soubiran et al. (2005, Cat. J/A+A/43/8139)
* RaveDR2 = from Rave DR2
* Lit = from literature
C=============================================================================
C Loading file 'arcs_a.dat' ! Parameters of 207 faint red clump stars
C Format for file interpretation
1 format(
+ I6,2X,I6,2X,I5,1X,I5,1X,I1,1X,A11,1X,F12.8,1X,F12.8,1X,F8.4,
+ 1X,F8.4,1X,F6.4,1X,F6.4,1X,F4.2,1X,F4.2,1X,F7.2,1X,F6.2,1X,I5,
+ 1X,I5,1X,F7.2,1X,F6.2,1X,I6,1X,I6,A1,F6.3,1X,F5.3,1X,F6.3,1X,
+ F5.3,1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.3,1X,F5.3,1X,F6.3,
+ 1X,F5.3,1X,F6.3,1X,F5.3,1X,A3,1X,F6.3,1X,F4.2,1X,I3,1X,F11.3,
+ 1X,F6.1,1X,F3.1,1X,F4.2,1X,I4,1X,I2,1X,F4.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F5.1,1X,F4.1,1X,I5,1X,I3,1X,F9.7,1X,I5,1X,I5,1X,I5,
+ 1X,I5,1X,I5,1X,I5,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'arcs_a.dat')
write(6,*) '....Loading file: arcs_a.dat'
do i__=1,273
read(1,'(A406)')ar__
read(ar__,1)
+ HD(i__),HIP(i__),TYC1(i__),TYC2(i__),TYC3(i__),SpType(i__),
+ RAdeg(i__),DEdeg(i__),Glon(i__),Glat(i__),Hpmag(i__),
+ e_Hpmag(i__),V_I(i__),e_V_I(i__),plx1(i__),e_plx1(i__),
+ dmax(i__),dmin(i__),plx2(i__),e_plx2(i__),dmax2(i__),
+ dmin2(i__),n_dmin2(i__),BTmag(i__),e_BTmag(i__),VTmag(i__),
+ e_VTmag(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__),
+ Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),Kmag(i__),
+ e_Kmag(i__),q2M(i__),Imag(i__),e_Imag(i__),q_Imag(i__),
+ HJD(i__),HRV(i__),e_HRV(i__),TCor(i__),Teff(i__),e_Teff(i__),
+ logg(i__),e_logg(i__),v_M_H_(i__),e__M_H_(i__),Vrot(i__),
+ e_Vrot(i__),d(i__),e_d(i__),E_B_V(i__),X(i__),Y(i__),Z(i__),
+ Uvel(i__),Vvel(i__),Wvel(i__),Com(i__)
if(ar__(1:6) .EQ. '') HD(i__) = iNULL__
if(ar__(9:14) .EQ. '') HIP(i__) = iNULL__
if(ar__(87:92) .EQ. '') Hpmag(i__) = rNULL__
if(ar__(94:99) .EQ. '') e_Hpmag(i__) = rNULL__
if(ar__(101:104) .EQ. '') V_I(i__) = rNULL__
if(ar__(106:109) .EQ. '') e_V_I(i__) = rNULL__
if(ar__(111:117) .EQ. '') plx1(i__) = rNULL__
if(ar__(119:124) .EQ. '') e_plx1(i__) = rNULL__
if(ar__(126:130) .EQ. '') dmax(i__) = iNULL__
if(ar__(132:136) .EQ. '') dmin(i__) = iNULL__
if(ar__(138:144) .EQ. '') plx2(i__) = rNULL__
if(ar__(146:151) .EQ. '') e_plx2(i__) = rNULL__
if(ar__(153:158) .EQ. '') dmax2(i__) = iNULL__
if(ar__(160:165) .EQ. '') dmin2(i__) = iNULL__
if(ar__(193:198) .EQ. '') pmRA(i__) = rNULL__
if(ar__(200:203) .EQ. '') e_pmRA(i__) = rNULL__
if(ar__(205:210) .EQ. '') pmDE(i__) = rNULL__
if(ar__(212:215) .EQ. '') e_pmDE(i__) = rNULL__
if(ar__(260:265) .EQ. '') Imag(i__) = rNULL__
if(ar__(267:270) .EQ. '') e_Imag(i__) = rNULL__
if(ar__(272:274) .EQ. '') q_Imag(i__) = iNULL__
if(ar__(333:337) .EQ. '') Vrot(i__) = rNULL__
if(ar__(339:342) .EQ. '') e_Vrot(i__) = rNULL__
if(ar__(344:348) .EQ. '') d(i__) = iNULL__
if(ar__(350:352) .EQ. '') e_d(i__) = iNULL__
if(ar__(364:368) .EQ. '') X(i__) = iNULL__
if(ar__(370:374) .EQ. '') Y(i__) = iNULL__
if(ar__(376:380) .EQ. '') Z(i__) = iNULL__
if(ar__(382:386) .EQ. '') Uvel(i__) = iNULL__
if(ar__(388:392) .EQ. '') Vvel(i__) = iNULL__
if(ar__(394:398) .EQ. '') Wvel(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ HD(i__),HIP(i__),TYC1(i__),TYC2(i__),TYC3(i__),SpType(i__),
+ RAdeg(i__),DEdeg(i__),Glon(i__),Glat(i__),Hpmag(i__),
+ e_Hpmag(i__),V_I(i__),e_V_I(i__),plx1(i__),e_plx1(i__),
+ dmax(i__),dmin(i__),plx2(i__),e_plx2(i__),dmax2(i__),
+ dmin2(i__),n_dmin2(i__),BTmag(i__),e_BTmag(i__),VTmag(i__),
+ e_VTmag(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__),
+ Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),Kmag(i__),
+ e_Kmag(i__),q2M(i__),Imag(i__),e_Imag(i__),q_Imag(i__),
+ HJD(i__),HRV(i__),e_HRV(i__),TCor(i__),Teff(i__),e_Teff(i__),
+ logg(i__),e_logg(i__),v_M_H_(i__),e__M_H_(i__),Vrot(i__),
+ e_Vrot(i__),d(i__),e_d(i__),E_B_V(i__),X(i__),Y(i__),Z(i__),
+ Uvel(i__),Vvel(i__),Wvel(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'arcs_b.dat' ! Parameters of 70 brighter red clump stars
C Format for file interpretation
2 format(
+ I6,2X,I6,2X,I5,1X,I5,1X,I1,1X,A11,1X,F12.8,1X,F12.8,1X,F8.4,
+ 1X,F8.4,1X,F6.4,1X,F6.4,1X,F4.2,1X,F4.2,1X,F7.2,1X,F6.2,1X,I5,
+ 1X,I5,1X,F7.2,1X,F6.2,1X,I6,1X,I6,A1,F6.3,1X,F5.3,1X,F6.3,1X,
+ F5.3,1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.3,1X,F5.3,1X,F6.3,
+ 1X,F5.3,1X,F6.3,1X,F5.3,1X,A3,1X,F6.3,1X,F4.2,1X,I3,1X,F11.3,
+ 1X,F6.1,1X,F3.1,1X,F4.2,1X,I4,1X,I2,1X,F4.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F5.1,1X,F4.1,1X,I5,1X,I3,1X,F9.7,1X,I5,1X,I5,1X,I5,
+ 1X,I5,1X,I5,1X,I5,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'arcs_b.dat')
write(6,*) '....Loading file: arcs_b.dat'
do i__=1,103
read(1,'(A406)')ar__1
read(ar__1,2)
+ HD_1(i__),HIP_1(i__),TYC1_1(i__),TYC2_1(i__),TYC3_1(i__),
+ SpType_1(i__),RAdeg_1(i__),DEdeg_1(i__),Glon_1(i__),
+ Glat_1(i__),Hpmag_1(i__),e_Hpmag_1(i__),V_I_1(i__),
+ e_V_I_1(i__),plx1_1(i__),e_plx1_1(i__),dmax_1(i__),
+ dmin_1(i__),plx2_1(i__),e_plx2_1(i__),dmax2_1(i__),
+ dmin2_1(i__),n_dmin2_1(i__),BTmag_1(i__),e_BTmag_1(i__),
+ VTmag_1(i__),e_VTmag_1(i__),pmRA_1(i__),e_pmRA_1(i__),
+ pmDE_1(i__),e_pmDE_1(i__),Jmag_1(i__),e_Jmag_1(i__),
+ Hmag_1(i__),e_Hmag_1(i__),Kmag_1(i__),e_Kmag_1(i__),
+ q2M_1(i__),Imag_1(i__),e_Imag_1(i__),q_Imag_1(i__),HJD_1(i__),
+ HRV_1(i__),e_HRV_1(i__),TCor_1(i__),Teff_1(i__),e_Teff_1(i__),
+ logg_1(i__),e_logg_1(i__),v_M_H__1(i__),e__M_H__1(i__),
+ Vrot_1(i__),e_Vrot_1(i__),d_1(i__),e_d_1(i__),E_B_V_1(i__),
+ X_1(i__),Y_1(i__),Z_1(i__),Uvel_1(i__),Vvel_1(i__),
+ Wvel_1(i__),Com_1(i__)
if(ar__1(1:6) .EQ. '') HD_1(i__) = iNULL__
if(ar__1(9:14) .EQ. '') HIP_1(i__) = iNULL__
if(ar__1(87:92) .EQ. '') Hpmag_1(i__) = rNULL__
if(ar__1(94:99) .EQ. '') e_Hpmag_1(i__) = rNULL__
if(ar__1(101:104) .EQ. '') V_I_1(i__) = rNULL__
if(ar__1(106:109) .EQ. '') e_V_I_1(i__) = rNULL__
if(ar__1(111:117) .EQ. '') plx1_1(i__) = rNULL__
if(ar__1(119:124) .EQ. '') e_plx1_1(i__) = rNULL__
if(ar__1(126:130) .EQ. '') dmax_1(i__) = iNULL__
if(ar__1(132:136) .EQ. '') dmin_1(i__) = iNULL__
if(ar__1(138:144) .EQ. '') plx2_1(i__) = rNULL__
if(ar__1(146:151) .EQ. '') e_plx2_1(i__) = rNULL__
if(ar__1(153:158) .EQ. '') dmax2_1(i__) = iNULL__
if(ar__1(160:165) .EQ. '') dmin2_1(i__) = iNULL__
if(ar__1(193:198) .EQ. '') pmRA_1(i__) = rNULL__
if(ar__1(200:203) .EQ. '') e_pmRA_1(i__) = rNULL__
if(ar__1(205:210) .EQ. '') pmDE_1(i__) = rNULL__
if(ar__1(212:215) .EQ. '') e_pmDE_1(i__) = rNULL__
if(ar__1(260:265) .EQ. '') Imag_1(i__) = rNULL__
if(ar__1(267:270) .EQ. '') e_Imag_1(i__) = rNULL__
if(ar__1(272:274) .EQ. '') q_Imag_1(i__) = iNULL__
if(ar__1(333:337) .EQ. '') Vrot_1(i__) = rNULL__
if(ar__1(339:342) .EQ. '') e_Vrot_1(i__) = rNULL__
if(ar__1(344:348) .EQ. '') d_1(i__) = iNULL__
if(ar__1(350:352) .EQ. '') e_d_1(i__) = iNULL__
if(ar__1(364:368) .EQ. '') X_1(i__) = iNULL__
if(ar__1(370:374) .EQ. '') Y_1(i__) = iNULL__
if(ar__1(376:380) .EQ. '') Z_1(i__) = iNULL__
if(ar__1(382:386) .EQ. '') Uvel_1(i__) = iNULL__
if(ar__1(388:392) .EQ. '') Vvel_1(i__) = iNULL__
if(ar__1(394:398) .EQ. '') Wvel_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ HD_1(i__),HIP_1(i__),TYC1_1(i__),TYC2_1(i__),TYC3_1(i__),
+ SpType_1(i__),RAdeg_1(i__),DEdeg_1(i__),Glon_1(i__),
+ Glat_1(i__),Hpmag_1(i__),e_Hpmag_1(i__),V_I_1(i__),
+ e_V_I_1(i__),plx1_1(i__),e_plx1_1(i__),dmax_1(i__),
+ dmin_1(i__),plx2_1(i__),e_plx2_1(i__),dmax2_1(i__),
+ dmin2_1(i__),n_dmin2_1(i__),BTmag_1(i__),e_BTmag_1(i__),
+ VTmag_1(i__),e_VTmag_1(i__),pmRA_1(i__),e_pmRA_1(i__),
+ pmDE_1(i__),e_pmDE_1(i__),Jmag_1(i__),e_Jmag_1(i__),
+ Hmag_1(i__),e_Hmag_1(i__),Kmag_1(i__),e_Kmag_1(i__),
+ q2M_1(i__),Imag_1(i__),e_Imag_1(i__),q_Imag_1(i__),HJD_1(i__),
+ HRV_1(i__),e_HRV_1(i__),TCor_1(i__),Teff_1(i__),e_Teff_1(i__),
+ logg_1(i__),e_logg_1(i__),v_M_H__1(i__),e__M_H__1(i__),
+ Vrot_1(i__),e_Vrot_1(i__),d_1(i__),e_d_1(i__),E_B_V_1(i__),
+ X_1(i__),Y_1(i__),Z_1(i__),Uvel_1(i__),Vvel_1(i__),
+ Wvel_1(i__),Com_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'arcs_c.dat' ! Parameters of 55 test stars
C Format for file interpretation
3 format(
+ I6,2X,I6,2X,I5,1X,I5,1X,I1,1X,A11,1X,F12.8,1X,F12.8,1X,F8.4,
+ 1X,F8.4,1X,F6.4,1X,F6.4,1X,F4.2,1X,F4.2,1X,F7.2,1X,F6.2,1X,I5,
+ 1X,I5,1X,F7.2,1X,F6.2,1X,I6,1X,I6,A1,F6.3,1X,F5.3,1X,F6.3,1X,
+ F5.3,1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.3,1X,F5.3,1X,F6.3,
+ 1X,F5.3,1X,F6.3,1X,F5.3,1X,A3,1X,F6.3,1X,F4.2,1X,I3,1X,F11.3,
+ 1X,F6.1,1X,F3.1,1X,F4.2,1X,I4,1X,I2,1X,F4.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F5.1,1X,F4.1,1X,I5,1X,I3,1X,F9.7,1X,I5,1X,I5,1X,I5,
+ 1X,I5,1X,I5,1X,I5,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'arcs_c.dat')
write(6,*) '....Loading file: arcs_c.dat'
do i__=1,55
read(1,'(A406)')ar__2
read(ar__2,3)
+ HD_2(i__),HIP_2(i__),TYC1_2(i__),TYC2_2(i__),TYC3_2(i__),
+ SpType_2(i__),RAdeg_2(i__),DEdeg_2(i__),Glon_2(i__),
+ Glat_2(i__),Hpmag_2(i__),e_Hpmag_2(i__),V_I_2(i__),
+ e_V_I_2(i__),plx1_2(i__),e_plx1_2(i__),dmax_2(i__),
+ dmin_2(i__),plx2_2(i__),e_plx2_2(i__),dmax2_2(i__),
+ dmin2_2(i__),n_dmin2_2(i__),BTmag_2(i__),e_BTmag_2(i__),
+ VTmag_2(i__),e_VTmag_2(i__),pmRA_2(i__),e_pmRA_2(i__),
+ pmDE_2(i__),e_pmDE_2(i__),Jmag_2(i__),e_Jmag_2(i__),
+ Hmag_2(i__),e_Hmag_2(i__),Kmag_2(i__),e_Kmag_2(i__),
+ q2M_2(i__),Imag_2(i__),e_Imag_2(i__),q_Imag_2(i__),HJD_2(i__),
+ HRV_2(i__),e_HRV_2(i__),TCor_2(i__),Teff_2(i__),e_Teff_2(i__),
+ logg_2(i__),e_logg_2(i__),v_M_H__2(i__),e__M_H__2(i__),
+ Vrot_2(i__),e_Vrot_2(i__),d_2(i__),e_d_2(i__),E_B_V_1(i__),
+ X_2(i__),Y_2(i__),Z_2(i__),Uvel_2(i__),Vvel_2(i__),
+ Wvel_2(i__),Com_2(i__)
if(ar__2(1:6) .EQ. '') HD_2(i__) = iNULL__
if(ar__2(9:14) .EQ. '') HIP_2(i__) = iNULL__
if(ar__2(87:92) .EQ. '') Hpmag_2(i__) = rNULL__
if(ar__2(94:99) .EQ. '') e_Hpmag_2(i__) = rNULL__
if(ar__2(101:104) .EQ. '') V_I_2(i__) = rNULL__
if(ar__2(106:109) .EQ. '') e_V_I_2(i__) = rNULL__
if(ar__2(111:117) .EQ. '') plx1_2(i__) = rNULL__
if(ar__2(119:124) .EQ. '') e_plx1_2(i__) = rNULL__
if(ar__2(126:130) .EQ. '') dmax_2(i__) = iNULL__
if(ar__2(132:136) .EQ. '') dmin_2(i__) = iNULL__
if(ar__2(138:144) .EQ. '') plx2_2(i__) = rNULL__
if(ar__2(146:151) .EQ. '') e_plx2_2(i__) = rNULL__
if(ar__2(153:158) .EQ. '') dmax2_2(i__) = iNULL__
if(ar__2(160:165) .EQ. '') dmin2_2(i__) = iNULL__
if(ar__2(193:198) .EQ. '') pmRA_2(i__) = rNULL__
if(ar__2(200:203) .EQ. '') e_pmRA_2(i__) = rNULL__
if(ar__2(205:210) .EQ. '') pmDE_2(i__) = rNULL__
if(ar__2(212:215) .EQ. '') e_pmDE_2(i__) = rNULL__
if(ar__2(260:265) .EQ. '') Imag_2(i__) = rNULL__
if(ar__2(267:270) .EQ. '') e_Imag_2(i__) = rNULL__
if(ar__2(272:274) .EQ. '') q_Imag_2(i__) = iNULL__
if(ar__2(333:337) .EQ. '') Vrot_2(i__) = rNULL__
if(ar__2(339:342) .EQ. '') e_Vrot_2(i__) = rNULL__
if(ar__2(344:348) .EQ. '') d_2(i__) = iNULL__
if(ar__2(350:352) .EQ. '') e_d_2(i__) = iNULL__
if(ar__2(364:368) .EQ. '') X_2(i__) = iNULL__
if(ar__2(370:374) .EQ. '') Y_2(i__) = iNULL__
if(ar__2(376:380) .EQ. '') Z_2(i__) = iNULL__
if(ar__2(382:386) .EQ. '') Uvel_2(i__) = iNULL__
if(ar__2(388:392) .EQ. '') Vvel_2(i__) = iNULL__
if(ar__2(394:398) .EQ. '') Wvel_2(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ HD_2(i__),HIP_2(i__),TYC1_2(i__),TYC2_2(i__),TYC3_2(i__),
+ SpType_2(i__),RAdeg_2(i__),DEdeg_2(i__),Glon_2(i__),
+ Glat_2(i__),Hpmag_2(i__),e_Hpmag_2(i__),V_I_2(i__),
+ e_V_I_2(i__),plx1_2(i__),e_plx1_2(i__),dmax_2(i__),
+ dmin_2(i__),plx2_2(i__),e_plx2_2(i__),dmax2_2(i__),
+ dmin2_2(i__),n_dmin2_2(i__),BTmag_2(i__),e_BTmag_2(i__),
+ VTmag_2(i__),e_VTmag_2(i__),pmRA_2(i__),e_pmRA_2(i__),
+ pmDE_2(i__),e_pmDE_2(i__),Jmag_2(i__),e_Jmag_2(i__),
+ Hmag_2(i__),e_Hmag_2(i__),Kmag_2(i__),e_Kmag_2(i__),
+ q2M_2(i__),Imag_2(i__),e_Imag_2(i__),q_Imag_2(i__),HJD_2(i__),
+ HRV_2(i__),e_HRV_2(i__),TCor_2(i__),Teff_2(i__),e_Teff_2(i__),
+ logg_2(i__),e_logg_2(i__),v_M_H__2(i__),e__M_H__2(i__),
+ Vrot_2(i__),e_Vrot_2(i__),d_2(i__),e_d_2(i__),E_B_V_1(i__),
+ X_2(i__),Y_2(i__),Z_2(i__),Uvel_2(i__),Vvel_2(i__),
+ Wvel_2(i__),Com_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end