Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/138/1116 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-29 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/AJ/138/1116 A Spitzer view of NGC 2264 (Sung+, 2009) *================================================================================ *A Spitzer view of the young open cluster NGC 2264. * Sung H., Stauffer J.R., Bessell M.S. * <Astron. J., 138, 1116-1136 (2009)> * =2009AJ....138.1116S C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Catalog of sources detected with Spitzer Space Telescope IRAC and MIPS 24 micron integer*4 nr__ parameter (nr__=21991) ! Number of records character*165 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 SST ! Running identification number (G1) integer*4 RAh ! (h) Hour of Right Ascension (J2000) integer*4 RAm ! (min) Minute of Right Ascension (J2000) real*4 RAs ! (s) Second of Right Ascension (J2000) character*1 DE_ ! Sign of Declination (J2000) integer*4 DEd ! (deg) Degree of Declination (J2000) integer*4 DEm ! (arcmin) Arcminute of Declination (J2000) real*4 DEs ! (arcsec) Arcsecond of Declination (J2000) real*4 v_3_6_ ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude real*4 v_4_5_ ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude real*4 v_5_8_ ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude real*4 v_8_0_ ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude real*4 v_24_ ! (mag) ? Spitzer/MIPS 24 micron band magnitude real*4 e__3_6_ ! (mag) ? Uncertainty in [3.6] real*4 e__4_5_ ! (mag) ? Uncertainty in [4.5] real*4 e__5_8_ ! (mag) ? Uncertainty in [5.8] real*4 e__8_0_ ! (mag) ? Uncertainty in [8.0] real*4 e__24_ ! (mag) ? Uncertainty in [24] integer*4 o__3_6_ ! Number of [3.6] observations integer*4 o__4_5_ ! Number of [4.5] observations integer*4 o__5_8_ ! Number of [5.8] observations integer*4 o__8_0_ ! Number of [8.0] observations integer*4 o__24_ ! Number of [24] observations character*1 Dup ! [D/T] Duplicity code (1) character*8 Mem ! Membership (2) character*16 v2MASS ! 2MASS identification (HHMMSSss+DDMMSSs) character*20 Sung ! Sung et al. identification(s) (G2) character*13 Note ! Counterparts of weak X-ray source (3) *Note (1): Code as follows: * D = SST source having two optical counterparts within 2" searching radii; * T = SST source having three optical counterparts within 2" searching radii. *Note (2): IR classification (Class I, II, II/III, pre-TD (pre-transition disk), * and TD(transition disk)) and other membership listed in * Sung et al. (2008, Cat. J/AJ/135/441) as follows: * H = H{alpha} emission star, * h = H{alpha} emission candidate; * X = X-ray emission star, * x = X-ray emission candidate; * + = H{alpha} emission star with X-ray emission, * - = H{alpha} emission candidate with X-ray emission, * E = strong H{alpha} emission star from spectroscopy (EWHa>=10{AA} in Dahm & * Simon (2005, Cat. J/AJ/129/829) or H{alpha} index >=2.5 in Reipurth et al. * (2004, Cat. J/AJ/127/1117)), * e = weak H{alpha} emission star from spectroscopy (EWHa<10{AA} in Dahm & Simon * (2005, Cat. J/AJ/129/829) or H{alpha} index <=2 in Reipurth et al. (2004, * Cat. J/AJ/127/1117)); * M = X-ray emission star with membership class "e", i.e. "M"="e"+"X". * P = X-ray emission star with membership class "E", i.e. "P"="E"+"X"; * p = X-ray candidate with membership class "E", i.e. "p"="E"+"x". *Note (3): 30 IRAC sources as counterparts of weak X-ray sources with * significance between 3.3 and 5.0. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Spitzer counterparts of X-ray sources without optical or 2MASS counterparts integer*4 nr__1 parameter (nr__1=43) ! Number of records character*117 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 SST_1 ! Running identification number (G1) integer*4 RAh_1 ! (h) Hour of Right Ascension (J2000) integer*4 RAm_1 ! (min) Minute of Right Ascension (J2000) real*4 RAs_1 ! (s) Second of Right Ascension (J2000) character*1 DE__1 ! Sign of declination (J2000) integer*4 DEd_1 ! (deg) Degree of Declination (J2000) integer*4 DEm_1 ! (arcmin) Arcminute of Declination (J2000) real*4 DEs_1 ! (arcsec) Arcsecond of Declination (J2000) real*4 v_3_6__1 ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude real*4 v_4_5__1 ! (mag) Spitzer/IRAC 4.5 micron band magnitude real*4 v_5_8__1 ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude real*4 v_8_0__1 ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude real*4 v_24__1 ! (mag) ? Spitzer/MIPS 24 micron band magnitude real*4 e__3_6__1 ! (mag) ? Uncertainty in [3.6] real*4 e__4_5__1 ! (mag) Uncertainty in [4.5] real*4 e__5_8__1 ! (mag) ? Uncertainty in [5.8] real*4 e__8_0__1 ! (mag) ? Uncertainty in [8.0] real*4 e__24__1 ! (mag) ? Uncertainty in [24] integer*4 o__3_6__1 ! Number of [3.6] observations integer*4 o__4_5__1 ! Number of [4.5] observations integer*4 o__5_8__1 ! Number of [5.8] observations integer*4 o__8_0__1 ! Number of [8.0] observations integer*4 o__24__1 ! Number of [24] observations character*6 Class ! Object classification (Galaxy, I, II or II/III) real*4 Signi ! X-ray detection significance (1) real*4 Dist ! (arcsec) Distance Spitzer/X-ray *Note (1): X-ray detection significance from "pwdetect" (see Sung et al. 2004, * Cat. J/AJ/128/1684; Flaccomio et al. 2006, Cat. J/A+A/455/903). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Variables detected in IRAC [3.6] and [4.5] integer*4 nr__2 parameter (nr__2=28) ! Number of records character*102 ar__2 ! Full-size record integer*4 SST_2 ! Running identification number (G1) character*1 f_SST ! [a] long exposure images used (1) real*4 v_3_6_6 ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude on * 2004 Mar 6 real*4 e__3_6_6 ! (mag) ? Uncertainty in [3.6] real*4 v_3_6_8 ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude on * 2004 Oct 8 real*4 e__3_6_8 ! (mag) ? Uncertainty in [3.6] real*4 d_3_6_ ! (mag) ? [3.6] difference real*4 v_4_5_6 ! (mag) ? Spitzer/IRAC [4.5 micron band magnitude on * 2004 Mar 6 real*4 e__4_5_6 ! (mag) ? Uncertainty in [4.5] real*4 v_4_5_8 ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude on * 2004 Oct 8 real*4 e__4_5_8 ! (mag) ? Uncertainty in [4.5] real*4 d_4_5_ ! (mag) ? [4.5] difference character*6 Class_1 ! Infra-red classification (I, II or II/III) character*8 Other ! Other classification (X or X+Halpha) character*13 OName ! Sung et al. identification(s) (G2) *Note (1): * a = Used only data from long exposure images because short and long exposure * images show large systematic differences probably due to the nearby * bright star SST 13808. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! New BMS star candidates (stars below the pre-main-sequence locus) integer*4 nr__3 parameter (nr__3=10) ! Number of records character*82 ar__3 ! Full-size record integer*4 SST_3 ! Running identification number (G1) character*1 n_SST ! [bc] Note on source (1) character*13 OName_1 ! Sung et al. identification(s) (G2) real*4 v_3_6__2 ! (mag) Spitzer/IRAC 3.6 micron band magnitude real*4 v_4_5__2 ! (mag) Spitzer/IRAC 4.5 micron band magnitude real*4 v_5_8__2 ! (mag) Spitzer/IRAC 5.8 micron band magnitude real*4 v_8_0__2 ! (mag) Spitzer/IRAC 8.0 micron band magnitude real*4 v_24__2 ! (mag) ? Spitzer/MIPS 24 micron band magnitude character*3 Ha_ ! [Yes/No/? ] H{alpha} emission? character*2 Class_2 ! IR classification (I or II) real*4 SED ! SED slope real*4 e_SED ! SED uncertainty real*4 Qcc ! Mean value of a YSO class Q_CC_ (2) real*4 D_I ! Magnitude difference from the faint limit of * the PMS locus in NGC 2264 *Note (1): Flag as follows: * b = If the actual optical counterpart of SST 10184 is C30962, then SST 10184 * may be a SFG with PAH emission. See Section 2.5. * c = Ogura (1984PASJ...36..139O) noted as an H{alpha} emission star, but no * signature of H{alpha} emission from CCD photometry was detected. *Note (2): The mean value of a YSO class (Cl_i_) (call this quantity Q_CC_) is: * Q_CC_={Sigma}_i_Cl_i_xq_i_xw_i_/{Sigma}_i_q_i_xw_i_ * where q_i_ and w_i_ are different weights. See section 3.1. C============================================================================= C Loading file 'table1.dat' ! Catalog of sources detected with Spitzer Space * Telescope IRAC and MIPS 24 micron C Format for file interpretation 1 format( + I5,1X,I1,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F6.3,1X,F6.3, + 1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X, + F5.3,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,A1,1X,A8,1X,A16,1X,A20, + 1X,A13) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,21991 read(1,'(A165)')ar__ read(ar__,1) + SST,RAh,RAm,RAs,DE_,DEd,DEm,DEs,v_3_6_,v_4_5_,v_5_8_,v_8_0_, + v_24_,e__3_6_,e__4_5_,e__5_8_,e__8_0_,e__24_,o__3_6_,o__4_5_, + o__5_8_,o__8_0_,o__24_,Dup,Mem,v2MASS,Sung,Note if(ar__(30:35) .EQ. '') v_3_6_ = rNULL__ if(ar__(37:42) .EQ. '') v_4_5_ = rNULL__ if(ar__(44:49) .EQ. '') v_5_8_ = rNULL__ if(ar__(51:56) .EQ. '') v_8_0_ = rNULL__ if(ar__(58:62) .EQ. '') v_24_ = rNULL__ if(ar__(64:68) .EQ. '') e__3_6_ = rNULL__ if(ar__(70:74) .EQ. '') e__4_5_ = rNULL__ if(ar__(76:80) .EQ. '') e__5_8_ = rNULL__ if(ar__(82:86) .EQ. '') e__8_0_ = rNULL__ if(ar__(88:92) .EQ. '') e__24_ = 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) + SST,RAh,RAm,RAs,DE_,DEd,DEm,DEs,v_3_6_,v_4_5_,v_5_8_,v_8_0_, + v_24_,e__3_6_,e__4_5_,e__5_8_,e__8_0_,e__24_,o__3_6_,o__4_5_, + o__5_8_,o__8_0_,o__24_,Dup,Mem,v2MASS,Sung,Note write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Spitzer counterparts of X-ray sources without * optical or 2MASS counterparts C Format for file interpretation 2 format( + I5,1X,I1,1X,I2,1X,F5.2,1X,A1,I1,1X,I2,1X,F4.1,1X,F6.3,1X,F6.3, + 1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X, + F5.3,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,A6,1X,F4.1,1X,F3.1) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,43 read(1,'(A117)')ar__1 read(ar__1,2) + SST_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,v_3_6__1, + v_4_5__1,v_5_8__1,v_8_0__1,v_24__1,e__3_6__1,e__4_5__1, + e__5_8__1,e__8_0__1,e__24__1,o__3_6__1,o__4_5__1,o__5_8__1, + o__8_0__1,o__24__1,Class,Signi,Dist if(ar__1(29:34) .EQ. '') v_3_6__1 = rNULL__ if(ar__1(43:48) .EQ. '') v_5_8__1 = rNULL__ if(ar__1(50:55) .EQ. '') v_8_0__1 = rNULL__ if(ar__1(57:61) .EQ. '') v_24__1 = rNULL__ if(ar__1(63:67) .EQ. '') e__3_6__1 = rNULL__ if(ar__1(75:79) .EQ. '') e__5_8__1 = rNULL__ if(ar__1(81:85) .EQ. '') e__8_0__1 = rNULL__ if(ar__1(87:91) .EQ. '') e__24__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) + SST_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,v_3_6__1, + v_4_5__1,v_5_8__1,v_8_0__1,v_24__1,e__3_6__1,e__4_5__1, + e__5_8__1,e__8_0__1,e__24__1,o__3_6__1,o__4_5__1,o__5_8__1, + o__8_0__1,o__24__1,Class,Signi,Dist write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Variables detected in IRAC [3.6] and [4.5] C Format for file interpretation 3 format( + I5,A1,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F6.3,1X,F5.3, + 1X,F6.3,1X,F5.3,1X,F6.3,1X,A6,1X,A8,1X,A13) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,28 read(1,'(A102)')ar__2 read(ar__2,3) + SST_2,f_SST,v_3_6_6,e__3_6_6,v_3_6_8,e__3_6_8,d_3_6_,v_4_5_6, + e__4_5_6,v_4_5_8,e__4_5_8,d_4_5_,Class_1,Other,OName if(ar__2(8:13) .EQ. '') v_3_6_6 = rNULL__ if(ar__2(15:19) .EQ. '') e__3_6_6 = rNULL__ if(ar__2(21:26) .EQ. '') v_3_6_8 = rNULL__ if(ar__2(28:32) .EQ. '') e__3_6_8 = rNULL__ if(ar__2(34:39) .EQ. '') d_3_6_ = rNULL__ if(ar__2(41:46) .EQ. '') v_4_5_6 = rNULL__ if(ar__2(48:52) .EQ. '') e__4_5_6 = rNULL__ if(ar__2(54:59) .EQ. '') v_4_5_8 = rNULL__ if(ar__2(61:65) .EQ. '') e__4_5_8 = rNULL__ if(ar__2(67:72) .EQ. '') d_4_5_ = rNULL__ c ..............Just test output........... write(6,3) + SST_2,f_SST,v_3_6_6,e__3_6_6,v_3_6_8,e__3_6_8,d_3_6_,v_4_5_6, + e__4_5_6,v_4_5_8,e__4_5_8,d_4_5_,Class_1,Other,OName c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! New BMS star candidates (stars below the * pre-main-sequence locus) C Format for file interpretation 4 format( + I5,A1,1X,A13,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,A3,1X, + A2,1X,F5.2,1X,F4.2,1X,F3.1,1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,10 read(1,'(A82)')ar__3 read(ar__3,4) + SST_3,n_SST,OName_1,v_3_6__2,v_4_5__2,v_5_8__2,v_8_0__2, + v_24__2,Ha_,Class_2,SED,e_SED,Qcc,D_I if(ar__3(50:54) .EQ. '') v_24__2 = rNULL__ c ..............Just test output........... write(6,4) + SST_3,n_SST,OName_1,v_3_6__2,v_4_5__2,v_5_8__2,v_8_0__2, + v_24__2,Ha_,Class_2,SED,e_SED,Qcc,D_I c .......End.of.Just test output........... end do close(1) C============================================================================= stop end