Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/609/A130 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 2024-Apr-20
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/609/A130 [CII] luminosities of galaxies in G.A.S.+Cloudy (Lagache+, 2018)
*================================================================================
*The [CII] 158 micron line emission in high-redshift galaxies.
* Lagache G., Cousin M., Chatzikos M.
* <Astron. Astrophys. 609, A130 (2018)>
* =2018A&A...609A.130L (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gcciiz4.dat' ! Outputs from the model for z=4
integer*4 nr__
parameter (nr__=10803) ! Number of records
character*126 ar__ ! Full-size record
real*8 logMstar ! ([Msun]) Stellar mass (log_M_star)
real*8 logMgas ! ([Msun]) Gas mass (log_Mgas)
real*8 ZgOHindex ! Gaz metallicity in OH index unit
* (Zg_OH_index)
real*8 ZgmZ ! (Sun) Gaz metallicity (Zg_mZ)
real*8 logISRF ! Interstellar radiation field in G0 unit
* (log_ISRF)
real*8 lognHI ! ([cm-3]) HI density (log_nHI)
real*8 logSFR ! ([Msun/yr]) Star formation rate from the model (log_SFR)
real*8 logSFRobs ! ([Msun/yr]) Star formation rate computed from the UV and
* IR luminosities (log_SFR_obs)
real*8 logLIR ! ([Lsun]) 8-1000 micron luminosity (log_L_IR)
real*8 logLCII ! ([Lsun]) CII luminosity (log_L_CII)
real*8 logLCIIw ! ([Lsun]) CII luminosity without CMB attenuation
* (log_L_CII_woa)
real*8 logLCIIww ! ([Lsun]) CII luminosity without CMB attenuation and
* heating (log_L_CII_woa_woh)
real*8 logLCIIVa ! ([Lsun]) CII luminosity using N=2 (Eq. 11)
* (log_L_CII_Vallini)
real*8 logLAGN ! ([Lsun]) ? AGN bolometric luminosity (log_L_agn)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gcciiz5.dat' ! Outputs from the model for z=4.7
integer*4 nr__1
parameter (nr__1=6175) ! Number of records
character*126 ar__1 ! Full-size record
real*8 logMstar_1 ! ([Msun]) Stellar mass (log_M_star)
real*8 logMgas_1 ! ([Msun]) Gas mass (log_Mgas)
real*8 ZgOHindex_1 ! Gaz metallicity in OH index unit
* (Zg_OH_index)
real*8 ZgmZ_1 ! (Sun) Gaz metallicity (Zg_mZ)
real*8 logISRF_1 ! Interstellar radiation field in G0 unit
* (log_ISRF)
real*8 lognHI_1 ! ([cm-3]) HI density (log_nHI)
real*8 logSFR_1 ! ([Msun/yr]) Star formation rate from the model (log_SFR)
real*8 logSFRobs_1 ! ([Msun/yr]) Star formation rate computed from the UV and
* IR luminosities (log_SFR_obs)
real*8 logLIR_1 ! ([Lsun]) 8-1000 micron luminosity (log_L_IR)
real*8 logLCII_1 ! ([Lsun]) CII luminosity (log_L_CII)
real*8 logLCIIw_1 ! ([Lsun]) CII luminosity without CMB attenuation
* (log_L_CII_woa)
real*8 logLCIIww_1 ! ([Lsun]) CII luminosity without CMB attenuation and
* heating (log_L_CII_woa_woh)
real*8 logLCIIVa_1 ! ([Lsun]) CII luminosity using N=2 (Eq. 11)
* (log_L_CII_Vallini)
real*8 logLAGN_1 ! ([Lsun]) ? AGN bolometric luminosity (log_L_agn)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gcciiz6.dat' ! Outputs from the model for z=5.9
integer*4 nr__2
parameter (nr__2=2374) ! Number of records
character*126 ar__2 ! Full-size record
real*8 logMstar_2 ! ([Msun]) Stellar mass (log_M_star)
real*8 logMgas_2 ! ([Msun]) Gas mass (log_Mgas)
real*8 ZgOHindex_2 ! Gaz metallicity in OH index unit
* (Zg_OH_index)
real*8 ZgmZ_2 ! (Sun) Gaz metallicity (Zg_mZ)
real*8 logISRF_2 ! Interstellar radiation field in G0 unit
* (log_ISRF)
real*8 lognHI_2 ! ([cm-3]) HI density (log_nHI)
real*8 logSFR_2 ! ([Msun/yr]) Star formation rate from the model (log_SFR)
real*8 logSFRobs_2 ! ([Msun/yr]) Star formation rate computed from the UV and
* IR luminosities (log_SFR_obs)
real*8 logLIR_2 ! ([Lsun]) 8-1000 micron luminosity (log_L_IR)
real*8 logLCII_2 ! ([Lsun]) CII luminosity (log_L_CII)
real*8 logLCIIw_2 ! ([Lsun]) CII luminosity without CMB attenuation
* (log_L_CII_woa)
real*8 logLCIIww_2 ! ([Lsun]) CII luminosity without CMB attenuation and
* heating (log_L_CII_woa_woh)
real*8 logLCIIVa_2 ! ([Lsun]) CII luminosity using N=2 (Eq. 11)
* (log_L_CII_Vallini)
real*8 logLAGN_2 ! ([Lsun]) ? AGN bolometric luminosity (log_L_agn)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gcciiz7.dat' ! Outputs from the model for z=6.7
integer*4 nr__3
parameter (nr__3=1074) ! Number of records
character*126 ar__3 ! Full-size record
real*8 logMstar_3 ! ([Msun]) Stellar mass (log_M_star)
real*8 logMgas_3 ! ([Msun]) Gas mass (log_Mgas)
real*8 ZgOHindex_3 ! Gaz metallicity in OH index unit
* (Zg_OH_index)
real*8 ZgmZ_3 ! (Sun) Gaz metallicity (Zg_mZ)
real*8 logISRF_3 ! Interstellar radiation field in G0 unit
* (log_ISRF)
real*8 lognHI_3 ! ([cm-3]) HI density (log_nHI)
real*8 logSFR_3 ! ([Msun/yr]) Star formation rate from the model (log_SFR)
real*8 logSFRobs_3 ! ([Msun/yr]) Star formation rate computed from the UV and
* IR luminosities (log_SFR_obs)
real*8 logLIR_3 ! ([Lsun]) 8-1000 micron luminosity (log_L_IR)
real*8 logLCII_3 ! ([Lsun]) CII luminosity (log_L_CII)
real*8 logLCIIw_3 ! ([Lsun]) CII luminosity without CMB attenuation
* (log_L_CII_woa)
real*8 logLCIIww_3 ! ([Lsun]) CII luminosity without CMB attenuation and
* heating (log_L_CII_woa_woh)
real*8 logLCIIVa_3 ! ([Lsun]) CII luminosity using N=2 (Eq. 11)
* (log_L_CII_Vallini)
real*8 logLAGN_3 ! ([Lsun]) ? AGN bolometric luminosity (log_L_agn)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gcciiz8.dat' ! Outputs from the model for z=7.6
integer*4 nr__4
parameter (nr__4=388) ! Number of records
character*126 ar__4 ! Full-size record
real*8 logMstar_4 ! ([Msun]) Stellar mass (log_M_star)
real*8 logMgas_4 ! ([Msun]) Gas mass (log_Mgas)
real*8 ZgOHindex_4 ! Gaz metallicity in OH index unit
* (Zg_OH_index)
real*8 ZgmZ_4 ! (Sun) Gaz metallicity (Zg_mZ)
real*8 logISRF_4 ! Interstellar radiation field in G0 unit
* (log_ISRF)
real*8 lognHI_4 ! ([cm-3]) HI density (log_nHI)
real*8 logSFR_4 ! ([Msun/yr]) Star formation rate from the model (log_SFR)
real*8 logSFRobs_4 ! ([Msun/yr]) Star formation rate computed from the UV and
* IR luminosities (log_SFR_obs)
real*8 logLIR_4 ! ([Lsun]) 8-1000 micron luminosity (log_L_IR)
real*8 logLCII_4 ! ([Lsun]) CII luminosity (log_L_CII)
real*8 logLCIIw_4 ! ([Lsun]) CII luminosity without CMB attenuation
* (log_L_CII_woa)
real*8 logLCIIww_4 ! ([Lsun]) CII luminosity without CMB attenuation and
* heating (log_L_CII_woa_woh)
real*8 logLCIIVa_4 ! ([Lsun]) CII luminosity using N=2 (Eq. 11)
* (log_L_CII_Vallini)
real*8 logLAGN_4 ! ([Lsun]) ? AGN bolometric luminosity (log_L_agn)
C=============================================================================
C Loading file 'gcciiz4.dat' ! Outputs from the model for z=4
C Format for file interpretation
1 format(
+ F8.5,1X,F8.5,1X,F8.5,1X,F9.6,1X,F7.5,1X,F7.5,1X,F9.6,1X,F9.6,
+ 1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.6)
C Effective file loading
open(unit=1,status='old',file=
+'gcciiz4.dat')
write(6,*) '....Loading file: gcciiz4.dat'
do i__=1,10803
read(1,'(A126)')ar__
read(ar__,1)
+ logMstar,logMgas,ZgOHindex,ZgmZ,logISRF,lognHI,logSFR,
+ logSFRobs,logLIR,logLCII,logLCIIw,logLCIIww,logLCIIVa,logLAGN
if(ar__(119:126) .EQ. '') logLAGN = rNULL__
c ..............Just test output...........
write(6,1)
+ logMstar,logMgas,ZgOHindex,ZgmZ,logISRF,lognHI,logSFR,
+ logSFRobs,logLIR,logLCII,logLCIIw,logLCIIww,logLCIIVa,logLAGN
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'gcciiz5.dat' ! Outputs from the model for z=4.7
C Format for file interpretation
2 format(
+ F8.5,1X,F8.5,1X,F8.5,1X,F9.6,1X,F7.5,1X,F7.5,1X,F9.6,1X,F9.6,
+ 1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.6)
C Effective file loading
open(unit=1,status='old',file=
+'gcciiz5.dat')
write(6,*) '....Loading file: gcciiz5.dat'
do i__=1,6175
read(1,'(A126)')ar__1
read(ar__1,2)
+ logMstar_1,logMgas_1,ZgOHindex_1,ZgmZ_1,logISRF_1,lognHI_1,
+ logSFR_1,logSFRobs_1,logLIR_1,logLCII_1,logLCIIw_1,
+ logLCIIww_1,logLCIIVa_1,logLAGN_1
if(ar__1(119:126) .EQ. '') logLAGN_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ logMstar_1,logMgas_1,ZgOHindex_1,ZgmZ_1,logISRF_1,lognHI_1,
+ logSFR_1,logSFRobs_1,logLIR_1,logLCII_1,logLCIIw_1,
+ logLCIIww_1,logLCIIVa_1,logLAGN_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'gcciiz6.dat' ! Outputs from the model for z=5.9
C Format for file interpretation
3 format(
+ F8.5,1X,F8.5,1X,F8.5,1X,F9.6,1X,F7.5,1X,F7.5,1X,F9.6,1X,F9.6,
+ 1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.6)
C Effective file loading
open(unit=1,status='old',file=
+'gcciiz6.dat')
write(6,*) '....Loading file: gcciiz6.dat'
do i__=1,2374
read(1,'(A126)')ar__2
read(ar__2,3)
+ logMstar_2,logMgas_2,ZgOHindex_2,ZgmZ_2,logISRF_2,lognHI_2,
+ logSFR_2,logSFRobs_2,logLIR_2,logLCII_2,logLCIIw_2,
+ logLCIIww_2,logLCIIVa_2,logLAGN_2
if(ar__2(119:126) .EQ. '') logLAGN_2 = rNULL__
c ..............Just test output...........
write(6,3)
+ logMstar_2,logMgas_2,ZgOHindex_2,ZgmZ_2,logISRF_2,lognHI_2,
+ logSFR_2,logSFRobs_2,logLIR_2,logLCII_2,logLCIIw_2,
+ logLCIIww_2,logLCIIVa_2,logLAGN_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'gcciiz7.dat' ! Outputs from the model for z=6.7
C Format for file interpretation
4 format(
+ F8.5,1X,F8.5,1X,F8.5,1X,F9.6,1X,F7.5,1X,F7.5,1X,F9.6,1X,F9.6,
+ 1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.6)
C Effective file loading
open(unit=1,status='old',file=
+'gcciiz7.dat')
write(6,*) '....Loading file: gcciiz7.dat'
do i__=1,1074
read(1,'(A126)')ar__3
read(ar__3,4)
+ logMstar_3,logMgas_3,ZgOHindex_3,ZgmZ_3,logISRF_3,lognHI_3,
+ logSFR_3,logSFRobs_3,logLIR_3,logLCII_3,logLCIIw_3,
+ logLCIIww_3,logLCIIVa_3,logLAGN_3
if(ar__3(119:126) .EQ. '') logLAGN_3 = rNULL__
c ..............Just test output...........
write(6,4)
+ logMstar_3,logMgas_3,ZgOHindex_3,ZgmZ_3,logISRF_3,lognHI_3,
+ logSFR_3,logSFRobs_3,logLIR_3,logLCII_3,logLCIIw_3,
+ logLCIIww_3,logLCIIVa_3,logLAGN_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'gcciiz8.dat' ! Outputs from the model for z=7.6
C Format for file interpretation
5 format(
+ F8.5,1X,F8.5,1X,F8.5,1X,F9.6,1X,F7.5,1X,F7.5,1X,F9.6,1X,F9.6,
+ 1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.6)
C Effective file loading
open(unit=1,status='old',file=
+'gcciiz8.dat')
write(6,*) '....Loading file: gcciiz8.dat'
do i__=1,388
read(1,'(A126)')ar__4
read(ar__4,5)
+ logMstar_4,logMgas_4,ZgOHindex_4,ZgmZ_4,logISRF_4,lognHI_4,
+ logSFR_4,logSFRobs_4,logLIR_4,logLCII_4,logLCIIw_4,
+ logLCIIww_4,logLCIIVa_4,logLAGN_4
if(ar__4(119:126) .EQ. '') logLAGN_4 = rNULL__
c ..............Just test output...........
write(6,5)
+ logMstar_4,logMgas_4,ZgOHindex_4,ZgmZ_4,logISRF_4,lognHI_4,
+ logSFR_4,logSFRobs_4,logLIR_4,logLCII_4,logLCIIw_4,
+ logLCIIww_4,logLCIIVa_4,logLAGN_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end