Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/414/2602 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-19
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/414/2602 Automated classification of HIP variables (Dubath+, 2011)
*================================================================================
*Random forest automated supervised classification of Hipparcos periodic variable
*stars.
* Dubath P., Rimoldini L., Suveges M., Blomme J., Lopez M., Sarro L.M.,
* De Ridder J., Cuypers J., Guy L., Lecoeur I., Nienartowicz K., Jan A.,
* Beck M., Mowlavi N., De Cat P., Lebzelter T., Eyer L.
* <Mon. Not. R. Astron. Soc., 414, 2602-2617 (2011)>
* =2011MNRAS.414.2602D
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! The Hipparcos training set star list with
literature types and attribute values
integer*4 nr__
parameter (nr__=1661) ! Number of records
character*110 ar__ ! Full-size record
integer*4 HIP (nr__) ! HIP number
character*8 Type (nr__) ! Type of variable
character*7 PType (nr__) ! Predicted type (only in table4)
real*8 logP (nr__) ! ([d]) Period extracted with the Lomb-Scargle method
* [Log(Period)]
real*4 logA (nr__) ! ([---]) Amplitude of the light-curve model
* [Log(Amplitude)]
real*4 V_I (nr__) ! (mag) Mean V-I colour index [V-I]
real*4 Mhip (nr__) ! (mag) Hipparcos absolute magnitude derived from the
* parallaxes neglecting interstellar
* absorption [M_Hipparcos_] (1)
real*4 res_raw (nr__) ! Median absolute of the residuals (obtained by
* subtracting model values from the raw light
* curve) divided by the Median Absolute
* Deviation (MAD) of the raw light-curve values
* around the median [Scatter:res/raw]
real*4 Skew (nr__) ! Unbiased weighted skewness of the magnitude
* distribution [Skewness]
real*4 log_1_A2_A1(nr__) ! ([-]) Amplitude ratio between the second harmonic
* and the fundamental (plus one, to avoid
* negative values) [Log(1+A2/A1)]
real*4 P2p_2P (nr__) ! Sum of the squares of the magnitude
* differences between pairs of successive data
* points in the light curve folded around twice
* the period divided by the same quantity
* derived from the raw light curve
* [P2p scatter:2P/raw]
real*4 P2p (nr__) ! Median of the absolute values of the
* differences between successive magnitudes in
* the raw light curve normalized by the MAD
* around the median [P2p_scatter]
real*4 P90 (nr__) ! The 90th percentile of the absolute residual
* values around the 2P model divided by the
* same quantity for the residuals around the P
* model. The 2P model is a model recomputed
* using twice the period value
* [Percentile90:2P/P]
real*4 Res (nr__) ! (%) Mean of the squared residuals around the
* model [Residual_scatter]
real*4 Phase2 (nr__) ! (rad) Phase of the second harmonic after setting the
* phase of the fundamental to zero [Phase2] (2)
real*4 P2P_P (nr__) ! Median of the absolute values of the
* differences between successive magnitudes in
* the folded light curve normalized by the MAD
* around the median of the raw light curve
* [P2p_scatter:P/raw]
real*4 Slope (nr__) ! Sum of the square of the slopes of lines
* joining the data points before and after a
* number of selected outliers towards faint
* magnitude (e.g., data points during eclipses)
* [P2p_slope] (3)
*Note (1): Because of measurement uncertainties, some stars have negative
* parallax values. Each of these values is replaced by a positive value taken
* randomly from a Gaussian distribution with zero mean and a standard deviation
* equal to the measurement uncertainty. In many cases, the derived absolute
* magnitudes represent lower limits as the parallax measurements are not
* significant.
*Note (2): by an appropriate transformation
* Phase2=arctan(sin({phi}_2_-2{phi}_1_),cos({phi}_2_-2{phi}_1_))
* (Debosscher et al., 2007, Cat. J/A+A/475/1159).
*Note (3): This is set to zero if there are no such outliers in the light curve.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Results obtained for the Hipparcos stars
excluded from the training set
integer*4 nr__1
parameter (nr__1=882) ! Number of records
character*110 ar__1 ! Full-size record
integer*4 HIP_1 (nr__1) ! HIP number
character*8 Type_1 (nr__1) ! Type of variable
character*7 PType_1 (nr__1) ! Predicted type (only in table4)
real*8 logP_1 (nr__1) ! ([d]) Period extracted with the Lomb-Scargle method
* [Log(Period)]
real*4 logA_1 (nr__1) ! ([---]) Amplitude of the light-curve model
* [Log(Amplitude)]
real*4 V_I_1 (nr__1) ! (mag) Mean V-I colour index [V-I]
real*4 Mhip_1 (nr__1) ! (mag) Hipparcos absolute magnitude derived from the
* parallaxes neglecting interstellar
* absorption [M_Hipparcos_] (1)
real*4 res_raw_1 (nr__1) ! Median absolute of the residuals (obtained by
* subtracting model values from the raw light
* curve) divided by the Median Absolute
* Deviation (MAD) of the raw light-curve values
* around the median [Scatter:res/raw]
real*4 Skew_1 (nr__1) ! Unbiased weighted skewness of the magnitude
* distribution [Skewness]
real*4 log_1_A2_A1_1(nr__1) ! ([-]) Amplitude ratio between the second harmonic
* and the fundamental (plus one, to avoid
* negative values) [Log(1+A2/A1)]
real*4 P2p_2P_1 (nr__1) ! Sum of the squares of the magnitude
* differences between pairs of successive data
* points in the light curve folded around twice
* the period divided by the same quantity
* derived from the raw light curve
* [P2p scatter:2P/raw]
real*4 P2p_1 (nr__1) ! Median of the absolute values of the
* differences between successive magnitudes in
* the raw light curve normalized by the MAD
* around the median [P2p_scatter]
real*4 P90_1 (nr__1) ! The 90th percentile of the absolute residual
* values around the 2P model divided by the
* same quantity for the residuals around the P
* model. The 2P model is a model recomputed
* using twice the period value
* [Percentile90:2P/P]
real*4 Res_1 (nr__1) ! (%) Mean of the squared residuals around the
* model [Residual_scatter]
real*4 Phase2_1 (nr__1) ! (rad) Phase of the second harmonic after setting the
* phase of the fundamental to zero [Phase2] (2)
real*4 P2P_P_1 (nr__1) ! Median of the absolute values of the
* differences between successive magnitudes in
* the folded light curve normalized by the MAD
* around the median of the raw light curve
* [P2p_scatter:P/raw]
real*4 Slope_1 (nr__1) ! Sum of the square of the slopes of lines
* joining the data points before and after a
* number of selected outliers towards faint
* magnitude (e.g., data points during eclipses)
* [P2p_slope] (3)
*Note (1): Because of measurement uncertainties, some stars have negative
* parallax values. Each of these values is replaced by a positive value taken
* randomly from a Gaussian distribution with zero mean and a standard deviation
* equal to the measurement uncertainty. In many cases, the derived absolute
* magnitudes represent lower limits as the parallax measurements are not
* significant.
*Note (2): by an appropriate transformation
* Phase2=arctan(sin({phi}_2_-2{phi}_1_),cos({phi}_2_-2{phi}_1_))
* (Debosscher et al., 2007, Cat. J/A+A/475/1159).
*Note (3): This is set to zero if there are no such outliers in the light curve.
C=============================================================================
C Loading file 'table3.dat' ! The Hipparcos training set star list with
* literature types and attribute values
C Format for file interpretation
1 format(
+ I6,1X,A8,1X,A7,1X,F7.4,1X,F5.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,1661
read(1,'(A110)')ar__
read(ar__,1)
+ HIP(i__),Type(i__),PType(i__),logP(i__),logA(i__),V_I(i__),
+ Mhip(i__),res_raw(i__),Skew(i__),log_1_A2_A1(i__),P2p_2P(i__),
+ P2p(i__),P90(i__),Res(i__),Phase2(i__),P2P_P(i__),Slope(i__)
c ..............Just test output...........
write(6,1)
+ HIP(i__),Type(i__),PType(i__),logP(i__),logA(i__),V_I(i__),
+ Mhip(i__),res_raw(i__),Skew(i__),log_1_A2_A1(i__),P2p_2P(i__),
+ P2p(i__),P90(i__),Res(i__),Phase2(i__),P2P_P(i__),Slope(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Results obtained for the Hipparcos stars
* excluded from the training set
C Format for file interpretation
2 format(
+ I6,1X,A8,1X,A7,1X,F7.4,1X,F5.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,882
read(1,'(A110)')ar__1
read(ar__1,2)
+ HIP_1(i__),Type_1(i__),PType_1(i__),logP_1(i__),logA_1(i__),
+ V_I_1(i__),Mhip_1(i__),res_raw_1(i__),Skew_1(i__),
+ log_1_A2_A1_1(i__),P2p_2P_1(i__),P2p_1(i__),P90_1(i__),
+ Res_1(i__),Phase2_1(i__),P2P_P_1(i__),Slope_1(i__)
c ..............Just test output...........
write(6,2)
+ HIP_1(i__),Type_1(i__),PType_1(i__),logP_1(i__),logA_1(i__),
+ V_I_1(i__),Mhip_1(i__),res_raw_1(i__),Skew_1(i__),
+ log_1_A2_A1_1(i__),P2p_2P_1(i__),P2p_1(i__),P90_1(i__),
+ Res_1(i__),Phase2_1(i__),P2P_P_1(i__),Slope_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end