Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/405/2570 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/MNRAS/405/2570 Linear polarization of nearby bright stars (Bailey+, 2010) *================================================================================ *The linear polarization of nearby bright stars measured at the parts per *million level. * Bailey J., Lucas P.W., Hough J.H. * <Mon. Not. R. Astron. Soc., 405, 2570-2578 (2010)> * =2010MNRAS.405.2570B C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Properties of sample stars integer*4 nr__ parameter (nr__=49) ! Number of records character*73 ar__ ! Full-size record C J2000 position composed of: RAh RAm DE- DEd DEm real*8 RAdeg ! (deg) Right Ascension J2000 real*8 DEdeg ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 HR ! HR number character*20 OName ! Other names real*4 Vmag ! (mag) V magnitude character*8 SpType ! MK spectral type real*4 Dist ! (pc) Heliocentric distance integer*4 RAh ! (h) Right ascension (J2000) integer*4 RAm ! (min) Right ascension (J2000) character*1 DE_ ! Declination sign (J2000) integer*4 DEd ! (deg) Declination (J2000) integer*4 DEm ! (arcmin) Declination (J2000) real*4 GLON ! (deg) Galactic longitude real*4 GLAT ! (deg) Galactic latitude character*1 l_vsini ! Limit flag on vsini integer*4 vsini ! (km/s) ? Rotational velocity c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Previous polarization measurements integer*4 nr__1 parameter (nr__1=35) ! Number of records character*42 ar__1 ! Full-size record integer*4 HR_1 ! HR number real*4 Pol ! (%) ? Polarization (Heiles, 2000, Cat. II/226) real*4 e_Pol ! (%) ? rms uncertainty on PH integer*4 Q_I_T ! (10-5) ? Q/I Stokes parameter * (Tinbergen, 1982A&A...105...53T) integer*4 e_Q_I_T ! (10-5) ? rms uncertainty on Q/IT integer*4 U_I_T ! (10-5) ? U/I Stokes parameter * (Tinbergen, 1982A&A...105...53T) integer*4 e_U_I_T ! (10-5) ? rms uncertainty on U/IT integer*4 Q_I_P ! (10-5) ? Q/I Stokes parameter * (Piirola, 1977A&AS...30..213P) integer*4 e_Q_I_P ! (10-5) ? rms uncertainty on Q/IP integer*4 U_I_P ! (10-5) ? U/I Stokes parameter * (Piirola, 1977A&AS...30..213P) integer*4 e_U_I_P ! (10-5) ? rms uncertainty on U/IP c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Planetpol linear polarization measurements integer*4 nr__2 parameter (nr__2=49) ! Number of records character*66 ar__2 ! Full-size record integer*4 HR_2 ! HR number character*17 Dates ! Date(s) of observation real*4 Q_I ! (10-6) Q/I Stokes parameter real*4 e_Q_I ! (10-6) rms uncertainty on Q/I real*4 U_I ! (10-6) U/I Stokes parameter real*4 e_U_I ! (10-6) rms uncertainty on U/I real*4 Pol_1 ! (10-6) Polarisation (is calculated from * Pol=sqrt[(Q/I)^2^+(U/I)^2^]) real*4 e_Pol_1 ! (10-6) rms uncertainty on P real*4 theta ! (deg) Polarisation position angle real*4 e_theta ! (deg) rms uncertainty on theta C============================================================================= C Loading file 'table1.dat' ! Properties of sample stars C Format for file interpretation 1 format( + I4,1X,A20,F5.2,1X,A8,1X,F4.1,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X, + F5.1,1X,F4.1,1X,A1,I3) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,49 read(1,'(A73)')ar__ read(ar__,1) + HR,OName,Vmag,SpType,Dist,RAh,RAm,DE_,DEd,DEm,GLON,GLAT, + l_vsini,vsini if(ar__(71:73) .EQ. '') vsini = iNULL__ RAdeg = rNULL__ DEdeg = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh .GT. -180) RAdeg=RAh*15. if(RAm .GT. -180) RAdeg=RAdeg+RAm/4. if(DEd .GE. 0) DEdeg=DEd if(DEm .GE. 0) DEdeg=DEdeg+DEm/60. if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg c ..............Just test output........... write(6,1) + HR,OName,Vmag,SpType,Dist,RAh,RAm,DE_,DEd,DEm,GLON,GLAT, + l_vsini,vsini write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Previous polarization measurements C Format for file interpretation 2 format( + I4,1X,F5.3,1X,F5.3,1X,I3,1X,I1,1X,I3,1X,I1,1X,I3,1X,I2,1X,I3, + 1X,I2) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,35 read(1,'(A42)')ar__1 read(ar__1,2) + HR_1,Pol,e_Pol,Q_I_T,e_Q_I_T,U_I_T,e_U_I_T,Q_I_P,e_Q_I_P, + U_I_P,e_U_I_P if(ar__1(6:10) .EQ. '') Pol = rNULL__ if(ar__1(12:16) .EQ. '') e_Pol = rNULL__ if(ar__1(18:20) .EQ. '') Q_I_T = iNULL__ if(ar__1(22:22) .EQ. '') e_Q_I_T = iNULL__ if(ar__1(24:26) .EQ. '') U_I_T = iNULL__ if(ar__1(28:28) .EQ. '') e_U_I_T = iNULL__ if(ar__1(30:32) .EQ. '') Q_I_P = iNULL__ if(ar__1(34:35) .EQ. '') e_Q_I_P = iNULL__ if(ar__1(37:39) .EQ. '') U_I_P = iNULL__ if(ar__1(41:42) .EQ. '') e_U_I_P = iNULL__ c ..............Just test output........... write(6,2) + HR_1,Pol,e_Pol,Q_I_T,e_Q_I_T,U_I_T,e_U_I_T,Q_I_P,e_Q_I_P, + U_I_P,e_U_I_P c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Planetpol linear polarization measurements C Format for file interpretation 3 format( + I4,1X,A17,F6.1,1X,F3.1,1X,F6.1,1X,F3.1,1X,F6.1,1X,F3.1,1X, + F5.1,1X,F5.1) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,49 read(1,'(A66)')ar__2 read(ar__2,3) + HR_2,Dates,Q_I,e_Q_I,U_I,e_U_I,Pol_1,e_Pol_1,theta,e_theta c ..............Just test output........... write(6,3) + HR_2,Dates,Q_I,e_Q_I,U_I,e_U_I,Pol_1,e_Pol_1,theta,e_theta c .......End.of.Just test output........... end do close(1) C============================================================================= stop end