Conversion of standardized ReadMe file for
file /./home/cats/J/A_A/287/990 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.
****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-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/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 ! code name
character*22 Name ! common name of the stream
integer*4 lo ! (deg) solar longitude for which radiant position
* is valid
integer*4 RAdeg ! (deg) right ascension of radiant
integer*4 DEdeg ! (deg) declination of radiant
real*4 dRA ! radiant drift in right ascension (degree
* of arc per degree of solar longitude)
real*4 dDE ! radiant drift in declination
integer*4 Vinf ! (km/s) geocentric entry velocity of meteoroids in the
* atmosphere, which includes the Earth's
* gravitational acceleration
character*1 n_Chi ! [*] when *, Chi varies with solar longitude
real*4 Chi ! magnitude distribution index
* (sporadic meteors Chi = 3.4)
integer*4 N ! total number of meteors observed of this stream
integer*4 Neff ! (h) total number of effective observing time
integer*4 Nobs ! ? 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 ! abbreviation name
real*4 lo_max ! (deg) ? solar longitude of peak position (from
* intersection point of slopes fitted to
* the data). Equinox 1950.0 (1)
character*1 u_lo_max ! [)] when ')', the values are for an assumed
* semi-major axis
real*4 e_lo_max ! (deg) ? rms uncertainty on lo(max)
character*1 l_ZHRmax ! limit flag on ZHRmax
real*4 ZHRmax ! (h-1) peak rate (Zenith Hourly Rates) (1), (2)
character*1 u_ZHRmax ! [)] when ')', the values are for an assumed
* semi-major axis
real*4 e_ZHRmax ! (h-1) ? rms uncertainty on ZHRmax
real*4 B ! (deg-1) ? slope of ascending and descending
* branches (1)
character*1 u_B ! [)] when ')', the values are for an assumed
* semi-major axis
character*1 n_B ! note on the slope B (4)
real*4 e_B ! (deg-1) ? rms uncertainty on the slope B
real*4 P ! (yr) ? mean orbital period of meteoroids. (3)
character*1 n_P ! [ :*i] : is an uncertainty flag on P
* * means that P is orbital period of
* parent comet
* i means infinity
real*4 q ! (AU) perihelion distance
integer*4 i ! (deg) inclination of the orbit with respect
* to the ecliptic
integer*4 Pi ! (deg) longitude of perihelion (sum of argument of
* perihelion (small omega) and the ascending
* node (large omega))
real*4 M_0 ! (g) ? mass of zero magnitude meteor
character*1 l_Rho ! limit flag on Rho
real*4 Rho ! (10-24g/cm3) density of matter in the peak of the
* meteoroid stream
character*1 u_Rho ! [)] when ')', the values are for an assumed
* semi-major axis
real*4 e_Rho ! (10-24g/cm3) ? rms uncertainty on the density Rho
character*1 l_M1yr ! limit flag on the mass M1yr
real*8 M1yr ! (10+14g) ? total mass in a flux tube of length 1
* year times heliocentric velocity
character*1 u_M1yr ! [)] when ')', the values are for an assumed
* semi-major axis
real*8 e_M1yr ! (10+14g) ? rms uncertainty on the mass M1yr
character*1 l_Mtot ! limit flag on the mass Mtot
real*4 Mtot ! (10+15g) ? total mass estimate (Mtot = M1yr x P)
character*1 n_Mtot ! when *, value of Mtot is for an assumed
* period P = 20 yr. These streams lack
* photographically determined orbits
real*4 e_Mtot ! (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 ! abbreviation name
real*4 lo_max_1 ! (deg) solar longitude of peak position 1950.0
real*4 ZHRmax_p ! (h-1) peak rate for the main peak
real*4 e_ZHRmax_p ! (h-1) rms uncertainty on ZHRmax
character*1 l_Bp ! [ ~] limit flag on the slope Bp
real*4 Bp ! (deg-1) ? slope of ascending and descending
* branches for the main peak
character*1 n_Bp ! note on the slope Bp (1)
real*4 e_Bp ! (deg-1) ? rms uncertainty on the slope Bp
character*1 l_Mtot_p ! [ ~] limit flag on the mass Mtot(p)
real*8 Mtot_p ! (10+15g) total mass estimate for the main peak
real*4 e_Mtot_p ! (10+15g) ? rms uncertainty on the mass Mtot(p)
real*4 ZHRmax_b ! (h-1) peak rate for the background
real*4 e_ZHRmax_b ! (h-1) rms uncertainty on ZHRmax(b)
real*4 Bb_ ! (deg-1) slope of ascending branch
real*4 e_Bb_ ! (deg-1) rms uncertainty on the slope Bb+
character*1 l_Bb_ ! [ ~>] limit flag on the slope Bd-
real*4 Bb__1 ! (deg-1) slope of descending branch
real*4 e_Bb__1 ! (deg-1) ? rms uncertainty on the slope Bb-
real*4 Mtot_b ! (10+15g) total mass estimate for the background
real*4 e_Mtot_b ! (10+15g) rms uncertainty on the mass Mtot(b)
real*4 Mb_Mp ! background to main peak mass ratio
real*4 e_Mb_Mp ! rms uncertainty on Mb/Mp
real*4 Mtot_1 ! (10+15g) total mass estimate
real*4 e_Mtot_1 ! (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,Name,lo,RAdeg,DEdeg,dRA,dDE,Vinf,n_Chi,Chi,N,Neff,Nobs
if(ar__(74:75) .EQ. '') Nobs = iNULL__
c ..............Just test output...........
write(6,1)
+ Code,Name,lo,RAdeg,DEdeg,dRA,dDE,Vinf,n_Chi,Chi,N,Neff,Nobs
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,lo_max,u_lo_max,e_lo_max,l_ZHRmax,ZHRmax,u_ZHRmax,
+ e_ZHRmax,B,u_B,n_B,e_B,P,n_P,q,i,Pi,M_0,l_Rho,Rho,u_Rho,e_Rho,
+ l_M1yr,M1yr,u_M1yr,e_M1yr,l_Mtot,Mtot,n_Mtot,e_Mtot
if(ar__1(6:11) .EQ. '') lo_max = rNULL__
if(ar__1(13:16) .EQ. '') e_lo_max = rNULL__
if(ar__1(29:32) .EQ. '') e_ZHRmax = rNULL__
if(ar__1(35:39) .EQ. '') B = rNULL__
if(ar__1(44:48) .EQ. '') e_B = rNULL__
if(ar__1(53:58) .EQ. '') P = rNULL__
if(ar__1(76:81) .EQ. '') M_0 = rNULL__
if(ar__1(94:98) .EQ. '') e_Rho = rNULL__
if(ar__1(102:108) .EQ. '') M1yr = rNULL__
if(ar__1(113:119) .EQ. '') e_M1yr = rNULL__
if(ar__1(123:128) .EQ. '') Mtot = rNULL__
if(ar__1(133:137) .EQ. '') e_Mtot = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,lo_max,u_lo_max,e_lo_max,l_ZHRmax,ZHRmax,u_ZHRmax,
+ e_ZHRmax,B,u_B,n_B,e_B,P,n_P,q,i,Pi,M_0,l_Rho,Rho,u_Rho,e_Rho,
+ l_M1yr,M1yr,u_M1yr,e_M1yr,l_Mtot,Mtot,n_Mtot,e_Mtot
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,lo_max_1,ZHRmax_p,e_ZHRmax_p,l_Bp,Bp,n_Bp,e_Bp,
+ l_Mtot_p,Mtot_p,e_Mtot_p,ZHRmax_b,e_ZHRmax_b,Bb_,e_Bb_,l_Bb_,
+ Bb__1,e_Bb__1,Mtot_b,e_Mtot_b,Mb_Mp,e_Mb_Mp,Mtot_1,e_Mtot_1
if(ar__2(26:29) .EQ. '') Bp = rNULL__
if(ar__2(38:41) .EQ. '') e_Bp = rNULL__
if(ar__2(53:58) .EQ. '') e_Mtot_p = rNULL__
if(ar__2(92:96) .EQ. '') e_Bb__1 = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,lo_max_1,ZHRmax_p,e_ZHRmax_p,l_Bp,Bp,n_Bp,e_Bp,
+ l_Mtot_p,Mtot_p,e_Mtot_p,ZHRmax_b,e_ZHRmax_b,Bb_,e_Bb_,l_Bb_,
+ Bb__1,e_Bb__1,Mtot_b,e_Mtot_b,Mb_Mp,e_Mb_Mp,Mtot_1,e_Mtot_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end