FORTRAN Generation
(/./ftp/cats/J/A_A/541/A142)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/541/A142 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.



++++Column 'RAdeg' at line 214: Right ascension in decimal degrees (J2000)
                                 at mean epoch
    position has no associated Epoch (last analysed word: '')
++++Column 'RAh' at line 294: Right ascension (J2000) of star at epoch
    position has no associated Epoch (last analysed word: '')
      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Mar-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/A+A/541/A142    Prediction of stellar occultations 2008-2015  (Assafin+, 2012)
*================================================================================
*Candidate stellar occultations by large TNOs up to 2015.
*    Assafin M., Camargo J.I.B., Vieira Martins R., Braga-Ribas F., Sicardy B.,
*    Andrei A.H., da Silva Neto D.N.
*   <Astron. Astrophys. 541, A142 (2012)>
*   =2012A&A...541A.142A
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'list.dat'	! List of objects with catalog star positions sky path
                             and prediction stellar occultations

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

      integer*4     Year       (nr__) ! (yr) Year of occultation
      character*10  Object     (nr__) ! Object name
      character*12  CatFile    (nr__) ! Name of the file with catalog star positions
*                                  sky path, in "catal" subdirectory
      character*12  TabFile    (nr__) ! Name of the file with prediction stellar
*                                  occultations, in "tables" subdirectory
      character*24  Com        (nr__) ! Comment

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

C  Declarations for 'catal/*'	! Individual files of catalog star positions sky path

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg      (nr__1) ! (deg) Right ascension in decimal degrees (J2000)
*                                 at mean epoch
      real*8        DEdeg      (nr__1) ! (deg) Declination in decimal degrees (J2000)
*                                 at mean epoch
      real*8        JD         (nr__1) ! (d) Mean epoch (Julian Date) of coordinates
      integer*4     e_RAs      (nr__1) ! (mas) ?=999999 Right ascension error at mean epoch (1)
      integer*4     e_DEs      (nr__1) ! (mas) ?=999999 Declination error at mean epoch (1)
      integer*4     pmRA       (nr__1) ! (mas/yr) ?=999999 Proper motion in Right Ascension
*                                          (RA*cosDE) (2)
      integer*4     pmDE       (nr__1) ! (mas/yr) ?=999999 Proper motion in Declination (2)
      real*4        Rmag       (nr__1) ! (mag) Magnitude (apparent) in R band
      real*4        Jmag       (nr__1) ! (mag) Magnitude (apparent) in J band (3)
      real*4        Hmag       (nr__1) ! (mag) Magnitude (apparent) in H band (3)
      real*4        Kmag       (nr__1) ! (mag) Magnitude (apparent) in K band (3)
      integer*4     e_x        (nr__1) ! (mas) Measured error in CCD frames (x) coordinates
      integer*4     e_y        (nr__1) ! (mas) Measured error in CCD frames (y) coordinates
      integer*4     o_Obs      (nr__1) ! Number of observations per position
      integer*4     f_cat      (nr__1) ! Flag on catalog cross-identification (4)
      integer*4     f_Pos      (nr__1) ! Astrometric multiplicity flag (5)
