FORTRAN Generation
(/./ftp/cats/J/A_A/506/1277)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/506/1277 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-Mar-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/A+A/506/1277   Modelling of Magellanic Cloud C/O stars (Groenewegen+, 2009)
*================================================================================
*Luminosities and mass-loss rates of SMC and LMC AGB stars and red supergiants.
*    Groenewegen M.A.T., Sloan G.C., Soszynski I., Petersen E.A.
*   <Astron. Astrophys. 506, 1277 (2009)>
*   =2009A&A...506.1277G
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1c.dat'	! C-star sample: identifiers and pulsation periods

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

C  J2000 position composed of: RAdeg DEdeg
      character*62  Names      (nr__) ! Various identifiers
      real*8        RAdeg      (nr__) ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg      (nr__) ! (deg) Declination, in decimal degrees (J2000)
      character*11  Ident      (nr__) ! Identifier in Tables and Plots
      character*19  Per        (nr__) ! (d) Pulsation Period(s)
      character*14  Ref        (nr__) ! Reference code
      real*4        Amp1       (nr__) ! (mag) ? Pulsation Amplitude
      character*1   n_Amp1     (nr__) ! [RIKVB] Filter for Amp1
      real*4        Amp2       (nr__) ! (mag) ? Pulsation Amplitude
      character*1   n_Amp2     (nr__) ! [RK] Filter for Amp2

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

C  Declarations for 'table2o.dat'	! O-star sample: identifiers and pulsation periods

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

C  J2000 position composed of: RAdeg DEdeg
      character*62  Names_1    (nr__1) ! Various identifiers
      real*8        RAdeg_1    (nr__1) ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Declination, in decimal degrees (J2000)
      character*11  Ident_1    (nr__1) ! Identifier in Tables and Plots
      character*19  Per_1      (nr__1) ! (d) Pulsation Period(s)
      character*14  Ref_1      (nr__1) ! Reference code
      real*4        Amp1_1     (nr__1) ! (mag) ? Pulsation Amplitude
      character*1   n_Amp1_1   (nr__1) ! [RIKVB] Filter for Amp1
      real*4        Amp2_1     (nr__1) ! (mag) ? Pulsation Amplitude
      character*1   n_Amp2_1   (nr__1) ! [RK] Filter for Amp2

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

C  Declarations for 'table3c.dat'	! Fit results for the C-star sample

      integer*4 nr__2
      parameter (nr__2=101)	! Number of records
      character*91 ar__2  	! Full-size record

      character*11  Ident_2    (nr__2) ! Identifier
      integer*4     Teff       (nr__2) ! (K) Effective temperature
      character*5   Type       (nr__2) ! Identification of type of dust (1)
      integer*4     Lum        (nr__2) ! (solLum) Luminosity in solar units
      real*4        dM_dt      (nr__2) ! (solMass/yr) Mass-loss rate
      integer*4     n_dM_dt    (nr__2) ! [0/1] Code (1=free param; 0=fixed value)
      integer*4     Tc         (nr__2) ! (K) Temperature at inner radius
      integer*4     n_Tc       (nr__2) ! [0/1] Code (1=free param; 0=fixed value)
      real*8        Router     (nr__2) ! Outer radius (relative to inner radius)
      real*4        tau1       (nr__2) ! Optical depth at 0.55 micron
      real*4        tau2       (nr__2) ! Optical depth at 11.3 micron
*Note (1): combination of amorphous carbon (AMC) and silicon carbide (SiC)
*     ({alpha}-SiC and {beta}-SiC). The number refers to the percentage of SiC.

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

C  Declarations for 'table4m.dat'	! Fit results for the O-star sample

      integer*4 nr__3
      parameter (nr__3=172)	! Number of records
      character*109 ar__3  	! Full-size record

      character*11  Ident_3    (nr__3) ! Identifier
      character*3   SpType     (nr__3) ! MK spectral type
      character*21  Type_1     (nr__3) ! Identification of type of dust (1)
      integer*4     Lum_1      (nr__3) ! (solLum) Luminosity in solar units
      real*4        dM_dt_1    (nr__3) ! (solMass/yr) Mass-loss rate
      integer*4     n_dM_dt_1  (nr__3) ! [0/1] Code (1=free param; 0=fixed value)
      integer*4     Tc_1       (nr__3) ! (K) Temperature at inner radius
      integer*4     n_Tc_1     (nr__3) ! [0/1] Code (1=free param; 0=fixed value)
      real*8        Router_1   (nr__3) ! Outer radius (relative to inner radius)
      real*4        tau1_1     (nr__3) ! Optical depth at 0.55 micron
      real*4        tau2_1     (nr__3) ! Optical depth at 10 micron
