Conversion of standardized ReadMe file for
catalog 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.5, on 2013-May-25
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/A+A/502/695 Observations of HD 80606 planetary system (Pont+, 2009)
*================================================================================
*Spin-orbit misalignment in the HD 80606 planetary system.
* Pont F., Hebrard G., Irwin J.M., Bouchy F., Moutou C., Ehrenreich D.,
* Guillot T., Aigrain S., Bonfils X., Berta Z., Boisse I., Burke C.,
* Charbonneau D., Delfosse X., Desort M., Eggenberger A., Forveille T.,
* Lagrange A.-M., Lovis C., Nutzman P., Pepe F., Perrier C., Queloz D.,
* Santos N.C., Segransan D., Udry S., Vidal-Madjar A.
* <Astron. Astrophys., 502, 695-703 (2009)>
* =2009A&A...502..695P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'lc.dat' ! Light curve data for HD 80606 from OHP
integer*4 nr__
parameter (nr__=563) ! Number of records
character*36 ar__ ! Full-size record
real*8 HJD (nr__) ! (d) Heliocentric Julian date of observation
real*8 Flux (nr__) ! Normalized flux from 715nm long-pass filter
real*8 e_Flux (nr__) ! rms uncertainty on Flux
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rv.dat' ! Radial velocity data for HD 80606 with SOPHIE
integer*4 nr__1
parameter (nr__1=48) ! Number of records
character*29 ar__1 ! Full-size record
real*8 HJD_1 (nr__1) ! (d) Heliocentric Julian date of observation
real*4 RV (nr__1) ! (km/s) Radial velocity
real*4 e_RV (nr__1) ! (km/s) rms uncertainty on RV
C=============================================================================
C Loading file 'lc.dat' ! Light curve data for HD 80606 from OHP
C Format for file interpretation
1 format(F13.5,2X,F8.6,2X,F11.9)
C Effective file loading
open(unit=1,file='lc.dat', status='old')
write(6,*) '....Loading file: lc.dat'
do i__=1,563
read(1,'(A36)')ar__
read(ar__,1)HJD(i__),Flux(i__),e_Flux(i__)
c ..............Just test output...........
write(6,1)HJD(i__),Flux(i__),e_Flux(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'rv.dat' ! Radial velocity data for HD 80606 with SOPHIE
C Format for file interpretation
2 format(F13.5,3X,F6.4,1X,F6.4)
C Effective file loading
open(unit=1,file='rv.dat', status='old')
write(6,*) '....Loading file: rv.dat'
do i__=1,48
read(1,'(A29)')ar__1
read(ar__1,2)HJD_1(i__),RV(i__),e_RV(i__)
c ..............Just test output...........
write(6,2)HJD_1(i__),RV(i__),e_RV(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end