Conversion of standardized ReadMe file for
file /./ftp/cats/IX/30 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 2026-Oct-11
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. IX/30 Second ROSAT PSPC Catalog (ROSAT, 2000)
*================================================================================
*The Second ROSAT Source Catalog of Pointed Observations
* ROSAT Consortium
* <ROSAT News 72, 25-May-2000 (2000)>
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '2rxp.dat' ! The 2RXP Catalog
integer*4 nr__
parameter (nr__=95331) ! Number of records
character*393 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*16 Seq (nr__) ! ROSAT detection sequence number
character*1 S (nr__) ! *[S] for source in "short" catalog
character*1 m_2RXP (nr__) ! *[a-l] Multiplicity index
character*16 v2RXP (nr__) ! *2RXP designation (Jhhmmss.s+ddmmss)
integer*4 ErrPos (nr__) ! (arcsec) ? Error in position (G1)
real*4 Crate (nr__) ! (ct/s) Estimate of net count rate
real*4 e_Crate (nr__) ! (ct/s) ? Error of Crate (G1)
real*4 LhMdet (nr__) ! Source likelihood from the L/MDETEC algorithm
* (MD_LIKE, truncated at 999.9)
integer*4 oax (nr__) ! (arcmin) Off-axis angle averaged over observation
integer*4 ExpTime (nr__) ! (s) Exposure time
character*1 Var (nr__) ! [NTP] Source variability flag:
* Negative / Too_few / Possible
character*13 SrcFlags (nr__) ! *[cdefislmnpavw.] Source Flags
character*8 SeqFlags (nr__) ! *[DABEFMSX.] Sequence flags
integer*4 QCdate (nr__) ! ("YYMMDD") Date of sequence and source checking
character*4 QCsite (nr__) ! Checker Site
character*4 QDname (nr__) ! Checker ID
character*5 Com (nr__) ! Source comment
integer*4 Tbeg (nr__) ! ("YYMMDD") Observation start date
integer*4 Tend (nr__) ! ("YYMMDD") Observation finish date
character*15 ObsID (nr__) ! Observation ID (e.g. US400623H.N1)
character*4 Site (nr__) ! Processing Site (USRS / MPE / PE / SRS)
real*8 RAdeg (nr__) ! (deg) Source Right ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Source Declination (J2000)
real*8 GLON (nr__) ! (deg) Galactic longitude
real*8 GLAT (nr__) ! (deg) Galactic latitude
real*4 imX (nr__) ! (15arcsec) Blocked x source position
real*4 imY (nr__) ! (15arcsec) Blocked Y source position
real*4 e_imX (nr__) ! (15arcsec) ? Error on imX (G1)
real*4 e_imY (nr__) ! (15arcsec) ? Error on imY (G1)
integer*4 Drb (nr__) ! (arcmin) Corrected source-rib distance
* (truncated at 99)
integer*4 Dsr (nr__) ! (arcmin) [0,99]?=999 Corrected source-source distance
* (truncated at 99)
real*4 Vig (nr__) ! *? Vignetting factor (G1)
real*4 bgCt (nr__) ! (ct/s/arcmin2) Background countrate
* (truncated to -9.99 ... 99.99)
character*6 DetFlags (nr__) ! *[01] Source detection string (3)
character*1 DetML (nr__) ! *[ABCDb ] band in which ML detection occurred
integer*4 LhMdet2 (nr__) ! Source likelihood from the L/MDETEC algorithm
* (truncated at 9999)
integer*4 Dcs (nr__) ! (arcsec) Width of detect square
integer*4 MSX (nr__) ! Source number from merged source list (MPLSX)
character*1 Mdet (nr__) ! [M.] set to 'M' if Crate was taken from MDETEC
* (as well as some other parameters)
integer*4 LhML (nr__) ! ? Source likelihood from the ML algorithm
* (ML_LIKE, truncated to 999) (G1)
integer*4 rML (nr__) ! (arcsec) Extraction radius used in ML algorithm
integer*4 Ext (nr__) ! (arcsec) sigma of Gaussion distribution fitted to data
integer*4 ExtL (nr__) ! Likelihood of source extent
* (truncated to 9999)
character*4 Vflags (nr__) ! *[ckuf.] variability flags
real*4 HR1 (nr__) ! *Hardness ratio 1 (truncated to -9.99...99.99)
real*4 e_HR1 (nr__) ! Error on HR1
real*4 HR2 (nr__) ! *hardness ratio 2 (truncated to -9.99...99.99)
real*4 e_HR2 (nr__) ! Error on HR2
integer*4 CtA (nr__) ! (ct) Source net counts PHAch. 11-41 (~0.1-0.4 keV)
integer*4 CtB (nr__) ! (ct) Source net counts PHAch. 52-201 (~0.5-2.0 keV)
integer*4 CtC (nr__) ! (ct) Source net counts PHAch. 52-90 (~0.5-0.9 keV)
integer*4 CtD (nr__) ! (ct) Source net counts PHAch. 91-201 (~0.9-2.0 keV)
real*4 e_CtA (nr__) ! (ct) Mean error on CtA
real*4 e_CtB (nr__) ! (ct) Mean error on CtB
real*4 e_CtC (nr__) ! (ct) Mean error on CtC
real*4 e_CtD (nr__) ! (ct) Mean error on CtD
real*4 cSpl (nr__) ! Reduced chi-square of power law fit
* to net spectrum (truncated at 99.9)
real*4 cSRS (nr__) ! reduced chi-square of Raymond-Smith fit
* to net spectrum (truncated at 99.9)
*Note on S: this flag was added at CDS.
* The "short" version file, 2rxs_s.dat, can exactly be reconstructed
* from the "full" version file 2rxs.dat by a filter on this column,
* i.e. the file 2rxs_s.dat can be recreated by the Unix shell command
* grep '^....................S' 2rxs.dat > 2rxs_s.dat
*Note on 2RXP, m_2RXP: the name is NOT UNIQUE.
* Detections with the same position will have the same detection name;
* such detections may or may not represent the same astrophysical
* source. A letter (a,b,etc...) was added in the m_2RXP column for
* such duplicated duplicated names.
*Note on SrcFlags: flags with a dot (.) are turned off (false).
* c=CONFUSED : close sources confused detection algorithm
* d=DET_STRUCT : source near detector structure
* e=EXTENDED : source is probably resolved
* f=FALSE_DET : detection judged to be spurious
* i=INTEN_SUSP : source intensity suspect
* l=LTC_SUSP : source lightcurve may be suspect
* s=SPEC_SUSP : source spectrum suspect
* m=MULT_DET : complex emission split into individual sources
* n=NOT_CHECKED : source not checked visually by operator
* p=POS_SUSP : source position suspect
* a=ASSOCIATED : another source nearby may influence derived source
* v=VARFLG_SUSP : SASS variability flag suspect
* w=WITHIN_EXT : source within extended emission
*Note on SeqFlags: flags with a dot (.) are turned off (false).
* D=DEFERRED : Field not visually examined by operator [always off]
* A=ASP_SUSP : Aspect solution suspect [always off]
* B=BKG_SUSP : Background suspect
* E=EXT_EMIS : Extended emission in field
* F=FIELD_FDET : False detections in field
* M=SRC_MISS : Sources in field missed
* S=SLIS_MIS : Source list missing
* X=EXPO_MAP : Exposure map is not correct
*Note on Vig: the vignetting factor increases from 1 for sources detected
* on the axis to higher values when the source was detected off-axis.
*Note on DetFlags: the flags indicate the sliding window detection history
* using either the background map (M) or the local background (L),
* in the hard / soft and broad bands, and take the values
* 0 = no detection, 1 = detection
* The order of the flags is:
* M-broad, L-broad, M-hard, L-hard, M-soft, L-soft
*Note on DetML: the value indicates in which band the source was detected:
* A = PHA range 11-41 (~0.1-0.4 keV)
* B = PHA range 52-201 (~0.5-2.0 keV)
* C = PHA range 52-90 (~0.5-0.9 keV)
* D = PHA range 91-201 (~0.9-2.0 keV)
* b or blank means 'broad'
*Note on Vflags: flags with a dot (.) are turned off (false).
* c = variability flag from chi-square test (CSQ_FLAG)
* k = variability flag from KS test (KS_FLAG )
* u = variability flag from U test (U_FLAG )
* f = variability flag from FFT test (FFT_FLAG)
*Note on HR1, HR2: The hardness ratios are defined by:
* HR1 = ( B - A ) / ( B + A ) , where
* A = countrate in PHA range 11 - 41 (~0.1-0.4 keV)
* B = countrate in PHA range 52 - 201 (~0.5-2.0 keV)
* HR2 = ( D - C ) / ( D + C ) , where
* C = countrate in PHA range 52 - 90 (~0.5-0.9 keV)
* D = countrate in PHA range 91 - 201 (~0.9-2.0 keV)
* (PHA = Pulse Height Amplitude)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '2rxp_s.dat' ! The "short" 2RXP Catalog
(restricted to high confidence detections)
integer*4 nr__1
parameter (nr__1=43156) ! Number of records
character*393 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*16 Seq_1 (nr__1) ! ROSAT detection sequence number
character*1 S_1 (nr__1) ! *[S] for source in "short" catalog
character*1 m_2RXP_1 (nr__1) ! *[a-l] Multiplicity index
character*16 v2RXP_1 (nr__1) ! *2RXP designation (Jhhmmss.s+ddmmss)
integer*4 ErrPos_1 (nr__1) ! (arcsec) ? Error in position (G1)
real*4 Crate_1 (nr__1) ! (ct/s) Estimate of net count rate
real*4 e_Crate_1 (nr__1) ! (ct/s) ? Error of Crate (G1)
real*4 LhMdet_1 (nr__1) ! Source likelihood from the L/MDETEC algorithm
* (MD_LIKE, truncated at 999.9)
integer*4 oax_1 (nr__1) ! (arcmin) Off-axis angle averaged over observation
integer*4 ExpTime_1 (nr__1) ! (s) Exposure time
character*1 Var_1 (nr__1) ! [NTP] Source variability flag:
* Negative / Too_few / Possible
character*13 SrcFlags_1 (nr__1) ! *[cdefislmnpavw.] Source Flags
character*8 SeqFlags_1 (nr__1) ! *[DABEFMSX.] Sequence flags
integer*4 QCdate_1 (nr__1) ! ("YYMMDD") Date of sequence and source checking
character*4 QCsite_1 (nr__1) ! Checker Site
character*4 QDname_1 (nr__1) ! Checker ID
character*5 Com_1 (nr__1) ! Source comment
integer*4 Tbeg_1 (nr__1) ! ("YYMMDD") Observation start date
integer*4 Tend_1 (nr__1) ! ("YYMMDD") Observation finish date
character*15 ObsID_1 (nr__1) ! Observation ID (e.g. US400623H.N1)
character*4 Site_1 (nr__1) ! Processing Site (USRS / MPE / PE / SRS)
real*8 RAdeg_1 (nr__1) ! (deg) Source Right ascension (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Source Declination (J2000)
real*8 GLON_1 (nr__1) ! (deg) Galactic longitude
real*8 GLAT_1 (nr__1) ! (deg) Galactic latitude
real*4 imX_1 (nr__1) ! (15arcsec) Blocked x source position
real*4 imY_1 (nr__1) ! (15arcsec) Blocked Y source position
real*4 e_imX_1 (nr__1) ! (15arcsec) ? Error on imX (G1)
real*4 e_imY_1 (nr__1) ! (15arcsec) ? Error on imY (G1)
integer*4 Drb_1 (nr__1) ! (arcmin) Corrected source-rib distance
* (truncated at 99)
integer*4 Dsr_1 (nr__1) ! (arcmin) [0,99]?=999 Corrected source-source distance
* (truncated at 99)
real*4 Vig_1 (nr__1) ! *? Vignetting factor (G1)
real*4 bgCt_1 (nr__1) ! (ct/s/arcmin2) Background countrate
* (truncated to -9.99 ... 99.99)
character*6 DetFlags_1 (nr__1) ! *[01] Source detection string (3)
character*1 DetML_1 (nr__1) ! *[ABCDb ] band in which ML detection occurred
integer*4 LhMdet2_1 (nr__1) ! Source likelihood from the L/MDETEC algorithm
* (truncated at 9999)
integer*4 Dcs_1 (nr__1) ! (arcsec) Width of detect square
integer*4 MSX_1 (nr__1) ! Source number from merged source list (MPLSX)
character*1 Mdet_1 (nr__1) ! [M.] set to 'M' if Crate was taken from MDETEC
* (as well as some other parameters)
integer*4 LhML_1 (nr__1) ! ? Source likelihood from the ML algorithm
* (ML_LIKE, truncated to 999) (G1)
integer*4 rML_1 (nr__1) ! (arcsec) Extraction radius used in ML algorithm
integer*4 Ext_1 (nr__1) ! (arcsec) sigma of Gaussion distribution fitted to data
integer*4 ExtL_1 (nr__1) ! Likelihood of source extent
* (truncated to 9999)
character*4 Vflags_1 (nr__1) ! *[ckuf.] variability flags
real*4 HR1_1 (nr__1) ! *Hardness ratio 1 (truncated to -9.99...99.99)
real*4 e_HR1_1 (nr__1) ! Error on HR1
real*4 HR2_1 (nr__1) ! *hardness ratio 2 (truncated to -9.99...99.99)
real*4 e_HR2_1 (nr__1) ! Error on HR2
integer*4 CtA_1 (nr__1) ! (ct) Source net counts PHAch. 11-41 (~0.1-0.4 keV)
integer*4 CtB_1 (nr__1) ! (ct) Source net counts PHAch. 52-201 (~0.5-2.0 keV)
integer*4 CtC_1 (nr__1) ! (ct) Source net counts PHAch. 52-90 (~0.5-0.9 keV)
integer*4 CtD_1 (nr__1) ! (ct) Source net counts PHAch. 91-201 (~0.9-2.0 keV)
real*4 e_CtA_1 (nr__1) ! (ct) Mean error on CtA
real*4 e_CtB_1 (nr__1) ! (ct) Mean error on CtB
real*4 e_CtC_1 (nr__1) ! (ct) Mean error on CtC
real*4 e_CtD_1 (nr__1) ! (ct) Mean error on CtD
real*4 cSpl_1 (nr__1) ! Reduced chi-square of power law fit
* to net spectrum (truncated at 99.9)
real*4 cSRS_1 (nr__1) ! reduced chi-square of Raymond-Smith fit
* to net spectrum (truncated at 99.9)
*Note on S: this flag was added at CDS.
* The "short" version file, 2rxs_s.dat, can exactly be reconstructed
* from the "full" version file 2rxs.dat by a filter on this column,
* i.e. the file 2rxs_s.dat can be recreated by the Unix shell command
* grep '^....................S' 2rxs.dat > 2rxs_s.dat
*Note on 2RXP, m_2RXP: the name is NOT UNIQUE.
* Detections with the same position will have the same detection name;
* such detections may or may not represent the same astrophysical
* source. A letter (a,b,etc...) was added in the m_2RXP column for
* such duplicated duplicated names.
*Note on SrcFlags: flags with a dot (.) are turned off (false).
* c=CONFUSED : close sources confused detection algorithm
* d=DET_STRUCT : source near detector structure
* e=EXTENDED : source is probably resolved
* f=FALSE_DET : detection judged to be spurious
* i=INTEN_SUSP : source intensity suspect
* l=LTC_SUSP : source lightcurve may be suspect
* s=SPEC_SUSP : source spectrum suspect
* m=MULT_DET : complex emission split into individual sources
* n=NOT_CHECKED : source not checked visually by operator
* p=POS_SUSP : source position suspect
* a=ASSOCIATED : another source nearby may influence derived source
* v=VARFLG_SUSP : SASS variability flag suspect
* w=WITHIN_EXT : source within extended emission
*Note on SeqFlags: flags with a dot (.) are turned off (false).
* D=DEFERRED : Field not visually examined by operator [always off]
* A=ASP_SUSP : Aspect solution suspect [always off]
* B=BKG_SUSP : Background suspect
* E=EXT_EMIS : Extended emission in field
* F=FIELD_FDET : False detections in field
* M=SRC_MISS : Sources in field missed
* S=SLIS_MIS : Source list missing
* X=EXPO_MAP : Exposure map is not correct
*Note on Vig: the vignetting factor increases from 1 for sources detected
* on the axis to higher values when the source was detected off-axis.
*Note on DetFlags: the flags indicate the sliding window detection history
* using either the background map (M) or the local background (L),
* in the hard / soft and broad bands, and take the values
* 0 = no detection, 1 = detection
* The order of the flags is:
* M-broad, L-broad, M-hard, L-hard, M-soft, L-soft
*Note on DetML: the value indicates in which band the source was detected:
* A = PHA range 11-41 (~0.1-0.4 keV)
* B = PHA range 52-201 (~0.5-2.0 keV)
* C = PHA range 52-90 (~0.5-0.9 keV)
* D = PHA range 91-201 (~0.9-2.0 keV)
* b or blank means 'broad'
*Note on Vflags: flags with a dot (.) are turned off (false).
* c = variability flag from chi-square test (CSQ_FLAG)
* k = variability flag from KS test (KS_FLAG )
* u = variability flag from U test (U_FLAG )
* f = variability flag from FFT test (FFT_FLAG)
*Note on HR1, HR2: The hardness ratios are defined by:
* HR1 = ( B - A ) / ( B + A ) , where
* A = countrate in PHA range 11 - 41 (~0.1-0.4 keV)
* B = countrate in PHA range 52 - 201 (~0.5-2.0 keV)
* HR2 = ( D - C ) / ( D + C ) , where
* C = countrate in PHA range 52 - 90 (~0.5-0.9 keV)
* D = countrate in PHA range 91 - 201 (~0.9-2.0 keV)
* (PHA = Pulse Height Amplitude)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'seqp.dat' ! The sequences of observations
integer*4 nr__2
parameter (nr__2=4093) ! Number of records
character*174 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*24 RDF (nr__2) ! Result Data directory and File
character*4 QCSite_2 (nr__2) ! Checker Site
character*4 QDname_2 (nr__2) ! Checker ID
integer*4 Cdate (nr__2) ! ("YYMMDD") Checking date
character*8 SeqFlags_2 (nr__2) ! [DABEFMSX.] Sequence flags, as in 2rxp.dat
integer*4 ExpTime_2 (nr__2) ! (s) Exposure Time
integer*4 Tbeg_2 (nr__2) ! ("YYMMDD") Observation begin date
integer*4 Tend_2 (nr__2) ! ("YYMMDD") Observation end date
integer*4 Nsrc (nr__2) ! Number of sources
integer*4 Mis (nr__2) ! Number of missed sources
integer*4 Nc (nr__2) ! Statistics of confused sources
integer*4 Nd (nr__2) ! Statistics of sources near detector
integer*4 Ne (nr__2) ! Statistics of probably resolved sources
integer*4 Nf (nr__2) ! Statistics of spurious detections
integer*4 Ni (nr__2) ! Statistics of suspect source intensity
integer*4 Ns (nr__2) ! Statistics of suspect lightcurve sources
integer*4 Nl (nr__2) ! Statistics of suspect spectrum sources
integer*4 Nm (nr__2) ! Statistics of split complex emission sources
integer*4 Nn (nr__2) ! Statistics of sources not checked visually
integer*4 Np (nr__2) ! Statistics of suspect source positions
integer*4 Na (nr__2) ! Statistics of nearby associated sources
integer*4 Nv (nr__2) ! Statistics of suspected variability
integer*4 Nw (nr__2) ! Statistics of sources within extended emission
real*8 RAdeg_2 (nr__2) ! (deg) J2000 Right Ascension of observation
real*8 DEdeg_2 (nr__2) ! (deg) J2000 Declination of observation
character*4 Site_2 (nr__2) ! Processing Site
C=============================================================================
C Loading file '2rxp.dat' ! The 2RXP Catalog
C Format for file interpretation
1 format(
+ A16,4X,A1,A1,5X,A16,1X,I3,2X,E9.4,1X,E8.3,1X,F5.1,1X,I2,2X,I5,
+ 1X,A1,2X,A13,1X,A8,1X,I6,1X,A4,1X,A4,1X,A5,3X,I6,1X,I6,1X,A15,
+ 2X,A4,3X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,2X,F6.2,1X,F6.2,2X,F4.2,
+ 1X,F5.2,2X,I2,1X,I3,1X,F5.3,1X,F5.2,1X,A6,A1,1X,I4,1X,I3,2X,
+ I3,1X,A1,1X,I4,1X,I3,2X,I3,1X,I4,2X,A4,2X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F4.2,1X,I6,1X,I6,1X,I6,1X,I6,1X,F6.1,1X,F6.1,1X,F6.1,
+ 1X,F6.1,1X,F4.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'2rxp.dat')
write(6,*) '....Loading file: 2rxp.dat'
do i__=1,95331
read(1,'(A393)')ar__
read(ar__,1)
+ Seq(i__),S(i__),m_2RXP(i__),v2RXP(i__),ErrPos(i__),Crate(i__),
+ e_Crate(i__),LhMdet(i__),oax(i__),ExpTime(i__),Var(i__),
+ SrcFlags(i__),SeqFlags(i__),QCdate(i__),QCsite(i__),
+ QDname(i__),Com(i__),Tbeg(i__),Tend(i__),ObsID(i__),Site(i__),
+ RAdeg(i__),DEdeg(i__),GLON(i__),GLAT(i__),imX(i__),imY(i__),
+ e_imX(i__),e_imY(i__),Drb(i__),Dsr(i__),Vig(i__),bgCt(i__),
+ DetFlags(i__),DetML(i__),LhMdet2(i__),Dcs(i__),MSX(i__),
+ Mdet(i__),LhML(i__),rML(i__),Ext(i__),ExtL(i__),Vflags(i__),
+ HR1(i__),e_HR1(i__),HR2(i__),e_HR2(i__),CtA(i__),CtB(i__),
+ CtC(i__),CtD(i__),e_CtA(i__),e_CtB(i__),e_CtC(i__),e_CtD(i__),
+ cSpl(i__),cSRS(i__)
if(ar__(45:47) .EQ. '') ErrPos(i__) = iNULL__
if(ar__(60:67) .EQ. '') e_Crate(i__) = rNULL__
if(ar__(226:229) .EQ. '') e_imX(i__) = rNULL__
if(ar__(231:235) .EQ. '') e_imY(i__) = rNULL__
if(ar__(245:249) .EQ. '') Vig(i__) = rNULL__
if(ar__(281:284) .EQ. '') LhML(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Seq(i__),S(i__),m_2RXP(i__),v2RXP(i__),ErrPos(i__),Crate(i__),
+ e_Crate(i__),LhMdet(i__),oax(i__),ExpTime(i__),Var(i__),
+ SrcFlags(i__),SeqFlags(i__),QCdate(i__),QCsite(i__),
+ QDname(i__),Com(i__),Tbeg(i__),Tend(i__),ObsID(i__),Site(i__),
+ RAdeg(i__),DEdeg(i__),GLON(i__),GLAT(i__),imX(i__),imY(i__),
+ e_imX(i__),e_imY(i__),Drb(i__),Dsr(i__),Vig(i__),bgCt(i__),
+ DetFlags(i__),DetML(i__),LhMdet2(i__),Dcs(i__),MSX(i__),
+ Mdet(i__),LhML(i__),rML(i__),Ext(i__),ExtL(i__),Vflags(i__),
+ HR1(i__),e_HR1(i__),HR2(i__),e_HR2(i__),CtA(i__),CtB(i__),
+ CtC(i__),CtD(i__),e_CtA(i__),e_CtB(i__),e_CtC(i__),e_CtD(i__),
+ cSpl(i__),cSRS(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file '2rxp_s.dat' ! The "short" 2RXP Catalog
* (restricted to high confidence detections)
C Format for file interpretation
2 format(
+ A16,4X,A1,A1,5X,A16,1X,I3,2X,E9.4,1X,E8.3,1X,F5.1,1X,I2,2X,I5,
+ 1X,A1,2X,A13,1X,A8,1X,I6,1X,A4,1X,A4,1X,A5,3X,I6,1X,I6,1X,A15,
+ 2X,A4,3X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,2X,F6.2,1X,F6.2,2X,F4.2,
+ 1X,F5.2,2X,I2,1X,I3,1X,F5.3,1X,F5.2,1X,A6,A1,1X,I4,1X,I3,2X,
+ I3,1X,A1,1X,I4,1X,I3,2X,I3,1X,I4,2X,A4,2X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F4.2,1X,I6,1X,I6,1X,I6,1X,I6,1X,F6.1,1X,F6.1,1X,F6.1,
+ 1X,F6.1,1X,F4.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'2rxp_s.dat')
write(6,*) '....Loading file: 2rxp_s.dat'
do i__=1,43156
read(1,'(A393)')ar__1
read(ar__1,2)
+ Seq_1(i__),S_1(i__),m_2RXP_1(i__),v2RXP_1(i__),ErrPos_1(i__),
+ Crate_1(i__),e_Crate_1(i__),LhMdet_1(i__),oax_1(i__),
+ ExpTime_1(i__),Var_1(i__),SrcFlags_1(i__),SeqFlags_1(i__),
+ QCdate_1(i__),QCsite_1(i__),QDname_1(i__),Com_1(i__),
+ Tbeg_1(i__),Tend_1(i__),ObsID_1(i__),Site_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),GLON_1(i__),GLAT_1(i__),imX_1(i__),imY_1(i__),
+ e_imX_1(i__),e_imY_1(i__),Drb_1(i__),Dsr_1(i__),Vig_1(i__),
+ bgCt_1(i__),DetFlags_1(i__),DetML_1(i__),LhMdet2_1(i__),
+ Dcs_1(i__),MSX_1(i__),Mdet_1(i__),LhML_1(i__),rML_1(i__),
+ Ext_1(i__),ExtL_1(i__),Vflags_1(i__),HR1_1(i__),e_HR1_1(i__),
+ HR2_1(i__),e_HR2_1(i__),CtA_1(i__),CtB_1(i__),CtC_1(i__),
+ CtD_1(i__),e_CtA_1(i__),e_CtB_1(i__),e_CtC_1(i__),
+ e_CtD_1(i__),cSpl_1(i__),cSRS_1(i__)
if(ar__1(45:47) .EQ. '') ErrPos_1(i__) = iNULL__
if(ar__1(60:67) .EQ. '') e_Crate_1(i__) = rNULL__
if(ar__1(226:229) .EQ. '') e_imX_1(i__) = rNULL__
if(ar__1(231:235) .EQ. '') e_imY_1(i__) = rNULL__
if(ar__1(245:249) .EQ. '') Vig_1(i__) = rNULL__
if(ar__1(281:284) .EQ. '') LhML_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),S_1(i__),m_2RXP_1(i__),v2RXP_1(i__),ErrPos_1(i__),
+ Crate_1(i__),e_Crate_1(i__),LhMdet_1(i__),oax_1(i__),
+ ExpTime_1(i__),Var_1(i__),SrcFlags_1(i__),SeqFlags_1(i__),
+ QCdate_1(i__),QCsite_1(i__),QDname_1(i__),Com_1(i__),
+ Tbeg_1(i__),Tend_1(i__),ObsID_1(i__),Site_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),GLON_1(i__),GLAT_1(i__),imX_1(i__),imY_1(i__),
+ e_imX_1(i__),e_imY_1(i__),Drb_1(i__),Dsr_1(i__),Vig_1(i__),
+ bgCt_1(i__),DetFlags_1(i__),DetML_1(i__),LhMdet2_1(i__),
+ Dcs_1(i__),MSX_1(i__),Mdet_1(i__),LhML_1(i__),rML_1(i__),
+ Ext_1(i__),ExtL_1(i__),Vflags_1(i__),HR1_1(i__),e_HR1_1(i__),
+ HR2_1(i__),e_HR2_1(i__),CtA_1(i__),CtB_1(i__),CtC_1(i__),
+ CtD_1(i__),e_CtA_1(i__),e_CtB_1(i__),e_CtC_1(i__),
+ e_CtD_1(i__),cSpl_1(i__),cSRS_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'seqp.dat' ! The sequences of observations
C Format for file interpretation
3 format(
+ A24,1X,A4,1X,A4,1X,I6,10X,A8,3X,I5,2X,I6,1X,I6,2X,I3,2X,I2,
+ 23X,I1,1X,I2,1X,I2,I3,I3,I3,I3,I3,I3,I3,I3,I3,I3,1X,F8.4,1X,
+ F8.4,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'seqp.dat')
write(6,*) '....Loading file: seqp.dat'
do i__=1,4093
read(1,'(A174)')ar__2
read(ar__2,3)
+ RDF(i__),QCSite_2(i__),QDname_2(i__),Cdate(i__),
+ SeqFlags_2(i__),ExpTime_2(i__),Tbeg_2(i__),Tend_2(i__),
+ Nsrc(i__),Mis(i__),Nc(i__),Nd(i__),Ne(i__),Nf(i__),Ni(i__),
+ Ns(i__),Nl(i__),Nm(i__),Nn(i__),Np(i__),Na(i__),Nv(i__),
+ Nw(i__),RAdeg_2(i__),DEdeg_2(i__),Site_2(i__)
c ..............Just test output...........
write(6,3)
+ RDF(i__),QCSite_2(i__),QDname_2(i__),Cdate(i__),
+ SeqFlags_2(i__),ExpTime_2(i__),Tbeg_2(i__),Tend_2(i__),
+ Nsrc(i__),Mis(i__),Nc(i__),Nd(i__),Ne(i__),Nf(i__),Ni(i__),
+ Ns(i__),Nl(i__),Nm(i__),Nn(i__),Np(i__),Na(i__),Nv(i__),
+ Nw(i__),RAdeg_2(i__),DEdeg_2(i__),Site_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end