Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/498/313 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 2024-Apr-20 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/498/313 Photometry of 3 main belt asteroids (Marciniak+, 2009) *================================================================================ *Photometry and models of selected main belt asteroids. *VI. 160 Una, 747 Winchester, and 849 Ara. * Marciniak A., Michalowski T., Hirsch R., Polinska M., * Kaminski K., Kwiatkowski T., Kryszczynska A., Behrend R., * Bernasconi L., Michalowski J., Starczewski S., Fagas M., Sobkowiak K. * <Astron. Astrophys. 498, 313 (2009)> * =2009A&A...498..313M C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'objects.dat' ! Asteroid fundamental parameters integer*4 nr__ parameter (nr__=3) ! Number of records character*68 ar__ ! Full-size record integer*4 Ast (nr__) ! Asteroid number character*10 Name (nr__) ! Asteroid name real*4 H (nr__) ! (mag) Absolute magnitude H parameter (1) real*4 Diam (nr__) ! (km) IRAS diameter (1) real*8 i (nr__) ! (deg) Inclination (1) real*8 e (nr__) ! Eccentricity (1) real*8 a (nr__) ! (AU) Semimajor axis (1) character*10 FileName (nr__) ! Name of the file containing the light curve *Note (1): Data from The Asteroid Orbital Elements Database, Cat. <B/astorb>. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'aspdat.dat' ! Aspect data of 160 Una, 747 Winchester and 849 Ara integer*4 nr__1 parameter (nr__1=61) ! Number of records character*67 ar__1 ! Full-size record C J2000 position composed of: Elon Elat integer*4 Ast_1 (nr__1) ! Asteroid number character*10 Name_1 (nr__1) ! Asteroid name integer*4 Obs_Y (nr__1) ! (yr) Year of an observation integer*4 Obs_M (nr__1) ! (m) Month of an observation real*4 Obs_D (nr__1) ! (d) Mid-time of an observation in days real*4 r (nr__1) ! (AU) Distance from object to the Sun real*4 delta (nr__1) ! (AU) Distance from object to the Earth real*4 phAngle (nr__1) ! (deg) Sun-asteroid-Earth phase angle real*4 Elon (nr__1) ! (deg) J2000 ecliptic longitude of object ({lambda}) real*4 Elat (nr__1) ! (deg) J2000 ecliptic latitude of object ({beta}) character*4 Site (nr__1) ! Observation site (1) *Note (1): Observatory codes as follows: * Bor = Borowiec * EnO = Les Engarouines Observatory * SAAO = South African Astronomical Observatory c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for '160lcs.dat' ! Asteroid 160 Una individual lightcurves integer*4 nr__2 parameter (nr__2=933) ! Number of records character*111 ar__2 ! Full-size record real*8 JD (nr__2) ! (d) JD epoch corrected for the light-time * corresponding to the Earth-asteroid distance real*4 br (nr__2) ! Relative brightness in intensity units, * mean brightness of each lightcurve is unity real*4 Sx (nr__2) ! (AU) x component of the vector from the asteroid to * the Sun in J2000 ecliptic Cartesian coordinates real*4 Sy (nr__2) ! (AU) y component of the vector from the asteroid to * the Sun in J2000 ecliptic Cartesian coordinates real*4 Sz (nr__2) ! (AU) z component of the vector from the asteroid to * the Sun in J2000 ecliptic Cartesian coordinates real*4 Ex (nr__2) ! (AU) x component of the vector from the asteroid to * the Earth in J2000 ecliptic Cartesian coordinates real*4 Ey (nr__2) ! (AU) y component of the vector from the asteroid to * the Earth in J2000 ecliptic Cartesian coordinates real*4 Ez (nr__2) ! (AU) z component of the vector from the asteroid to * the Earth in J2000 ecliptic Cartesian coordinates c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for '747lcs.dat' ! Asteroid 747 Winchester individual lightcurves integer*4 nr__3 parameter (nr__3=2173) ! Number of records character*111 ar__3 ! Full-size record real*8 JD_1 (nr__3) ! (d) JD epoch corrected for the light-time * corresponding to the Earth-asteroid distance real*4 br_1 (nr__3) ! Relative brightness in intensity units, * mean brightness of each lightcurve is unity real*4 Sx_1 (nr__3) ! (AU) x component of the vector from the asteroid to * the Sun in J2000 ecliptic Cartesian coordinates real*4 Sy_1 (nr__3) ! (AU) y component of the vector from the asteroid to * the Sun in J2000 ecliptic Cartesian coordinates real*4 Sz_1 (nr__3) ! (AU) z component of the vector from the asteroid to * the Sun in J2000 ecliptic Cartesian coordinates real*4 Ex_1 (nr__3) ! (AU) x component of the vector from the asteroid to * the Earth in J2000 ecliptic Cartesian coordinates real*4 Ey_1 (nr__3) ! (AU) y component of the vector from the asteroid to * the Earth in J2000 ecliptic Cartesian coordinates real*4 Ez_1 (nr__3) ! (AU) z component of the vector from the asteroid to * the Earth in J2000 ecliptic Cartesian coordinates c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for '849lcs.dat' ! Asteroid 849 Ara individual lightcurves integer*4 nr__4 parameter (nr__4=480) ! Number of records character*111 ar__4 ! Full-size record real*8 JD_2 (nr__4) ! (d) JD epoch corrected for the light-time * corresponding to the Earth-asteroid distance real*4 br_2 (nr__4) ! Relative brightness in intensity units, * mean brightness of each lightcurve is unity real*4 Sx_2 (nr__4) ! (AU) x component of the vector from the asteroid to * the Sun in J2000 ecliptic Cartesian coordinates real*4 Sy_2 (nr__4) ! (AU) y component of the vector from the asteroid to * the Sun in J2000 ecliptic Cartesian coordinates real*4 Sz_2 (nr__4) ! (AU) z component of the vector from the asteroid to * the Sun in J2000 ecliptic Cartesian coordinates real*4 Ex_2 (nr__4) ! (AU) x component of the vector from the asteroid to * the Earth in J2000 ecliptic Cartesian coordinates real*4 Ey_2 (nr__4) ! (AU) y component of the vector from the asteroid to * the Earth in J2000 ecliptic Cartesian coordinates real*4 Ez_2 (nr__4) ! (AU) z component of the vector from the asteroid to * the Earth in J2000 ecliptic Cartesian coordinates C============================================================================= C Loading file 'objects.dat' ! Asteroid fundamental parameters C Format for file interpretation 1 format( + I3,1X,A10,1X,F4.2,1X,F5.1,1X,F9.6,1X,F10.8,1X,F10.8,1X,A10) C Effective file loading open(unit=1,status='old',file= +'objects.dat') write(6,*) '....Loading file: objects.dat' do i__=1,3 read(1,'(A68)')ar__ read(ar__,1) + Ast(i__),Name(i__),H(i__),Diam(i__),i(i__),e(i__),a(i__), + FileName(i__) c ..............Just test output........... write(6,1) + Ast(i__),Name(i__),H(i__),Diam(i__),i(i__),e(i__),a(i__), + FileName(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'aspdat.dat' ! Aspect data of 160 Una, 747 Winchester and * 849 Ara C Format for file interpretation 2 format( + I3,1X,A10,1X,I4,1X,I2,1X,F4.1,2X,F6.4,1X,F6.4,1X,F5.2,1X,F6.2, + 1X,F6.2,1X,A4) C Effective file loading open(unit=1,status='old',file= +'aspdat.dat') write(6,*) '....Loading file: aspdat.dat' do i__=1,61 read(1,'(A67)')ar__1 read(ar__1,2) + Ast_1(i__),Name_1(i__),Obs_Y(i__),Obs_M(i__),Obs_D(i__), + r(i__),delta(i__),phAngle(i__),Elon(i__),Elat(i__),Site(i__) c ..............Just test output........... write(6,2) + Ast_1(i__),Name_1(i__),Obs_Y(i__),Obs_M(i__),Obs_D(i__), + r(i__),delta(i__),phAngle(i__),Elon(i__),Elat(i__),Site(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file '160lcs.dat' ! Asteroid 160 Una individual lightcurves C Format for file interpretation 3 format( + F14.6,1X,E12.7,1X,E13.8,1X,E13.8,1X,E13.6,1X,E13.8,1X,E13.8, + 1X,E13.8) C Effective file loading open(unit=1,status='old',file= +'160lcs.dat') write(6,*) '....Loading file: 160lcs.dat' do i__=1,933 read(1,'(A111)')ar__2 read(ar__2,3) + JD(i__),br(i__),Sx(i__),Sy(i__),Sz(i__),Ex(i__),Ey(i__), + Ez(i__) c ..............Just test output........... write(6,3) + JD(i__),br(i__),Sx(i__),Sy(i__),Sz(i__),Ex(i__),Ey(i__), + Ez(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file '747lcs.dat' ! Asteroid 747 Winchester individual lightcurves C Format for file interpretation 4 format( + F14.6,1X,E12.7,1X,E13.8,1X,E13.8,1X,E13.6,1X,E13.8,1X,E13.8, + 1X,E13.8) C Effective file loading open(unit=1,status='old',file= +'747lcs.dat') write(6,*) '....Loading file: 747lcs.dat' do i__=1,2173 read(1,'(A111)')ar__3 read(ar__3,4) + JD_1(i__),br_1(i__),Sx_1(i__),Sy_1(i__),Sz_1(i__),Ex_1(i__), + Ey_1(i__),Ez_1(i__) c ..............Just test output........... write(6,4) + JD_1(i__),br_1(i__),Sx_1(i__),Sy_1(i__),Sz_1(i__),Ex_1(i__), + Ey_1(i__),Ez_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file '849lcs.dat' ! Asteroid 849 Ara individual lightcurves C Format for file interpretation 5 format( + F14.6,1X,E12.7,1X,E13.8,1X,E13.8,1X,E13.6,1X,E13.8,1X,E13.8, + 1X,E13.8) C Effective file loading open(unit=1,status='old',file= +'849lcs.dat') write(6,*) '....Loading file: 849lcs.dat' do i__=1,480 read(1,'(A111)')ar__4 read(ar__4,5) + JD_2(i__),br_2(i__),Sx_2(i__),Sy_2(i__),Sz_2(i__),Ex_2(i__), + Ey_2(i__),Ez_2(i__) c ..............Just test output........... write(6,5) + JD_2(i__),br_2(i__),Sx_2(i__),Sy_2(i__),Sz_2(i__),Ex_2(i__), + Ey_2(i__),Ez_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end