Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/835/152 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-16
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/ApJ/835/152 HST FGS-1r parallaxes for 8 metal-poor stars (Chaboyer+, 2017)
*================================================================================
*Testing metal-poor stellar models and isochrones with HST parallaxes of
*metal-poor stars.
* Chaboyer B., McArthur B.E., O'Malley E., Benedict G.F., Feiden G.A.,
* Harrison T.E., McWilliam A., Nelan E.P., Patterson R.J., Sarajedini A.
* <Astrophys. J., 835, 152-152 (2017)>
* =2017ApJ...835..152C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Metal-poor target star parallaxes
integer*4 nr__
parameter (nr__=8) ! Number of records
character*43 ar__ ! Full-size record
integer*4 HIP (nr__) ! [46120/108200] HIP number
real*4 Plx (nr__) ! (mas) [3.4/15.1] HST parallax
real*4 e_Plx (nr__) ! (mas) [0.08/0.2] Plx uncertainty
real*4 PlxH07 (nr__) ! (mas) [9.5/16] HIP parallax (Cat. I/311)
real*4 e_PlxH07 (nr__) ! (mas) [0.9/3] PlxH07 uncertainty
real*4 PlxT17 (nr__) ! (mas) [2.4/15]? TGAS parallax (Cat. I/337/tgas)
real*4 e_PlxT17 (nr__) ! (mas) [0.2/0.9]? PlxT17 uncertainty
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Metal-poor target star proper motions
integer*4 nr__1
parameter (nr__1=8) ! Number of records
character*84 ar__1 ! Full-size record
integer*4 HIP_1 (nr__1) ! [46120/108200] HIP number
real*8 pmRA (nr__1) ! (mas/yr) [-567.1/760.7] {mu}_R.A._ from HST
real*4 e_pmRA (nr__1) ! (mas/yr) [0.07/0.2] pmRA uncertainty
real*8 pmDE (nr__1) ! (mas/yr) [-625.6/1236.4] {mu}_Decl._ from HST
real*4 e_pmDE (nr__1) ! (mas/yr) [0.08/0.2] pmDE uncertainty
real*8 pmRAH07 (nr__1) ! (mas/yr) [-568/763] HIP proper motion along RA
real*4 e_pmRAH07 (nr__1) ! (mas/yr) [1/4] pmRAH07 uncertainty
real*8 pmDEH07 (nr__1) ! (mas/yr) [-627/1237] HIP proper motion along DEC
real*4 e_pmDEH07 (nr__1) ! (mas/yr) [1/3] pmDEH07 uncertainty
real*4 pmRAPP (nr__1) ! (mas/yr) [-382/762] PPMXL {mu}_R.A._ (Cat. I/317)
real*4 e_pmRAPP (nr__1) ! (mas/yr) [1/9] pmRAPP uncertainty
real*4 pmDEPP (nr__1) ! (mas/yr) [-630/1238] PPMXL {mu}_Decl._
real*4 e_pmDEPP (nr__1) ! (mas/yr) [1/7] pmDEPP uncertainty
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Summary of parallax star properties
integer*4 nr__2
parameter (nr__2=8) ! Number of records
character*55 ar__2 ! Full-size record
integer*4 HIP_2 (nr__2) ! [46120/108200] HIP number
real*4 Vmag (nr__2) ! (mag) [10/11.6] V-band magnitude (Vega system)
real*4 V_I (nr__2) ! (mag) [0.7/1] (V-I) color index
real*4 E_B_V (nr__2) ! (mag) [0/0.2] Color excess (1)
real*4 v_V_I_0 (nr__2) ! [0.6/1] Dereddened (V-I) color index
real*4 VMag_1 (nr__2) ! (mag) [3.9/6.7] Absolute V-band magnitude (M_V_)
real*4 e_VMag (nr__2) ! (mag) [0.01/0.1] VMag uncertainty
real*4 F606W (nr__2) ! (mag) [3.7/6.4] F606W absolute magnitude (M_F606W_) (2)
real*4 v_Fe_H_ (nr__2) ! ([-]) [-2.7/-1.4] [Fe/H] abundance (1)
real*4 v_a_Fe_ (nr__2) ! ([-]) [0.01/0.7] [{alpha}/Fe] abundance (1)(3)
*Note (1): High-resolution spectroscopic abundances and reddening estimates are
* from O'Malley+ (2017ApJ...838...90O). The reddening estimates are
* based upon the strength of the interstellar NaI line in the
* high-resolution spectra of the target stars. See section 3.
*Note (2): Uncertainty is the same as the VMag values.
*Note (3): Based upon the abundance of calcium, an {alpha}-capture element.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! Astrometric data for the 8 targets and
reference stars
integer*4 nr__3
parameter (nr__3=2345) ! Number of records
character*765 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*9 obsID (nr__3) ! HST Observation ID (Dataset)
character*8 Mode (nr__3) ! HST observing mode; POSITION or TRANSFER
integer*4 PID (nr__3) ! Raw target ID from proposal (Proposal ID)
integer*4 Plate (nr__3) ! Proposal number of observations
integer*4 Star (nr__3) ! [1/580] Star number of observation from proposal
real*4 pVmag (nr__3) ! (mag) [10/16] Predicted V band magnitude
real*4 Vmag_2 (nr__3) ! (mag) [9.5/16.2] Actual V band magnitude
real*8 RApdeg (nr__3) ! (deg) Predicted Right Ascension (J2000)
real*8 DEpdeg (nr__3) ! (deg) Predicted Declination (J2000)
real*8 RAdeg (nr__3) ! (deg) Right Ascension of V1 telescope axis (J2000)
real*8 DEdeg (nr__3) ! (deg) Declination of V1 telescope axis (J2000)
real*8 rollV3 (nr__3) ! (deg) Roll about V3 of telescope
real*8 Xpos (nr__3) ! (km) X position of telescope
real*8 Ypos (nr__3) ! (km) Y position of telescope
real*8 Zpos (nr__3) ! (km) Z position of telescope
real*8 Xvel (nr__3) ! (km/s) X velocity of telescope
real*8 Yvel (nr__3) ! (km/s) Y velocity of telescope
real*8 Zvel (nr__3) ! (km/s) Z velocity of telescope
real*8 JD (nr__3) ! (d) Observation Julian Date; corrected
real*8 MJD (nr__3) ! (d) Modified Julian observation time
integer*4 Obs_Y (nr__3) ! (yr) [2008/2011] Year of observation
integer*4 Obs_D (nr__3) ! (d) [2/356] Day of observation
integer*4 Obs_h (nr__3) ! (h) [0/23] The hour of observation
integer*4 Astro (nr__3) ! [1] FGS used as the astrometer (always "1")
character*5 Filt (nr__3) ! [F583W] Filter used for observation
* (always F583W)
integer*4 Set (nr__3) ! [1/12] Set number of observation group
integer*4 E_Xvel (nr__3) ! (km/s) [0] Earth plus HST X velocity
integer*4 E_Yvel (nr__3) ! (km/s) [0] Earth plus HST Y velocity
integer*4 E_Zvel (nr__3) ! (km/s) [0] Earth plus HST Z velocity
real*8 Xave (nr__3) ! (arcsec) Uncorrected X position calculated by average
real*8 Yave (nr__3) ! (arcsec) Uncorrected Y position calculated by average
real*8 e_Xave (nr__3) ! (arcsec) [0.001/0.6] Standard deviation in Xave
real*8 e_Yave (nr__3) ! (arcsec) [0.001/0.02] Standard deviation in Yave
real*4 XYcor (nr__3) ! (arcsec) XY correlation
real*4 XYcov (nr__3) ! (arcsec) XY covariance
real*8 Xmed (nr__3) ! (arcsec) Uncorrected X position calculated by median
real*8 Ymed (nr__3) ! (arcsec) Uncorrected Y position calculated by median
real*8 Xmad (nr__3) ! (arcsec) [0.001/0.2] X median average deviation
real*8 Ymad (nr__3) ! (arcsec) [0.0009/0.007] Y median average deviation
real*4 exptm (nr__3) ! (s) [20/124] Seconds of FINELOCK data
integer*4 npix (nr__3) ! [829/4969] Number of samples
* from FINELOCK interval
character*5 VLabSi (nr__3) ! Instrument velocity aberration measured from
* HST/FGS-1r (always FGS1r)
integer*4 VLabFGS (nr__3) ! [1] FGS velocity aberration measured from
real*8 XVLab (nr__3) ! (arcsec) X position of velocity aberration
real*8 YVLab (nr__3) ! (arcsec) Y position of velocity aberration
integer*4 Dom_GS (nr__3) ! [2/3] Dominant FGS guider
real*8 pmtsum (nr__3) ! [40/11825] Sum of pmts
real*8 pmtxa (nr__3) ! pmt from channel xa
real*8 pmtxb (nr__3) ! pmt from channel xb
real*8 pmtya (nr__3) ! pmt from channel ya
real*8 pmtyb (nr__3) ! pmt from channel yb
integer*4 S_bckgd (nr__3) ! [100] Number of samples for background
integer*4 NSamX (nr__3) ! [0/1733] Number of samples for walkdown
* pmt diff/sum in X axis
real*8 xfld_s (nr__3) ! [-0.6/0.3] X pmt diff/sum from fine lock dv
* interval
real*8 xwdd_s (nr__3) ! [-1/0.2] X pmt diff/sum from walkdown interval
real*4 k1x (nr__3) ! [0.01/0.04] Reference X axis S_curve inverse
* slope of Upgren69
real*8 xdiff (nr__3) ! (arcsec) [-1/0.004] X centroid shift in arcseconds;
* add to Xmed
integer*4 NSamY (nr__3) ! [87/4322] Number of samples for walkdown
* pmt diff/sum in Y axis
real*8 yfld_s (nr__3) ! [-0.4/0.02] Y pmt diff/sum from fine lock dv
* interval
real*8 ywdd_s (nr__3) ! [-1/-0.08] Y pmt diff/sum from walkdown interval
real*4 k1y (nr__3) ! [0.01/0.05] Reference Y axis S_curve inverse
* slope of Upgren69
real*8 ydiff (nr__3) ! (arcsec) [-1/0.005] Y centroid shift in arcseconds;
* add to Ymed
real*8 ctX (nr__3) ! (arcsec) Course track X position
real*8 ctY (nr__3) ! (arcsec) Course track Y position
integer*4 ssAmed (nr__3) ! (arcsec) Star selector position; theta A
integer*4 ssBmed (nr__3) ! (arcsec) Star selector position; theta B
real*4 Xvar (nr__3) ! (arcsec) Variance of X
real*4 Yvar (nr__3) ! (arcsec) Variance of Y
real*4 Xmvar (nr__3) ! (arcsec) [/0.02] Median average deviation variance of X
real*4 Ymvar (nr__3) ! (arcsec) Median average deviation variance of Y
real*8 ssA1 (nr__3) ! Predicted star selector position; theta A
real*8 ssB1 (nr__3) ! Predicted star selector position; theta B
real*8 intrhoA (nr__3) ! Interpolated rhoA lever arm
real*8 intkA (nr__3) ! Interpolated kA lever arm
real*8 Q1 (nr__3) ! Spacecraft quaternion
real*8 Q2 (nr__3) ! Spacecraft quaternion
real*8 Q3 (nr__3) ! Spacecraft quaternion
real*8 Q4 (nr__3) ! Spacecraft quaternion
real*8 Palpha (nr__3) ! [-3.5/2.3] Parallax factor {alpha}
real*8 Pdelta (nr__3) ! [-1/1.1] Parallax factor {delta}
real*4 Xsinfit (nr__3) ! (arcsec) [-0.003/0.002] Sinfit position correction for X
real*8 Xfcqc (nr__3) ! (arcsec) Final X corrected position
* with polynomial drift correction
real*8 Yfcqc (nr__3) ! (arcsec) Final Y corrected position
integer*4 HIP_3 (nr__3) ! Target name (HIP number)
C=============================================================================
C Loading file 'table2.dat' ! Metal-poor target star parallaxes
C Format for file interpretation
1 format(I6,1X,F6.3,1X,F5.3,1X,F5.2,1X,F4.2,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,8
read(1,'(A43)')ar__
read(ar__,1)
+ HIP(i__),Plx(i__),e_Plx(i__),PlxH07(i__),e_PlxH07(i__),
+ PlxT17(i__),e_PlxT17(i__)
if(ar__(32:37) .EQ. '') PlxT17(i__) = rNULL__
if(ar__(39:43) .EQ. '') e_PlxT17(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ HIP(i__),Plx(i__),e_Plx(i__),PlxH07(i__),e_PlxH07(i__),
+ PlxT17(i__),e_PlxT17(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Metal-poor target star proper motions
C Format for file interpretation
2 format(
+ I6,1X,F8.3,1X,F5.3,1X,F8.3,1X,F5.3,1X,F7.2,1X,F4.2,1X,F7.2,1X,
+ F4.2,1X,F6.1,1X,F3.1,1X,F6.1,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,8
read(1,'(A84)')ar__1
read(ar__1,2)
+ HIP_1(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__),
+ pmRAH07(i__),e_pmRAH07(i__),pmDEH07(i__),e_pmDEH07(i__),
+ pmRAPP(i__),e_pmRAPP(i__),pmDEPP(i__),e_pmDEPP(i__)
c ..............Just test output...........
write(6,2)
+ HIP_1(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__),
+ pmRAH07(i__),e_pmRAH07(i__),pmDEH07(i__),e_pmDEH07(i__),
+ pmRAPP(i__),e_pmRAPP(i__),pmDEPP(i__),e_pmDEPP(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Summary of parallax star properties
C Format for file interpretation
3 format(
+ I6,1X,F5.2,1X,F5.3,1X,F4.2,1X,F5.3,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,8
read(1,'(A55)')ar__2
read(ar__2,3)
+ HIP_2(i__),Vmag(i__),V_I(i__),E_B_V(i__),v_V_I_0(i__),
+ VMag_1(i__),e_VMag(i__),F606W(i__),v_Fe_H_(i__),v_a_Fe_(i__)
c ..............Just test output...........
write(6,3)
+ HIP_2(i__),Vmag(i__),V_I(i__),E_B_V(i__),v_V_I_0(i__),
+ VMag_1(i__),e_VMag(i__),F606W(i__),v_Fe_H_(i__),v_a_Fe_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! Astrometric data for the 8 targets and
* reference stars
C Format for file interpretation
4 format(
+ A9,1X,A8,1X,I9,1X,I5,1X,I3,1X,F5.2,1X,F5.2,1X,F9.5,1X,F9.5,1X,
+ F8.4,1X,F10.6,1X,F9.5,1X,F12.6,1X,F12.6,1X,F12.6,1X,F9.6,1X,
+ F9.6,1X,F9.6,1X,F11.3,1X,F11.5,1X,I4,1X,I3,7X,I2,1X,I1,1X,A5,
+ 1X,I2,1X,I1,1X,I1,1X,I1,1X,F10.5,1X,F9.5,1X,F7.5,1X,F7.5,1X,
+ E9.2,1X,E9.2,1X,F14.9,1X,F11.7,1X,F7.5,1X,F7.5,1X,E8.2,1X,I4,
+ 1X,A5,1X,I1,1X,F14.9,1X,F11.7,1X,I1,1X,F10.4,1X,F9.4,1X,F9.4,
+ 1X,F9.4,1X,F9.4,1X,I3,1X,I4,1X,F7.4,1X,F7.4,1X,F6.4,1X,F9.6,
+ 1X,I4,1X,F7.4,1X,F7.4,1X,F6.4,1X,F9.6,1X,F8.3,1X,F7.3,1X,I6,
+ 1X,I6,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,1X,F11.4,1X,F11.4,
+ 1X,F11.9,1X,F11.9,1X,F12.9,1X,F12.9,1X,F12.9,1X,F12.9,1X,
+ F12.9,1X,F12.9,1X,E12.5,1X,F14.9,1X,F11.7,4X,I6)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,2345
read(1,'(A765)')ar__3
read(ar__3,4)
+ obsID(i__),Mode(i__),PID(i__),Plate(i__),Star(i__),pVmag(i__),
+ Vmag_2(i__),RApdeg(i__),DEpdeg(i__),RAdeg(i__),DEdeg(i__),
+ rollV3(i__),Xpos(i__),Ypos(i__),Zpos(i__),Xvel(i__),Yvel(i__),
+ Zvel(i__),JD(i__),MJD(i__),Obs_Y(i__),Obs_D(i__),Obs_h(i__),
+ Astro(i__),Filt(i__),Set(i__),E_Xvel(i__),E_Yvel(i__),
+ E_Zvel(i__),Xave(i__),Yave(i__),e_Xave(i__),e_Yave(i__),
+ XYcor(i__),XYcov(i__),Xmed(i__),Ymed(i__),Xmad(i__),Ymad(i__),
+ exptm(i__),npix(i__),VLabSi(i__),VLabFGS(i__),XVLab(i__),
+ YVLab(i__),Dom_GS(i__),pmtsum(i__),pmtxa(i__),pmtxb(i__),
+ pmtya(i__),pmtyb(i__),S_bckgd(i__),NSamX(i__),xfld_s(i__),
+ xwdd_s(i__),k1x(i__),xdiff(i__),NSamY(i__),yfld_s(i__),
+ ywdd_s(i__),k1y(i__),ydiff(i__),ctX(i__),ctY(i__),ssAmed(i__),
+ ssBmed(i__),Xvar(i__),Yvar(i__),Xmvar(i__),Ymvar(i__),
+ ssA1(i__),ssB1(i__),intrhoA(i__),intkA(i__),Q1(i__),Q2(i__),
+ Q3(i__),Q4(i__),Palpha(i__),Pdelta(i__),Xsinfit(i__),
+ Xfcqc(i__),Yfcqc(i__),HIP_3(i__)
c ..............Just test output...........
write(6,4)
+ obsID(i__),Mode(i__),PID(i__),Plate(i__),Star(i__),pVmag(i__),
+ Vmag_2(i__),RApdeg(i__),DEpdeg(i__),RAdeg(i__),DEdeg(i__),
+ rollV3(i__),Xpos(i__),Ypos(i__),Zpos(i__),Xvel(i__),Yvel(i__),
+ Zvel(i__),JD(i__),MJD(i__),Obs_Y(i__),Obs_D(i__),Obs_h(i__),
+ Astro(i__),Filt(i__),Set(i__),E_Xvel(i__),E_Yvel(i__),
+ E_Zvel(i__),Xave(i__),Yave(i__),e_Xave(i__),e_Yave(i__),
+ XYcor(i__),XYcov(i__),Xmed(i__),Ymed(i__),Xmad(i__),Ymad(i__),
+ exptm(i__),npix(i__),VLabSi(i__),VLabFGS(i__),XVLab(i__),
+ YVLab(i__),Dom_GS(i__),pmtsum(i__),pmtxa(i__),pmtxb(i__),
+ pmtya(i__),pmtyb(i__),S_bckgd(i__),NSamX(i__),xfld_s(i__),
+ xwdd_s(i__),k1x(i__),xdiff(i__),NSamY(i__),yfld_s(i__),
+ ywdd_s(i__),k1y(i__),ydiff(i__),ctX(i__),ctY(i__),ssAmed(i__),
+ ssBmed(i__),Xvar(i__),Yvar(i__),Xmvar(i__),Ymvar(i__),
+ ssA1(i__),ssB1(i__),intrhoA(i__),intkA(i__),Q1(i__),Q2(i__),
+ Q3(i__),Q4(i__),Palpha(i__),Pdelta(i__),Xsinfit(i__),
+ Xfcqc(i__),Yfcqc(i__),HIP_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end