Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/593/1093 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-18 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/593/1093 IC 348 membership (Luhman+, 2003) *================================================================================ *A census of the young cluster IC 348. * Luhman K.L., Stauffer J.R., Muench A.A., Rieke G.H., Lada E.A., Bouvier J., * Lada C.J. * <Astrophys. J. 593, 1093 (2003)> * =2003ApJ...593.1093L C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Data for Members of IC 348 integer*4 nr__ parameter (nr__=288) ! Number of records character*146 ar__ ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Seq (nr__) ! Identification number, * Cl* IC 348 LRL NNNNNA in Simbad character*1 m_Seq (nr__) ! Multiplicity on Seq character*1 f_Seq (nr__) ! [hijkl] Flag on Seq (1) integer*4 RAh (nr__) ! (h) Hour of Right Ascension (J2000) (2) integer*4 RAm (nr__) ! (min) Minute of Right Ascension (J2000) (2) real*4 RAs (nr__) ! (s) Second of Right Ascension (J2000) (2) character*1 DE_ (nr__) ! Declination sign (J2000) integer*4 DEd (nr__) ! (deg) Degree of Declination (J2000) (2) integer*4 DEm (nr__) ! (arcmin) Arcminute of Declination (J2000) (2) real*4 DEs (nr__) ! (arcsec) Arcsecond of Declination (J2000) (2) character*32 SpType (nr__) ! Spectral Type/H{alpha} equivalent width (3) character*5 r_SpType (nr__) ! References on SpType (4) character*5 Adopt (nr__) ! Adopted spectral type character*10 Mem (nr__) ! Membership in IC 348 code (5) integer*4 Teff (nr__) ! (K) ? Effective temperature (6) real*4 AJ (nr__) ! (mag) ? Extinction in the J band real*8 LBol (nr__) ! (solLum) ? Bolometric luminosity real*4 R_I (nr__) ! (mag) ? The (R-I) color (7) real*4 Imag (nr__) ! (mag) ? The I band magnitude (8) real*4 I_Z (nr__) ! (mag) ? The (I-Z) color from this work real*4 J_H (nr__) ! (mag) ? The (J-H) color real*4 H_Ks (nr__) ! (mag) ? The (H-Ks) color real*4 Kmag (nr__) ! (mag) ? The K_s_ band magnitude character*1 r_J_H (nr__) ! JHKs reference (9) character*3 IMF (nr__) ! [yes no] Included in the IMF sample? *Note (1): the flag has the following meaning: * h = Sources 1, 6, 9, 52, 86, 92, 221, and 226 are resolved as binaries with * separations of 0.61, 0.559, 0.409, 1.028, 0.249, 0.26, 0.13, and 0.70 * arcseconds and magnitude differences of {Delta}H = 0.14, 1.49, 1.84, * 1.48, 0.67, 0.67, 1.35, and 0.75, respectively (Duchene et al., * 1999A&A...343..831D). * i = R-I and JHK_s_ are for A+B. * j = JHK_s_ are for A+B. * k = Spectral type and JHK_s_ are for A+B. LBol is for A. * l = Spectral type and JHK_s_ are for A+B. *Note (2): From the 2MASS Point Source Catalog for 1, 2, 3, 4, 6, 7, 8, * 10, 16, 19, 63, 88, 90, 122, 242, 286, and 367. The coordinates for 91 * are derived in the manner described in Section 4.1. All other * coordinates are from the I-band images in this work. *Note (3): Measurement uncertainties for the optical spectral types from * Luhman (1999, Cat. <J/ApJ/525/466>) and this work are +/-0.5 and 0.25 * subclass for K and M types, respectively, unless noted otherwise. For * H{alpha} equivalent widths from Luhman et al. (1998ApJ...508..347L), * Luhman (1999, Cat. <J/ApJ/525/466>), and this work, uncertainties are * 0.5-1{AA} unless noted otherwise. *Note (4): References are numbered as follows: * 1 = Harris et al. (1954ApJ...119..622H); * 2 = Luhman et al. (1998ApJ...508..347L) * 4 = Strom et al. (1974PASP...86..798S); * 5 = Herbig (1998, Cat. <J/ApJ/497/736>); * 6 = Measurements described by Luhman (1999, Cat. <J/ApJ/525/466>) and * presented here for the first time; * 7 = this work; * 8 = Luhman (1999, Cat. <J/ApJ/525/466>). *Note (5): Membership in IC 348 code use the following notation: * A_V_=A_V_>~1 and a position above the main sequence for the distance * of IC 348; * e = emission lines; * NaK = spectral features such as Na I and K I; * ex = Ks excess emission; * pm = proper motion measurements (Fredrick, 1956AJ.....61..437F). *Note (6): From Schmidt-Kaler (1982, in Landolt-Bornstein, Group VI, Vol. 2, * ed. K.-H. Hellwege (Berlin: Springer), 454) (<=M0) and Table 8 (>M0). *Note (7): From Luhman (1999, Cat. <J/ApJ/525/466>) for stars with R-I>=1.5 * and for 435 and from Herbig (1998, Cat. <J/ApJ/497/736>) for the * remaining objects. *Note (8): From Herbig (1998, Cat. <J/ApJ/497/736>) for 16, 33, 63, 71, * 88, 90, and 91, from Luhman (1999, Cat. <J/ApJ/525/466>) for 367, and * from this work for the remaining objects. *Note (9): JHKs reference: * 3 = 2MASS Point Source Catalog; Cat. <II/246> * 9 = Muench et al. (2003AJ....125.2029M). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Foreground Stars integer*4 nr__1 parameter (nr__1=20) ! Number of records character*106 ar__1 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000 real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Seq_1 (nr__1) ! Sequential number integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000) (1) integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000) real*4 RAs_1 (nr__1) ! (s) Right ascension (J2000) character*1 DE__1 (nr__1) ! Declination sign (J2000) integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000) integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000) real*4 DEs_1 (nr__1) ! (arcsec) Declination (J2000) character*24 SpType_1 (nr__1) ! Spectral type character*5 r_SpType_1 (nr__1) ! Reference for spectral type (6) character*9 FE (nr__1) ! Foregroung or Background evidence (2) real*4 R_I_1 (nr__1) ! (mag) ? R-I colour index (3) real*4 Imag_1 (nr__1) ! (mag) ? I magnitude (4) real*4 I_Z_1 (nr__1) ! (mag) ? I-Z colour index (4) real*4 J_H_1 (nr__1) ! (mag) ? I-Z colour index (5) real*4 H_Ks_1 (nr__1) ! (mag) ? H-Ks colour index (5) real*4 Ksmag (nr__1) ! (mag) ? Ks magnitude (5) *Note (1): The positions are from the 2MASS Point Source Catalog for * 18, 28, 206 and 370. All other coordinates are from the I-band images * in this work. *Note (2): Foreground nature is indicated by proper-motion * measurements ("pm"; Blaauw, 1952BAN....11..405B; Scholz et al., 1999, * Cat. <J/A+AS/137/305>) or a combination of no Li absorption ("Li"), * A_V_~0 ("A_V_"), and spectral features such as strong Na I and K I * ("NaK"; Luhman et al., 1998ApJ...508..347L; this work). *Note (3): From Herbig, 1998, Cat. <J/ApJ/497/736> for 77, 106 and 106 and * from Luhman 1999, Cat. <J/ApJ/525/466> for the remaining objects. *Note (4): From this work. *Note (5): From 2MASS for 14, 18, 27, 28, 34, 43, 54, 57, 73, 77, 80, 131, * 189, 320, 400, 1788, 131, 1932, 1941, 10226, 10252, and 10338 and from * Muench et al. (2003AJ....125.2029M) for the remaining objects. *Note (6): References are numbered as follows: * (1): Luhman et al., 1998ApJ...508..347L * (2): this work. * (3): measurements described by Luhman 1999, Cat. <J/ApJ/525/466> and * presented here for the first time; * (4): Fredrick, 1956AJ.....61..437F * (5): Herbig, 1998, Cat. <J/ApJ/497/736> c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Background Stars integer*4 nr__2 parameter (nr__2=103) ! Number of records character*107 ar__2 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension J2000 real*8 DEdeg_2 (nr__2) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Seq_2 (nr__2) ! Sequential number integer*4 RAh_2 (nr__2) ! (h) Right ascension (J2000) (1) integer*4 RAm_2 (nr__2) ! (min) Right ascension (J2000) real*4 RAs_2 (nr__2) ! (s) Right ascension (J2000) character*1 DE__2 (nr__2) ! Declination sign (J2000) integer*4 DEd_2 (nr__2) ! (deg) Declination (J2000) integer*4 DEm_2 (nr__2) ! (arcmin) Declination (J2000) real*4 DEs_2 (nr__2) ! (arcsec) Declination (J2000) character*24 SpType_2 (nr__2) ! Spectral type character*5 r_SpType_2 (nr__2) ! Reference for spectral type (6) character*9 FE_1 (nr__2) ! Foregroung or Background evidence (2) real*4 R_I_2 (nr__2) ! (mag) ? R-I colour index (3) real*4 Imag_2 (nr__2) ! (mag) ? I magnitude (4) real*4 I_Z_2 (nr__2) ! (mag) ? I-Z colour index (4) real*4 J_H_2 (nr__2) ! (mag) ? I-Z colour index (5) real*4 H_Ks_2 (nr__2) ! (mag) ? H-Ks colour index (5) real*4 Ksmag_1 (nr__2) ! (mag) ? Ks magnitude (5) *Note (1): The positions are from the 2MASS Point Source Catalog for * 18, 28, 206 and 370. All other coordinates are from the I-band images * in this work. *Note (2): Foreground nature is indicated by proper-motion * measurements ("pm"; Blaauw, 1952BAN....11..405B; Scholz et al., 1999, * Cat. <J/A+AS/137/305>) or a combination of no Li absorption ("Li"), * A_V_~0 ("A_V_"), and spectral features such as strong Na I and K I * ("NaK"; Luhman et al., 1998ApJ...508..347L; this work). *Note (3): From Herbig, 1998, Cat. <J/ApJ/497/736> for 77, 106 and 106 and * from Luhman 1999, Cat. <J/ApJ/525/466> for the remaining objects. *Note (4): From this work. *Note (5): From 2MASS for 14, 18, 27, 28, 34, 43, 54, 57, 73, 77, 80, 131, * 189, 320, 400, 1788, 131, 1932, 1941, 10226, 10252, and 10338 and from * Muench et al. (2003AJ....125.2029M) for the remaining objects. *Note (6): References are numbered as follows: * (1): Luhman et al., 1998ApJ...508..347L * (2): this work. * (3): measurements described by Luhman 1999, Cat. <J/ApJ/525/466> and * presented here for the first time; * (4): Fredrick, 1956AJ.....61..437F * (5): Herbig, 1998, Cat. <J/ApJ/497/736> c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Recent Identifications for Sources in IC 348 integer*4 nr__3 parameter (nr__3=360) ! Number of records character*40 ar__3 ! Full-size record integer*4 Seq_3 (nr__3) ! Identification number (1) character*1 m_Seq_1 (nr__3) ! Multiplicity indicator integer*4 H98 (nr__3) ! ? Herbig (1998, Cat. J/ApJ/497/736) ID number character*6 Naj00 (nr__3) ! Najita et al. (2000, Cat. J/ApJ/541/977) ID integer*4 Pre02 (nr__3) ! ? Preibisch & Zinnecker (2002AJ....123.1613P) ID character*17 v2MASS (nr__3) ! The 2MASS ID number *Note (1): From Luhman et al. 1998ApJ...508..347L) * Luhman (1999ApJ...525..466L), and this work. C============================================================================= C Loading file 'table2.dat' ! Data for Members of IC 348 C Format for file interpretation 1 format( + I5,A1,A1,1X,I2,1X,I2,1X,F5.2,A1,I2,1X,I2,1X,F4.1,1X,A32,1X,A5, + 1X,A5,1X,A10,1X,I5,1X,F4.2,1X,F9.5,1X,F4.2,1X,F5.2,1X,F4.2,1X, + F4.2,1X,F4.2,1X,F5.2,1X,A1,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,288 read(1,'(A146)')ar__ read(ar__,1) + Seq(i__),m_Seq(i__),f_Seq(i__),RAh(i__),RAm(i__),RAs(i__), + DE_(i__),DEd(i__),DEm(i__),DEs(i__),SpType(i__),r_SpType(i__), + Adopt(i__),Mem(i__),Teff(i__),AJ(i__),LBol(i__),R_I(i__), + Imag(i__),I_Z(i__),J_H(i__),H_Ks(i__),Kmag(i__),r_J_H(i__), + IMF(i__) if(ar__(88:92) .EQ. '') Teff(i__) = iNULL__ if(ar__(94:97) .EQ. '') AJ(i__) = rNULL__ if(ar__(99:107) .EQ. '') LBol(i__) = rNULL__ if(ar__(109:112) .EQ. '') R_I(i__) = rNULL__ if(ar__(114:118) .EQ. '') Imag(i__) = rNULL__ if(ar__(120:123) .EQ. '') I_Z(i__) = rNULL__ if(ar__(125:128) .EQ. '') J_H(i__) = rNULL__ if(ar__(130:133) .EQ. '') H_Ks(i__) = rNULL__ if(ar__(135:139) .EQ. '') Kmag(i__) = rNULL__ RAdeg(i__) = rNULL__ DEdeg(i__) = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15. if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4. if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240. if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__) if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60. if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,1) + Seq(i__),m_Seq(i__),f_Seq(i__),RAh(i__),RAm(i__),RAs(i__), + DE_(i__),DEd(i__),DEm(i__),DEs(i__),SpType(i__),r_SpType(i__), + Adopt(i__),Mem(i__),Teff(i__),AJ(i__),LBol(i__),R_I(i__), + Imag(i__),I_Z(i__),J_H(i__),H_Ks(i__),Kmag(i__),r_J_H(i__), + IMF(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Foreground Stars C Format for file interpretation 2 format( + I5,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A24,1X,A5,2X, + A9,2X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,20 read(1,'(A106)')ar__1 read(ar__1,2) + Seq_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__), + DEd_1(i__),DEm_1(i__),DEs_1(i__),SpType_1(i__), + r_SpType_1(i__),FE(i__),R_I_1(i__),Imag_1(i__),I_Z_1(i__), + J_H_1(i__),H_Ks_1(i__),Ksmag(i__) if(ar__1(74:77) .EQ. '') R_I_1(i__) = rNULL__ if(ar__1(79:83) .EQ. '') Imag_1(i__) = rNULL__ if(ar__1(85:88) .EQ. '') I_Z_1(i__) = rNULL__ if(ar__1(90:94) .EQ. '') J_H_1(i__) = rNULL__ if(ar__1(96:100) .EQ. '') H_Ks_1(i__) = rNULL__ if(ar__1(102:106) .EQ. '') Ksmag(i__) = rNULL__ RAdeg_1(i__) = rNULL__ DEdeg_1(i__) = 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(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15. if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4. if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240. if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__) if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60. if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600. if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__) c ..............Just test output........... write(6,2) + Seq_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__), + DEd_1(i__),DEm_1(i__),DEs_1(i__),SpType_1(i__), + r_SpType_1(i__),FE(i__),R_I_1(i__),Imag_1(i__),I_Z_1(i__), + J_H_1(i__),H_Ks_1(i__),Ksmag(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Background Stars C Format for file interpretation 3 format( + I5,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A24,1X,A5,2X, + A9,2X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,103 read(1,'(A107)')ar__2 read(ar__2,3) + Seq_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__), + DEd_2(i__),DEm_2(i__),DEs_2(i__),SpType_2(i__), + r_SpType_2(i__),FE_1(i__),R_I_2(i__),Imag_2(i__),I_Z_2(i__), + J_H_2(i__),H_Ks_2(i__),Ksmag_1(i__) if(ar__2(74:77) .EQ. '') R_I_2(i__) = rNULL__ if(ar__2(79:83) .EQ. '') Imag_2(i__) = rNULL__ if(ar__2(85:88) .EQ. '') I_Z_2(i__) = rNULL__ if(ar__2(90:94) .EQ. '') J_H_2(i__) = rNULL__ if(ar__2(96:100) .EQ. '') H_Ks_2(i__) = rNULL__ if(ar__2(102:106) .EQ. '') Ksmag_1(i__) = rNULL__ RAdeg_2(i__) = rNULL__ DEdeg_2(i__) = 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(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15. if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4. if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240. if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__) if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60. if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600. if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__) c ..............Just test output........... write(6,3) + Seq_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__), + DEd_2(i__),DEm_2(i__),DEs_2(i__),SpType_2(i__), + r_SpType_2(i__),FE_1(i__),R_I_2(i__),Imag_2(i__),I_Z_2(i__), + J_H_2(i__),H_Ks_2(i__),Ksmag_1(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table6.dat' ! Recent Identifications for Sources in IC 348 C Format for file interpretation 4 format(I5,A1,2X,I3,1X,A6,1X,I3,1X,A17) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,360 read(1,'(A40)')ar__3 read(ar__3,4) + Seq_3(i__),m_Seq_1(i__),H98(i__),Naj00(i__),Pre02(i__), + v2MASS(i__) if(ar__3(9:11) .EQ. '') H98(i__) = iNULL__ if(ar__3(20:22) .EQ. '') Pre02(i__) = iNULL__ c ..............Just test output........... write(6,4) + Seq_3(i__),m_Seq_1(i__),H98(i__),Naj00(i__),Pre02(i__), + v2MASS(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end