FORTRAN Generation
(/./ftp/cats/J/MNRAS/411/2099)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/411/2099 into FORTRAN code for reading data files line by line.

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-20
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/411/2099    Fast-rotating M-dwarf stars in NGC 2547  (Jeffries+, 2011)
*================================================================================
*Investigating coronal saturation and supersaturation in fast-rotating
*M-dwarf stars.
*    Jeffries R.D., Jackson R.J., Briggs K.R., Evans P.A., Pye J.P.
*   <Mon. Not. R. Astron. Soc., 411, 2099-2112 (2011)>
*   =2011MNRAS.411.2099J
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table1.dat'	! X-ray detections of sources in the Irwin et al.
                           (2008, Cat. J/MNRAS/383/1588) catalogue of NGC 2547
                           members with rotation periods

      integer*4 nr__
      parameter (nr__=69)	! Number of records
      character*267 ar__   	! Full-size record

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension J2000
      real*8        DEdeg       ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*14  v_IHA2008_  ! Name (N2547-N-N-NNNN) (G1)
      integer*4     RAh         ! (h) Right ascension (J2000) (G1)
      integer*4     RAm         ! (min) Right ascension (J2000) (G1)
      real*4        RAs         ! (s) Right ascension (J2000) (G1)
      character*1   DE_         ! Declination sign (J2000) (G1)
      integer*4     DEd         ! (deg) Declination (J2000) (G1)
      integer*4     DEm         ! (arcmin) Declination (J2000) (G1)
      real*4        DEs         ! (arcsec) Declination (J2000) (G1)
      real*4        Vmag        ! (mag) Johnson V magnitude (G2)
      real*4        V_Ic        ! (mag) Johnson-Cousins V-I colour index (G2)
      real*4        Per         ! (d) Rotation period (G1)
      real*4        BC          ! (mag) V-band bolometric correction
      real*4        logLbol     ! ([Lsun]) Bolometric luminosity
      real*4        logtc       ! ([d]) Convective turnover time
      real*4        logNr       ! ([-]) Rossby number
      real*4        M           ! (Msun) Stellar mass
      real*4        R           ! (Rsun) Stellar radius
      real*4        Rkep_R      ! Keplerian corotation radius as a multiple of
*                                   the stellar radius
      real*8        RAXdeg      ! (deg) X-ray source Right ascension
*                                   in decimal degrees (J2000)
      real*8        DEXdeg      ! (deg) X-ray source Declination
*                                   in decimal degrees (J2000)
      real*8        ML          ! ?=0 Maximum likelihood of the detection
      real*4        Sep         ! (arcsec) ?=0 Separation from optical counterpart
      real*4        PN0         ! (ct/s) ?=-9. Total PN count rate (0.3-3keV)
      real*4        e_PN0       ! (ct/s) ?=-9. rms uncertainty on PN0
      real*4        PN1         ! (ct/s) ?=-9. PN count rate (0.3-1keV, S-band)
      real*4        e_PN1       ! (ct/s) ?=-9. rms uncertainty on PN1
      real*4        PN2         ! (ct/s) ?=-9. PN count rate (1-3keV, H-band)
      real*4        e_PN2       ! (ct/s) ?=-9. rms uncertainty on PN2
      real*4        HR          ! ?=-9. HPN hardness ratio [(H-S)/(H+S)]
      real*4        e_HR        ! ?=-9. rms uncertainty on HR
      real*4        M1          ! (ct/s) ?=-9. Total MOS1 count rate (0.3-3keV)
      real*4        e_M1        ! (ct/s) ?=-9. rms uncertainty on M1
      real*4        M2          ! (ct/s) ?=-9. Total MOS2 count rate (0.3-3keV)
      real*4        e_M2        ! (ct/s) ?=-9. rms uncertainty on M2
      real*4        Flux        ! (mW/m2) X-ray flux (0.3-3keV)
      real*4        e_Flux      ! (mW/m2) rms uncertainty on Flux
      real*4        logLx_L     ! X-ray-to-bolometric luminosity ratio
      real*4        e_logLx_L   ! rms uncertainty on logLx/Lbol
      character*7   NTJ         ! Naylor et al., 2002, Cat. J/MNRAS/335/291
