FORTRAN Generation
(/./ftp/cats/J/ApJ/702/1507)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/702/1507 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 2026-Oct-14
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/702/1507      Mid-IR photometry in IC 1396A      (Morales-Calderon+, 2009)
*================================================================================
*Mid-infrared variability of protostars in IC 1396A.
*    Morales-Calderon M., Stauffer J.R., Rebull L., Whitney B.A.,
*    Barrado y Navascues D., Ardila D.R., Song I., Brooke T.Y., Hartmann L.,
*    Calvet N.
*   <Astrophys. J., 702, 1507-1529 (2009)>
*   =2009ApJ...702.1507M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! New candidates photometry from the co-added
                             deep mosaics

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

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

      integer*4     ID         (nr__) ! [60/74] YSO candidate name
*                                 (<[MSR2009] IC1396A-NN> in Simbad) (1)
      character*1   f_ID       (nr__) ! [d] source close to a brighter one (2)
      integer*4     RAh        (nr__) ! (h) Hour of right ascension (J2000)
      integer*4     RAm        (nr__) ! (min) Minute of right ascension (J2000)
      real*4        RAs        (nr__) ! (s) Second of right ascension (J2000)
      character*1   DE_        (nr__) ! Sign of declination (J2000)
      integer*4     DEd        (nr__) ! (deg) Degree of declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs        (nr__) ! (arcsec) Arcsecond of declination (J2000)
      real*4        Jmag       (nr__) ! (mag) ? 2MASS J-band magnitude
      character*1   l_Jmag     (nr__) ! [c] c: upper limit
      real*4        e_Jmag     (nr__) ! (mag) ? Jmag uncertainty
      real*4        Hmag       (nr__) ! (mag) ? 2MASS H-band magnitude
      character*1   l_Hmag     (nr__) ! [c] c: upper limit
      real*4        e_Hmag     (nr__) ! (mag) ? Hmag uncertainty
      real*4        Ksmag      (nr__) ! (mag) ? 2MASS Ks-band magnitude
      character*1   l_Ksmag    (nr__) ! [c] c: upper limit
      real*4        e_Ksmag    (nr__) ! (mag) ? Ksmag uncertainty
      real*4        v_3_6_     (nr__) ! (mag) Spitzer/IRAC 3.6um-band magnitude
      real*4        e__3_6_    (nr__) ! (mag) 3.6um-band magnitude uncertainty
      real*4        v_4_5_     (nr__) ! (mag) Spitzer/IRAC 4.5um-band magnitude
      real*4        e__4_5_    (nr__) ! (mag) 4.5um-band magnitude uncertainty
      real*4        v_5_8_     (nr__) ! (mag) Spitzer/IRAC 5.8um-band magnitude
      real*4        e__5_8_    (nr__) ! (mag) 5.8um-band magnitude uncertainty
      real*4        v_8_0_     (nr__) ! (mag) ? Spitzer/IRAC 8.0um-band magnitude
      real*4        e__8_0_    (nr__) ! (mag) ? 8.0um-band magnitude uncertainty
      real*4        v_24_      (nr__) ! (mag) ? Spitzer 24um-band magnitude
      real*4        e__24_     (nr__) ! (mag) ? 24um-band magnitude uncertainty
      character*2   Class      (nr__) ! Candidate class (I or II)
*Note (1): Sicilia-Aguilar et al. (2006, Cat. J/ApJ/638/897)
*     identified 57 YSOs in this region and there are another two objects
*     that were identified by Reach et al. (2005PASP..117..978R), thus we
*     name our new candidate members beginning with number 60.
*Note (2): d = This object is very faint at 24um and it is located too close
*     to a brighter object to derive a reliable magnitude.

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

C  Declarations for 'table2.dat'	! Long-term variability: main results from mapping
                             data

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

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

      integer*4     ID_1       (nr__1) ! [1/74] YSO candidate identification
*                               (<[MSR2009] IC1396A-NN> in Simbad)(1)
      character*1   f_ID_1     (nr__1) ! [f] Uncertain member (2)
      integer*4     RAh_1      (nr__1) ! (h) Hour of right ascension (J2000)
      integer*4     RAm_1      (nr__1) ! (min) Minute of right ascension (J2000)
      real*4        RAs_1      (nr__1) ! (s) Second of right ascension (J2000)
      character*1   DE__1      (nr__1) ! Sign of declination (J2000)
      integer*4     DEd_1      (nr__1) ! (deg) Degree of declination (J2000)
      integer*4     DEm_1      (nr__1) ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs_1      (nr__1) ! (arcsec) Arcsecond of declination (J2000)
      character*2   SA         (nr__1) ! YSO class as in Sicilia-Aguilar et al. (2006,
*                               Cat. J/ApJ/638/897)
      character*12  R04        (nr__1) ! ID and YSO Class as in Reach et al.
