Conversion of standardized ReadMe file for
file /./ftp/cats/VII/108C 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-18
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. VII/108C Asteroids II Machine-Readable Data Base (Binzel+ 1987)
*================================================================================
*Asteroids II
* Binzel R. P., Gehrels T., Matthews M. S. (ed)
* <The University Arizona Press, Tucson, p. 997 (1989)>
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'discover.dat' ! Names and discovery circumstances
integer*4 nr__
parameter (nr__=6213) ! Number of records
character*125 ar__ ! Full-size record
integer*4 Planet (nr__) ! Number of minor planet
character*17 name (nr__) ! Name of minor planet
character*9 pre_name (nr__) ! Preliminary name
integer*4 disc_Y (nr__) ! (yr) Year of discovery
character*3 disc_M (nr__) ! ("month") Month of discovery
integer*4 disc_D (nr__) ! (d) Day of discovery (in month)
character*29 dis_name (nr__) ! Name of discoverer
character*21 dis_pl (nr__) ! Place of discovery
integer*4 note (nr__) ! *?Reference to note
character*23 MPC (nr__) ! *Reference to Minor Planet Cir.
*Note on note:
* Notes referencing two or more discoverers with names of combined
* length too great to fit in the discoverer column, to give a more
* complete description of programs involving several persons, and to
* reference cases in which two numbered planets were subsequently
* discovered to be identical and the number and name of one of these was
* re-assigned to a newly discovered planet. Notes have also been used to
* reference conflicting discovery claims and list important independent
* discoveries which are no longer regarded as official.
*Note on MPC:
* References the Minor Planet Circular or Circulars in which information
* concerning the discovery circumstances and/or name assignment, and
* occasional error corrections, was first published. The author first
* listed these references only on 1988 Feb 2 with the publication of MPC
* 12781. Previously this table had been maintained without references to
* the original literature. Adding references from earlier publications
* remains a task for the future (if ever). Would anyone else like to
* contribute to this important assignment?
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'dis_note.dat' ! The notes for discover.dat
integer*4 nr__1
parameter (nr__1=315) ! Number of records
character*80 ar__1 ! Full-size record
integer*4 Note_1 (nr__1) ! *Note number
character*76 Text (nr__1) ! Text of Note
*Note on Note:
* The note number is repeated when several lines are necessary.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'lightcrv.dat' ! Lightcurve parameters
integer*4 nr__2
parameter (nr__2=528) ! Number of records
character*118 ar__2 ! Full-size record
integer*4 Planet_1 (nr__2) ! ?Asteroid number
character*16 name_1 (nr__2) ! Asteroid name
character*1 l_per (nr__2) ! [>] Limit sign for period
real*8 per (nr__2) ! (h) *?Period
character*1 l_mag1 (nr__2) ! [><]Limit sign for mag1
real*4 mag1 (nr__2) ! (mag) ?First mag., amplitude or variation
character*1 l_mag2 (nr__2) ! [-><]Separator or limit sign
real*4 mag2 (nr__2) ! (mag) ?Second mag., amplitude or variation
character*1 u_mag2 (nr__2) ! [? ] Uncertainty note on mag2
integer*4 n_mag (nr__2) ! *[1,4]? Reliability code
character*2 rem (nr__2) ! *Remarks code
character*69 refs (nr__2) ! *Reference codes, as a comma-separated list.
*Note on per:
* The period for 288 Glauke is 1150. and is in the bytes 28-32.
*Note on n_mag:
* The Reliability Codes are as follows:
* 1 = Very tentative result, may be completely wrong.
* 2 = Reasonably secure result, based on over half coverage of the
* lightcurve.
* 3 = Secure result, full lightcurve coverage, no ambiguity of period.
* 4 = Multiple apparition coverage, pole position reported.
*Note on rem:
* The Remarks Codes are as follows:
* 3 = Number of extrema per rotation cycle (e.g., 1, 3).
* Unless otherwise noted, two per cycle is assumed.
* A = Ambiguous period. The "most likely" period is listed, with other
* possibilities listed in a footnote ordered by asteroid number, below.
* D = Rotation period "determined" from published data,
* but not given by author(s) of original data.
* F = Footnote, below, ordered by asteroid number, containing
* additional information.
* N = No lightcurve published.
* P = Photographic photometry.
* V = Visual photometry.
*Note on refs:
* All references to the asteroid are cited by number in the reference
* file, in chronological order. The last few citations are the most
* recent, and should be the most useful for evaluating the current state
* of knowledge of a given asteroid.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ltcv_ref.dat' ! References file for lightcrv.dat
integer*4 nr__3
parameter (nr__3=497) ! Number of records
character*80 ar__3 ! Full-size record
integer*4 Ref (nr__3) ! *Reference number
character*75 Text_1 (nr__3) ! Reference
*Note on Ref:
* The reference number is repeated when several lines are necessary.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'proper.dat' ! Proper elements and family identities
integer*4 nr__4
parameter (nr__4=1796) ! Number of records
character*70 ar__4 ! Full-size record
integer*4 Planet_2 (nr__4) ! Asteroid number
real*4 a (nr__4) ! (AU) Semimajor axis of orbit
real*4 e (nr__4) ! ?Eccentricity
real*4 sin_i (nr__4) ! ?Sine of inclination
real*4 wbar (nr__4) ! (deg) ?Proper longitude of perihelion (1950.0)
real*4 anode (nr__4) ! (deg) ?Proper longitude, ascending node (1950.0)
real*4 dwbar (nr__4) ! (arcsec/a) ?Rate of longitude of perihelion
real*8 dnode (nr__4) ! (arcsec/a) ?Rate of longitude of node
character*3 res (nr__4) ! *Note indicating resonance or special case
real*4 dMars (nr__4) ! (AU) *?Closest colinear distance to Mars
real*4 dJup (nr__4) ! (AU) ?Closest colinear distance to Jupiter
character*1 q_Planet (nr__4) ! *Quality code for Palomar-Leiden astr.
character*3 family (nr__4) ! Family identity if appropriate
integer*4 Ref_1 (nr__4) ! [1,2] Reference number
*Note on res:
* Values 1-16 identify a secular resonance number. Commensurabilities
* are indicated as 1:1 for trojans, 3:2 for Hildas, etc.
* W L denotes an argument of perihelion librator.
* JXR = Jupiter crosser,
* SXR = Saturn crosser.
* This flag, when set, also means that rates and distances will have
* degraded accuracy (see notes in file "prop_doc.txt")
*Note on dMars:
* A negative value indicates a Mars crosser.
*Note on q_Planet:
* Quality code reserved for Palomar-Leiden asteroids. Values 1, 2,
* 3, 4, X, 0 are in decreasing order of accuracy. see
* Palomar-Leiden survey for details. A blank, not a zero, is used for
* numbered asteroids.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pole.dat' ! Pole determinations
integer*4 nr__5
parameter (nr__5=248) ! Number of records
character*75 ar__5 ! Full-size record
integer*4 Planet_3 (nr__5) ! Asteroid number
character*4 n_Planet (nr__5) ! *Source of basic data
integer*4 ELpole (nr__5) ! (deg) *?Ecliptic longitude of pole
character*1 u_pole (nr__5) ! [* ] indicates discordant positions
integer*4 EBpole (nr__5) ! (deg) *?Ecliptic latitude of pole
character*1 n_pole (nr__5) ! [* ] Possible parenthesis or discord
character*1 spin (nr__5) ! *[-NMS] Spin direction
character*1 u_spin (nr__5) ! [) ] Possible parenthesis
integer*4 EL2pole (nr__5) ! (deg) *? Ecliptic longitude of pole
integer*4 EB2pole (nr__5) ! (deg) *? Ecliptic latitude of pole
character*1 spin2 (nr__5) ! *[-?NS] Spin direction
character*1 u_spin2 (nr__5) ! [) ] Possible parenthesis
character*1 n_pole2 (nr__5) ! *[A-C ] Note on pole position
character*1 rot (nr__5) ! *Sense of rotation
real*8 per_1 (nr__5) ! (d) ?Period of rotation
character*1 n_per (nr__5) ! Mean of discordant periods = *
character*1 u_per (nr__5) ! [* ] indicates discordant periods
character*1 l_a_b (nr__5) ! [(< ] Possible parenthesis or limit
real*4 a_b (nr__5) ! ?Model axial ratio
character*1 u_a_b (nr__5) ! [) ] Possible parenthesis
real*4 b_c (nr__5) ! ?Model axial ratio
character*1 u_b_c (nr__5) ! [) ] Possible parenthesis
character*1 n_shape (nr__5) ! *[*?E-J ] Note on shape
character*7 refs_1 (nr__5) ! Reference codes
*Note on n_Planet:
* A = Amplitudes of lightcurves;
* D = Individual data-points of photometric lightcurves;
* E = Epochs (e.g., times of lightcurve extreme);
* F = Fourier coefficients of photometric lightcurves;
* I = Infrared pre- and post-opposition differences;
* M = Magnitudes (usually at maximum light);
* O = 0ccultation observations;
* P = Infrared polarimetry;
* R = Radar observations;
* S = Speckle interferometry;
* V = Visual position angles;
* Z = Zero and nonzero amplitude apparitions implying pole-on view in
* former case.
*Note on ELpole, EBpole, EL2pole, EB2pole:
* Pole Coordinates are given in the ecliptical reference frame of
* equinox 1950.
*Note on spin and spin2:
* The sense of rotation is termed indeterminate (1) when the spin
* direction is known, but the accuracy of the pole coordinates is
* insufficient to give an unambiguous sense of rotation, or the pole is
* so close to the ecliptic and/or orbital plane that forced precession
* will cause the sense of rotation to alternate with time.
*Note on n_pole2:
* A = concentric ring region
* B = aspect circle
* C = at least one pole position rejected
*Note on rot:
* - = Sense of rotation not determined;
* P = Prograde rotation;
* R = Retrograde rotation;
* I = Indeterminate case.
* M = Moving
*Note on n_shape:
* * = Mean value of two significantly different solutions (or discordant
* shape?)
* ? = One or both shapes questionable
* E = flat region
* F = see ref.
* G = "kettle"
* H = rnd. cyl.
* J = complex
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pole_ref.dat' ! References for pole.dat
integer*4 nr__6
parameter (nr__6=140) ! Number of records
character*80 ar__6 ! Full-size record
character*6 Ref_2 (nr__6) ! *Reference identification
character*72 Text_2 (nr__6) ! Reference
*Note on Ref:
* The reference number is repeated when several lines are necessary.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'taxonomy.dat' ! Taxonomic classifications
integer*4 nr__7
parameter (nr__7=982) ! Number of records
character*46 ar__7 ! Full-size record
integer*4 Planet_4 (nr__7) ! ?Asteroid number
character*13 name_2 (nr__7) ! Asteroid name
character*4 ThCl (nr__7) ! *Tholen class
character*2 n_ThCl (nr__7) ! Uncertainty note on ThCl
character*2 BarCl (nr__7) ! Barucci class
character*3 group (nr__7) ! *Asteroid group
integer*4 n_Planet_1 (nr__7) ! *? Note on asteroid
*Note on ThCl:
* The following notation appears in the classifications:
* U suffix indicating an unusual spectrum; falls far from cluster center
* : suffix indicating noisy data
* :: suffix indicating very noisy data
* --- indicates data that are too noisy to permit classification
* (essentially all types would be allowed)
*Note on group:
* Due to popular demand, orbital group designations have been included in this
* table. The 2- or 3-letter abbreviations stand for the following groups:
* ATE Aten
* APO Apollo
* AMO Amor
* MC Mars crosser
* HUN Hungaria
* PHO Phocaea
* GRI Griqua
* CYB Cybele
* HIL Hilda
* TRO Trojan
*Note on n_Planet:
* 1: 13 Egeria. Tholen (1984) listed classification as CG.
* C eliminated on the basis of 24-color data.
* 2: 192 Nausikaa. Tholen (1984) listed classification as RS.
* R eliminated on the basis of 24-color data.
* 3: 344 Desiderata. Tholen (1984) listed classification as CSU.
* SU eliminated on the basis of 24-color data.
* 4: 515 Athalia. Tholen (1984) listed classification as U.
* Changed to I as explained in the introduction to the printed article
* on the basis of 24-color data.
* 5: 3200 Phaethon. From Tholen (1985IAUC.4034....1H)
* 6: 3551 1983 RD, 3352 1983 SA, and 1984 BC. Unpublished data of
* Tholen cited in Hartmann et al. (1987Icar...69...33H).
* 7: 3908 1980 PA and 4055 1985 DO2. From Tholen et al. (1988)
* 8: 1975 U2. Not a proper provisional designation. Data taken from TRIAD
* UBV table as published by Bowell et al. (1979). In turn, they took the
* data from unpublished observations by Tedesco. It is not known what this
* object really is.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'mag.dat' ! Magnitudes, colors, albedos, and diameters
integer*4 nr__8
parameter (nr__8=3317) ! Number of records
character*63 ar__8 ! Full-size record
integer*4 Planet_5 (nr__8) ! *Minor planet number
real*4 Mag (nr__8) ! (mag) *Mean visual absolute magnitude
real*4 slope (nr__8) ! (mag) *Slope parameter
integer*4 q_slope (nr__8) ! *Quality of slope (3 highest)
real*4 U_B (nr__8) ! (mag) *?U-B color
integer*4 q_U_B (nr__8) ! *?Quality of U-B
real*4 B_V (nr__8) ! (mag) *?B-V color
integer*4 q_B_V (nr__8) ! *?Quality of B-V
integer*4 r_B_V (nr__8) ! *?Source of colors
real*4 albedo (nr__8) ! *?Visual geometric albedo
real*4 e_albedo (nr__8) ! *?Standard deviation of albedo
real*4 diam (nr__8) ! (km) *?Diameter
real*4 e_diam (nr__8) ! (km) *?Standard deviation of diameter
integer*4 sight (nr__8) ! *?Number of sightings
*Note on Planet:
* Asteroids numbered through 3318 are included. There is one record for
* each asteroid.
*Note on Mag:
* The visual absolute magnitude from the IRAS ground-based input data
* set. These differ from the blue absolute magnitudes (HB) given in
* final data product No. 4. In many cases they also differ from those
* published in the 1988 Russian ephemeris.
*Note on slope:
* The slope parameter from the IRAS ground-based input data set. In many
* cases these differ from those published in the 1988 Russian ephemeris.
*Note on q_slope:
* The quality code for the absolute magnitude and slope parameter. 3 is
* best and 1 worst.
*Note on U-B, q_U-B, B-V, and q_B-V:
* The U-B color index from IRAS Asteroid and Comet Survey final data
* product 13.
* The quality codes for the UBV colors were obtained as follows:
* Quality Code For sigma
* 0 > 0.05 magnitude
* 1 <= 0.05 but > 0.03 magnitude
* 2 <= 0.03 but > 0.02 magnitude, and
* 3 <= 0.02 magnitude
* 4 <= 0.02 magnitude and good agreement between
* results from references 1 and 2.
*Note on r_B-V:
* The reference for the UBV color(s):
* 1 = Bowell et al. (1979) In Asteroids (T. Gehrels, ed.), pp. 1108-1129.
* ("TRIAD")
* 2 = Zellner et al. (1985Icar...61..355Z) Icarus 61, 355-416. ("ECAS")
* 3 = Weighted mean of values from references 1 and 2.
*Note on albedo:
* Visual geometric albedo from IRAS Asteroid and Comet Survey final data
* product No. 4
*Note on e_albedo:
* One sigma formal uncertainty in the visual geometric albedo from IRAS
* Asteroid and Comet Survey final data product No. 4
*Note on diam and e_diam:
* Diameter (in km) and one sigma uncertainty from IRAS Asteroid and
* Comet Survey final data product No. 4
*Note on sight:
* Number of IRAS observations used in computing the albedo and diameter
* (from IRAS Asteroid and Comet Survey final data product No. 4)
C=============================================================================
C Loading file 'discover.dat' ! Names and discovery circumstances
C Format for file interpretation
1 format(
+ I4,1X,A17,1X,A9,1X,I4,1X,A3,1X,I2,2X,A29,1X,A21,1X,I3,1X,A23)
C Effective file loading
open(unit=1,status='old',file=
+'discover.dat')
write(6,*) '....Loading file: discover.dat'
do i__=1,6213
read(1,'(A125)')ar__
read(ar__,1)
+ Planet(i__),name(i__),pre_name(i__),disc_Y(i__),disc_M(i__),
+ disc_D(i__),dis_name(i__),dis_pl(i__),note(i__),MPC(i__)
if(ar__(99:101) .EQ. '') note(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Planet(i__),name(i__),pre_name(i__),disc_Y(i__),disc_M(i__),
+ disc_D(i__),dis_name(i__),dis_pl(i__),note(i__),MPC(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'dis_note.dat' ! The notes for discover.dat
C Format for file interpretation
2 format(I3,1X,A76)
C Effective file loading
open(unit=1,status='old',file=
+'dis_note.dat')
write(6,*) '....Loading file: dis_note.dat'
do i__=1,315
read(1,'(A80)')ar__1
read(ar__1,2)Note_1(i__),Text(i__)
c ..............Just test output...........
write(6,2)Note_1(i__),Text(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'lightcrv.dat' ! Lightcurve parameters
C Format for file interpretation
3 format(I4,1X,A16,A1,F8.3,1X,A1,F4.2,A1,F4.2,A1,1X,I1,A2,3X,A69)
C Effective file loading
open(unit=1,status='old',file=
+'lightcrv.dat')
write(6,*) '....Loading file: lightcrv.dat'
do i__=1,528
read(1,'(A118)')ar__2
read(ar__2,3)
+ Planet_1(i__),name_1(i__),l_per(i__),per(i__),l_mag1(i__),
+ mag1(i__),l_mag2(i__),mag2(i__),u_mag2(i__),n_mag(i__),
+ rem(i__),refs(i__)
if(ar__2(1:4) .EQ. '') Planet_1(i__) = iNULL__
if(ar__2(23:30) .EQ. '') per(i__) = rNULL__
if(ar__2(33:36) .EQ. '') mag1(i__) = rNULL__
if(ar__2(38:41) .EQ. '') mag2(i__) = rNULL__
if(ar__2(44:44) .EQ. '') n_mag(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Planet_1(i__),name_1(i__),l_per(i__),per(i__),l_mag1(i__),
+ mag1(i__),l_mag2(i__),mag2(i__),u_mag2(i__),n_mag(i__),
+ rem(i__),refs(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ltcv_ref.dat' ! References file for lightcrv.dat
C Format for file interpretation
4 format(I3,2X,A75)
C Effective file loading
open(unit=1,status='old',file=
+'ltcv_ref.dat')
write(6,*) '....Loading file: ltcv_ref.dat'
do i__=1,497
read(1,'(A80)')ar__3
read(ar__3,4)Ref(i__),Text_1(i__)
c ..............Just test output...........
write(6,4)Ref(i__),Text_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'proper.dat' ! Proper elements and family identities
C Format for file interpretation
5 format(
+ I4,1X,F6.3,F5.3,F5.3,F6.1,F6.1,F6.1,F7.1,1X,A3,F6.3,F6.3,1X,
+ A1,1X,A3,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'proper.dat')
write(6,*) '....Loading file: proper.dat'
do i__=1,1796
read(1,'(A70)')ar__4
read(ar__4,5)
+ Planet_2(i__),a(i__),e(i__),sin_i(i__),wbar(i__),anode(i__),
+ dwbar(i__),dnode(i__),res(i__),dMars(i__),dJup(i__),
+ q_Planet(i__),family(i__),Ref_1(i__)
if(ar__4(12:16) .EQ. '') e(i__) = rNULL__
if(ar__4(17:21) .EQ. '') sin_i(i__) = rNULL__
if(ar__4(22:27) .EQ. '') wbar(i__) = rNULL__
if(ar__4(28:33) .EQ. '') anode(i__) = rNULL__
if(ar__4(34:39) .EQ. '') dwbar(i__) = rNULL__
if(ar__4(40:46) .EQ. '') dnode(i__) = rNULL__
if(ar__4(51:56) .EQ. '') dMars(i__) = rNULL__
if(ar__4(57:62) .EQ. '') dJup(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Planet_2(i__),a(i__),e(i__),sin_i(i__),wbar(i__),anode(i__),
+ dwbar(i__),dnode(i__),res(i__),dMars(i__),dJup(i__),
+ q_Planet(i__),family(i__),Ref_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pole.dat' ! Pole determinations
C Format for file interpretation
6 format(
+ I4,1X,A4,2X,I3,A1,I3,A1,A1,A1,2X,I3,1X,I3,1X,A1,A1,1X,A1,1X,
+ A1,4X,F10.8,A1,A1,A1,F4.2,A1,1X,F4.2,A1,1X,A1,A7)
C Effective file loading
open(unit=1,status='old',file=
+'pole.dat')
write(6,*) '....Loading file: pole.dat'
do i__=1,248
read(1,'(A75)')ar__5
read(ar__5,6)
+ Planet_3(i__),n_Planet(i__),ELpole(i__),u_pole(i__),
+ EBpole(i__),n_pole(i__),spin(i__),u_spin(i__),EL2pole(i__),
+ EB2pole(i__),spin2(i__),u_spin2(i__),n_pole2(i__),rot(i__),
+ per_1(i__),n_per(i__),u_per(i__),l_a_b(i__),a_b(i__),
+ u_a_b(i__),b_c(i__),u_b_c(i__),n_shape(i__),refs_1(i__)
if(ar__5(12:14) .EQ. '') ELpole(i__) = iNULL__
if(ar__5(16:18) .EQ. '') EBpole(i__) = iNULL__
if(ar__5(24:26) .EQ. '') EL2pole(i__) = iNULL__
if(ar__5(28:30) .EQ. '') EB2pole(i__) = iNULL__
if(ar__5(42:51) .EQ. '') per_1(i__) = rNULL__
if(ar__5(55:58) .EQ. '') a_b(i__) = rNULL__
if(ar__5(61:64) .EQ. '') b_c(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ Planet_3(i__),n_Planet(i__),ELpole(i__),u_pole(i__),
+ EBpole(i__),n_pole(i__),spin(i__),u_spin(i__),EL2pole(i__),
+ EB2pole(i__),spin2(i__),u_spin2(i__),n_pole2(i__),rot(i__),
+ per_1(i__),n_per(i__),u_per(i__),l_a_b(i__),a_b(i__),
+ u_a_b(i__),b_c(i__),u_b_c(i__),n_shape(i__),refs_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pole_ref.dat' ! References for pole.dat
C Format for file interpretation
7 format(A6,2X,A72)
C Effective file loading
open(unit=1,status='old',file=
+'pole_ref.dat')
write(6,*) '....Loading file: pole_ref.dat'
do i__=1,140
read(1,'(A80)')ar__6
read(ar__6,7)Ref_2(i__),Text_2(i__)
c ..............Just test output...........
write(6,7)Ref_2(i__),Text_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'taxonomy.dat' ! Taxonomic classifications
C Format for file interpretation
8 format(I4,1X,A13,1X,A4,A2,1X,A2,1X,A3,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'taxonomy.dat')
write(6,*) '....Loading file: taxonomy.dat'
do i__=1,982
read(1,'(A46)')ar__7
read(ar__7,8)
+ Planet_4(i__),name_2(i__),ThCl(i__),n_ThCl(i__),BarCl(i__),
+ group(i__),n_Planet_1(i__)
if(ar__7(1:4) .EQ. '') Planet_4(i__) = iNULL__
if(ar__7(34:34) .EQ. '') n_Planet_1(i__) = iNULL__
c ..............Just test output...........
write(6,8)
+ Planet_4(i__),name_2(i__),ThCl(i__),n_ThCl(i__),BarCl(i__),
+ group(i__),n_Planet_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'mag.dat' ! Magnitudes, colors, albedos, and diameters
C Format for file interpretation
9 format(
+ I4,1X,F5.2,F6.3,1X,I1,1X,F4.2,1X,I1,1X,F4.2,1X,I1,1X,I1,2X,
+ F5.3,1X,F5.3,1X,F6.2,1X,F5.2,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'mag.dat')
write(6,*) '....Loading file: mag.dat'
do i__=1,3317
read(1,'(A63)')ar__8
read(ar__8,9)
+ Planet_5(i__),Mag(i__),slope(i__),q_slope(i__),U_B(i__),
+ q_U_B(i__),B_V(i__),q_B_V(i__),r_B_V(i__),albedo(i__),
+ e_albedo(i__),diam(i__),e_diam(i__),sight(i__)
if(ar__8(20:23) .EQ. '') U_B(i__) = rNULL__
if(ar__8(25:25) .EQ. '') q_U_B(i__) = iNULL__
if(ar__8(27:30) .EQ. '') B_V(i__) = rNULL__
if(ar__8(32:32) .EQ. '') q_B_V(i__) = iNULL__
if(ar__8(34:34) .EQ. '') r_B_V(i__) = iNULL__
if(ar__8(37:41) .EQ. '') albedo(i__) = rNULL__
if(ar__8(43:47) .EQ. '') e_albedo(i__) = rNULL__
if(ar__8(49:54) .EQ. '') diam(i__) = rNULL__
if(ar__8(56:60) .EQ. '') e_diam(i__) = rNULL__
if(ar__8(62:63) .EQ. '') sight(i__) = iNULL__
c ..............Just test output...........
write(6,9)
+ Planet_5(i__),Mag(i__),slope(i__),q_slope(i__),U_B(i__),
+ q_U_B(i__),B_V(i__),q_B_V(i__),r_B_V(i__),albedo(i__),
+ e_albedo(i__),diam(i__),e_diam(i__),sight(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end