Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/384/879 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-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/384/879 Galactic orbits of stars with planets (Barbieri+, 2002) *================================================================================ *Galactic orbits of stars with planets * Barbieri M., Gratton R.G. * <Astron. Astrophys. 384, 879 (2002)> * =2002A&A...384..879B C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Spatial properties of stars with planets (SWP) integer*4 nr__ parameter (nr__=58) ! Number of records character*69 ar__ ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*10 Name (nr__) ! Name of the star real*4 plx (nr__) ! (arcsec) ? Parallax real*8 HRV (nr__) ! (km/s) ? Heliocentric radial velocity real*8 pmRA (nr__) ! (mas/yr) Proper motion in right ascension real*8 pmDE (nr__) ! (mas/yr) Proper motion declination integer*4 RAh (nr__) ! (h) Right ascension (J2000) integer*4 RAm (nr__) ! (min) Right ascension (J2000) real*4 RAs (nr__) ! (s) Right ascension (J2000) character*1 DE_ (nr__) ! Declination sign (J2000) integer*4 DEd (nr__) ! (deg) Declination (J2000) integer*4 DEm (nr__) ! (arcmin) Declination (J2000) real*4 DEs (nr__) ! (arcsec) Declination (J2000) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Physical properties of stars with planets (SWP) integer*4 nr__1 parameter (nr__1=58) ! Number of records character*61 ar__1 ! Full-size record character*10 Name_1 (nr__1) ! Name of the star character*5 SpType (nr__1) ! Spectral type character*5 LClass (nr__1) ! Luminosity class real*4 Vmag (nr__1) ! (mag) ? Visual magnitude real*4 v_Fe_H_ (nr__1) ! ([Sun]) ? Metallicity real*4 e__Fe_H_ (nr__1) ! ([Sun]) ? Metallicity error integer*4 r__Fe_H_ (nr__1) ! [1/11]? Reference for metallicities (1) real*4 Age (nr__1) ! (Gyr) ? Age real*4 e_Age (nr__1) ! (Gyr) ? Age error real*4 Mass (nr__1) ! (Sun) Mass real*4 e_Mass (nr__1) ! (Sun) ? Mass error *Note (1): References for metallicities of SWP * 1 = Santos, Israelian & Mayor, 2001A&A...373.1019S * 2 = CORAVEL web site: http://obsunige.ch/~udry/planet/planet.html * 3 = Gonzalez & Laws, 2000AJ....119..390G * 4 = Gonzales, Laws, Tyagi & Reddy, 2001AJ....121..432G * 5 = Gonzalez & Vanture, 1998A&A...339L..29G * 6 = Naef et al., 2001, Cat. <J/A+A/375/L27> * 7 = Gonzalez, 1998A&A...334..221G * 8 = Gonzalez, 1999MNRAS.308..447G * 9 = Cayrel de Strobel, Soubiran & Ralite, 2001, Cat. <III/221> * 10 = Marsakov & Shevelev, 1988, Cat. <V/58> * 11 = Laws & Gonzalez, 2001ApJ...553..405L c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Orbital data of stars with planets (SWP) integer*4 nr__2 parameter (nr__2=56) ! Number of records character*48 ar__2 ! Full-size record character*10 Name_2 (nr__2) ! Name of the star real*4 Distpg (nr__2) ! (kpc) Perigalactic distance real*4 Distag (nr__2) ! (kpc) Apogalactic distance real*4 Zmax (nr__2) ! (kpc) Maximum height of the orbit on Galactic Plane real*4 e (nr__2) ! Orbital eccentricity real*8 Energy (nr__2) ! (km2/s2) Orbital energy real*4 Jz (nr__2) ! (kpc.km/s) Orbital angular moment C============================================================================= C Loading file 'table1.dat' ! Spatial properties of stars with planets (SWP) C Format for file interpretation 1 format( + A10,2X,F6.2,1X,F7.3,2X,F7.2,2X,F7.2,2X,I2,1X,I2,1X,F5.2,1X,A1, + I2,1X,I2,1X,F4.1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,58 read(1,'(A69)')ar__ read(ar__,1) + Name(i__),plx(i__),HRV(i__),pmRA(i__),pmDE(i__),RAh(i__), + RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__) if(ar__(13:18) .EQ. '') plx(i__) = rNULL__ if(ar__(20:26) .EQ. '') HRV(i__) = rNULL__ RAdeg(i__) = rNULL__ DEdeg(i__) = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15. if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4. if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240. if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__) if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60. if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,1) + Name(i__),plx(i__),HRV(i__),pmRA(i__),pmDE(i__),RAh(i__), + RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Physical properties of stars with planets (SWP) C Format for file interpretation 2 format( + A10,1X,A5,1X,A5,1X,F4.2,1X,F5.2,1X,F4.2,2X,I2,1X,F4.1,1X,F3.1, + 1X,F4.2,1X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,58 read(1,'(A61)')ar__1 read(ar__1,2) + Name_1(i__),SpType(i__),LClass(i__),Vmag(i__),v_Fe_H_(i__), + e__Fe_H_(i__),r__Fe_H_(i__),Age(i__),e_Age(i__),Mass(i__), + e_Mass(i__) if(ar__1(24:27) .EQ. '') Vmag(i__) = rNULL__ if(ar__1(29:33) .EQ. '') v_Fe_H_(i__) = rNULL__ if(ar__1(35:38) .EQ. '') e__Fe_H_(i__) = rNULL__ if(ar__1(41:42) .EQ. '') r__Fe_H_(i__) = iNULL__ if(ar__1(44:47) .EQ. '') Age(i__) = rNULL__ if(ar__1(49:51) .EQ. '') e_Age(i__) = rNULL__ if(ar__1(58:61) .EQ. '') e_Mass(i__) = rNULL__ c ..............Just test output........... write(6,2) + Name_1(i__),SpType(i__),LClass(i__),Vmag(i__),v_Fe_H_(i__), + e__Fe_H_(i__),r__Fe_H_(i__),Age(i__),e_Age(i__),Mass(i__), + e_Mass(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Orbital data of stars with planets (SWP) C Format for file interpretation 3 format(A10,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F9.1,1X,F6.1) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,56 read(1,'(A48)')ar__2 read(ar__2,3) + Name_2(i__),Distpg(i__),Distag(i__),Zmax(i__),e(i__), + Energy(i__),Jz(i__) c ..............Just test output........... write(6,3) + Name_2(i__),Distpg(i__),Distag(i__),Zmax(i__),e(i__), + Energy(i__),Jz(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end