FORTRAN Generation
(/./ftp/cats/J/A_A/331/633)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/331/633 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/A+A/331/633       IUE spectra of lambda Boo stars (Solano+ 1998)
*================================================================================
*Identification of Lambda Bootis stars using IUE spectra. I. Low Resolution Data
*       Solano E., Paunzen E.
*      <Astron. Astrophys. 331, 633 (1998)>
*      =1998A&A...331..633S      (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Lambda Bootis stars included in Paper I
                              (Pauzen et al., 1997, Cat. <J/A+AS/123/93>)
                              with ultraviolet observations

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension J2000
      real*8        DEdeg       ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*9   Name        ! Star identification
      integer*4     Teff        ! (K) ? Effective temperature (this work)
      integer*4     Teff2       ! (K) ? Effective temperature from other work (1)
      integer*4     r_Teff2     ! ? Reference for Teff2
      integer*4     Teff3       ! (K) ? Effective temperature from other work (1)
      integer*4     r_Teff3     ! ? Reference for Teff3
      integer*4     Teff4       ! (K) ? Effective temperature from other work (1)
      integer*4     r_Teff4     ! ? Reference for Teff4
      real*4        log_g       ! ([cm/s2]) ? Surface gravity
      integer*4     r_log_g     ! ? Reference for log(g) (2)
      character*9   Images      ! IUE Images (SWP/LWP). (3)
      real*4        EW_1657_2800 ! (0.1nm) ? Equivalent width at 1657/2800A
      real*4        EW_1670     ! (0.1nm) ? Equivalent width at 1670A
      real*4        EW_1741     ! (0.1nm) ? Equivalent width at 1741A
      real*4        EW_1931     ! (0.1nm) ? Equivalent width at 1931A
      integer*4     RAh         ! (h) Right Ascension J2000 (hours)
      integer*4     RAm         ! (min) Right Ascension J2000 (minutes)
      real*4        RAs         ! (s) Right Ascension J2000 (seconds)
      character*1   DE_         ! Declination J2000 (sign)
      integer*4     DEd         ! (deg) Declination J2000 (degrees)
      integer*4     DEm         ! (arcmin) Declination J2000 (minutes)
      integer*4     DEs         ! (arcsec) ?Declination J2000 (seconds)
*Note (1): Sources for T_eff_:
*       1: Iliev & Barzova (1995A&A...302..735I)
*       2: Stuerenburg (1993A&A...277..139S)
*       3: Venn & Lambert (1990ApJ...363..234V)
*       4: North et al. (1994A&A...281..775N)
*       5: Cacciari et al. (1987A&A...183..314C)
*       6: Huenemoerder et al. (1984AJ.....89..851H)
*       7: De Boer et al. (1997A&A...317L..23)
*       8: Adelman & Philip (1990MNRAS.247..132A)
*       9: Hayes & Philip (1988PASP..100..801H)
*      11: Holweger & Rentzsch-Holm (1995A&A...303..819H)
*      12: Smalley & Dworetsky (1995A&A...293..446S)
*      13: Blackwell & Lynas-Gray (1994, Cat. <J/A+A/282/899>)
*      14: Alonso et al. (1996A&A...313..873A)
*      15: Malagnini & Morossi (1990A&AS...85.1015M)
*      16: Smalley (1993MNRAS.265.1035S)
*      17: Smalley & Dworetsky (1993A&A...271..515S)
*      18: Napiwotzki et al. (1993A&A...268..653N)
*      19: Adelman et al. (1980ApJS...43..491A)
*      20: Nissen (1981A&A....97..145N)
*      21: Van Veer-Menneret et al. (1989A&A...224..171V)
*      22: Hill & Landstreet (1993, Cat. <J/A+A/276/142>)
*      23: Gerbaldi et al. (1989A&AS...81..127G)
*      24: Adelman (1978ApJ...222..547A)
*Note (2): Sources for log(g):
*       6: Huenemoerder et al. (1984AJ.....89..851H)
*       9: Hayes & Philip (1988PASP..100..801H)
* nothing: this work
*Note (3): Camera and image number. The L and S following the image number
*           indicates the aperture used (Large or Small)

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

C  Declarations for 'table2.dat'	! Same as table1.dat but for the field horizontal
                              branch stars

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_1     ! (deg) Right Ascension J2000
      real*8        DEdeg_1     ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*9   Name_1      ! Star identification
      integer*4     Teff_1      ! (K) ? Effective temperature (this work)
      integer*4     Teff2_1     ! (K) ? Effective temperature from other work (1)
      integer*4     r_Teff2_1   ! ? Reference for Teff2
      integer*4     Teff3_1     ! (K) ? Effective temperature from other work (1)
      integer*4     r_Teff3_1   ! ? Reference for Teff3
      integer*4     Teff4_1     ! (K) ? Effective temperature from other work (1)
      integer*4     r_Teff4_1   ! ? Reference for Teff4
      real*4        log_g_1     ! ([cm/s2]) ? Surface gravity
      integer*4     r_log_g_1   ! ? Reference for log(g) (2)
      character*9   Images_1    ! IUE Images (SWP/LWP). (3)
      real*4        EW_1657_2800_1 ! (0.1nm) ? Equivalent width at 1657/2800A
      real*4        EW_1670_1   ! (0.1nm) ? Equivalent width at 1670A
      real*4        EW_1741_1   ! (0.1nm) ? Equivalent width at 1741A
      real*4        EW_1931_1   ! (0.1nm) ? Equivalent width at 1931A
      integer*4     RAh_1       ! (h) Right Ascension J2000 (hours)
      integer*4     RAm_1       ! (min) Right Ascension J2000 (minutes)
      real*4        RAs_1       ! (s) Right Ascension J2000 (seconds)
      character*1   DE__1       ! Declination J2000 (sign)
      integer*4     DEd_1       ! (deg) Declination J2000 (degrees)
      integer*4     DEm_1       ! (arcmin) Declination J2000 (minutes)
      integer*4     DEs_1       ! (arcsec) ?Declination J2000 (seconds)
*Note (1): Sources for T_eff_:
*       1: Iliev & Barzova (1995A&A...302..735I)
*       2: Stuerenburg (1993A&A...277..139S)
*       3: Venn & Lambert (1990ApJ...363..234V)
*       4: North et al. (1994A&A...281..775N)
*       5: Cacciari et al. (1987A&A...183..314C)
*       6: Huenemoerder et al. (1984AJ.....89..851H)
*       7: De Boer et al. (1997A&A...317L..23)
*       8: Adelman & Philip (1990MNRAS.247..132A)
*       9: Hayes & Philip (1988PASP..100..801H)
*      11: Holweger & Rentzsch-Holm (1995A&A...303..819H)
*      12: Smalley & Dworetsky (1995A&A...293..446S)
*      13: Blackwell & Lynas-Gray (1994, Cat. <J/A+A/282/899>)
*      14: Alonso et al. (1996A&A...313..873A)
*      15: Malagnini & Morossi (1990A&AS...85.1015M)
*      16: Smalley (1993MNRAS.265.1035S)
*      17: Smalley & Dworetsky (1993A&A...271..515S)
*      18: Napiwotzki et al. (1993A&A...268..653N)
*      19: Adelman et al. (1980ApJS...43..491A)
*      20: Nissen (1981A&A....97..145N)
*      21: Van Veer-Menneret et al. (1989A&A...224..171V)
*      22: Hill & Landstreet (1993, Cat. <J/A+A/276/142>)
*      23: Gerbaldi et al. (1989A&AS...81..127G)
*      24: Adelman (1978ApJ...222..547A)
*Note (2): Sources for log(g):
*       6: Huenemoerder et al. (1984AJ.....89..851H)
*       9: Hayes & Philip (1988PASP..100..801H)
* nothing: this work
*Note (3): Camera and image number. The L and S following the image number
*           indicates the aperture used (Large or Small)

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

C  Declarations for 'table3.dat'	! Same as table1.dat but for the standard stars

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_2     ! (deg) Right Ascension J2000
      real*8        DEdeg_2     ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*9   Name_2      ! Star identification
      integer*4     Teff_2      ! (K) ? Effective temperature (this work)
      integer*4     Teff2_2     ! (K) ? Effective temperature from other work (1)
      integer*4     r_Teff2_2   ! ? Reference for Teff2
      integer*4     Teff3_2     ! (K) ? Effective temperature from other work (1)
      integer*4     r_Teff3_2   ! ? Reference for Teff3
      integer*4     Teff4_2     ! (K) ? Effective temperature from other work (1)
      integer*4     r_Teff4_2   ! ? Reference for Teff4
      real*4        log_g_1     ! ([cm/s2]) ? Surface gravity
      integer*4     r_log_g_1   ! ? Reference for log(g) (2)
      character*9   Images_2    ! IUE Images (SWP/LWP). (3)
      real*4        EW_1657_2800_1 ! (0.1nm) ? Equivalent width at 1657/2800A
      real*4        EW_1670_1   ! (0.1nm) ? Equivalent width at 1670A
      real*4        EW_1741_1   ! (0.1nm) ? Equivalent width at 1741A
      real*4        EW_1931_1   ! (0.1nm) ? Equivalent width at 1931A
      integer*4     RAh_2       ! (h) Right Ascension J2000 (hours)
      integer*4     RAm_2       ! (min) Right Ascension J2000 (minutes)
      real*4        RAs_2       ! (s) Right Ascension J2000 (seconds)
      character*1   DE__2       ! Declination J2000 (sign)
      integer*4     DEd_2       ! (deg) Declination J2000 (degrees)
      integer*4     DEm_2       ! (arcmin) Declination J2000 (minutes)
      integer*4     DEs_2       ! (arcsec) ?Declination J2000 (seconds)
*Note (1): Sources for T_eff_:
*       1: Iliev & Barzova (1995A&A...302..735I)
*       2: Stuerenburg (1993A&A...277..139S)
*       3: Venn & Lambert (1990ApJ...363..234V)
*       4: North et al. (1994A&A...281..775N)
*       5: Cacciari et al. (1987A&A...183..314C)
*       6: Huenemoerder et al. (1984AJ.....89..851H)
*       7: De Boer et al. (1997A&A...317L..23)
*       8: Adelman & Philip (1990MNRAS.247..132A)
*       9: Hayes & Philip (1988PASP..100..801H)
*      11: Holweger & Rentzsch-Holm (1995A&A...303..819H)
*      12: Smalley & Dworetsky (1995A&A...293..446S)
*      13: Blackwell & Lynas-Gray (1994, Cat. <J/A+A/282/899>)
*      14: Alonso et al. (1996A&A...313..873A)
*      15: Malagnini & Morossi (1990A&AS...85.1015M)
*      16: Smalley (1993MNRAS.265.1035S)
*      17: Smalley & Dworetsky (1993A&A...271..515S)
*      18: Napiwotzki et al. (1993A&A...268..653N)
*      19: Adelman et al. (1980ApJS...43..491A)
*      20: Nissen (1981A&A....97..145N)
*      21: Van Veer-Menneret et al. (1989A&A...224..171V)
*      22: Hill & Landstreet (1993, Cat. <J/A+A/276/142>)
*      23: Gerbaldi et al. (1989A&AS...81..127G)
*      24: Adelman (1978ApJ...222..547A)
*Note (2): Sources for log(g):
*       6: Huenemoerder et al. (1984AJ.....89..851H)
*       9: Hayes & Philip (1988PASP..100..801H)
* nothing: this work
*Note (3): Camera and image number. The L and S following the image number
*           indicates the aperture used (Large or Small)

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

C  Loading file 'table1.dat'	! Lambda Bootis stars included in Paper I
*                              (Pauzen et al., 1997, Cat. <J/A+AS/123/93>)
*                              with ultraviolet observations

C  Format for file interpretation

    1 format(
     +  A9,1X,I4,1X,I4,1X,I2,2X,I4,1X,I2,2X,I4,1X,I2,1X,F4.2,2X,I2,2X,
     +  A9,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,I2,1X,I2,1X,F4.1,1X,A1,
     +  I2,1X,I2,1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,41
        read(1,'(A102)')ar__
        read(ar__,1)
     +  Name,Teff,Teff2,r_Teff2,Teff3,r_Teff3,Teff4,r_Teff4,log_g,
     +  r_log_g,Images,EW_1657_2800,EW_1670,EW_1741,EW_1931,RAh,RAm,
     +  RAs,DE_,DEd,DEm,DEs
        if(ar__(11:14) .EQ. '') Teff = iNULL__
        if(ar__(16:19) .EQ. '') Teff2 = iNULL__
        if(ar__(21:22) .EQ. '') r_Teff2 = iNULL__
        if(ar__(25:28) .EQ. '') Teff3 = iNULL__
        if(ar__(30:31) .EQ. '') r_Teff3 = iNULL__
        if(ar__(34:37) .EQ. '') Teff4 = iNULL__
        if(ar__(39:40) .EQ. '') r_Teff4 = iNULL__
        if(ar__(42:45) .EQ. '') log_g = rNULL__
        if(ar__(48:49) .EQ. '') r_log_g = iNULL__
        if(ar__(62:66) .EQ. '') EW_1657_2800 = rNULL__
        if(ar__(68:71) .EQ. '') EW_1670 = rNULL__
        if(ar__(73:76) .EQ. '') EW_1741 = rNULL__
        if(ar__(78:81) .EQ. '') EW_1931 = rNULL__
        if(ar__(101:102) .EQ. '') DEs = 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(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  Name,Teff,Teff2,r_Teff2,Teff3,r_Teff3,Teff4,r_Teff4,log_g,
     +  r_log_g,Images,EW_1657_2800,EW_1670,EW_1741,EW_1931,RAh,RAm,
     +  RAs,DE_,DEd,DEm,DEs
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Same as table1.dat but for the field horizontal
*                              branch stars

C  Format for file interpretation

    2 format(
     +  A9,1X,I4,1X,I4,1X,I2,2X,I4,1X,I2,2X,I4,1X,I2,1X,F4.2,2X,I2,2X,
     +  A9,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,I2,1X,I2,1X,F4.1,1X,A1,
     +  I2,1X,I2,1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,12
        read(1,'(A102)')ar__1
        read(ar__1,2)
     +  Name_1,Teff_1,Teff2_1,r_Teff2_1,Teff3_1,r_Teff3_1,Teff4_1,
     +  r_Teff4_1,log_g_1,r_log_g_1,Images_1,EW_1657_2800_1,EW_1670_1,
     +  EW_1741_1,EW_1931_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1
        if(ar__1(11:14) .EQ. '') Teff_1 = iNULL__
        if(ar__1(16:19) .EQ. '') Teff2_1 = iNULL__
        if(ar__1(21:22) .EQ. '') r_Teff2_1 = iNULL__
        if(ar__1(25:28) .EQ. '') Teff3_1 = iNULL__
        if(ar__1(30:31) .EQ. '') r_Teff3_1 = iNULL__
        if(ar__1(34:37) .EQ. '') Teff4_1 = iNULL__
        if(ar__1(39:40) .EQ. '') r_Teff4_1 = iNULL__
        if(ar__1(42:45) .EQ. '') log_g_1 = rNULL__
        if(ar__1(48:49) .EQ. '') r_log_g_1 = iNULL__
        if(ar__1(62:66) .EQ. '') EW_1657_2800_1 = rNULL__
        if(ar__1(68:71) .EQ. '') EW_1670_1 = rNULL__
        if(ar__1(73:76) .EQ. '') EW_1741_1 = rNULL__
        if(ar__1(78:81) .EQ. '') EW_1931_1 = rNULL__
        if(ar__1(101:102) .EQ. '') DEs_1 = iNULL__
        RAdeg_1 = rNULL__
        DEdeg_1 = rNULL__
c  Derive coordinates RAdeg_1 and DEdeg_1 from input data
c  (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
        if(RAh_1 .GT. -180) RAdeg_1=RAh_1*15.
        if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
        if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
        if(DEd_1 .GE. 0) DEdeg_1=DEd_1
        if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
        if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
        if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c    ..............Just test output...........
        write(6,2)
     +  Name_1,Teff_1,Teff2_1,r_Teff2_1,Teff3_1,r_Teff3_1,Teff4_1,
     +  r_Teff4_1,log_g_1,r_log_g_1,Images_1,EW_1657_2800_1,EW_1670_1,
     +  EW_1741_1,EW_1931_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1
        write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Same as table1.dat but for the standard stars

C  Format for file interpretation

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

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,48
        read(1,'(A102)')ar__2
        read(ar__2,3)
     +  Name_2,Teff_2,Teff2_2,r_Teff2_2,Teff3_2,r_Teff3_2,Teff4_2,
     +  r_Teff4_2,log_g_1,r_log_g_1,Images_2,EW_1657_2800_1,EW_1670_1,
     +  EW_1741_1,EW_1931_1,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2
        if(ar__2(11:14) .EQ. '') Teff_2 = iNULL__
        if(ar__2(16:19) .EQ. '') Teff2_2 = iNULL__
        if(ar__2(21:22) .EQ. '') r_Teff2_2 = iNULL__
        if(ar__2(25:28) .EQ. '') Teff3_2 = iNULL__
        if(ar__2(30:31) .EQ. '') r_Teff3_2 = iNULL__
        if(ar__2(34:37) .EQ. '') Teff4_2 = iNULL__
        if(ar__2(39:40) .EQ. '') r_Teff4_2 = iNULL__
        if(ar__2(42:45) .EQ. '') log_g_1 = rNULL__
        if(ar__2(48:49) .EQ. '') r_log_g_1 = iNULL__
        if(ar__2(62:66) .EQ. '') EW_1657_2800_1 = rNULL__
        if(ar__2(68:71) .EQ. '') EW_1670_1 = rNULL__
        if(ar__2(73:76) .EQ. '') EW_1741_1 = rNULL__
        if(ar__2(78:81) .EQ. '') EW_1931_1 = rNULL__
        if(ar__2(101:102) .EQ. '') DEs_2 = iNULL__
        RAdeg_2 = rNULL__
        DEdeg_2 = rNULL__
c  Derive coordinates RAdeg_2 and DEdeg_2 from input data
c  (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
        if(RAh_2 .GT. -180) RAdeg_2=RAh_2*15.
        if(RAm_2 .GT. -180) RAdeg_2=RAdeg_2+RAm_2/4.
        if(RAs_2 .GT. -180) RAdeg_2=RAdeg_2+RAs_2/240.
        if(DEd_2 .GE. 0) DEdeg_2=DEd_2
        if(DEm_2 .GE. 0) DEdeg_2=DEdeg_2+DEm_2/60.
        if(DEs_2 .GE. 0) DEdeg_2=DEdeg_2+DEs_2/3600.
        if(DE__2.EQ.'-'.AND.DEdeg_2.GE.0) DEdeg_2=-DEdeg_2
c    ..............Just test output...........
        write(6,3)
     +  Name_2,Teff_2,Teff2_2,r_Teff2_2,Teff3_2,r_Teff3_2,Teff4_2,
     +  r_Teff4_2,log_g_1,r_log_g_1,Images_2,EW_1657_2800_1,EW_1670_1,
     +  EW_1741_1,EW_1931_1,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2
        write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_2
c    .......End.of.Just test output...........
      end do
      close(1)

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