FORTRAN Generation
(/./ftp/cats/J/ApJ/731/120)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/731/120 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-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/ApJ/731/120     Intrinsic SN Ia light curves               (Mandel+, 2011)
*================================================================================
*Type Ia supernova light curve inference: hierarchical models in the optical and
*near-infrared.
*    Mandel K.S., Narayan G., Kirshner R.P.
*   <Astrophys. J., 731, 120 (2011)>
*   =2011ApJ...731..120M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table4.dat'	! Apparent light curve and dust estimates for
                             individual SN Ia

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

      character*6   SN          ! SN name (YYYYaa)
      real*4        B0mag       ! (mag) ?=99.99 Rest-frame B band maximum magnitude (1)
      real*4        e_B0mag     ! (mag) ?=9.99 Uncertainty in B0mag
      real*4        Dm15        ! (mag) ?=9.99 B band decline rate {Delta}m_15_(B)
*                                        from peak to 15 day past maximum.
      real*4        e_Dm15      ! (mag) ?=9.99 Uncertainty in Dm15
      real*4        V0mag       ! (mag) ?=99.99 Rest-frame V band magnitude (2)
      real*4        e_V0mag     ! (mag) ?=9.99 Uncertainty in V0mag
      real*4        R0mag       ! (mag) ?=99.99 Rest-frame R band magnitude (2)
      real*4        e_R0mag     ! (mag) ?=9.99 Uncertainty in R0mag
      real*4        I0mag       ! (mag) ?=99.99 Rest-frame I band magnitude (2)
      real*4        e_I0mag     ! (mag) ?=9.99 Uncertainty in I0mag
      real*4        J0mag       ! (mag) ?=99.99 Rest-frame J band magnitude (2)
      real*4        e_J0mag     ! (mag) ?=9.99 Uncertainty in J0mag
      real*4        H0mag       ! (mag) ?=99.99 Rest-frame H band magnitude (2)
      real*4        e_H0mag     ! (mag) Uncertainty in H0mag
      real*4        AV          ! (mag) Marginal posterior mode of V band extinction
      real*4        e_AV        ! (mag) Lower 68% interval in AV
      real*4        E_AV_1      ! (mag) Upper 68% interval in AV
      real*4        Rv          ! Marginal posterior mean ratio of
*                                 total-to-selective dust absorption
*                                 (R_V_=A_V_/(A_B_-A_V_))
      real*4        e_Rv        ! Standard Deviation uncertainty in Rv
      character*7   Ref         ! Reference code (3)
*Note (1): Corrected for Milky Way extinction and K-corrections.
*     Estimates only listed if SN was observed in the filter.
*Note (2): At time of maximum in B, corrected for Milky Way extinction and
*     K-corrections.
*Note (3): Reference as follows:
*    CfA3 = Hicken et al. (2009, Cat. J/ApJ/700/331);
*    WV08 = Wood-Vasey et al. (PAIRITEL; 2008, Cat. J/ApJ/689/377);
*     WC3 = WV08+CfA3;
*     J99 = Jha et al. (1999ApJS..125...73J);
*     H00 = Hernandez et al. (2000MNRAS.319..223H);
*     K00 = Krisciunas et al. (2000ApJ...539..658K);
*     K01 = Krisciunas et al. (2001AJ....122.1616K);
*    DP02 = Di Paola et al. (2002A&A...393L..21D);
*     V03 = Valentini et al. (2003ApJ...595..779V);
*     K03 = Krisciunas et al. (2003AJ....125..166K);
*    K04b = Krisciunas et al. (2004AJ....127.1664K);
*    K04c = Krisciunas et al. (2004AJ....128.3034K);
*     K07 = Krisciunas et al. (2007AJ....133...58K);
*    ER06 = Elias-Rosa et al. (2006, Cat. J/MNRAS/369/1880);
*    ER07 = Elias-Rosa et al. (2008, Cat. J/MNRAS/384/107);
*    Pa07 = Pastorello et al. (2007, Cat. J/MNRAS/376/1301);
*    St07 = Stanishev et al. (2007A&A...469..645S);
*     P08 = Pignata et al. (2008MNRAS.388..971P).

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

C  Declarations for 'table5.dat'	! Distance modulus predictions for SN Ia

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

      character*6   SN_1        ! SN name (YYYYaa)
      real*8        cz          ! (km/s) Redshift
      real*4        muz         ! (mag) Distance modulus expected from redshift (4)
      real*4        e_muz       ! (mag) Distance modulus standard deviation (5)
      real*4        mu1         ! (mag) Distance modulus estimated from resubstitution
      real*4        mu2         ! (mag) Mean predicted distance modulus over bootstrap
      real*4        e_mu2       ! (mag) Standard deviation of mu2 (6)
      real*4        sig2        ! (mag) Average standard deviation of uncertainty in mu2
      real*4        AVm         ! (mag) Marginal posterior mode V band extinction (7)
*Note (4): Assuming h=0.72, {Omega}_M_=0.27, {Omega}_{Lambda}_=0.73, w=-1.
*Note (5): Assuming peculiar velocity dispersion {sigma}_pec_=150km/s.
*Note (6): Zero values are less than 0.005mag.
*Note (7): Averaged over the prediction sets for each SN.

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

C  Loading file 'table4.dat'	! Apparent light curve and dust estimates for
*                             individual SN Ia

C  Format for file interpretation

    1 format(
     +  2X,A6,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,
     +  1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,
     +  F4.2,1X,F4.2,1X,F4.2,1X,F3.1,1X,F3.1,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,110
        read(1,'(A115)')ar__
        read(ar__,1)
     +  SN,B0mag,e_B0mag,Dm15,e_Dm15,V0mag,e_V0mag,R0mag,e_R0mag,
     +  I0mag,e_I0mag,J0mag,e_J0mag,H0mag,e_H0mag,AV,e_AV,E_AV_1,Rv,
     +  e_Rv,Ref
c    ..............Just test output...........
        write(6,1)
     +  SN,B0mag,e_B0mag,Dm15,e_Dm15,V0mag,e_V0mag,R0mag,e_R0mag,
     +  I0mag,e_I0mag,J0mag,e_J0mag,H0mag,e_H0mag,AV,e_AV,E_AV_1,Rv,
     +  e_Rv,Ref
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Distance modulus predictions for SN Ia

C  Format for file interpretation

    2 format(
     +  2X,A6,1X,F8.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,110
        read(1,'(A55)')ar__1
        read(ar__1,2)SN_1,cz,muz,e_muz,mu1,mu2,e_mu2,sig2,AVm
c    ..............Just test output...........
        write(6,2)SN_1,cz,muz,e_muz,mu1,mu2,e_mu2,sig2,AVm
c    .......End.of.Just test output...........
      end do
      close(1)

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