*Note (1): Error values marked "999999" mean no estimated errors; this regards
*     to astrometric flagged positions or to positions from one observation
*Note (2): proper motion values marked "999999" mean no computed proper motions;
*     this is because no first epoch position counterpart could be found in
*     the 2MASS or USNO B1.0 catalogs
*Note (3): J, H or K magnitudes marked "99.999" mean no infrared magnitudes
*     available; this means that it is not a 2MASS star
*Note (4): Flag on catalog cross-identification as follows:
*   1 = UCAC2 star
*   2 = 2MASS star
*   9 = not a UCAC2/2MASS star
*Note (5): Astrometric multiplicity flags range from 0 to 6 as follows:
*   0 = final WFI mosaic star position came from all contributing
*       individual CCD positions within 0.2arcsec from each other,
*       after the final WFI mosaic solutions; no nearby individual CCD
*       position was identified within 1.5arcsec of the final catalog
*       star position, indicating good astrometry
*   1 = a number of individual CCD positions was identified within
*       1.5arcsec of each other; more than one of these individual
*       positions was labeled as a UCAC2 star; the final catalog
*       position assigned for this star came from the average over these
*       individual CCD positions labeled as UCAC2; the same flag is
*       assigned if it refers to the 2MASS catalog (one can distinguish
*       between both cases looking at the catalog flag)
*   2 = a number of individual CCD positions was identified within
*       1.5arcsec of each other; but only one single of these positions
*       was labeled UCAC2/2MASS; this single position was assigned as
*       the final catalog star position (the catalog flag indicates
*       which catalog is involved - UCAC2 or 2MASS)
*   3 = some contributing individual CCD positions were within 0.2arcsec
*       from each other after the final WFI mosaic solutions, but the 
*       resulting positions were nearby others within 1.5arcsec; none of
*       these positions were labeled UCAC2 or 2MASS; here, the final 
*       catalog star position assigned for this star came from the entry
*       with highest number (N not equal to 1) of individual position
*       contributions from the WFI mosaic solution
*   4 = same situation as in flag "3", but more than one entry presented
*       the same highest number of individual position contributions;
*       in this case, the final catalog star position assigned for this
*       star came from the single entry with least (x,y) measured errors
*   5 = same situation as in flag "4", but more than one entry presented
*       the same best (x,y) measurement errors; in this case, the final
*       catalog star position assigned for this star came from the
*       single entry with brightest R magnitude
*   6 = same situation as in flag "5", but more than one entry presented
*       the same brightest R magnitudes; in this case, the final catalog
*       star position assigned for this star came from the average over
*       all detected multiple entries

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

C  Declarations for 'tables/*'	! Individual files for prediction stellar occultations

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

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

      integer*4     App_Y      (nr__2) ! (yr) UTC instant at closest approach (year)
      integer*4     App_M      (nr__2) ! ("month") UTC instant at closest approach (month)
      integer*4     App_D      (nr__2) ! (d) UTC instant at closest approach (day)
      integer*4     App_h      (nr__2) ! (h) UTC instant at closest approach (hours)
      integer*4     App_m_1    (nr__2) ! (min) UTC instant at closest approach (minutes)
      real*4        App_s      (nr__2) ! (s) UTC instant at closest approach (seconds)
      integer*4     RAh        (nr__2) ! (h) Right ascension (J2000) of star at epoch
      integer*4     RAm        (nr__2) ! (min) Right ascension (J2000) of star at epoch
      real*8        RAs        (nr__2) ! (s) Right ascension (J2000) of star at epoch
      character*1   DE_        (nr__2) ! Declination sign (J2000) of star at epoch
      integer*4     DEd        (nr__2) ! (deg) Declination (J2000) of star at epoch
      integer*4     DEm        (nr__2) ! (arcmin) Declination (J2000) of star at epoch
      real*4        DEs        (nr__2) ! (arcsec) Declination (J2000) of star at epoch
      integer*4     RAoh       (nr__2) ! (h) TNO geocentric right ascension (J2000)
*                                 at occultation
      integer*4     RAom       (nr__2) ! (min) TNO geocentric right ascension (J2000)
*                                 at occultation
      real*8        RAos       (nr__2) ! (s) TNO geocentric right ascension (J2000)
*                                 at occultation
      character*1   DEo_       (nr__2) ! TNO declination sign (J2000) at occultation
      integer*4     DEod       (nr__2) ! (deg) TNO declination (J2000) at occultation
      integer*4     DEom       (nr__2) ! (arcmin) TNO declination (J2000) at occultation
      real*4        DEos       (nr__2) ! (arcsec) TNO declination (J2000) at occultation
      real*4        Sep        (nr__2) ! (arcsec) Minimum Separation at closest approach (CA)
      real*4        PA         (nr__2) ! (deg) Position angle (PA) at closest approach (1)
      real*8        Vshad      (nr__2) ! (km/s) Occultation shadow velocity (2)
      real*4        Dist       (nr__2) ! (AU) Apparent geocentric distance at occultation
      real*4        Rmag_1     (nr__2) ! (mag) Normalized magnitude (apparent) in R band
      real*4        Jmag_1     (nr__2) ! (mag) Normalized magnitude (apparent) in J band (3)
      real*4        Hmag_1     (nr__2) ! (mag) Normalized magnitude (apparent) in H band (3)
      real*4        Kmag_1     (nr__2) ! (mag) Normalized magnitude (apparent) in K band (3)
      real*4        Long       (nr__2) ! (deg) Longitude of the sub-solar point at occultation
      integer*4     LST_h      (nr__2) ! (h) Local solar time at closest approach (hours)
      integer*4     LST_m      (nr__2) ! (min) Local solar time at closest approach (minutes)
      real*8        oRA        (nr__2) ! (mas) Right ascension ephemeris offset correction
      real*8        oDE        (nr__2) ! (mas) Declination ephemeris offset correction
      character*2   f_pm       (nr__2) ! [ok/no] Proper motion flag (4)
      character*2   f_cat_1    (nr__2) ! Catalog cross-identification flag (5)
      integer*4     f_Pos_1    (nr__2) ! Astrometric multiplicity flag (6)
      integer*4     e_RAs_1    (nr__2) ! (mas) ?=9999 Right ascension error at catalog
