Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/715/671 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-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/ApJ/715/671 X-ray and IR emission from YSOs (Wolk+, 2010) *================================================================================ *X-ray and infrared emission from young stellar objects near LkH{alpha} 101. * Wolk S.J., Winston E., Bourke T.L., Gutermuth R., Megeath S.T., * Spitzbart B.D., Osten R. * <Astrophys. J., 715, 671-695 (2010)> * =2010ApJ...715..671W C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! X-ray sources detected in the LkH{alpha} 101 field integer*4 nr__ parameter (nr__=213) ! Number of records character*41 ar__ ! Full-size record integer*4 Seq (nr__) ! [2/231] Running source number (skips indicate * the spurious detections) character*18 ID (nr__) ! Source identification; (HHMMSS.ss+DDMMSS.s; * <LkHa101CXO JHHMMSS.ss+DDMMSS.s> in Simbad) real*4 Off (nr__) ! (arcsec) Offset from axis integer*4 RCts (nr__) ! (ct) Raw counts real*4 NCts (nr__) ! (ct) Net counts c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! 2MASS photometry of X-ray sources integer*4 nr__1 parameter (nr__1=140) ! Number of records character*51 ar__1 ! Full-size record integer*4 Seq_1 (nr__1) ! [2/231] Source index number real*4 X_2M (nr__1) ! (arcsec) X-ray/2MASS positional offset real*4 Off2M (nr__1) ! (arcsec) Offset from axis real*4 Jmag (nr__1) ! (mag) 2MASS J band magnitude real*4 e_Jmag (nr__1) ! (mag) ? Uncertainty in Jmag real*4 Hmag (nr__1) ! (mag) 2MASS H band magnitude real*4 e_Hmag (nr__1) ! (mag) ? Uncertainty in Hmag real*4 Kmag (nr__1) ! (mag) 2MASS K band magnitude real*4 e_Kmag (nr__1) ! (mag) ? Uncertainty in Kmag character*3 Qual (nr__1) ! 2MASS photometry quality codes c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Spitzer photometry of X-ray sources integer*4 nr__2 parameter (nr__2=158) ! Number of records character*89 ar__2 ! Full-size record integer*4 Seq_2 (nr__2) ! [2/231] Source index number (1) real*4 X_Sptz (nr__2) ! (arcsec) X-ray/Spitzer positional offset real*4 Off_1 (nr__2) ! (arcsec) Offset from axis real*4 v_3_6_ (nr__2) ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude real*4 v_4_5_ (nr__2) ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude real*4 v_5_8_ (nr__2) ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude real*4 v_8_0_ (nr__2) ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude real*4 v_24_ (nr__2) ! (mag) ? Spitzer/MIPS 24 micron band magnitude real*4 e__3_6_ (nr__2) ! (mag) ? Uncertainty in [3.6] real*4 e__4_5_ (nr__2) ! (mag) ? Uncertainty in [4.5] real*4 e__5_8_ (nr__2) ! (mag) ? Uncertainty in [5.8] real*4 e__8_0_ (nr__2) ! (mag) ? Uncertainty in [8.0] real*4 e__24_ (nr__2) ! (mag) ? Uncertainty in [24] real*4 AK (nr__2) ! (mag) The K band extinction character*14 Class (nr__2) ! Source classification (1, 2, III, T(transition), * gal (background galaxy), or LkH{alpha} 101) character*1 f_Class (nr__2) ! [b] b: star within 90" of LkH{alpha} 101 (2) *Note (1): In addition to their IAU approved Names: * <LkHa101CXO JHHMMSS.ss+DDMMSS.s> given in Table 1, these sources are * equally identified as: <LkHa101SST JHHMMSS.ss+DDMMSS.s>. *Note (2): Stars with a notation "b" are within 90 arcsec of LkH{alpha} 101. * Their proximity, led to the conclusion that they were cluster * members, despite missing some IRAC data. For these stars 2MASS data * played an enhanced role in determining Class identification. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! LkH{alpha} 101 cluster spectral classification integer*4 nr__3 parameter (nr__3=87) ! Number of records character*57 ar__3 ! Full-size record integer*4 Seq_3 (nr__3) ! [2/231]? Chandra identification number character*18 ID_1 (nr__3) ! LkHaSST identification (HHMMSS.ss+DDMMSS.s; * <LkHa101SST JHHMMSS.ss+DDMMSS.s> in Simbad) character*4 SpT (nr__3) ! MK spectral type real*4 e_SpT (nr__3) ! SpT uncertainty integer*4 Teff (nr__3) ! (K) Effective temperature real*4 logL (nr__3) ! ([-]) Overall luminosity log(L_bol_/L_0_) real*8 EWHa (nr__3) ! (0.1nm) H{alpha} (6563{AA}) equivalent width real*4 EWLi (nr__3) ! (0.1nm) ? LiI (6707{AA}) equivalent width c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Spectral fits to X-ray sources integer*4 nr__4 parameter (nr__4=88) ! Number of records character*64 ar__4 ! Full-size record integer*4 Seq_4 (nr__4) ! [2/231] Source index number real*4 NCts_1 (nr__4) ! (ct) Net counts real*4 NH (nr__4) ! (10+22/cm2) Hydrogen column density real*4 e_NH (nr__4) ! (10+22/cm2) ? The 1{sigma} uncertainty in NH character*1 l_kT (nr__4) ! Limit flag on kT real*4 kT (nr__4) ! (keV) Effective temperature real*4 e_kT (nr__4) ! (keV) ? The 1{sigma} uncertainty in kT real*4 logF_A (nr__4) ! ([mW/m2]) Log of Absorbed flux; erg/s/cm^2^ real*4 logF_U (nr__4) ! ([mW/m2]) Log of Unabsorbed flux; erg/s/cm^2^ real*4 chi2 (nr__4) ! The reduced {chi}^2^ of the fit (1) character*9 Model (nr__4) ! Model used in fit (2) *Note (1): Included for "cstat" models for completeness. One model was * calculated with C statistics, it is unclear how this relates to * goodness of fit. The formal errors are unbiased however. *Note (2): The model used for the fit is coded in two parts. The left hand side * indicates the type of model used. The right hand side indicates the * data set. The overall model is written: LHS_RHS. * LHS as follows: * cstat = C-statistics used for the fit with an initial guess of kT=1keV and * NH=10^21^cm^-2^. * crs = c-stat used for initial guess with those input parameters and then * refit using Raymond-Smith plasma for final fit. * Mekal = again c-stat used for initial guess and then Mekal for final fit, * APEC = cstat used for initial guess and then APEC for final fit. * RHS as follows: * 29 = only data from first observation (Obsid = 5429) used in fit, * 28 = only data from second observation (Obsid = 5428) used in fit, * all = data from both observations used in fit. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Spectral fits to X-ray sources (lower quality) integer*4 nr__5 parameter (nr__5=30) ! Number of records character*64 ar__5 ! Full-size record integer*4 Seq_5 (nr__5) ! [2/231] Source index number real*4 NCts_2 (nr__5) ! (ct) Net counts real*4 NH_1 (nr__5) ! (10+22/cm2) Hydrogen column density real*4 e_NH_1 (nr__5) ! (10+22/cm2) ? The 1{sigma} uncertainty in NH character*1 l_kT_1 (nr__5) ! Limit flag on kT real*4 kT_1 (nr__5) ! (keV) Effective temperature real*4 e_kT_1 (nr__5) ! (keV) ? The 1{sigma} uncertainty in kT real*4 logF_A_1 (nr__5) ! ([mW/m2]) Log of Absorbed flux; erg/s/cm^2^ real*4 logF_U_1 (nr__5) ! ([mW/m2]) Log of Unabsorbed flux; erg/s/cm^2^ real*4 chi2_1 (nr__5) ! The reduced {chi}^2^ of the fit (1) character*9 Model_1 (nr__5) ! Model used in fit (2) *Note (1): Included for "cstat" models for completeness. One model was * calculated with C statistics, it is unclear how this relates to * goodness of fit. The formal errors are unbiased however. *Note (2): The model used for the fit is coded in two parts. The left hand side * indicates the type of model used. The right hand side indicates the * data set. The overall model is written: LHS_RHS. * LHS as follows: * cstat = C-statistics used for the fit with an initial guess of kT=1keV and * NH=10^21^cm^-2^. * crs = c-stat used for initial guess with those input parameters and then * refit using Raymond-Smith plasma for final fit. * Mekal = again c-stat used for initial guess and then Mekal for final fit, * APEC = cstat used for initial guess and then APEC for final fit. * RHS as follows: * 29 = only data from first observation (Obsid = 5429) used in fit, * 28 = only data from second observation (Obsid = 5428) used in fit, * all = data from both observations used in fit. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table10.dat' ! PMS stars detected by Spitzer but not detected in X-rays (2MASS colors) integer*4 nr__6 parameter (nr__6=69) ! Number of records character*55 ar__6 ! Full-size record character*18 ID_2 (nr__6) ! Source identification; (HHMMSS.ss+DDMMSS.s; * <LkHa101SST JHHMMSS.ss+DDMMSS.s> in Simbad) (1) real*4 Jmag_1 (nr__6) ! (mag) ? 2MASS J band magnitude real*4 e_Jmag_1 (nr__6) ! (mag) ? Uncertainty in Jmag real*4 Hmag_1 (nr__6) ! (mag) ? 2MASS H band magnitude real*4 e_Hmag_1 (nr__6) ! (mag) ? Uncertainty in Hmag real*4 Ksmag (nr__6) ! (mag) ? 2MASS Ks band magnitude real*4 e_Ksmag (nr__6) ! (mag) ? Uncertainty in Kmag character*3 Qual_1 (nr__6) ! 2MASS photometry quality codes *Note (1): Source 042955.68+352021.2 has been modified at CDS into * 042955.68+352021.0 as in table 11 and 2MASS catalog. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table11.dat' ! PMS stars not detected in X-rays (Spitzer) integer*4 nr__7 parameter (nr__7=69) ! Number of records character*81 ar__7 ! Full-size record character*18 ID_3 (nr__7) ! Source identification; (HHMMSS.ss+DDMMSS.s; * <LkHa101SST JHHMMSS.ss+DDMMSS.s> in Simbad) character*1 f_ID (nr__7) ! [*] outside Chandra field of view (1) real*4 v_3_6__1 (nr__7) ! (mag) Spitzer/IRAC 3.6 micron band magnitude real*4 e__3_6__1 (nr__7) ! (mag) Uncertainty in [3.6] real*4 v_4_5__1 (nr__7) ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude real*4 e__4_5__1 (nr__7) ! (mag) ? Uncertainty in [4.5] real*4 v_5_6_ (nr__7) ! (mag) ? Spitzer/IRAC 5.6 micron band magnitude real*4 e__5_6_ (nr__7) ! (mag) ? Uncertainty in [5.6] real*4 v_8_0__1 (nr__7) ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude real*4 e__8_0__1 (nr__7) ! (mag) ? Uncertainty in [8.0] real*4 v_24__1 (nr__7) ! (mag) ? Spitzer/MIPS 24 micron band magnitude real*4 e__24__1 (nr__7) ! (mag) ? Uncertainty in [24] real*4 AK_1 (nr__7) ! (mag) The K band extinction character*1 Cl (nr__7) ! [012T] Classification (Class 0, I, II or * T=possible transition objects) *Note (1): * * = These Spitzer identified PMS stars are off the Chandra field of view. C============================================================================= C Loading file 'table1.dat' ! X-ray sources detected in the LkH{alpha} 101 field C Format for file interpretation 1 format(I3,1X,A18,1X,F6.1,1X,I4,1X,F6.1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,213 read(1,'(A41)')ar__ read(ar__,1)Seq(i__),ID(i__),Off(i__),RCts(i__),NCts(i__) c ..............Just test output........... write(6,1)Seq(i__),ID(i__),Off(i__),RCts(i__),NCts(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! 2MASS photometry of X-ray sources C Format for file interpretation 2 format( + I3,1X,F4.2,1X,F5.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X, + F4.2,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,140 read(1,'(A51)')ar__1 read(ar__1,2) + Seq_1(i__),X_2M(i__),Off2M(i__),Jmag(i__),e_Jmag(i__), + Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),Qual(i__) if(ar__1(22:25) .EQ. '') e_Jmag(i__) = rNULL__ if(ar__1(33:36) .EQ. '') e_Hmag(i__) = rNULL__ if(ar__1(44:47) .EQ. '') e_Kmag(i__) = rNULL__ c ..............Just test output........... write(6,2) + Seq_1(i__),X_2M(i__),Off2M(i__),Jmag(i__),e_Jmag(i__), + Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),Qual(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Spitzer photometry of X-ray sources C Format for file interpretation 3 format( + I3,1X,F4.2,1X,F5.1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X, + F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A14,A1) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,158 read(1,'(A89)')ar__2 read(ar__2,3) + Seq_2(i__),X_Sptz(i__),Off_1(i__),v_3_6_(i__),v_4_5_(i__), + v_5_8_(i__),v_8_0_(i__),v_24_(i__),e__3_6_(i__),e__4_5_(i__), + e__5_8_(i__),e__8_0_(i__),e__24_(i__),AK(i__),Class(i__), + f_Class(i__) if(ar__2(16:20) .EQ. '') v_3_6_(i__) = rNULL__ if(ar__2(22:26) .EQ. '') v_4_5_(i__) = rNULL__ if(ar__2(28:32) .EQ. '') v_5_8_(i__) = rNULL__ if(ar__2(34:38) .EQ. '') v_8_0_(i__) = rNULL__ if(ar__2(40:43) .EQ. '') v_24_(i__) = rNULL__ if(ar__2(45:48) .EQ. '') e__3_6_(i__) = rNULL__ if(ar__2(50:53) .EQ. '') e__4_5_(i__) = rNULL__ if(ar__2(55:58) .EQ. '') e__5_8_(i__) = rNULL__ if(ar__2(60:63) .EQ. '') e__8_0_(i__) = rNULL__ if(ar__2(65:68) .EQ. '') e__24_(i__) = rNULL__ c ..............Just test output........... write(6,3) + Seq_2(i__),X_Sptz(i__),Off_1(i__),v_3_6_(i__),v_4_5_(i__), + v_5_8_(i__),v_8_0_(i__),v_24_(i__),e__3_6_(i__),e__4_5_(i__), + e__5_8_(i__),e__8_0_(i__),e__24_(i__),AK(i__),Class(i__), + f_Class(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! LkH{alpha} 101 cluster spectral classification C Format for file interpretation 4 format(I3,1X,A18,1X,A4,1X,F3.1,1X,I5,1X,F6.3,1X,F8.2,1X,F3.1) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,87 read(1,'(A57)')ar__3 read(ar__3,4) + Seq_3(i__),ID_1(i__),SpT(i__),e_SpT(i__),Teff(i__),logL(i__), + EWHa(i__),EWLi(i__) if(ar__3(1:3) .EQ. '') Seq_3(i__) = iNULL__ if(ar__3(55:57) .EQ. '') EWLi(i__) = rNULL__ c ..............Just test output........... write(6,4) + Seq_3(i__),ID_1(i__),SpT(i__),e_SpT(i__),Teff(i__),logL(i__), + EWHa(i__),EWLi(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! Spectral fits to X-ray sources C Format for file interpretation 5 format( + I3,1X,F6.1,1X,F5.2,1X,F4.2,1X,A1,1X,F6.2,1X,F4.2,1X,F6.2,1X, + F6.2,1X,F4.2,1X,A9) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,88 read(1,'(A64)')ar__4 read(ar__4,5) + Seq_4(i__),NCts_1(i__),NH(i__),e_NH(i__),l_kT(i__),kT(i__), + e_kT(i__),logF_A(i__),logF_U(i__),chi2(i__),Model(i__) if(ar__4(18:21) .EQ. '') e_NH(i__) = rNULL__ if(ar__4(32:35) .EQ. '') e_kT(i__) = rNULL__ c ..............Just test output........... write(6,5) + Seq_4(i__),NCts_1(i__),NH(i__),e_NH(i__),l_kT(i__),kT(i__), + e_kT(i__),logF_A(i__),logF_U(i__),chi2(i__),Model(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! Spectral fits to X-ray sources (lower quality) C Format for file interpretation 6 format( + I3,1X,F6.1,1X,F5.2,1X,F4.2,1X,A1,1X,F6.2,1X,F4.2,1X,F6.2,1X, + F6.2,1X,F4.2,1X,A9) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,30 read(1,'(A64)')ar__5 read(ar__5,6) + Seq_5(i__),NCts_2(i__),NH_1(i__),e_NH_1(i__),l_kT_1(i__), + kT_1(i__),e_kT_1(i__),logF_A_1(i__),logF_U_1(i__),chi2_1(i__), + Model_1(i__) if(ar__5(18:21) .EQ. '') e_NH_1(i__) = rNULL__ if(ar__5(32:35) .EQ. '') e_kT_1(i__) = rNULL__ c ..............Just test output........... write(6,6) + Seq_5(i__),NCts_2(i__),NH_1(i__),e_NH_1(i__),l_kT_1(i__), + kT_1(i__),e_kT_1(i__),logF_A_1(i__),logF_U_1(i__),chi2_1(i__), + Model_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table10.dat' ! PMS stars detected by Spitzer but not detected in * X-rays (2MASS colors) C Format for file interpretation 7 format( + A18,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table10.dat') write(6,*) '....Loading file: table10.dat' do i__=1,69 read(1,'(A55)')ar__6 read(ar__6,7) + ID_2(i__),Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__), + Ksmag(i__),e_Ksmag(i__),Qual_1(i__) if(ar__6(20:24) .EQ. '') Jmag_1(i__) = rNULL__ if(ar__6(26:29) .EQ. '') e_Jmag_1(i__) = rNULL__ if(ar__6(31:35) .EQ. '') Hmag_1(i__) = rNULL__ if(ar__6(37:40) .EQ. '') e_Hmag_1(i__) = rNULL__ if(ar__6(42:46) .EQ. '') Ksmag(i__) = rNULL__ if(ar__6(48:51) .EQ. '') e_Ksmag(i__) = rNULL__ c ..............Just test output........... write(6,7) + ID_2(i__),Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__), + Ksmag(i__),e_Ksmag(i__),Qual_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table11.dat' ! PMS stars not detected in X-rays (Spitzer) C Format for file interpretation 8 format( + A18,1X,A1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X, + F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A1) C Effective file loading open(unit=1,status='old',file= +'table11.dat') write(6,*) '....Loading file: table11.dat' do i__=1,69 read(1,'(A81)')ar__7 read(ar__7,8) + ID_3(i__),f_ID(i__),v_3_6__1(i__),e__3_6__1(i__), + v_4_5__1(i__),e__4_5__1(i__),v_5_6_(i__),e__5_6_(i__), + v_8_0__1(i__),e__8_0__1(i__),v_24__1(i__),e__24__1(i__), + AK_1(i__),Cl(i__) if(ar__7(33:37) .EQ. '') v_4_5__1(i__) = rNULL__ if(ar__7(39:42) .EQ. '') e__4_5__1(i__) = rNULL__ if(ar__7(44:48) .EQ. '') v_5_6_(i__) = rNULL__ if(ar__7(50:53) .EQ. '') e__5_6_(i__) = rNULL__ if(ar__7(55:59) .EQ. '') v_8_0__1(i__) = rNULL__ if(ar__7(61:64) .EQ. '') e__8_0__1(i__) = rNULL__ if(ar__7(66:69) .EQ. '') v_24__1(i__) = rNULL__ if(ar__7(71:74) .EQ. '') e__24__1(i__) = rNULL__ c ..............Just test output........... write(6,8) + ID_3(i__),f_ID(i__),v_3_6__1(i__),e__3_6__1(i__), + v_4_5__1(i__),e__4_5__1(i__),v_5_6_(i__),e__5_6_(i__), + v_8_0__1(i__),e__8_0__1(i__),v_24__1(i__),e__24__1(i__), + AK_1(i__),Cl(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end