FORTRAN Generation
(/./ftp/cats/J/PASP/106/949)

Conversion of standardized ReadMe file for file /./ftp/cats/J/PASP/106/949 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-Mar-29
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/PASP/106/949      Multicolor polarimetry of Be stars      (McDavid, 1994)
*================================================================================
*Multicolor polarimetry of selected Be stars: 1990-93
*    McDavid D.
*   <Publ. Astron. Soc. Pac. 106, 949 (1994)>
*   =1994PASP..106..949M      (SIMBAD/NED Reference)
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table1'	! Program Be stars

      integer*4 nr__
      parameter (nr__=8)	! Number of records
      character*44 ar__   	! Full-size record

      character*7   Name       (nr__) ! Star name
      integer*4     HD         (nr__) ! HD (Cat. <III/135>) catalog number
      integer*4     HR         (nr__) ! BS (Cat. <V/50>) catalog number
      real*4        Vmag       (nr__) ! (mag) Visual magnitude (BSC4, See Cat. <V/50>)
      character*15  Sp         (nr__) ! Spectral type (1)
      integer*4     vsini      (nr__) ! (km/s) vsini (1)
*Note (1): According to Slettebak, A. 1982ApJS...50...55S

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2'	! Filter system parameters

      integer*4 nr__1
      parameter (nr__1=5)	! Number of records
      character*11 ar__1  	! Full-size record

      character*1   Passband   (nr__1) ! Filter system passband
      integer*4     Wavelength (nr__1) ! (0.1nm) Effective wavelength
      integer*4     FWHM       (nr__1) ! (0.1nm) Bandpass

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table3'	! Polarization of 2H Cam (standard star)

      integer*4 nr__2
      parameter (nr__2=40)	! Number of records
      character*62 ar__2  	! Full-size record

      integer*4     Month      (nr__2) ! Month of observations
      integer*4     Year_1900  (nr__2) ! (yr) Year (offset 1900) of observation
      character*1   Passband_1 (nr__2) ! [UBVRI] Filter system passband
      integer*4     nobs       (nr__2) ! []? Number of observations
      real*4        q          (nr__2) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q        (nr__2) ! (%) []? rms uncertainty on q
      real*4        u          (nr__2) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u        (nr__2) ! (%) []? rms uncertainty on u
      real*4        p          (nr__2) ! (%) Mean degree of polarization
      real*4        e_p        (nr__2) ! (%) []? rms uncertainty on p
      real*4        theta      (nr__2) ! (deg) Mean equatorial position angle
      real*4        e_theta    (nr__2) ! (deg) []? rms uncertainty on theta
      real*4        dpi        (nr__2) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai    (nr__2) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table4'	! Polarization of {omicron} Sco (standard star)

      integer*4 nr__3
      parameter (nr__3=40)	! Number of records
      character*62 ar__3  	! Full-size record

      integer*4     Month_1    (nr__3) ! Month of observations
      integer*4     Year_1900_1(nr__3) ! (yr) Year (offset 1900) of observation
      character*1   Passband_2 (nr__3) ! [UBVRI] Filter system passband
      integer*4     nobs_1     (nr__3) ! []? Number of observations
      real*4        q_1        (nr__3) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q_1      (nr__3) ! (%) []? rms uncertainty on q
      real*4        u_1        (nr__3) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u_1      (nr__3) ! (%) []? rms uncertainty on u
      real*4        p_1        (nr__3) ! (%) Mean degree of polarization
      real*4        e_p_1      (nr__3) ! (%) []? rms uncertainty on p
      real*4        theta_1    (nr__3) ! (deg) Mean equatorial position angle
      real*4        e_theta_1  (nr__3) ! (deg) []? rms uncertainty on theta
      real*4        dpi_1      (nr__3) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai_1  (nr__3) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table5'	! Polarization of {gamma} Cas

      integer*4 nr__4
      parameter (nr__4=40)	! Number of records
      character*62 ar__4  	! Full-size record

      integer*4     Month_2    (nr__4) ! Month of observations
      integer*4     Year_1900_2(nr__4) ! (yr) Year (offset 1900) of observation
      character*1   Passband_3 (nr__4) ! [UBVRI] Filter system passband
      integer*4     nobs_2     (nr__4) ! []? Number of observations
      real*4        q_2        (nr__4) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q_2      (nr__4) ! (%) []? rms uncertainty on q
      real*4        u_2        (nr__4) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u_2      (nr__4) ! (%) []? rms uncertainty on u
      real*4        p_2        (nr__4) ! (%) Mean degree of polarization
      real*4        e_p_2      (nr__4) ! (%) []? rms uncertainty on p
      real*4        theta_2    (nr__4) ! (deg) Mean equatorial position angle
      real*4        e_theta_2  (nr__4) ! (deg) []? rms uncertainty on theta
      real*4        dpi_2      (nr__4) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai_2  (nr__4) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table6'	! Polarization of {phi} Per

      integer*4 nr__5
      parameter (nr__5=45)	! Number of records
      character*62 ar__5  	! Full-size record

      integer*4     Month_3    (nr__5) ! Month of observations
      integer*4     Year_1900_3(nr__5) ! (yr) Year (offset 1900) of observation
      character*1   Passband_4 (nr__5) ! [UBVRI] Filter system passband
      integer*4     nobs_3     (nr__5) ! []? Number of observations
      real*4        q_3        (nr__5) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q_3      (nr__5) ! (%) []? rms uncertainty on q
      real*4        u_3        (nr__5) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u_3      (nr__5) ! (%) []? rms uncertainty on u
      real*4        p_3        (nr__5) ! (%) Mean degree of polarization
      real*4        e_p_3      (nr__5) ! (%) []? rms uncertainty on p
      real*4        theta_3    (nr__5) ! (deg) Mean equatorial position angle
      real*4        e_theta_3  (nr__5) ! (deg) []? rms uncertainty on theta
      real*4        dpi_3      (nr__5) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai_3  (nr__5) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table7'	! Polarization of 48 Per

      integer*4 nr__6
      parameter (nr__6=45)	! Number of records
      character*62 ar__6  	! Full-size record

      integer*4     Month_4    (nr__6) ! Month of observations
      integer*4     Year_1900_4(nr__6) ! (yr) Year (offset 1900) of observation
      character*1   Passband_5 (nr__6) ! [UBVRI] Filter system passband
      integer*4     nobs_4     (nr__6) ! []? Number of observations
      real*4        q_4        (nr__6) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q_4      (nr__6) ! (%) []? rms uncertainty on q
      real*4        u_4        (nr__6) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u_4      (nr__6) ! (%) []? rms uncertainty on u
      real*4        p_4        (nr__6) ! (%) Mean degree of polarization
      real*4        e_p_4      (nr__6) ! (%) []? rms uncertainty on p
      real*4        theta_4    (nr__6) ! (deg) Mean equatorial position angle
      real*4        e_theta_4  (nr__6) ! (deg) []? rms uncertainty on theta
      real*4        dpi_4      (nr__6) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai_4  (nr__6) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table8'	! Polarization of {zeta} Tau

      integer*4 nr__7
      parameter (nr__7=45)	! Number of records
      character*62 ar__7  	! Full-size record

      integer*4     Month_5    (nr__7) ! Month of observations
      integer*4     Year_1900_5(nr__7) ! (yr) Year (offset 1900) of observation
      character*1   Passband_6 (nr__7) ! [UBVRI] Filter system passband
      integer*4     nobs_5     (nr__7) ! []? Number of observations
      real*4        q_5        (nr__7) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q_5      (nr__7) ! (%) []? rms uncertainty on q
      real*4        u_5        (nr__7) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u_5      (nr__7) ! (%) []? rms uncertainty on u
      real*4        p_5        (nr__7) ! (%) Mean degree of polarization
      real*4        e_p_5      (nr__7) ! (%) []? rms uncertainty on p
      real*4        theta_5    (nr__7) ! (deg) Mean equatorial position angle
      real*4        e_theta_5  (nr__7) ! (deg) []? rms uncertainty on theta
      real*4        dpi_5      (nr__7) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai_5  (nr__7) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table9'	! Polarization of 48 Lib

      integer*4 nr__8
      parameter (nr__8=45)	! Number of records
      character*62 ar__8  	! Full-size record

      integer*4     Month_6    (nr__8) ! Month of observations
      integer*4     Year_1900_6(nr__8) ! (yr) Year (offset 1900) of observation
      character*1   Passband_7 (nr__8) ! [UBVRI] Filter system passband
      integer*4     nobs_6     (nr__8) ! []? Number of observations
      real*4        q_6        (nr__8) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q_6      (nr__8) ! (%) []? rms uncertainty on q
      real*4        u_6        (nr__8) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u_6      (nr__8) ! (%) []? rms uncertainty on u
      real*4        p_6        (nr__8) ! (%) Mean degree of polarization
      real*4        e_p_6      (nr__8) ! (%) []? rms uncertainty on p
      real*4        theta_6    (nr__8) ! (deg) Mean equatorial position angle
      real*4        e_theta_6  (nr__8) ! (deg) []? rms uncertainty on theta
      real*4        dpi_6      (nr__8) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai_6  (nr__8) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table10'	! Polarization of {chi} Oph

      integer*4 nr__9
      parameter (nr__9=45)	! Number of records
      character*62 ar__9  	! Full-size record

      integer*4     Month_7    (nr__9) ! Month of observations
      integer*4     Year_1900_7(nr__9) ! (yr) Year (offset 1900) of observation
      character*1   Passband_8 (nr__9) ! [UBVRI] Filter system passband
      integer*4     nobs_7     (nr__9) ! []? Number of observations
      real*4        q_7        (nr__9) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q_7      (nr__9) ! (%) []? rms uncertainty on q
      real*4        u_7        (nr__9) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u_7      (nr__9) ! (%) []? rms uncertainty on u
      real*4        p_7        (nr__9) ! (%) Mean degree of polarization
      real*4        e_p_7      (nr__9) ! (%) []? rms uncertainty on p
      real*4        theta_7    (nr__9) ! (deg) Mean equatorial position angle
      real*4        e_theta_7  (nr__9) ! (deg) []? rms uncertainty on theta
      real*4        dpi_7      (nr__9) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai_7  (nr__9) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table11'	! Polarization of {pi} Aqr

      integer*4 nr__10
      parameter (nr__10=45)	! Number of records
      character*62 ar__10 	! Full-size record

      integer*4     Month_8    (nr__10) ! Month of observations
      integer*4     Year_1900_8(nr__10) ! (yr) Year (offset 1900) of observation
      character*1   Passband_9 (nr__10) ! [UBVRI] Filter system passband
      integer*4     nobs_8     (nr__10) ! []? Number of observations
      real*4        q_8        (nr__10) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q_8      (nr__10) ! (%) []? rms uncertainty on q
      real*4        u_8        (nr__10) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u_8      (nr__10) ! (%) []? rms uncertainty on u
      real*4        p_8        (nr__10) ! (%) Mean degree of polarization
      real*4        e_p_8      (nr__10) ! (%) []? rms uncertainty on p
      real*4        theta_8    (nr__10) ! (deg) Mean equatorial position angle
      real*4        e_theta_8  (nr__10) ! (deg) []? rms uncertainty on theta
      real*4        dpi_8      (nr__10) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai_8  (nr__10) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table12'	! Polarization of {omicron} And

      integer*4 nr__11
      parameter (nr__11=40)	! Number of records
      character*62 ar__11 	! Full-size record

      integer*4     Month_9    (nr__11) ! Month of observations
      integer*4     Year_1900_9(nr__11) ! (yr) Year (offset 1900) of observation
      character*1   Passband_10(nr__11) ! [UBVRI] Filter system passband
      integer*4     nobs_9     (nr__11) ! []? Number of observations
      real*4        q_9        (nr__11) ! (%) []? Mean normalized Stokes parameter q
      real*4        e_q_9      (nr__11) ! (%) []? rms uncertainty on q
      real*4        u_9        (nr__11) ! (%) []? Mean normalized Stokes parameter u
      real*4        e_u_9      (nr__11) ! (%) []? rms uncertainty on u
      real*4        p_9        (nr__11) ! (%) Mean degree of polarization
      real*4        e_p_9      (nr__11) ! (%) []? rms uncertainty on p
      real*4        theta_9    (nr__11) ! (deg) Mean equatorial position angle
      real*4        e_theta_9  (nr__11) ! (deg) []? rms uncertainty on theta
      real*4        dpi_9      (nr__11) ! (%) []? Instrumental precision of q, u, or p (1)
      real*4        dthetai_9  (nr__11) ! (deg) []? Instrumental precision of theta (2)