*                                       mean epoch (7)
      integer*4     e_DEs_1    (nr__2) ! (mas) ?=9999 Declination error at catalog
*                                       mean epoch (7)
      integer*4     pmRA_1     (nr__2) ! (mas/yr) ?=9999 Proper motion in Right Ascension
*                                       (RA*cosDE) (8)
      integer*4     pmDE_1     (nr__2) ! (mas/yr) ?=9999 Proper motion in Declination (8)
*Note (1): This is the position angle of the body with respect to the star
*     at closest approach; PA is zero when the body is north of the star and
*     is counted clockwise
*Note (2): This is the occultation shadow velocity across the Earth at the
*     apparent geocentric distance of the body; it has a positive/negative
*     sign for prograde/retrograde velocities, that is, when the TNO's
*     geocentric right ascension is respectively increasing/decreasing
*Note (3): J, H or K magnitudes marked "50.0" mean no infrared magnitudes
*     available; this means that it is not a 2MASS star
*Note (4): Star has proper motion? If yes, flag = "ok", otherwise flag = "no"
*Note (5): Flag as follows:
*  uc = UCAC2 star
*  2m = 2MASS star
*  fs = field star (not a UCAC2 or 2MASS catalog star)
*Note (6): Astrometric multiplicity flags range from 0 to 6 as follows:
*   0 = final WFI mosaic star position ok; all contributing individual
*       CCD positions were within 0.2arcsec from each other after the
*       final WFI mosaic solutions; no nearby individual CCD position
*       was identified within 1.5arcsec of the final catalog star
*       position, indicating good astrometry
*   1 = a number of individual CCD positions was identified within
*       1.5arcsec of each other; more than one of these individual
*       positions was labeled as a UCAC2 star; the final catalog
*       position assigned for this star came from the average over these
*       individual CCD positions labeled as UCAC2; the same flag is
*       assigned if it refers to the 2MASS catalog (one can distinguish
*       between both cases looking at the catalog flag)
*   2 = a number of individual CCD positions was identified within
*       1.5arcsec of each other; but only one single of these positions
*       was labeled UCAC2/2MASS; this single position was assigned as
*       the final catalog star position (the catalog flag indicates
*       which catalog is involved - UCAC2 or 2MASS)
*   3 = some contributing individual CCD positions were within 0.2arcsec
*       from each other after the final WFI mosaic solutions, but the
*       resulting positions were nearby others within 1.5arcsec; none of
*       these positions were labeled UCAC2 or 2MASS; here, the final
*       catalog star position assigned for this star came from the entry
*       with highest number (N not equal to 1) of individual position
*       contributions from the WFI mosaic solution
*   4 = same situation as in flag "3", but more than one entry presented
*       the same highest number of individual position contributions; in
*       this case, the final catalog star position assigned for this star
*       came from the single entry with least (x,y) measured errors
*   5 = same situation as in flag "4", but more than one entry presented
*       the same best (x,y) measurement errors; in this case, the final
*       catalog star position assigned for this star came from the single
*       entry with brightest R magnitude
*   6 = same situation as in flag "5", but more than one entry presented
*       the same brightest R magnitudes; in this case, the final catalog
*       star position assigned for this star came from the average over
*       all detected multiple entries
*Note (7): Error values marked "9999" mean no estimated errors; this regards
*     to astrometric flagged positions or to positions from one observation.
*Note (8): Proper motion values marked "9999" mean no computed proper motions;
*     this is because no first epoch position counterpart could be found in
*     the 2MASS or USNO B1.0 catalogs.

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

C  Loading file 'list.dat'	! List of objects with catalog star positions sky path
*                             and prediction stellar occultations

