Conversion of standardized ReadMe file for
file /./ftp/cats/IX/30 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 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 ! ROSAT detection sequence number
character*1 S ! *[S] for source in "short" catalog
character*1 m_2RXP ! *[a-l] Multiplicity index
character*16 v2RXP ! *2RXP designation (Jhhmmss.s+ddmmss)
integer*4 ErrPos ! (arcsec) ? Error in position (G1)
real*4 Crate ! (ct/s) Estimate of net count rate
real*4 e_Crate ! (ct/s) ? Error of Crate (G1)
real*4 LhMdet ! Source likelihood from the L/MDETEC algorithm
* (MD_LIKE, truncated at 999.9)
integer*4 oax ! (arcmin) Off-axis angle averaged over observation
integer*4 ExpTime ! (s) Exposure time
character*1 Var ! [NTP] Source variability flag:
* Negative / Too_few / Possible
character*13 SrcFlags ! *[cdefislmnpavw.] Source Flags
character*8 SeqFlags ! *[DABEFMSX.] Sequence flags
integer*4 QCdate ! ("YYMMDD") Date of sequence and source checking
character*4 QCsite ! Checker Site
character*4 QDname ! Checker ID
character*5 Com ! Source comment
integer*4 Tbeg ! ("YYMMDD") Observation start date
integer*4 Tend ! ("YYMMDD") Observation finish date
character*15 ObsID ! Observation ID (e.g. US400623H.N1)
character*4 Site ! Processing Site (USRS / MPE / PE / SRS)
real*8 RAdeg ! (deg) Source Right ascension (J2000)
real*8 DEdeg ! (deg) Source Declination (J2000)
real*8 GLON ! (deg) Galactic longitude
real*8 GLAT ! (deg) Galactic latitude
real*4 imX ! (15arcsec) Blocked x source position
real*4 imY ! (15arcsec) Blocked Y source position
real*4 e_imX ! (15arcsec) ? Error on imX (G1)
real*4 e_imY ! (15arcsec) ? Error on imY (G1)
integer*4 Drb ! (arcmin) Corrected source-rib distance
* (truncated at 99)
integer*4 Dsr ! (arcmin) [0,99]?=999 Corrected source-source distance
* (truncated at 99)
real*4 Vig ! *? Vignetting factor (G1)
real*4 bgCt ! (ct/s/arcmin2) Background countrate
* (truncated to -9.99 ... 99.99)
character*6 DetFlags ! *[01] Source detection string (3)
character*1 DetML ! *[ABCDb ] band in which ML detection occurred
integer*4 LhMdet2 ! Source likelihood from the L/MDETEC algorithm
* (truncated at 9999)
integer*4 Dcs ! (arcsec) Width of detect square
integer*4 MSX ! Source number from merged source list (MPLSX)
character*1 Mdet ! [M.] set to 'M' if Crate was taken from MDETEC
* (as well as some other parameters)
integer*4 LhML ! ? Source likelihood from the ML algorithm
* (ML_LIKE, truncated to 999) (G1)
integer*4 rML ! (arcsec) Extraction radius used in ML algorithm
integer*4 Ext ! (arcsec) sigma of Gaussion distribution fitted to data
integer*4 ExtL ! Likelihood of source extent
* (truncated to 9999)
character*4 Vflags ! *[ckuf.] variability flags
real*4 HR1 ! *Hardness ratio 1 (truncated to -9.99...99.99)
real*4 e_HR1 ! Error on HR1
real*4 HR2 ! *hardness ratio 2 (truncated to -9.99...99.99)
real*4 e_HR2 ! Error on HR2
integer*4 CtA ! (ct) Source net counts PHAch. 11-41 (~0.1-0.4 keV)
integer*4 CtB ! (ct) Source net counts PHAch. 52-201 (~0.5-2.0 keV)
integer*4 CtC ! (ct) Source net counts PHAch. 52-90 (~0.5-0.9 keV)
integer*4 CtD ! (ct) Source net counts PHAch. 91-201 (~0.9-2.0 keV)
real*4 e_CtA ! (ct) Mean error on CtA
real*4 e_CtB ! (ct) Mean error on CtB
real*4 e_CtC ! (ct) Mean error on CtC
real*4 e_CtD ! (ct) Mean error on CtD
real*4 cSpl ! Reduced chi-square of power law fit
* to net spectrum (truncated at 99.9)
real*4 cSRS ! 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 ! ROSAT detection sequence number
character*1 S_1 ! *[S] for source in "short" catalog
character*1 m_2RXP_1 ! *[a-l] Multiplicity index
character*16 v2RXP_1 ! *2RXP designation (Jhhmmss.s+ddmmss)
integer*4 ErrPos_1 ! (arcsec) ? Error in position (G1)
real*4 Crate_1 ! (ct/s) Estimate of net count rate
real*4 e_Crate_1 ! (ct/s) ? Error of Crate (G1)
real*4 LhMdet_1 ! Source likelihood from the L/MDETEC algorithm
* (MD_LIKE, truncated at 999.9)
integer*4 oax_1 ! (arcmin) Off-axis angle averaged over observation
integer*4 ExpTime_1 ! (s) Exposure time
character*1 Var_1 ! [NTP] Source variability flag:
* Negative / Too_few / Possible
character*13 SrcFlags_1 ! *[cdefislmnpavw.] Source Flags
character*8 SeqFlags_1 ! *[DABEFMSX.] Sequence flags
integer*4 QCdate_1 ! ("YYMMDD") Date of sequence and source checking
character*4 QCsite_1 ! Checker Site
character*4 QDname_1 ! Checker ID
character*5 Com_1 ! Source comment
integer*4 Tbeg_1 ! ("YYMMDD") Observation start date
integer*4 Tend_1 ! ("YYMMDD") Observation finish date
character*15 ObsID_1 ! Observation ID (e.g. US400623H.N1)
character*4 Site_1 ! Processing Site (USRS / MPE / PE / SRS)
real*8 RAdeg_1 ! (deg) Source Right ascension (J2000)
real*8 DEdeg_1 ! (deg) Source Declination (J2000)
real*8 GLON_1 ! (deg) Galactic longitude
real*8 GLAT_1 ! (deg) Galactic latitude
real*4 imX_1 ! (15arcsec) Blocked x source position
real*4 imY_1 ! (15arcsec) Blocked Y source position
real*4 e_imX_1 ! (15arcsec) ? Error on imX (G1)
real*4 e_imY_1 ! (15arcsec) ? Error on imY (G1)
integer*4 Drb_1 ! (arcmin) Corrected source-rib distance
* (truncated at 99)
integer*4 Dsr_1 ! (arcmin) [0,99]?=999 Corrected source-source distance
* (truncated at 99)
real*4 Vig_1 ! *? Vignetting factor (G1)
real*4 bgCt_1 ! (ct/s/arcmin2) Background countrate
* (truncated to -9.99 ... 99.99)
character*6 DetFlags_1 ! *[01] Source detection string (3)
character*1 DetML_1 ! *[ABCDb ] band in which ML detection occurred
integer*4 LhMdet2_1 ! Source likelihood from the L/MDETEC algorithm
* (truncated at 9999)
integer*4 Dcs_1 ! (arcsec) Width of detect square
integer*4 MSX_1 ! Source number from merged source list (MPLSX)
character*1 Mdet_1 ! [M.] set to 'M' if Crate was taken from MDETEC
* (as well as some other parameters)
integer*4 LhML_1 ! ? Source likelihood from the ML algorithm
* (ML_LIKE, truncated to 999) (G1)
integer*4 rML_1 ! (arcsec) Extraction radius used in ML algorithm
integer*4 Ext_1 ! (arcsec) sigma of Gaussion distribution fitted to data
integer*4 ExtL_1 ! Likelihood of source extent
* (truncated to 9999)
character*4 Vflags_1 ! *[ckuf.] variability flags
real*4 HR1_1 ! *Hardness ratio 1 (truncated to -9.99...99.99)
real*4 e_HR1_1 ! Error on HR1
real*4 HR2_1 ! *hardness ratio 2 (truncated to -9.99...99.99)
real*4 e_HR2_1 ! Error on HR2
integer*4 CtA_1 ! (ct) Source net counts PHAch. 11-41 (~0.1-0.4 keV)
integer*4 CtB_1 ! (ct) Source net counts PHAch. 52-201 (~0.5-2.0 keV)
integer*4 CtC_1 ! (ct) Source net counts PHAch. 52-90 (~0.5-0.9 keV)
integer*4 CtD_1 ! (ct) Source net counts PHAch. 91-201 (~0.9-2.0 keV)
real*4 e_CtA_1 ! (ct) Mean error on CtA
real*4 e_CtB_1 ! (ct) Mean error on CtB
real*4 e_CtC_1 ! (ct) Mean error on CtC
real*4 e_CtD_1 ! (ct) Mean error on CtD
real*4 cSpl_1 ! Reduced chi-square of power law fit
* to net spectrum (truncated at 99.9)
real*4 cSRS_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 ! Result Data directory and File
character*4 QCSite_2 ! Checker Site
character*4 QDname_2 ! Checker ID
integer*4 Cdate ! ("YYMMDD") Checking date
character*8 SeqFlags_2 ! [DABEFMSX.] Sequence flags, as in 2rxp.dat
integer*4 ExpTime_2 ! (s) Exposure Time
integer*4 Tbeg_2 ! ("YYMMDD") Observation begin date
integer*4 Tend_2 ! ("YYMMDD") Observation end date
integer*4 Nsrc ! Number of sources
integer*4 Mis ! Number of missed sources
integer*4 Nc ! Statistics of confused sources
integer*4 Nd ! Statistics of sources near detector
integer*4 Ne ! Statistics of probably resolved sources
integer*4 Nf ! Statistics of spurious detections
integer*4 Ni ! Statistics of suspect source intensity
integer*4 Ns ! Statistics of suspect lightcurve sources
integer*4 Nl ! Statistics of suspect spectrum sources
integer*4 Nm ! Statistics of split complex emission sources
integer*4 Nn ! Statistics of sources not checked visually
integer*4 Np ! Statistics of suspect source positions
integer*4 Na ! Statistics of nearby associated sources
integer*4 Nv ! Statistics of suspected variability
integer*4 Nw ! Statistics of sources within extended emission
real*8 RAdeg_2 ! (deg) J2000 Right Ascension of observation
real*8 DEdeg_2 ! (deg) J2000 Declination of observation
character*4 Site_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,S,m_2RXP,v2RXP,ErrPos,Crate,e_Crate,LhMdet,oax,ExpTime,
+ Var,SrcFlags,SeqFlags,QCdate,QCsite,QDname,Com,Tbeg,Tend,
+ ObsID,Site,RAdeg,DEdeg,GLON,GLAT,imX,imY,e_imX,e_imY,Drb,Dsr,
+ Vig,bgCt,DetFlags,DetML,LhMdet2,Dcs,MSX,Mdet,LhML,rML,Ext,
+ ExtL,Vflags,HR1,e_HR1,HR2,e_HR2,CtA,CtB,CtC,CtD,e_CtA,e_CtB,
+ e_CtC,e_CtD,cSpl,cSRS
if(ar__(45:47) .EQ. '') ErrPos = iNULL__
if(ar__(60:67) .EQ. '') e_Crate = rNULL__
if(ar__(226:229) .EQ. '') e_imX = rNULL__
if(ar__(231:235) .EQ. '') e_imY = rNULL__
if(ar__(245:249) .EQ. '') Vig = rNULL__
if(ar__(281:284) .EQ. '') LhML = iNULL__
c ..............Just test output...........
write(6,1)
+ Seq,S,m_2RXP,v2RXP,ErrPos,Crate,e_Crate,LhMdet,oax,ExpTime,
+ Var,SrcFlags,SeqFlags,QCdate,QCsite,QDname,Com,Tbeg,Tend,
+ ObsID,Site,RAdeg,DEdeg,GLON,GLAT,imX,imY,e_imX,e_imY,Drb,Dsr,
+ Vig,bgCt,DetFlags,DetML,LhMdet2,Dcs,MSX,Mdet,LhML,rML,Ext,
+ ExtL,Vflags,HR1,e_HR1,HR2,e_HR2,CtA,CtB,CtC,CtD,e_CtA,e_CtB,
+ e_CtC,e_CtD,cSpl,cSRS
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,S_1,m_2RXP_1,v2RXP_1,ErrPos_1,Crate_1,e_Crate_1,
+ LhMdet_1,oax_1,ExpTime_1,Var_1,SrcFlags_1,SeqFlags_1,QCdate_1,
+ QCsite_1,QDname_1,Com_1,Tbeg_1,Tend_1,ObsID_1,Site_1,RAdeg_1,
+ DEdeg_1,GLON_1,GLAT_1,imX_1,imY_1,e_imX_1,e_imY_1,Drb_1,Dsr_1,
+ Vig_1,bgCt_1,DetFlags_1,DetML_1,LhMdet2_1,Dcs_1,MSX_1,Mdet_1,
+ LhML_1,rML_1,Ext_1,ExtL_1,Vflags_1,HR1_1,e_HR1_1,HR2_1,
+ e_HR2_1,CtA_1,CtB_1,CtC_1,CtD_1,e_CtA_1,e_CtB_1,e_CtC_1,
+ e_CtD_1,cSpl_1,cSRS_1
if(ar__1(45:47) .EQ. '') ErrPos_1 = iNULL__
if(ar__1(60:67) .EQ. '') e_Crate_1 = rNULL__
if(ar__1(226:229) .EQ. '') e_imX_1 = rNULL__
if(ar__1(231:235) .EQ. '') e_imY_1 = rNULL__
if(ar__1(245:249) .EQ. '') Vig_1 = rNULL__
if(ar__1(281:284) .EQ. '') LhML_1 = iNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1,S_1,m_2RXP_1,v2RXP_1,ErrPos_1,Crate_1,e_Crate_1,
+ LhMdet_1,oax_1,ExpTime_1,Var_1,SrcFlags_1,SeqFlags_1,QCdate_1,
+ QCsite_1,QDname_1,Com_1,Tbeg_1,Tend_1,ObsID_1,Site_1,RAdeg_1,
+ DEdeg_1,GLON_1,GLAT_1,imX_1,imY_1,e_imX_1,e_imY_1,Drb_1,Dsr_1,
+ Vig_1,bgCt_1,DetFlags_1,DetML_1,LhMdet2_1,Dcs_1,MSX_1,Mdet_1,
+ LhML_1,rML_1,Ext_1,ExtL_1,Vflags_1,HR1_1,e_HR1_1,HR2_1,
+ e_HR2_1,CtA_1,CtB_1,CtC_1,CtD_1,e_CtA_1,e_CtB_1,e_CtC_1,
+ e_CtD_1,cSpl_1,cSRS_1
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,QCSite_2,QDname_2,Cdate,SeqFlags_2,ExpTime_2,Tbeg_2,
+ Tend_2,Nsrc,Mis,Nc,Nd,Ne,Nf,Ni,Ns,Nl,Nm,Nn,Np,Na,Nv,Nw,
+ RAdeg_2,DEdeg_2,Site_2
c ..............Just test output...........
write(6,3)
+ RDF,QCSite_2,QDname_2,Cdate,SeqFlags_2,ExpTime_2,Tbeg_2,
+ Tend_2,Nsrc,Mis,Nc,Nd,Ne,Nf,Ni,Ns,Nl,Nm,Nn,Np,Na,Nv,Nw,
+ RAdeg_2,DEdeg_2,Site_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end