*Note (1): dpi = mean instrumental error in the degree of polarization for each
*           set of measurements according to photon statistics, unless photon
*           statistics gives a value less than 0.02%, in which case the error is
*           approximated as 0.02% due to atmospheric scintillation for the
*           200-second integration time
*Note (2): dthetai = mean instrumental error in equatorial position angle of the
*           polarization, based on the relation dthetai=28.65(dpi/p)

C=============================================================================

C  Loading file 'table1'	! Program Be stars

C  Format for file interpretation

    1 format(A7,1X,I6,1X,I4,1X,F4.2,1X,A15,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1')
      write(6,*) '....Loading file: table1'
      do i__=1,8
        read(1,'(A44)')ar__
        read(ar__,1)
     +  Name(i__),HD(i__),HR(i__),Vmag(i__),Sp(i__),vsini(i__)
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),HD(i__),HR(i__),Vmag(i__),Sp(i__),vsini(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table2'	! Filter system parameters

C  Format for file interpretation

    2 format(A1,1X,I4,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2')
      write(6,*) '....Loading file: table2'
      do i__=1,5
        read(1,'(A11)')ar__1
        read(ar__1,2)Passband(i__),Wavelength(i__),FWHM(i__)
c    ..............Just test output...........
        write(6,2)Passband(i__),Wavelength(i__),FWHM(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table3'	! Polarization of 2H Cam (standard star)

C  Format for file interpretation

    3 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3')
      write(6,*) '....Loading file: table3'
      do i__=1,40
        read(1,'(A62)')ar__2
        read(ar__2,3)
     +  Month(i__),Year_1900(i__),Passband_1(i__),nobs(i__),q(i__),
     +  e_q(i__),u(i__),e_u(i__),p(i__),e_p(i__),theta(i__),
     +  e_theta(i__),dpi(i__),dthetai(i__)
        if(ar__2(9:10) .EQ. '') nobs(i__) = iNULL__
        if(ar__2(12:16) .EQ. '') q(i__) = rNULL__
        if(ar__2(18:21) .EQ. '') e_q(i__) = rNULL__
        if(ar__2(23:27) .EQ. '') u(i__) = rNULL__
        if(ar__2(29:32) .EQ. '') e_u(i__) = rNULL__
        if(ar__2(39:42) .EQ. '') e_p(i__) = rNULL__
        if(ar__2(50:53) .EQ. '') e_theta(i__) = rNULL__
        if(ar__2(55:58) .EQ. '') dpi(i__) = rNULL__
        if(ar__2(60:62) .EQ. '') dthetai(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Month(i__),Year_1900(i__),Passband_1(i__),nobs(i__),q(i__),
     +  e_q(i__),u(i__),e_u(i__),p(i__),e_p(i__),theta(i__),
     +  e_theta(i__),dpi(i__),dthetai(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table4'	! Polarization of {omicron} Sco (standard star)

C  Format for file interpretation

    4 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4')
      write(6,*) '....Loading file: table4'
      do i__=1,40
        read(1,'(A62)')ar__3
        read(ar__3,4)
     +  Month_1(i__),Year_1900_1(i__),Passband_2(i__),nobs_1(i__),
     +  q_1(i__),e_q_1(i__),u_1(i__),e_u_1(i__),p_1(i__),e_p_1(i__),
     +  theta_1(i__),e_theta_1(i__),dpi_1(i__),dthetai_1(i__)
        if(ar__3(9:10) .EQ. '') nobs_1(i__) = iNULL__
        if(ar__3(12:16) .EQ. '') q_1(i__) = rNULL__
        if(ar__3(18:21) .EQ. '') e_q_1(i__) = rNULL__
        if(ar__3(23:27) .EQ. '') u_1(i__) = rNULL__
        if(ar__3(29:32) .EQ. '') e_u_1(i__) = rNULL__
        if(ar__3(39:42) .EQ. '') e_p_1(i__) = rNULL__
        if(ar__3(50:53) .EQ. '') e_theta_1(i__) = rNULL__
        if(ar__3(55:58) .EQ. '') dpi_1(i__) = rNULL__
        if(ar__3(60:62) .EQ. '') dthetai_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Month_1(i__),Year_1900_1(i__),Passband_2(i__),nobs_1(i__),
     +  q_1(i__),e_q_1(i__),u_1(i__),e_u_1(i__),p_1(i__),e_p_1(i__),
     +  theta_1(i__),e_theta_1(i__),dpi_1(i__),dthetai_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table5'	! Polarization of {gamma} Cas

C  Format for file interpretation

    5 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5')
      write(6,*) '....Loading file: table5'
      do i__=1,40
        read(1,'(A62)')ar__4
        read(ar__4,5)
     +  Month_2(i__),Year_1900_2(i__),Passband_3(i__),nobs_2(i__),
     +  q_2(i__),e_q_2(i__),u_2(i__),e_u_2(i__),p_2(i__),e_p_2(i__),
     +  theta_2(i__),e_theta_2(i__),dpi_2(i__),dthetai_2(i__)
        if(ar__4(9:10) .EQ. '') nobs_2(i__) = iNULL__
        if(ar__4(12:16) .EQ. '') q_2(i__) = rNULL__
        if(ar__4(18:21) .EQ. '') e_q_2(i__) = rNULL__
        if(ar__4(23:27) .EQ. '') u_2(i__) = rNULL__
        if(ar__4(29:32) .EQ. '') e_u_2(i__) = rNULL__
        if(ar__4(39:42) .EQ. '') e_p_2(i__) = rNULL__
        if(ar__4(50:53) .EQ. '') e_theta_2(i__) = rNULL__
        if(ar__4(55:58) .EQ. '') dpi_2(i__) = rNULL__
        if(ar__4(60:62) .EQ. '') dthetai_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Month_2(i__),Year_1900_2(i__),Passband_3(i__),nobs_2(i__),
     +  q_2(i__),e_q_2(i__),u_2(i__),e_u_2(i__),p_2(i__),e_p_2(i__),
     +  theta_2(i__),e_theta_2(i__),dpi_2(i__),dthetai_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table6'	! Polarization of {phi} Per

C  Format for file interpretation

    6 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6')
      write(6,*) '....Loading file: table6'
      do i__=1,45
        read(1,'(A62)')ar__5
        read(ar__5,6)
     +  Month_3(i__),Year_1900_3(i__),Passband_4(i__),nobs_3(i__),
     +  q_3(i__),e_q_3(i__),u_3(i__),e_u_3(i__),p_3(i__),e_p_3(i__),
     +  theta_3(i__),e_theta_3(i__),dpi_3(i__),dthetai_3(i__)
        if(ar__5(9:10) .EQ. '') nobs_3(i__) = iNULL__
        if(ar__5(12:16) .EQ. '') q_3(i__) = rNULL__
        if(ar__5(18:21) .EQ. '') e_q_3(i__) = rNULL__
        if(ar__5(23:27) .EQ. '') u_3(i__) = rNULL__
        if(ar__5(29:32) .EQ. '') e_u_3(i__) = rNULL__
        if(ar__5(39:42) .EQ. '') e_p_3(i__) = rNULL__
        if(ar__5(50:53) .EQ. '') e_theta_3(i__) = rNULL__
        if(ar__5(55:58) .EQ. '') dpi_3(i__) = rNULL__
        if(ar__5(60:62) .EQ. '') dthetai_3(i__) = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  Month_3(i__),Year_1900_3(i__),Passband_4(i__),nobs_3(i__),
     +  q_3(i__),e_q_3(i__),u_3(i__),e_u_3(i__),p_3(i__),e_p_3(i__),
     +  theta_3(i__),e_theta_3(i__),dpi_3(i__),dthetai_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table7'	! Polarization of 48 Per

C  Format for file interpretation

    7 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7')
      write(6,*) '....Loading file: table7'
      do i__=1,45
        read(1,'(A62)')ar__6
        read(ar__6,7)
     +  Month_4(i__),Year_1900_4(i__),Passband_5(i__),nobs_4(i__),
     +  q_4(i__),e_q_4(i__),u_4(i__),e_u_4(i__),p_4(i__),e_p_4(i__),
     +  theta_4(i__),e_theta_4(i__),dpi_4(i__),dthetai_4(i__)
        if(ar__6(9:10) .EQ. '') nobs_4(i__) = iNULL__
        if(ar__6(12:16) .EQ. '') q_4(i__) = rNULL__
        if(ar__6(18:21) .EQ. '') e_q_4(i__) = rNULL__
        if(ar__6(23:27) .EQ. '') u_4(i__) = rNULL__
        if(ar__6(29:32) .EQ. '') e_u_4(i__) = rNULL__
        if(ar__6(39:42) .EQ. '') e_p_4(i__) = rNULL__
        if(ar__6(50:53) .EQ. '') e_theta_4(i__) = rNULL__
        if(ar__6(55:58) .EQ. '') dpi_4(i__) = rNULL__
        if(ar__6(60:62) .EQ. '') dthetai_4(i__) = rNULL__
c    ..............Just test output...........
        write(6,7)
     +  Month_4(i__),Year_1900_4(i__),Passband_5(i__),nobs_4(i__),
     +  q_4(i__),e_q_4(i__),u_4(i__),e_u_4(i__),p_4(i__),e_p_4(i__),
     +  theta_4(i__),e_theta_4(i__),dpi_4(i__),dthetai_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table8'	! Polarization of {zeta} Tau

C  Format for file interpretation

    8 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8')
      write(6,*) '....Loading file: table8'
      do i__=1,45
        read(1,'(A62)')ar__7
        read(ar__7,8)
     +  Month_5(i__),Year_1900_5(i__),Passband_6(i__),nobs_5(i__),
     +  q_5(i__),e_q_5(i__),u_5(i__),e_u_5(i__),p_5(i__),e_p_5(i__),
     +  theta_5(i__),e_theta_5(i__),dpi_5(i__),dthetai_5(i__)
        if(ar__7(9:10) .EQ. '') nobs_5(i__) = iNULL__
        if(ar__7(12:16) .EQ. '') q_5(i__) = rNULL__
        if(ar__7(18:21) .EQ. '') e_q_5(i__) = rNULL__
        if(ar__7(23:27) .EQ. '') u_5(i__) = rNULL__
        if(ar__7(29:32) .EQ. '') e_u_5(i__) = rNULL__
        if(ar__7(39:42) .EQ. '') e_p_5(i__) = rNULL__
        if(ar__7(50:53) .EQ. '') e_theta_5(i__) = rNULL__
        if(ar__7(55:58) .EQ. '') dpi_5(i__) = rNULL__
        if(ar__7(60:62) .EQ. '') dthetai_5(i__) = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  Month_5(i__),Year_1900_5(i__),Passband_6(i__),nobs_5(i__),
     +  q_5(i__),e_q_5(i__),u_5(i__),e_u_5(i__),p_5(i__),e_p_5(i__),
     +  theta_5(i__),e_theta_5(i__),dpi_5(i__),dthetai_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table9'	! Polarization of 48 Lib

C  Format for file interpretation

    9 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9')
      write(6,*) '....Loading file: table9'
      do i__=1,45
        read(1,'(A62)')ar__8
        read(ar__8,9)
     +  Month_6(i__),Year_1900_6(i__),Passband_7(i__),nobs_6(i__),
     +  q_6(i__),e_q_6(i__),u_6(i__),e_u_6(i__),p_6(i__),e_p_6(i__),
     +  theta_6(i__),e_theta_6(i__),dpi_6(i__),dthetai_6(i__)
        if(ar__8(9:10) .EQ. '') nobs_6(i__) = iNULL__
        if(ar__8(12:16) .EQ. '') q_6(i__) = rNULL__
        if(ar__8(18:21) .EQ. '') e_q_6(i__) = rNULL__
        if(ar__8(23:27) .EQ. '') u_6(i__) = rNULL__
        if(ar__8(29:32) .EQ. '') e_u_6(i__) = rNULL__
        if(ar__8(39:42) .EQ. '') e_p_6(i__) = rNULL__
        if(ar__8(50:53) .EQ. '') e_theta_6(i__) = rNULL__
        if(ar__8(55:58) .EQ. '') dpi_6(i__) = rNULL__
        if(ar__8(60:62) .EQ. '') dthetai_6(i__) = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  Month_6(i__),Year_1900_6(i__),Passband_7(i__),nobs_6(i__),
     +  q_6(i__),e_q_6(i__),u_6(i__),e_u_6(i__),p_6(i__),e_p_6(i__),
     +  theta_6(i__),e_theta_6(i__),dpi_6(i__),dthetai_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table10'	! Polarization of {chi} Oph

C  Format for file interpretation

   10 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10')
      write(6,*) '....Loading file: table10'
      do i__=1,45
        read(1,'(A62)')ar__9
        read(ar__9,10)
     +  Month_7(i__),Year_1900_7(i__),Passband_8(i__),nobs_7(i__),
     +  q_7(i__),e_q_7(i__),u_7(i__),e_u_7(i__),p_7(i__),e_p_7(i__),
     +  theta_7(i__),e_theta_7(i__),dpi_7(i__),dthetai_7(i__)
        if(ar__9(9:10) .EQ. '') nobs_7(i__) = iNULL__
        if(ar__9(12:16) .EQ. '') q_7(i__) = rNULL__
        if(ar__9(18:21) .EQ. '') e_q_7(i__) = rNULL__
        if(ar__9(23:27) .EQ. '') u_7(i__) = rNULL__
        if(ar__9(29:32) .EQ. '') e_u_7(i__) = rNULL__
        if(ar__9(39:42) .EQ. '') e_p_7(i__) = rNULL__
        if(ar__9(50:53) .EQ. '') e_theta_7(i__) = rNULL__
        if(ar__9(55:58) .EQ. '') dpi_7(i__) = rNULL__
        if(ar__9(60:62) .EQ. '') dthetai_7(i__) = rNULL__
c    ..............Just test output...........
        write(6,10)
     +  Month_7(i__),Year_1900_7(i__),Passband_8(i__),nobs_7(i__),
     +  q_7(i__),e_q_7(i__),u_7(i__),e_u_7(i__),p_7(i__),e_p_7(i__),
     +  theta_7(i__),e_theta_7(i__),dpi_7(i__),dthetai_7(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table11'	! Polarization of {pi} Aqr

C  Format for file interpretation

   11 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table11')
      write(6,*) '....Loading file: table11'
      do i__=1,45
        read(1,'(A62)')ar__10
        read(ar__10,11)
     +  Month_8(i__),Year_1900_8(i__),Passband_9(i__),nobs_8(i__),
     +  q_8(i__),e_q_8(i__),u_8(i__),e_u_8(i__),p_8(i__),e_p_8(i__),
     +  theta_8(i__),e_theta_8(i__),dpi_8(i__),dthetai_8(i__)
        if(ar__10(9:10) .EQ. '') nobs_8(i__) = iNULL__
        if(ar__10(12:16) .EQ. '') q_8(i__) = rNULL__
        if(ar__10(18:21) .EQ. '') e_q_8(i__) = rNULL__
        if(ar__10(23:27) .EQ. '') u_8(i__) = rNULL__
        if(ar__10(29:32) .EQ. '') e_u_8(i__) = rNULL__
        if(ar__10(39:42) .EQ. '') e_p_8(i__) = rNULL__
        if(ar__10(50:53) .EQ. '') e_theta_8(i__) = rNULL__
        if(ar__10(55:58) .EQ. '') dpi_8(i__) = rNULL__
        if(ar__10(60:62) .EQ. '') dthetai_8(i__) = rNULL__
c    ..............Just test output...........
        write(6,11)
     +  Month_8(i__),Year_1900_8(i__),Passband_9(i__),nobs_8(i__),
     +  q_8(i__),e_q_8(i__),u_8(i__),e_u_8(i__),p_8(i__),e_p_8(i__),
     +  theta_8(i__),e_theta_8(i__),dpi_8(i__),dthetai_8(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table12'	! Polarization of {omicron} And

C  Format for file interpretation

   12 format(
     +  I2,1X,I2,1X,A1,1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F4.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table12')
      write(6,*) '....Loading file: table12'
      do i__=1,40
        read(1,'(A62)')ar__11
        read(ar__11,12)
     +  Month_9(i__),Year_1900_9(i__),Passband_10(i__),nobs_9(i__),
     +  q_9(i__),e_q_9(i__),u_9(i__),e_u_9(i__),p_9(i__),e_p_9(i__),
     +  theta_9(i__),e_theta_9(i__),dpi_9(i__),dthetai_9(i__)
        if(ar__11(9:10) .EQ. '') nobs_9(i__) = iNULL__
        if(ar__11(12:16) .EQ. '') q_9(i__) = rNULL__
        if(ar__11(18:21) .EQ. '') e_q_9(i__) = rNULL__
        if(ar__11(23:27) .EQ. '') u_9(i__) = rNULL__
        if(ar__11(29:32) .EQ. '') e_u_9(i__) = rNULL__
        if(ar__11(39:42) .EQ. '') e_p_9(i__) = rNULL__
        if(ar__11(50:53) .EQ. '') e_theta_9(i__) = rNULL__
        if(ar__11(55:58) .EQ. '') dpi_9(i__) = rNULL__
        if(ar__11(60:62) .EQ. '') dthetai_9(i__) = rNULL__
c    ..............Just test output...........
        write(6,12)
     +  Month_9(i__),Year_1900_9(i__),Passband_10(i__),nobs_9(i__),
     +  q_9(i__),e_q_9(i__),u_9(i__),e_u_9(i__),p_9(i__),e_p_9(i__),
     +  theta_9(i__),e_theta_9(i__),dpi_9(i__),dthetai_9(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end