*                               (2004ApJS..154..385R)
      character*3   Class_1    (nr__1) ! YSO class in this paper (3)
      character*1   f_Class    (nr__1) ! [e] e: classification in disagreement with
*                               Sicilia-Aguilar et al. (2006, Cat. J/ApJ/638/897)
      real*4        v_3_6__1   (nr__1) ! (mag) Spitzer/IRAC 3.6um-band magnitude
      real*4        e__3_6__1  (nr__1) ! (mag) 3.6mag uncertainty
      real*4        rms1       (nr__1) ! 3.6um rms amplitude
      character*1   V1         (nr__1) ! [N/V] {chi}^2^ test variability result
      real*4        T1         (nr__1) ! (d) ? Variability timescale for the IRAC Ch. 1 (4)
      real*4        v_4_5__1   (nr__1) ! (mag) ? Spitzer/IRAC 4.5um-band magnitude
      real*4        e__4_5__1  (nr__1) ! (mag) ? 4.5mag uncertainty
      real*4        rms2       (nr__1) ! ? 4.5um rms amplitude
      character*1   V2         (nr__1) ! [N/V] {chi}^2^ test variability result
      real*4        T2         (nr__1) ! (d) ? Variability timescale for the IRAC Ch. 2 (4)
      real*4        Js         (nr__1) ! Stetson J statistic (Js>0.55 for a variable)
      character*1   Var        (nr__1) ! [N/V] Final variability tag
*Note (1): Objects are named following the order in Table 6 of
*     Sicilia-Aguilar et al. (2006, Cat. J/ApJ/638/897; note that not all
*     their members are included), objects 58 and 59 are YSOs from Reach et
*     al. (2004ApJS..154..385R) not in common with Sicilia-Aguilar et al.
*     (2006, Cat. J/ApJ/638/897), and objects from 60 to 74 are our new
*     candidate members from Table 1.
*Note (2): f = Uncertain cluster members following Sicilia-Aguilar et al.
*     (2006, Cat. J/ApJ/638/897).
*Note (3): Members in Sicilia-Aguilar et al. (2006, Cat. J/ApJ/638/897) that
*     do not show IR excess in our data have been classified as Class III.
*Note (4): Variability timescales rather than periods since the timespan of
*     our observations are not long enough to confirm the periodicity.

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

C  Declarations for 'table3.dat'	! Time series in the four IRAC bands for the
                             variable YSOs of our sample

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

      integer*4     ID_2       (nr__2) ! [1/74] Object ID
*                                 (<[MSR2009] IC1396A-NN> in Simbad)
      real*8        MJD        (nr__2) ! (d) Modified Julian Date
      real*4        v_3_6__2   (nr__2) ! (mag) ?=99.999 Spitzer/IRAC 3.6 micron band magnitude
      real*4        e__3_6__2  (nr__2) ! (mag) ?=9.999 Uncertainty in 3.6mag
      real*4        v_4_5__2   (nr__2) ! (mag) ?=99.999 Spitzer/IRAC 4.5 micron band magnitude
      real*4        e__4_5__2  (nr__2) ! (mag) ?=9.999 Uncertainty in 4.5mag
      real*4        v_5_8__1   (nr__2) ! (mag) ?=99.999 Spitzer/IRAC 5.8 micron band magnitude
      real*4        e__5_8__1  (nr__2) ! (mag) ?=9.999 Uncertainty in 5.8mag
      real*4        v_8_0__1   (nr__2) ! (mag) ?=99.999 Spitzer/IRAC 8.0 micron band magnitude
      real*4        e__8_0__1  (nr__2) ! (mag) ?=9.999 Uncertainty in 8.0mag

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

C  Declarations for 'table6.dat'	! Photometry for the variable objects found in the
                             field

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

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

      integer*4     Seq        (nr__3) ! Object internal identification
