Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/676/286 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-Jun-08
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/ApJ/676/286 Spectral templates for galaxies from 0.2 to 10um (Assef+, 2008)
*================================================================================
*Low-resolution spectral templates for galaxies from 0.2 to 10 {mu}m.
* Assef R.J., Kochanek C.S., Brodwin M., Brown M.J.I., Caldwell N.,
* Cool R.J., Eisenhardt P., Eisenstein D., Gonzalez A.H., Jannuzi B.T.,
* Jones C., McKenzie E., Murray S.S., Stern D.
* <Astrophys. J., 676, 286-303 (2008)>
* =2008ApJ...676..286A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Three template model absolute magnitudes
integer*4 nr__
parameter (nr__=453) ! Number of records
character*158 ar__ ! Full-size record
real*4 z ! Redshift
integer*4 t ! [1/4] Template number: 1=E, 2=Sbc, 3=Im,
* 4=E+A (post-starburst galaxy)
real*8 BWMAG ! (mag) Absolute B_W_ band magnitude (2)
real*8 BMAG ! (mag) Absolute B band magnitude (2)
real*8 VMAG ! (mag) Absolute V band magnitude (2)
real*8 RMAG ! (mag) Absolute R band magnitude (2)
real*8 IMAG ! (mag) Absolute I band magnitude (2)
real*8 uMAG ! (mag) Absolute SDSS u' band AB magnitude (2)
real*8 gMAG ! (mag) Absolute SDSS g' band AB magnitude (2)
real*8 rMAG_1 ! (mag) Absolute SDSS r' band AB magnitude (2)
real*8 iMAG_1 ! (mag) Absolute SDSS i' band AB magnitude (2)
real*8 zMAG ! (mag) Absolute SDSS z' band AB magnitude (2)
real*8 JMAG ! (mag) Absolute 2MASS J band magnitude (2)
real*8 HMAG ! (mag) Absolute 2MASS H band magnitude (2)
real*8 KSMAG ! (mag) Absolute 2MASS K_S_ band magnitude (2)
real*8 KMAG ! (mag) Absolute K band magnitude (2)
real*8 v3_6MAG ! (mag) Absolute Spitzer/IRAC 3.6 micron band magnitude (2)
real*8 v4_5MAG ! (mag) Absolute Spitzer/IRAC 4.5 micron band magnitude (2)
real*8 v5_8MAG ! (mag) Absolute Spitzer/IRAC 5.8 micron band magnitude (2)
real*8 v8_0MAG ! (mag) Absolute Spitzer/IRAC 8.0 micron band magnitude (2)
real*4 DM ! (mag) ? Distance Modulus
*Note (2): The absolute magnitude we present here corresponds to the
* canonical definition of the absolute magnitude (as in, for example,
* eq. 26 of Hogg 1999, astro-ph/9905116) plus the K correction term.
* This allows the calculation of photometric redshifts and K corrections
* from the table. To determine photometric redshifts, colors should be
* calculated and matched to the data by varying the a_k_ coefficients
* (see Section 3.1) and the redshift. For a galaxy at redshift z with
* template coefficients a_k_, the model magnitude in band b is given by
* M_b(z)_=-2.5log[sum_k_(a_k_ 10^-0.4M_b,k_(z)^)]. Apparent magnitudes
* can be determined by adding the distance modulus to the absolute ones.
* To determine K corrections for a galaxy at redshift z, coefficients
* a_k_ should also be determined to match the observed colors as above.
* With the same coefficients, redshift z and redshift zero model
* absolute magnitudes can be determined, and the difference between them
* will correspond to the desired K correction.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Four template model absolute magnitudes
integer*4 nr__1
parameter (nr__1=604) ! Number of records
character*158 ar__1 ! Full-size record
real*4 z_1 ! Redshift
integer*4 t_1 ! [1/4] Template number: 1=E, 2=Sbc, 3=Im,
* 4=E+A (post-starburst galaxy)
real*8 BWMAG_1 ! (mag) Absolute B_W_ band magnitude (2)
real*8 BMAG_1 ! (mag) Absolute B band magnitude (2)
real*8 VMAG_1 ! (mag) Absolute V band magnitude (2)
real*8 RMAG_2 ! (mag) Absolute R band magnitude (2)
real*8 IMAG_2 ! (mag) Absolute I band magnitude (2)
real*8 uMAG_1 ! (mag) Absolute SDSS u' band AB magnitude (2)
real*8 gMAG_1 ! (mag) Absolute SDSS g' band AB magnitude (2)
real*8 rMAG_3 ! (mag) Absolute SDSS r' band AB magnitude (2)
real*8 iMAG_3 ! (mag) Absolute SDSS i' band AB magnitude (2)
real*8 zMAG_1 ! (mag) Absolute SDSS z' band AB magnitude (2)
real*8 JMAG_1 ! (mag) Absolute 2MASS J band magnitude (2)
real*8 HMAG_1 ! (mag) Absolute 2MASS H band magnitude (2)
real*8 KSMAG_1 ! (mag) Absolute 2MASS K_S_ band magnitude (2)
real*8 KMAG_1 ! (mag) Absolute K band magnitude (2)
real*8 v3_6MAG_1 ! (mag) Absolute Spitzer/IRAC 3.6 micron band magnitude (2)
real*8 v4_5MAG_1 ! (mag) Absolute Spitzer/IRAC 4.5 micron band magnitude (2)
real*8 v5_8MAG_1 ! (mag) Absolute Spitzer/IRAC 5.8 micron band magnitude (2)
real*8 v8_0MAG_1 ! (mag) Absolute Spitzer/IRAC 8.0 micron band magnitude (2)
real*4 DM_1 ! (mag) ? Distance Modulus
*Note (2): The absolute magnitude we present here corresponds to the
* canonical definition of the absolute magnitude (as in, for example,
* eq. 26 of Hogg 1999, astro-ph/9905116) plus the K correction term.
* This allows the calculation of photometric redshifts and K corrections
* from the table. To determine photometric redshifts, colors should be
* calculated and matched to the data by varying the a_k_ coefficients
* (see Section 3.1) and the redshift. For a galaxy at redshift z with
* template coefficients a_k_, the model magnitude in band b is given by
* M_b(z)_=-2.5log[sum_k_(a_k_ 10^-0.4M_b,k_(z)^)]. Apparent magnitudes
* can be determined by adding the distance modulus to the absolute ones.
* To determine K corrections for a galaxy at redshift z, coefficients
* a_k_ should also be determined to match the observed colors as above.
* With the same coefficients, redshift z and redshift zero model
* absolute magnitudes can be determined, and the difference between them
* will correspond to the desired K correction.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Three component spectral templates
integer*4 nr__2
parameter (nr__2=160) ! Number of records
character*52 ar__2 ! Full-size record
real*8 lambda ! (um) Wavelength
real*4 E ! (mW/m2/Hz) Best fit Elliptical template F_nu_ (1)
real*4 Sbc ! (mW/m2/Hz) Best fit Sbc Spiral Sbc template F_nu_ (1)
real*4 Im ! (mW/m2/Hz) Best fit Irregular Im template F_nu_ (1)
real*4 E_A ! (mW/m2/Hz) ? Best fit E+A template F_nu_ (table4 only )(1)
*Note (1): In units of erg/s/cm^2^/Hz (=mW/m^2^/Hz=10^23^Jy). Templates
* are normalized to be at a distance of 10pc and to have an integrated
* luminosity between the wavelength boundaries of 10^10^L_{sun}_.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Four component spectral templates
integer*4 nr__3
parameter (nr__3=160) ! Number of records
character*52 ar__3 ! Full-size record
real*8 lambda_1 ! (um) Wavelength
real*4 E_1 ! (mW/m2/Hz) Best fit Elliptical template F_nu_ (1)
real*4 Sbc_1 ! (mW/m2/Hz) Best fit Sbc Spiral Sbc template F_nu_ (1)
real*4 Im_1 ! (mW/m2/Hz) Best fit Irregular Im template F_nu_ (1)
real*4 E_A_1 ! (mW/m2/Hz) ? Best fit E+A template F_nu_ (table4 only )(1)
*Note (1): In units of erg/s/cm^2^/Hz (=mW/m^2^/Hz=10^23^Jy). Templates
* are normalized to be at a distance of 10pc and to have an integrated
* luminosity between the wavelength boundaries of 10^10^L_{sun}_.
C=============================================================================
C Loading file 'table1.dat' ! Three template model absolute magnitudes
C Format for file interpretation
1 format(
+ F5.3,1X,I1,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,
+ F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,
+ 1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,453
read(1,'(A158)')ar__
read(ar__,1)
+ z,t,BWMAG,BMAG,VMAG,RMAG,IMAG,uMAG,gMAG,rMAG_1,iMAG_1,zMAG,
+ JMAG,HMAG,KSMAG,KMAG,v3_6MAG,v4_5MAG,v5_8MAG,v8_0MAG,DM
if(ar__(153:158) .EQ. '') DM = rNULL__
c ..............Just test output...........
write(6,1)
+ z,t,BWMAG,BMAG,VMAG,RMAG,IMAG,uMAG,gMAG,rMAG_1,iMAG_1,zMAG,
+ JMAG,HMAG,KSMAG,KMAG,v3_6MAG,v4_5MAG,v5_8MAG,v8_0MAG,DM
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Four template model absolute magnitudes
C Format for file interpretation
2 format(
+ F5.3,1X,I1,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,
+ F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,
+ 1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,604
read(1,'(A158)')ar__1
read(ar__1,2)
+ z_1,t_1,BWMAG_1,BMAG_1,VMAG_1,RMAG_2,IMAG_2,uMAG_1,gMAG_1,
+ rMAG_3,iMAG_3,zMAG_1,JMAG_1,HMAG_1,KSMAG_1,KMAG_1,v3_6MAG_1,
+ v4_5MAG_1,v5_8MAG_1,v8_0MAG_1,DM_1
if(ar__1(153:158) .EQ. '') DM_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ z_1,t_1,BWMAG_1,BMAG_1,VMAG_1,RMAG_2,IMAG_2,uMAG_1,gMAG_1,
+ rMAG_3,iMAG_3,zMAG_1,JMAG_1,HMAG_1,KSMAG_1,KMAG_1,v3_6MAG_1,
+ v4_5MAG_1,v5_8MAG_1,v8_0MAG_1,DM_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Three component spectral templates
C Format for file interpretation
3 format(F8.6,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,160
read(1,'(A52)')ar__2
read(ar__2,3)lambda,E,Sbc,Im,E_A
if(ar__2(43:52) .EQ. '') E_A = rNULL__
c ..............Just test output...........
write(6,3)lambda,E,Sbc,Im,E_A
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Four component spectral templates
C Format for file interpretation
4 format(F8.6,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,160
read(1,'(A52)')ar__3
read(ar__3,4)lambda_1,E_1,Sbc_1,Im_1,E_A_1
if(ar__3(43:52) .EQ. '') E_A_1 = rNULL__
c ..............Just test output...........
write(6,4)lambda_1,E_1,Sbc_1,Im_1,E_A_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end