Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/162/184 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-Aug-10
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/AJ/162/184 APOGEE: Properties & radial velocities of binaries (Kounkel+, 2021)
*================================================================================
*Double-lined spectroscopic binaries in the APOGEE DR16 and DR17 data.
* Kounkel M., Covey K.R., Stassun K.G., Price-Whelan A.M., Holtzman J.,
* Chojnowski D., Longa-Pena P., Roman-Zuniga C.G., Hernandez J., Serna J.,
* Badenes C., De Lee N., Majewski S., Stringfellow G.S., Kratter K.M.,
* Moe M., Frinchaboy P.M., Beaton R.L., Fernandez-Trincado J.G.,
* Mahadevan S., Minniti D., Beers T.C., Schneider D.P., Barba R.,
* Brownstein J.R., Garcia-Hernandez D.A., Pan K., Bizyaev D.
* <Astron. J., 162, 184 (2021)>
* =2021AJ....162..184K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Properties of the identified SB2s and higher order
multiples
integer*4 nr__
parameter (nr__=8105) ! Number of records
character*582 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*18 ID (nr__) ! APOGEE identifier
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
integer*4 SBn (nr__) ! [2/4] Number of deconvolved components for a
* source (1)
integer*4 Nepoch (nr__) ! [1/47] Number of epochs available for the
* source
real*8 qW (nr__) ! [-373.1/79.2]? Mass ratio from a Wilson plot
real*4 e_qW (nr__) ! [0/65947.1]? qW uncertainty
real*8 RVelW (nr__) ! (km/s) [-3748/291]? Barycenter radial velocity from a
* Wilson plot
real*4 e_RVelW (nr__) ! (km/s) [-58300/2200]? RVelW uncertainty
real*8 Per (nr__) ! (d) [0.8/348]? Period (2)
real*4 e_Per (nr__) ! (d) [3.1e-7/52]? Per uncertainty
real*8 T0 (nr__) ! (d) [2.4e+6/2.5e+6]? Time of periastron passage (2)
real*4 e_T0 (nr__) ! (d) [8.8e-5/52]? T0 uncertainty
real*4 e (nr__) ! [1.07e-5/0.92]? Eccentricity (2)
real*4 e_e (nr__) ! [9.5e-6/0.25]? e uncertainty
real*8 omega (nr__) ! (deg) [-46.5/608]? Longitude of periastron (2)
real*8 e_omega (nr__) ! (deg) [0.01/271]? omega uncertainty
real*8 RVel (nr__) ! (km/s) [-117.2/243]? Barycenter radial velocity (2)
real*4 e_RVel (nr__) ! (km/s) [3.6e-4/6.4]? RVel uncertainty
real*8 K1 (nr__) ! (km/s) [7.1/173.1]? Semiamplitude of the velocity for
* component 1 (2)
real*4 e_K1 (nr__) ! (km/s) [6.7e-4/71]? K1 uncertainty
real*8 K2 (nr__) ! (km/s) [10.8/169]? Semiamplitude of the velocity for
* component 2 (2)
real*4 e_K2 (nr__) ! (km/s) [7.8e-4/40.1]? K2 uncertainty
real*4 M1sin3i (nr__) ! (Msun) [2.7e-4/4.8]? Inclination-dependent mass of
* component 1 (2)
real*4 e_M1sin3i (nr__) ! (Msun) [6.3e-6/5.3]? M1sin3i uncertainty
real*4 M2sin3i (nr__) ! (Msun) [1.5e-4/6.4]? Inclination-dependent mass of
* component 2 (2)
real*4 e_M2sin3i (nr__) ! (Msun) [5.4e-6/5.5]? M2sin3i uncertainty
real*4 asini (nr__) ! (km) [220359/1.9e+8]? Inclination-dependent
* semi-major axis (2)
real*4 e_asini (nr__) ! (km) [46.6/4.6e+7]? asini uncertainty
real*8 MaxdRV (nr__) ! (km/s) [18.3/407.6] Maximum observed separation in RV
* between the primary and the secondary in any
* epoch
real*4 AmpRV1 (nr__) ! (km/s) [0/354.4]? Observed amplitude of variation in
* radial velocity of component 1 across all of
* the available data
real*4 AmpRV2 (nr__) ! (km/s) [0/387.5]? Observed amplitude of variation in
* radial velocity of component 2 across all of
* the available data
real*8 MaxT (nr__) ! (d) [0/2979]? Maximum temporal baseline
character*5 fLOS (nr__) ! Line-of-sight coincidence flag (3)
character*3 fVar (nr__) ! Flag for a variable light curve in TESS data
* (4)
real*8 PerTESS (nr__) ! (d) [0.22/172.1]? Period derived from TESS light
* curves for periodic variables
*Note (1): Number of components for a source as follows:
* 2 = 7273 occurrences
* 3 = 813 occurrences
* 4 = 19 occurrences
*Note (2): Derived from the orbital fitting.
*Note (3): Flags as follows:
* true = Likely line-of-sight coincidence, defined as having the total radial
* velocity change for each component of less than 5 km/s; 286 occurrences
* false = Unlikely to be a line-of-sight coincidence; 7819 occurrences
*Note (4): Flags as follows:
* f = in TESS footprint (3981 occurrences)
* v = variable (1135 occurrences)
* d = detached eclipsing (369 occurrences)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Vetted radial velocities of the individual
components
integer*4 nr__1
parameter (nr__1=32642) ! Number of records
character*229 ar__1 ! Full-size record
character*18 ID_1 (nr__1) ! APOGEE identifier
real*4 HJD (nr__1) ! (d) [55804/58933] Heliocentric Julian Date;
* HJD-2400000
real*4 RVel1 (nr__1) ! (km/s) [-317.9/327.4]? Radial velocity of component 1
real*4 e_RVel1 (nr__1) ! (km/s) [0/170]? RV1 uncertainty
real*4 RVel2 (nr__1) ! (km/s) [-302.4/329.7]? Radial velocity of component 2
real*4 e_RVel2 (nr__1) ! (km/s) [0/180]? RV2 uncertainty
real*4 RVel3 (nr__1) ! (km/s) [-3.3e+7/6.8e+9]? Radial velocity of component 3
real*4 e_RVel3 (nr__1) ! (km/s) [0/2.5e+10]? RV3 uncertainty
real*4 RVel4 (nr__1) ! (km/s) [-98/157.6]? Radial velocity of component 4
real*4 e_RVel4 (nr__1) ! (km/s) [0/20]? RV4 uncertainty
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'dr17samp.dat' ! Sample of systems from Table A1 (Table added by CDS)
integer*4 nr__2
parameter (nr__2=7204) ! Number of records
character*49 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*18 ID_2 (nr__2) ! APOGEE identifier
real*8 RAdeg_1 (nr__2) ! (deg) Right ascension (J2000)
real*8 DEdeg_1 (nr__2) ! (deg) Declination (J2000)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Parameters of the cross-correlation function (CCF)
components extracted by the pipeline in DR17 data
integer*4 nr__3
parameter (nr__3=31313) ! Number of records
character*284 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*18 ID_3 (nr__3) ! APOGEE identifier
real*8 HJD_1 (nr__3) ! (d) [55804/59236] Heliocentric Julian Date;
* HJD-2400000
real*8 RAdeg_2 (nr__3) ! (deg) Right ascension (J2000)
real*8 DEdeg_2 (nr__3) ! (deg) Declination (J2000)
real*4 Amp1 (nr__3) ! [0/0.65]? Amplitude of component 1
real*8 e_Amp1 (nr__3) ! [2e-3/4.7e+5]? Amp1 uncertainty
real*4 Amp2 (nr__3) ! [0/3.5]? Amplitude of component 2
real*4 e_Amp2 (nr__3) ! [1.8e-3/2.8e+8]? Amp2 uncertainty
real*4 Amp3 (nr__3) ! [0/1.03]? Amplitude of component 3
real*4 e_Amp3 (nr__3) ! [6e-3/2.2e+11]? Amp3 uncertainty
real*4 Amp4 (nr__3) ! [0/0.3]? Amplitude of component 4
real*4 e_Amp4 (nr__3) ! [0.016/6.6]? Amp4 uncertainty
real*8 RVel1_1 (nr__3) ! (km/s) [-318/328.7]? Radial velocity of component 1
real*8 e_RVel1_1 (nr__3) ! (km/s) [0.6/1e+6]? RVel1 uncertainty
real*8 RVel2_1 (nr__3) ! (km/s) [-298/329]? Radial velocity of component 2
real*4 e_RVel2_1 (nr__3) ! (km/s) [0.6/7.7e+9]? RVel2 uncertainty
real*8 RVel3_1 (nr__3) ! (km/s) [-179.1/321]? Radial velocity of component 3
real*4 e_RVel3_1 (nr__3) ! (km/s) [0.9/1.7e+10]? RVel3 uncertainty
real*4 e_RVel4_1 (nr__3) ! (km/s) [1.5/100]? RVel4 uncertainty
real*8 RVel4_1 (nr__3) ! (km/s) [-144/200.8]? Radial velocity of component 4
real*8 FWHM1 (nr__3) ! [0/1861.1]? Full width at half maximum of
* component 1
real*8 e_FWHM1 (nr__3) ! [1.4/8.7e+4]? FWHM1 uncertainty
real*8 FWHM2 (nr__3) ! [0/1310.3]? Full width at half maximum of
* component 2
real*4 e_FWHM2 (nr__3) ! [1e-3/1.5e+11]? FWHM2 uncertainty
real*4 FWHM3 (nr__3) ! [0/648]? Full width at half maximum of comp
*onent 3
real*4 e_FWHM3 (nr__3) ! [2.3/6.8e+10]? FWHM3 uncertainty
real*4 FWHM4 (nr__3) ! [0.11/55]? Full width at half maximum of comp
*onent 4
real*4 e_FWHM4 (nr__3) ! [3.5/130]? FWHM4 uncertainty
integer*4 Flag1 (nr__3) ! [0/4]? Classification flag for component 1 (1)
integer*4 Flag2 (nr__3) ! [0/4]? Classification flag for component 2 (1)
integer*4 Flag3 (nr__3) ! [0/4]? Classification flag for component 3 (1)
integer*4 Flag4 (nr__3) ! [0/2]? Classification flag for component 4 (1)
integer*4 N (nr__3) ! [0/4] Number of components (2)
character*8 DR16comp (nr__3) ! Comparison to DR16 (3)
*Note (1): Classification flags as follows:
*
* 0 = No component.
* Comp 1: 308 occurrences;
* Comp 2: 7979 occurrences;
* Comp 3: 30429 occurrences;
* Comp 4: 31286 occurrences
*
* 1 = Failed FWHM/amplitude test; could be noise.
* Comp 1: 901 occurrences;
* Comp 2: 3659 occurrences;
* Comp 3: 701 occurrences;
* Comp 4: 25 occurrences
*
* 2 = Failed the symmetry test; could be falsely multiply deconvolved.
* Comp 1: 0 occurrence;
* Comp 2: 3659 occurrences;
* Comp 3: 61 occurrences;
* Comp 4: 2 occurrences
*
* 3 = Blended asymmetric peaks; generally a robust detection of a companion,
* but may be a signature of spots in strongly magnetic stars.
* Comp 1: 0 occurrence;
* Comp 2: 2449 occurrences;
* Comp 3: 36 occurrences;
* Comp 4: 0 occurrence
*
* 4 = The primary peaks in all stars, and the secondary peaks of bona fide
* SB2s that are not blended with the primary.
* Comp 1: 30104 occurrences
* Comp 2: 9842 occurrences
* Comp 3: 86 occurrences
* Comp 4: 0 occurrence
*Note (2): Number of components as follows:
* 0 = 302 occurrences
* 1 = 7668 occurrences
* 2 = 22446 occurrences
* 3 = 870 occurrences
* 4 = 27 occurrences
*Note (3): Comparison as follows:
* New obs = Epoch not a part of Table 3 (10887 occurrences)
* New comp = Epoch in Table 3 is unresolved with a single component;
* component with Flag>=3 is detected in DR17 (242 occurrences)
C=============================================================================
C Loading file 'table1.dat' ! Properties of the identified SB2s and higher order
* multiples
C Format for file interpretation
1 format(
+ A18,1X,F15.11,1X,F17.13,1X,I1,1X,I2,1X,F24.19,1X,E29.23,1X,
+ F24.18,1X,E23.19,1X,F11.7,1X,E13.10,1X,F9.1,1X,E12.10,1X,
+ E10.8,1X,E12.10,1X,F10.6,1X,F13.9,1X,F11.6,1X,E12.10,1X,F10.6,
+ 1X,E12.9,1X,F10.6,1X,E11.9,1X,E21.19,1X,E21.19,1X,E21.19,1X,
+ E21.19,1X,E20.17,1X,E20.17,1X,F19.15,1X,E22.19,1X,E21.19,1X,
+ F23.18,1X,A5,1X,A3,1X,F21.17)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,8105
read(1,'(A582)')ar__
read(ar__,1)
+ ID(i__),RAdeg(i__),DEdeg(i__),SBn(i__),Nepoch(i__),qW(i__),
+ e_qW(i__),RVelW(i__),e_RVelW(i__),Per(i__),e_Per(i__),T0(i__),
+ e_T0(i__),e(i__),e_e(i__),omega(i__),e_omega(i__),RVel(i__),
+ e_RVel(i__),K1(i__),e_K1(i__),K2(i__),e_K2(i__),M1sin3i(i__),
+ e_M1sin3i(i__),M2sin3i(i__),e_M2sin3i(i__),asini(i__),
+ e_asini(i__),MaxdRV(i__),AmpRV1(i__),AmpRV2(i__),MaxT(i__),
+ fLOS(i__),fVar(i__),PerTESS(i__)
if(ar__(59:82) .EQ. '') qW(i__) = rNULL__
if(ar__(84:112) .EQ. '') e_qW(i__) = rNULL__
if(ar__(114:137) .EQ. '') RVelW(i__) = rNULL__
if(ar__(139:161) .EQ. '') e_RVelW(i__) = rNULL__
if(ar__(163:173) .EQ. '') Per(i__) = rNULL__
if(ar__(175:187) .EQ. '') e_Per(i__) = rNULL__
if(ar__(189:197) .EQ. '') T0(i__) = rNULL__
if(ar__(199:210) .EQ. '') e_T0(i__) = rNULL__
if(ar__(212:221) .EQ. '') e(i__) = rNULL__
if(ar__(223:234) .EQ. '') e_e(i__) = rNULL__
if(ar__(236:245) .EQ. '') omega(i__) = rNULL__
if(ar__(247:259) .EQ. '') e_omega(i__) = rNULL__
if(ar__(261:271) .EQ. '') RVel(i__) = rNULL__
if(ar__(273:284) .EQ. '') e_RVel(i__) = rNULL__
if(ar__(286:295) .EQ. '') K1(i__) = rNULL__
if(ar__(297:308) .EQ. '') e_K1(i__) = rNULL__
if(ar__(310:319) .EQ. '') K2(i__) = rNULL__
if(ar__(321:331) .EQ. '') e_K2(i__) = rNULL__
if(ar__(333:353) .EQ. '') M1sin3i(i__) = rNULL__
if(ar__(355:375) .EQ. '') e_M1sin3i(i__) = rNULL__
if(ar__(377:397) .EQ. '') M2sin3i(i__) = rNULL__
if(ar__(399:419) .EQ. '') e_M2sin3i(i__) = rNULL__
if(ar__(421:440) .EQ. '') asini(i__) = rNULL__
if(ar__(442:461) .EQ. '') e_asini(i__) = rNULL__
if(ar__(483:504) .EQ. '') AmpRV1(i__) = rNULL__
if(ar__(506:526) .EQ. '') AmpRV2(i__) = rNULL__
if(ar__(528:550) .EQ. '') MaxT(i__) = rNULL__
if(ar__(562:582) .EQ. '') PerTESS(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),RAdeg(i__),DEdeg(i__),SBn(i__),Nepoch(i__),qW(i__),
+ e_qW(i__),RVelW(i__),e_RVelW(i__),Per(i__),e_Per(i__),T0(i__),
+ e_T0(i__),e(i__),e_e(i__),omega(i__),e_omega(i__),RVel(i__),
+ e_RVel(i__),K1(i__),e_K1(i__),K2(i__),e_K2(i__),M1sin3i(i__),
+ e_M1sin3i(i__),M2sin3i(i__),e_M2sin3i(i__),asini(i__),
+ e_asini(i__),MaxdRV(i__),AmpRV1(i__),AmpRV2(i__),MaxT(i__),
+ fLOS(i__),fVar(i__),PerTESS(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Vetted radial velocities of the individual
* components
C Format for file interpretation
2 format(
+ A18,1X,E22.17,1X,E23.17,1X,E22.17,1X,E23.17,1X,E22.17,1X,
+ E23.17,1X,E22.17,1X,E23.17,1X,E22.17)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,32642
read(1,'(A229)')ar__1
read(ar__1,2)
+ ID_1(i__),HJD(i__),RVel1(i__),e_RVel1(i__),RVel2(i__),
+ e_RVel2(i__),RVel3(i__),e_RVel3(i__),RVel4(i__),e_RVel4(i__)
if(ar__1(43:65) .EQ. '') RVel1(i__) = rNULL__
if(ar__1(67:88) .EQ. '') e_RVel1(i__) = rNULL__
if(ar__1(90:112) .EQ. '') RVel2(i__) = rNULL__
if(ar__1(114:135) .EQ. '') e_RVel2(i__) = rNULL__
if(ar__1(137:159) .EQ. '') RVel3(i__) = rNULL__
if(ar__1(161:182) .EQ. '') e_RVel3(i__) = rNULL__
if(ar__1(184:206) .EQ. '') RVel4(i__) = rNULL__
if(ar__1(208:229) .EQ. '') e_RVel4(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1(i__),HJD(i__),RVel1(i__),e_RVel1(i__),RVel2(i__),
+ e_RVel2(i__),RVel3(i__),e_RVel3(i__),RVel4(i__),e_RVel4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'dr17samp.dat' ! Sample of systems from Table A1 (Table added by CDS)
C Format for file interpretation
3 format(A18,1X,F14.10,1X,F15.11)
C Effective file loading
open(unit=1,status='old',file=
+'dr17samp.dat')
write(6,*) '....Loading file: dr17samp.dat'
do i__=1,7204
read(1,'(A49)')ar__2
read(ar__2,3)ID_2(i__),RAdeg_1(i__),DEdeg_1(i__)
c ..............Just test output...........
write(6,3)ID_2(i__),RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea1.dat' ! Parameters of the cross-correlation function (CCF)
* components extracted by the pipeline in DR17 data
C Format for file interpretation
4 format(
+ A18,1X,F18.12,1X,F14.10,1X,F15.11,1X,F6.4,1X,F11.4,1X,F6.4,1X,
+ E9.4,1X,F6.4,1X,E8.4,1X,F5.3,1X,F5.3,1X,F7.2,1X,F10.2,1X,F7.2,
+ 1X,E9.2,1X,F7.2,1X,E9.2,1X,F5.1,1X,F7.2,1X,F7.2,1X,F7.1,1X,
+ F8.3,1X,E9.3,1X,F6.2,1X,E9.2,1X,F5.2,1X,F5.1,1X,I1,1X,I1,1X,
+ I1,1X,I1,1X,I1,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,31313
read(1,'(A284)')ar__3
read(ar__3,4)
+ ID_3(i__),HJD_1(i__),RAdeg_2(i__),DEdeg_2(i__),Amp1(i__),
+ e_Amp1(i__),Amp2(i__),e_Amp2(i__),Amp3(i__),e_Amp3(i__),
+ Amp4(i__),e_Amp4(i__),RVel1_1(i__),e_RVel1_1(i__),
+ RVel2_1(i__),e_RVel2_1(i__),RVel3_1(i__),e_RVel3_1(i__),
+ e_RVel4_1(i__),RVel4_1(i__),FWHM1(i__),e_FWHM1(i__),
+ FWHM2(i__),e_FWHM2(i__),FWHM3(i__),e_FWHM3(i__),FWHM4(i__),
+ e_FWHM4(i__),Flag1(i__),Flag2(i__),Flag3(i__),Flag4(i__),
+ N(i__),DR16comp(i__)
if(ar__3(70:75) .EQ. '') Amp1(i__) = rNULL__
if(ar__3(77:87) .EQ. '') e_Amp1(i__) = rNULL__
if(ar__3(89:94) .EQ. '') Amp2(i__) = rNULL__
if(ar__3(96:104) .EQ. '') e_Amp2(i__) = rNULL__
if(ar__3(106:111) .EQ. '') Amp3(i__) = rNULL__
if(ar__3(113:120) .EQ. '') e_Amp3(i__) = rNULL__
if(ar__3(122:126) .EQ. '') Amp4(i__) = rNULL__
if(ar__3(128:132) .EQ. '') e_Amp4(i__) = rNULL__
if(ar__3(134:140) .EQ. '') RVel1_1(i__) = rNULL__
if(ar__3(142:151) .EQ. '') e_RVel1_1(i__) = rNULL__
if(ar__3(153:159) .EQ. '') RVel2_1(i__) = rNULL__
if(ar__3(161:169) .EQ. '') e_RVel2_1(i__) = rNULL__
if(ar__3(171:177) .EQ. '') RVel3_1(i__) = rNULL__
if(ar__3(179:187) .EQ. '') e_RVel3_1(i__) = rNULL__
if(ar__3(189:193) .EQ. '') e_RVel4_1(i__) = rNULL__
if(ar__3(195:201) .EQ. '') RVel4_1(i__) = rNULL__
if(ar__3(203:209) .EQ. '') FWHM1(i__) = rNULL__
if(ar__3(211:217) .EQ. '') e_FWHM1(i__) = rNULL__
if(ar__3(219:226) .EQ. '') FWHM2(i__) = rNULL__
if(ar__3(228:236) .EQ. '') e_FWHM2(i__) = rNULL__
if(ar__3(238:243) .EQ. '') FWHM3(i__) = rNULL__
if(ar__3(245:253) .EQ. '') e_FWHM3(i__) = rNULL__
if(ar__3(255:259) .EQ. '') FWHM4(i__) = rNULL__
if(ar__3(261:265) .EQ. '') e_FWHM4(i__) = rNULL__
if(ar__3(267:267) .EQ. '') Flag1(i__) = iNULL__
if(ar__3(269:269) .EQ. '') Flag2(i__) = iNULL__
if(ar__3(271:271) .EQ. '') Flag3(i__) = iNULL__
if(ar__3(273:273) .EQ. '') Flag4(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ ID_3(i__),HJD_1(i__),RAdeg_2(i__),DEdeg_2(i__),Amp1(i__),
+ e_Amp1(i__),Amp2(i__),e_Amp2(i__),Amp3(i__),e_Amp3(i__),
+ Amp4(i__),e_Amp4(i__),RVel1_1(i__),e_RVel1_1(i__),
+ RVel2_1(i__),e_RVel2_1(i__),RVel3_1(i__),e_RVel3_1(i__),
+ e_RVel4_1(i__),RVel4_1(i__),FWHM1(i__),e_FWHM1(i__),
+ FWHM2(i__),e_FWHM2(i__),FWHM3(i__),e_FWHM3(i__),FWHM4(i__),
+ e_FWHM4(i__),Flag1(i__),Flag2(i__),Flag3(i__),Flag4(i__),
+ N(i__),DR16comp(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end