Conversion of standardized ReadMe file for
file /./home/cats/J/A_A/287/990 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.
****Equinox Unknown in Position Column 'RAdeg' for: table3a
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-May-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. J/A+A/287/990 Meteor stream activity. (Jenniskens 1994)
*================================================================================
*Meteor stream activity. 1. The annual streams.
* Jenniskens P.
* <Astron. Astrophys. 287, 990 (1994)>
* =1994A&A...287..990J
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3a' ! Catalogue of meteor streams
integer*4 nr__
parameter (nr__=50) ! Number of records
character*77 ar__ ! Full-size record
character*4 Code (nr__) ! code name
character*22 Name (nr__) ! common name of the stream
integer*4 lo (nr__) ! (deg) solar longitude for which radiant position
* is valid
integer*4 RAdeg (nr__) ! (deg) right ascension of radiant
integer*4 DEdeg (nr__) ! (deg) declination of radiant
real*4 dRA (nr__) ! radiant drift in right ascension (degree
* of arc per degree of solar longitude)
real*4 dDE (nr__) ! radiant drift in declination
integer*4 Vinf (nr__) ! (km/s) geocentric entry velocity of meteoroids in the
* atmosphere, which includes the Earth's
* gravitational acceleration
character*1 n_Chi (nr__) ! [*] when *, Chi varies with solar longitude
real*4 Chi (nr__) ! magnitude distribution index
* (sporadic meteors Chi = 3.4)
integer*4 N (nr__) ! total number of meteors observed of this stream
integer*4 Neff (nr__) ! (h) total number of effective observing time
integer*4 Nobs (nr__) ! ? number of observers that contributed to
* data of this stream
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3b' ! Properties related to the activity curve of
the meteor streams
integer*4 nr__1
parameter (nr__1=50) ! Number of records
character*138 ar__1 ! Full-size record
character*4 Name_1 (nr__1) ! abbreviation name
real*4 lo_max (nr__1) ! (deg) ? solar longitude of peak position (from
* intersection point of slopes fitted to
* the data). Equinox 1950.0 (1)
character*1 u_lo_max (nr__1) ! [)] when ')', the values are for an assumed
* semi-major axis
real*4 e_lo_max (nr__1) ! (deg) ? rms uncertainty on lo(max)
character*1 l_ZHRmax (nr__1) ! limit flag on ZHRmax
real*4 ZHRmax (nr__1) ! (h-1) peak rate (Zenith Hourly Rates) (1), (2)
character*1 u_ZHRmax (nr__1) ! [)] when ')', the values are for an assumed
* semi-major axis
real*4 e_ZHRmax (nr__1) ! (h-1) ? rms uncertainty on ZHRmax
real*4 B (nr__1) ! (deg-1) ? slope of ascending and descending
* branches (1)
character*1 u_B (nr__1) ! [)] when ')', the values are for an assumed
* semi-major axis
character*1 n_B (nr__1) ! note on the slope B (4)
real*4 e_B (nr__1) ! (deg-1) ? rms uncertainty on the slope B
real*4 P (nr__1) ! (yr) ? mean orbital period of meteoroids. (3)
character*1 n_P (nr__1) ! [ :*i] : is an uncertainty flag on P
* * means that P is orbital period of
* parent comet
* i means infinity
real*4 q (nr__1) ! (AU) perihelion distance
integer*4 i (nr__1) ! (deg) inclination of the orbit with respect
* to the ecliptic
integer*4 Pi (nr__1) ! (deg) longitude of perihelion (sum of argument of
* perihelion (small omega) and the ascending
* node (large omega))
real*4 M_0 (nr__1) ! (g) ? mass of zero magnitude meteor
character*1 l_Rho (nr__1) ! limit flag on Rho
real*4 Rho (nr__1) ! (10-24g/cm3) density of matter in the peak of the
* meteoroid stream
character*1 u_Rho (nr__1) ! [)] when ')', the values are for an assumed
* semi-major axis
real*4 e_Rho (nr__1) ! (10-24g/cm3) ? rms uncertainty on the density Rho
character*1 l_M1yr (nr__1) ! limit flag on the mass M1yr
real*8 M1yr (nr__1) ! (10+14g) ? total mass in a flux tube of length 1
* year times heliocentric velocity
character*1 u_M1yr (nr__1) ! [)] when ')', the values are for an assumed
* semi-major axis
real*8 e_M1yr (nr__1) ! (10+14g) ? rms uncertainty on the mass M1yr
character*1 l_Mtot (nr__1) ! limit flag on the mass Mtot
real*4 Mtot (nr__1) ! (10+15g) ? total mass estimate (Mtot = M1yr x P)
character*1 n_Mtot (nr__1) ! when *, value of Mtot is for an assumed
* period P = 20 yr. These streams lack
* photographically determined orbits
real*4 e_Mtot (nr__1) ! (10+15g) ? rms uncertainty on the mass Mtot
*Note (1): ZHR = ZHRmax 10^( -B |lo - lo(max)|)
*Note (2): The peak rate is in number per hour that can be seen by a standard
* observing in perfectly clear sky (star limiting magnitude = 6.5) and
* a radiant position in the zenith.
*Note (3): The average of the mean orbital period of meteoroids is over 1/a,
* where a is the semi-major axis of photographically determined orbits
* in the literature.
*Note (4): a: for Boo, B = 1.8 (4)
* b: for Gem, B = 0.39(4)/0.72(10)
* Values between brackets are for assumed semi-major axis
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3c' ! Results for activity profiles with a
background component
integer*4 nr__2
parameter (nr__2=6) ! Number of records
character*135 ar__2 ! Full-size record
character*4 Name_2 (nr__2) ! abbreviation name
real*4 lo_max_1 (nr__2) ! (deg) solar longitude of peak position 1950.0
real*4 ZHRmax_p (nr__2) ! (h-1) peak rate for the main peak
real*4 e_ZHRmax_p (nr__2) ! (h-1) rms uncertainty on ZHRmax
character*1 l_Bp (nr__2) ! [ ~] limit flag on the slope Bp
real*4 Bp (nr__2) ! (deg-1) ? slope of ascending and descending
* branches for the main peak
character*1 n_Bp (nr__2) ! note on the slope Bp (1)
real*4 e_Bp (nr__2) ! (deg-1) ? rms uncertainty on the slope Bp
character*1 l_Mtot_p (nr__2) ! [ ~] limit flag on the mass Mtot(p)
real*8 Mtot_p (nr__2) ! (10+15g) total mass estimate for the main peak
real*4 e_Mtot_p (nr__2) ! (10+15g) ? rms uncertainty on the mass Mtot(p)
real*4 ZHRmax_b (nr__2) ! (h-1) peak rate for the background
real*4 e_ZHRmax_b (nr__2) ! (h-1) rms uncertainty on ZHRmax(b)
real*4 Bb_ (nr__2) ! (deg-1) slope of ascending branch
real*4 e_Bb_ (nr__2) ! (deg-1) rms uncertainty on the slope Bb+
character*1 l_Bb_ (nr__2) ! [ ~>] limit flag on the slope Bd-
real*4 Bb__1 (nr__2) ! (deg-1) slope of descending branch
real*4 e_Bb__1 (nr__2) ! (deg-1) ? rms uncertainty on the slope Bb-
real*4 Mtot_b (nr__2) ! (10+15g) total mass estimate for the background
real*4 e_Mtot_b (nr__2) ! (10+15g) rms uncertainty on the mass Mtot(b)
real*4 Mb_Mp (nr__2) ! background to main peak mass ratio
real*4 e_Mb_Mp (nr__2) ! rms uncertainty on Mb/Mp
real*4 Mtot_1 (nr__2) ! (10+15g) total mass estimate
real*4 e_Mtot_1 (nr__2) ! (10+15g) rms uncertainty on the mass Mtot
*Note (1): a: For Gem, Bp = (0.59/0.81) +/-0.07
* Values between brackets are for assumed semi-major axis
C=============================================================================
C Loading file 'table3a' ! Catalogue of meteor streams
C Format for file interpretation
1 format(
+ A4,1X,A22,I3,1X,I3,1X,I3,1X,F4.1,1X,F4.1,1X,I2,A1,F3.1,1X,I5,
+ 1X,I4,7X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table3a')
write(6,*) '....Loading file: table3a'
do i__=1,50
read(1,'(A77)')ar__
read(ar__,1)
+ Code(i__),Name(i__),lo(i__),RAdeg(i__),DEdeg(i__),dRA(i__),
+ dDE(i__),Vinf(i__),n_Chi(i__),Chi(i__),N(i__),Neff(i__),
+ Nobs(i__)
if(ar__(74:75) .EQ. '') Nobs(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Code(i__),Name(i__),lo(i__),RAdeg(i__),DEdeg(i__),dRA(i__),
+ dDE(i__),Vinf(i__),n_Chi(i__),Chi(i__),N(i__),Neff(i__),
+ Nobs(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3b' ! Properties related to the activity curve of
* the meteor streams
C Format for file interpretation
2 format(
+ A4,1X,F6.2,A1,F4.2,2X,A1,F5.1,A1,3X,F4.1,2X,F5.3,A1,A1,2X,
+ F5.3,4X,F6.2,A1,F4.2,2X,I3,2X,I3,2X,F6.3,1X,A1,F6.3,A1,3X,
+ F5.3,2X,A1,F7.4,A1,3X,F7.4,2X,A1,F6.3,A1,3X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3b')
write(6,*) '....Loading file: table3b'
do i__=1,50
read(1,'(A138)')ar__1
read(ar__1,2)
+ Name_1(i__),lo_max(i__),u_lo_max(i__),e_lo_max(i__),
+ l_ZHRmax(i__),ZHRmax(i__),u_ZHRmax(i__),e_ZHRmax(i__),B(i__),
+ u_B(i__),n_B(i__),e_B(i__),P(i__),n_P(i__),q(i__),i(i__),
+ Pi(i__),M_0(i__),l_Rho(i__),Rho(i__),u_Rho(i__),e_Rho(i__),
+ l_M1yr(i__),M1yr(i__),u_M1yr(i__),e_M1yr(i__),l_Mtot(i__),
+ Mtot(i__),n_Mtot(i__),e_Mtot(i__)
if(ar__1(6:11) .EQ. '') lo_max(i__) = rNULL__
if(ar__1(13:16) .EQ. '') e_lo_max(i__) = rNULL__
if(ar__1(29:32) .EQ. '') e_ZHRmax(i__) = rNULL__
if(ar__1(35:39) .EQ. '') B(i__) = rNULL__
if(ar__1(44:48) .EQ. '') e_B(i__) = rNULL__
if(ar__1(53:58) .EQ. '') P(i__) = rNULL__
if(ar__1(76:81) .EQ. '') M_0(i__) = rNULL__
if(ar__1(94:98) .EQ. '') e_Rho(i__) = rNULL__
if(ar__1(102:108) .EQ. '') M1yr(i__) = rNULL__
if(ar__1(113:119) .EQ. '') e_M1yr(i__) = rNULL__
if(ar__1(123:128) .EQ. '') Mtot(i__) = rNULL__
if(ar__1(133:137) .EQ. '') e_Mtot(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),lo_max(i__),u_lo_max(i__),e_lo_max(i__),
+ l_ZHRmax(i__),ZHRmax(i__),u_ZHRmax(i__),e_ZHRmax(i__),B(i__),
+ u_B(i__),n_B(i__),e_B(i__),P(i__),n_P(i__),q(i__),i(i__),
+ Pi(i__),M_0(i__),l_Rho(i__),Rho(i__),u_Rho(i__),e_Rho(i__),
+ l_M1yr(i__),M1yr(i__),u_M1yr(i__),e_M1yr(i__),l_Mtot(i__),
+ Mtot(i__),n_Mtot(i__),e_Mtot(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3c' ! Results for activity profiles with a
* background component
C Format for file interpretation
3 format(
+ A4,1X,F6.2,2X,F5.1,1X,F4.1,1X,A1,F4.2,A1,7X,F4.2,1X,A1,F7.4,
+ 2X,F6.4,1X,F4.1,2X,F3.1,1X,F5.3,2X,F5.3,2X,A1,F5.3,2X,F5.3,3X,
+ F6.3,2X,F5.3,1X,F3.1,2X,F3.1,1X,F6.3,2X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3c')
write(6,*) '....Loading file: table3c'
do i__=1,6
read(1,'(A135)')ar__2
read(ar__2,3)
+ Name_2(i__),lo_max_1(i__),ZHRmax_p(i__),e_ZHRmax_p(i__),
+ l_Bp(i__),Bp(i__),n_Bp(i__),e_Bp(i__),l_Mtot_p(i__),
+ Mtot_p(i__),e_Mtot_p(i__),ZHRmax_b(i__),e_ZHRmax_b(i__),
+ Bb_(i__),e_Bb_(i__),l_Bb_(i__),Bb__1(i__),e_Bb__1(i__),
+ Mtot_b(i__),e_Mtot_b(i__),Mb_Mp(i__),e_Mb_Mp(i__),Mtot_1(i__),
+ e_Mtot_1(i__)
if(ar__2(26:29) .EQ. '') Bp(i__) = rNULL__
if(ar__2(38:41) .EQ. '') e_Bp(i__) = rNULL__
if(ar__2(53:58) .EQ. '') e_Mtot_p(i__) = rNULL__
if(ar__2(92:96) .EQ. '') e_Bb__1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),lo_max_1(i__),ZHRmax_p(i__),e_ZHRmax_p(i__),
+ l_Bp(i__),Bp(i__),n_Bp(i__),e_Bp(i__),l_Mtot_p(i__),
+ Mtot_p(i__),e_Mtot_p(i__),ZHRmax_b(i__),e_ZHRmax_b(i__),
+ Bb_(i__),e_Bb_(i__),l_Bb_(i__),Bb__1(i__),e_Bb__1(i__),
+ Mtot_b(i__),e_Mtot_b(i__),Mb_Mp(i__),e_Mb_Mp(i__),Mtot_1(i__),
+ e_Mtot_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end