*                                   name <Cl* NGC 2547 NTJ NN-NNNN>

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2.dat'	! The properties of stars from the Irwin et al.
                           (2008, Cat. J/MNRAS/383/1588) catalogue that have
                           known rotation periods but were not detected
                           within the XMM-Newton field

      integer*4 nr__1
      parameter (nr__1=28)	! Number of records
      character*130 ar__1  	! Full-size record

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_1     ! (deg) Right Ascension J2000
      real*8        DEdeg_1     ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*14  v_IHA2008__1 ! Name (N2547-N-N-NNNN) (G1)
      integer*4     RAh_1       ! (h) Right ascension (J2000) (G1)
      integer*4     RAm_1       ! (min) Right ascension (J2000) (G1)
      real*4        RAs_1       ! (s) Right ascension (J2000) (G1)
      character*1   DE__1       ! Declination sign (J2000) (G1)
      integer*4     DEd_1       ! (deg) Declination (J2000) (G1)
      integer*4     DEm_1       ! (arcmin) Declination (J2000) (G1)
      real*4        DEs_1       ! (arcsec) Declination (J2000) (G1)
      real*4        Vmag_1      ! (mag) Johnson V magnitude (G2)
      real*4        V_Ic_1      ! (mag) Johnson-Cousins V-I colour index (G2)
      real*4        Per_1       ! (d) Rotation period (G1)
      real*4        BC_1        ! (mag) V-band bolometric correction
      real*4        logLbol_1   ! ([Lsun]) Bolometric luminosity
      real*4        logtc_1     ! ([d]) Convective turnover time
      real*4        logNr_1     ! ([-]) Rossby number
      real*4        M_1         ! (Msun) Stellar mass
      real*4        R_1         ! (Rsun) Stellar radius
      real*4        Rkep_R_1    ! Keplerian corotation radius
*                                   (relative to the stellar radius)
      character*3   Inst        ! Instrument from which X-ray count-rate upper
*                                    limit was derived (P=PN, M1=MOS1, M2=MOS2,
*                                    M12 = average from MOS1 and MOS2)
      real*4        CR          ! (ct/s) Upper limit count-rate in instrument Inst
      real*4        Flux_1      ! (mW/m2) Upper limit X-ray flux (0.3-3keV)
      real*4        logLx_L_1   ! Upper limit X-ray-to-bolometric luminosity
*                                    ratio
      character*7   NTJ_1       ! Naylor et al., 2002, Cat. J/MNRAS/335/291
*                                    name <Cl* NGC 2547 NTJ NN-NNNN>

C=============================================================================

C  Loading file 'table1.dat'	! X-ray detections of sources in the Irwin et al.
*                           (2008, Cat. J/MNRAS/383/1588) catalogue of NGC 2547
*                           members with rotation periods

C  Format for file interpretation

    1 format(
     +  A14,1X,I1,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X,
     +  F4.2,1X,F6.3,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F5.2,1X,F9.5,1X,F8.4,1X,F7.1,1X,F4.2,1X,E8.3,1X,E8.3,1X,
     +  E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,F5.2,1X,F5.2,1X,E8.3,1X,E8.3,
     +  1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,F5.2,2X,F4.2,2X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,69
        read(1,'(A267)')ar__
        read(ar__,1)
     +  v_IHA2008_,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vmag,V_Ic,Per,BC,
     +  logLbol,logtc,logNr,M,R,Rkep_R,RAXdeg,DEXdeg,ML,Sep,PN0,e_PN0,
     +  PN1,e_PN1,PN2,e_PN2,HR,e_HR,M1,e_M1,M2,e_M2,Flux,e_Flux,
     +  logLx_L,e_logLx_L,NTJ
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  v_IHA2008_,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vmag,V_Ic,Per,BC,
     +  logLbol,logtc,logNr,M,R,Rkep_R,RAXdeg,DEXdeg,ML,Sep,PN0,e_PN0,
     +  PN1,e_PN1,PN2,e_PN2,HR,e_HR,M1,e_M1,M2,e_M2,Flux,e_Flux,
     +  logLx_L,e_logLx_L,NTJ
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table2.dat'	! The properties of stars from the Irwin et al.
*                           (2008, Cat. J/MNRAS/383/1588) catalogue that have
*                           known rotation periods but were not detected
*                           within the XMM-Newton field

C  Format for file interpretation

    2 format(
     +  A14,1X,I1,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X,
     +  F4.2,2X,F5.3,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F5.2,1X,A3,1X,E8.3,1X,E8.3,1X,F5.2,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,28
        read(1,'(A130)')ar__1
        read(ar__1,2)
     +  v_IHA2008__1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Vmag_1,
     +  V_Ic_1,Per_1,BC_1,logLbol_1,logtc_1,logNr_1,M_1,R_1,Rkep_R_1,
     +  Inst,CR,Flux_1,logLx_L_1,NTJ_1
        RAdeg_1 = rNULL__
        DEdeg_1 = rNULL__
c  Derive coordinates RAdeg_1 and DEdeg_1 from input data
c  (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
        if(RAh_1 .GT. -180) RAdeg_1=RAh_1*15.
        if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
        if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
        if(DEd_1 .GE. 0) DEdeg_1=DEd_1
        if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
        if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
        if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c    ..............Just test output...........
        write(6,2)
     +  v_IHA2008__1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Vmag_1,
     +  V_Ic_1,Per_1,BC_1,logLbol_1,logtc_1,logNr_1,M_1,R_1,Rkep_R_1,
     +  Inst,CR,Flux_1,logLx_L_1,NTJ_1
        write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end