Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/193/25 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-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/ApJS/193/25 The NAN complex. II. MIPS observations (Rebull+, 2011) *================================================================================ *The North American and Pelican nebulae. *II. MIPS observations and analysis. * Rebull L.M., Guieu S., Stauffer J.R., Hillenbrand L.A., Noriega-Crespo A., * Stapelfeldt K.R., Carey S.J., Carpenter J.M., Cole D.M., Padgett D.L., * Strom S.E., Wolff S.C. * <Astrophys. J. Suppl. Ser., 193, 25 (2011)> * =2011ApJS..193...25R C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Complete list of YSOs I. Spitzer measurements integer*4 nr__ parameter (nr__=1286) ! Number of records character*187 ar__ ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 v_RGS2011_ ! Source name (HHMMSS.ss+DDMMSS.s on J2000) real*8 RAdeg ! (deg) Right Ascension in decimal degrees (J2000) real*8 DEdeg ! (deg) Declination in decimal degrees (J2000) character*16 v2MASS ! 2MASS name (HHMMSSss+DDMMSSs) character*23 AName ! Alternate name character*12 SpType ! MK spectral type character*22 r_SpType ! Reference for SpType (1) character*1 l__3_6_ ! Limit flag on 3.6] real*4 v_3_6_ ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude real*4 e__3_6_ ! (mag) ? Uncertainty in 3.6mag character*1 l__4_5_ ! Limit flag on 4.5] real*4 v_4_5_ ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude real*4 e__4_5_ ! (mag) ? Uncertainty in 4.5mag character*1 l__5_8_ ! Limit flag on 5.8] real*4 v_5_8_ ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude real*4 e__5_8_ ! (mag) ? Uncertainty in 5.8mag character*1 l__8_0_ ! Limit flag on 8.0] real*4 v_8_0_ ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude real*4 e__8_0_ ! (mag) ? Uncertainty in 8.0mag character*1 l__24_ ! Limit flag on 24] real*4 v_24_ ! (mag) ? Spitzer/MIPS 24 micron band magnitude real*4 e__24_ ! (mag) ? Uncertainty in 24mag real*4 v_70_ ! (mag) ? Spitzer/MIPS 70 micron band magnitude real*4 e__70_ ! (mag) ? Uncertainty in 70mag *Note (1): References as follows: * BallyScoville = 1980ApJ...239..121B * CohenKuhi79 = 1979ApJS...41..743C * Comeron05 = 2005A&A...430..541C; <[CP2005] NN> in Simbad * Corbally09 = 2009BaltA..18..111C; <[CSL2009] CSL N> in Simbad * Laugalys02 = 2002BaltA..11..205L * Laugalys06a = 2006BaltA..15..483L (Cat. J/BaltA/15/483; * <[LSV2006] Area II NNN> in Simbad) * Laugalys06b = 2006BaltA..15..327L * MerrillBurwell49 = 1949ApJ...110..387M * Ogura2002 = 2002AJ....123.2597O (Cat. J/AJ/123/2597; * <[OSP2002] BRC 31 NN> in Simbad) * Shevchenko93 = 1993Ap&SS.202..121S * SIMBAD = http://simbad.u-strasbg.fr/simbad/ * Straizys08a = 2008BaltA..17....1S (Cat. J/BaltA/17/1) * Straizys08b = 2008BaltA..17..143S; <[SL2008c] NN> in Simbad * Straizys09 = 2009BaltA..18..141S (Cat. J/BaltA/18/141) * Straizys89 = 1989VilOB..83....3S * Straizys93 = 1993BaltA...2..171S (Cat. II/196) * Straizys99 = 1999BaltA...8..355S * Terranegra94 = 1994A&AS..104..557T (Cat. J/A+AS/104/557) * Vill04 = 2004A&A...419..149V (Cat. J/A+A/419/149; * <Cl* NGC 6996 VBC NNNN> in Simbad) * Welin73 = 1973A&AS....9..183W * Zdan90 = 1990Ap&SS.173..309Z c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Complete list of YSOs II. YSO criteria integer*4 nr__1 parameter (nr__1=1286) ! Number of records character*82 ar__1 ! Full-size record character*18 v_RGS2011__1 ! Source name (HHMMSS.ss+DDMMSS.s) character*8 Clu ! Cluster Name (1) character*1 Q ! [ABC] YSO quality (A for good confidence) character*4 cl ! YSO class (2) character*3 cr1 ! [yes/no] Selected from [3.6] vs [3.6]-[24] CMD (3) character*3 cr2 ! [yes/no] Not too faint in [3.6] vs [3.6]-[24] * CMD (3) character*3 cr3 ! [yes/no] Selected from Ks vs Ks-[24] CMD (3) character*3 cr4 ! [yes/no] Detected at 70 microns (3) character*3 cr5 ! [yes/no] Found in Guieu et al. 2009 * Cat. J/ApJ/697/787 ([GRS2009] in Simbad) character*3 cr6 ! [yes/no] Close to the middle of the nebula (3) character*3 cr7 ! [yes/no] In a cluster (3) character*3 cr8 ! [yes/no] Previously identified YSO (3) character*3 cr9 ! [yes/no] Selected from Ic vs. V-Ic CMD (3) character*3 cr10 ! [yes/no] Selected from IPHAS r' vs. r'-i' CMD (3) character*3 cr11 ! [yes/no] A bright IPHAS H{alpha} source (3) character*3 cr12 ! [yes/no] Selected from IPHAS r'-H{alpha} vs * r'-i' diagram (3) *Note (1): Clusters as follows: * Pelican = IC 5070. See section 6.3 for further details. * PelHat = Pelican's Hat. North of the Pelican's head, in roughly the same * location as would be found a hat on its head, is a previously * unidentified dark lane at 24um, similar to but much smaller than * the Gulf of Mexico. See section 6.4 for further details. * Gulf = The "Gulf of Mexico" region is the area of highest extinction in * the NAN (The North American (NGC 7000) and Pelican (IC 5070) * Nebulae) complex. See section 6.2 for further details. *Note (2): In the spirit of Wilking et al. (2001ApJ...551..357W), we define the * near- to mid-IR slope of the SED, * {alpha}=dlog{lambda}F_{lambda}_/dlog{lambda}, where {alpha}>0.3 for a * Class I, 0.3 to -0.3 for a flat-spectrum source, -0.3 to -1.6 for a * Class II, and <-1.6 for a Class III. *Note (3): For each individual criterion, a "yes" means the object met that * given criterion for YSO selection, "no" means it did not, and blank * means that it was unable to be assessed for that criterion. * For example, column S3624 indicates whether each object was selected * as a YSO candidate in the [3.6] vs. [3.6]-[24] diagram ("yes"), or * it was not ("no"), or, it was missing a 3.6 or 24 micron * measurement, and thus could not be assessed in the [3.6] vs. * [3.6]-[24] diagram (blank). * See section 4.2 for further details on each criterion. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! IRAC-selected YSOs not recovered by our MIPS-based search integer*4 nr__2 parameter (nr__2=796) ! Number of records character*187 ar__2 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 v_RGS2011__2 ! Source name (HHMMSS.ss+DDMMSS.s on J2000) real*8 RAdeg_1 ! (deg) Right Ascension in decimal degrees (J2000) real*8 DEdeg_1 ! (deg) Declination in decimal degrees (J2000) character*16 v2MASS_1 ! 2MASS name (HHMMSSss+DDMMSSs) character*23 AName_1 ! Alternate name character*12 SpType_1 ! MK spectral type character*22 r_SpType_1 ! Reference for SpType (1) character*1 l__3_6__1 ! Limit flag on 3.6] real*4 v_3_6__1 ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude real*4 e__3_6__1 ! (mag) ? Uncertainty in 3.6mag character*1 l__4_5__1 ! Limit flag on 4.5] real*4 v_4_5__1 ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude real*4 e__4_5__1 ! (mag) ? Uncertainty in 4.5mag character*1 l__5_8__1 ! Limit flag on 5.8] real*4 v_5_8__1 ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude real*4 e__5_8__1 ! (mag) ? Uncertainty in 5.8mag character*1 l__8_0__1 ! Limit flag on 8.0] real*4 v_8_0__1 ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude real*4 e__8_0__1 ! (mag) ? Uncertainty in 8.0mag character*1 l__24__1 ! Limit flag on 24] real*4 v_24__1 ! (mag) ? Spitzer/MIPS 24 micron band magnitude real*4 e__24__1 ! (mag) ? Uncertainty in 24mag real*4 v_70__1 ! (mag) ? Spitzer/MIPS 70 micron band magnitude real*4 e__70__1 ! (mag) ? Uncertainty in 70mag *Note (1): References as follows: * BallyScoville = 1980ApJ...239..121B * CohenKuhi79 = 1979ApJS...41..743C * Comeron05 = 2005A&A...430..541C; <[CP2005] NN> in Simbad * Corbally09 = 2009BaltA..18..111C; <[CSL2009] CSL N> in Simbad * Laugalys02 = 2002BaltA..11..205L * Laugalys06a = 2006BaltA..15..483L (Cat. J/BaltA/15/483; * <[LSV2006] Area II NNN> in Simbad) * Laugalys06b = 2006BaltA..15..327L * MerrillBurwell49 = 1949ApJ...110..387M * Ogura2002 = 2002AJ....123.2597O (Cat. J/AJ/123/2597; * <[OSP2002] BRC 31 NN> in Simbad) * Shevchenko93 = 1993Ap&SS.202..121S * SIMBAD = http://simbad.u-strasbg.fr/simbad/ * Straizys08a = 2008BaltA..17....1S (Cat. J/BaltA/17/1) * Straizys08b = 2008BaltA..17..143S; <[SL2008c] NN> in Simbad * Straizys09 = 2009BaltA..18..141S (Cat. J/BaltA/18/141) * Straizys89 = 1989VilOB..83....3S * Straizys93 = 1993BaltA...2..171S (Cat. II/196) * Straizys99 = 1999BaltA...8..355S * Terranegra94 = 1994A&AS..104..557T (Cat. J/A+AS/104/557) * Vill04 = 2004A&A...419..149V (Cat. J/A+A/419/149; * <Cl* NGC 6996 VBC NNNN> in Simbad) * Welin73 = 1973A&AS....9..183W * Zdan90 = 1990Ap&SS.173..309Z c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table9.dat' ! Previously suggested YSOs not recovered by our MIPS or IRAC searches integer*4 nr__3 parameter (nr__3=114) ! Number of records character*187 ar__3 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 v_RGS2011__3 ! Source name (HHMMSS.ss+DDMMSS.s on J2000) real*8 RAdeg_2 ! (deg) Right Ascension in decimal degrees (J2000) real*8 DEdeg_2 ! (deg) Declination in decimal degrees (J2000) character*16 v2MASS_2 ! 2MASS name (HHMMSSss+DDMMSSs) character*23 AName_2 ! Alternate name character*12 SpType_2 ! MK spectral type character*22 r_SpType_2 ! Reference for SpType (1) character*1 l__3_6__2 ! Limit flag on 3.6] real*4 v_3_6__2 ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude real*4 e__3_6__2 ! (mag) ? Uncertainty in 3.6mag character*1 l__4_5__2 ! Limit flag on 4.5] real*4 v_4_5__2 ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude real*4 e__4_5__2 ! (mag) ? Uncertainty in 4.5mag character*1 l__5_8__2 ! Limit flag on 5.8] real*4 v_5_8__2 ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude real*4 e__5_8__2 ! (mag) ? Uncertainty in 5.8mag character*1 l__8_0__2 ! Limit flag on 8.0] real*4 v_8_0__2 ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude real*4 e__8_0__2 ! (mag) ? Uncertainty in 8.0mag character*1 l__24__2 ! Limit flag on 24] real*4 v_24__2 ! (mag) ? Spitzer/MIPS 24 micron band magnitude real*4 e__24__2 ! (mag) ? Uncertainty in 24mag real*4 v_70__2 ! (mag) ? Spitzer/MIPS 70 micron band magnitude real*4 e__70__2 ! (mag) ? Uncertainty in 70mag *Note (1): References as follows: * BallyScoville = 1980ApJ...239..121B * CohenKuhi79 = 1979ApJS...41..743C * Comeron05 = 2005A&A...430..541C; <[CP2005] NN> in Simbad * Corbally09 = 2009BaltA..18..111C; <[CSL2009] CSL N> in Simbad * Laugalys02 = 2002BaltA..11..205L * Laugalys06a = 2006BaltA..15..483L (Cat. J/BaltA/15/483; * <[LSV2006] Area II NNN> in Simbad) * Laugalys06b = 2006BaltA..15..327L * MerrillBurwell49 = 1949ApJ...110..387M * Ogura2002 = 2002AJ....123.2597O (Cat. J/AJ/123/2597; * <[OSP2002] BRC 31 NN> in Simbad) * Shevchenko93 = 1993Ap&SS.202..121S * SIMBAD = http://simbad.u-strasbg.fr/simbad/ * Straizys08a = 2008BaltA..17....1S (Cat. J/BaltA/17/1) * Straizys08b = 2008BaltA..17..143S; <[SL2008c] NN> in Simbad * Straizys09 = 2009BaltA..18..141S (Cat. J/BaltA/18/141) * Straizys89 = 1989VilOB..83....3S * Straizys93 = 1993BaltA...2..171S (Cat. II/196) * Straizys99 = 1999BaltA...8..355S * Terranegra94 = 1994A&AS..104..557T (Cat. J/A+AS/104/557) * Vill04 = 2004A&A...419..149V (Cat. J/A+A/419/149; * <Cl* NGC 6996 VBC NNNN> in Simbad) * Welin73 = 1973A&AS....9..183W * Zdan90 = 1990Ap&SS.173..309Z C============================================================================= C Loading file 'table3.dat' ! Complete list of YSOs I. Spitzer measurements C Format for file interpretation 1 format( + A18,1X,F10.6,1X,F9.6,1X,A16,1X,A23,1X,A12,1X,A22,1X,A1,F5.2, + 1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X, + F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,1286 read(1,'(A187)')ar__ read(ar__,1) + v_RGS2011_,RAdeg,DEdeg,v2MASS,AName,SpType,r_SpType,l__3_6_, + v_3_6_,e__3_6_,l__4_5_,v_4_5_,e__4_5_,l__5_8_,v_5_8_,e__5_8_, + l__8_0_,v_8_0_,e__8_0_,l__24_,v_24_,e__24_,v_70_,e__70_ if(ar__(119:123) .EQ. '') v_3_6_ = rNULL__ if(ar__(125:128) .EQ. '') e__3_6_ = rNULL__ if(ar__(131:135) .EQ. '') v_4_5_ = rNULL__ if(ar__(137:140) .EQ. '') e__4_5_ = rNULL__ if(ar__(143:147) .EQ. '') v_5_8_ = rNULL__ if(ar__(149:152) .EQ. '') e__5_8_ = rNULL__ if(ar__(155:159) .EQ. '') v_8_0_ = rNULL__ if(ar__(161:164) .EQ. '') e__8_0_ = rNULL__ if(ar__(167:171) .EQ. '') v_24_ = rNULL__ if(ar__(173:176) .EQ. '') e__24_ = rNULL__ if(ar__(178:182) .EQ. '') v_70_ = rNULL__ if(ar__(184:187) .EQ. '') e__70_ = rNULL__ c ..............Just test output........... write(6,1) + v_RGS2011_,RAdeg,DEdeg,v2MASS,AName,SpType,r_SpType,l__3_6_, + v_3_6_,e__3_6_,l__4_5_,v_4_5_,e__4_5_,l__5_8_,v_5_8_,e__5_8_, + l__8_0_,v_8_0_,e__8_0_,l__24_,v_24_,e__24_,v_70_,e__70_ c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Complete list of YSOs II. YSO criteria C Format for file interpretation 2 format( + A18,1X,A8,1X,A1,1X,A4,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X, + A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,1286 read(1,'(A82)')ar__1 read(ar__1,2) + v_RGS2011__1,Clu,Q,cl,cr1,cr2,cr3,cr4,cr5,cr6,cr7,cr8,cr9, + cr10,cr11,cr12 c ..............Just test output........... write(6,2) + v_RGS2011__1,Clu,Q,cl,cr1,cr2,cr3,cr4,cr5,cr6,cr7,cr8,cr9, + cr10,cr11,cr12 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! IRAC-selected YSOs not recovered by our * MIPS-based search C Format for file interpretation 3 format( + A18,1X,F10.6,1X,F9.6,1X,A16,1X,A23,1X,A12,1X,A22,1X,A1,F5.2, + 1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X, + F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,796 read(1,'(A187)')ar__2 read(ar__2,3) + v_RGS2011__2,RAdeg_1,DEdeg_1,v2MASS_1,AName_1,SpType_1, + r_SpType_1,l__3_6__1,v_3_6__1,e__3_6__1,l__4_5__1,v_4_5__1, + e__4_5__1,l__5_8__1,v_5_8__1,e__5_8__1,l__8_0__1,v_8_0__1, + e__8_0__1,l__24__1,v_24__1,e__24__1,v_70__1,e__70__1 if(ar__2(119:123) .EQ. '') v_3_6__1 = rNULL__ if(ar__2(125:128) .EQ. '') e__3_6__1 = rNULL__ if(ar__2(131:135) .EQ. '') v_4_5__1 = rNULL__ if(ar__2(137:140) .EQ. '') e__4_5__1 = rNULL__ if(ar__2(143:147) .EQ. '') v_5_8__1 = rNULL__ if(ar__2(149:152) .EQ. '') e__5_8__1 = rNULL__ if(ar__2(155:159) .EQ. '') v_8_0__1 = rNULL__ if(ar__2(161:164) .EQ. '') e__8_0__1 = rNULL__ if(ar__2(167:171) .EQ. '') v_24__1 = rNULL__ if(ar__2(173:176) .EQ. '') e__24__1 = rNULL__ if(ar__2(178:182) .EQ. '') v_70__1 = rNULL__ if(ar__2(184:187) .EQ. '') e__70__1 = rNULL__ c ..............Just test output........... write(6,3) + v_RGS2011__2,RAdeg_1,DEdeg_1,v2MASS_1,AName_1,SpType_1, + r_SpType_1,l__3_6__1,v_3_6__1,e__3_6__1,l__4_5__1,v_4_5__1, + e__4_5__1,l__5_8__1,v_5_8__1,e__5_8__1,l__8_0__1,v_8_0__1, + e__8_0__1,l__24__1,v_24__1,e__24__1,v_70__1,e__70__1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table9.dat' ! Previously suggested YSOs not recovered by our * MIPS or IRAC searches C Format for file interpretation 4 format( + A18,1X,F10.6,1X,F9.6,1X,A16,1X,A23,1X,A12,1X,A22,1X,A1,F5.2, + 1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X, + F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table9.dat') write(6,*) '....Loading file: table9.dat' do i__=1,114 read(1,'(A187)')ar__3 read(ar__3,4) + v_RGS2011__3,RAdeg_2,DEdeg_2,v2MASS_2,AName_2,SpType_2, + r_SpType_2,l__3_6__2,v_3_6__2,e__3_6__2,l__4_5__2,v_4_5__2, + e__4_5__2,l__5_8__2,v_5_8__2,e__5_8__2,l__8_0__2,v_8_0__2, + e__8_0__2,l__24__2,v_24__2,e__24__2,v_70__2,e__70__2 if(ar__3(119:123) .EQ. '') v_3_6__2 = rNULL__ if(ar__3(125:128) .EQ. '') e__3_6__2 = rNULL__ if(ar__3(131:135) .EQ. '') v_4_5__2 = rNULL__ if(ar__3(137:140) .EQ. '') e__4_5__2 = rNULL__ if(ar__3(143:147) .EQ. '') v_5_8__2 = rNULL__ if(ar__3(149:152) .EQ. '') e__5_8__2 = rNULL__ if(ar__3(155:159) .EQ. '') v_8_0__2 = rNULL__ if(ar__3(161:164) .EQ. '') e__8_0__2 = rNULL__ if(ar__3(167:171) .EQ. '') v_24__2 = rNULL__ if(ar__3(173:176) .EQ. '') e__24__2 = rNULL__ if(ar__3(178:182) .EQ. '') v_70__2 = rNULL__ if(ar__3(184:187) .EQ. '') e__70__2 = rNULL__ c ..............Just test output........... write(6,4) + v_RGS2011__3,RAdeg_2,DEdeg_2,v2MASS_2,AName_2,SpType_2, + r_SpType_2,l__3_6__2,v_3_6__2,e__3_6__2,l__4_5__2,v_4_5__2, + e__4_5__2,l__5_8__2,v_5_8__2,e__5_8__2,l__8_0__2,v_8_0__2, + e__8_0__2,l__24__2,v_24__2,e__24__2,v_70__2,e__70__2 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end