Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/527/A145 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/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 ! (deg) Right Ascension J2000 real*8 DEdeg ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Seq ! [1/84] Source sequential number (1) integer*4 RAh ! (h) Right ascension (J2000) integer*4 RAm ! (min) Right ascension (J2000) real*4 RAs ! (s) Right ascension (J2000) character*1 DE_ ! Declination sign (J2000) integer*4 DEd ! (deg) Declination (J2000) integer*4 DEm ! (arcmin) Declination (J2000) real*4 DEs ! (arcsec) Declination (J2000) real*4 Fp ! (Jy) Peak flux density at 870um (in Jy/21.2"-beam) real*4 Ftot ! (Jy) Total flux density at 870um real*4 amaj ! (arcsec) Major axis (FWHM) real*4 amin ! (arcsec) Minor axis (FWHM) real*4 PA ! (deg) [-90/90] Position angle (east from north) real*4 Size ! (arcsec) ? Mean source size * (geometrical mean of amaj and amin) character*2 Type ! Type of source based on the possible * associations found in the SIMBAD database (2) character*23 SIMBAD ! Source found in the SIMBAD database within the * FWHM ellipse character*3 Stype ! SIMBAD type of the SIMBAD source real*4 Dist ! (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 ! (deg) Right Ascension J2000 real*8 DEdeg_1 ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Seq_1 ! [1/84] Source sequential number (1) integer*4 RAh_1 ! (h) Right ascension (J2000) integer*4 RAm_1 ! (min) Right ascension (J2000) real*4 RAs_1 ! (s) Right ascension (J2000) character*1 DE__1 ! Declination sign (J2000) integer*4 DEd_1 ! (deg) Declination (J2000) integer*4 DEm_1 ! (arcmin) Declination (J2000) real*4 DEs_1 ! (arcsec) Declination (J2000) real*4 Fp_1 ! (Jy) Peak flux density at 870um (in Jy/21.2"-beam) real*4 Ftot_1 ! (Jy) Total flux density at 870um real*4 amaj_1 ! (arcsec) Major axis (FWHM) real*4 amin_1 ! (arcsec) Minor axis (FWHM) real*4 PA_1 ! (deg) [-90/90] Position angle (east from north) real*4 Size_1 ! (arcsec) ? Mean source size * (geometrical mean of amaj and amin) character*2 Type_1 ! Type of source based on the possible * associations found in the SIMBAD database (2) character*23 SIMBAD_1 ! Source found in the SIMBAD database within the * FWHM ellipse character*3 Stype_1 ! SIMBAD type of the SIMBAD source real*4 Dist_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 ! Name, Cha1-CNN integer*4 Seq_2 ! [1/84] Numbering of Gaussclumps sources like * in Table 2 real*4 FWHM1 ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHM2 ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHMa1 ! (kAU) Size of the fitted Gaussian (1) real*4 FWHMa2 ! (kAU) Size of the fitted Gaussian (1) real*4 Ra ! Aspect ratio (1) real*4 Npeak ! (10+21cm-2) Peak H2 column density (2) real*4 AV ! (mag) Visual extinction derived from 2MASS real*4 Mpeak ! (Msun) Mass in the central beam (HPBW = 21.2") (3) real*4 Mtot ! (Msun) Total mass derived from the Gaussian fit (3) real*4 M50_ ! (Msun) Mass computed from the flux measured in an * aperture of 50" in diameter (3) integer*4 CM ! (%) Mass concentration Mpeak/M50" integer*4 e_CM ! (%) Statistical rms uncertainty on CM real*4 npeak_1 ! (10+5cm-3) Beam-averaged free-particle density within * the central beam (4) real*4 n50_ ! (10+5cm-3) Mean free-particle density computed for the * mass M50" in the aperture of diameter * 50" (4) real*4 cn ! Density contrast npeak/n50" real*4 e_cn ! 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 ! Name, Cha1-SNN integer*4 Seq_3 ! [1/87] Numbering of Gaussclumps sources like * in table2 (table8) or table5 (table9) real*4 FWHM1_1 ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHM2_1 ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHMa1_1 ! (kAU) Deconvolved source size (1) real*4 FWHMa2_1 ! (kAU) Deconvolved source size (1) real*4 SNR ! Peak signal-to-noise ratio of the fitted * Gaussian real*4 Npeak_2 ! (10+21cm-2) Peak H_2_ column density (2) real*4 AV_1 ! (mag) Visual extinction derived from 2MASS real*4 Mpeak_1 ! (Msun) Mass in the central beam (HPBW=21.2") (3) real*4 Mtot_1 ! (Msun) Total mass derived from the Gaussian fit (3) real*4 M50__1 ! (Msun) Mass computed from the flux measured in an * aperture of 50" in diameter (3) integer*4 CM_1 ! (%) Mass concentration Mpeak/M50" integer*4 e_CM_1 ! (%) Statistical rms uncertainty on CM character*4 Class ! Infrared class of the nearest YSO listed in * column SIMBAD of table 2 (table8) or * table5 (table9) integer*4 Ref ! 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 ! Name, Cha1-SNN integer*4 Seq_4 ! [1/87] Numbering of Gaussclumps sources like * in table2 (table8) or table5 (table9) real*4 FWHM1_2 ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHM2_2 ! (arcsec) Size of the fitted Gaussian (1) real*4 FWHMa1_2 ! (kAU) Deconvolved source size (1) real*4 FWHMa2_2 ! (kAU) Deconvolved source size (1) real*4 SNR_1 ! Peak signal-to-noise ratio of the fitted * Gaussian real*4 Npeak_3 ! (10+21cm-2) Peak H_2_ column density (2) real*4 AV_2 ! (mag) Visual extinction derived from 2MASS real*4 Mpeak_2 ! (Msun) Mass in the central beam (HPBW=21.2") (3) real*4 Mtot_2 ! (Msun) Total mass derived from the Gaussian fit (3) real*4 M50__2 ! (Msun) Mass computed from the flux measured in an * aperture of 50" in diameter (3) integer*4 CM_2 ! (%) Mass concentration Mpeak/M50" integer*4 e_CM_2 ! (%) Statistical rms uncertainty on CM character*4 Class_1 ! Infrared class of the nearest YSO listed in * column SIMBAD of table 2 (table8) or * table5 (table9) integer*4 Ref_1 ! 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,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Fp,Ftot,amaj,amin,PA,Size, + Type,SIMBAD,Stype,Dist if(ar__(57:60) .EQ. '') Size = rNULL__ if(ar__(93:96) .EQ. '') Dist = rNULL__ RAdeg = rNULL__ DEdeg = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh .GT. -180) RAdeg=RAh*15. if(RAm .GT. -180) RAdeg=RAdeg+RAm/4. if(RAs .GT. -180) RAdeg=RAdeg+RAs/240. if(DEd .GE. 0) DEdeg=DEd if(DEm .GE. 0) DEdeg=DEdeg+DEm/60. if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600. if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg c ..............Just test output........... write(6,1) + Seq,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Fp,Ftot,amaj,amin,PA,Size, + Type,SIMBAD,Stype,Dist write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg 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,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Fp_1,Ftot_1, + amaj_1,amin_1,PA_1,Size_1,Type_1,SIMBAD_1,Stype_1,Dist_1 if(ar__1(57:60) .EQ. '') Size_1 = rNULL__ if(ar__1(93:96) .EQ. '') Dist_1 = rNULL__ RAdeg_1 = rNULL__ DEdeg_1 = 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 .GT. -180) RAdeg_1=RAh_1*15. if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4. if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240. if(DEd_1 .GE. 0) DEdeg_1=DEd_1 if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60. if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600. if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1 c ..............Just test output........... write(6,2) + Seq_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Fp_1,Ftot_1, + amaj_1,amin_1,PA_1,Size_1,Type_1,SIMBAD_1,Stype_1,Dist_1 write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1 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,Seq_2,FWHM1,FWHM2,FWHMa1,FWHMa2,Ra,Npeak,AV,Mpeak,Mtot, + M50_,CM,e_CM,npeak_1,n50_,cn,e_cn c ..............Just test output........... write(6,3) + Name,Seq_2,FWHM1,FWHM2,FWHMa1,FWHMa2,Ra,Npeak,AV,Mpeak,Mtot, + M50_,CM,e_CM,npeak_1,n50_,cn,e_cn 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,Seq_3,FWHM1_1,FWHM2_1,FWHMa1_1,FWHMa2_1,SNR,Npeak_2, + AV_1,Mpeak_1,Mtot_1,M50__1,CM_1,e_CM_1,Class,Ref c ..............Just test output........... write(6,4) + Name_1,Seq_3,FWHM1_1,FWHM2_1,FWHMa1_1,FWHMa2_1,SNR,Npeak_2, + AV_1,Mpeak_1,Mtot_1,M50__1,CM_1,e_CM_1,Class,Ref 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,Seq_4,FWHM1_2,FWHM2_2,FWHMa1_2,FWHMa2_2,SNR_1,Npeak_3, + AV_2,Mpeak_2,Mtot_2,M50__2,CM_2,e_CM_2,Class_1,Ref_1 c ..............Just test output........... write(6,5) + Name_2,Seq_4,FWHM1_2,FWHM2_2,FWHMa1_2,FWHMa2_2,SNR_1,Npeak_3, + AV_2,Mpeak_2,Mtot_2,M50__2,CM_2,e_CM_2,Class_1,Ref_1 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end