Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/715/671 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/ApJ/715/671 X-ray and IR emission from YSOs (Wolk+, 2010) *================================================================================ *X-ray and infrared emission from young stellar objects near LkH{alpha} 101. * Wolk S.J., Winston E., Bourke T.L., Gutermuth R., Megeath S.T., * Spitzbart B.D., Osten R. * <Astrophys. J., 715, 671-695 (2010)> * =2010ApJ...715..671W C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! X-ray sources detected in the LkH{alpha} 101 field integer*4 nr__ parameter (nr__=213) ! Number of records character*41 ar__ ! Full-size record integer*4 Seq ! [2/231] Running source number (skips indicate * the spurious detections) character*18 ID ! Source identification; (HHMMSS.ss+DDMMSS.s; * <LkHa101CXO JHHMMSS.ss+DDMMSS.s> in Simbad) real*4 Off ! (arcsec) Offset from axis integer*4 RCts ! (ct) Raw counts real*4 NCts ! (ct) Net counts c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! 2MASS photometry of X-ray sources integer*4 nr__1 parameter (nr__1=140) ! Number of records character*51 ar__1 ! Full-size record integer*4 Seq_1 ! [2/231] Source index number real*4 X_2M ! (arcsec) X-ray/2MASS positional offset real*4 Off2M ! (arcsec) Offset from axis real*4 Jmag ! (mag) 2MASS J band magnitude real*4 e_Jmag ! (mag) ? Uncertainty in Jmag real*4 Hmag ! (mag) 2MASS H band magnitude real*4 e_Hmag ! (mag) ? Uncertainty in Hmag real*4 Kmag ! (mag) 2MASS K band magnitude real*4 e_Kmag ! (mag) ? Uncertainty in Kmag character*3 Qual ! 2MASS photometry quality codes c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Spitzer photometry of X-ray sources integer*4 nr__2 parameter (nr__2=158) ! Number of records character*89 ar__2 ! Full-size record integer*4 Seq_2 ! [2/231] Source index number (1) real*4 X_Sptz ! (arcsec) X-ray/Spitzer positional offset real*4 Off_1 ! (arcsec) Offset from axis 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] real*4 AK ! (mag) The K band extinction character*14 Class ! Source classification (1, 2, III, T(transition), * gal (background galaxy), or LkH{alpha} 101) character*1 f_Class ! [b] b: star within 90" of LkH{alpha} 101 (2) *Note (1): In addition to their IAU approved Names: * <LkHa101CXO JHHMMSS.ss+DDMMSS.s> given in Table 1, these sources are * equally identified as: <LkHa101SST JHHMMSS.ss+DDMMSS.s>. *Note (2): Stars with a notation "b" are within 90 arcsec of LkH{alpha} 101. * Their proximity, led to the conclusion that they were cluster * members, despite missing some IRAC data. For these stars 2MASS data * played an enhanced role in determining Class identification. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! LkH{alpha} 101 cluster spectral classification integer*4 nr__3 parameter (nr__3=87) ! Number of records character*57 ar__3 ! Full-size record integer*4 Seq_3 ! [2/231]? Chandra identification number character*18 ID_1 ! LkHaSST identification (HHMMSS.ss+DDMMSS.s; * <LkHa101SST JHHMMSS.ss+DDMMSS.s> in Simbad) character*4 SpT ! MK spectral type real*4 e_SpT ! SpT uncertainty integer*4 Teff ! (K) Effective temperature real*4 logL ! ([-]) Overall luminosity log(L_bol_/L_0_) real*8 EWHa ! (0.1nm) H{alpha} (6563{AA}) equivalent width real*4 EWLi ! (0.1nm) ? LiI (6707{AA}) equivalent width c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Spectral fits to X-ray sources integer*4 nr__4 parameter (nr__4=88) ! Number of records character*64 ar__4 ! Full-size record integer*4 Seq_4 ! [2/231] Source index number real*4 NCts_1 ! (ct) Net counts real*4 NH ! (10+22/cm2) Hydrogen column density real*4 e_NH ! (10+22/cm2) ? The 1{sigma} uncertainty in NH character*1 l_kT ! Limit flag on kT real*4 kT ! (keV) Effective temperature real*4 e_kT ! (keV) ? The 1{sigma} uncertainty in kT real*4 logF_A ! ([mW/m2]) Log of Absorbed flux; erg/s/cm^2^ real*4 logF_U ! ([mW/m2]) Log of Unabsorbed flux; erg/s/cm^2^ real*4 chi2 ! The reduced {chi}^2^ of the fit (1) character*9 Model ! Model used in fit (2) *Note (1): Included for "cstat" models for completeness. One model was * calculated with C statistics, it is unclear how this relates to * goodness of fit. The formal errors are unbiased however. *Note (2): The model used for the fit is coded in two parts. The left hand side * indicates the type of model used. The right hand side indicates the * data set. The overall model is written: LHS_RHS. * LHS as follows: * cstat = C-statistics used for the fit with an initial guess of kT=1keV and * NH=10^21^cm^-2^. * crs = c-stat used for initial guess with those input parameters and then * refit using Raymond-Smith plasma for final fit. * Mekal = again c-stat used for initial guess and then Mekal for final fit, * APEC = cstat used for initial guess and then APEC for final fit. * RHS as follows: * 29 = only data from first observation (Obsid = 5429) used in fit, * 28 = only data from second observation (Obsid = 5428) used in fit, * all = data from both observations used in fit. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Spectral fits to X-ray sources (lower quality) integer*4 nr__5 parameter (nr__5=30) ! Number of records character*64 ar__5 ! Full-size record integer*4 Seq_5 ! [2/231] Source index number real*4 NCts_2 ! (ct) Net counts real*4 NH_1 ! (10+22/cm2) Hydrogen column density real*4 e_NH_1 ! (10+22/cm2) ? The 1{sigma} uncertainty in NH character*1 l_kT_1 ! Limit flag on kT real*4 kT_1 ! (keV) Effective temperature real*4 e_kT_1 ! (keV) ? The 1{sigma} uncertainty in kT real*4 logF_A_1 ! ([mW/m2]) Log of Absorbed flux; erg/s/cm^2^ real*4 logF_U_1 ! ([mW/m2]) Log of Unabsorbed flux; erg/s/cm^2^ real*4 chi2_1 ! The reduced {chi}^2^ of the fit (1) character*9 Model_1 ! Model used in fit (2) *Note (1): Included for "cstat" models for completeness. One model was * calculated with C statistics, it is unclear how this relates to * goodness of fit. The formal errors are unbiased however. *Note (2): The model used for the fit is coded in two parts. The left hand side * indicates the type of model used. The right hand side indicates the * data set. The overall model is written: LHS_RHS. * LHS as follows: * cstat = C-statistics used for the fit with an initial guess of kT=1keV and * NH=10^21^cm^-2^. * crs = c-stat used for initial guess with those input parameters and then * refit using Raymond-Smith plasma for final fit. * Mekal = again c-stat used for initial guess and then Mekal for final fit, * APEC = cstat used for initial guess and then APEC for final fit. * RHS as follows: * 29 = only data from first observation (Obsid = 5429) used in fit, * 28 = only data from second observation (Obsid = 5428) used in fit, * all = data from both observations used in fit. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table10.dat' ! PMS stars detected by Spitzer but not detected in X-rays (2MASS colors) integer*4 nr__6 parameter (nr__6=69) ! Number of records character*55 ar__6 ! Full-size record character*18 ID_2 ! Source identification; (HHMMSS.ss+DDMMSS.s; * <LkHa101SST JHHMMSS.ss+DDMMSS.s> in Simbad) (1) real*4 Jmag_1 ! (mag) ? 2MASS J band magnitude real*4 e_Jmag_1 ! (mag) ? Uncertainty in Jmag real*4 Hmag_1 ! (mag) ? 2MASS H band magnitude real*4 e_Hmag_1 ! (mag) ? Uncertainty in Hmag real*4 Ksmag ! (mag) ? 2MASS Ks band magnitude real*4 e_Ksmag ! (mag) ? Uncertainty in Kmag character*3 Qual_1 ! 2MASS photometry quality codes *Note (1): Source 042955.68+352021.2 has been modified at CDS into * 042955.68+352021.0 as in table 11 and 2MASS catalog. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table11.dat' ! PMS stars not detected in X-rays (Spitzer) integer*4 nr__7 parameter (nr__7=69) ! Number of records character*81 ar__7 ! Full-size record character*18 ID_3 ! Source identification; (HHMMSS.ss+DDMMSS.s; * <LkHa101SST JHHMMSS.ss+DDMMSS.s> in Simbad) character*1 f_ID ! [*] outside Chandra field of view (1) real*4 v_3_6__1 ! (mag) Spitzer/IRAC 3.6 micron band magnitude real*4 e__3_6__1 ! (mag) Uncertainty in [3.6] real*4 v_4_5__1 ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude real*4 e__4_5__1 ! (mag) ? Uncertainty in [4.5] real*4 v_5_6_ ! (mag) ? Spitzer/IRAC 5.6 micron band magnitude real*4 e__5_6_ ! (mag) ? Uncertainty in [5.6] real*4 v_8_0__1 ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude real*4 e__8_0__1 ! (mag) ? Uncertainty in [8.0] real*4 v_24__1 ! (mag) ? Spitzer/MIPS 24 micron band magnitude real*4 e__24__1 ! (mag) ? Uncertainty in [24] real*4 AK_1 ! (mag) The K band extinction character*1 Cl ! [012T] Classification (Class 0, I, II or * T=possible transition objects) *Note (1): * * = These Spitzer identified PMS stars are off the Chandra field of view. C============================================================================= C Loading file 'table1.dat' ! X-ray sources detected in the LkH{alpha} 101 field C Format for file interpretation 1 format(I3,1X,A18,1X,F6.1,1X,I4,1X,F6.1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,213 read(1,'(A41)')ar__ read(ar__,1)Seq,ID,Off,RCts,NCts c ..............Just test output........... write(6,1)Seq,ID,Off,RCts,NCts c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! 2MASS photometry of X-ray sources C Format for file interpretation 2 format( + I3,1X,F4.2,1X,F5.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X, + F4.2,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,140 read(1,'(A51)')ar__1 read(ar__1,2) + Seq_1,X_2M,Off2M,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,e_Kmag,Qual if(ar__1(22:25) .EQ. '') e_Jmag = rNULL__ if(ar__1(33:36) .EQ. '') e_Hmag = rNULL__ if(ar__1(44:47) .EQ. '') e_Kmag = rNULL__ c ..............Just test output........... write(6,2) + Seq_1,X_2M,Off2M,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,e_Kmag,Qual c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Spitzer photometry of X-ray sources C Format for file interpretation 3 format( + I3,1X,F4.2,1X,F5.1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X, + F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A14,A1) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,158 read(1,'(A89)')ar__2 read(ar__2,3) + Seq_2,X_Sptz,Off_1,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_,AK,Class,f_Class if(ar__2(16:20) .EQ. '') v_3_6_ = rNULL__ if(ar__2(22:26) .EQ. '') v_4_5_ = rNULL__ if(ar__2(28:32) .EQ. '') v_5_8_ = rNULL__ if(ar__2(34:38) .EQ. '') v_8_0_ = rNULL__ if(ar__2(40:43) .EQ. '') v_24_ = rNULL__ if(ar__2(45:48) .EQ. '') e__3_6_ = rNULL__ if(ar__2(50:53) .EQ. '') e__4_5_ = rNULL__ if(ar__2(55:58) .EQ. '') e__5_8_ = rNULL__ if(ar__2(60:63) .EQ. '') e__8_0_ = rNULL__ if(ar__2(65:68) .EQ. '') e__24_ = rNULL__ c ..............Just test output........... write(6,3) + Seq_2,X_Sptz,Off_1,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_,AK,Class,f_Class c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! LkH{alpha} 101 cluster spectral classification C Format for file interpretation 4 format(I3,1X,A18,1X,A4,1X,F3.1,1X,I5,1X,F6.3,1X,F8.2,1X,F3.1) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,87 read(1,'(A57)')ar__3 read(ar__3,4)Seq_3,ID_1,SpT,e_SpT,Teff,logL,EWHa,EWLi if(ar__3(1:3) .EQ. '') Seq_3 = iNULL__ if(ar__3(55:57) .EQ. '') EWLi = rNULL__ c ..............Just test output........... write(6,4)Seq_3,ID_1,SpT,e_SpT,Teff,logL,EWHa,EWLi c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! Spectral fits to X-ray sources C Format for file interpretation 5 format( + I3,1X,F6.1,1X,F5.2,1X,F4.2,1X,A1,1X,F6.2,1X,F4.2,1X,F6.2,1X, + F6.2,1X,F4.2,1X,A9) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,88 read(1,'(A64)')ar__4 read(ar__4,5) + Seq_4,NCts_1,NH,e_NH,l_kT,kT,e_kT,logF_A,logF_U,chi2,Model if(ar__4(18:21) .EQ. '') e_NH = rNULL__ if(ar__4(32:35) .EQ. '') e_kT = rNULL__ c ..............Just test output........... write(6,5) + Seq_4,NCts_1,NH,e_NH,l_kT,kT,e_kT,logF_A,logF_U,chi2,Model c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! Spectral fits to X-ray sources (lower quality) C Format for file interpretation 6 format( + I3,1X,F6.1,1X,F5.2,1X,F4.2,1X,A1,1X,F6.2,1X,F4.2,1X,F6.2,1X, + F6.2,1X,F4.2,1X,A9) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,30 read(1,'(A64)')ar__5 read(ar__5,6) + Seq_5,NCts_2,NH_1,e_NH_1,l_kT_1,kT_1,e_kT_1,logF_A_1,logF_U_1, + chi2_1,Model_1 if(ar__5(18:21) .EQ. '') e_NH_1 = rNULL__ if(ar__5(32:35) .EQ. '') e_kT_1 = rNULL__ c ..............Just test output........... write(6,6) + Seq_5,NCts_2,NH_1,e_NH_1,l_kT_1,kT_1,e_kT_1,logF_A_1,logF_U_1, + chi2_1,Model_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table10.dat' ! PMS stars detected by Spitzer but not detected in * X-rays (2MASS colors) C Format for file interpretation 7 format( + A18,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table10.dat') write(6,*) '....Loading file: table10.dat' do i__=1,69 read(1,'(A55)')ar__6 read(ar__6,7) + ID_2,Jmag_1,e_Jmag_1,Hmag_1,e_Hmag_1,Ksmag,e_Ksmag,Qual_1 if(ar__6(20:24) .EQ. '') Jmag_1 = rNULL__ if(ar__6(26:29) .EQ. '') e_Jmag_1 = rNULL__ if(ar__6(31:35) .EQ. '') Hmag_1 = rNULL__ if(ar__6(37:40) .EQ. '') e_Hmag_1 = rNULL__ if(ar__6(42:46) .EQ. '') Ksmag = rNULL__ if(ar__6(48:51) .EQ. '') e_Ksmag = rNULL__ c ..............Just test output........... write(6,7) + ID_2,Jmag_1,e_Jmag_1,Hmag_1,e_Hmag_1,Ksmag,e_Ksmag,Qual_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table11.dat' ! PMS stars not detected in X-rays (Spitzer) C Format for file interpretation 8 format( + A18,1X,A1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X, + F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A1) C Effective file loading open(unit=1,status='old',file= +'table11.dat') write(6,*) '....Loading file: table11.dat' do i__=1,69 read(1,'(A81)')ar__7 read(ar__7,8) + ID_3,f_ID,v_3_6__1,e__3_6__1,v_4_5__1,e__4_5__1,v_5_6_, + e__5_6_,v_8_0__1,e__8_0__1,v_24__1,e__24__1,AK_1,Cl if(ar__7(33:37) .EQ. '') v_4_5__1 = rNULL__ if(ar__7(39:42) .EQ. '') e__4_5__1 = rNULL__ if(ar__7(44:48) .EQ. '') v_5_6_ = rNULL__ if(ar__7(50:53) .EQ. '') e__5_6_ = rNULL__ if(ar__7(55:59) .EQ. '') v_8_0__1 = rNULL__ if(ar__7(61:64) .EQ. '') e__8_0__1 = rNULL__ if(ar__7(66:69) .EQ. '') v_24__1 = rNULL__ if(ar__7(71:74) .EQ. '') e__24__1 = rNULL__ c ..............Just test output........... write(6,8) + ID_3,f_ID,v_3_6__1,e__3_6__1,v_4_5__1,e__4_5__1,v_5_6_, + e__5_6_,v_8_0__1,e__8_0__1,v_24__1,e__24__1,AK_1,Cl c .......End.of.Just test output........... end do close(1) C============================================================================= stop end