FORTRAN Generation
(/./ftp/cats/J/A_A/630/A75)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/630/A75 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-Apr-23
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/630/A75       Six luminous red novae photometry        (Pastorello+, 2019)
*================================================================================
*Luminous red novae: Stellar mergers or giant eruptions?
*    Pastorello A., Mason E., Taubenberger S., Fraser M., Cortini G.,
*    Tomasella L., Botticella M.T., Elias-Rosa N., Kotak R., Smartt S.J.,
*    Benetti S., Cappellaro E., Turatto M., Tartaglia L., Djorgovski S.G.,
*    Drake A.J., Berton M., Briganti F., Brimacombe J., Bufano F., Cai Y.-Z.,
*    Chen S., Christensen E.J., Ciabattari F., Congiu E., Dimai A., Inserra C.,
*    Kankare E., Magill L., Maguire K., Martinelli F., Morales-Garoffolo A.,
*    Ochner P., Pignata G., Reguitti A., Sollerman J., Spiro S., Terreran G.,
*    Wright D.E.
*    <Astron. Astrophys. 630, A75 (2019)>
*    =2019A&A...630A..75P        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'objects.dat'	! List of studied novae

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

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

      character*21  Name        ! Name
      integer*4     RAh         ! (h) Simbad Hour of Right Ascension (J2000.0)
      integer*4     RAm         ! (min) Simbad Minute of Right Ascension (J2000.0)
      real*4        RAs         ! (s) Simbad Second of Right Ascension (J2000.0)
      character*1   DE_         ! Simbad Sign of the Declination (J2000.0)
      integer*4     DEd         ! (deg) Simbad Degree of Declination (J2000.0)
      integer*4     DEm         ! (arcmin) Simbad Arcminute of Declination (J2000.0)
      real*4        DEs         ! (arcsec) Simbad Arcsecond of Declination (J2000.0)
      character*7   FileName    ! Name of the table with photometry
      character*16  OName       ! Other name

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

C  Declarations for 'tablea1.dat'	! Johnson-Cousins UBVRI (Vega mag) photometry of
                              NGC4490-2011OT1 (PSN J12304185+4137498)

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

      character*8   Date        ! ("date") Date of observation (dd/mm/yy)
      real*8        JD          ! (d) Julian Date (JD-2400000)
      character*1   l_Umag      ! Limit flag on Umag
      real*4        Umag        ! (mag) ?=- Johnson-Cousins U band magnitude (Vega)
      real*4        e_Umag      ! (mag) ? Johnson-Cousins U band magnitude error
      character*1   l_Bmag      ! Limit flag on Bmag
      real*4        Bmag        ! (mag) ?=- Johnson-Cousins B band magnitude (Vega)
      real*4        e_Bmag      ! (mag) ? Johnson-Cousins B band magnitude error
      character*1   l_Vmag      ! Limit flag on Vmag
      real*4        Vmag        ! (mag) ?=- Johnson-Cousins V band magnitude (Vega)
      real*4        e_Vmag      ! (mag) ? Johnson-Cousins V band magnitude error
      character*1   l_Rmag      ! Limit flag on Rmag
      real*4        Rmag        ! (mag) ?=- Johnson-Cousins R band magnitude (Vega)
      real*4        e_Rmag      ! (mag) ? Johnson-Cousins R band magnitude error
      character*1   l_Imag      ! Limit flag on Imag
      real*4        Imag        ! (mag) ?=- Johnson-Cousins I band magnitude (Vega)
      real*4        e_Imag      ! (mag) ? Johnson-Cousins I band magnitude error
      integer*4     Tel         ! [1/20] Telescope numerical code (1)
      character*1   Note        ! [abc] Notes (2)
