Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/604/A119 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/A+A/604/A119 SHARDS S0 and E/S0 with anti-truncated disc (Borlaff+, 2017)
*================================================================================
*Evolution of the anti-truncated stellar profiles of S0 galaxies since
*z = 0.6 in the SHARDS survey. I. Sample and methods.
* Borlaff A., Eliche-Moral M.C., Beckman J.E., Ciambur B.C.,
* Perez-Gonzalez P.G., Barro G., Cava A., Cardiel N.
* <Astron. Astrophys., 604, A119 (2017)>
* =2017A&A...604A.119B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Properties of the initial red sample and
morphological classification
integer*4 nr__
parameter (nr__=150) ! Number of records
character*114 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq (nr__) ! [1/150] Sequential number
integer*4 SHARDS (nr__) ! SHARDS identification number
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
character*5 MType (nr__) ! Morphological type (1)
real*4 zphot (nr__) ! Photometric redshift
real*4 zspec (nr__) ! ?=- Spectroscopic redshift (when available)
real*4 logM_ (nr__) ! ([Msun]) Decimal logarithm of the stellar mass
real*4 e_logM_ (nr__) ! ([Msun]) rms uncertainty on logM*
character*1 l_SFR (nr__) ! [-] Upper limit flag (2)
real*4 SFR (nr__) ! (Msun/yr) Total SFR (2)
real*4 VMAG (nr__) ! (mag) Absolute synthetic rest-frame magnitude in
* Johnson V band from the Rainbow database
* (Sect. 2.5)
real*4 KsMAG (nr__) ! (mag) Absolute synthetic rest-frame magnitude in
* Ks band from the Rainbow database
* (Sect. 2.5)
real*4 AF775W (nr__) ! (10-2mag) Milky Way extinction in the F775W band
real*4 E_AF775W (nr__) ! (10-2mag) Error on AF775W (upper value)
real*4 e_AF775W_1 (nr__) ! (10-2mag) Error on AF775W (lower value)
real*4 R_F775W (nr__) ! (mag) ?=- Estimated K-correction for the F775W
* filters to the Steidel R band (3)
*Note (1): Morphological classification as follows:
* E = elliptical
* E/S0 = E/S0
* S0 = S0
* Sp = spiral
* CPSTB = compact post-starburst
* OM = ongoing merger
* GP = green pea (compact extremely star-forming galaxies)
* DF = diffuse galaxy
*Note (2): Total SFR (from the emission in 24um and 2800{AA} provided in the
* Rainbow database, following the procedure described in Barro et al.,
* 2013ApJ...765..104B). Negative values (in column l_SFR) correspond to upper
* limits following the notation on the Rainbow Database (see Sect. 3).
*Note (3): for the redshift of each galaxy considering the fit in Eq. (8)
* (see Sect. 2.5).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! Profile classification and photometric
parameters of S0 and E/S0 galaxies at 0.2<z<0.6
in the GOODS-N field
integer*4 nr__1
parameter (nr__1=44) ! Number of records
character*172 ar__1 ! Full-size record
integer*4 Seq_1 (nr__1) ! [1/44] Numerical ID in our final sample
* of S0 and E/S0 galaxies at 0.2<z<0.6
character*1 n_Seq (nr__1) ! [*] Note on Seq (1)
integer*4 SHARDS_1 (nr__1) ! SHARDS identification number
* (release DR2 beta)
character*3 Type (nr__1) ! Profile type
integer*4 i (nr__1) ! (deg) Galaxy inclination
real*4 mub (nr__1) ! (mag/arcsec+2) ?=- Surface brightness of the profile at
* the break radius {mu}_break_
real*4 E_mub (nr__1) ! (mag/arcsec+2) ? Error on mub (upper value)
real*4 e_mub_1 (nr__1) ! (mag/arcsec+2) ? Error on mub (lower value)
real*4 Rb (nr__1) ! (kpc) ?=- Break radius
real*4 E_Rb (nr__1) ! (kpc) ? Error on Rb (upper value)
real*4 e_Rb_1 (nr__1) ! (kpc) ? Error on Rb (lower value)
real*4 mu0i (nr__1) ! (mag/arcsec+2) Central surface brightness of the inner
* profile {mu}_0,i_
real*4 E_mu0i (nr__1) ! (mag/arcsec+2) ? Error on mu0i (upper value)
real*4 e_mu0i_1 (nr__1) ! (mag/arcsec+2) ? Error on mu0i (lower value)
real*4 hi (nr__1) ! (kpc) Scale-length of the inner profile
real*4 E_hi (nr__1) ! (kpc) Error on hi (upper value)
real*4 e_hi_1 (nr__1) ! (kpc) Error on hi (lower value)
real*4 mu0o (nr__1) ! (mag/arcsec+2) ?=- Central surface brightness of the
* outer profile {mu}_0,o_
real*4 E_mu0o (nr__1) ! (mag/arcsec+2) ? Error on mu0o (upper value)
real*4 e_mu0o_1 (nr__1) ! (mag/arcsec+2) ? Error on mu0o (lower value)
real*4 ho (nr__1) ! (kpc) ?=- Scale-length of the outer profile
real*4 E_ho (nr__1) ! (kpc) ? Error on ho (upper value)
real*4 e_ho_1 (nr__1) ! (kpc) ? Error on ho (lower value)
character*1 l_Ph (nr__1) ! Limit flag on Ph
real*4 Ph (nr__1) ! Elbow p-value for the likelihood that
* hi and ho are equal
character*1 l_Pmu0 (nr__1) ! Limit flag on Pmu0
real*4 Pmu0 (nr__1) ! Elbow p-value for for the likelihood
* that mu0o and mu0o are equal
real*4 R23 (nr__1) ! (kpc) Isophotal radius for 23mag/arcsec^2^ of
* the object in the corrected Steidel R
* band
real*4 E_R23 (nr__1) ! (kpc) Error on R23 (upper value)
real*4 e_R23_1 (nr__1) ! (kpc) Error on R23 (lower value)
real*4 Rlim (nr__1) ! (kpc) Limiting radius of the profile
*Note (1): * indicates objects in the AGN subsample (see Sect. 2.8).
C=============================================================================
C Loading file 'tablea1.dat' ! Properties of the initial red sample and
* morphological classification
C Format for file interpretation
1 format(
+ I3,7X,I8,1X,F8.4,1X,F7.4,1X,A5,1X,F4.2,1X,F6.4,1X,F6.3,1X,
+ F5.3,1X,A1,F6.3,1X,F6.2,1X,F6.2,1X,F6.3,1X,F5.3,1X,F5.3,1X,
+ F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,150
read(1,'(A114)')ar__
read(ar__,1)
+ Seq(i__),SHARDS(i__),RAdeg(i__),DEdeg(i__),MType(i__),
+ zphot(i__),zspec(i__),logM_(i__),e_logM_(i__),l_SFR(i__),
+ SFR(i__),VMAG(i__),KsMAG(i__),AF775W(i__),E_AF775W(i__),
+ e_AF775W_1(i__),R_F775W(i__)
if (idig(ar__(48:53)).EQ.0) zspec(i__) = rNULL__
if (idig(ar__(109:114)).EQ.0) R_F775W(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq(i__),SHARDS(i__),RAdeg(i__),DEdeg(i__),MType(i__),
+ zphot(i__),zspec(i__),logM_(i__),e_logM_(i__),l_SFR(i__),
+ SFR(i__),VMAG(i__),KsMAG(i__),AF775W(i__),E_AF775W(i__),
+ e_AF775W_1(i__),R_F775W(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb1.dat' ! Profile classification and photometric
* parameters of S0 and E/S0 galaxies at 0.2<z<0.6
* in the GOODS-N field
C Format for file interpretation
2 format(
+ I2,A1,7X,I8,1X,A3,1X,I2,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.3,1X,F5.3,1X,F5.3,A1,E6.2,1X,A1,
+ E6.2,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,44
read(1,'(A172)')ar__1
read(ar__1,2)
+ Seq_1(i__),n_Seq(i__),SHARDS_1(i__),Type(i__),i(i__),mub(i__),
+ E_mub(i__),e_mub_1(i__),Rb(i__),E_Rb(i__),e_Rb_1(i__),
+ mu0i(i__),E_mu0i(i__),e_mu0i_1(i__),hi(i__),E_hi(i__),
+ e_hi_1(i__),mu0o(i__),E_mu0o(i__),e_mu0o_1(i__),ho(i__),
+ E_ho(i__),e_ho_1(i__),l_Ph(i__),Ph(i__),l_Pmu0(i__),Pmu0(i__),
+ R23(i__),E_R23(i__),e_R23_1(i__),Rlim(i__)
if (idig(ar__1(27:32)).EQ.0) mub(i__) = rNULL__
if(ar__1(34:38) .EQ. '') E_mub(i__) = rNULL__
if(ar__1(40:44) .EQ. '') e_mub_1(i__) = rNULL__
if (idig(ar__1(46:50)).EQ.0) Rb(i__) = rNULL__
if(ar__1(52:55) .EQ. '') E_Rb(i__) = rNULL__
if(ar__1(57:60) .EQ. '') e_Rb_1(i__) = rNULL__
if(ar__1(69:73) .EQ. '') E_mu0i(i__) = rNULL__
if(ar__1(75:79) .EQ. '') e_mu0i_1(i__) = rNULL__
if (idig(ar__1(99:103)).EQ.0) mu0o(i__) = rNULL__
if(ar__1(105:108) .EQ. '') E_mu0o(i__) = rNULL__
if(ar__1(110:113) .EQ. '') e_mu0o_1(i__) = rNULL__
if (idig(ar__1(115:119)).EQ.0) ho(i__) = rNULL__
if(ar__1(121:125) .EQ. '') E_ho(i__) = rNULL__
if(ar__1(127:131) .EQ. '') e_ho_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),n_Seq(i__),SHARDS_1(i__),Type(i__),i(i__),mub(i__),
+ E_mub(i__),e_mub_1(i__),Rb(i__),E_Rb(i__),e_Rb_1(i__),
+ mu0i(i__),E_mu0i(i__),e_mu0i_1(i__),hi(i__),E_hi(i__),
+ e_hi_1(i__),mu0o(i__),E_mu0o(i__),e_mu0o_1(i__),ho(i__),
+ E_ho(i__),e_ho_1(i__),l_Ph(i__),Ph(i__),l_Pmu0(i__),Pmu0(i__),
+ R23(i__),E_R23(i__),e_R23_1(i__),Rlim(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