Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/490/3508 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-14
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/MNRAS/490/3508 The Southern Stellar Stream Sp. Survey (S^5^) (Li+, 2019)
*================================================================================
*The Southern Stellar Stream Spectroscopic Survey (S^5^): overview, target
*selection, data reduction, validation, and early science.
* Li T.S., Koposov S.E., Zucker D.B., Lewis G.F., Kuehn K., Simpson J.D.,
* Ji A.P., Shipp N., Mao Y.-Y., Geha M., Pace A.B., Mackey A.D., Allam S.,
* Tucker D.L., Da Costa G.S., Erkal D., Simon J.D., Mould J.R., Martell S.L.,
* Wan Z., De Silva G.M., Bechtol K., Balbinot E., Belokurov V.,
* Bland-Hawthorn J., Casey A.R., Cullinane L., Drlica-Wagner A., Sharma S.,
* Vivas A.K., Wechsler R.H., Yanny B., (the S5 Collaboration)
* <Mon. Not. R. Astron. Soc., 490, 3508-3531 (2019)>
* =2019MNRAS.490.3508L (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 's5dr1.dat' ! First Public Data Release of the Southern Stellar
Stream Spectroscopic Survey (S^5^ DR1)
integer*4 nr__
parameter (nr__=31777) ! Number of records
character*799 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 vel1 (nr__) ! (km/s) [-2000/1999] Radial velocity
* statistic from rvspecfit 1st
* percentile of the RV posterior
* distribution (1)
real*8 vel99 (nr__) ! (km/s) [-1990/2000] Radial velocity
* statistic from rvspecfit 99th
* percentile of the RV posterior
* distribution (1)
real*8 vel50 (nr__) ! (km/s) [-1998/2000] Radial velocity
* statistic from rvspecfit 50th
* percentile of the RV posterior
* distribution (1)
real*8 vel16 (nr__) ! (km/s) [-2000/2000] Radial velocity
* statistic from rvspecfit 16th
* percentile of the RV posterior
* distribution (1)
real*8 vel84 (nr__) ! (km/s) [-1996/2000] Radial velocity
* statistic from rvspecfit 84th
* percentile of the RV posterior
* distribution (1)
real*8 vel_std (nr__) ! (km/s) [0.027/228] Standard deviation of the
* distribution
real*4 vel_skew (nr__) ! (km/s) [-7.3/5.7] Vel skewness
real*8 vel_kurt (nr__) ! (km/s) [-1.9/120] Vel kurtosis
real*8 feh1 (nr__) ! ([-]) [-4.9/1.1] Metallicity statistic from
* rvspecfit, 1st percentile of the
* distribution
real*8 feh99 (nr__) ! ([-]) [-4.3/1.9] Metallicity statistic from
* rvspecfit, 99th percentile of the
* distribution
real*8 feh50 (nr__) ! ([-]) [-4.4/1.5] Metallicity statistic from
* rvspecfit, 50th percentile of the
* distribution (2)
real*8 feh16 (nr__) ! ([-]) [-4.5/1.4] Metallicity statistic from
* rvspecfit, 16th percentile of the
* distribution
real*8 feh84 (nr__) ! ([-]) [-4.4/1.9] Metallicity statistic from
* rvspecfit, 84th percentile of the
* distribution
real*4 feh_std (nr__) ! ([-]) [1.9e-5/2.2] Metallicity standard
* deviation
real*8 feh_skew (nr__) ! ([-]) [-9.2/8.7] feh skewness
real*8 feh_kurt (nr__) ! ([-]) [-2/95] feh kurtosis
real*8 logg1 (nr__) ! ([cm/s2]) [-1.5/6.4] Log of surface gravity
* from rvspecfit, 1st percentile
real*8 logg99 (nr__) ! ([cm/s2]) [-0.9/7.6] Log of surface gravity
* from rvspecfit, 99th percentile
real*8 logg50 (nr__) ! ([cm/s2]) [-1.3/7.3] Log of surface gravity
* from rvspecfit, 50th percentile
real*8 logg16 (nr__) ! ([cm/s2]) [-1.4/7] Log of surface gravity from
* rvspecfit, 16th percentile
real*8 logg84 (nr__) ! ([cm/s2]) [-1/7.3] Log of surface gravity from
* rvspecfit, 84th percentile
real*8 logg_std (nr__) ! ([cm/s2]) [0.0001/3.1] logg standard deviation
real*8 logg_skew (nr__) ! ([cm/s2]) [-15.4/12.3] logg skewness
real*8 logg_kurt (nr__) ! ([cm/s2]) [-2/336] logg kurtosis
real*8 teff1 (nr__) ! (K) [1655/9893] Effective temperature
* from rvspecfit, 1st percentile
real*8 teff99 (nr__) ! (K) [2070/18973] Effective temperature
* from rvspecfit, 99th percentile
real*8 teff50 (nr__) ! (K) [1781/12514] Effective temperature
* from rvspecfit, 50th percentile
real*8 teff16 (nr__) ! (K) [1774/11514] Effective temperature
* from rvspecfit, 16th percentile
real*8 teff84 (nr__) ! (K) [1790/14340] Effective temperature
* from rvspecfit, 84th percentile
real*8 teff_std (nr__) ! (K) [0.07/3238] teff standard deviation
real*4 teff_skew (nr__) ! (K) [-7.4/11.2] teff skewness
real*8 teff_kurt (nr__) ! (K) [-2/161] teff kurtosis
real*8 velcalib (nr__) ! (km/s) [-1999/1999] Radial velocity
* calibrated to have no zeropoint
* offset (vel_calib) (3)
real*4 velcalib_std(nr__) ! (km/s) [0.6/294] Calibrated radial velocity
* uncertainty using repeated
* observations (vel_calib_std)
real*4 fehcalib_std(nr__) ! [2.5e-5/2.8] Calibrated [Fe/H]
* uncertainty using repeated
* observations (feh_calib_std)
real*8 rchi2_1700d(nr__) ! [0.8/481] Reduced chi square from the
* best fit template model from
* rvspecfit (redchisq_1700d)
real*8 chi2_1700d (nr__) ! [1579/956425] Chi square from the
* best fit template model from
* rvspecfit (chisq_1700d)
real*4 absdev_1700d(nr__) ! [0.006/0.5] Median absolute deviation
* from the best fit template model from
* rvspecfit (4)
real*8 RAdeg (nr__) ! (deg) Right ascension in decimal degrees
* (J2000) (ra)
real*8 DEdeg (nr__) ! (deg) [-77/13.2] Declination in decimal
* degrees (J2000) (dec)
character*19 Name (nr__) ! Object name used for AAT observation
* (5)
character*19 Field (nr__) ! Name of the stream fields
* as in Table2 (6)
character*50 FileName (nr__) ! Filename of the reduced AAT spectra
* from each field - files not released
* in DR1 (fileinfo)
real*8 MJD (nr__) ! [58141.4/58579.8] UT Modified Julian
* Date of the observation (utmjd)
real*4 sn_1700d (nr__) ! [3/228] Median Signal-to-Noise ratio
* of the 1700d spectra (7)
real*8 MoonDist (nr__) ! (deg) [12/179] Moon distance (moon_dist)
real*4 MoonPhase (nr__) ! [0.04/100] Moon phase
* (0-100% fraction of surface
* illumination; 0%=new moon -
* 100%=full moon) (moon_phase)
real*8 MoonAlt (nr__) ! (deg) [-78/61] Moon altitude (moon_alt)
integer*4 Flag (nr__) ! [0] QSO flag based on WISE
* photometry; always 0 = not a QSO
* (see Section 4.4) (qso_flag_wise)
integer*4 P (nr__) ! [-1/9] Initial target selection
* priorities (priority) (8)
integer*4 Prim (nr__) ! [0/1] For some object, multiple
* observations were taken (primary) (9)
real*4 Prob (nr__) ! [0.5/1] Good star probability based
* on random forest (good_star_pb) (10)
integer*4 good_star (nr__) ! [1] A flag to select stellar spectra
* with good fitting quality
* (always 1=good) (11)
real*4 EBV (nr__) ! (mag) [0.006/1.4]? E(B-V) from
* Schlegel+ (1998ApJ...500..525S)
integer*8 Gaia (nr__) ! ? Gaia EDR3 identifier
* (gaia_source_id)
real*4 plx (nr__) ! (mas) [-5.4/9.8]? Gaia EDR3 Parallax
* (parallax)
real*4 e_plx (nr__) ! (mas) [0.01/1.3]? Parallax uncertainty
* (parallax_error)
real*8 pmRA (nr__) ! (mas/yr) [-60/131]? Gaia EDR3 Proper motion in
* right ascension
real*8 pmDE (nr__) ! (mas/yr) [-113/41]? Gaia EDR3 Proper motion in
* declination (pmdec)
real*4 e_pmRA (nr__) ! (mas/yr) [0.01/1.4]? pmRA uncertainty
* (pmra_error)
real*4 e_pmDE (nr__) ! (mas/yr) [0.009/1.3]? pmDE uncertainty
* (pmdec_error)
real*8 Gmag (nr__) ! (mag) [11.99/21.61]? Gaia EDR3 G-band
* magnitude (phot_g_mean_mag)
real*8 BPmag (nr__) ! (mag) [12.2/20.8]? Gaia EDR3 BP-band
* magnitude (phot_bp_mean_mag)
real*8 RPmag (nr__) ! (mag) [10.87/19.9]? Gaia EDR3 RP-band
* magnitude (phot_rp_mean_mag)
real*8 E_BP_RP (nr__) ! [0.86/25.5]? Gaia EDR3 BP-RP excess
* factor (phot_bp_rp_excess_factor)
real*8 gmag_1 (nr__) ! (mag) [11.6/21.3]? DECam g-band magnitude
* (decam_g) (12)
real*8 rmag (nr__) ! (mag) [11.7/20.3]? DECam r-band magnitude
* (decam_r) (12)
real*8 imag (nr__) ! (mag) [11.3/20.2]? DECam i-band magnitude
* (decam_i) (12)
real*8 zmag (nr__) ! (mag) [10.7/20.6]? DECam z-band magnitude
* (decam_z) (12)
integer*4 r_DECam (nr__) ! [0/1755] The source of the decam
* photometry (decam_phot_src_bits) (13)
real*4 FW1 (nr__) ! (309.05nJy) [33.5/1.9e+6]? UnWISE W1 (3.4um) band
* flux density in "Vega nMgy" (flux_w1)
real*4 FW2 (nr__) ! (167.66nJy) [93.6/2.4e+6]? UnWISE W2 (4.6um) band
* flux density in "Vega nMgy" (flux_w2)
real*4 e_FW1 (nr__) ! (309.05nJy) [4.5/1233]? FW1 uncertainty (dflux_w1)
real*4 e_FW2 (nr__) ! (167.66nJy) [17/2856]? FW2 uncertainty (dflux_w2)
real*8 umagS (nr__) ! (mag) [13.3/20]? Skymapper DR2 u-band
* magnitude (u_skm)
real*8 gmagS (nr__) ! (mag) [12/20.8]? Skymapper DR2 g-band
* magnitude (g_skm)
real*8 rmagS (nr__) ! (mag) [11.8/20]? Skymapper DR2 r-band
* magnitude (r_skm)
real*8 imagS (nr__) ! (mag) [11/20.1]? Skymapper DR2 i-band
* magnitude (i_skm)
real*8 zmagS (nr__) ! (mag) [10.7/19.7]? Skymapper DR2 z-band
* magnitude (z_skm)
real*8 Dist05 (nr__) ! (kpc) [0.1/54]? Starhorse distance at 5th
* percentile
* from Queiroz+ 2018MNRAS.476.2556Q
* (starhorse_dist05)
real*8 Dist16 (nr__) ! (kpc) [0.1/67]? Starhorse distance at 16th
* percentile
* from Queiroz+ 2018MNRAS.476.2556Q
* (starhorse_dist16)
real*8 Dist50 (nr__) ! (kpc) [0.1/78]? Starhorse distance at 50th
* percentile
* from Queiroz+ 2018MNRAS.476.2556Q
* (starhorse_dist50)
real*8 Dist84 (nr__) ! (kpc) [0.1/97]? Starhorse distance at 84th
* percentile
* from Queiroz+ 2018MNRAS.476.2556Q
* (starhorse_dist84)
real*8 Dist95 (nr__) ! (kpc) [0.1/120]? Starhorse distance at 95th
* percentile
* from Queiroz+ 2018MNRAS.476.2556Q
* (starhorse_dist95)
*Note (1): For radial velocity, we recommend the velcalib column.
*Note (2): feh50 is recommended to be used as mean metallicity for each star.
*Note (3): Offset is determined by cross-matching with APOGEE data.
* See details in Section 4.1
*Note (4): Median absolute deviation from the best fit template model from
* rvspecfit, defined as Median(Data-Model)/Median(Model)
*Note (5): For stellar targets, name is usually source_id from Gaia DR2, for
* galaxy targets (S5-lowz), the name is usually object_id from DES DR1.
* Some observations prior S5 were also included and for those, names
* could follow different conventions.
*Note (6): The field names given in S^5^ DR1 were made early on and can
* differ from those given in Table 2 of the paper. They have been
* homogenized by the CDS.
*Note (7): For DR1, we decide to only release objects with sn_1700d>3 to ensure
* the high quality of this dataset.
*Note (8): Priority target are:
* -1 = few fields were observed prior to the official S5 observations
* (2528 occurrences)
* 7-9 = stream targets
* 3-6 = halo star targets
*
* For details, refer to Section 2.
*Note (9): Primary = True ("1") if this is the observation with the highest S/N
* (29,584 occurrences). For objects that only one observation was
* taken, Primary is always True. This column is useful if you want to
* get only one measurement per target to remove duplicates.
*Note (10): Good star probability based on random forest trained on the following
* parameters as features: chi-squared values of the fit, radial
* velocity error, radial velocity posterior skewness and kurtosis,
* effective temperature, median signal to noise in the spectrum and
* median absolute deviation from the best fit template model.
* See details in Section 4.5.
*Note (11): A flag to select stellar spectra with good fitting quality.
* good_star=0 if it meets one of the 3 following criteria:
* 1. good_star_pb<0.5
* 2. priority<=2
* 3. qso_flag_wise==1
* See details in Section 4.5.
*
* For DR1, we decide to only release objects with good_star=1 to
* ensure the high quality of this dataset.
*Note (12): DECam photometry (without reddening correction) from DES DR1,
* DECaLS DR9 or NSC DR2.
* If the object is in multiple catalogs, then DES is first considered,
* then DECaLS and finally All-sky NOAO Source Catalog (NSC).
* See r_DECam for the source of the photometry
*Note (13): The source of the decam photometry. It is a 12 bit integer
* "ABC DEF GHI JKL" where ABC bits describe the source of g photometry,
* DEF bits the source of r photometry, etc. The value of 1 corresponds
* to DES DR2, 2 to DECALS, 3 to NSC. For example, r_DECam=585
* corresponds to 001 001 001 001, which means that griz photometry all
* comes from DES DR2; r_DECam=83 corresponds to 000 001 010 011, which
* means that the griz photometry comes from NONE, DES, DECALS, NSC.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! S^5^ fields observed with the Anglo-Australian
Telescope (AAT) as of 2019 June
integer*4 nr__1
parameter (nr__1=115) ! Number of records
character*79 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*19 Field_1 (nr__1) ! S^5^ field name (1)
real*8 RAdeg_1 (nr__1) ! (deg) Centre of the field right ascension (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) [-76/12.2] Centre of the field
* declination (J2000)
character*10 Date (nr__1) ! ("Y/M/D") UT date of the observation (with highest S/N
* if observed multiple times)
real*8 MJD_1 (nr__1) ! (d) [58141.4/58579.7] Modified Julian Date (start)
* of the observation
integer*4 Texp (nr__1) ! (s) [600/9600] Total exposure time
real*4 Gmag_2 (nr__1) ! (mag) [14.7/19.3] Average Gaia DR2 G-band magnitude
* at S/N=5 per pixel (from red arm spectra)
integer*4 Nt (nr__1) ! [96/367] Total number of targets
integer*4 Ns (nr__1) ! [67/351] Number of stars (N_goodstar_) with
* good measurements (2)
*Note (1): All fields are grouped into four categories, which are fields in
* the DES footprint, fields outside the DES footprint (see Section 2.6),
* calibration fields for survey validation (see Section 4), and fields
* observed prior to S5 with the same instrument setup but previously
* unpublished.
*Note (2): i.e. good_star=1; see definition in Section 4.5.
* We note that for fields in some streams such as ATLAS, Elqui, Phoenix,
* etc., N_goodstar_ is usually much lower than Ntargets, because these
* streams are at high Galactic latitude; therefore, we used ~100 spare
* fibres to target low redshift galaxies (see Section 2.4), and all
* galaxy targets are assigned good_star=0 regardless of the quality of
* the spectra.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! A list of 674 QSOs identified by S^5^ with
robust redshift measurements
integer*4 nr__2
parameter (nr__2=674) ! Number of records
character*51 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*8 Gaia_1 (nr__2) ! Gaia DR2 source ID
real*8 RAdeg_2 (nr__2) ! (deg) [4.69/358.13] Right ascension (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) [-64.76/-18.59] Declination (J2000)
real*4 Gmag_3 (nr__2) ! (mag) [15.1/20.3] Gaia DR2 G-band magnitude
real*4 z (nr__2) ! [0.04/3.4] Measured redshift
C=============================================================================
C Loading file 's5dr1.dat' ! First Public Data Release of the Southern Stellar
* Stream Spectroscopic Survey (S^5^ DR1)
C Format for file interpretation
1 format(
+ F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F7.3,1X,F6.3,1X,F7.3,
+ 1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,E9.4,1X,F7.4,1X,
+ F7.4,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F7.5,1X,F9.5,
+ 1X,F9.5,1X,F8.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F8.3,1X,
+ F6.3,1X,F7.3,1X,F8.2,1X,F6.2,1X,E10.5,1X,F7.3,1X,F10.3,1X,
+ F6.4,1X,F11.7,1X,F11.7,1X,A19,1X,A19,1X,A50,1X,F11.5,1X,F6.2,
+ 1X,F10.6,1X,F6.3,1X,F9.5,1X,I1,1X,I2,1X,I1,1X,F5.3,1X,I1,1X,
+ F6.4,1X,I19,1X,F6.3,1X,F5.3,1X,F7.3,1X,F8.3,1X,F5.3,1X,F5.3,
+ 1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,I4,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F7.4,1X,
+ F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,
+ 1X,F10.6)
C Effective file loading
open(unit=1,status='old',file=
+'s5dr1.dat')
write(6,*) '....Loading file: s5dr1.dat'
do i__=1,31777
read(1,'(A799)')ar__
read(ar__,1)
+ vel1(i__),vel99(i__),vel50(i__),vel16(i__),vel84(i__),
+ vel_std(i__),vel_skew(i__),vel_kurt(i__),feh1(i__),feh99(i__),
+ feh50(i__),feh16(i__),feh84(i__),feh_std(i__),feh_skew(i__),
+ feh_kurt(i__),logg1(i__),logg99(i__),logg50(i__),logg16(i__),
+ logg84(i__),logg_std(i__),logg_skew(i__),logg_kurt(i__),
+ teff1(i__),teff99(i__),teff50(i__),teff16(i__),teff84(i__),
+ teff_std(i__),teff_skew(i__),teff_kurt(i__),velcalib(i__),
+ velcalib_std(i__),fehcalib_std(i__),rchi2_1700d(i__),
+ chi2_1700d(i__),absdev_1700d(i__),RAdeg(i__),DEdeg(i__),
+ Name(i__),Field(i__),FileName(i__),MJD(i__),sn_1700d(i__),
+ MoonDist(i__),MoonPhase(i__),MoonAlt(i__),Flag(i__),P(i__),
+ Prim(i__),Prob(i__),good_star(i__),EBV(i__),Gaia(i__),
+ plx(i__),e_plx(i__),pmRA(i__),pmDE(i__),e_pmRA(i__),
+ e_pmDE(i__),Gmag(i__),BPmag(i__),RPmag(i__),E_BP_RP(i__),
+ gmag_1(i__),rmag(i__),imag(i__),zmag(i__),r_DECam(i__),
+ FW1(i__),FW2(i__),e_FW1(i__),e_FW2(i__),umagS(i__),gmagS(i__),
+ rmagS(i__),imagS(i__),zmagS(i__),Dist05(i__),Dist16(i__),
+ Dist50(i__),Dist84(i__),Dist95(i__)
if(ar__(516:521) .EQ. '') EBV(i__) = rNULL__
if(ar__(523:541) .EQ. '') Gaia(i__) = iNULL__
if(ar__(543:548) .EQ. '') plx(i__) = rNULL__
if(ar__(550:554) .EQ. '') e_plx(i__) = rNULL__
if(ar__(556:562) .EQ. '') pmRA(i__) = rNULL__
if(ar__(564:571) .EQ. '') pmDE(i__) = rNULL__
if(ar__(573:577) .EQ. '') e_pmRA(i__) = rNULL__
if(ar__(579:583) .EQ. '') e_pmDE(i__) = rNULL__
if(ar__(585:591) .EQ. '') Gmag(i__) = rNULL__
if(ar__(593:599) .EQ. '') BPmag(i__) = rNULL__
if(ar__(601:607) .EQ. '') RPmag(i__) = rNULL__
if(ar__(609:615) .EQ. '') E_BP_RP(i__) = rNULL__
if(ar__(617:624) .EQ. '') gmag_1(i__) = rNULL__
if(ar__(626:633) .EQ. '') rmag(i__) = rNULL__
if(ar__(635:642) .EQ. '') imag(i__) = rNULL__
if(ar__(644:651) .EQ. '') zmag(i__) = rNULL__
if(ar__(658:669) .EQ. '') FW1(i__) = rNULL__
if(ar__(671:682) .EQ. '') FW2(i__) = rNULL__
if(ar__(684:695) .EQ. '') e_FW1(i__) = rNULL__
if(ar__(697:708) .EQ. '') e_FW2(i__) = rNULL__
if(ar__(710:716) .EQ. '') umagS(i__) = rNULL__
if(ar__(718:724) .EQ. '') gmagS(i__) = rNULL__
if(ar__(726:732) .EQ. '') rmagS(i__) = rNULL__
if(ar__(734:740) .EQ. '') imagS(i__) = rNULL__
if(ar__(742:748) .EQ. '') zmagS(i__) = rNULL__
if(ar__(750:758) .EQ. '') Dist05(i__) = rNULL__
if(ar__(760:768) .EQ. '') Dist16(i__) = rNULL__
if(ar__(770:778) .EQ. '') Dist50(i__) = rNULL__
if(ar__(780:788) .EQ. '') Dist84(i__) = rNULL__
if(ar__(790:799) .EQ. '') Dist95(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ vel1(i__),vel99(i__),vel50(i__),vel16(i__),vel84(i__),
+ vel_std(i__),vel_skew(i__),vel_kurt(i__),feh1(i__),feh99(i__),
+ feh50(i__),feh16(i__),feh84(i__),feh_std(i__),feh_skew(i__),
+ feh_kurt(i__),logg1(i__),logg99(i__),logg50(i__),logg16(i__),
+ logg84(i__),logg_std(i__),logg_skew(i__),logg_kurt(i__),
+ teff1(i__),teff99(i__),teff50(i__),teff16(i__),teff84(i__),
+ teff_std(i__),teff_skew(i__),teff_kurt(i__),velcalib(i__),
+ velcalib_std(i__),fehcalib_std(i__),rchi2_1700d(i__),
+ chi2_1700d(i__),absdev_1700d(i__),RAdeg(i__),DEdeg(i__),
+ Name(i__),Field(i__),FileName(i__),MJD(i__),sn_1700d(i__),
+ MoonDist(i__),MoonPhase(i__),MoonAlt(i__),Flag(i__),P(i__),
+ Prim(i__),Prob(i__),good_star(i__),EBV(i__),Gaia(i__),
+ plx(i__),e_plx(i__),pmRA(i__),pmDE(i__),e_pmRA(i__),
+ e_pmDE(i__),Gmag(i__),BPmag(i__),RPmag(i__),E_BP_RP(i__),
+ gmag_1(i__),rmag(i__),imag(i__),zmag(i__),r_DECam(i__),
+ FW1(i__),FW2(i__),e_FW1(i__),e_FW2(i__),umagS(i__),gmagS(i__),
+ rmagS(i__),imagS(i__),zmagS(i__),Dist05(i__),Dist16(i__),
+ Dist50(i__),Dist84(i__),Dist95(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! S^5^ fields observed with the Anglo-Australian
* Telescope (AAT) as of 2019 June
C Format for file interpretation
2 format(
+ A19,1X,F10.6,1X,F10.6,1X,A10,1X,F8.2,1X,I4,1X,F4.1,1X,I3,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,115
read(1,'(A79)')ar__1
read(ar__1,2)
+ Field_1(i__),RAdeg_1(i__),DEdeg_1(i__),Date(i__),MJD_1(i__),
+ Texp(i__),Gmag_2(i__),Nt(i__),Ns(i__)
c ..............Just test output...........
write(6,2)
+ Field_1(i__),RAdeg_1(i__),DEdeg_1(i__),Date(i__),MJD_1(i__),
+ Texp(i__),Gmag_2(i__),Nt(i__),Ns(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb1.dat' ! A list of 674 QSOs identified by S^5^ with
* robust redshift measurements
C Format for file interpretation
3 format(I19,1X,F10.6,1X,F10.6,1X,F4.1,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,674
read(1,'(A51)')ar__2
read(ar__2,3)
+ Gaia_1(i__),RAdeg_2(i__),DEdeg_2(i__),Gmag_3(i__),z(i__)
c ..............Just test output...........
write(6,3)
+ Gaia_1(i__),RAdeg_2(i__),DEdeg_2(i__),Gmag_3(i__),z(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end