Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/471/165 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-Mar-28 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/471/165 UBVRI polarisation in NGC 5606 (Orsatti+, 2007) *================================================================================ *Revisiting the region of the open cluster NGC 5606. * Orsatti A.M., Feinstein C., Vega E.I., Vergne M.M. * <Astron. Astrophys., 471, 165-172 (2007)> * =2007A&A...471..165O C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Polarimetric observations towards NGC 5606 integer*4 nr__ parameter (nr__=54) ! Number of records character*118 ar__ ! Full-size record integer*4 VF ! Star identification from Vazquez & Feinstein * (1991A&AS...87..383V), * Cl* NGC 5606 VF NN in Simbad character*1 n_VF ! [*] * NGC 5606 members for Vazquez & Feinstein character*4 MV ! Identifications from Moffat & Vogt * (1973A&AS...10..135M), MVNN (1) real*4 PolU ! (%) ?=- Polarisation U filter real*4 e_PolU ! (%) ?=- rms uncertainty on PolU real*4 PAU ! (deg) ?=- Polarisation position angle in U filter real*4 e_PAU ! (deg) ?=- rms uncertainty on PAU integer*4 n ! Number of iterations real*4 PolB ! (%) ?=- Polarisation B filter real*4 e_PolB ! (%) ?=- rms uncertainty on PolB real*4 PAB ! (deg) ?=- Polarisation position angle in B filter real*4 e_PAB ! (deg) ?=- rms uncertainty on PAB real*4 PolV ! (%) Polarisation in V filter real*4 e_PolV ! (%) rms uncertainty on PolV real*4 PAV ! (deg) Polarisation position angle in V filter real*4 e_PAV ! (deg) rms uncertainty on PAV real*4 PolR ! (%) Polarisation in R filter real*4 e_PolR ! (%) rms uncertainty on PolR real*4 PAR ! (deg) Polarisation position angle in Cousins R filter real*4 e_PAR ! (deg) rms uncertainty on PAR real*4 PolI ! (%) Polarisation in Cousins I filter real*4 e_PolI ! (%) rms uncertainty on PolI real*4 PAI ! (deg) Polarisation position angle in Cousins I filter real*4 e_PAI ! (deg) rms uncertainty on PAI *Note (1): MVNN = NGC 5606 NN in Simbad c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Polarization results of stars towards NGC 5606 integer*4 nr__1 parameter (nr__1=54) ! Number of records character*54 ar__1 ! Full-size record integer*4 VF_1 ! Star identifications from Vazquez & Feinstein * (1991A&AS...87..383V), * Cl* NGC 5606 VF NN in Simbad character*1 n_VF_1 ! [*] * NGC 5606 members for Vazquez & Feinstein real*4 Pmax ! (%) Maximum polarisation real*4 e_Pmax ! (%) rms uncertainty on Pmax real*4 sigma1 ! Unit weight error of the fit (1) real*4 lam_max ! (um) Wavelength of maximum polarisation real*4 e_lam_max ! (um) rms uncertainty on lambdamax character*5 SpType ! MK Spectral type from the WEBDA database character*19 Memb ! Membership revision *Note (1): (sigma1)^2^={Sum}[(r_{lambda}_/P_{lambda}_)^2^]/(m-2) * where m is the number of colors and * r_{lambda}_=P_{lambda}_-P_max_*exp(-K*[ln({lambda}_max_/{lambda})]^2^) C============================================================================= C Loading file 'table1.dat' ! Polarimetric observations towards NGC 5606 C Format for file interpretation 1 format( + 3X,I2,A1,1X,A4,1X,F4.2,1X,F4.2,1X,F5.1,1X,F4.1,1X,I1,1X,F4.2, + 1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F4.2,1X,F5.1,1X,F4.1,1X, + F4.2,1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F4.2,1X,F5.1,1X,F4.1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,54 read(1,'(A118)')ar__ read(ar__,1) + VF,n_VF,MV,PolU,e_PolU,PAU,e_PAU,n,PolB,e_PolB,PAB,e_PAB,PolV, + e_PolV,PAV,e_PAV,PolR,e_PolR,PAR,e_PAR,PolI,e_PolI,PAI,e_PAI if (idig(ar__(13:16)).EQ.0) PolU = rNULL__ if (idig(ar__(18:21)).EQ.0) e_PolU = rNULL__ if (idig(ar__(23:27)).EQ.0) PAU = rNULL__ if (idig(ar__(29:32)).EQ.0) e_PAU = rNULL__ if (idig(ar__(36:39)).EQ.0) PolB = rNULL__ if (idig(ar__(41:44)).EQ.0) e_PolB = rNULL__ if (idig(ar__(46:50)).EQ.0) PAB = rNULL__ if (idig(ar__(52:55)).EQ.0) e_PAB = rNULL__ c ..............Just test output........... write(6,1) + VF,n_VF,MV,PolU,e_PolU,PAU,e_PAU,n,PolB,e_PolB,PAB,e_PAB,PolV, + e_PolV,PAV,e_PAV,PolR,e_PolR,PAR,e_PAR,PolI,e_PolI,PAI,e_PAI c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Polarization results of stars towards NGC 5606 C Format for file interpretation 2 format( + I2,A1,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A5,1X,A19) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,54 read(1,'(A54)')ar__1 read(ar__1,2) + VF_1,n_VF_1,Pmax,e_Pmax,sigma1,lam_max,e_lam_max,SpType,Memb c ..............Just test output........... write(6,2) + VF_1,n_VF_1,Pmax,e_Pmax,sigma1,lam_max,e_lam_max,SpType,Memb c .......End.of.Just test output........... end do close(1) C============================================================================= stop end C Locate position of first digit in string; or return 0 integer function idig(c) character*(*) c character*1 c1 integer lc,i lc=len(c) idig=0 do i=1,lc if(c(i:i).ne.' ') go to 1 end do 1 if(i.gt.lc) return c1=c(i:i) if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1 if(i.gt.lc) return c1=c(i:i) if(c1.ge.'0'.and.c1.le.'9') idig=i return end