Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/456/4488 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-Apr-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/MNRAS/456/4488 XQ-100 survey neutral gas (Sanchez-Ramirez+, 2016)
*================================================================================
*The evolution of neutral gas in damped Lyman {alpha} systems from the
*XQ-100 survey.
* Sanchez-Ramirez R., Ellison S.L., Prochaska J.X., Berg T.A.M., Lopez S.,
* D'odorico V., Becker G.D., Christensen L., Cupani G., Denney K.D.,
* Paris I., Worseck G., Gorosabel J.
* <Mon. Not. R. Astron. Soc., 456, 4488-4505 (2016)>
* =2016MNRAS.456.4488S (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! XQ-100 DLA catalog
integer*4 nr__
parameter (nr__=108) ! Number of records
character*116 ar__ ! Full-size record
character*2 n_Name (nr__) ! [ABC ] Note on Name (1)
character*10 Name (nr__) ! XQ-100 name (JHHMM+DDMM)
real*4 zem (nr__) ! ? Emission redshift
real*4 zmin (nr__) ! ? Minimum redshift range
real*4 zmax (nr__) ! ? Maximum redshift range
real*4 zabs (nr__) ! ? Absorption redshift
real*4 logN_HI (nr__) ! ([cm-2]) ? HI column density
real*4 e_logN_HI (nr__) ! ([cm-2]) ? rms uncertainty on logN(HI)
character*54 Lines (nr__) ! Lines covered
*Note (1): Note as follows:
* A = DLA already identified in a previous survey
* B = PDLA
* C = color biased sight-line
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Combined sample catalogue
integer*4 nr__1
parameter (nr__1=5874) ! Number of records
character*84 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*25 QSO (nr__1) ! QSO name
real*8 RAdeg (nr__1) ! (deg) Rigth ascension (J2000)
real*8 DEdeg (nr__1) ! (deg) Declination (J2000)
real*4 zem_1 (nr__1) ! ?=- Emission redshift
real*4 zmin_1 (nr__1) ! ?=- Minimum redshift range
real*4 zmax_1 (nr__1) ! ?=- Maximum redshift range
real*4 zabs_1 (nr__1) ! ?=- Absorption redshift
real*4 logN_HI_1 (nr__1) ! ([cm-2]) ?=- HI column density
real*4 e_logN_HI_1(nr__1) ! ([cm-2]) ?=- rms uncertainty on logN
C=============================================================================
C Loading file 'table1.dat' ! XQ-100 DLA catalog
C Format for file interpretation
1 format(
+ A2,1X,A10,1X,F6.4,2X,F6.4,2X,F6.4,2X,F6.4,2X,F5.2,3X,F4.2,4X,
+ A54)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,108
read(1,'(A116)')ar__
read(ar__,1)
+ n_Name(i__),Name(i__),zem(i__),zmin(i__),zmax(i__),zabs(i__),
+ logN_HI(i__),e_logN_HI(i__),Lines(i__)
if(ar__(15:20) .EQ. '') zem(i__) = rNULL__
if(ar__(23:28) .EQ. '') zmin(i__) = rNULL__
if(ar__(31:36) .EQ. '') zmax(i__) = rNULL__
if(ar__(39:44) .EQ. '') zabs(i__) = rNULL__
if(ar__(47:51) .EQ. '') logN_HI(i__) = rNULL__
if(ar__(55:58) .EQ. '') e_logN_HI(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ n_Name(i__),Name(i__),zem(i__),zmin(i__),zmax(i__),zabs(i__),
+ logN_HI(i__),e_logN_HI(i__),Lines(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Combined sample catalogue
C Format for file interpretation
2 format(
+ A25,5X,F9.5,1X,F8.4,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.4,1X,F5.2,
+ 1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,5874
read(1,'(A84)')ar__1
read(ar__1,2)
+ QSO(i__),RAdeg(i__),DEdeg(i__),zem_1(i__),zmin_1(i__),
+ zmax_1(i__),zabs_1(i__),logN_HI_1(i__),e_logN_HI_1(i__)
if (idig(ar__1(50:54)).EQ.0) zem_1(i__) = rNULL__
if (idig(ar__1(56:60)).EQ.0) zmin_1(i__) = rNULL__
if (idig(ar__1(62:66)).EQ.0) zmax_1(i__) = rNULL__
if (idig(ar__1(68:73)).EQ.0) zabs_1(i__) = rNULL__
if (idig(ar__1(75:79)).EQ.0) logN_HI_1(i__) = rNULL__
if (idig(ar__1(81:84)).EQ.0) e_logN_HI_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ QSO(i__),RAdeg(i__),DEdeg(i__),zem_1(i__),zmin_1(i__),
+ zmax_1(i__),zabs_1(i__),logN_HI_1(i__),e_logN_HI_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