FORTRAN Generation
(/./ftp/cats/J/A_A/634/A37)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/634/A37 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 2026-Sep-12
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/634/A37       TURTLS Light curves of ^56^Ni distributions   (Magee+, 2020)
*================================================================================
*Determining the ^56^Ni distribution of type Ia supernovae from observations
*within days of explosion.
*    Magee M.R., Maguire K., Kotak R., Sim S.A., Gillanders J.H., Prentice S.J.,
*    Skillen K.
*    <Astron. Astrophys. 634, A37 (2020)>
*    =2020A&A...634A..37M        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'list.dat'	! List of files

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

      real*4        NiMass      ! (Msun) ^56^Ni mass  (0.4, 0.6 or 0.8)
      real*4        KE          ! Kinetic energy (0.50, 0.60, 0.65, 0.78, 1.10,
*                                   1.40, 1.68, 1.81 or 2.18)
      real*4        P           ! Scaling parameter (3, 4.4, 9.7, 21 or 100)
      character*36  dp          ! Name of the file with model density profiles
*                                   in subdirectory dp
      character*28  lc          ! Name of the file with model light curves
*                                   in subdirectory lc

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

C  Declarations for 'lc/*'	! Individual model light curves

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

      real*4        Time        ! (d) Days since explosion
      real*4        UMAG        ! (mag) ? Absolute U band magnitude
      real*4        BMAG        ! (mag) ? Absolute B band magnitude
      real*4        VMAG        ! (mag) Absolute V band magnitude
      real*4        RMAG        ! (mag) Absolute R band magnitude
      real*4        IMAG        ! (mag) Absolute I band magnitude
      real*4        gMAG        ! (mag) ? Absolute g band magnitude
      character*1   n_gMAG      ! [i] i for infinity
      real*4        rMAG_1      ! (mag) Absolute r band magnitude
      real*4        iMAG_1      ! (mag) Absolute i band magnitude

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

C  Declarations for 'dp/*'	! Individual model density profiles

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

      real*4        Index       ! Zone index
      real*4        Velocity    ! (km/s) Zone outer velocity
      real*4        Density     ! (g/cm3) Zone density at 20s after explosion
      real*4        Fraction    ! ^56^Ni mass fraction at explosion

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

C  Loading file 'list.dat'	! List of files

C  Format for file interpretation

    1 format(2X,F3.1,3X,F4.2,2X,F5.1,1X,A36,4X,A28)

C  Effective file loading

      open(unit=1,status='old',file=
     +'list.dat')
      write(6,*) '....Loading file: list.dat'
      do i__=1,255
        read(1,'(A88)')ar__
        read(ar__,1)NiMass,KE,P,dp,lc
c    ..............Just test output...........
        write(6,1)NiMass,KE,P,dp,lc
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'lc/*'	! Individual model light curves

C  Format for file interpretation

    2 format(
     +  E12.7,4X,E13.7,3X,E13.7,3X,E13.7,3X,E13.7,3X,E13.7,3X,E13.7,
     +  A1,2X,E13.7,3X,E13.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'lc/*')
      write(6,*) '....Loading file: lc/*'
      do i__=1,255
        read(1,'(A141)')ar__1
        read(ar__1,2)
     +  Time,UMAG,BMAG,VMAG,RMAG,IMAG,gMAG,n_gMAG,rMAG_1,iMAG_1
        if(ar__1(17:29) .EQ. '') UMAG = rNULL__
        if(ar__1(33:45) .EQ. '') BMAG = rNULL__
        if(ar__1(97:109) .EQ. '') gMAG = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Time,UMAG,BMAG,VMAG,RMAG,IMAG,gMAG,n_gMAG,rMAG_1,iMAG_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'dp/*'	! Individual model density profiles

C  Format for file interpretation

    3 format(E12.7,2X,E12.7,7X,E12.7,3X,E12.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'dp/*')
      write(6,*) '....Loading file: dp/*'
      do i__=1,255
        read(1,'(A60)')ar__2
        read(ar__2,3)Index,Velocity,Density,Fraction
c    ..............Just test output...........
        write(6,3)Index,Velocity,Density,Fraction
c    .......End.of.Just test output...........
      end do
      close(1)

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