*Note (1): For every star there are two lines, for "laboratory silicates" and
*     "astronomical silicates". The first line provides the best fitting one.

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

C  Loading file 'table1c.dat'	! C-star sample: identifiers and pulsation periods

C  Format for file interpretation

    1 format(
     +  A62,1X,F10.6,1X,F10.6,2X,A11,4X,A19,3X,A14,2X,F5.3,1X,A1,3X,
     +  F4.2,2X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1c.dat')
      write(6,*) '....Loading file: table1c.dat'
      do i__=1,101
        read(1,'(A157)')ar__
        read(ar__,1)
     +  Names(i__),RAdeg(i__),DEdeg(i__),Ident(i__),Per(i__),Ref(i__),
     +  Amp1(i__),n_Amp1(i__),Amp2(i__),n_Amp2(i__)
        if(ar__(140:144) .EQ. '') Amp1(i__) = rNULL__
        if(ar__(150:153) .EQ. '') Amp2(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Names(i__),RAdeg(i__),DEdeg(i__),Ident(i__),Per(i__),Ref(i__),
     +  Amp1(i__),n_Amp1(i__),Amp2(i__),n_Amp2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2o.dat'	! O-star sample: identifiers and pulsation periods

C  Format for file interpretation

    2 format(
     +  A62,1X,F10.6,1X,F10.6,2X,A11,4X,A19,3X,A14,2X,F5.3,1X,A1,3X,
     +  F4.2,2X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2o.dat')
      write(6,*) '....Loading file: table2o.dat'
      do i__=1,86
        read(1,'(A157)')ar__1
        read(ar__1,2)
     +  Names_1(i__),RAdeg_1(i__),DEdeg_1(i__),Ident_1(i__),
     +  Per_1(i__),Ref_1(i__),Amp1_1(i__),n_Amp1_1(i__),Amp2_1(i__),
     +  n_Amp2_1(i__)
        if(ar__1(140:144) .EQ. '') Amp1_1(i__) = rNULL__
        if(ar__1(150:153) .EQ. '') Amp2_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Names_1(i__),RAdeg_1(i__),DEdeg_1(i__),Ident_1(i__),
     +  Per_1(i__),Ref_1(i__),Amp1_1(i__),n_Amp1_1(i__),Amp2_1(i__),
     +  n_Amp2_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3c.dat'	! Fit results for the C-star sample

C  Format for file interpretation

    3 format(
     +  A11,3X,I4,4X,A5,5X,I5,3X,E8.2,3X,I1,3X,I4,3X,I1,3X,F7.1,3X,
     +  F6.3,4X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3c.dat')
      write(6,*) '....Loading file: table3c.dat'
      do i__=1,101
        read(1,'(A91)')ar__2
        read(ar__2,3)
     +  Ident_2(i__),Teff(i__),Type(i__),Lum(i__),dM_dt(i__),
     +  n_dM_dt(i__),Tc(i__),n_Tc(i__),Router(i__),tau1(i__),tau2(i__)
c    ..............Just test output...........
        write(6,3)
     +  Ident_2(i__),Teff(i__),Type(i__),Lum(i__),dM_dt(i__),
     +  n_dM_dt(i__),Tc(i__),n_Tc(i__),Router(i__),tau1(i__),tau2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4m.dat'	! Fit results for the O-star sample

C  Format for file interpretation

    4 format(
     +  A11,3X,A3,7X,A21,3X,I7,3X,E8.2,3X,I1,3X,I4,3X,I1,3X,F7.1,3X,
     +  F6.3,4X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4m.dat')
      write(6,*) '....Loading file: table4m.dat'
      do i__=1,172
        read(1,'(A109)')ar__3
        read(ar__3,4)
     +  Ident_3(i__),SpType(i__),Type_1(i__),Lum_1(i__),dM_dt_1(i__),
     +  n_dM_dt_1(i__),Tc_1(i__),n_Tc_1(i__),Router_1(i__),
     +  tau1_1(i__),tau2_1(i__)
c    ..............Just test output...........
        write(6,4)
     +  Ident_3(i__),SpType(i__),Type_1(i__),Lum_1(i__),dM_dt_1(i__),
     +  n_dM_dt_1(i__),Tc_1(i__),n_Tc_1(i__),Router_1(i__),
     +  tau1_1(i__),tau2_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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