Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/474/876 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/474/876 WASP-80b wavelength-binned light curves (Kirk+, 2018)
*================================================================================
*LRG-BEASTS III: ground-based transmission spectrum of the gas giant orbiting
*the cool dwarf WASP-80.
* Kirk J., Wheatley P.J., Louden T., Skillen I., King G.W., Mccormac J.,
* Irwin P.G.J.
* <Mon. Not. R. Astron. Soc., 474, 876-885 (2018)>
* =2018MNRAS.474..876K (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Reduced light curves from night 1 with
time and airmass
integer*4 nr__
parameter (nr__=230) ! Number of records
character*579 ar__ ! Full-size record
real*8 MJD ! (d) Modified Julian Date
real*4 Airmass ! Mean effective airmass
real*8 FWl ! Normalised white light flux
real*8 e_FWl ! Normalised white light error
real*8 F4790A ! Normalised flux in 300A-wide bin
* centred on 4790{AA}
real*8 e_F4790A ! Normalised error in 300A-wide bin
* centred on 4790{AA}
real*8 F5090A ! Normalised flux in 300A-wide bin
* centred on 5090{AA}
real*8 e_F5090A ! Normalised error in 300A-wide bin
* centred on 5090{AA}
real*8 F5340A ! Normalised flux in 200A-wide bin
* centred on 5340{AA}
real*8 e_F5340A ! Normalised error in 200A-wide bin
* centred on 5340{AA}
real*8 F5540A ! Normalised flux in 200A-wide bin
* centred on 5540{AA}
real*8 e_F5540A ! Normalised error in 200A-wide bin
* centred on 5540{AA}
real*8 F5690A ! Normalised flux in 100A-wide bin
* centred on 5690{AA}
real*8 e_F5690A ! Normalised error in 100A-wide bin
* centred on 5690{AA}
real*8 F5790A ! Normalised flux in 100A-wide bin
* centred on 5790{AA}
real*8 e_F5790A ! Normalised error in 100A-wide bin
* centred on 5790{AA}
real*8 F5890A ! Normalised flux in 100A-wide bin
* centred on 5890{AA}
real*8 e_F5890A ! Normalised error in 100A-wide bin
* centred on 5890{AA}
real*8 F5990A ! Normalised flux in 100A-wide bin
* centred on 5990{AA}
real*8 e_F5990A ! Normalised error in 100A-wide bin
* centred on 5990{AA}
real*8 F6090A ! Normalised flux in 100A-wide bin
* centred on 6090{AA}
real*8 e_F6090A ! Normalised error in 100A-wide bin
* centred on 6090{AA}
real*8 F6190A ! Normalised flux in 100A-wide bin
* centred on 6190{AA}
real*8 e_F6190A ! Normalised error in 100A-wide bin
* centred on 6190{AA}
real*8 F6290A ! Normalised flux in 100A-wide bin
* centred on 6290{AA}
real*8 e_F6290A ! Normalised error in 100A-wide bin
* centred on 6290{AA}
real*8 F6390A ! Normalised flux in 100A-wide bin
* centred on 6390{AA}
real*8 e_F6390A ! Normalised error in 100A-wide bin
* centred on 6390{AA}
real*8 F6490A ! Normalised flux in 100A-wide bin
* centred on 6490{AA}
real*8 e_F6490A ! Normalised error in 100A-wide bin
* centred on 6490{AA}
real*8 F6590A ! Normalised flux in 100A-wide bin
* centred on 6590{AA}
real*8 e_F6590A ! Normalised error in 100A-wide bin
* centred on 6590{AA}
real*8 F6690A ! Normalised flux in 100A-wide bin
* centred on 6690{AA}
real*8 e_F6690A ! Normalised error in 100A-wide bin
* centred on 6690{AA}
real*8 F6790A ! Normalised flux in 100A-wide bin
* centred on 6790{AA}
real*8 e_F6790A ! Normalised error in 100A-wide bin
* centred on 6790{AA}
real*8 F6890A ! Normalised flux in 100A-wide bin
* centred on 6890{AA}
real*8 e_F6890A ! Normalised error in 100A-wide bin
* centred on 6890{AA}
real*8 F6990A ! Normalised flux in 100A-wide bin
* centred on 6990{AA}
real*8 e_F6990A ! Normalised error in 100A-wide bin
* centred on 6990{AA}
real*8 F7090A ! Normalised flux in 100A-wide bin
* centred on 7090{AA}
real*8 e_F7090A ! Normalised error in 100A-wide bin
* centred on 7090{AA}
real*8 F7190A ! Normalised flux in 100A-wide bin
* centred on 7190{AA}
real*8 e_F7190A ! Normalised error in 100A-wide bin
* centred on 7190{AA}
real*8 F7290A ! Normalised flux in 100A-wide bin
* centred on 7290{AA}
real*8 e_F7290A ! Normalised error in 100A-wide bin
* centred on 7290{AA}
real*8 F7390A ! Normalised flux in 100A-wide bin
* centred on 7390{AA}
real*8 e_F7390A ! Normalised error in 100A-wide bin
* centred on 7390{AA}
real*8 F7490A ! Normalised flux in 100A-wide bin
* centred on 7490{AA}
real*8 e_F7490A ! Normalised error in 100A-wide bin
* centred on 7490{AA}
real*8 F7790A ! Normalised flux in 100A-wide bin
* centred on 7790{AA}
real*8 e_F7790A ! Normalised error in 100A-wide bin
* centred on 7790{AA}
real*8 F7890A ! Normalised flux in 100A-wide bin
* centred on 7890{AA}
real*8 e_F7890A ! Normalised error in 100A-wide bin
* centred on 7890{AA}
real*8 F7990A ! Normalised flux in 100A-wide bin
* centred on 7990{AA}
real*8 e_F7990A ! Normalised error in 100A-wide bin
* centred on 7990{AA}
real*8 F8090A ! Normalised flux in 100A-wide bin
* centred on 8090{AA}
real*8 e_F8090A ! Normalised error in 100A-wide bin
* centred on 8090{AA}
real*8 F8190A ! Normalised flux in 100A-wide bin
* centred on 8190{AA}
real*8 e_F8190A ! Normalised error in 100A-wide bin
* centred on 8190{AA}
real*8 F8290A ! Normalised flux in 100A-wide bin
* centred on 8290{AA}
real*8 e_F8290A ! Normalised error in 100A-wide bin
* centred on 8290{AA}
real*8 F8390A ! Normalised flux in 100A-wide bin
* centred on 8390{AA}
real*8 e_F8390A ! Normalised error in 100A-wide bin
* centred on 8390{AA}
real*8 F8490A ! Normalised flux in 100A-wide bin
* centred on 8490{AA}
real*8 e_F8490A ! Normalised error in 100A-wide bin
* centred on 8490{AA}
real*8 F8590A ! Normalised flux in 100A-wide bin
* centred on 8590{AA}
real*8 e_F8590A ! Normalised error in 100A-wide bin
* centred on 8590{AA}
real*8 F8690A ! Normalised flux in 100A-wide bin
* centred on 8690{AA}
real*8 e_F8690A ! Normalised error in 100A-wide bin
* centred on 8690{AA}
real*8 F8790A ! Normalised flux in 100A-wide bin
* centred on 8790{AA}
real*8 e_F8790A ! Normalised error in 100A-wide bin
* centred on 8790{AA}
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Reduced light curves from night 2 with
time and airmass
integer*4 nr__1
parameter (nr__1=173) ! Number of records
character*579 ar__1 ! Full-size record
real*8 MJD_1 ! (d) Modified Julian Date
real*4 Airmass_1 ! Mean effective airmass
real*8 FWl_1 ! Normalised white light flux
real*8 e_FWl_1 ! Normalised white light error
real*8 F4790A_1 ! Normalised flux in 300A-wide bin
* centred on 4790{AA}
real*8 e_F4790A_1 ! Normalised error in 300A-wide bin
* centred on 4790{AA}
real*8 F5090A_1 ! Normalised flux in 300A-wide bin
* centred on 5090{AA}
real*8 e_F5090A_1 ! Normalised error in 300A-wide bin
* centred on 5090{AA}
real*8 F5340A_1 ! Normalised flux in 200A-wide bin
* centred on 5340{AA}
real*8 e_F5340A_1 ! Normalised error in 200A-wide bin
* centred on 5340{AA}
real*8 F5540A_1 ! Normalised flux in 200A-wide bin
* centred on 5540{AA}
real*8 e_F5540A_1 ! Normalised error in 200A-wide bin
* centred on 5540{AA}
real*8 F5690A_1 ! Normalised flux in 100A-wide bin
* centred on 5690{AA}
real*8 e_F5690A_1 ! Normalised error in 100A-wide bin
* centred on 5690{AA}
real*8 F5790A_1 ! Normalised flux in 100A-wide bin
* centred on 5790{AA}
real*8 e_F5790A_1 ! Normalised error in 100A-wide bin
* centred on 5790{AA}
real*8 F5890A_1 ! Normalised flux in 100A-wide bin
* centred on 5890{AA}
real*8 e_F5890A_1 ! Normalised error in 100A-wide bin
* centred on 5890{AA}
real*8 F5990A_1 ! Normalised flux in 100A-wide bin
* centred on 5990{AA}
real*8 e_F5990A_1 ! Normalised error in 100A-wide bin
* centred on 5990{AA}
real*8 F6090A_1 ! Normalised flux in 100A-wide bin
* centred on 6090{AA}
real*8 e_F6090A_1 ! Normalised error in 100A-wide bin
* centred on 6090{AA}
real*8 F6190A_1 ! Normalised flux in 100A-wide bin
* centred on 6190{AA}
real*8 e_F6190A_1 ! Normalised error in 100A-wide bin
* centred on 6190{AA}
real*8 F6290A_1 ! Normalised flux in 100A-wide bin
* centred on 6290{AA}
real*8 e_F6290A_1 ! Normalised error in 100A-wide bin
* centred on 6290{AA}
real*8 F6390A_1 ! Normalised flux in 100A-wide bin
* centred on 6390{AA}
real*8 e_F6390A_1 ! Normalised error in 100A-wide bin
* centred on 6390{AA}
real*8 F6490A_1 ! Normalised flux in 100A-wide bin
* centred on 6490{AA}
real*8 e_F6490A_1 ! Normalised error in 100A-wide bin
* centred on 6490{AA}
real*8 F6590A_1 ! Normalised flux in 100A-wide bin
* centred on 6590{AA}
real*8 e_F6590A_1 ! Normalised error in 100A-wide bin
* centred on 6590{AA}
real*8 F6690A_1 ! Normalised flux in 100A-wide bin
* centred on 6690{AA}
real*8 e_F6690A_1 ! Normalised error in 100A-wide bin
* centred on 6690{AA}
real*8 F6790A_1 ! Normalised flux in 100A-wide bin
* centred on 6790{AA}
real*8 e_F6790A_1 ! Normalised error in 100A-wide bin
* centred on 6790{AA}
real*8 F6890A_1 ! Normalised flux in 100A-wide bin
* centred on 6890{AA}
real*8 e_F6890A_1 ! Normalised error in 100A-wide bin
* centred on 6890{AA}
real*8 F6990A_1 ! Normalised flux in 100A-wide bin
* centred on 6990{AA}
real*8 e_F6990A_1 ! Normalised error in 100A-wide bin
* centred on 6990{AA}
real*8 F7090A_1 ! Normalised flux in 100A-wide bin
* centred on 7090{AA}
real*8 e_F7090A_1 ! Normalised error in 100A-wide bin
* centred on 7090{AA}
real*8 F7190A_1 ! Normalised flux in 100A-wide bin
* centred on 7190{AA}
real*8 e_F7190A_1 ! Normalised error in 100A-wide bin
* centred on 7190{AA}
real*8 F7290A_1 ! Normalised flux in 100A-wide bin
* centred on 7290{AA}
real*8 e_F7290A_1 ! Normalised error in 100A-wide bin
* centred on 7290{AA}
real*8 F7390A_1 ! Normalised flux in 100A-wide bin
* centred on 7390{AA}
real*8 e_F7390A_1 ! Normalised error in 100A-wide bin
* centred on 7390{AA}
real*8 F7490A_1 ! Normalised flux in 100A-wide bin
* centred on 7490{AA}
real*8 e_F7490A_1 ! Normalised error in 100A-wide bin
* centred on 7490{AA}
real*8 F7790A_1 ! Normalised flux in 100A-wide bin
* centred on 7790{AA}
real*8 e_F7790A_1 ! Normalised error in 100A-wide bin
* centred on 7790{AA}
real*8 F7890A_1 ! Normalised flux in 100A-wide bin
* centred on 7890{AA}
real*8 e_F7890A_1 ! Normalised error in 100A-wide bin
* centred on 7890{AA}
real*8 F7990A_1 ! Normalised flux in 100A-wide bin
* centred on 7990{AA}
real*8 e_F7990A_1 ! Normalised error in 100A-wide bin
* centred on 7990{AA}
real*8 F8090A_1 ! Normalised flux in 100A-wide bin
* centred on 8090{AA}
real*8 e_F8090A_1 ! Normalised error in 100A-wide bin
* centred on 8090{AA}
real*8 F8190A_1 ! Normalised flux in 100A-wide bin
* centred on 8190{AA}
real*8 e_F8190A_1 ! Normalised error in 100A-wide bin
* centred on 8190{AA}
real*8 F8290A_1 ! Normalised flux in 100A-wide bin
* centred on 8290{AA}
real*8 e_F8290A_1 ! Normalised error in 100A-wide bin
* centred on 8290{AA}
real*8 F8390A_1 ! Normalised flux in 100A-wide bin
* centred on 8390{AA}
real*8 e_F8390A_1 ! Normalised error in 100A-wide bin
* centred on 8390{AA}
real*8 F8490A_1 ! Normalised flux in 100A-wide bin
* centred on 8490{AA}
real*8 e_F8490A_1 ! Normalised error in 100A-wide bin
* centred on 8490{AA}
real*8 F8590A_1 ! Normalised flux in 100A-wide bin
* centred on 8590{AA}
real*8 e_F8590A_1 ! Normalised error in 100A-wide bin
* centred on 8590{AA}
real*8 F8690A_1 ! Normalised flux in 100A-wide bin
* centred on 8690{AA}
real*8 e_F8690A_1 ! Normalised error in 100A-wide bin
* centred on 8690{AA}
real*8 F8790A_1 ! Normalised flux in 100A-wide bin
* centred on 8790{AA}
real*8 e_F8790A_1 ! Normalised error in 100A-wide bin
* centred on 8790{AA}
C=============================================================================
C Loading file 'table1.dat' ! Reduced light curves from night 1 with
* time and airmass
C Format for file interpretation
1 format(
+ F12.6,1X,F6.4,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,
+ F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,
+ F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,
+ F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,
+ F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,230
read(1,'(A579)')ar__
read(ar__,1)
+ MJD,Airmass,FWl,e_FWl,F4790A,e_F4790A,F5090A,e_F5090A,F5340A,
+ e_F5340A,F5540A,e_F5540A,F5690A,e_F5690A,F5790A,e_F5790A,
+ F5890A,e_F5890A,F5990A,e_F5990A,F6090A,e_F6090A,F6190A,
+ e_F6190A,F6290A,e_F6290A,F6390A,e_F6390A,F6490A,e_F6490A,
+ F6590A,e_F6590A,F6690A,e_F6690A,F6790A,e_F6790A,F6890A,
+ e_F6890A,F6990A,e_F6990A,F7090A,e_F7090A,F7190A,e_F7190A,
+ F7290A,e_F7290A,F7390A,e_F7390A,F7490A,e_F7490A,F7790A,
+ e_F7790A,F7890A,e_F7890A,F7990A,e_F7990A,F8090A,e_F8090A,
+ F8190A,e_F8190A,F8290A,e_F8290A,F8390A,e_F8390A,F8490A,
+ e_F8490A,F8590A,e_F8590A,F8690A,e_F8690A,F8790A,e_F8790A
c ..............Just test output...........
write(6,1)
+ MJD,Airmass,FWl,e_FWl,F4790A,e_F4790A,F5090A,e_F5090A,F5340A,
+ e_F5340A,F5540A,e_F5540A,F5690A,e_F5690A,F5790A,e_F5790A,
+ F5890A,e_F5890A,F5990A,e_F5990A,F6090A,e_F6090A,F6190A,
+ e_F6190A,F6290A,e_F6290A,F6390A,e_F6390A,F6490A,e_F6490A,
+ F6590A,e_F6590A,F6690A,e_F6690A,F6790A,e_F6790A,F6890A,
+ e_F6890A,F6990A,e_F6990A,F7090A,e_F7090A,F7190A,e_F7190A,
+ F7290A,e_F7290A,F7390A,e_F7390A,F7490A,e_F7490A,F7790A,
+ e_F7790A,F7890A,e_F7890A,F7990A,e_F7990A,F8090A,e_F8090A,
+ F8190A,e_F8190A,F8290A,e_F8290A,F8390A,e_F8390A,F8490A,
+ e_F8490A,F8590A,e_F8590A,F8690A,e_F8690A,F8790A,e_F8790A
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Reduced light curves from night 2 with
* time and airmass
C Format for file interpretation
2 format(
+ F12.6,1X,F6.4,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,
+ F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,
+ F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,
+ F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,
+ F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,
+ 1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,173
read(1,'(A579)')ar__1
read(ar__1,2)
+ MJD_1,Airmass_1,FWl_1,e_FWl_1,F4790A_1,e_F4790A_1,F5090A_1,
+ e_F5090A_1,F5340A_1,e_F5340A_1,F5540A_1,e_F5540A_1,F5690A_1,
+ e_F5690A_1,F5790A_1,e_F5790A_1,F5890A_1,e_F5890A_1,F5990A_1,
+ e_F5990A_1,F6090A_1,e_F6090A_1,F6190A_1,e_F6190A_1,F6290A_1,
+ e_F6290A_1,F6390A_1,e_F6390A_1,F6490A_1,e_F6490A_1,F6590A_1,
+ e_F6590A_1,F6690A_1,e_F6690A_1,F6790A_1,e_F6790A_1,F6890A_1,
+ e_F6890A_1,F6990A_1,e_F6990A_1,F7090A_1,e_F7090A_1,F7190A_1,
+ e_F7190A_1,F7290A_1,e_F7290A_1,F7390A_1,e_F7390A_1,F7490A_1,
+ e_F7490A_1,F7790A_1,e_F7790A_1,F7890A_1,e_F7890A_1,F7990A_1,
+ e_F7990A_1,F8090A_1,e_F8090A_1,F8190A_1,e_F8190A_1,F8290A_1,
+ e_F8290A_1,F8390A_1,e_F8390A_1,F8490A_1,e_F8490A_1,F8590A_1,
+ e_F8590A_1,F8690A_1,e_F8690A_1,F8790A_1,e_F8790A_1
c ..............Just test output...........
write(6,2)
+ MJD_1,Airmass_1,FWl_1,e_FWl_1,F4790A_1,e_F4790A_1,F5090A_1,
+ e_F5090A_1,F5340A_1,e_F5340A_1,F5540A_1,e_F5540A_1,F5690A_1,
+ e_F5690A_1,F5790A_1,e_F5790A_1,F5890A_1,e_F5890A_1,F5990A_1,
+ e_F5990A_1,F6090A_1,e_F6090A_1,F6190A_1,e_F6190A_1,F6290A_1,
+ e_F6290A_1,F6390A_1,e_F6390A_1,F6490A_1,e_F6490A_1,F6590A_1,
+ e_F6590A_1,F6690A_1,e_F6690A_1,F6790A_1,e_F6790A_1,F6890A_1,
+ e_F6890A_1,F6990A_1,e_F6990A_1,F7090A_1,e_F7090A_1,F7190A_1,
+ e_F7190A_1,F7290A_1,e_F7290A_1,F7390A_1,e_F7390A_1,F7490A_1,
+ e_F7490A_1,F7790A_1,e_F7790A_1,F7890A_1,e_F7890A_1,F7990A_1,
+ e_F7990A_1,F8090A_1,e_F8090A_1,F8190A_1,e_F8190A_1,F8290A_1,
+ e_F8290A_1,F8390A_1,e_F8390A_1,F8490A_1,e_F8490A_1,F8590A_1,
+ e_F8590A_1,F8690A_1,e_F8690A_1,F8790A_1,e_F8790A_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end