C  Format for file interpretation

    1 format(I4,1X,A10,1X,A12,1X,A12,1X,A24)

C  Effective file loading

      open(unit=1,status='old',file=
     +'list.dat')
      write(6,*) '....Loading file: list.dat'
      do i__=1,74
        read(1,'(A66)')ar__
        read(ar__,1)
     +  Year(i__),Object(i__),CatFile(i__),TabFile(i__),Com(i__)
c    ..............Just test output...........
        write(6,1)
     +  Year(i__),Object(i__),CatFile(i__),TabFile(i__),Com(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'catal/*'	! Individual files of catalog star positions sky path

C  Format for file interpretation

    2 format(
     +  1X,F13.9,1X,F13.9,1X,F16.8,1X,I6,1X,I6,1X,I6,1X,I6,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,F6.3,1X,I5,1X,I5,1X,I2,1X,I1,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catal/*')
      write(6,*) '....Loading file: catal/*'
      do i__=1,74
        read(1,'(A120)')ar__1
        read(ar__1,2)
     +  RAdeg(i__),DEdeg(i__),JD(i__),e_RAs(i__),e_DEs(i__),pmRA(i__),
     +  pmDE(i__),Rmag(i__),Jmag(i__),Hmag(i__),Kmag(i__),e_x(i__),
     +  e_y(i__),o_Obs(i__),f_cat(i__),f_Pos(i__)
c    ..............Just test output...........
        write(6,2)
     +  RAdeg(i__),DEdeg(i__),JD(i__),e_RAs(i__),e_DEs(i__),pmRA(i__),
     +  pmDE(i__),Rmag(i__),Jmag(i__),Hmag(i__),Kmag(i__),e_x(i__),
     +  e_y(i__),o_Obs(i__),f_cat(i__),f_Pos(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tables/*'	! Individual files for prediction stellar occultations

C  Format for file interpretation

    3 format(
     +  1X,I4,1X,I2,1X,I2,2X,I2,1X,I2,1X,F3.0,3X,I2,1X,I2,1X,F7.4,1X,
     +  A1,I2,1X,I2,2X,F6.3,3X,I2,1X,I2,1X,F7.4,1X,A1,I2,1X,I2,2X,
     +  F6.3,3X,F5.3,2X,F6.2,1X,F7.2,1X,F6.2,1X,F4.1,1X,F4.1,1X,F4.1,
     +  1X,F4.1,3X,F4.0,1X,I2,1X,I2,3X,F7.1,2X,F7.1,1X,A2,1X,A2,1X,I1,
     +  1X,I4,1X,I4,1X,I4,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tables/*')
      write(6,*) '....Loading file: tables/*'
      do i__=1,67
        read(1,'(A195)')ar__2
        read(ar__2,3)
     +  App_Y(i__),App_M(i__),App_D(i__),App_h(i__),App_m_1(i__),
     +  App_s(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__),RAoh(i__),RAom(i__),RAos(i__),DEo_(i__),
     +  DEod(i__),DEom(i__),DEos(i__),Sep(i__),PA(i__),Vshad(i__),
     +  Dist(i__),Rmag_1(i__),Jmag_1(i__),Hmag_1(i__),Kmag_1(i__),
     +  Long(i__),LST_h(i__),LST_m(i__),oRA(i__),oDE(i__),f_pm(i__),
     +  f_cat_1(i__),f_Pos_1(i__),e_RAs_1(i__),e_DEs_1(i__),
     +  pmRA_1(i__),pmDE_1(i__)
        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(i__) .GT. -180) RAdeg_1(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg_1(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c    ..............Just test output...........
        write(6,3)
     +  App_Y(i__),App_M(i__),App_D(i__),App_h(i__),App_m_1(i__),
     +  App_s(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__),RAoh(i__),RAom(i__),RAos(i__),DEo_(i__),
     +  DEod(i__),DEom(i__),DEos(i__),Sep(i__),PA(i__),Vshad(i__),
     +  Dist(i__),Rmag_1(i__),Jmag_1(i__),Hmag_1(i__),Kmag_1(i__),
     +  Long(i__),LST_h(i__),LST_m(i__),oRA(i__),oDE(i__),f_pm(i__),
     +  f_cat_1(i__),f_Pos_1(i__),e_RAs_1(i__),e_DEs_1(i__),
     +  pmRA_1(i__),pmDE_1(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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