Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/610/L12 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-Apr-20 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/610/L12 KIC 8462852 GTC spectra (Deeg+, 2018) *================================================================================ *Non-grey dimming events of KIC 8462852 from GTC spectrophotometry. * Deeg H.J., Alonso R., Nespral D., Boyajian T. * <Astron. Astrophys. 610, L12 (2018)> * =2018A&A...610L..12D (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Observed fluxes at each pointing in five wavelengths bands, indicated by their central wavelengths in nm integer*4 nr__ parameter (nr__=16) ! Number of records character*97 ar__ ! Full-size record integer*4 Nr (nr__) ! [1/16] Pointing sequential number character*1 n_Nr (nr__) ! [ab] Note (1) real*8 BJD (nr__) ! (d) Observation central time of each pointing * (BJD-2400000) real*4 F556 (nr__) ! Observed normalized flux in 556um band real*4 e_F556 (nr__) ! rms uncertainty on F556 real*4 F618 (nr__) ! Observed normalized flux in 618um band real*4 e_F618 (nr__) ! rms uncertainty on F618 real*4 F709 (nr__) ! Observed normalized flux in 709um band real*4 e_F709 (nr__) ! rms uncertainty on F709 real*4 F792 (nr__) ! Observed normalized flux in 792um band real*4 e_F792 (nr__) ! rms uncertainty on F792 real*4 F868 (nr__) ! Observed normalized flux in 868um band real*4 e_F868 (nr__) ! rms uncertainty on F868 character*12 FileName (nr__) ! Name of the table with times of mid-exposure * of each spectrum for this pointing *Note (1): Notes as follows: * a = Data not used in analysis due to large scatter from varying exposure times * b = Data not used due to poor sky transparency c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p01.dat' ! Pointing 01 times of mid-exposure of each spectrum integer*4 nr__1 parameter (nr__1=10) ! Number of records character*51 ar__1 ! Full-size record real*8 BJD_1 (nr__1) ! (d) Time (BJD-2400000) character*12 spTarget (nr__1) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp (nr__1) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p02.dat' ! Pointing 02 times of mid-exposure of each spectrum integer*4 nr__2 parameter (nr__2=45) ! Number of records character*51 ar__2 ! Full-size record real*8 BJD_2 (nr__2) ! (d) Time (BJD-2400000) character*12 spTarget_1 (nr__2) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_1 (nr__2) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p03.dat' ! Pointing 03 times of mid-exposure of each spectrum integer*4 nr__3 parameter (nr__3=45) ! Number of records character*51 ar__3 ! Full-size record real*8 BJD_3 (nr__3) ! (d) Time (BJD-2400000) character*12 spTarget_2 (nr__3) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_2 (nr__3) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p04.dat' ! Pointing 04 times of mid-exposure of each spectrum integer*4 nr__4 parameter (nr__4=33) ! Number of records character*51 ar__4 ! Full-size record real*8 BJD_4 (nr__4) ! (d) Time (BJD-2400000) character*12 spTarget_3 (nr__4) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_3 (nr__4) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p05.dat' ! Pointing 05 times of mid-exposure of each spectrum integer*4 nr__5 parameter (nr__5=30) ! Number of records character*51 ar__5 ! Full-size record real*8 BJD_5 (nr__5) ! (d) Time (BJD-2400000) character*12 spTarget_4 (nr__5) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_4 (nr__5) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p06.dat' ! Pointing 06 times of mid-exposure of each spectrum integer*4 nr__6 parameter (nr__6=35) ! Number of records character*51 ar__6 ! Full-size record real*8 BJD_6 (nr__6) ! (d) Time (BJD-2400000) character*12 spTarget_5 (nr__6) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_5 (nr__6) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p07.dat' ! Pointing 07 times of mid-exposure of each spectrum integer*4 nr__7 parameter (nr__7=30) ! Number of records character*51 ar__7 ! Full-size record real*8 BJD_7 (nr__7) ! (d) Time (BJD-2400000) character*12 spTarget_6 (nr__7) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_6 (nr__7) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p08.dat' ! Pointing 08 times of mid-exposure of each spectrum integer*4 nr__8 parameter (nr__8=31) ! Number of records character*51 ar__8 ! Full-size record real*8 BJD_8 (nr__8) ! (d) Time (BJD-2400000) character*12 spTarget_7 (nr__8) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_7 (nr__8) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p09.dat' ! Pointing 09 times of mid-exposure of each spectrum integer*4 nr__9 parameter (nr__9=44) ! Number of records character*51 ar__9 ! Full-size record real*8 BJD_9 (nr__9) ! (d) Time (BJD-2400000) character*12 spTarget_8 (nr__9) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_8 (nr__9) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p10.dat' ! Pointing 10 times of mid-exposure of each spectrum integer*4 nr__10 parameter (nr__10=31) ! Number of records character*51 ar__10 ! Full-size record real*8 BJD_10 (nr__10) ! (d) Time (BJD-2400000) character*12 spTarget_9 (nr__10) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_9 (nr__10) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p11.dat' ! Pointing 11 times of mid-exposure of each spectrum integer*4 nr__11 parameter (nr__11=19) ! Number of records character*51 ar__11 ! Full-size record real*8 BJD_11 (nr__11) ! (d) Time (BJD-2400000) character*12 spTarget_10(nr__11) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_10 (nr__11) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p12.dat' ! Pointing 12 times of mid-exposure of each spectrum integer*4 nr__12 parameter (nr__12=24) ! Number of records character*51 ar__12 ! Full-size record real*8 BJD_12 (nr__12) ! (d) Time (BJD-2400000) character*12 spTarget_11(nr__12) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_11 (nr__12) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p13.dat' ! Pointing 13 times of mid-exposure of each spectrum integer*4 nr__13 parameter (nr__13=30) ! Number of records character*51 ar__13 ! Full-size record real*8 BJD_13 (nr__13) ! (d) Time (BJD-2400000) character*12 spTarget_12(nr__13) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_12 (nr__13) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p14.dat' ! Pointing 14 times of mid-exposure of each spectrum integer*4 nr__14 parameter (nr__14=35) ! Number of records character*51 ar__14 ! Full-size record real*8 BJD_14 (nr__14) ! (d) Time (BJD-2400000) character*12 spTarget_13(nr__14) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_13 (nr__14) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p15.dat' ! Pointing 15 times of mid-exposure of each spectrum integer*4 nr__15 parameter (nr__15=31) ! Number of records character*51 ar__15 ! Full-size record real*8 BJD_15 (nr__15) ! (d) Time (BJD-2400000) character*12 spTarget_14(nr__15) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_14 (nr__15) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'time_p16.dat' ! Pointing 16 times of mid-exposure of each spectrum integer*4 nr__16 parameter (nr__16=43) ! Number of records character*51 ar__16 ! Full-size record real*8 BJD_16 (nr__16) ! (d) Time (BJD-2400000) character*12 spTarget_15(nr__16) ! Name of the target spectrum name in * subdirectory sp_target character*12 spComp_15 (nr__16) ! Name of the comparison spectrum name in * subdirectory sp_compar c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'sp_target/*' ! Individual target (KIC 8462852) spectra integer*4 nr__17 parameter (nr__17=516) ! Number of records character*30 ar__17 ! Full-size record real*8 lambda (nr__17) ! (nm) Wavelength real*8 Flux (nr__17) ! Flux in arbitrary unit c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'sp_compar/*' ! Individual comparison (KIC 8462763) spectra integer*4 nr__18 parameter (nr__18=516) ! Number of records character*30 ar__18 ! Full-size record real*8 lambda_1 (nr__18) ! (nm) Wavelength real*8 Flux_1 (nr__18) ! Flux in arbitrary unit C============================================================================= C Loading file 'table1.dat' ! Observed fluxes at each pointing in five * wavelengths bands, indicated by their central * wavelengths in nm C Format for file interpretation 1 format( + I2,A1,1X,F10.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F6.4,1X, + F6.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,A12) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,16 read(1,'(A97)')ar__ read(ar__,1) + Nr(i__),n_Nr(i__),BJD(i__),F556(i__),e_F556(i__),F618(i__), + e_F618(i__),F709(i__),e_F709(i__),F792(i__),e_F792(i__), + F868(i__),e_F868(i__),FileName(i__) c ..............Just test output........... write(6,1) + Nr(i__),n_Nr(i__),BJD(i__),F556(i__),e_F556(i__),F618(i__), + e_F618(i__),F709(i__),e_F709(i__),F792(i__),e_F792(i__), + F868(i__),e_F868(i__),FileName(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p01.dat' ! Pointing 01 times of mid-exposure of each spectrum C Format for file interpretation 2 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p01.dat') write(6,*) '....Loading file: time_p01.dat' do i__=1,10 read(1,'(A51)')ar__1 read(ar__1,2)BJD_1(i__),spTarget(i__),spComp(i__) c ..............Just test output........... write(6,2)BJD_1(i__),spTarget(i__),spComp(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p02.dat' ! Pointing 02 times of mid-exposure of each spectrum C Format for file interpretation 3 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p02.dat') write(6,*) '....Loading file: time_p02.dat' do i__=1,45 read(1,'(A51)')ar__2 read(ar__2,3)BJD_2(i__),spTarget_1(i__),spComp_1(i__) c ..............Just test output........... write(6,3)BJD_2(i__),spTarget_1(i__),spComp_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p03.dat' ! Pointing 03 times of mid-exposure of each spectrum C Format for file interpretation 4 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p03.dat') write(6,*) '....Loading file: time_p03.dat' do i__=1,45 read(1,'(A51)')ar__3 read(ar__3,4)BJD_3(i__),spTarget_2(i__),spComp_2(i__) c ..............Just test output........... write(6,4)BJD_3(i__),spTarget_2(i__),spComp_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p04.dat' ! Pointing 04 times of mid-exposure of each spectrum C Format for file interpretation 5 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p04.dat') write(6,*) '....Loading file: time_p04.dat' do i__=1,33 read(1,'(A51)')ar__4 read(ar__4,5)BJD_4(i__),spTarget_3(i__),spComp_3(i__) c ..............Just test output........... write(6,5)BJD_4(i__),spTarget_3(i__),spComp_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p05.dat' ! Pointing 05 times of mid-exposure of each spectrum C Format for file interpretation 6 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p05.dat') write(6,*) '....Loading file: time_p05.dat' do i__=1,30 read(1,'(A51)')ar__5 read(ar__5,6)BJD_5(i__),spTarget_4(i__),spComp_4(i__) c ..............Just test output........... write(6,6)BJD_5(i__),spTarget_4(i__),spComp_4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p06.dat' ! Pointing 06 times of mid-exposure of each spectrum C Format for file interpretation 7 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p06.dat') write(6,*) '....Loading file: time_p06.dat' do i__=1,35 read(1,'(A51)')ar__6 read(ar__6,7)BJD_6(i__),spTarget_5(i__),spComp_5(i__) c ..............Just test output........... write(6,7)BJD_6(i__),spTarget_5(i__),spComp_5(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p07.dat' ! Pointing 07 times of mid-exposure of each spectrum C Format for file interpretation 8 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p07.dat') write(6,*) '....Loading file: time_p07.dat' do i__=1,30 read(1,'(A51)')ar__7 read(ar__7,8)BJD_7(i__),spTarget_6(i__),spComp_6(i__) c ..............Just test output........... write(6,8)BJD_7(i__),spTarget_6(i__),spComp_6(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p08.dat' ! Pointing 08 times of mid-exposure of each spectrum C Format for file interpretation 9 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p08.dat') write(6,*) '....Loading file: time_p08.dat' do i__=1,31 read(1,'(A51)')ar__8 read(ar__8,9)BJD_8(i__),spTarget_7(i__),spComp_7(i__) c ..............Just test output........... write(6,9)BJD_8(i__),spTarget_7(i__),spComp_7(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p09.dat' ! Pointing 09 times of mid-exposure of each spectrum C Format for file interpretation 10 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p09.dat') write(6,*) '....Loading file: time_p09.dat' do i__=1,44 read(1,'(A51)')ar__9 read(ar__9,10)BJD_9(i__),spTarget_8(i__),spComp_8(i__) c ..............Just test output........... write(6,10)BJD_9(i__),spTarget_8(i__),spComp_8(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p10.dat' ! Pointing 10 times of mid-exposure of each spectrum C Format for file interpretation 11 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p10.dat') write(6,*) '....Loading file: time_p10.dat' do i__=1,31 read(1,'(A51)')ar__10 read(ar__10,11)BJD_10(i__),spTarget_9(i__),spComp_9(i__) c ..............Just test output........... write(6,11)BJD_10(i__),spTarget_9(i__),spComp_9(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p11.dat' ! Pointing 11 times of mid-exposure of each spectrum C Format for file interpretation 12 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p11.dat') write(6,*) '....Loading file: time_p11.dat' do i__=1,19 read(1,'(A51)')ar__11 read(ar__11,12)BJD_11(i__),spTarget_10(i__),spComp_10(i__) c ..............Just test output........... write(6,12)BJD_11(i__),spTarget_10(i__),spComp_10(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p12.dat' ! Pointing 12 times of mid-exposure of each spectrum C Format for file interpretation 13 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p12.dat') write(6,*) '....Loading file: time_p12.dat' do i__=1,24 read(1,'(A51)')ar__12 read(ar__12,13)BJD_12(i__),spTarget_11(i__),spComp_11(i__) c ..............Just test output........... write(6,13)BJD_12(i__),spTarget_11(i__),spComp_11(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p13.dat' ! Pointing 13 times of mid-exposure of each spectrum C Format for file interpretation 14 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p13.dat') write(6,*) '....Loading file: time_p13.dat' do i__=1,30 read(1,'(A51)')ar__13 read(ar__13,14)BJD_13(i__),spTarget_12(i__),spComp_12(i__) c ..............Just test output........... write(6,14)BJD_13(i__),spTarget_12(i__),spComp_12(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p14.dat' ! Pointing 14 times of mid-exposure of each spectrum C Format for file interpretation 15 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p14.dat') write(6,*) '....Loading file: time_p14.dat' do i__=1,35 read(1,'(A51)')ar__14 read(ar__14,15)BJD_14(i__),spTarget_13(i__),spComp_13(i__) c ..............Just test output........... write(6,15)BJD_14(i__),spTarget_13(i__),spComp_13(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p15.dat' ! Pointing 15 times of mid-exposure of each spectrum C Format for file interpretation 16 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p15.dat') write(6,*) '....Loading file: time_p15.dat' do i__=1,31 read(1,'(A51)')ar__15 read(ar__15,16)BJD_15(i__),spTarget_14(i__),spComp_14(i__) c ..............Just test output........... write(6,16)BJD_15(i__),spTarget_14(i__),spComp_14(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'time_p16.dat' ! Pointing 16 times of mid-exposure of each spectrum C Format for file interpretation 17 format(3X,F22.16,1X,A12,1X,A12) C Effective file loading open(unit=1,status='old',file= +'time_p16.dat') write(6,*) '....Loading file: time_p16.dat' do i__=1,43 read(1,'(A51)')ar__16 read(ar__16,17)BJD_16(i__),spTarget_15(i__),spComp_15(i__) c ..............Just test output........... write(6,17)BJD_16(i__),spTarget_15(i__),spComp_15(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'sp_target/*' ! Individual target (KIC 8462852) spectra C Format for file interpretation 18 format(2X,F13.7,1X,F14.7) C Effective file loading open(unit=1,status='old',file= +'sp_target/*') write(6,*) '....Loading file: sp_target/*' do i__=1,516 read(1,'(A30)')ar__17 read(ar__17,18)lambda(i__),Flux(i__) c ..............Just test output........... write(6,18)lambda(i__),Flux(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'sp_compar/*' ! Individual comparison (KIC 8462763) spectra C Format for file interpretation 19 format(2X,F13.7,1X,F14.7) C Effective file loading open(unit=1,status='old',file= +'sp_compar/*') write(6,*) '....Loading file: sp_compar/*' do i__=1,516 read(1,'(A30)')ar__18 read(ar__18,19)lambda_1(i__),Flux_1(i__) c ..............Just test output........... write(6,19)lambda_1(i__),Flux_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end