Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/152/213 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-17 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/152/213 Interferometry and spectroscopy of sigma Orionis (Schaefer+, 2016) *================================================================================ *Orbits, distance, and stellar masses of the massive triple star {sigma} Orionis. * Schaefer G.H., Hummel C.A., Gies D.R., Zavala R.T., Monnier J.D., * Walter F.M., Turner N.H., Baron F., ten Brummelaar T., Che X., * Farrington C.D., Kraus S., Sturmann J., Sturmann L. * <Astron. J., 152, 213-213 (2016)> * =2016AJ....152..213S (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Positions of {sigma} Orionis Aa,Ab measured with CHARA Michigan Infrared Combiner (MIRC) integer*4 nr__ parameter (nr__=9) ! Number of records character*99 ar__ ! Full-size record character*11 Date (nr__) ! ("Y:M:D") UT date of the Michigan Infrared Combiner (MIRC) * observation obtained at the CHARA Array real*8 HJD (nr__) ! (d) Median Heliocentric Julian Date for the time * of observation (HJD-2400000) real*4 Sep (nr__) ! (mas) [1/7.4] Separation {rho} real*8 PA (nr__) ! (deg) [74.5/250.4] Position Angle {theta} (measured * east of north) real*4 emaj (nr__) ! (mas) [0.0046/0.32] Semimajor axis of the error * ellipse {sigma}_maj_ (1) real*4 emin (nr__) ! (mas) [0.0019/0.053] Semiminor axis of the error * ellipse {sigma}_min_ (1) real*4 ePA (nr__) ! (deg) [41.1/175.7] Position angle of the major axis * of the error ellipse real*4 fAa (nr__) ! [0.41/0.6] Fractional flux contribution for Aa * component real*4 e_fAa (nr__) ! [0.0024/0.05] Error in fAa real*4 fAb (nr__) ! [0.19/0.39] Fractional flux contribution for Ab * component real*4 e_fAb (nr__) ! [0.002/0.045] Error in fAb real*4 fB (nr__) ! [0.09/0.35] Fractional flux contribution for B * component real*4 e_fB (nr__) ! [0.003/0.07] Error in fB *Note (1): The size of the major and minor axes of the error ellipse for the * binary positions have been scaled by a factor of 2.24 to force * {chi}^2^_{nu}_=1 for an orbit fit to only the MIRC derived positions. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Navy Prototype Optical Interferometer (NPOI) calibrators integer*4 nr__1 parameter (nr__1=33) ! Number of records character*48 ar__1 ! Full-size record integer*4 HD (nr__1) ! [886/222603] HD identification number character*10 SpT (nr__1) ! Spectral type real*4 Vmag (nr__1) ! (mag) [1.7/5.12] V-band magnitude real*4 V_K (nr__1) ! (mag) [-0.94/1.31] The V-K color index real*4 E_B_V (nr__1) ! (mag) [0/0.17] The (B-V) color excess real*4 Diam (nr__1) ! (mas) [0.32/1.01] Angular diameter {theta}_V-K_ * corrected for extinction real*4 V2min (nr__1) ! [0.45/1] Minimum estimated calibrator visibility * V^2^_min_ based on {theta}_V-K_ integer*4 Nn (nr__1) ! [1/16] Number of nights c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Navy Prototype Optical Interferometer (NPOI) positions for {sigma} Orionis A,B integer*4 nr__2 parameter (nr__2=25) ! Number of records character*67 ar__2 ! Full-size record character*11 Date_1 (nr__2) ! ("Y:M:D") UT date of the observation real*8 Year (nr__2) ! (yr) Julian year at 7UT on the date of observation integer*4 Nv (nr__2) ! [30/1479] Number of measured visibilities real*4 Sep_1 (nr__2) ! (mas) [254.39/260.14] Derived separation {rho} (1) real*4 PA_1 (nr__2) ! (deg) [84.17/109.8] Position Angle {theta} real*4 emaj_1 (nr__2) ! (mas) [0.59/7.59] Semimajor axis of the uncertainty * ellipse {sigma}_maj_ real*4 emin_1 (nr__2) ! (mas) [0.4/2.25] Semiminor axis of the uncertainty * ellipse {sigma}_min_ real*4 ePA_1 (nr__2) ! (deg) [7.4/175.3] Position angle of the major axis * of the uncertainty ellipse real*4 O_Cr (nr__2) ! (mas) [-1.7/3.9] Deviation of the fitted relative binary * position from the model values (O-C)_{rho}_ real*4 O_Ct (nr__2) ! (deg) [-1.6/1] Deviation of the fitted relative binary * position from the model values (O-C)_{theta}_ *Note (1): Relative to the center of mass of the close pair. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Navy Prototype Optical Interferometer (NPOI) positions for {sigma} Orionis Aa,Ab integer*4 nr__3 parameter (nr__3=17) ! Number of records character*65 ar__3 ! Full-size record character*11 Date_2 (nr__3) ! ("Y:M:D") UT date of the observation real*8 HJD_1 (nr__3) ! (yr) Heliocentric Julian Date (HJD-2400000) (1) integer*4 Nv_1 (nr__3) ! [75/1479] Number of measured visibilities real*4 Sep_2 (nr__3) ! (mas) [0.49/7.4] Derived separation {rho} real*4 PA_2 (nr__3) ! (deg) [123.7/326.3] Position Angle {theta} real*4 emaj_2 (nr__3) ! (mas) [0.08/0.9] Semimajor axis of the uncertainty * ellipse {sigma}_maj_ real*4 emin_2 (nr__3) ! (mas) [0.06/0.26] Semiminor axis of the uncertainty * ellipse {sigma}_min_ real*4 ePA_2 (nr__3) ! (deg) [0.1/179.4] Position angle of the major axis * of the uncertainty ellipse real*4 O_Cr_1 (nr__3) ! (mas) [-0.6/1] Deviation of the fitted relative binary * position from the model values (O-C)_{rho}_ real*4 O_Ct_1 (nr__3) ! (deg) [-34/9.8] Deviation of the fitted relative binary * position from the model values (O-C)_{theta}_ *Note (1): Computed at local midnight (7UT) on the date of observation. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Cerro Tololo Inter-American Observatory (CTIO) radial velocities measured for {sigma} Orionis Aa,Ab integer*4 nr__4 parameter (nr__4=61) ! Number of records character*60 ar__4 ! Full-size record character*11 Date_3 (nr__4) ! ("Y:M:D") UT date real*8 HJD_2 (nr__4) ! (d) Heliocentric Julian Date (HJD-2400000) real*4 Phase (nr__4) ! [0/1] Orbital phase real*4 RVAa (nr__4) ! (km/s) [-87.9/59.3] Radial velocity for Aa component real*4 e_RVAa (nr__4) ! (km/s) [2.6/50.6] Error in RVAa real*4 RVAb (nr__4) ! (km/s) [13.3/204] Radial velocity for Ab component real*4 e_RVAb (nr__4) ! (km/s) [1.8/17] Error in RVAb real*4 O_CAa (nr__4) ! (km/s) [-6/24.1] The radial velocity residual for * Aa component real*4 O_CAb (nr__4) ! (km/s) [-6.7/8.4] The radial velocity residual for * Ab component C============================================================================= C Loading file 'table3.dat' ! Positions of {sigma} Orionis Aa,Ab measured * with CHARA Michigan Infrared Combiner (MIRC) C Format for file interpretation 1 format( + A11,1X,F9.3,1X,F6.4,1X,F7.3,1X,F6.4,1X,F6.4,1X,F6.2,1X,F6.4, + 1X,F6.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F6.4) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,9 read(1,'(A99)')ar__ read(ar__,1) + Date(i__),HJD(i__),Sep(i__),PA(i__),emaj(i__),emin(i__), + ePA(i__),fAa(i__),e_fAa(i__),fAb(i__),e_fAb(i__),fB(i__), + e_fB(i__) c ..............Just test output........... write(6,1) + Date(i__),HJD(i__),Sep(i__),PA(i__),emaj(i__),emin(i__), + ePA(i__),fAa(i__),e_fAa(i__),fAb(i__),e_fAb(i__),fB(i__), + e_fB(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! Navy Prototype Optical Interferometer (NPOI) * calibrators C Format for file interpretation 2 format(3X,I6,1X,A10,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,I2) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,33 read(1,'(A48)')ar__1 read(ar__1,2) + HD(i__),SpT(i__),Vmag(i__),V_K(i__),E_B_V(i__),Diam(i__), + V2min(i__),Nn(i__) c ..............Just test output........... write(6,2) + HD(i__),SpT(i__),Vmag(i__),V_K(i__),E_B_V(i__),Diam(i__), + V2min(i__),Nn(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table6.dat' ! Navy Prototype Optical Interferometer (NPOI) * positions for {sigma} Orionis A,B C Format for file interpretation 3 format( + A11,1X,F9.4,1X,I4,1X,F6.2,1X,F6.2,1X,F4.2,1X,F4.2,1X,F5.1,1X, + F5.2,1X,F4.1) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,25 read(1,'(A67)')ar__2 read(ar__2,3) + Date_1(i__),Year(i__),Nv(i__),Sep_1(i__),PA_1(i__), + emaj_1(i__),emin_1(i__),ePA_1(i__),O_Cr(i__),O_Ct(i__) c ..............Just test output........... write(6,3) + Date_1(i__),Year(i__),Nv(i__),Sep_1(i__),PA_1(i__), + emaj_1(i__),emin_1(i__),ePA_1(i__),O_Cr(i__),O_Ct(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! Navy Prototype Optical Interferometer (NPOI) * positions for {sigma} Orionis Aa,Ab C Format for file interpretation 4 format( + A11,1X,F8.2,1X,I4,1X,F4.2,1X,F6.2,1X,F4.2,1X,F4.2,1X,F5.1,1X, + F5.2,1X,F5.1) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,17 read(1,'(A65)')ar__3 read(ar__3,4) + Date_2(i__),HJD_1(i__),Nv_1(i__),Sep_2(i__),PA_2(i__), + emaj_2(i__),emin_2(i__),ePA_2(i__),O_Cr_1(i__),O_Ct_1(i__) c ..............Just test output........... write(6,4) + Date_2(i__),HJD_1(i__),Nv_1(i__),Sep_2(i__),PA_2(i__), + emaj_2(i__),emin_2(i__),ePA_2(i__),O_Cr_1(i__),O_Ct_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! Cerro Tololo Inter-American Observatory (CTIO) * radial velocities measured for {sigma} Orionis * Aa,Ab C Format for file interpretation 5 format( + A11,1X,F9.3,1X,F5.3,1X,F6.2,1X,F4.1,1X,F5.1,1X,F4.1,1X,F4.1, + 1X,F4.1) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,61 read(1,'(A60)')ar__4 read(ar__4,5) + Date_3(i__),HJD_2(i__),Phase(i__),RVAa(i__),e_RVAa(i__), + RVAb(i__),e_RVAb(i__),O_CAa(i__),O_CAb(i__) c ..............Just test output........... write(6,5) + Date_3(i__),HJD_2(i__),Phase(i__),RVAa(i__),e_RVAa(i__), + RVAb(i__),e_RVAb(i__),O_CAa(i__),O_CAb(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end