*Note (1): Telescope numerical code as follows:
*  1 = 68/76cm CSS Schmidt Telescope + SI 600-277 CCD
*  2 = 1.2m S. Oschin Telescope (PTF) + CHF12K CCD Mosaic
*  3 = 2.0m Faulkes Telescope North + EM01
*  4 = 0.356m C14 Telescope (Monte Maggiore Observatory) + Starlight MX 916
*  5 = 0.53m Telescope (Montarrenti Observatory) + Apogee Alta U47 CCD
*  6 = 2.54m Isaac Newton Telescope + WFC (chip 4)
*  7 = 0.356m C14 Telescope (Lajatico Astronomical Centre) + kaf0260 CCD
*  8 = 0.51m Lotti Telescope (Monte Agliale Observatory) + FLI 4710 Proline CCD
*  9 = 0.30m Maioni Telescope (Cortina) + SBIG ST-8 CCD
* 10 = 1.82m Ekar Copernico Telescope + AFOSC
* 11 = 4.2m William Herschel Telescope + ACAM
* 12 = Ekar 67/92cm Schmidt Telecope + Sbig STL-11000M + KAI-11000M CCD
* 13 = 3.58m Telescopio Nazionale Galileo + Dolores
* 14 = 2.56m Nordic Optical Telescope + ALFOSC
* 15 = Unfiltered observation taken from the TOCP page
*       (http://www.cbat.eps.harvard.edu/unconf/followups/J12304185+4137492.html)
* 16 = 2.0m Liverpool Telescope + RATCam
* 17 = 2.2m Calar Alto Telescope + CAFOS
* 18 = 2x8.4m Large Binocular Telescope + MODS
* 19 = 10.4m Gran Telescopio Canarias + OSIRIS
* 20 = 67/92cm Schmidt Telescope + Moravian G4-16000LC + KAF-16803 CCD
*Note (2): Notes as follows:
*   a = Unfiltered data, calibrated to Johnson-Cousins R
*   b = Sloan g, r, i filters, magnitudes converted to Johnson-Cousins B, R and I
*   c = Unfiltered, template-subtracted CRTS photometry. The template was
*        obtained by stacking good-quality CRTS images collected during 2014 and
*        2015. he unfiltered data were then reported to Johnson-Cousins V
*        adopting a colour correction as in Drake et al. (2013ApJ...763...32D)

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

C  Declarations for 'tablea2.dat'	! *Johnson-Cousins UBVRI (Vega mag) photometry of
                              NGC3437-2011OT1 (PSN J10523453+2256052)

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

      character*17  Date_1      ! Date of observation
*                                   (dd/mm/yy or Stack MMMDD-MMMDD)
      real*8        JD_1        ! (d) Julian Date (JD-2400000)
      character*1   l_Umag_1    ! Limit flag on Umag
      real*4        Umag_1      ! (mag) ?=- Johnson-Cousins U band magnitude
      real*4        e_Umag_1    ! (mag) ? Johnson-Cousins U band magnitude error
      character*1   l_Bmag_1    ! Limit flag on Bmag
      real*4        Bmag_1      ! (mag) ?=- Johnson-Cousins B band magnitude
      real*4        e_Bmag_1    ! (mag) ? Johnson-Cousins B band magnitude error
      character*1   l_Vmag_1    ! Limit flag on Vmag
      real*4        Vmag_1      ! (mag) ?=- Johnson-Cousins V band magnitude
      real*4        e_Vmag_1    ! (mag) ? Johnson-Cousins V band magnitude error
      character*1   l_Rmag_1    ! Limit flag on Rmag
      real*4        Rmag_1      ! (mag) ?=- Johnson-Cousins R band magnitude
      real*4        e_Rmag_1    ! (mag) ? Johnson-Cousins R band magnitude error
      character*1   l_Imag_1    ! Limit flag on Imag
      real*4        Imag_1      ! (mag) ?=- Johnson-Cousins I band magnitude
      real*4        e_Imag_1    ! (mag) ? Johnson-Cousins I band magnitude error
      integer*4     Tel_1       ! [1/15] Telescope numerical code (1)
      character*3   Note_1      ! Notes (2)
*Note (1): Telescope numerical code as follows:
*  1 = 68/76cm CSS Schmidt Telescope + SI 600-277 CCD
*  2 = 2.54m Isaac Newton Telescope  + WFC (chip 4)
*  3 = 1.2m S. Oschin Telescope (PTF) + CHF12K CCD Mosaic
*  4 = 0.51m Lotti Telescope (Monte Agliale Observatory) + FLI 4710 Proline CCD
*  5 = CBET 2637
*  6 = 20-inch Telescope (New Mexico Skies)
*  7 = Liverpool Telescope + RATCam
*  8 = 2.2m Calar Alto Telescope + CAFOS
*  9 = 0.35m Telescope (Lajatico Astronomical Centre) + DTA KAF 0260 CCD
* 10 = 3.58m New Technology Telescope + EFOSC2
* 11 = 0.203m Meade Telescope (Clair Observatory) + ST10-XME CCD
* 12 = 0.356m C14 Telescope (Monte Maggiore Observatory) + Starlight MX 916
* 13 = 1.82m Ekar Copernico Telescope + AFOSC
* 14 = Ekar~67/92cm Schmidt Telescope + SBIG CCD
* 15 = Ekar~67/92cm Schmidt Telescope + Moravian G4-16000LC + KAF-16803 CCD
*Note (2): Notes as follows:
*  a = CRTS/CSS Y filter data, converted to Johnson-Cousins R band
*  b = Measures obtained with images after the template subtraction
*  c = CRTS/CSS Johnson-Cousins R-band data

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

C  Declarations for 'tablea3.dat'	! *Sloan g (AB mag) and Johnson-Cousins BVRI
                              (Vega mag) photometry of UGC12307-2013OT1
                              (PSN J23011153+1243218)

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

      character*8   Date_2      ! ("date") Date of observation (dd/mm/yy)
      real*8        JD_2        ! (d) Julian Date (JD-2400000)
      character*1   l_gmag      ! Limit flag on gmag
      real*4        gmag        ! (mag) ?=- Sloan g band magnitude (AB)
      character*1   l_Bmag_2    ! Limit flag on Bmag
      real*4        Bmag_2      ! (mag) ?=- Johnson-Cousins B band magnitude (Vega)
      real*4        e_Bmag_2    ! (mag) ? Johnson-Cousins B band magnitude error
      character*1   l_Vmag_2    ! Limit flag on Vmag
      real*4        Vmag_2      ! (mag) ?=- Johnson-Cousins V band magnitude (Vega)
      real*4        e_Vmag_2    ! (mag) ? Johnson-Cousins V band magnitude error
      character*1   l_Rmag_2    ! Limit flag on Rmag
      real*4        Rmag_2      ! (mag) ?=- Johnson-Cousins R band magnitude (Vega)
      real*4        e_Rmag_2    ! (mag) ? Johnson-Cousins R band magnitude error
      character*1   l_Imag_2    ! Limit flag on Imag
      real*4        Imag_2      ! (mag) ?=- Johnson-Cousins I band magnitude (Vega)
      real*4        e_Imag_2    ! (mag) ? Johnson-Cousins I band magnitude error
      integer*4     Tel_2       ! [1/11] Telescope numerical code (1)
      character*1   Note_2      ! Notes (2)
*Note (1): Telescope numerical code as follows:
*   1 = 68/76cm CSS Schmidt Telescope + SI 600-277 CCD
*   2 = 1.2m S. Oschin Telescope (PTF) + CHF12K CCD Mosaic
*   3 = 0.51m Lotti Telescope (Monte Agliale Observatory) + SBIG ST-9E CCD
*   4 = 0.51m Lotti Telescope (Monte Agliale Observatory) + FLI 4710 Proline CCD
*   5 = 0.41m Coral Towers Observatory + SBIG STL 6K CCD
*   6 = 1.82m Ekar Copernico Telescope + AFOSC
*   7 = 3.58m Telescopio Nazionale Galileo + Dolores
*   8 = 2m Liverpool Telescope + RATCam
*   9 = 10.4m Gran Telescopio Canarias + OSIRIS
*  10 = 4.2m William Herschel Telescope + ACAM
*  11 = Ekar 67/92cm Schmidt Telescope + Moravian G4-16000LC + KAF-16803 CCD
*Note (2): Notes as follows:
*  a = Sloan g, r and i observations, converted to Johnson-Cousins B, V, R, I
*        magnitudes using the relations from Jester et al. (2005AJ....130..873J,
*        Cat. J/AJ/130/873)

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

C  Declarations for 'tablea4.dat'	! Johnson-Cousins BVR (Vega mag) and Sloan ugriz
                              (AB mag) photometry of M101-2015OT1
                              (PSN J14021678+5426205)

      integer*4 nr__4
      parameter (nr__4=75)	! Number of records
      character*116 ar__4  	! Full-size record

      character*8   Date_3      ! ("date") Date of observation (dd/mm/yy)
      real*8        JD_3        ! (d) Julian Date (JD-2400000)
      character*1   l_Bmag_3    ! Limit flag on Bmag
      real*4        Bmag_3      ! (mag) ?=- Johnson-Cousins B band magnitude (Vega)
      real*4        e_Bmag_3    ! (mag) ? Johnson-Cousins B band magnitude error
      character*1   l_Vmag_3    ! Limit flag on Vmag
      real*4        Vmag_3      ! (mag) ?=- Johnson-Cousins V band magnitude (Vega) (2)
      real*4        e_Vmag_3    ! (mag) ? Johnson-Cousins V band magnitude error
      character*1   l_Rmag_3    ! Limit flag on Rmag
      real*4        Rmag_3      ! (mag) ?=- Johnson-Cousins R band magnitude (Vega) (2)
      real*4        e_Rmag_3    ! (mag) ? Johnson-Cousins R band magnitude error
      character*1   l_umag_2    ! Limit flag on umag
      real*4        umag_2      ! (mag) ?=- Sloan u band magnitude (AB)
      real*4        e_umag_2    ! (mag) ? Sloan u band magnitude error
      character*1   l_gmag_1    ! Limit flag on gmag
      real*4        gmag_1      ! (mag) ?=- Sloan g band magnitude (AB)
      real*4        e_gmag      ! (mag) ? Sloan g band magnitude error
      character*1   l_rmag_4    ! Limit flag on rmag
      real*4        rmag_4      ! (mag) ?=- Sloan r band magnitude (AB)
      real*4        e_rmag_4    ! (mag) ? Sloan r band magnitude error
      character*1   l_imag_3    ! Limit flag on imag
      real*4        imag_3      ! (mag) ?=- Sloan i band magnitude (AB)
      real*4        e_imag_3    ! (mag) ? Sloan i band magnitude error
      character*1   l_zmag      ! Limit flag on zmag
      real*4        zmag        ! (mag) ?=- Sloan z band magnitude (AB)
      real*4        e_zmag      ! (mag) ? Sloan z band magnitude error
      integer*4     Tel_3       ! [1/13] Telescope numerical code (1)
*Note (1): Telescope numerical code as follows:
*   1 = 0.51m Lotti Telescope (Monte Agliale Observatory) + FLI 4710 Proline CCD
*   2 = 0.356m C14 Telescope (Lajatico Astronomical Centre) + DTA KAF 0261 CCD
*   3 = 0.356m C14 Telescope (Lajatico Astronomical Centre) + DTA KAF 1001E CCD
*   4 = 0.356m C14 Telescope (Monte Maggiore Observatory) + Starlight MX 916
*   5 = 0.50m Helmut Ullrich Telescope (Col Driuscie' Observatory)
*        + SBIG ST-8 3 CCD
*   6 = 0.30m Vittorio Maioni Telescope (Col Driusce' Observatory)
*        + SBIG ST-8 3 CCD
*   7 = 0.254 Meade Telescope (Lajatico Astronomical Centre) + DTA KAF 1600 CCD
*   8 = 0.302 Meade Telescope (Lajatico Astronomical Centre)
*        + SBIG ST-7 KAF 402 CCD
*   9 = 0.30m Telescope at the Itagaki Astronomical Observatory + KAF-1001E CCD
*  10 = 1.82m Ekar Copernico Telescope + AFOSC
*  11 = 2.56m Nordic Optical Telescope + ALFOSC
*  12 = 10.4m Gran telescopio Canarias + OSIRIS
*  13 = Ekar 67/92cm Schmidt Telescope + Moravian G4-16000LC + KAF-16803 CCD
*Note (2): Unfiltered amateur observations are reported to Johnson-Cousins V
* or R band, depending on the wavelengths of the CCD quantum efficiency peaks.

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

C  Declarations for 'tablea5.dat'	! Johnson-Cousins UBVRI (Vega mag) photometry of
                              SN 1997bs

      integer*4 nr__5
      parameter (nr__5=3)	! Number of records
      character*75 ar__5  	! Full-size record

      character*8   Date_4      ! ("date") Date of observation (dd/mm/yy)
      real*8        JD_4        ! (d) Julian Date (JD-2400000)
      real*4        Umag_3      ! (mag) ?=- Johnson-Cousins U band magnitude (Vega)
      real*4        e_Umag_3    ! (mag) ? Johnson-Cousins U band magnitude error
      real*4        Bmag_4      ! (mag) ?=- Johnson-Cousins B band magnitude (Vega)
      real*4        e_Bmag_4    ! (mag) ? Johnson-Cousins B band magnitude error
      real*4        Vmag_4      ! (mag) ?=- Johnson-Cousins V band magnitude (Vega)
      real*4        e_Vmag_4    ! (mag) ? Johnson-Cousins V band magnitude error
      character*1   l_Rmag_5    ! Limit flag on Rmag
      real*4        Rmag_5      ! (mag) Johnson-Cousins R band magnitude (Vega)
      real*4        e_Rmag_5    ! (mag) ? Johnson-Cousins R band magnitude error
      real*4        Imag_4      ! (mag) ?=- Johnson-Cousins I band magnitude (Vega)
      real*4        e_Imag_4    ! (mag) ? Johnson-Cousins I band magnitude error
      integer*4     Tel_4       ! Telescope numerical code (1)
*Note (1): Telescope numerical code as follows:
*    1 = Dutch 0.92m Telescope (ESO-La Silla)

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

C  Declarations for 'tablea6.dat'	! Johnson-Cousins VR (Vega mag) and Sloan gr (AB mag)
                              photometry of NGC5806-2014OT1 (SNhunt248)

      integer*4 nr__6
      parameter (nr__6=261)	! Number of records
      character*57 ar__6  	! Full-size record

      character*8   Date_5      ! ("date") Date of observation (dd/mm/yy)
      real*8        JD_5        ! (d) Julian Date (JD-2400000)
      character*1   l_Vmag_4    ! Limit flag on Vmag
      real*4        Vmag_5      ! (mag) ?=- Johnson-Cousins V band magnitude (Vega)
      character*1   l_Rmag_6    ! Limit flag on Rmag
      real*4        Rmag_6      ! (mag) ?=- Johnson-Cousins R band magnitude (Vega)
      real*4        e_Rmag_6    ! (mag) ? Johnson-Cousins R band magnitude error
      character*1   l_gmag_2    ! Limit flag on gmag (AB)
      real*4        gmag_2      ! (mag) ?=- Sloan g band magnitude
      real*4        e_gmag_1    ! (mag) ? Sloan g band magnitude error
      character*1   l_rmag_7    ! Limit flag on rmag
      real*4        rmag_7      ! (mag) ?=- Sloan r band magnitude (AB)
      integer*4     Tel_5       ! Telescope numerical code (1)
*Note (1): Telescope numerical code as follows:
*  1 = 1.2m S.Oschin Telescope (PTF) + CHF12K CCD Mosaic
*  2 = 1.82m Ekar Copernico Telescope + AFOSC

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

C  Loading file 'objects.dat'	! List of studied novae

C  Format for file interpretation

    1 format(
     +  A21,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A7,1X,A16)

C  Effective file loading

      open(unit=1,status='old',file=
     +'objects.dat')
      write(6,*) '....Loading file: objects.dat'
      do i__=1,6
        read(1,'(A70)')ar__
        read(ar__,1)Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,FileName,OName
        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)Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,FileName,OName
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea1.dat'	! Johnson-Cousins UBVRI (Vega mag) photometry of
*                              NGC4490-2011OT1 (PSN J12304185+4137498)

C  Format for file interpretation

    2 format(
     +  A8,1X,F8.2,1X,A1,F6.3,1X,F5.3,1X,A1,F6.3,1X,F5.3,1X,A1,F6.3,
     +  1X,F5.3,1X,A1,F6.3,1X,F5.3,1X,A1,F6.3,1X,F5.3,1X,I2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea1.dat')
      write(6,*) '....Loading file: tablea1.dat'
      do i__=1,346
        read(1,'(A92)')ar__1
        read(ar__1,2)
     +  Date,JD,l_Umag,Umag,e_Umag,l_Bmag,Bmag,e_Bmag,l_Vmag,Vmag,
     +  e_Vmag,l_Rmag,Rmag,e_Rmag,l_Imag,Imag,e_Imag,Tel,Note
        if (idig(ar__1(20:25)).EQ.0) Umag =  rNULL__
        if(ar__1(27:31) .EQ. '') e_Umag = rNULL__
        if (idig(ar__1(34:39)).EQ.0) Bmag =  rNULL__
        if(ar__1(41:45) .EQ. '') e_Bmag = rNULL__
        if (idig(ar__1(48:53)).EQ.0) Vmag =  rNULL__
        if(ar__1(55:59) .EQ. '') e_Vmag = rNULL__
        if (idig(ar__1(62:67)).EQ.0) Rmag =  rNULL__
        if(ar__1(69:73) .EQ. '') e_Rmag = rNULL__
        if (idig(ar__1(76:81)).EQ.0) Imag =  rNULL__
        if(ar__1(83:87) .EQ. '') e_Imag = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Date,JD,l_Umag,Umag,e_Umag,l_Bmag,Bmag,e_Bmag,l_Vmag,Vmag,
     +  e_Vmag,l_Rmag,Rmag,e_Rmag,l_Imag,Imag,e_Imag,Tel,Note
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea2.dat'	! *Johnson-Cousins UBVRI (Vega mag) photometry of
*                              NGC3437-2011OT1 (PSN J10523453+2256052)

C  Format for file interpretation

    3 format(
     +  A17,1X,F8.2,1X,A1,F6.3,1X,F5.3,1X,A1,F6.3,1X,F5.3,1X,A1,F6.3,
     +  1X,F5.3,1X,A1,F6.3,1X,F5.3,1X,A1,F6.3,1X,F5.3,1X,I2,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea2.dat')
      write(6,*) '....Loading file: tablea2.dat'
      do i__=1,218
        read(1,'(A103)')ar__2
        read(ar__2,3)
     +  Date_1,JD_1,l_Umag_1,Umag_1,e_Umag_1,l_Bmag_1,Bmag_1,e_Bmag_1,
     +  l_Vmag_1,Vmag_1,e_Vmag_1,l_Rmag_1,Rmag_1,e_Rmag_1,l_Imag_1,
     +  Imag_1,e_Imag_1,Tel_1,Note_1
        if (idig(ar__2(29:34)).EQ.0) Umag_1 =  rNULL__
        if(ar__2(36:40) .EQ. '') e_Umag_1 = rNULL__
        if (idig(ar__2(43:48)).EQ.0) Bmag_1 =  rNULL__
        if(ar__2(50:54) .EQ. '') e_Bmag_1 = rNULL__
        if (idig(ar__2(57:62)).EQ.0) Vmag_1 =  rNULL__
        if(ar__2(64:68) .EQ. '') e_Vmag_1 = rNULL__
        if (idig(ar__2(71:76)).EQ.0) Rmag_1 =  rNULL__
        if(ar__2(78:82) .EQ. '') e_Rmag_1 = rNULL__
        if (idig(ar__2(85:90)).EQ.0) Imag_1 =  rNULL__
        if(ar__2(92:96) .EQ. '') e_Imag_1 = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Date_1,JD_1,l_Umag_1,Umag_1,e_Umag_1,l_Bmag_1,Bmag_1,e_Bmag_1,
     +  l_Vmag_1,Vmag_1,e_Vmag_1,l_Rmag_1,Rmag_1,e_Rmag_1,l_Imag_1,
     +  Imag_1,e_Imag_1,Tel_1,Note_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea3.dat'	! *Sloan g (AB mag) and Johnson-Cousins BVRI
*                              (Vega mag) photometry of UGC12307-2013OT1
*                              (PSN J23011153+1243218)

C  Format for file interpretation

    4 format(
     +  A8,1X,F8.2,1X,A1,F5.2,1X,A1,F6.3,1X,F5.3,1X,A1,F6.3,1X,F5.3,
     +  1X,A1,F6.3,1X,F5.3,1X,A1,F6.3,1X,F5.3,1X,I2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea3.dat')
      write(6,*) '....Loading file: tablea3.dat'
      do i__=1,293
        read(1,'(A85)')ar__3
        read(ar__3,4)
     +  Date_2,JD_2,l_gmag,gmag,l_Bmag_2,Bmag_2,e_Bmag_2,l_Vmag_2,
     +  Vmag_2,e_Vmag_2,l_Rmag_2,Rmag_2,e_Rmag_2,l_Imag_2,Imag_2,
     +  e_Imag_2,Tel_2,Note_2
        if (idig(ar__3(20:24)).EQ.0) gmag =  rNULL__
        if (idig(ar__3(27:32)).EQ.0) Bmag_2 =  rNULL__
        if(ar__3(34:38) .EQ. '') e_Bmag_2 = rNULL__
        if (idig(ar__3(41:46)).EQ.0) Vmag_2 =  rNULL__
        if(ar__3(48:52) .EQ. '') e_Vmag_2 = rNULL__
        if (idig(ar__3(55:60)).EQ.0) Rmag_2 =  rNULL__
        if(ar__3(62:66) .EQ. '') e_Rmag_2 = rNULL__
        if (idig(ar__3(69:74)).EQ.0) Imag_2 =  rNULL__
        if(ar__3(76:80) .EQ. '') e_Imag_2 = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Date_2,JD_2,l_gmag,gmag,l_Bmag_2,Bmag_2,e_Bmag_2,l_Vmag_2,
     +  Vmag_2,e_Vmag_2,l_Rmag_2,Rmag_2,e_Rmag_2,l_Imag_2,Imag_2,
     +  e_Imag_2,Tel_2,Note_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea4.dat'	! Johnson-Cousins BVR (Vega mag) and Sloan ugriz
*                              (AB mag) photometry of M101-2015OT1
*                              (PSN J14021678+5426205)

C  Format for file interpretation

    5 format(
     +  A8,1X,F8.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,
     +  1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
     +  F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea4.dat')
      write(6,*) '....Loading file: tablea4.dat'
      do i__=1,75
        read(1,'(A116)')ar__4
        read(ar__4,5)
     +  Date_3,JD_3,l_Bmag_3,Bmag_3,e_Bmag_3,l_Vmag_3,Vmag_3,e_Vmag_3,
     +  l_Rmag_3,Rmag_3,e_Rmag_3,l_umag_2,umag_2,e_umag_2,l_gmag_1,
     +  gmag_1,e_gmag,l_rmag_4,rmag_4,e_rmag_4,l_imag_3,imag_3,
     +  e_imag_3,l_zmag,zmag,e_zmag,Tel_3
        if (idig(ar__4(20:24)).EQ.0) Bmag_3 =  rNULL__
        if(ar__4(26:29) .EQ. '') e_Bmag_3 = rNULL__
        if (idig(ar__4(32:36)).EQ.0) Vmag_3 =  rNULL__
        if(ar__4(38:41) .EQ. '') e_Vmag_3 = rNULL__
        if (idig(ar__4(44:48)).EQ.0) Rmag_3 =  rNULL__
        if(ar__4(50:53) .EQ. '') e_Rmag_3 = rNULL__
        if (idig(ar__4(56:60)).EQ.0) umag_2 =  rNULL__
        if(ar__4(62:65) .EQ. '') e_umag_2 = rNULL__
        if (idig(ar__4(68:72)).EQ.0) gmag_1 =  rNULL__
        if(ar__4(74:77) .EQ. '') e_gmag = rNULL__
        if (idig(ar__4(80:84)).EQ.0) rmag_4 =  rNULL__
        if(ar__4(86:89) .EQ. '') e_rmag_4 = rNULL__
        if (idig(ar__4(92:96)).EQ.0) imag_3 =  rNULL__
        if(ar__4(98:101) .EQ. '') e_imag_3 = rNULL__
        if (idig(ar__4(104:108)).EQ.0) zmag =  rNULL__
        if(ar__4(110:113) .EQ. '') e_zmag = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Date_3,JD_3,l_Bmag_3,Bmag_3,e_Bmag_3,l_Vmag_3,Vmag_3,e_Vmag_3,
     +  l_Rmag_3,Rmag_3,e_Rmag_3,l_umag_2,umag_2,e_umag_2,l_gmag_1,
     +  gmag_1,e_gmag,l_rmag_4,rmag_4,e_rmag_4,l_imag_3,imag_3,
     +  e_imag_3,l_zmag,zmag,e_zmag,Tel_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea5.dat'	! Johnson-Cousins UBVRI (Vega mag) photometry of
*                              SN 1997bs

C  Format for file interpretation

    6 format(
     +  A8,1X,F8.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,
     +  A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea5.dat')
      write(6,*) '....Loading file: tablea5.dat'
      do i__=1,3
        read(1,'(A75)')ar__5
        read(ar__5,6)
     +  Date_4,JD_4,Umag_3,e_Umag_3,Bmag_4,e_Bmag_4,Vmag_4,e_Vmag_4,
     +  l_Rmag_5,Rmag_5,e_Rmag_5,Imag_4,e_Imag_4,Tel_4
        if (idig(ar__5(19:23)).EQ.0) Umag_3 =  rNULL__
        if(ar__5(25:28) .EQ. '') e_Umag_3 = rNULL__
        if (idig(ar__5(30:34)).EQ.0) Bmag_4 =  rNULL__
        if(ar__5(36:39) .EQ. '') e_Bmag_4 = rNULL__
        if (idig(ar__5(41:45)).EQ.0) Vmag_4 =  rNULL__
        if(ar__5(47:50) .EQ. '') e_Vmag_4 = rNULL__
        if(ar__5(59:62) .EQ. '') e_Rmag_5 = rNULL__
        if (idig(ar__5(64:68)).EQ.0) Imag_4 =  rNULL__
        if(ar__5(70:73) .EQ. '') e_Imag_4 = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  Date_4,JD_4,Umag_3,e_Umag_3,Bmag_4,e_Bmag_4,Vmag_4,e_Vmag_4,
     +  l_Rmag_5,Rmag_5,e_Rmag_5,Imag_4,e_Imag_4,Tel_4
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea6.dat'	! Johnson-Cousins VR (Vega mag) and Sloan gr (AB mag)
*                              photometry of NGC5806-2014OT1 (SNhunt248)

C  Format for file interpretation

    7 format(
     +  A8,1X,F8.2,1X,A1,F5.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,
     +  1X,A1,F5.2,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea6.dat')
      write(6,*) '....Loading file: tablea6.dat'
      do i__=1,261
        read(1,'(A57)')ar__6
        read(ar__6,7)
     +  Date_5,JD_5,l_Vmag_4,Vmag_5,l_Rmag_6,Rmag_6,e_Rmag_6,l_gmag_2,
     +  gmag_2,e_gmag_1,l_rmag_7,rmag_7,Tel_5
        if (idig(ar__6(20:24)).EQ.0) Vmag_5 =  rNULL__
        if (idig(ar__6(27:31)).EQ.0) Rmag_6 =  rNULL__
        if(ar__6(33:36) .EQ. '') e_Rmag_6 = rNULL__
        if (idig(ar__6(39:43)).EQ.0) gmag_2 =  rNULL__
        if(ar__6(45:48) .EQ. '') e_gmag_1 = rNULL__
        if (idig(ar__6(51:55)).EQ.0) rmag_7 =  rNULL__
c    ..............Just test output...........
        write(6,7)
     +  Date_5,JD_5,l_Vmag_4,Vmag_5,l_Rmag_6,Rmag_6,e_Rmag_6,l_gmag_2,
     +  gmag_2,e_gmag_1,l_rmag_7,rmag_7,Tel_5
c    .......End.of.Just test output...........
      end do
      close(1)

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

C Locate position of first digit in string; or return 0
      integer function idig(c)
      character*(*) c
      character*1 c1
      integer lc,i
      lc=len(c)
      idig=0
      do i=1,lc
         if(c(i:i).ne.' ') go to 1
      end do
    1 if(i.gt.lc) return
      c1=c(i:i)
      if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
      if(i.gt.lc) return
      c1=c(i:i)
      if(c1.ge.'0'.and.c1.le.'9') idig=i
      return
      end