Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/567/A126 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.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-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/567/A126 Warsaw Catalogue of cometary orbits (Krolikowska, 2014)
*================================================================================
*Warsaw Catalogue of cometary orbits: 119 near-parabolic comets.
* Krolikowska M.
* <Astron. Astrophys. 567, A126 (2014)>
* =2014A&A...567A.126K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Description of observational material and the new
orbital quality assessment for 119 near-parabolic
comets
integer*4 nr__
parameter (nr__=135) ! Number of records
character*184 ar__ ! Full-size record
character*2 Sample ! [A1 A2 B] Sample name (G1)
character*12 Desig ! Comet's designation
character*24 Cname ! Comet's name
real*4 qosc ! (AU) Occulating perihelion distance
integer*4 Time ! ("YYYYMMDD") Perihelion time
integer*4 arc1 ! ("YYYYMMDD") Date of the first observation
integer*4 arc2 ! ("YYYYMMDD") Date of the last observation
integer*4 Nobs ! [7/5337] number of observations
* taken into account
real*4 arcy ! (yr) interval of data (arc in years)
real*4 dh1 ! (AU) heliocentric distance of the first observation
real*4 dh2 ! (AU) heliocentric distance of the last observation
character*10 datat ! Data type in the sense of time-distribution
* around perihelion: pre, post, full
character*8 Model ! Model type (GR=pure gravitational,
* NG=non-gravitational)
real*4 Q_GR ! ? Accuracy of gravitational orbit
real*4 Q_NG ! ? Accuracy of non-gravitational orbit when
* determinable
character*3 Qnew ! New quality orbit assessment (1a+ is best)
real*4 RMS ! (arcsec) Root-mean-square error
integer*4 Nres ! [14/10589] Number of residuals taken for orbit
* determination for a given model
character*7 Ref ! References to the solution (see the original
* article for the remaining symbols) (1)
*Note (1): References to the solution are:
* P1 = Krolikowska and Dybczynski, 2010MNRAS.404.1886K
* P2 = Dybczynski and Krolikowska, 2011MNRAS.416...51D
* P3 = Krolikowska, Dybczynski and Sitarski, 2012A&A...544A.119K
* P4 = Krolikowska and Dybczynski, 2013MNRAS.435..440K
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb.dat' ! Heliocentric orbital elements of osculating
orbits for comets from Samples A1, A2 and B
integer*4 nr__1
parameter (nr__1=134) ! Number of records
character*198 ar__1 ! Full-size record
character*2 Sample_1 ! [A1 A2 B] Sample name (G1)
character*1 Com ! [a-l] Comments (1)
character*12 Desig_1 ! Comet's designation
integer*4 Epoch ! ("YYYYMMDD") Epoch, year
integer*4 T_Y ! (yr) Perihelion time, year
integer*4 T_M ! ("month") Perihelion time, month
real*8 T_D ! (d) Perihelion time, day
real*8 q ! (AU) Perihelion distance
real*8 e ! [0.96/1.06] Eccentricity
real*8 ARG ! (deg) [0/360] Argument of perihelion {omega} (J2000)
real*8 LONG ! (deg) [0/360] Longitude of the ascending node {Omega}
* (J2000)
real*8 i ! (deg) Inclination of the orbit in degrees (J2000)
real*8 v1_a ! (10-6/AU) Inverse of semimajor axis
real*8 e_T_D ! (d) Uncertainty of Perihelion time in days
real*8 e_q ! (AU) Uncertainty of perihelion distance
real*8 e_e ! Uncertainty of eccentricity
real*8 e_ARG ! (deg) Uncertainty of argument of perihelion
real*8 e_LONG ! (deg) Uncertainty of longitude of the ascending node
real*8 e_i ! (deg) Uncertainty of inclination
real*8 e_1_a ! (10-6/AU) Uncertainty of inverse of semimajor axis
*Note (1): Comments as follows:
* a = Osculating barycentric orbits based on the entire data sets (tableb)
* or Original barycentric orbits based on the entire data sets (tablec)
* b = Osculating orbits based on DIST or pre-perihelion data sets
* c = Osculating orbits based on post-perihelion data sets
* d = Osculating orbits from the entire data sets
* e = Osculating orbits from pre-perihelion data sets for investigating the
* past dynamical evolution or Original barycentric orbits based on
* pre-perihelion data sets
* f = Osculating orbits from post-perihelion data sets for investigating the
* future dynamical evolution
* g = Osculating orbits based on the entire data sets
* h = Original orbits based on the osculating orbits determined from the
* entire data sets
* i = Future barycentric orbits based on the entire data sets
* j = Future barycentric orbits based on DIST or post-perihelion data sets
* k = Future barycentric orbits based on post-perihelion data sets
* l = Future orbits based on the osculating orbits determined from the
* entire data sets
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablec.dat' ! Barycentric orbital elements of original orbits,
i.e. before entering the planetary zone for
near-parabolic comets of Samples A1, A2 and B
integer*4 nr__2
parameter (nr__2=128) ! Number of records
character*198 ar__2 ! Full-size record
character*2 Sample_2 ! [A1 A2 B] Sample name (G1)
character*1 Com_1 ! [a-l] Comments (1)
character*12 Desig_2 ! Comet's designation
integer*4 Epoch_1 ! ("YYYYMMDD") Epoch, year
integer*4 T_Y_1 ! (yr) Perihelion time, year
integer*4 T_M_1 ! ("month") Perihelion time, month
real*8 T_D_1 ! (d) Perihelion time, day
real*8 q_1 ! (AU) Perihelion distance
real*8 e_1 ! [0.96/1.06] Eccentricity
real*8 ARG_1 ! (deg) [0/360] Argument of perihelion {omega} (J2000)
real*8 LONG_1 ! (deg) [0/360] Longitude of the ascending node {Omega}
* (J2000)
real*8 i_1 ! (deg) Inclination of the orbit in degrees (J2000)
real*8 v1_a_1 ! (10-6/AU) Inverse of semimajor axis
real*8 e_T_D_1 ! (d) Uncertainty of Perihelion time in days
real*8 e_q_1 ! (AU) Uncertainty of perihelion distance
real*8 e_e_1 ! Uncertainty of eccentricity
real*8 e_ARG_1 ! (deg) Uncertainty of argument of perihelion
real*8 e_LONG_1 ! (deg) Uncertainty of longitude of the ascending node
real*8 e_i_1 ! (deg) Uncertainty of inclination
real*8 e_1_a_1 ! (10-6/AU) Uncertainty of inverse of semimajor axis
*Note (1): Comments as follows:
* a = Osculating barycentric orbits based on the entire data sets (tableb)
* or Original barycentric orbits based on the entire data sets (tablec)
* b = Osculating orbits based on DIST or pre-perihelion data sets
* c = Osculating orbits based on post-perihelion data sets
* d = Osculating orbits from the entire data sets
* e = Osculating orbits from pre-perihelion data sets for investigating the
* past dynamical evolution or Original barycentric orbits based on
* pre-perihelion data sets
* f = Osculating orbits from post-perihelion data sets for investigating the
* future dynamical evolution
* g = Osculating orbits based on the entire data sets
* h = Original orbits based on the osculating orbits determined from the
* entire data sets
* i = Future barycentric orbits based on the entire data sets
* j = Future barycentric orbits based on DIST or post-perihelion data sets
* k = Future barycentric orbits based on post-perihelion data sets
* l = Future orbits based on the osculating orbits determined from the
* entire data sets
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled.dat' ! Barycentric orbital elements of future orbits,
i.e. after leaving the planetary zone for
near-parabolic comets of Samples A1, A2 and B
integer*4 nr__3
parameter (nr__3=126) ! Number of records
character*198 ar__3 ! Full-size record
character*2 Sample_3 ! [A1 A2 B] Sample name (G1)
character*1 Com_2 ! [a-l] Comments (1)
character*12 Desig_3 ! Comet's designation
integer*4 Epoch_2 ! ("YYYYMMDD") Epoch, year
integer*4 T_Y_2 ! (yr) Perihelion time, year
integer*4 T_M_2 ! ("month") Perihelion time, month
real*8 T_D_2 ! (d) Perihelion time, day
real*8 q_2 ! (AU) Perihelion distance
real*8 e_2 ! [0.96/1.06] Eccentricity
real*8 ARG_2 ! (deg) [0/360] Argument of perihelion {omega} (J2000)
real*8 LONG_2 ! (deg) [0/360] Longitude of the ascending node {Omega}
* (J2000)
real*8 i_2 ! (deg) Inclination of the orbit in degrees (J2000)
real*8 v1_a_2 ! (10-6/AU) Inverse of semimajor axis
real*8 e_T_D_2 ! (d) Uncertainty of Perihelion time in days
real*8 e_q_2 ! (AU) Uncertainty of perihelion distance
real*8 e_e_2 ! Uncertainty of eccentricity
real*8 e_ARG_2 ! (deg) Uncertainty of argument of perihelion
real*8 e_LONG_2 ! (deg) Uncertainty of longitude of the ascending node
real*8 e_i_2 ! (deg) Uncertainty of inclination
real*8 e_1_a_2 ! (10-6/AU) Uncertainty of inverse of semimajor axis
*Note (1): Comments as follows:
* a = Osculating barycentric orbits based on the entire data sets (tableb)
* or Original barycentric orbits based on the entire data sets (tablec)
* b = Osculating orbits based on DIST or pre-perihelion data sets
* c = Osculating orbits based on post-perihelion data sets
* d = Osculating orbits from the entire data sets
* e = Osculating orbits from pre-perihelion data sets for investigating the
* past dynamical evolution or Original barycentric orbits based on
* pre-perihelion data sets
* f = Osculating orbits from post-perihelion data sets for investigating the
* future dynamical evolution
* g = Osculating orbits based on the entire data sets
* h = Original orbits based on the osculating orbits determined from the
* entire data sets
* i = Future barycentric orbits based on the entire data sets
* j = Future barycentric orbits based on DIST or post-perihelion data sets
* k = Future barycentric orbits based on post-perihelion data sets
* l = Future orbits based on the osculating orbits determined from the
* entire data sets
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb4.dat' ! NG parameters for comets with NG heliocentric
orbital elements of osculating orbits in tableb.dat
integer*4 nr__4
parameter (nr__4=54) ! Number of records
character*115 ar__4 ! Full-size record
character*2 Sample_4 ! [A1 A2 B] Sample name (G1)
character*12 Desig_4 ! Comet's designation
real*8 A1 ! (10-8AU/d2) Radial non-gravitational (NG) parameter
real*8 e_A1 ! (10-8AU/d2) Uncertainy of radial NG parameter
real*8 A2 ! (10-8AU/d2) Transverse non-gravitational (NG) parameter
real*8 e_A2 ! (10-8AU/d2) Uncertainy of transverse NG parameter
real*8 A3 ! (10-8AU/d2) Normal non-gravitational (NG) parameter
real*8 e_A3 ! (10-8AU/d2) Uncertainy of normal NG parameter
real*8 tau ! (d) Time shift of maximum of g(r) relative to
* perihelion time (2)
real*8 e_tau ! (d) ? Uncertainty of perihelion time (3)
character*4 Model_1 ! Type of model (STD, PRE, POST, or DIST)
*Note (2): negative value means that maximum of g(r) is before perihelion passage
*Note (3): for C/1959 Y1 tau was assumed
C=============================================================================
C Loading file 'tablea1.dat' ! Description of observational material and the new
* orbital quality assessment for 119 near-parabolic
* comets
C Format for file interpretation
1 format(
+ A2,1X,A12,1X,A24,F6.3,6X,I8,6X,I8,2X,I8,4X,I5,4X,F5.2,4X,F5.2,
+ 1X,F5.2,5X,A10,A8,F3.1,5X,F3.1,5X,A3,4X,F4.2,I6,4X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,135
read(1,'(A184)')ar__
read(ar__,1)
+ Sample,Desig,Cname,qosc,Time,arc1,arc2,Nobs,arcy,dh1,dh2,
+ datat,Model,Q_GR,Q_NG,Qnew,RMS,Nres,Ref
if(ar__(141:143) .EQ. '') Q_GR = rNULL__
if(ar__(149:151) .EQ. '') Q_NG = rNULL__
c ..............Just test output...........
write(6,1)
+ Sample,Desig,Cname,qosc,Time,arc1,arc2,Nobs,arcy,dh1,dh2,
+ datat,Model,Q_GR,Q_NG,Qnew,RMS,Nres,Ref
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb.dat' ! Heliocentric orbital elements of osculating
* orbits for comets from Samples A1, A2 and B
C Format for file interpretation
2 format(
+ A2,1X,A1,1X,A12,I8,2X,I4,I2,F9.6,F14.8,F14.8,F12.6,F12.6,
+ F12.6,F9.2,F10.6,F14.8,F14.8,F12.6,F12.6,F12.6,F9.2)
C Effective file loading
open(unit=1,status='old',file=
+'tableb.dat')
write(6,*) '....Loading file: tableb.dat'
do i__=1,134
read(1,'(A198)')ar__1
read(ar__1,2)
+ Sample_1,Com,Desig_1,Epoch,T_Y,T_M,T_D,q,e,ARG,LONG,i,v1_a,
+ e_T_D,e_q,e_e,e_ARG,e_LONG,e_i,e_1_a
c ..............Just test output...........
write(6,2)
+ Sample_1,Com,Desig_1,Epoch,T_Y,T_M,T_D,q,e,ARG,LONG,i,v1_a,
+ e_T_D,e_q,e_e,e_ARG,e_LONG,e_i,e_1_a
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablec.dat' ! Barycentric orbital elements of original orbits,
* i.e. before entering the planetary zone for
* near-parabolic comets of Samples A1, A2 and B
C Format for file interpretation
3 format(
+ A2,1X,A1,1X,A12,I8,2X,I4,I2,F9.6,F14.8,F14.8,F12.6,F12.6,
+ F12.6,F9.2,F10.6,F14.8,F14.8,F12.6,F12.6,F12.6,F9.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablec.dat')
write(6,*) '....Loading file: tablec.dat'
do i__=1,128
read(1,'(A198)')ar__2
read(ar__2,3)
+ Sample_2,Com_1,Desig_2,Epoch_1,T_Y_1,T_M_1,T_D_1,q_1,e_1,
+ ARG_1,LONG_1,i_1,v1_a_1,e_T_D_1,e_q_1,e_e_1,e_ARG_1,e_LONG_1,
+ e_i_1,e_1_a_1
c ..............Just test output...........
write(6,3)
+ Sample_2,Com_1,Desig_2,Epoch_1,T_Y_1,T_M_1,T_D_1,q_1,e_1,
+ ARG_1,LONG_1,i_1,v1_a_1,e_T_D_1,e_q_1,e_e_1,e_ARG_1,e_LONG_1,
+ e_i_1,e_1_a_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled.dat' ! Barycentric orbital elements of future orbits,
* i.e. after leaving the planetary zone for
* near-parabolic comets of Samples A1, A2 and B
C Format for file interpretation
4 format(
+ A2,1X,A1,1X,A12,I8,2X,I4,I2,F9.6,F14.8,F14.8,F12.6,F12.6,
+ F12.6,F9.2,F10.6,F14.8,F14.8,F12.6,F12.6,F12.6,F9.2)
C Effective file loading
open(unit=1,status='old',file=
+'tabled.dat')
write(6,*) '....Loading file: tabled.dat'
do i__=1,126
read(1,'(A198)')ar__3
read(ar__3,4)
+ Sample_3,Com_2,Desig_3,Epoch_2,T_Y_2,T_M_2,T_D_2,q_2,e_2,
+ ARG_2,LONG_2,i_2,v1_a_2,e_T_D_2,e_q_2,e_e_2,e_ARG_2,e_LONG_2,
+ e_i_2,e_1_a_2
c ..............Just test output...........
write(6,4)
+ Sample_3,Com_2,Desig_3,Epoch_2,T_Y_2,T_M_2,T_D_2,q_2,e_2,
+ ARG_2,LONG_2,i_2,v1_a_2,e_T_D_2,e_q_2,e_e_2,e_ARG_2,e_LONG_2,
+ e_i_2,e_1_a_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb4.dat' ! NG parameters for comets with NG heliocentric
* orbital elements of osculating orbits in tableb.dat
C Format for file interpretation
5 format(
+ A2,1X,A12,2X,F12.5,F12.5,F11.5,1X,F11.5,2X,F11.5,F9.5,3X,F7.2,
+ F7.2,8X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'tableb4.dat')
write(6,*) '....Loading file: tableb4.dat'
do i__=1,54
read(1,'(A115)')ar__4
read(ar__4,5)
+ Sample_4,Desig_4,A1,e_A1,A2,e_A2,A3,e_A3,tau,e_tau,Model_1
if(ar__4(97:103) .EQ. '') e_tau = rNULL__
c ..............Just test output...........
write(6,5)
+ Sample_4,Desig_4,A1,e_A1,A2,e_A2,A3,e_A3,tau,e_tau,Model_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end