*                                 (<[MSR2009] NNNN> in Simbad)
      integer*4     RAh_2      (nr__3) ! (h) Hour of right ascension (J2000)
      integer*4     RAm_2      (nr__3) ! (min) Minute of right ascension (J2000)
      real*4        RAs_2      (nr__3) ! (s) Second of right ascension (J2000)
      character*1   DE__2      (nr__3) ! Sign of declination (J2000)
      integer*4     DEd_2      (nr__3) ! (deg) Degree of declination (J2000)
      integer*4     DEm_2      (nr__3) ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs_2      (nr__3) ! (arcsec) Arcsecond of declination (J2000)
      real*4        Jmag_1     (nr__3) ! (mag) ? 2MASS J-band magnitude
      character*1   f_Jmag     (nr__3) ! [c] c: upper limit
      real*4        e_Jmag_1   (nr__3) ! (mag) ? Jmag uncertainty
      real*4        Hmag_1     (nr__3) ! (mag) ? 2MASS H-band magnitude
      real*4        e_Hmag_1   (nr__3) ! (mag) ? Hmag uncertainty
      real*4        Ksmag_1    (nr__3) ! (mag) ? 2MASS Ks-band magnitude
      character*1   f_Ksmag    (nr__3) ! [c] c: upper limit
      real*4        e_Ksmag_1  (nr__3) ! (mag) ? Ksmag uncertainty
      real*4        v_3_6__3   (nr__3) ! (mag) Spitzer/IRAC 3.6um-band magnitude
      real*4        e__3_6__3  (nr__3) ! (mag) Spitzer/IRAC 3.6um-band magnitude uncertainty
      real*4        v_4_5__3   (nr__3) ! (mag) Spitzer/IRAC 4.5um-band magnitude
      real*4        e__4_5__3  (nr__3) ! (mag) Spitzer/IRAC 4.5um-band magnitude uncertainty
      real*4        v_5_8__2   (nr__3) ! (mag) ? Spitzer/IRAC 5.8um-band magnitude
      real*4        e__5_8__2  (nr__3) ! (mag) ? Spitzer/IRAC 5.8um-band magnitude uncertainty
      real*4        v_8_0__2   (nr__3) ! (mag) ? Spitzer/IRAC 8.0um-band magnitude
      real*4        e__8_0__2  (nr__3) ! (mag) ? Spitzer/IRAC 8.0um-band magnitude uncertainty

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

C  Loading file 'table1.dat'	! New candidates photometry from the co-added
*                             deep mosaics

