Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/505/5567 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-Oct-07
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/505/5567 Red Be Stars in the LMC and the SMC (Vieira+, 2021)
*================================================================================
*Red Be stars in the Magellanic Clouds?
* Vieira K., Garcia-Varela A., Sabogal B., Rimulo L.R., Hernandez J.
* <Mon. Not. R. Astron. Soc., 505, 5567-5583 (2021)>
* =2021MNRAS.505.5567V (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Sample of the LMC Be stars data used
in this study
integer*4 nr__
parameter (nr__=2393) ! Number of records
character*273 ar__ ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=2015.5)
integer*4 ID (nr__) ! Identifier number index from this work (idBe)
character*16 OGLE (nr__) ! Designation HHMMSSss-DDMMSSs (J2000)
* (OGLE) (G1)
real*4 Vmag (nr__) ! (mag) Apparent magnitude in V band (Vmag) (G1)
real*4 B_V (nr__) ! (mag) B-V colour index (B-V) (G1)
real*4 V_I (nr__) ! (mag) V-I colour index (V-I) (G1)
real*4 Type (nr__) ! [1/4] Type of the star (Type_num) (G2)
integer*8 GaiaDR2 (nr__) ! Gaia DR2 source identifier
* (GaiaDR2_source_id) (G3)
real*8 RAdeg (nr__) ! (deg) Barycentric right ascension (ICRS)
* at Ep=2015.5 (ra) (G3)
real*8 DEdeg (nr__) ! (deg) Barycentric declination (ICRS)
* at Ep=2015.5 (dec) (G3)
real*8 plx (nr__) ! (mas) ?=- Absolute stellar parallax (parallax) (G3)
real*4 e_plx (nr__) ! (mas) ?=- Mean error of plx (parallax_error) (G3)
real*8 pmRA (nr__) ! (mas/yr) ?=- Proper motion in right ascension direction
* pmRA*cosDE (pmra) (G3)
real*4 e_pmRA (nr__) ! (mas/yr) ?=- Mean error of pmRA (pmra_error) (G3)
real*8 pmDE (nr__) ! (mas/yr) ?=- Proper motion in declination direction
* (pmdec) (G3)
real*4 e_pmDE (nr__) ! (mas/yr) ?=- Mean error of pmDE (pmdec_error) (G3)
real*8 Gmag (nr__) ! (mag) G-band mean magnitude (Vega)
* (phot_g_mean_mag) (G3)
real*8 BPmag (nr__) ! (mag) ?=- Integrated BP mean magnitude (Vega)
* (phot_bp_mean_mag) (G3)
real*8 RPmag (nr__) ! (mag) ?=- Integrated RP mean magnitude (Vega)
* (phot_rp_mean_mag) (G3)
real*8 Rank (nr__) ! Position mag ranking value (rank) (G4)
real*8 Chi (nr__) ! ?=- {Chi}-value for membership to the Clouds
* ({Chi}) (G5)
character*13 ONames (nr__) ! BEL crossmatch confirmed with linked catalogues
* (other_names) (G6)
real*4 Jmag (nr__) ! (mag) ?=- Apparent magnitude in J band (Jmag) (G7)
real*4 Hmag (nr__) ! (mag) ?=- Apparent magnitude in H band (Hmag) (G7)
real*4 Ksmag (nr__) ! (mag) ?=- Apparent magnitude in Ks band (Ksmag) (G7)
real*4 v3_6mag (nr__) ! (mag) ?=- Apparent magnitude Spitzer/IRAC at 3.6
* micron band magnitude ([3.6]) (G8)
real*4 v4_5mag (nr__) ! (mag) ?=- Apparent magnitude Spitzer/IRAC
* at 4.5 micron band magnitude ([4.5]) (G8)
real*4 v5_8mag (nr__) ! (mag) ?=- Apparent magnitude Spitzer/IRAC
* at 5.8 micron band magnitude ([5.8]) (G8)
real*4 v8_0mag (nr__) ! (mag) ?=- Apparent magnitude Spitzer/IRAC
* at 8.0 micron band magnitude ([8.0]) (G8)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Sample of the four exceptionally red Be stars,
ID = 909 is in the LMC and ID = 75, 184,
and 547 are in the SMC
integer*4 nr__1
parameter (nr__1=4) ! Number of records
character*269 ar__1 ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=2015.5)
integer*4 ID_1 (nr__1) ! Identifier number index from this work (idBe)
character*18 OGLE_1 (nr__1) ! Designation HHMMSSss-DDMMSSs (J2000)
* (OGLE) (G1)
real*4 Vmag_1 (nr__1) ! (mag) Apparent magnitude in V band (Vmag) (G1)
real*4 B_V_1 (nr__1) ! (mag) B-V colour index (B-V) (G1)
real*4 V_I_1 (nr__1) ! (mag) V-I colour index (V-I) (G1)
real*4 Type_1 (nr__1) ! [1/4] Type of the star (Type_num) (G2)
integer*8 GaiaDR2_1 (nr__1) ! Gaia DR2 source identifier
* (GaiaDR2_source_id) (G3)
real*8 RAdeg_1 (nr__1) ! (deg) Barycentric right ascension (ICRS)
* at Ep=2015.5 (ra) (G3)
real*8 DEdeg_1 (nr__1) ! (deg) Barycentric declination (ICRS)
* at Ep=2015.5 (dec) (G3)
real*8 plx_1 (nr__1) ! (mas) Absolute stellar parallax (parallax) (G3)
real*4 e_plx_1 (nr__1) ! (mas) Mean error of plx (parallax_error) (G3)
real*4 pmRA_1 (nr__1) ! (mas/yr) Proper motion in right ascension direction
* pmRA*cosDE (pmra) (G3)
real*4 e_pmRA_1 (nr__1) ! (mas/yr) Mean error of pmRA (pmra_error) (G3)
real*4 pmDE_1 (nr__1) ! (mas/yr) Proper motion in declination direction
* (pmdec) (G3)
real*4 e_pmDE_1 (nr__1) ! (mas/yr) Mean error of pmDE (pmdec_error) (G3)
real*8 Gmag_1 (nr__1) ! (mag) G-band mean magnitude (Vega)
* (phot_g_mean_mag) (G3)
real*8 BPmag_1 (nr__1) ! (mag) Integrated BP mean magnitude (Vega)
* (phot_bp_mean_mag) (G3)
real*8 RPmag_1 (nr__1) ! (mag) Integrated RP mean magnitude (Vega)
* (phot_rp_mean_mag) (G3)
real*8 Rank_1 (nr__1) ! Position mag ranking value (rank) (G4)
real*8 Chi_1 (nr__1) ! {Chi}-value for membership to the Clouds
* ({Chi}) (G5)
character*14 ONames_1 (nr__1) ! BEL crossmatch confirmed with linked
* catalogues (other_names) (G6)
real*4 Jmag_1 (nr__1) ! (mag) Apparent magnitude in J band (Jmag) (G7)
real*4 Hmag_1 (nr__1) ! (mag) Apparent magnitude in H band (Hmag) (G7)
real*4 Ksmag_1 (nr__1) ! (mag) Apparent magnitude in Ks band (Ksmag) (G7)
real*4 v3_6mag_1 (nr__1) ! (mag) Apparent magnitude Spitzer/IRAC
* at 3.6 micron band magnitude ([3.6]) (G8)
real*4 v4_5mag_1 (nr__1) ! (mag) Apparent magnitude Spitzer/IRAC
* at 4.5 micron band magnitude ([4.5]) (G8)
real*4 v5_8mag_1 (nr__1) ! (mag) Apparent magnitude Spitzer/IRAC
* at 5.8 micron band magnitude ([5.8]) (G8)
real*4 v8_0mag_1 (nr__1) ! (mag) Apparent magnitude Spitzer/IRAC
* at 8.0 micron band magnitude ([8.0]) (G8)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea3.dat' ! Sample of the SMC Be stars data used
in this study
integer*4 nr__2
parameter (nr__2=1004) ! Number of records
character*317 ar__2 ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=2015.5)
integer*4 ID_2 (nr__2) ! Identifier number index from this work (idBe)
character*18 OGLE_2 (nr__2) ! Designation HHMMSSss-DDMMSSs (J2000)
* (OGLE) (G1)
real*4 Vmag_2 (nr__2) ! (mag) Apparent magnitude in V band (Vmag) (G1)
real*4 B_V_2 (nr__2) ! (mag) B-V colour index (B-V) (G1)
real*4 V_I_2 (nr__2) ! (mag) V-I colour index (V-I) (G1)
real*4 Type_2 (nr__2) ! [1/4] Type of the star (Type_num) (G2)
integer*8 GaiaDR2_2 (nr__2) ! Gaia DR2 source identifier
* (GaiaDR2_source_id) (G3)
real*8 RAdeg_2 (nr__2) ! (deg) Barycentric right ascension (ICRS)
* at Ep=2015.5 (ra) (G3)
real*8 DEdeg_2 (nr__2) ! (deg) Barycentric declination (ICRS)
* at Ep=2015.5 (dec) (G3)
real*8 plx_2 (nr__2) ! (mas) ?=- Absolute stellar parallax (parallax) (G3)
real*4 e_plx_2 (nr__2) ! (mas) ?=- Mean error of plx (parallax_error) (G3)
real*8 pmRA_2 (nr__2) ! (mas/yr) ?=- Proper motion in right ascension direction
* pmRA*cosDE (pmra) (G3)
real*4 e_pmRA_2 (nr__2) ! (mas/yr) ?=- Mean error of pmRA (pmra_error) (G3)
real*8 pmDE_2 (nr__2) ! (mas/yr) ?=- Proper motion in declination
* direction (pmdec) (G3)
real*4 e_pmDE_2 (nr__2) ! (mas/yr) ?=- Mean error of pmDE (pmdec_error) (G3)
real*8 Gmag_2 (nr__2) ! (mag) G-band mean magnitude (Vega)
* (phot_g_mean_mag) (G3)
real*8 BPmag_2 (nr__2) ! (mag) ?=- Integrated BP mean magnitude (Vega)
* (phot_bp_mean_mag) (G3)
real*8 RPmag_2 (nr__2) ! (mag) ?=- Integrated RP mean magnitude (Vega)
* (phot_rp_mean_mag) (G3)
real*8 Rank_2 (nr__2) ! Position mag ranking value (rank) (G4)
real*8 Chi_2 (nr__2) ! ?=- {Chi}-value for membership to the Clouds
* ({Chi}) (G5)
character*56 ONames_2 (nr__2) ! BEL crossmatch confirmed with linked
* catalogues (other_names) (G6)
real*4 Jmag_2 (nr__2) ! (mag) ?=- Apparent magnitude in J band (Jmag) (G7)
real*4 Hmag_2 (nr__2) ! (mag) ?=- Apparent magnitude in H band (Hmag) (G7)
real*4 Ksmag_2 (nr__2) ! (mag) ?=- Apparent magnitude in Ks band (Ksmag) (G7)
real*4 v3_6mag_2 (nr__2) ! (mag) ?=- Apparent magnitude Spitzer/IRAC
* at 3.6 micron band magnitude ([3.6]) (G8)
real*4 v4_5mag_2 (nr__2) ! (mag) ?=- Apparent magnitude Spitzer/IRAC
* at 4.5 micron band magnitude ([4.5]) (G8)
real*4 v5_8mag_2 (nr__2) ! (mag) ?=- Apparent magnitude Spitzer/IRAC
* at 5.8 micron band magnitude ([5.8]) (G8)
real*4 v8_0mag_2 (nr__2) ! (mag) ?=- Apparent magnitude Spitzer/IRAC
* at 8.0 micron band magnitude ([8.0]) (G8)
C=============================================================================
C Loading file 'tablea1.dat' ! Sample of the LMC Be stars data used
* in this study
C Format for file interpretation
1 format(
+ I4,1X,A16,1X,F6.3,1X,F6.3,1X,F6.3,1X,F3.1,1X,I19,1X,F14.11,1X,
+ F15.11,1X,F7.4,1X,F6.4,1X,F7.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,
+ F9.6,1X,F9.6,1X,F9.6,1X,F20.18,1X,F21.17,1X,A13,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,2393
read(1,'(A273)')ar__
read(ar__,1)
+ ID(i__),OGLE(i__),Vmag(i__),B_V(i__),V_I(i__),Type(i__),
+ GaiaDR2(i__),RAdeg(i__),DEdeg(i__),plx(i__),e_plx(i__),
+ pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__),Gmag(i__),
+ BPmag(i__),RPmag(i__),Rank(i__),Chi(i__),ONames(i__),
+ Jmag(i__),Hmag(i__),Ksmag(i__),v3_6mag(i__),v4_5mag(i__),
+ v5_8mag(i__),v8_0mag(i__)
if (idig(ar__(99:105)).EQ.0) plx(i__) = rNULL__
if (idig(ar__(107:112)).EQ.0) e_plx(i__) = rNULL__
if (idig(ar__(114:120)).EQ.0) pmRA(i__) = rNULL__
if (idig(ar__(122:126)).EQ.0) e_pmRA(i__) = rNULL__
if (idig(ar__(128:134)).EQ.0) pmDE(i__) = rNULL__
if (idig(ar__(136:140)).EQ.0) e_pmDE(i__) = rNULL__
if (idig(ar__(152:160)).EQ.0) BPmag(i__) = rNULL__
if (idig(ar__(162:170)).EQ.0) RPmag(i__) = rNULL__
if (idig(ar__(193:213)).EQ.0) Chi(i__) = rNULL__
if (idig(ar__(229:233)).EQ.0) Jmag(i__) = rNULL__
if (idig(ar__(235:239)).EQ.0) Hmag(i__) = rNULL__
if (idig(ar__(241:245)).EQ.0) Ksmag(i__) = rNULL__
if (idig(ar__(247:252)).EQ.0) v3_6mag(i__) = rNULL__
if (idig(ar__(254:259)).EQ.0) v4_5mag(i__) = rNULL__
if (idig(ar__(261:266)).EQ.0) v5_8mag(i__) = rNULL__
if (idig(ar__(268:273)).EQ.0) v8_0mag(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),OGLE(i__),Vmag(i__),B_V(i__),V_I(i__),Type(i__),
+ GaiaDR2(i__),RAdeg(i__),DEdeg(i__),plx(i__),e_plx(i__),
+ pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__),Gmag(i__),
+ BPmag(i__),RPmag(i__),Rank(i__),Chi(i__),ONames(i__),
+ Jmag(i__),Hmag(i__),Ksmag(i__),v3_6mag(i__),v4_5mag(i__),
+ v5_8mag(i__),v8_0mag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Sample of the four exceptionally red Be stars,
* ID = 909 is in the LMC and ID = 75, 184,
* and 547 are in the SMC
C Format for file interpretation
2 format(
+ I3,1X,A18,1X,F6.3,1X,F6.3,1X,F6.3,1X,F3.1,1X,I19,1X,F14.11,1X,
+ F15.11,1X,F7.4,1X,F6.4,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,
+ F9.6,1X,F10.6,1X,F9.6,1X,F19.17,1X,F18.16,1X,A14,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,4
read(1,'(A269)')ar__1
read(ar__1,2)
+ ID_1(i__),OGLE_1(i__),Vmag_1(i__),B_V_1(i__),V_I_1(i__),
+ Type_1(i__),GaiaDR2_1(i__),RAdeg_1(i__),DEdeg_1(i__),
+ plx_1(i__),e_plx_1(i__),pmRA_1(i__),e_pmRA_1(i__),pmDE_1(i__),
+ e_pmDE_1(i__),Gmag_1(i__),BPmag_1(i__),RPmag_1(i__),
+ Rank_1(i__),Chi_1(i__),ONames_1(i__),Jmag_1(i__),Hmag_1(i__),
+ Ksmag_1(i__),v3_6mag_1(i__),v4_5mag_1(i__),v5_8mag_1(i__),
+ v8_0mag_1(i__)
c ..............Just test output...........
write(6,2)
+ ID_1(i__),OGLE_1(i__),Vmag_1(i__),B_V_1(i__),V_I_1(i__),
+ Type_1(i__),GaiaDR2_1(i__),RAdeg_1(i__),DEdeg_1(i__),
+ plx_1(i__),e_plx_1(i__),pmRA_1(i__),e_pmRA_1(i__),pmDE_1(i__),
+ e_pmDE_1(i__),Gmag_1(i__),BPmag_1(i__),RPmag_1(i__),
+ Rank_1(i__),Chi_1(i__),ONames_1(i__),Jmag_1(i__),Hmag_1(i__),
+ Ksmag_1(i__),v3_6mag_1(i__),v4_5mag_1(i__),v5_8mag_1(i__),
+ v8_0mag_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea3.dat' ! Sample of the SMC Be stars data used
* in this study
C Format for file interpretation
3 format(
+ I4,1X,A18,1X,F6.3,1X,F6.3,1X,F6.3,1X,F3.1,1X,I19,1X,F14.11,1X,
+ F15.11,1X,F7.4,1X,F6.4,1X,F7.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,
+ F9.6,1X,F9.5,1X,F9.6,1X,F19.17,1X,F21.17,1X,A56,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablea3.dat')
write(6,*) '....Loading file: tablea3.dat'
do i__=1,1004
read(1,'(A317)')ar__2
read(ar__2,3)
+ ID_2(i__),OGLE_2(i__),Vmag_2(i__),B_V_2(i__),V_I_2(i__),
+ Type_2(i__),GaiaDR2_2(i__),RAdeg_2(i__),DEdeg_2(i__),
+ plx_2(i__),e_plx_2(i__),pmRA_2(i__),e_pmRA_2(i__),pmDE_2(i__),
+ e_pmDE_2(i__),Gmag_2(i__),BPmag_2(i__),RPmag_2(i__),
+ Rank_2(i__),Chi_2(i__),ONames_2(i__),Jmag_2(i__),Hmag_2(i__),
+ Ksmag_2(i__),v3_6mag_2(i__),v4_5mag_2(i__),v5_8mag_2(i__),
+ v8_0mag_2(i__)
if (idig(ar__2(101:107)).EQ.0) plx_2(i__) = rNULL__
if (idig(ar__2(109:114)).EQ.0) e_plx_2(i__) = rNULL__
if (idig(ar__2(116:122)).EQ.0) pmRA_2(i__) = rNULL__
if (idig(ar__2(124:128)).EQ.0) e_pmRA_2(i__) = rNULL__
if (idig(ar__2(130:136)).EQ.0) pmDE_2(i__) = rNULL__
if (idig(ar__2(138:142)).EQ.0) e_pmDE_2(i__) = rNULL__
if (idig(ar__2(154:162)).EQ.0) BPmag_2(i__) = rNULL__
if (idig(ar__2(164:172)).EQ.0) RPmag_2(i__) = rNULL__
if (idig(ar__2(194:214)).EQ.0) Chi_2(i__) = rNULL__
if (idig(ar__2(273:277)).EQ.0) Jmag_2(i__) = rNULL__
if (idig(ar__2(279:283)).EQ.0) Hmag_2(i__) = rNULL__
if (idig(ar__2(285:289)).EQ.0) Ksmag_2(i__) = rNULL__
if (idig(ar__2(291:296)).EQ.0) v3_6mag_2(i__) = rNULL__
if (idig(ar__2(298:303)).EQ.0) v4_5mag_2(i__) = rNULL__
if (idig(ar__2(305:310)).EQ.0) v5_8mag_2(i__) = rNULL__
if (idig(ar__2(312:317)).EQ.0) v8_0mag_2(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ ID_2(i__),OGLE_2(i__),Vmag_2(i__),B_V_2(i__),V_I_2(i__),
+ Type_2(i__),GaiaDR2_2(i__),RAdeg_2(i__),DEdeg_2(i__),
+ plx_2(i__),e_plx_2(i__),pmRA_2(i__),e_pmRA_2(i__),pmDE_2(i__),
+ e_pmDE_2(i__),Gmag_2(i__),BPmag_2(i__),RPmag_2(i__),
+ Rank_2(i__),Chi_2(i__),ONames_2(i__),Jmag_2(i__),Hmag_2(i__),
+ Ksmag_2(i__),v3_6mag_2(i__),v4_5mag_2(i__),v5_8mag_2(i__),
+ v8_0mag_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end