Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/527/A145 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/A+A/527/A145 Chamaeleon I 870um sources (Belloche+, 2011) *================================================================================ *The end of star formation in Chamaeleon I ? *A LABOCA census of starless and protostellar cores. * Belloche A., Schuller F., Parise B., Andre P., Hatchell J., Jorgensen J.K., * Bontemps S., Weiss A., Menten , K.M., Muders D. * <Astron. Astrophys. 527, A145 (2011)> * =2011A&A...527A.145B C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Sources extracted with Gaussclumps in the 870um continuum sum map of Cha I at scale 5, and possible associations found in the SIMBAD database integer*4 nr__ parameter (nr__=118) ! Number of records character*96 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__) ! [1/84] Source sequential number (1) integer*4 RAh (nr__) ! (h) Right ascension (J2000) integer*4 RAm (nr__) ! (min) Right ascension (J2000) real*4 RAs (nr__) ! (s) Right ascension (J2000) character*1 DE_ (nr__) ! Declination sign (J2000) integer*4 DEd (nr__) ! (deg) Declination (J2000) integer*4 DEm (nr__) ! (arcmin) Declination (J2000) real*4 DEs (nr__) ! (arcsec) Declination (J2000) real*4 Fp (nr__) ! (Jy) Peak flux density at 870um (in Jy/21.2"-beam) real*4 Ftot (nr__) ! (Jy) Total flux density at 870um real*4 amaj (nr__) ! (arcsec) Major axis (FWHM) real*4 amin (nr__) ! (arcsec) Minor axis (FWHM) real*4 PA (nr__) ! (deg) [-90/90] Position angle (east from north) real*4 Size (nr__) ! (arcsec) ? Mean source size * (geometrical mean of amaj and amin) character*2 Type (nr__) ! Type of source based on the possible * associations found in the SIMBAD database (2) character*23 SIMBAD (nr__) ! Source found in the SIMBAD database within the * FWHM ellipse character*3 Stype (nr__) ! SIMBAD type of the SIMBAD source real*4 Dist (nr__) ! (arcsec) ? Distance of the possibly associated SIMBAD * source to the fitted peak position of the * LABOCA source *Note (1): Sequential number 1-84 in table2, and * Additional sequential number 1-5 in table5. *Note (2): Type of source based on the associations found in the SIMBAD * database as follows: * C = starless core or Class 0 protostar * S = young stellar object (Class I or more evolved) * Sc = candidate association with a YSO * R = residual of a stronger nearby source * A = artefact c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Additional compact sources with SIMBAD association with a formal peak flux density above 3.5{sigma} in the 870um continuum map of Cha I filtered up to scale 5 integer*4 nr__1 parameter (nr__1=9) ! Number of records character*96 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) ! [1/84] Source sequential number (1) integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000) 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) real*4 Fp_1 (nr__1) ! (Jy) Peak flux density at 870um (in Jy/21.2"-beam) real*4 Ftot_1 (nr__1) ! (Jy) Total flux density at 870um real*4 amaj_1 (nr__1) ! (arcsec) Major axis (FWHM) real*4 amin_1 (nr__1) ! (arcsec) Minor axis (FWHM) real*4 PA_1 (nr__1) ! (deg) [-90/90] Position angle (east from north) real*4 Size_1 (nr__1) ! (arcsec) ? Mean source size * (geometrical mean of amaj and amin) character*2 Type_1 (nr__1) ! Type of source based on the possible * associations found in the SIMBAD database (2) character*23 SIMBAD_1 (nr__1) ! Source found in the SIMBAD database within the * FWHM ellipse character*3 Stype_1 (nr__1) ! SIMBAD type of the SIMBAD source real*4 Dist_1 (nr__1) ! (arcsec) ? Distance of the possibly associated SIMBAD * source to the fitted peak position of the * LABOCA source *Note (1): Sequential number 1-84 in table2, and * Additional sequential number 1-5 in table5. *Note (2): Type of source based on the associations found in the SIMBAD * database as follows: * C = starless core or Class 0 protostar * S = young stellar object (Class I or more evolved) * Sc = candidate association with a YSO * R = residual of a stronger nearby source * A = artefact c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Characteristics of starless (or Class 0) sources extracted with Gaussclumps in the 870um continuum map of Cha I filtered up to scale 5 integer*4 nr__2 parameter (nr__2=60) ! Number of records character*86 ar__2 ! Full-size record character*8 Name (nr__2) ! Name, Cha1-CNN integer*4 Seq_2 (nr__2) ! [1/84] Numbering of Gaussclumps sources like * in Table 2 real*4 FWHM1 (nr__2) ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHM2 (nr__2) ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHMa1 (nr__2) ! (kAU) Size of the fitted Gaussian (1) real*4 FWHMa2 (nr__2) ! (kAU) Size of the fitted Gaussian (1) real*4 Ra (nr__2) ! Aspect ratio (1) real*4 Npeak (nr__2) ! (10+21cm-2) Peak H2 column density (2) real*4 AV (nr__2) ! (mag) Visual extinction derived from 2MASS real*4 Mpeak (nr__2) ! (Msun) Mass in the central beam (HPBW = 21.2") (3) real*4 Mtot (nr__2) ! (Msun) Total mass derived from the Gaussian fit (3) real*4 M50_ (nr__2) ! (Msun) Mass computed from the flux measured in an * aperture of 50" in diameter (3) integer*4 CM (nr__2) ! (%) Mass concentration Mpeak/M50" integer*4 e_CM (nr__2) ! (%) Statistical rms uncertainty on CM real*4 npeak_1 (nr__2) ! (10+5cm-3) Beam-averaged free-particle density within * the central beam (4) real*4 n50_ (nr__2) ! (10+5cm-3) Mean free-particle density computed for the * mass M50" in the aperture of diameter * 50" (4) real*4 cn (nr__2) ! Density contrast npeak/n50" real*4 e_cn (nr__2) ! Statistical rms uncertainty on cn *Note (1): Sizes smaller than 25.4" were set to 25.4" to compute the deconvolved * sizes, in order to account for a fit inaccuracy corresponding to a 5{sigma} * detection in peak flux density. As a result, the minimum size that can be * measured is about 2100AU. The aspect ratio is the ratio of the deconvolved * sizes along the major and minor axes. *Note (2): Peak H_2_ column density computed assuming a dust opacity of * 0.01cm^2^/g. The statistical rms uncertainty is 1.1x10^21^cm^-2^. *Note (3): The statistical rms uncertainties of Mpeak and M50" are 0.006 and * 0.011M_{sun}_, respectively. *Note (4): The statistical rms uncertainties of npeak and n50" are * 5.6x10^4^ and 7.2x10^3^cm^-3^, respectively. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Characteristics of YSOs extracted with Gaussclumps in the 870um continuum map of Cha I filtered up to scale 5 integer*4 nr__3 parameter (nr__3=16) ! Number of records character*83 ar__3 ! Full-size record character*8 Name_1 (nr__3) ! Name, Cha1-SNN integer*4 Seq_3 (nr__3) ! [1/87] Numbering of Gaussclumps sources like * in table2 (table8) or table5 (table9) real*4 FWHM1_1 (nr__3) ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHM2_1 (nr__3) ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHMa1_1 (nr__3) ! (kAU) Deconvolved source size (1) real*4 FWHMa2_1 (nr__3) ! (kAU) Deconvolved source size (1) real*4 SNR (nr__3) ! Peak signal-to-noise ratio of the fitted * Gaussian real*4 Npeak_2 (nr__3) ! (10+21cm-2) Peak H_2_ column density (2) real*4 AV_1 (nr__3) ! (mag) Visual extinction derived from 2MASS real*4 Mpeak_1 (nr__3) ! (Msun) Mass in the central beam (HPBW=21.2") (3) real*4 Mtot_1 (nr__3) ! (Msun) Total mass derived from the Gaussian fit (3) real*4 M50__1 (nr__3) ! (Msun) Mass computed from the flux measured in an * aperture of 50" in diameter (3) integer*4 CM_1 (nr__3) ! (%) Mass concentration Mpeak/M50" integer*4 e_CM_1 (nr__3) ! (%) Statistical rms uncertainty on CM character*4 Class (nr__3) ! Infrared class of the nearest YSO listed in * column SIMBAD of table 2 (table8) or * table5 (table9) integer*4 Ref (nr__3) ! Reference (4) *Note (1): The deconvolved size is computed from the fitted size multiplied by * (1+1/SNR), with SNR the peak signal-to-noise ratio, and therefore is an * upper limit in most cases. A zero value means that the uncertainty was * still underestimated. *Note (2): Peak H_2_ column density computed assuming a dust temperature of 20K * and a dust opacity of 0.03cm^2^/g. The statistical rms uncertainty is * 1.5x10^20^cm^-2^. *Note (3): The statistical rms uncertainties of Mpeak and M50" are 0.0009 and * 0.0015M_{sun}_, respectively. *Note (4): References as follows: * 2 = Luhman et al. (2008, Cat. J/ApJ/675/1375) * 3 = Furlan et al. (2009, Cat. J/ApJ/703/1964) * 4 = Henning et al. (1993A&A...276..129H) * 5 = Luhman & Muench (2008, Cat. J/ApJ/684/654) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table9.dat' ! Characteristics of the additional compact sources listed in Table 5 integer*4 nr__4 parameter (nr__4=5) ! Number of records character*83 ar__4 ! Full-size record character*8 Name_2 (nr__4) ! Name, Cha1-SNN integer*4 Seq_4 (nr__4) ! [1/87] Numbering of Gaussclumps sources like * in table2 (table8) or table5 (table9) real*4 FWHM1_2 (nr__4) ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHM2_2 (nr__4) ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHMa1_2 (nr__4) ! (kAU) Deconvolved source size (1) real*4 FWHMa2_2 (nr__4) ! (kAU) Deconvolved source size (1) real*4 SNR_1 (nr__4) ! Peak signal-to-noise ratio of the fitted * Gaussian real*4 Npeak_3 (nr__4) ! (10+21cm-2) Peak H_2_ column density (2) real*4 AV_2 (nr__4) ! (mag) Visual extinction derived from 2MASS real*4 Mpeak_2 (nr__4) ! (Msun) Mass in the central beam (HPBW=21.2") (3) real*4 Mtot_2 (nr__4) ! (Msun) Total mass derived from the Gaussian fit (3) real*4 M50__2 (nr__4) ! (Msun) Mass computed from the flux measured in an * aperture of 50" in diameter (3) integer*4 CM_2 (nr__4) ! (%) Mass concentration Mpeak/M50" integer*4 e_CM_2 (nr__4) ! (%) Statistical rms uncertainty on CM character*4 Class_1 (nr__4) ! Infrared class of the nearest YSO listed in * column SIMBAD of table 2 (table8) or * table5 (table9) integer*4 Ref_1 (nr__4) ! Reference (4) *Note (1): The deconvolved size is computed from the fitted size multiplied by * (1+1/SNR), with SNR the peak signal-to-noise ratio, and therefore is an * upper limit in most cases. A zero value means that the uncertainty was * still underestimated. *Note (2): Peak H_2_ column density computed assuming a dust temperature of 20K * and a dust opacity of 0.03cm^2^/g. The statistical rms uncertainty is * 1.5x10^20^cm^-2^. *Note (3): The statistical rms uncertainties of Mpeak and M50" are 0.0009 and * 0.0015M_{sun}_, respectively. *Note (4): References as follows: * 2 = Luhman et al. (2008, Cat. J/ApJ/675/1375) * 3 = Furlan et al. (2009, Cat. J/ApJ/703/1964) * 4 = Henning et al. (1993A&A...276..129H) * 5 = Luhman & Muench (2008, Cat. J/ApJ/684/654) C============================================================================= C Loading file 'table2.dat' ! Sources extracted with Gaussclumps in the 870um * continuum sum map of Cha I at scale 5, and * possible associations found in the SIMBAD database C Format for file interpretation 1 format( + I2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.3,1X,F5.3, + 1X,F5.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,A2,1X,A23,1X,A3,1X,F4.1) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,118 read(1,'(A96)')ar__ read(ar__,1) + Seq(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),Fp(i__),Ftot(i__),amaj(i__),amin(i__), + PA(i__),Size(i__),Type(i__),SIMBAD(i__),Stype(i__),Dist(i__) if(ar__(57:60) .EQ. '') Size(i__) = rNULL__ if(ar__(93:96) .EQ. '') Dist(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__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),Fp(i__),Ftot(i__),amaj(i__),amin(i__), + PA(i__),Size(i__),Type(i__),SIMBAD(i__),Stype(i__),Dist(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 'table5.dat' ! Additional compact sources with SIMBAD association * with a formal peak flux density above 3.5{sigma} * in the 870um continuum map of Cha I filtered up * to scale 5 C Format for file interpretation 2 format( + I2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.3,1X,F5.3, + 1X,F5.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,A2,1X,A23,1X,A3,1X,F4.1) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,9 read(1,'(A96)')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__),Fp_1(i__),Ftot_1(i__), + amaj_1(i__),amin_1(i__),PA_1(i__),Size_1(i__),Type_1(i__), + SIMBAD_1(i__),Stype_1(i__),Dist_1(i__) if(ar__1(57:60) .EQ. '') Size_1(i__) = rNULL__ if(ar__1(93:96) .EQ. '') Dist_1(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__),Fp_1(i__),Ftot_1(i__), + amaj_1(i__),amin_1(i__),PA_1(i__),Size_1(i__),Type_1(i__), + SIMBAD_1(i__),Stype_1(i__),Dist_1(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 'table6.dat' ! Characteristics of starless (or Class 0) sources * extracted with Gaussclumps in the 870um * continuum map of Cha I filtered up to scale 5 C Format for file interpretation 3 format( + A8,1X,I2,1X,F5.1,1X,F4.1,1X,F4.1,1X,F3.1,1X,F3.1,1X,F4.1,1X, + F4.1,1X,F5.3,1X,F5.3,F5.3,1X,I2,1X,I2,1X,F4.1,1X,F4.1,1X,F3.1, + 1X,F3.1) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,60 read(1,'(A86)')ar__2 read(ar__2,3) + Name(i__),Seq_2(i__),FWHM1(i__),FWHM2(i__),FWHMa1(i__), + FWHMa2(i__),Ra(i__),Npeak(i__),AV(i__),Mpeak(i__),Mtot(i__), + M50_(i__),CM(i__),e_CM(i__),npeak_1(i__),n50_(i__),cn(i__), + e_cn(i__) c ..............Just test output........... write(6,3) + Name(i__),Seq_2(i__),FWHM1(i__),FWHM2(i__),FWHMa1(i__), + FWHMa2(i__),Ra(i__),Npeak(i__),AV(i__),Mpeak(i__),Mtot(i__), + M50_(i__),CM(i__),e_CM(i__),npeak_1(i__),n50_(i__),cn(i__), + e_cn(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! Characteristics of YSOs extracted with Gaussclumps * in the 870um continuum map of Cha I filtered * up to scale 5 C Format for file interpretation 4 format( + A8,1X,I2,1X,F4.1,1X,F4.1,1X,F4.1,1X,F3.1,1X,F5.1,1X,F5.2,1X, + F4.1,1X,F6.4,1X,F6.4,1X,F6.4,1X,I3,1X,I3,1X,A4,1X,I1) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,16 read(1,'(A83)')ar__3 read(ar__3,4) + Name_1(i__),Seq_3(i__),FWHM1_1(i__),FWHM2_1(i__), + FWHMa1_1(i__),FWHMa2_1(i__),SNR(i__),Npeak_2(i__),AV_1(i__), + Mpeak_1(i__),Mtot_1(i__),M50__1(i__),CM_1(i__),e_CM_1(i__), + Class(i__),Ref(i__) c ..............Just test output........... write(6,4) + Name_1(i__),Seq_3(i__),FWHM1_1(i__),FWHM2_1(i__), + FWHMa1_1(i__),FWHMa2_1(i__),SNR(i__),Npeak_2(i__),AV_1(i__), + Mpeak_1(i__),Mtot_1(i__),M50__1(i__),CM_1(i__),e_CM_1(i__), + Class(i__),Ref(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table9.dat' ! Characteristics of the additional compact sources * listed in Table 5 C Format for file interpretation 5 format( + A8,1X,I2,1X,F4.1,1X,F4.1,1X,F4.1,1X,F3.1,1X,F5.1,1X,F5.2,1X, + F4.1,1X,F6.4,1X,F6.4,1X,F6.4,1X,I3,1X,I3,1X,A4,1X,I1) C Effective file loading open(unit=1,status='old',file= +'table9.dat') write(6,*) '....Loading file: table9.dat' do i__=1,5 read(1,'(A83)')ar__4 read(ar__4,5) + Name_2(i__),Seq_4(i__),FWHM1_2(i__),FWHM2_2(i__), + FWHMa1_2(i__),FWHMa2_2(i__),SNR_1(i__),Npeak_3(i__),AV_2(i__), + Mpeak_2(i__),Mtot_2(i__),M50__2(i__),CM_2(i__),e_CM_2(i__), + Class_1(i__),Ref_1(i__) c ..............Just test output........... write(6,5) + Name_2(i__),Seq_4(i__),FWHM1_2(i__),FWHM2_2(i__), + FWHMa1_2(i__),FWHMa2_2(i__),SNR_1(i__),Npeak_3(i__),AV_2(i__), + Mpeak_2(i__),Mtot_2(i__),M50__2(i__),CM_2(i__),e_CM_2(i__), + Class_1(i__),Ref_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end