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