FORTRAN Generation
(/./ftp/cats/J/ApJ/678/200)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/678/200 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/678/200             YSOs in N63 and N180                   (Caulet+, 2008)
*================================================================================
*Young stellar objects in the Large Magellanic Cloud: N63 and N180 HII regions.
*    Caulet A., Gruendl R.A., Chu Y.-H.
*   <Astrophys. J., 678, 200-218 (2008)>
*   =2008ApJ...678..200C
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! Galaxies and AGB stars

      integer*4 nr__
      parameter (nr__=18)	! Number of records
      character*141 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)

      character*9   Field      (nr__) ! Field of view
      character*1   f_Field    (nr__) ! [cd] Flag on Field (1)
      character*3   cl         (nr__) ! AGB or Galaxy type (2)
      character*1   f_cl       (nr__) ! [e] superimposed sources (3)
      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__) ! Declination sign (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
      real*4        e_Jmag     (nr__) ! (mag) ? Jmag uncertainty
      real*4        Hmag       (nr__) ! (mag) ? 2MASS H-band magnitude
      real*4        e_Hmag     (nr__) ! (mag) ? Hmag uncertainty
      real*4        Ksmag      (nr__) ! (mag) ? 2MASS Ks-band magnitude
      real*4        e_Ksmag    (nr__) ! (mag) ? Ksmag uncertainty
      real*4        v3_6mag    (nr__) ! (mag) Spitzer/IRAC 3.6um magnitude (G1)
      real*4        e_3_6mag   (nr__) ! (mag) 3.6mag uncertainty (G1)
      real*4        v4_5mag    (nr__) ! (mag) Spitzer/IRAC 4.5um magnitude (G1)
      real*4        e_4_5mag   (nr__) ! (mag) 4.5mag uncertainty (G1)
      real*4        v5_8mag    (nr__) ! (mag) Spitzer/IRAC 5.8um magnitude (G1)
      real*4        e_5_8mag   (nr__) ! (mag) 5.8mag uncertainty (G1)
      real*4        v8_0mag    (nr__) ! (mag) Spitzer/IRAC 8.0um magnitude (G1)
      real*4        e_8_0mag   (nr__) ! (mag) 8.0mag uncertainty (G1)
      real*4        v24mag     (nr__) ! (mag) ? Spitzer/MIPS 24.0um magnitude (G1)
      real*4        e_24mag    (nr__) ! (mag) ? 24mag uncertainty (G1)
*Note (1): Note on Field of view (FOV) as follows:
*    c = FOV of 40'x40' shown in Figs. 4 and 8 for N63 and N180, respectively.
*        Two AGB stars are located within the boundaries of Fig. 9.
*    d = FOV of 10'x10' of ISPI images for N63 and N180; the galaxies are shown
*        in Figs. 5, 6, and 9.
*Note (2): These objects were removed from our lists of classified YSOs.
*     They are evolved stars (AGB) identified in the CMD diagrams or
*     background galaxies (G) identified on ISPI images (see text).
*Note (3):
*    e = There are two 2MASS superimposed sources at this position; the
*        magnitudes of the other source are [J]=13.858+/-0.026,
*        [H]=12.864+/-0.034, and [Ks]=11.841+/-0.271 .

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

C  Declarations for 'table5.dat'	! Properties of the YSOs of N63

      integer*4 nr__1
      parameter (nr__1=29)	! Number of records
      character*141 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)

      character*3   Name       (nr__1) ! YSO name (see Nomenclature notes) (1)
      character*2   cl_1       (nr__1) ! [1/2/fs ] YSO class (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) ! Declination sign (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*1   f_Jmag     (nr__1) ! [*] This paper ISPI Jmag; 2MASS otherwise
      real*4        Jmag_1     (nr__1) ! (mag) ? 2MASS J-band magnitude
      real*4        e_Jmag_1   (nr__1) ! (mag) ? Jmag uncertainty
      real*4        Hmag_1     (nr__1) ! (mag) ? 2MASS H-band magnitude
      real*4        e_Hmag_1   (nr__1) ! (mag) ? Hmag uncertainty
      character*1   f_Ksmag    (nr__1) ! [*] This paper ISPI Ksmag; 2MASS otherwise
      real*4        Ksmag_1    (nr__1) ! (mag) ? 2MASS Ks-band magnitude
      real*4        e_Ksmag_1  (nr__1) ! (mag) ? Ksmag uncertainty
      real*4        v3_6mag_1  (nr__1) ! (mag) Spitzer/IRAC 3.6um magnitude (G1)
      real*4        e_3_6mag_1 (nr__1) ! (mag) 3.6mag uncertainty (G1)
      real*4        v4_5mag_1  (nr__1) ! (mag) Spitzer/IRAC 4.5um magnitude (G1)
      real*4        e_4_5mag_1 (nr__1) ! (mag) 4.5mag uncertainty (G1)
      real*4        v5_8mag_1  (nr__1) ! (mag) Spitzer/IRAC 5.8um magnitude (G1)
      real*4        e_5_8mag_1 (nr__1) ! (mag) 5.8mag uncertainty (G1)
      real*4        v8_0mag_1  (nr__1) ! (mag) Spitzer/IRAC 8.0um magnitude (G1)
      real*4        e_8_0mag_1 (nr__1) ! (mag) 8.0mag uncertainty (G1)
      real*4        v24mag_1   (nr__1) ! (mag) ? Spitzer/MIPS 24.0um magnitude (G1) (3)
      real*4        e_24mag_1  (nr__1) ! (mag) ? 24mag uncertainty (G1) (3)
      integer*4     Sep        (nr__1) ! (arcsec) YSO separation from center of HII region (4)
      character*1   m          (nr__1) ! [1-3m-] Multiplicity (5)
      character*1   n_Name     (nr__1) ! [f-mF-W] Note on YSOs. See notes.dat file (6)
*Note (1): YSOs without 24um photometry have a "C" appended to their number in
*          these tables, including N63 YSOs 8C, 16C, and 31C, which have 24um
*          manual photometry.
*Note (2): The class "flat spectrum" is designated by "fs".
*Note (3): The 24um magnitudes listed here used the SSC aperture corrections,
*          which likely underestimate fluxes by a factor of ~20% (see Sect. 2.2).
*Note (4): separation of YSOs from the approximate center of
*    * NGC 2030 (=LHA 120-N 63) at 05:35:32.3-66:02:38 (J2000.0), 
*      i.e., on star 2 (Fig. 6) in table5;
*    * the brightest ionizing star of NGC 2122 at 05:48:49.6-70:03:58 (J2000.0) 
*      (=HD 270145), i.e. star 118 (Fig. 9) in table6
*Note (5): Number of stars seen on J, Ks ISPI images within 3" of the YSO;
*          "m" is for several stars.
*Note (6): Lower case letter flag for table 5, upper case letter for table 6.

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

C  Declarations for 'table6.dat'	! Properties of the YSOs of N180

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

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

      character*3   Name_1     (nr__2) ! YSO name (see Nomenclature notes) (1)
      character*2   cl_2       (nr__2) ! [1/2/fs ] YSO class (2)
      integer*4     RAh_2      (nr__2) ! (h) Hour of right ascension (J2000)
      integer*4     RAm_2      (nr__2) ! (min) Minute of right ascension (J2000)
      real*4        RAs_2      (nr__2) ! (s) Second of right ascension (J2000)
      character*1   DE__2      (nr__2) ! Declination sign (J2000)
      integer*4     DEd_2      (nr__2) ! (deg) Degree of declination (J2000)
      integer*4     DEm_2      (nr__2) ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs_2      (nr__2) ! (arcsec) Arcsecond of declination (J2000)
      character*1   f_Jmag_1   (nr__2) ! [*] This paper ISPI Jmag; 2MASS otherwise
      real*4        Jmag_2     (nr__2) ! (mag) ? 2MASS J-band magnitude
      real*4        e_Jmag_2   (nr__2) ! (mag) ? Jmag uncertainty
      real*4        Hmag_2     (nr__2) ! (mag) ? 2MASS H-band magnitude
      real*4        e_Hmag_2   (nr__2) ! (mag) ? Hmag uncertainty
      character*1   f_Ksmag_1  (nr__2) ! [*] This paper ISPI Ksmag; 2MASS otherwise
      real*4        Ksmag_2    (nr__2) ! (mag) ? 2MASS Ks-band magnitude
      real*4        e_Ksmag_2  (nr__2) ! (mag) ? Ksmag uncertainty
      real*4        v3_6mag_2  (nr__2) ! (mag) Spitzer/IRAC 3.6um magnitude (G1)
      real*4        e_3_6mag_2 (nr__2) ! (mag) 3.6mag uncertainty (G1)
      real*4        v4_5mag_2  (nr__2) ! (mag) Spitzer/IRAC 4.5um magnitude (G1)
      real*4        e_4_5mag_2 (nr__2) ! (mag) 4.5mag uncertainty (G1)
      real*4        v5_8mag_2  (nr__2) ! (mag) Spitzer/IRAC 5.8um magnitude (G1)
      real*4        e_5_8mag_2 (nr__2) ! (mag) 5.8mag uncertainty (G1)
      real*4        v8_0mag_2  (nr__2) ! (mag) Spitzer/IRAC 8.0um magnitude (G1)
      real*4        e_8_0mag_2 (nr__2) ! (mag) 8.0mag uncertainty (G1)
      real*4        v24mag_2   (nr__2) ! (mag) ? Spitzer/MIPS 24.0um magnitude (G1) (3)
      real*4        e_24mag_2  (nr__2) ! (mag) ? 24mag uncertainty (G1) (3)
      integer*4     Sep_1      (nr__2) ! (arcsec) YSO separation from center of HII region (4)
      character*1   m_1        (nr__2) ! [1-3m-] Multiplicity (5)
      character*1   n_Name_1   (nr__2) ! [f-mF-W] Note on YSOs. See notes.dat file (6)
*Note (1): YSOs without 24um photometry have a "C" appended to their number in
*          these tables, including N63 YSOs 8C, 16C, and 31C, which have 24um
*          manual photometry.
*Note (2): The class "flat spectrum" is designated by "fs".
*Note (3): The 24um magnitudes listed here used the SSC aperture corrections,
*          which likely underestimate fluxes by a factor of ~20% (see Sect. 2.2).
*Note (4): separation of YSOs from the approximate center of
*    * NGC 2030 (=LHA 120-N 63) at 05:35:32.3-66:02:38 (J2000.0), 
*      i.e., on star 2 (Fig. 6) in table5;
*    * the brightest ionizing star of NGC 2122 at 05:48:49.6-70:03:58 (J2000.0) 
*      (=HD 270145), i.e. star 118 (Fig. 9) in table6
*Note (5): Number of stars seen on J, Ks ISPI images within 3" of the YSO;
*          "m" is for several stars.
*Note (6): Lower case letter flag for table 5, upper case letter for table 6.

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

C  Declarations for 'notes.dat'	! Notes on YSOs

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

      character*1   Flag       (nr__3) ! Note designation
      character*74  Note       (nr__3) ! Text of note

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

C  Loading file 'table3.dat'	! Galaxies and AGB stars

C  Format for file interpretation

    1 format(
     +  A9,A1,1X,A3,A1,1X,I1,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,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,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.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,18
        read(1,'(A141)')ar__
        read(ar__,1)
     +  Field(i__),f_Field(i__),cl(i__),f_cl(i__),RAh(i__),RAm(i__),
     +  RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),Jmag(i__),
     +  e_Jmag(i__),Hmag(i__),e_Hmag(i__),Ksmag(i__),e_Ksmag(i__),
     +  v3_6mag(i__),e_3_6mag(i__),v4_5mag(i__),e_4_5mag(i__),
     +  v5_8mag(i__),e_5_8mag(i__),v8_0mag(i__),e_8_0mag(i__),
     +  v24mag(i__),e_24mag(i__)
        if(ar__(40:45) .EQ. '') Jmag(i__) = rNULL__
        if(ar__(47:51) .EQ. '') e_Jmag(i__) = rNULL__
        if(ar__(53:58) .EQ. '') Hmag(i__) = rNULL__
        if(ar__(60:64) .EQ. '') e_Hmag(i__) = rNULL__
        if(ar__(66:71) .EQ. '') Ksmag(i__) = rNULL__
        if(ar__(73:77) .EQ. '') e_Ksmag(i__) = rNULL__
        if(ar__(131:135) .EQ. '') v24mag(i__) = rNULL__
        if(ar__(137:141) .EQ. '') e_24mag(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)
     +  Field(i__),f_Field(i__),cl(i__),f_cl(i__),RAh(i__),RAm(i__),
     +  RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),Jmag(i__),
     +  e_Jmag(i__),Hmag(i__),e_Hmag(i__),Ksmag(i__),e_Ksmag(i__),
     +  v3_6mag(i__),e_3_6mag(i__),v4_5mag(i__),e_4_5mag(i__),
     +  v5_8mag(i__),e_5_8mag(i__),v8_0mag(i__),e_8_0mag(i__),
     +  v24mag(i__),e_24mag(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 'table5.dat'	! Properties of the YSOs of N63

C  Format for file interpretation

    2 format(
     +  A3,1X,A2,1X,I1,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A1,
     +  F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,A1,F6.3,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,I3,1X,A1,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,29
        read(1,'(A141)')ar__1
        read(ar__1,2)
     +  Name(i__),cl_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
     +  DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),f_Jmag(i__),
     +  Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),
     +  f_Ksmag(i__),Ksmag_1(i__),e_Ksmag_1(i__),v3_6mag_1(i__),
     +  e_3_6mag_1(i__),v4_5mag_1(i__),e_4_5mag_1(i__),v5_8mag_1(i__),
     +  e_5_8mag_1(i__),v8_0mag_1(i__),e_8_0mag_1(i__),v24mag_1(i__),
     +  e_24mag_1(i__),Sep(i__),m(i__),n_Name(i__)
        if(ar__1(32:37) .EQ. '') Jmag_1(i__) = rNULL__
        if(ar__1(39:43) .EQ. '') e_Jmag_1(i__) = rNULL__
        if(ar__1(45:50) .EQ. '') Hmag_1(i__) = rNULL__
        if(ar__1(52:56) .EQ. '') e_Hmag_1(i__) = rNULL__
        if(ar__1(59:64) .EQ. '') Ksmag_1(i__) = rNULL__
        if(ar__1(66:70) .EQ. '') e_Ksmag_1(i__) = rNULL__
        if(ar__1(124:128) .EQ. '') v24mag_1(i__) = rNULL__
        if(ar__1(130:134) .EQ. '') e_24mag_1(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)
     +  Name(i__),cl_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
     +  DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),f_Jmag(i__),
     +  Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),
     +  f_Ksmag(i__),Ksmag_1(i__),e_Ksmag_1(i__),v3_6mag_1(i__),
     +  e_3_6mag_1(i__),v4_5mag_1(i__),e_4_5mag_1(i__),v5_8mag_1(i__),
     +  e_5_8mag_1(i__),v8_0mag_1(i__),e_8_0mag_1(i__),v24mag_1(i__),
     +  e_24mag_1(i__),Sep(i__),m(i__),n_Name(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 'table6.dat'	! Properties of the YSOs of N180

C  Format for file interpretation

    3 format(
     +  A3,1X,A2,1X,I1,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A1,
     +  F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,A1,F6.3,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,I3,1X,A1,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,38
        read(1,'(A141)')ar__2
        read(ar__2,3)
     +  Name_1(i__),cl_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),
     +  DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),f_Jmag_1(i__),
     +  Jmag_2(i__),e_Jmag_2(i__),Hmag_2(i__),e_Hmag_2(i__),
     +  f_Ksmag_1(i__),Ksmag_2(i__),e_Ksmag_2(i__),v3_6mag_2(i__),
     +  e_3_6mag_2(i__),v4_5mag_2(i__),e_4_5mag_2(i__),v5_8mag_2(i__),
     +  e_5_8mag_2(i__),v8_0mag_2(i__),e_8_0mag_2(i__),v24mag_2(i__),
     +  e_24mag_2(i__),Sep_1(i__),m_1(i__),n_Name_1(i__)
        if(ar__2(32:37) .EQ. '') Jmag_2(i__) = rNULL__
        if(ar__2(39:43) .EQ. '') e_Jmag_2(i__) = rNULL__
        if(ar__2(45:50) .EQ. '') Hmag_2(i__) = rNULL__
        if(ar__2(52:56) .EQ. '') e_Hmag_2(i__) = rNULL__
        if(ar__2(59:64) .EQ. '') Ksmag_2(i__) = rNULL__
        if(ar__2(66:70) .EQ. '') e_Ksmag_2(i__) = rNULL__
        if(ar__2(124:128) .EQ. '') v24mag_2(i__) = rNULL__
        if(ar__2(130:134) .EQ. '') e_24mag_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,3)
     +  Name_1(i__),cl_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),
     +  DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),f_Jmag_1(i__),
     +  Jmag_2(i__),e_Jmag_2(i__),Hmag_2(i__),e_Hmag_2(i__),
     +  f_Ksmag_1(i__),Ksmag_2(i__),e_Ksmag_2(i__),v3_6mag_2(i__),
     +  e_3_6mag_2(i__),v4_5mag_2(i__),e_4_5mag_2(i__),v5_8mag_2(i__),
     +  e_5_8mag_2(i__),v8_0mag_2(i__),e_8_0mag_2(i__),v24mag_2(i__),
     +  e_24mag_2(i__),Sep_1(i__),m_1(i__),n_Name_1(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'notes.dat'	! Notes on YSOs

C  Format for file interpretation

    4 format(A1,1X,A74)

C  Effective file loading

      open(unit=1,status='old',file=
     +'notes.dat')
      write(6,*) '....Loading file: notes.dat'
      do i__=1,42
        read(1,'(A76)')ar__3
        read(ar__3,4)Flag(i__),Note(i__)
c    ..............Just test output...........
        write(6,4)Flag(i__),Note(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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