Conversion of standardized ReadMe file for
file /./home/cats/J/A_A/480/611 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-May-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/A+A/480/611 XMM-Newton Slew Survey catalogue, XMMSL1 (Saxton+, 2008)
*================================================================================
*The first XMM-Newton Slew Survey catalogue: XMMSL1.
* Saxton R.D., Read A.M., Esquej P., Freyberg M.J., Altieri B., Bermejo D.
* <Astron. Astrophys., 480, 611-622 (2008)>
* =2008A&A...480..611S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'xmmsl1c.dat' ! XMMSL1 catalogue, clean catalog (V1.1, 07/04/08)
integer*4 nr__
parameter (nr__=4326) ! Number of records
character*1244 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*16 XMMSL1 (nr__) ! Name of the source (JHHMMSS.s+DDMMSS)
character*10 ObsID (nr__) ! Observation identification number
integer*4 SN (nr__) ! [1,139] Unique identification number
* for each detection in an observation
real*8 RAdeg (nr__) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
real*8 e_pos (nr__) ! (arcsec) rms uncertainty on position
real*8 GLON (nr__) ! (deg) Galactic longitude
real*8 GLAT (nr__) ! (deg) Galactic latitude
real*4 HR1 (nr__) ! ?=- Hardness ratio HR1=(H-S)/(H+S)
real*4 e_HR1 (nr__) ! ?=- rms uncertainty on HR1
character*19 Date_1 (nr__) ! ("date") Observation date start (YYYY-MM-DDThh:mm:ss)
character*19 Date_2 (nr__) ! ("date") Observation date, end (YYYY-MM-DDThh:mm:ss)
real*4 SCTST (nr__) ! (ct) ?=- Number of background subtracted
* counts, in the total energy band
* (0.2-12keV) (1)
real*4 e_SCTST (nr__) ! (ct) ?=- rms uncertainty on SCTST
real*4 SCTSH (nr__) ! (ct) ?=- Number of background subtracted counts
* in the hard energy band (2-12keV) (1)
real*4 e_SCTSH (nr__) ! (ct) ?=- rms uncertainty on SCTSH
real*4 SCTSS (nr__) ! (ct) ?=- Number of background subtracted counts
* in the soft energy band (0.2-2keV) (1)
real*4 e_SCTSS (nr__) ! (ct) ?=- rms uncertainty on SCTSS
real*8 XposT (nr__) ! (pix) ?=- X position on total band image (2)
real*8 YposT (nr__) ! (pix) ?=- Y position on total band image (2)
real*8 XposH (nr__) ! (pix) ?=- X position on hard band image (2)
real*8 YposH (nr__) ! (pix) ?=- Y position on hard band image (2)
real*8 XposS (nr__) ! (pix) ?=- X position on soft band image (2)
real*8 YposS (nr__) ! (pix) ?=- Y position on soft band image (2)
real*4 ExtT (nr__) ! (pix) ?=- Spatial extension of the source in
* the total energy band
real*4 e_ExtT (nr__) ! (pix) ?=- rms uncertainty on ExtT
real*4 ExtH (nr__) ! (pix) ?=- Spatial extension of the source in
* the hard energy band
real*4 e_ExtH (nr__) ! (pix) ?=- rms uncertainty on ExtH
real*4 ExtS (nr__) ! (pix) ?=- Spatial extension of the source in
* the soft energy band
real*4 e_ExtS (nr__) ! (pix) ?=- rms uncertainty on ExtS
real*4 DetMLT (nr__) ! ?=- Detection likelihood in the total band
real*4 DetMLH (nr__) ! ?=- Detection likelihood in the hard band
real*4 DetMLS (nr__) ! ?=- Detection likelihood in the soft band
real*4 ExtMLT (nr__) ! ?=- Likelihood of the source being
* extended in the total band
real*4 ExtMLH (nr__) ! ?=- Likelihood of the source being
* extended in the hard band
real*4 ExTMLS (nr__) ! ?=- Likelihood of the source being
* extended in the soft band
real*4 bgMapT (nr__) ! (ct/pix) ?=- Total band background value
real*4 bgMapH (nr__) ! (ct/pix) ?=- Hard band background value
real*4 bgMapS (nr__) ! (ct/pix) ?=- Soft band background value
real*4 ExpMapT (nr__) ! (s) ?=- Total band effective on-axis expo.time
real*4 ExpMapH (nr__) ! (s) ?=- Hard band effective on-axis expo.time
real*4 ExpMapS (nr__) ! (s) ?=- Soft band effective on-axis expo.time
real*4 FluxT (nr__) ! (fW/m2) ?=- Total band source flux
real*4 e_FluxT (nr__) ! (fW/m2) ?=- rms uncertainty on FluxT
real*4 FluxH (nr__) ! (fW/m2) ?=- Hard band source flux
real*4 e_FluxH (nr__) ! (fW/m2) ?=- rms uncertainty on FluxH
real*4 FluxS (nr__) ! (fW/m2) ?=- Soft band source flux
real*4 e_FluxS (nr__) ! (fW/m2) ?=- rms uncertainty on FluxS
real*4 RateT (nr__) ! (ct/s) ?=- Total band count rate
real*4 e_RateT (nr__) ! (ct/s) ?=- rms uncertainty on RateT
real*4 RateH (nr__) ! (ct/s) ?=- Hard band count rate
real*4 e_RateH (nr__) ! (ct/s) ?=- rms uncertainty on RateH
real*4 RateS (nr__) ! (ct/s) ?=- Soft band count rate
real*4 e_RateS (nr__) ! (ct/s) ?=- rms uncertainty on RateS
real*8 RAdegT (nr__) ! (deg) ?=- Total band right ascension (J2000)
real*8 DEdegT (nr__) ! (deg) ?=- Total band declination (J2000)
real*8 e_posT (nr__) ! (deg) ?=- Total band positional error
real*8 GLONT (nr__) ! (deg) ?=- Total band galactic longitude
real*8 GLATT (nr__) ! (deg) ?=- Total band galactic latitude
real*8 RAdegH (nr__) ! (deg) ?=- Hard band right ascension (J2000)
real*8 DEdegH (nr__) ! (deg) ?=- Hard band declination (J2000)
real*8 e_posH (nr__) ! (deg) ?=- Hard band positional error
real*8 GLONH (nr__) ! (deg) ?=- Hard band galactic longitude
real*8 GLATH (nr__) ! (deg) ?=- Hard band galactic latitude
real*8 RAdegS (nr__) ! (deg) ?=- Soft band right ascension (J2000)
real*8 DEdegS (nr__) ! (deg) ?=- Soft band declination (J2000)
real*8 e_posS (nr__) ! (deg) ?=- Soft band positional error
real*8 GLONS (nr__) ! (deg) ?=- Soft band galactic longitude
real*8 GLATS (nr__) ! (deg) ?=- Soft band galactic latitude
character*12 VarFlag (nr__) ! [CLEAN_SAMPLE XXXXXXXXXXXX]
character*37 SRCNameT (nr__) ! Source name defined in total band list (3)
character*37 SRCNameH (nr__) ! Source name defined in hard band list (3)
character*37 SRCNameS (nr__) ! Source name defined in soft band list (3)
character*40 XImNameT (nr__) ! Name of the image containing the total
* band detection
character*40 XImNameH (nr__) ! Name of the image containing the hard
* band detection
character*40 XImNameS (nr__) ! Name of the image containing the soft
* band detection
character*1 VerINEXT (nr__) ! [TF] Detection flag (4)
character*1 VerHALO (nr__) ! [TF] Detection flag (4)
character*1 VerHIBGND (nr__) ! [TF] Detection flag (4)
character*1 VerNREDG (nr__) ! [TF] Detection flag (4)
character*1 VerPSUSP (nr__) ! [TF] Detection flag (4)
character*1 VerFALSE (nr__) ! [TF] Detection flag (4)
character*40 Comment (nr__) ! Comment which explains why a particular
* source flag has been set to true
character*31 Name (nr__) ! Cross-correlations of the positions of
* the slew sources with astronomical
* databases and catalogues (best match)
character*30 OName (nr__) ! Alternative name for the best match source
character*21 RASS (nr__) ! Name of the closest Rosat All Sky Survey
* (RASS) source (5)
character*14 Type (nr__) ! Source type (6)
character*9 r_Name (nr__) ! Source of best match selection (7)
real*4 Dist (nr__) ! (arcmin) ? Distance between the best match
* candidate and the slew survey source
real*4 DistRASS (nr__) ! (arcsec) ?=- Distance from the best match ROSAT
* source and the slew source
character*3 ModeID (nr__) ! Observing mode of EPIC-pn camera (8)
real*4 bgRate (nr__) ! (ct/s) Peak background count rate in the image
* from which the source was extracted (9)
*Note (1): corrected for the PSF.
*Note (2): Units: 4.1"x4.1" pixels
*Note (3): The format is 'xs' followed by the revolution number, the observation
* ID and the source position, e.g. xs0841_9084100002_12:57:07.5+01:50:42
*Note (4): Detection flags signification as:
* VerINEXT: If set true, this flag notes that a detection has been
* found within an extended source and is probably false.
* This flag is used to ensure that only one source is quoted
* for the large supernova remnants.
* VerHALO: If set true, this flag indicates that a detection lies
* within the point spread function, or halo, of a very
* bright source and is probably spurious.
* VerHIBGND: If set true, this flag indicates that the source lies
* within a bright region, caused by high background, and
* is probably spurious.
* VerNREDG: If set true, this flag indicates that the source has been
* detected in two consecutive images of the same slew. This
* can occur if the source lies on the border between slews
* and implies that the position of the source is unlikely to
* be accurate.
* VerPSUSP: This flag indicates that the quoted position of the source
* is likely to be inaccurate. In a few slews the attitude
* reconstruction is poor and in these cases the actual
* position of the source is difficult to determine and can
* be wrong by an arcminute.
* VerFALSE: This is a catch-all for problems not included in the other
* flags. In the catalogue 29 sources have this flag set
* true, in two cases it is because the exposure time has
* been calculated as zero and in the other cases the
* detections are due to reflection from a close, but out of
* the field-of-view SCO X-1.
*Note (5): RASS name, 1RXS JHHMMSS.s+DDMMSS or JHHMMSS.s+DDMMSS
*Note (6): Source type as returned by SIMBAD, NED and the other resources
* used in the cross-matching process. This is directly taken from the
* catalogue in question and no attempt has been made to rationalise the
* values.
*Note (7): The astronomical database or catalogue from which the best match
* has been selected. e.g. SIMBAD, NED, etc.
*Note (8): Observing mode of the EPIC-pn camera ad follows:
* FF = Full frame mode
* eFF = extended full frame mode
* LW = large window mode
*Note (9): Peak background count rate in the image from which the source was
* extracted. This is measured as the count rate, over the whole image,
* for events with energy >10keV (PI>10000). It is quoted to a resolution
* of 0.1ct/s and is used to determine the overall background environment
* in which a source was detected. This value is used in the creation of
* a clean subsample of the catalogue (see top level description).
* Generally, the higher this value, the greater the probability that the
* source is spurious. This is especially true if the detection
* likelihood of the source is less than 14. The value has been tabulated
* for images with a background rate greater than 2ct/s. For quieter
* images the value has been artifically set to 1.9ct/s.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'xmmsl1t.dat' ! XMMSL1 catalogue, total catalog (V1.1, 07/04/08)
integer*4 nr__1
parameter (nr__1=8662) ! Number of records
character*1244 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*16 XMMSL1_1 (nr__1) ! Name of the source (JHHMMSS.s+DDMMSS)
character*10 ObsID_1 (nr__1) ! Observation identification number
integer*4 SN_1 (nr__1) ! [1,139] Unique identification number
* for each detection in an observation
real*8 RAdeg_1 (nr__1) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Declination in decimal degrees (J2000)
real*8 e_pos_1 (nr__1) ! (arcsec) rms uncertainty on position
real*8 GLON_1 (nr__1) ! (deg) Galactic longitude
real*8 GLAT_1 (nr__1) ! (deg) Galactic latitude
real*4 HR1_1 (nr__1) ! ?=- Hardness ratio HR1=(H-S)/(H+S)
real*4 e_HR1_1 (nr__1) ! ?=- rms uncertainty on HR1
character*19 Date_1_1 (nr__1) ! ("date") Observation date start (YYYY-MM-DDThh:mm:ss)
character*19 Date_2_1 (nr__1) ! ("date") Observation date, end (YYYY-MM-DDThh:mm:ss)
real*4 SCTST_1 (nr__1) ! (ct) ?=- Number of background subtracted
* counts, in the total energy band
* (0.2-12keV) (1)
real*4 e_SCTST_1 (nr__1) ! (ct) ?=- rms uncertainty on SCTST
real*4 SCTSH_1 (nr__1) ! (ct) ?=- Number of background subtracted counts
* in the hard energy band (2-12keV) (1)
real*4 e_SCTSH_1 (nr__1) ! (ct) ?=- rms uncertainty on SCTSH
real*4 SCTSS_1 (nr__1) ! (ct) ?=- Number of background subtracted counts
* in the soft energy band (0.2-2keV) (1)
real*4 e_SCTSS_1 (nr__1) ! (ct) ?=- rms uncertainty on SCTSS
real*8 XposT_1 (nr__1) ! (pix) ?=- X position on total band image (2)
real*8 YposT_1 (nr__1) ! (pix) ?=- Y position on total band image (2)
real*8 XposH_1 (nr__1) ! (pix) ?=- X position on hard band image (2)
real*8 YposH_1 (nr__1) ! (pix) ?=- Y position on hard band image (2)
real*8 XposS_1 (nr__1) ! (pix) ?=- X position on soft band image (2)
real*8 YposS_1 (nr__1) ! (pix) ?=- Y position on soft band image (2)
real*4 ExtT_1 (nr__1) ! (pix) ?=- Spatial extension of the source in
* the total energy band
real*4 e_ExtT_1 (nr__1) ! (pix) ?=- rms uncertainty on ExtT
real*4 ExtH_1 (nr__1) ! (pix) ?=- Spatial extension of the source in
* the hard energy band
real*4 e_ExtH_1 (nr__1) ! (pix) ?=- rms uncertainty on ExtH
real*4 ExtS_1 (nr__1) ! (pix) ?=- Spatial extension of the source in
* the soft energy band
real*4 e_ExtS_1 (nr__1) ! (pix) ?=- rms uncertainty on ExtS
real*4 DetMLT_1 (nr__1) ! ?=- Detection likelihood in the total band
real*4 DetMLH_1 (nr__1) ! ?=- Detection likelihood in the hard band
real*4 DetMLS_1 (nr__1) ! ?=- Detection likelihood in the soft band
real*4 ExtMLT_1 (nr__1) ! ?=- Likelihood of the source being
* extended in the total band
real*4 ExtMLH_1 (nr__1) ! ?=- Likelihood of the source being
* extended in the hard band
real*4 ExTMLS_1 (nr__1) ! ?=- Likelihood of the source being
* extended in the soft band
real*4 bgMapT_1 (nr__1) ! (ct/pix) ?=- Total band background value
real*4 bgMapH_1 (nr__1) ! (ct/pix) ?=- Hard band background value
real*4 bgMapS_1 (nr__1) ! (ct/pix) ?=- Soft band background value
real*4 ExpMapT_1 (nr__1) ! (s) ?=- Total band effective on-axis expo.time
real*4 ExpMapH_1 (nr__1) ! (s) ?=- Hard band effective on-axis expo.time
real*4 ExpMapS_1 (nr__1) ! (s) ?=- Soft band effective on-axis expo.time
real*4 FluxT_1 (nr__1) ! (fW/m2) ?=- Total band source flux
real*4 e_FluxT_1 (nr__1) ! (fW/m2) ?=- rms uncertainty on FluxT
real*4 FluxH_1 (nr__1) ! (fW/m2) ?=- Hard band source flux
real*4 e_FluxH_1 (nr__1) ! (fW/m2) ?=- rms uncertainty on FluxH
real*4 FluxS_1 (nr__1) ! (fW/m2) ?=- Soft band source flux
real*4 e_FluxS_1 (nr__1) ! (fW/m2) ?=- rms uncertainty on FluxS
real*4 RateT_1 (nr__1) ! (ct/s) ?=- Total band count rate
real*4 e_RateT_1 (nr__1) ! (ct/s) ?=- rms uncertainty on RateT
real*4 RateH_1 (nr__1) ! (ct/s) ?=- Hard band count rate
real*4 e_RateH_1 (nr__1) ! (ct/s) ?=- rms uncertainty on RateH
real*4 RateS_1 (nr__1) ! (ct/s) ?=- Soft band count rate
real*4 e_RateS_1 (nr__1) ! (ct/s) ?=- rms uncertainty on RateS
real*8 RAdegT_1 (nr__1) ! (deg) ?=- Total band right ascension (J2000)
real*8 DEdegT_1 (nr__1) ! (deg) ?=- Total band declination (J2000)
real*8 e_posT_1 (nr__1) ! (deg) ?=- Total band positional error
real*8 GLONT_1 (nr__1) ! (deg) ?=- Total band galactic longitude
real*8 GLATT_1 (nr__1) ! (deg) ?=- Total band galactic latitude
real*8 RAdegH_1 (nr__1) ! (deg) ?=- Hard band right ascension (J2000)
real*8 DEdegH_1 (nr__1) ! (deg) ?=- Hard band declination (J2000)
real*8 e_posH_1 (nr__1) ! (deg) ?=- Hard band positional error
real*8 GLONH_1 (nr__1) ! (deg) ?=- Hard band galactic longitude
real*8 GLATH_1 (nr__1) ! (deg) ?=- Hard band galactic latitude
real*8 RAdegS_1 (nr__1) ! (deg) ?=- Soft band right ascension (J2000)
real*8 DEdegS_1 (nr__1) ! (deg) ?=- Soft band declination (J2000)
real*8 e_posS_1 (nr__1) ! (deg) ?=- Soft band positional error
real*8 GLONS_1 (nr__1) ! (deg) ?=- Soft band galactic longitude
real*8 GLATS_1 (nr__1) ! (deg) ?=- Soft band galactic latitude
character*12 VarFlag_1 (nr__1) ! [CLEAN_SAMPLE XXXXXXXXXXXX]
character*37 SRCNameT_1 (nr__1) ! Source name defined in total band list (3)
character*37 SRCNameH_1 (nr__1) ! Source name defined in hard band list (3)
character*37 SRCNameS_1 (nr__1) ! Source name defined in soft band list (3)
character*40 XImNameT_1 (nr__1) ! Name of the image containing the total
* band detection
character*40 XImNameH_1 (nr__1) ! Name of the image containing the hard
* band detection
character*40 XImNameS_1 (nr__1) ! Name of the image containing the soft
* band detection
character*1 VerINEXT_1 (nr__1) ! [TF] Detection flag (4)
character*1 VerHALO_1 (nr__1) ! [TF] Detection flag (4)
character*1 VerHIBGND_1(nr__1) ! [TF] Detection flag (4)
character*1 VerNREDG_1 (nr__1) ! [TF] Detection flag (4)
character*1 VerPSUSP_1 (nr__1) ! [TF] Detection flag (4)
character*1 VerFALSE_1 (nr__1) ! [TF] Detection flag (4)
character*40 Comment_1 (nr__1) ! Comment which explains why a particular
* source flag has been set to true
character*31 Name_1 (nr__1) ! Cross-correlations of the positions of
* the slew sources with astronomical
* databases and catalogues (best match)
character*30 OName_1 (nr__1) ! Alternative name for the best match source
character*21 RASS_1 (nr__1) ! Name of the closest Rosat All Sky Survey
* (RASS) source (5)
character*14 Type_1 (nr__1) ! Source type (6)
character*9 r_Name_1 (nr__1) ! Source of best match selection (7)
real*4 Dist_1 (nr__1) ! (arcmin) ? Distance between the best match
* candidate and the slew survey source
real*4 DistRASS_1 (nr__1) ! (arcsec) ?=- Distance from the best match ROSAT
* source and the slew source
character*3 ModeID_1 (nr__1) ! Observing mode of EPIC-pn camera (8)
real*4 bgRate_1 (nr__1) ! (ct/s) Peak background count rate in the image
* from which the source was extracted (9)
*Note (1): corrected for the PSF.
*Note (2): Units: 4.1"x4.1" pixels
*Note (3): The format is 'xs' followed by the revolution number, the observation
* ID and the source position, e.g. xs0841_9084100002_12:57:07.5+01:50:42
*Note (4): Detection flags signification as:
* VerINEXT: If set true, this flag notes that a detection has been
* found within an extended source and is probably false.
* This flag is used to ensure that only one source is quoted
* for the large supernova remnants.
* VerHALO: If set true, this flag indicates that a detection lies
* within the point spread function, or halo, of a very
* bright source and is probably spurious.
* VerHIBGND: If set true, this flag indicates that the source lies
* within a bright region, caused by high background, and
* is probably spurious.
* VerNREDG: If set true, this flag indicates that the source has been
* detected in two consecutive images of the same slew. This
* can occur if the source lies on the border between slews
* and implies that the position of the source is unlikely to
* be accurate.
* VerPSUSP: This flag indicates that the quoted position of the source
* is likely to be inaccurate. In a few slews the attitude
* reconstruction is poor and in these cases the actual
* position of the source is difficult to determine and can
* be wrong by an arcminute.
* VerFALSE: This is a catch-all for problems not included in the other
* flags. In the catalogue 29 sources have this flag set
* true, in two cases it is because the exposure time has
* been calculated as zero and in the other cases the
* detections are due to reflection from a close, but out of
* the field-of-view SCO X-1.
*Note (5): RASS name, 1RXS JHHMMSS.s+DDMMSS or JHHMMSS.s+DDMMSS
*Note (6): Source type as returned by SIMBAD, NED and the other resources
* used in the cross-matching process. This is directly taken from the
* catalogue in question and no attempt has been made to rationalise the
* values.
*Note (7): The astronomical database or catalogue from which the best match
* has been selected. e.g. SIMBAD, NED, etc.
*Note (8): Observing mode of the EPIC-pn camera ad follows:
* FF = Full frame mode
* eFF = extended full frame mode
* LW = large window mode
*Note (9): Peak background count rate in the image from which the source was
* extracted. This is measured as the count rate, over the whole image,
* for events with energy >10keV (PI>10000). It is quoted to a resolution
* of 0.1ct/s and is used to determine the overall background environment
* in which a source was detected. This value is used in the creation of
* a clean subsample of the catalogue (see top level description).
* Generally, the higher this value, the greater the probability that the
* source is spurious. This is especially true if the detection
* likelihood of the source is less than 14. The value has been tabulated
* for images with a background rate greater than 2ct/s. For quieter
* images the value has been artifically set to 1.9ct/s.
C=============================================================================
C Loading file 'xmmsl1c.dat' ! XMMSL1 catalogue, clean catalog (V1.1, 07/04/08)
C Format for file interpretation
1 format(
+ 7X,A16,1X,A10,1X,I3,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F6.3,1X,F6.3,1X,A19,1X,A19,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,
+ F7.2,1X,F7.2,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,A12,1X,A37,1X,
+ A37,1X,A37,1X,A40,1X,A40,1X,A40,1X,A1,1X,A1,1X,A1,1X,A1,1X,A1,
+ 1X,A1,1X,A40,1X,A31,1X,A30,1X,A21,1X,A14,1X,A9,1X,E12.6,1X,
+ E12.6,1X,A3,1X,E12.7)
C Effective file loading
open(unit=1,status='old',file=
+'xmmsl1c.dat')
write(6,*) '....Loading file: xmmsl1c.dat'
do i__=1,4326
read(1,'(A1244)')ar__
read(ar__,1)
+ XMMSL1(i__),ObsID(i__),SN(i__),RAdeg(i__),DEdeg(i__),
+ e_pos(i__),GLON(i__),GLAT(i__),HR1(i__),e_HR1(i__),
+ Date_1(i__),Date_2(i__),SCTST(i__),e_SCTST(i__),SCTSH(i__),
+ e_SCTSH(i__),SCTSS(i__),e_SCTSS(i__),XposT(i__),YposT(i__),
+ XposH(i__),YposH(i__),XposS(i__),YposS(i__),ExtT(i__),
+ e_ExtT(i__),ExtH(i__),e_ExtH(i__),ExtS(i__),e_ExtS(i__),
+ DetMLT(i__),DetMLH(i__),DetMLS(i__),ExtMLT(i__),ExtMLH(i__),
+ ExTMLS(i__),bgMapT(i__),bgMapH(i__),bgMapS(i__),ExpMapT(i__),
+ ExpMapH(i__),ExpMapS(i__),FluxT(i__),e_FluxT(i__),FluxH(i__),
+ e_FluxH(i__),FluxS(i__),e_FluxS(i__),RateT(i__),e_RateT(i__),
+ RateH(i__),e_RateH(i__),RateS(i__),e_RateS(i__),RAdegT(i__),
+ DEdegT(i__),e_posT(i__),GLONT(i__),GLATT(i__),RAdegH(i__),
+ DEdegH(i__),e_posH(i__),GLONH(i__),GLATH(i__),RAdegS(i__),
+ DEdegS(i__),e_posS(i__),GLONS(i__),GLATS(i__),VarFlag(i__),
+ SRCNameT(i__),SRCNameH(i__),SRCNameS(i__),XImNameT(i__),
+ XImNameH(i__),XImNameS(i__),VerINEXT(i__),VerHALO(i__),
+ VerHIBGND(i__),VerNREDG(i__),VerPSUSP(i__),VerFALSE(i__),
+ Comment(i__),Name(i__),OName(i__),RASS(i__),Type(i__),
+ r_Name(i__),Dist(i__),DistRASS(i__),ModeID(i__),bgRate(i__)
if (idig(ar__(85:90)).EQ.0) HR1(i__) = rNULL__
if (idig(ar__(92:97)).EQ.0) e_HR1(i__) = rNULL__
if (idig(ar__(139:150)).EQ.0) SCTST(i__) = rNULL__
if (idig(ar__(152:163)).EQ.0) e_SCTST(i__) = rNULL__
if (idig(ar__(165:176)).EQ.0) SCTSH(i__) = rNULL__
if (idig(ar__(178:189)).EQ.0) e_SCTSH(i__) = rNULL__
if (idig(ar__(191:202)).EQ.0) SCTSS(i__) = rNULL__
if (idig(ar__(204:215)).EQ.0) e_SCTSS(i__) = rNULL__
if (idig(ar__(217:223)).EQ.0) XposT(i__) = rNULL__
if (idig(ar__(225:231)).EQ.0) YposT(i__) = rNULL__
if (idig(ar__(233:239)).EQ.0) XposH(i__) = rNULL__
if (idig(ar__(241:247)).EQ.0) YposH(i__) = rNULL__
if (idig(ar__(249:255)).EQ.0) XposS(i__) = rNULL__
if (idig(ar__(257:263)).EQ.0) YposS(i__) = rNULL__
if (idig(ar__(265:276)).EQ.0) ExtT(i__) = rNULL__
if (idig(ar__(278:289)).EQ.0) e_ExtT(i__) = rNULL__
if (idig(ar__(291:302)).EQ.0) ExtH(i__) = rNULL__
if (idig(ar__(304:315)).EQ.0) e_ExtH(i__) = rNULL__
if (idig(ar__(317:328)).EQ.0) ExtS(i__) = rNULL__
if (idig(ar__(330:341)).EQ.0) e_ExtS(i__) = rNULL__
if (idig(ar__(343:354)).EQ.0) DetMLT(i__) = rNULL__
if (idig(ar__(356:367)).EQ.0) DetMLH(i__) = rNULL__
if (idig(ar__(369:380)).EQ.0) DetMLS(i__) = rNULL__
if (idig(ar__(382:393)).EQ.0) ExtMLT(i__) = rNULL__
if (idig(ar__(395:406)).EQ.0) ExtMLH(i__) = rNULL__
if (idig(ar__(408:419)).EQ.0) ExTMLS(i__) = rNULL__
if (idig(ar__(421:432)).EQ.0) bgMapT(i__) = rNULL__
if (idig(ar__(434:445)).EQ.0) bgMapH(i__) = rNULL__
if (idig(ar__(447:458)).EQ.0) bgMapS(i__) = rNULL__
if (idig(ar__(460:471)).EQ.0) ExpMapT(i__) = rNULL__
if (idig(ar__(473:484)).EQ.0) ExpMapH(i__) = rNULL__
if (idig(ar__(486:497)).EQ.0) ExpMapS(i__) = rNULL__
if (idig(ar__(499:510)).EQ.0) FluxT(i__) = rNULL__
if (idig(ar__(512:523)).EQ.0) e_FluxT(i__) = rNULL__
if (idig(ar__(525:536)).EQ.0) FluxH(i__) = rNULL__
if (idig(ar__(538:549)).EQ.0) e_FluxH(i__) = rNULL__
if (idig(ar__(551:562)).EQ.0) FluxS(i__) = rNULL__
if (idig(ar__(564:575)).EQ.0) e_FluxS(i__) = rNULL__
if (idig(ar__(577:588)).EQ.0) RateT(i__) = rNULL__
if (idig(ar__(590:601)).EQ.0) e_RateT(i__) = rNULL__
if (idig(ar__(603:614)).EQ.0) RateH(i__) = rNULL__
if (idig(ar__(616:627)).EQ.0) e_RateH(i__) = rNULL__
if (idig(ar__(629:640)).EQ.0) RateS(i__) = rNULL__
if (idig(ar__(642:653)).EQ.0) e_RateS(i__) = rNULL__
if (idig(ar__(655:662)).EQ.0) RAdegT(i__) = rNULL__
if (idig(ar__(664:671)).EQ.0) DEdegT(i__) = rNULL__
if (idig(ar__(673:680)).EQ.0) e_posT(i__) = rNULL__
if (idig(ar__(682:689)).EQ.0) GLONT(i__) = rNULL__
if (idig(ar__(691:698)).EQ.0) GLATT(i__) = rNULL__
if (idig(ar__(700:707)).EQ.0) RAdegH(i__) = rNULL__
if (idig(ar__(709:716)).EQ.0) DEdegH(i__) = rNULL__
if (idig(ar__(718:725)).EQ.0) e_posH(i__) = rNULL__
if (idig(ar__(727:734)).EQ.0) GLONH(i__) = rNULL__
if (idig(ar__(736:743)).EQ.0) GLATH(i__) = rNULL__
if (idig(ar__(745:752)).EQ.0) RAdegS(i__) = rNULL__
if (idig(ar__(754:761)).EQ.0) DEdegS(i__) = rNULL__
if (idig(ar__(763:770)).EQ.0) e_posS(i__) = rNULL__
if (idig(ar__(772:779)).EQ.0) GLONS(i__) = rNULL__
if (idig(ar__(781:788)).EQ.0) GLATS(i__) = rNULL__
if(ar__(1203:1214) .EQ. '') Dist(i__) = rNULL__
if (idig(ar__(1216:1227)).EQ.0) DistRASS(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ XMMSL1(i__),ObsID(i__),SN(i__),RAdeg(i__),DEdeg(i__),
+ e_pos(i__),GLON(i__),GLAT(i__),HR1(i__),e_HR1(i__),
+ Date_1(i__),Date_2(i__),SCTST(i__),e_SCTST(i__),SCTSH(i__),
+ e_SCTSH(i__),SCTSS(i__),e_SCTSS(i__),XposT(i__),YposT(i__),
+ XposH(i__),YposH(i__),XposS(i__),YposS(i__),ExtT(i__),
+ e_ExtT(i__),ExtH(i__),e_ExtH(i__),ExtS(i__),e_ExtS(i__),
+ DetMLT(i__),DetMLH(i__),DetMLS(i__),ExtMLT(i__),ExtMLH(i__),
+ ExTMLS(i__),bgMapT(i__),bgMapH(i__),bgMapS(i__),ExpMapT(i__),
+ ExpMapH(i__),ExpMapS(i__),FluxT(i__),e_FluxT(i__),FluxH(i__),
+ e_FluxH(i__),FluxS(i__),e_FluxS(i__),RateT(i__),e_RateT(i__),
+ RateH(i__),e_RateH(i__),RateS(i__),e_RateS(i__),RAdegT(i__),
+ DEdegT(i__),e_posT(i__),GLONT(i__),GLATT(i__),RAdegH(i__),
+ DEdegH(i__),e_posH(i__),GLONH(i__),GLATH(i__),RAdegS(i__),
+ DEdegS(i__),e_posS(i__),GLONS(i__),GLATS(i__),VarFlag(i__),
+ SRCNameT(i__),SRCNameH(i__),SRCNameS(i__),XImNameT(i__),
+ XImNameH(i__),XImNameS(i__),VerINEXT(i__),VerHALO(i__),
+ VerHIBGND(i__),VerNREDG(i__),VerPSUSP(i__),VerFALSE(i__),
+ Comment(i__),Name(i__),OName(i__),RASS(i__),Type(i__),
+ r_Name(i__),Dist(i__),DistRASS(i__),ModeID(i__),bgRate(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'xmmsl1t.dat' ! XMMSL1 catalogue, total catalog (V1.1, 07/04/08)
C Format for file interpretation
2 format(
+ 7X,A16,1X,A10,1X,I3,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F6.3,1X,F6.3,1X,A19,1X,A19,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,
+ F7.2,1X,F7.2,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,A12,1X,A37,1X,
+ A37,1X,A37,1X,A40,1X,A40,1X,A40,1X,A1,1X,A1,1X,A1,1X,A1,1X,A1,
+ 1X,A1,1X,A40,1X,A31,1X,A30,1X,A21,1X,A14,1X,A9,1X,E12.6,1X,
+ E12.6,1X,A3,1X,E12.7)
C Effective file loading
open(unit=1,status='old',file=
+'xmmsl1t.dat')
write(6,*) '....Loading file: xmmsl1t.dat'
do i__=1,8662
read(1,'(A1244)')ar__1
read(ar__1,2)
+ XMMSL1_1(i__),ObsID_1(i__),SN_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),e_pos_1(i__),GLON_1(i__),GLAT_1(i__),HR1_1(i__),
+ e_HR1_1(i__),Date_1_1(i__),Date_2_1(i__),SCTST_1(i__),
+ e_SCTST_1(i__),SCTSH_1(i__),e_SCTSH_1(i__),SCTSS_1(i__),
+ e_SCTSS_1(i__),XposT_1(i__),YposT_1(i__),XposH_1(i__),
+ YposH_1(i__),XposS_1(i__),YposS_1(i__),ExtT_1(i__),
+ e_ExtT_1(i__),ExtH_1(i__),e_ExtH_1(i__),ExtS_1(i__),
+ e_ExtS_1(i__),DetMLT_1(i__),DetMLH_1(i__),DetMLS_1(i__),
+ ExtMLT_1(i__),ExtMLH_1(i__),ExTMLS_1(i__),bgMapT_1(i__),
+ bgMapH_1(i__),bgMapS_1(i__),ExpMapT_1(i__),ExpMapH_1(i__),
+ ExpMapS_1(i__),FluxT_1(i__),e_FluxT_1(i__),FluxH_1(i__),
+ e_FluxH_1(i__),FluxS_1(i__),e_FluxS_1(i__),RateT_1(i__),
+ e_RateT_1(i__),RateH_1(i__),e_RateH_1(i__),RateS_1(i__),
+ e_RateS_1(i__),RAdegT_1(i__),DEdegT_1(i__),e_posT_1(i__),
+ GLONT_1(i__),GLATT_1(i__),RAdegH_1(i__),DEdegH_1(i__),
+ e_posH_1(i__),GLONH_1(i__),GLATH_1(i__),RAdegS_1(i__),
+ DEdegS_1(i__),e_posS_1(i__),GLONS_1(i__),GLATS_1(i__),
+ VarFlag_1(i__),SRCNameT_1(i__),SRCNameH_1(i__),
+ SRCNameS_1(i__),XImNameT_1(i__),XImNameH_1(i__),
+ XImNameS_1(i__),VerINEXT_1(i__),VerHALO_1(i__),
+ VerHIBGND_1(i__),VerNREDG_1(i__),VerPSUSP_1(i__),
+ VerFALSE_1(i__),Comment_1(i__),Name_1(i__),OName_1(i__),
+ RASS_1(i__),Type_1(i__),r_Name_1(i__),Dist_1(i__),
+ DistRASS_1(i__),ModeID_1(i__),bgRate_1(i__)
if (idig(ar__1(85:90)).EQ.0) HR1_1(i__) = rNULL__
if (idig(ar__1(92:97)).EQ.0) e_HR1_1(i__) = rNULL__
if (idig(ar__1(139:150)).EQ.0) SCTST_1(i__) = rNULL__
if (idig(ar__1(152:163)).EQ.0) e_SCTST_1(i__) = rNULL__
if (idig(ar__1(165:176)).EQ.0) SCTSH_1(i__) = rNULL__
if (idig(ar__1(178:189)).EQ.0) e_SCTSH_1(i__) = rNULL__
if (idig(ar__1(191:202)).EQ.0) SCTSS_1(i__) = rNULL__
if (idig(ar__1(204:215)).EQ.0) e_SCTSS_1(i__) = rNULL__
if (idig(ar__1(217:223)).EQ.0) XposT_1(i__) = rNULL__
if (idig(ar__1(225:231)).EQ.0) YposT_1(i__) = rNULL__
if (idig(ar__1(233:239)).EQ.0) XposH_1(i__) = rNULL__
if (idig(ar__1(241:247)).EQ.0) YposH_1(i__) = rNULL__
if (idig(ar__1(249:255)).EQ.0) XposS_1(i__) = rNULL__
if (idig(ar__1(257:263)).EQ.0) YposS_1(i__) = rNULL__
if (idig(ar__1(265:276)).EQ.0) ExtT_1(i__) = rNULL__
if (idig(ar__1(278:289)).EQ.0) e_ExtT_1(i__) = rNULL__
if (idig(ar__1(291:302)).EQ.0) ExtH_1(i__) = rNULL__
if (idig(ar__1(304:315)).EQ.0) e_ExtH_1(i__) = rNULL__
if (idig(ar__1(317:328)).EQ.0) ExtS_1(i__) = rNULL__
if (idig(ar__1(330:341)).EQ.0) e_ExtS_1(i__) = rNULL__
if (idig(ar__1(343:354)).EQ.0) DetMLT_1(i__) = rNULL__
if (idig(ar__1(356:367)).EQ.0) DetMLH_1(i__) = rNULL__
if (idig(ar__1(369:380)).EQ.0) DetMLS_1(i__) = rNULL__
if (idig(ar__1(382:393)).EQ.0) ExtMLT_1(i__) = rNULL__
if (idig(ar__1(395:406)).EQ.0) ExtMLH_1(i__) = rNULL__
if (idig(ar__1(408:419)).EQ.0) ExTMLS_1(i__) = rNULL__
if (idig(ar__1(421:432)).EQ.0) bgMapT_1(i__) = rNULL__
if (idig(ar__1(434:445)).EQ.0) bgMapH_1(i__) = rNULL__
if (idig(ar__1(447:458)).EQ.0) bgMapS_1(i__) = rNULL__
if (idig(ar__1(460:471)).EQ.0) ExpMapT_1(i__) = rNULL__
if (idig(ar__1(473:484)).EQ.0) ExpMapH_1(i__) = rNULL__
if (idig(ar__1(486:497)).EQ.0) ExpMapS_1(i__) = rNULL__
if (idig(ar__1(499:510)).EQ.0) FluxT_1(i__) = rNULL__
if (idig(ar__1(512:523)).EQ.0) e_FluxT_1(i__) = rNULL__
if (idig(ar__1(525:536)).EQ.0) FluxH_1(i__) = rNULL__
if (idig(ar__1(538:549)).EQ.0) e_FluxH_1(i__) = rNULL__
if (idig(ar__1(551:562)).EQ.0) FluxS_1(i__) = rNULL__
if (idig(ar__1(564:575)).EQ.0) e_FluxS_1(i__) = rNULL__
if (idig(ar__1(577:588)).EQ.0) RateT_1(i__) = rNULL__
if (idig(ar__1(590:601)).EQ.0) e_RateT_1(i__) = rNULL__
if (idig(ar__1(603:614)).EQ.0) RateH_1(i__) = rNULL__
if (idig(ar__1(616:627)).EQ.0) e_RateH_1(i__) = rNULL__
if (idig(ar__1(629:640)).EQ.0) RateS_1(i__) = rNULL__
if (idig(ar__1(642:653)).EQ.0) e_RateS_1(i__) = rNULL__
if (idig(ar__1(655:662)).EQ.0) RAdegT_1(i__) = rNULL__
if (idig(ar__1(664:671)).EQ.0) DEdegT_1(i__) = rNULL__
if (idig(ar__1(673:680)).EQ.0) e_posT_1(i__) = rNULL__
if (idig(ar__1(682:689)).EQ.0) GLONT_1(i__) = rNULL__
if (idig(ar__1(691:698)).EQ.0) GLATT_1(i__) = rNULL__
if (idig(ar__1(700:707)).EQ.0) RAdegH_1(i__) = rNULL__
if (idig(ar__1(709:716)).EQ.0) DEdegH_1(i__) = rNULL__
if (idig(ar__1(718:725)).EQ.0) e_posH_1(i__) = rNULL__
if (idig(ar__1(727:734)).EQ.0) GLONH_1(i__) = rNULL__
if (idig(ar__1(736:743)).EQ.0) GLATH_1(i__) = rNULL__
if (idig(ar__1(745:752)).EQ.0) RAdegS_1(i__) = rNULL__
if (idig(ar__1(754:761)).EQ.0) DEdegS_1(i__) = rNULL__
if (idig(ar__1(763:770)).EQ.0) e_posS_1(i__) = rNULL__
if (idig(ar__1(772:779)).EQ.0) GLONS_1(i__) = rNULL__
if (idig(ar__1(781:788)).EQ.0) GLATS_1(i__) = rNULL__
if(ar__1(1203:1214) .EQ. '') Dist_1(i__) = rNULL__
if (idig(ar__1(1216:1227)).EQ.0) DistRASS_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ XMMSL1_1(i__),ObsID_1(i__),SN_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),e_pos_1(i__),GLON_1(i__),GLAT_1(i__),HR1_1(i__),
+ e_HR1_1(i__),Date_1_1(i__),Date_2_1(i__),SCTST_1(i__),
+ e_SCTST_1(i__),SCTSH_1(i__),e_SCTSH_1(i__),SCTSS_1(i__),
+ e_SCTSS_1(i__),XposT_1(i__),YposT_1(i__),XposH_1(i__),
+ YposH_1(i__),XposS_1(i__),YposS_1(i__),ExtT_1(i__),
+ e_ExtT_1(i__),ExtH_1(i__),e_ExtH_1(i__),ExtS_1(i__),
+ e_ExtS_1(i__),DetMLT_1(i__),DetMLH_1(i__),DetMLS_1(i__),
+ ExtMLT_1(i__),ExtMLH_1(i__),ExTMLS_1(i__),bgMapT_1(i__),
+ bgMapH_1(i__),bgMapS_1(i__),ExpMapT_1(i__),ExpMapH_1(i__),
+ ExpMapS_1(i__),FluxT_1(i__),e_FluxT_1(i__),FluxH_1(i__),
+ e_FluxH_1(i__),FluxS_1(i__),e_FluxS_1(i__),RateT_1(i__),
+ e_RateT_1(i__),RateH_1(i__),e_RateH_1(i__),RateS_1(i__),
+ e_RateS_1(i__),RAdegT_1(i__),DEdegT_1(i__),e_posT_1(i__),
+ GLONT_1(i__),GLATT_1(i__),RAdegH_1(i__),DEdegH_1(i__),
+ e_posH_1(i__),GLONH_1(i__),GLATH_1(i__),RAdegS_1(i__),
+ DEdegS_1(i__),e_posS_1(i__),GLONS_1(i__),GLATS_1(i__),
+ VarFlag_1(i__),SRCNameT_1(i__),SRCNameH_1(i__),
+ SRCNameS_1(i__),XImNameT_1(i__),XImNameH_1(i__),
+ XImNameS_1(i__),VerINEXT_1(i__),VerHALO_1(i__),
+ VerHIBGND_1(i__),VerNREDG_1(i__),VerPSUSP_1(i__),
+ VerFALSE_1(i__),Comment_1(i__),Name_1(i__),OName_1(i__),
+ RASS_1(i__),Type_1(i__),r_Name_1(i__),Dist_1(i__),
+ DistRASS_1(i__),ModeID_1(i__),bgRate_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end