C  Format for file interpretation

    1 format(
     +  8X,I2,A1,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F6.3,
     +  A1,1X,F5.3,1X,F6.3,A1,1X,F5.3,1X,F6.3,A1,1X,F5.3,1X,F6.3,1X,
     +  F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,
     +  1X,F5.3,1X,A2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,15
        read(1,'(A144)')ar__
        read(ar__,1)
     +  ID(i__),f_ID(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),Jmag(i__),l_Jmag(i__),e_Jmag(i__),
     +  Hmag(i__),l_Hmag(i__),e_Hmag(i__),Ksmag(i__),l_Ksmag(i__),
     +  e_Ksmag(i__),v_3_6_(i__),e__3_6_(i__),v_4_5_(i__),
     +  e__4_5_(i__),v_5_8_(i__),e__5_8_(i__),v_8_0_(i__),
     +  e__8_0_(i__),v_24_(i__),e__24_(i__),Class(i__)
        if(ar__(37:42) .EQ. '') Jmag(i__) = rNULL__
        if(ar__(45:49) .EQ. '') e_Jmag(i__) = rNULL__
        if(ar__(51:56) .EQ. '') Hmag(i__) = rNULL__
        if(ar__(59:63) .EQ. '') e_Hmag(i__) = rNULL__
        if(ar__(65:70) .EQ. '') Ksmag(i__) = rNULL__
        if(ar__(73:77) .EQ. '') e_Ksmag(i__) = rNULL__
        if(ar__(118:123) .EQ. '') v_8_0_(i__) = rNULL__
        if(ar__(125:129) .EQ. '') e__8_0_(i__) = rNULL__
        if(ar__(131:135) .EQ. '') v_24_(i__) = rNULL__
        if(ar__(137:141) .EQ. '') e__24_(i__) = rNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  ID(i__),f_ID(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),Jmag(i__),l_Jmag(i__),e_Jmag(i__),
     +  Hmag(i__),l_Hmag(i__),e_Hmag(i__),Ksmag(i__),l_Ksmag(i__),
     +  e_Ksmag(i__),v_3_6_(i__),e__3_6_(i__),v_4_5_(i__),
     +  e__4_5_(i__),v_5_8_(i__),e__5_8_(i__),v_8_0_(i__),
     +  e__8_0_(i__),v_24_(i__),e__24_(i__),Class(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Long-term variability: main results from mapping
*                             data

C  Format for file interpretation

    2 format(
     +  8X,I2,A1,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A2,1X,
     +  A12,1X,A3,A1,1X,F6.3,1X,F5.3,1X,F5.3,1X,A1,1X,F4.1,1X,F6.3,1X,
     +  F5.3,1X,F5.3,1X,A1,1X,F3.1,1X,E6.2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,69
        read(1,'(A116)')ar__1
        read(ar__1,2)
     +  ID_1(i__),f_ID_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
     +  DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),SA(i__),R04(i__),
     +  Class_1(i__),f_Class(i__),v_3_6__1(i__),e__3_6__1(i__),
     +  rms1(i__),V1(i__),T1(i__),v_4_5__1(i__),e__4_5__1(i__),
     +  rms2(i__),V2(i__),T2(i__),Js(i__),Var(i__)
        if(ar__1(79:82) .EQ. '') T1(i__) = rNULL__
        if(ar__1(84:89) .EQ. '') v_4_5__1(i__) = rNULL__
        if(ar__1(91:95) .EQ. '') e__4_5__1(i__) = rNULL__
        if(ar__1(97:101) .EQ. '') rms2(i__) = rNULL__
        if(ar__1(105:107) .EQ. '') T2(i__) = rNULL__
        RAdeg_1(i__) = rNULL__
        DEdeg_1(i__) = 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(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
        if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
        if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
        if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
        if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
        if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
        if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  ID_1(i__),f_ID_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
     +  DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),SA(i__),R04(i__),
     +  Class_1(i__),f_Class(i__),v_3_6__1(i__),e__3_6__1(i__),
     +  rms1(i__),V1(i__),T1(i__),v_4_5__1(i__),e__4_5__1(i__),
     +  rms2(i__),V2(i__),T2(i__),Js(i__),Var(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Time series in the four IRAC bands for the
*                             variable YSOs of our sample

C  Format for file interpretation

    3 format(
     +  8X,I2,1X,F13.7,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
     +  F5.3,1X,F6.3,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,1148
        read(1,'(A76)')ar__2
        read(ar__2,3)
     +  ID_2(i__),MJD(i__),v_3_6__2(i__),e__3_6__2(i__),v_4_5__2(i__),
     +  e__4_5__2(i__),v_5_8__1(i__),e__5_8__1(i__),v_8_0__1(i__),
     +  e__8_0__1(i__)
c    ..............Just test output...........
        write(6,3)
     +  ID_2(i__),MJD(i__),v_3_6__2(i__),e__3_6__2(i__),v_4_5__2(i__),
     +  e__4_5__2(i__),v_5_8__1(i__),e__5_8__1(i__),v_8_0__1(i__),
     +  e__8_0__1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Photometry for the variable objects found in the
*                             field

C  Format for file interpretation

    4 format(
     +  I4,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F6.3,A1,1X,
     +  F5.3,1X,F6.3,1X,F5.3,1X,F6.3,A1,1X,F5.3,1X,F6.3,1X,F5.3,1X,
     +  F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,46
        read(1,'(A121)')ar__3
        read(ar__3,4)
     +  Seq(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),
     +  DEd_2(i__),DEm_2(i__),DEs_2(i__),Jmag_1(i__),f_Jmag(i__),
     +  e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),Ksmag_1(i__),
     +  f_Ksmag(i__),e_Ksmag_1(i__),v_3_6__3(i__),e__3_6__3(i__),
     +  v_4_5__3(i__),e__4_5__3(i__),v_5_8__2(i__),e__5_8__2(i__),
     +  v_8_0__2(i__),e__8_0__2(i__)
        if(ar__3(30:35) .EQ. '') Jmag_1(i__) = rNULL__
        if(ar__3(38:42) .EQ. '') e_Jmag_1(i__) = rNULL__
        if(ar__3(44:49) .EQ. '') Hmag_1(i__) = rNULL__
        if(ar__3(51:55) .EQ. '') e_Hmag_1(i__) = rNULL__
        if(ar__3(57:62) .EQ. '') Ksmag_1(i__) = rNULL__
        if(ar__3(65:69) .EQ. '') e_Ksmag_1(i__) = rNULL__
        if(ar__3(97:102) .EQ. '') v_5_8__2(i__) = rNULL__
        if(ar__3(104:108) .EQ. '') e__5_8__2(i__) = rNULL__
        if(ar__3(110:115) .EQ. '') v_8_0__2(i__) = rNULL__
        if(ar__3(117:121) .EQ. '') e__8_0__2(i__) = rNULL__
        RAdeg_2(i__) = rNULL__
        DEdeg_2(i__) = rNULL__
c  Derive coordinates RAdeg_2 and DEdeg_2 from input data
c  (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
        if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
        if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
        if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
        if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
        if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
        if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
        if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c    ..............Just test output...........
        write(6,4)
     +  Seq(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),
     +  DEd_2(i__),DEm_2(i__),DEs_2(i__),Jmag_1(i__),f_Jmag(i__),
     +  e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),Ksmag_1(i__),
     +  f_Ksmag(i__),e_Ksmag_1(i__),v_3_6__3(i__),e__3_6__3(i__),
     +  v_4_5__3(i__),e__4_5__3(i__),v_5_8__2(i__),e__5_8__2(i__),
     +  v_8_0__2(i__),e__8_0__2(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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