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 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.



++++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        ! (yr) Year of occultation
      character*10  Object      ! Object name
      character*12  CatFile     ! Name of the file with catalog star positions
*                                  sky path, in "catal" subdirectory
      character*12  TabFile     ! Name of the file with prediction stellar
*                                  occultations, in "tables" subdirectory
      character*24  Com         ! 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       ! (deg) Right ascension in decimal degrees (J2000)
*                                 at mean epoch
      real*8        DEdeg       ! (deg) Declination in decimal degrees (J2000)
*                                 at mean epoch
      real*8        JD          ! (d) Mean epoch (Julian Date) of coordinates
      integer*4     e_RAs       ! (mas) ?=999999 Right ascension error at mean epoch (1)
      integer*4     e_DEs       ! (mas) ?=999999 Declination error at mean epoch (1)
      integer*4     pmRA        ! (mas/yr) ?=999999 Proper motion in Right Ascension
*                                          (RA*cosDE) (2)
      integer*4     pmDE        ! (mas/yr) ?=999999 Proper motion in Declination (2)
      real*4        Rmag        ! (mag) Magnitude (apparent) in R band
      real*4        Jmag        ! (mag) Magnitude (apparent) in J band (3)
      real*4        Hmag        ! (mag) Magnitude (apparent) in H band (3)
      real*4        Kmag        ! (mag) Magnitude (apparent) in K band (3)
      integer*4     e_x         ! (mas) Measured error in CCD frames (x) coordinates
      integer*4     e_y         ! (mas) Measured error in CCD frames (y) coordinates
      integer*4     o_Obs       ! Number of observations per position
      integer*4     f_cat       ! Flag on catalog cross-identification (4)
      integer*4     f_Pos       ! 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     ! (deg) Right Ascension J2000
      real*8        DEdeg_1     ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     App_Y       ! (yr) UTC instant at closest approach (year)
      integer*4     App_M       ! ("month") UTC instant at closest approach (month)
      integer*4     App_D       ! (d) UTC instant at closest approach (day)
      integer*4     App_h       ! (h) UTC instant at closest approach (hours)
      integer*4     App_m_1     ! (min) UTC instant at closest approach (minutes)
      real*4        App_s       ! (s) UTC instant at closest approach (seconds)
      integer*4     RAh         ! (h) Right ascension (J2000) of star at epoch
      integer*4     RAm         ! (min) Right ascension (J2000) of star at epoch
      real*8        RAs         ! (s) Right ascension (J2000) of star at epoch
      character*1   DE_         ! Declination sign (J2000) of star at epoch
      integer*4     DEd         ! (deg) Declination (J2000) of star at epoch
      integer*4     DEm         ! (arcmin) Declination (J2000) of star at epoch
      real*4        DEs         ! (arcsec) Declination (J2000) of star at epoch
      integer*4     RAoh        ! (h) TNO geocentric right ascension (J2000)
*                                 at occultation
      integer*4     RAom        ! (min) TNO geocentric right ascension (J2000)
*                                 at occultation
      real*8        RAos        ! (s) TNO geocentric right ascension (J2000)
*                                 at occultation
      character*1   DEo_        ! TNO declination sign (J2000) at occultation
      integer*4     DEod        ! (deg) TNO declination (J2000) at occultation
      integer*4     DEom        ! (arcmin) TNO declination (J2000) at occultation
      real*4        DEos        ! (arcsec) TNO declination (J2000) at occultation
      real*4        Sep         ! (arcsec) Minimum Separation at closest approach (CA)
      real*4        PA          ! (deg) Position angle (PA) at closest approach (1)
      real*8        Vshad       ! (km/s) Occultation shadow velocity (2)
      real*4        Dist        ! (AU) Apparent geocentric distance at occultation
      real*4        Rmag_1      ! (mag) Normalized magnitude (apparent) in R band
      real*4        Jmag_1      ! (mag) Normalized magnitude (apparent) in J band (3)
      real*4        Hmag_1      ! (mag) Normalized magnitude (apparent) in H band (3)
      real*4        Kmag_1      ! (mag) Normalized magnitude (apparent) in K band (3)
      real*4        Long        ! (deg) Longitude of the sub-solar point at occultation
      integer*4     LST_h       ! (h) Local solar time at closest approach (hours)
      integer*4     LST_m       ! (min) Local solar time at closest approach (minutes)
      real*8        oRA         ! (mas) Right ascension ephemeris offset correction
      real*8        oDE         ! (mas) Declination ephemeris offset correction
      character*2   f_pm        ! [ok/no] Proper motion flag (4)
      character*2   f_cat_1     ! Catalog cross-identification flag (5)
      integer*4     f_Pos_1     ! Astrometric multiplicity flag (6)
      integer*4     e_RAs_1     ! (mas) ?=9999 Right ascension error at catalog
*                                       mean epoch (7)
      integer*4     e_DEs_1     ! (mas) ?=9999 Declination error at catalog
*                                       mean epoch (7)
      integer*4     pmRA_1      ! (mas/yr) ?=9999 Proper motion in Right Ascension
*                                       (RA*cosDE) (8)
      integer*4     pmDE_1      ! (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,Object,CatFile,TabFile,Com
c    ..............Just test output...........
        write(6,1)Year,Object,CatFile,TabFile,Com
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,DEdeg,JD,e_RAs,e_DEs,pmRA,pmDE,Rmag,Jmag,Hmag,Kmag,e_x,
     +  e_y,o_Obs,f_cat,f_Pos
c    ..............Just test output...........
        write(6,2)
     +  RAdeg,DEdeg,JD,e_RAs,e_DEs,pmRA,pmDE,Rmag,Jmag,Hmag,Kmag,e_x,
     +  e_y,o_Obs,f_cat,f_Pos
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,App_M,App_D,App_h,App_m_1,App_s,RAh,RAm,RAs,DE_,DEd,DEm,
     +  DEs,RAoh,RAom,RAos,DEo_,DEod,DEom,DEos,Sep,PA,Vshad,Dist,
     +  Rmag_1,Jmag_1,Hmag_1,Kmag_1,Long,LST_h,LST_m,oRA,oDE,f_pm,
     +  f_cat_1,f_Pos_1,e_RAs_1,e_DEs_1,pmRA_1,pmDE_1
        RAdeg_1 = rNULL__
        DEdeg_1 = 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 .GT. -180) RAdeg_1=RAh*15.
        if(RAm .GT. -180) RAdeg_1=RAdeg_1+RAm/4.
        if(RAs .GT. -180) RAdeg_1=RAdeg_1+RAs/240.
        if(DEd .GE. 0) DEdeg_1=DEd
        if(DEm .GE. 0) DEdeg_1=DEdeg_1+DEm/60.
        if(DEs .GE. 0) DEdeg_1=DEdeg_1+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c    ..............Just test output...........
        write(6,3)
     +  App_Y,App_M,App_D,App_h,App_m_1,App_s,RAh,RAm,RAs,DE_,DEd,DEm,
     +  DEs,RAoh,RAom,RAos,DEo_,DEod,DEom,DEos,Sep,PA,Vshad,Dist,
     +  Rmag_1,Jmag_1,Hmag_1,Kmag_1,Long,LST_h,LST_m,oRA,oDE,f_pm,
     +  f_cat_1,f_Pos_1,e_RAs_1,e_DEs_1,pmRA_1,pmDE_1
        write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c    .......End.of.Just test output...........
      end do
      close(1)

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