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