Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/359/113 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-28 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/359/113 X-ray and IR study of Rho Oph dark cloud (Grosso+, 2000) *================================================================================ *X-rays and regions of star formation: *a combined ROSAT-HRI/near-to-mid IR study of the Rho Oph dark cloud * Grosso N., Montmerle T., Bontemps S., Andre P., Feigelson E.D. * <Astron. Astrophys. 359, 113 (2000)> * =2000A&A...359..113G C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablea1.dat' ! HRI X-ray sources integer*4 nr__ parameter (nr__=63) ! Number of records character*102 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) character*3 ROXR (nr__) ! ROXR designation integer*4 RAh (nr__) ! (h) Right ascension (J2000) (1) integer*4 RAm (nr__) ! (min) Right ascension (J2000) (1) real*4 RAs (nr__) ! (s) Right ascension (J2000) (1) character*1 DE_ (nr__) ! Declination sign (J2000) (1) integer*4 DEd (nr__) ! (deg) Declination (J2000) (1) integer*4 DEm (nr__) ! (arcmin) Declination (J2000) (1) integer*4 DEs (nr__) ! (arcsec) Declination (J2000) (1) integer*4 RA1950h (nr__) ! (h) Right ascension (B1950) (1) integer*4 RA1950m (nr__) ! (min) Right ascension (B1950) (1) real*4 RA1950s (nr__) ! (s) Right ascension (B1950) (1) character*1 DE1950_ (nr__) ! Declination sign (B1950) (1) integer*4 DE1950d (nr__) ! (deg) Declination (B1950) (1) integer*4 DE1950m (nr__) ! (arcmin) Declination (B1950) (1) integer*4 DE1950s (nr__) ! (arcsec) Declination (B1950) (1) integer*4 e_Pos (nr__) ! (arcsec) rms uncertainty on Position (1) real*4 LH (nr__) ! Likelihood of existence value (1) character*2 l_ct1 (nr__) ! [<= ] Limit flag on ct1 real*4 ct1 (nr__) ! (ct/ks) ? Count rate of observation 1 real*4 e_ct1 (nr__) ! (ct/ks) ? rms uncertainty on ct1 character*2 l_ct2 (nr__) ! [<= ] Limit flag on ct2 real*4 ct2 (nr__) ! (ct/ks) ? Count rate of observation 2 real*4 e_ct2 (nr__) ! (ct/ks) ? rms uncertainty on ct2 character*2 l_ct3 (nr__) ! [<= ] Limit flag on ct3 real*4 ct3 (nr__) ! (ct/ks) ? Count rate of observation 3 real*4 e_ct3 (nr__) ! (ct/ks) ? rms uncertainty on ct3 character*2 l_ct1_2_3 (nr__) ! [<= ] Limit flag on ct real*4 ct1_2_3 (nr__) ! (ct/ks) ? Count rate of observation (1+2+3) real*4 e_ct1_2_3 (nr__) ! (ct/ks) ? rms uncertainty on ct *Note (1): For the core F field, the indicated positions and the LH values * correspond to observation (1, 2 or 3) where LH and the position * accuracy are the best, i.e. when the count rate is highest. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tableb1.dat' ! HRI X-ray sources counterpart identifications in the {rho} Oph cloud cores A & F integer*4 nr__1 parameter (nr__1=63) ! Number of records character*102 ar__1 ! Full-size record character*3 ROXR_1 (nr__1) ! X-ray source number in this article * (ANN or FNN) integer*4 ROXR1 (nr__1) ! ? Casanova et al. (1995ApJ...439..752C) * X-ray source number (ROXR1) character*1 n_ROXR1 (nr__1) ! [* ou] Note (1) integer*4 ROXR2 (nr__1) ! ? Casanova (1994, Ph.D. thesis, Paris VI * University) X-ray source number (ROXR2) character*2 n_ROXR2 (nr__1) ! [* ou] Note (1) character*2 ASCA (nr__1) ! Kamata et al. (1997PASJ...49..461K) * X-ray source number character*1 n_ASCA (nr__1) ! [ou] Note (1) character*20 Name (nr__1) ! Counterpart name character*1 n_Name (nr__1) ! [* ou] Note (1) character*4 ISO (nr__1) ! IR excess from Bontemps et al. * (2000, A&A, submitted) (2) character*2 n_ISO (nr__1) ! [* ou] Note (1) character*7 ClassIR (nr__1) ! IR class (3) character*1 n_ClassIR (nr__1) ! [*] Note in notes.dat file real*4 AV (nr__1) ! (mag) ? Visual extinction character*1 n_AV (nr__1) ! [*] Note in notes.dat file real*8 Lum (nr__1) ! (solLum) ? Stellar luminosity character*1 n_Lum (nr__1) ! [*] Note in notes.dat file character*2 l_logLx1 (nr__1) ! [<= ] Limit flag on logLx1 real*4 logLx1 (nr__1) ! ([10-7W]) ? X-ray luminosity for d=140pc (obs. 1) character*1 u_logLx1 (nr__1) ! [?] Uncertainty flag on logLx1 character*2 l_logLx2 (nr__1) ! [<= ] Limit flag on logLx2 real*4 logLx2 (nr__1) ! ([10-7W]) ? X-ray luminosity for d=140pc (obs. 2) character*1 u_logLx2 (nr__1) ! [?] Uncertainty flag on logLx2 character*2 l_logLx3 (nr__1) ! [<= ] Limit flag on logLx3 real*4 logLx3 (nr__1) ! ([10-7W]) ? X-ray luminosity for d=140pc (obs. 3) character*1 u_logLx3 (nr__1) ! [?] Uncertainty flag on logLx3 character*2 l_logLxtot (nr__1) ! [<= ] Limit flag on logLxtot real*4 logLxtot (nr__1) ! ([10-7W]) ? X-ray luminosity for d=140pc (obs. 1+2+3) character*1 u_logLxtot (nr__1) ! [?] Uncertainty flag on logLxtot *Note (1): o = out of observation field * u = unobserved source * * = note in notes.dat file *Note (2): red = ISOCAM source with IR excess * blue = ISOCAM source without IR excess *Note (3): IR Class: * nII = new class II * nIII = new class III * ? = X-ray detected source for which intrinsic X-ray luminosity * cannot be determined * II-III? = class II or class III candidate c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'notes.dat' ! Individual notes integer*4 nr__2 parameter (nr__2=66) ! Number of records character*80 ar__2 ! Full-size record character*3 ROXR_2 (nr__2) ! ROXR number character*72 Note (nr__2) ! Text of the note (1) *Note (1): Acronyms used in the text: * B = Barsony et al. (1997, Cat. <J/ApJS/112/109>) * BBRCG = Barsony et al. (1989ApJ...346L..93B) * DoAr = Dolidze and Arakelyan (1959AZh....36..444D) * El = Elias (1978ApJ...224..453E) * GSS = Grasdalen et al. (1973ApJ...184L..53G) * GY = Greene and Young (1992ApJ...395..516G) * HD = The Henry Draper catalogue (Draper, 1918, Cat. <III/135>) * IRS = Wilking et al. (1989ApJ...340..823W) * ROXC = Montmerle et al. (1983ApJ...269..182M) * ROXs = Bouvier and Appenzeller (1992A&AS...92..481B) * S = abbreviation for ``Source'' in * Grasdalen et al. (1973ApJ...184L..53G) * SKS = Table 1 of Strom et al. (1995, Cat. <J/ApJ/438/813>) * SR = Struve and Rudkjoebing (1949ApJ...109...92S) * U = The PMM USNO-A1.0 Catalogue (Monet et al., 1996, Cat. <I/243>) * VSS = Vrba et al. (1976AJ.....81..958V) * VSSG = Vrba et al. (1975ApJ...197...77V) * WL = Wilking and Lada (1983ApJ...274..698W) * WSB = Wilking et al. (1987AJ.....94..106W) * YLW = Young et al. (1989ApJ...340..823W) C============================================================================= C Loading file 'tablea1.dat' ! HRI X-ray sources C Format for file interpretation 1 format( + A3,2X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,2X,I2,1X,I2,1X, + F4.1,1X,A1,I2,1X,I2,1X,I2,2X,I1,1X,F6.1,1X,A2,F4.1,2X,F3.1,1X, + A2,F4.1,1X,F3.1,1X,A2,F4.1,1X,F3.1,1X,A2,F4.1,1X,F3.1) C Effective file loading open(unit=1,status='old',file= +'tablea1.dat') write(6,*) '....Loading file: tablea1.dat' do i__=1,63 read(1,'(A102)')ar__ read(ar__,1) + ROXR(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),RA1950h(i__),RA1950m(i__),RA1950s(i__), + DE1950_(i__),DE1950d(i__),DE1950m(i__),DE1950s(i__), + e_Pos(i__),LH(i__),l_ct1(i__),ct1(i__),e_ct1(i__),l_ct2(i__), + ct2(i__),e_ct2(i__),l_ct3(i__),ct3(i__),e_ct3(i__), + l_ct1_2_3(i__),ct1_2_3(i__),e_ct1_2_3(i__) if(ar__(61:64) .EQ. '') ct1(i__) = rNULL__ if(ar__(67:69) .EQ. '') e_ct1(i__) = rNULL__ if(ar__(73:76) .EQ. '') ct2(i__) = rNULL__ if(ar__(78:80) .EQ. '') e_ct2(i__) = rNULL__ if(ar__(84:87) .EQ. '') ct3(i__) = rNULL__ if(ar__(89:91) .EQ. '') e_ct3(i__) = rNULL__ if(ar__(95:98) .EQ. '') ct1_2_3(i__) = rNULL__ if(ar__(100:102) .EQ. '') e_ct1_2_3(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) + ROXR(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),RA1950h(i__),RA1950m(i__),RA1950s(i__), + DE1950_(i__),DE1950d(i__),DE1950m(i__),DE1950s(i__), + e_Pos(i__),LH(i__),l_ct1(i__),ct1(i__),e_ct1(i__),l_ct2(i__), + ct2(i__),e_ct2(i__),l_ct3(i__),ct3(i__),e_ct3(i__), + l_ct1_2_3(i__),ct1_2_3(i__),e_ct1_2_3(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 'tableb1.dat' ! HRI X-ray sources counterpart identifications in * the {rho} Oph cloud cores A & F C Format for file interpretation 2 format( + A3,1X,I2,A1,1X,I2,A2,1X,A2,A1,1X,A20,A1,1X,A4,A2,1X,A7,A1,1X, + F4.1,A1,1X,F7.2,A1,1X,A2,F4.1,A1,1X,A2,F4.1,A1,1X,A2,F4.1,A1, + 1X,A2,1X,F4.1,A1) C Effective file loading open(unit=1,status='old',file= +'tableb1.dat') write(6,*) '....Loading file: tableb1.dat' do i__=1,63 read(1,'(A102)')ar__1 read(ar__1,2) + ROXR_1(i__),ROXR1(i__),n_ROXR1(i__),ROXR2(i__),n_ROXR2(i__), + ASCA(i__),n_ASCA(i__),Name(i__),n_Name(i__),ISO(i__), + n_ISO(i__),ClassIR(i__),n_ClassIR(i__),AV(i__),n_AV(i__), + Lum(i__),n_Lum(i__),l_logLx1(i__),logLx1(i__),u_logLx1(i__), + l_logLx2(i__),logLx2(i__),u_logLx2(i__),l_logLx3(i__), + logLx3(i__),u_logLx3(i__),l_logLxtot(i__),logLxtot(i__), + u_logLxtot(i__) if(ar__1(5:6) .EQ. '') ROXR1(i__) = iNULL__ if(ar__1(9:10) .EQ. '') ROXR2(i__) = iNULL__ if(ar__1(56:59) .EQ. '') AV(i__) = rNULL__ if(ar__1(62:68) .EQ. '') Lum(i__) = rNULL__ if(ar__1(73:76) .EQ. '') logLx1(i__) = rNULL__ if(ar__1(81:84) .EQ. '') logLx2(i__) = rNULL__ if(ar__1(89:92) .EQ. '') logLx3(i__) = rNULL__ if(ar__1(98:101) .EQ. '') logLxtot(i__) = rNULL__ c ..............Just test output........... write(6,2) + ROXR_1(i__),ROXR1(i__),n_ROXR1(i__),ROXR2(i__),n_ROXR2(i__), + ASCA(i__),n_ASCA(i__),Name(i__),n_Name(i__),ISO(i__), + n_ISO(i__),ClassIR(i__),n_ClassIR(i__),AV(i__),n_AV(i__), + Lum(i__),n_Lum(i__),l_logLx1(i__),logLx1(i__),u_logLx1(i__), + l_logLx2(i__),logLx2(i__),u_logLx2(i__),l_logLx3(i__), + logLx3(i__),u_logLx3(i__),l_logLxtot(i__),logLxtot(i__), + u_logLxtot(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'notes.dat' ! Individual notes C Format for file interpretation 3 format(4X,A3,1X,A72) C Effective file loading open(unit=1,status='old',file= +'notes.dat') write(6,*) '....Loading file: notes.dat' do i__=1,66 read(1,'(A80)')ar__2 read(ar__2,3)ROXR_2(i__),Note(i__) c ..............Just test output........... write(6,3)ROXR_2(i__),Note(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end