Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/472/4878 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-Sep-12
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/472/4878 Number counts predictions for surveys (Valentino+, 2017)
*================================================================================
*Predicting emission line fluxes and number counts of distant galaxies for
*cosmological surveys.
* Valentino F., Daddi E., Silverman J.D., Puglisi A., Kashino D., Renzini A.,
* Cimatti A., Pozzetti L., Rodighiero G., Pannella M., Gobat R., Zamorani G.
* <Mon. Not. R. Astron. Soc., 472, 4878-4899 (2017)>
* =2017MNRAS.472.4878V (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cosmos.dat' ! Photometric parameters and line flux predictions
for 31193 galaxies in COSMOS at 1.4<z<1.8
integer*4 nr__
parameter (nr__=31193) ! Number of records
character*171 ar__ ! Full-size record
C J2000.0 position composed of: RAdeg DEdeg
integer*4 ID ! Galaxy identification number (1)
real*8 RAdeg ! (deg) Right ascension (J2000.0)
real*8 DEdeg ! (deg) Declination (J2000.0)
real*4 logM ! ([Msun]) Logarithmic stellar mass (2)
real*4 logSFR ! ([Msun/yr]) Logarithmic star formation rate
real*4 zph ! Photometric redshift
real*4 AV ! (mag) V-band stellar extinction
real*4 FHa ! (10-19W/m2) Predicted flux for H{alpha} (G1)
real*4 F_OII_ ! (10-19W/m2) Predicted flux for [OII] (G1)
real*4 F_OIII_ ! (10-19W/m2) Predicted flux for [OIII] (G1)
real*4 FHb ! (10-19W/m2) Predicted flux for H{beta} (G1)
real*4 fabs ! H{beta} Balmer absorption factor
*Note (1): The galaxy identification number is from Laigle et al., 2016,
* Cat. J/ApJS/224/24 .
*Note (2): The logarithmic stellar mass is from Laigle et al., 2016,
* Cat. J/ApJS/224/24, converted to a Salpeter (1955ApJ...121..161S)
* Initial Mass Function adding 0.23dex.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'goodss.dat' ! Photometric parameters and line flux predictions
for 3858 galaxies in GOODS-S at 1.4<z<1.8
integer*4 nr__1
parameter (nr__1=3858) ! Number of records
character*120 ar__1 ! Full-size record
C J2000.0 position composed of: RAdeg DEdeg
integer*4 ID_1 ! Galaxy identification number
real*8 RAdeg_1 ! (deg) Right ascension (J2000.0)
real*8 DEdeg_1 ! (deg) Declination (J2000.0)
real*4 logM_1 ! ([Msun]) Logarithmic stellar mass (1)
real*4 logSFR_1 ! ([Msun/yr]) Logarithmic star formation rate
real*4 zphot ! Photometric redshift
real*4 AV_1 ! (mag) V-band stellar extinction
real*4 FHa_1 ! (10-19W/m2) Predicted flux for H{alpha} (G1)
real*4 F_OII__1 ! (10-19W/m2) Predicted flux for [OII] (G1)
real*4 F_OIII__1 ! (10-19W/m2) Predicted flux for [OIII] (G1)
real*4 FHb_1 ! (10-19W/m2) Predicted flux for H{beta} (G1)
real*4 fabs_1 ! H{beta} Balmer absorption factor
*Note (1): Salpeter (1955ApJ...121..161S) Initial Mass Function.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! H{alpha} number counts from COSMOS and GOODS-S
photometric samples
GOODS-S photometric samples at 1.4<z<1.8
integer*4 nr__2
parameter (nr__2=14) ! Number of records
character*133 ar__2 ! Full-size record
real*4 Flux ! (10-19W/m2) Limiting flux (G1)
integer*4 COSct1_4z1_8 ! (deg-2) Average counts at 1.4<z<1.8
* for the COSMOS sample (COS_aver_1.4z1.8)
integer*4 COSco1_4z1_8 ! (deg-2) Convolution error on counts at 1.4<z<1.8
* for the COSMOS sample (COS_conv_1.4z1.8)
integer*4 COSpo1_4z1_8 ! (deg-2) Poissonian error on counts at 1.4<z<1.8
* for the COSMOS sample (COS_pois_1.4z1.8)
integer*4 COSmc1_4z1_8 ! (deg-2) Monte Carlo error on counts at 1.4<z<1.8
* for the COSMOS sample (COS_mc_1.4z1.8)
integer*4 COSct0_9z1_8 ! (deg-2) Average counts at 0.9<z<1.8
* for the COSMOS sample (COS_aver_0.9z1.8)
integer*4 COSco0_9z1_8 ! (deg-2) Convolution error on counts at 0.9<z<1.8
* for the COSMOS sample (COS_conv_0.9z1.8)
integer*4 COSpo0_9z1_8 ! (deg-2) Poissonian error on counts at 0.9<z<1.8
* for the COSMOS sample (COS_pois_0.9z1.8)
integer*4 COSmc0_9z1_8 ! (deg-2) Monte Carlo error on counts at 0.9<z<1.8
* for the COSMOS sample (COS_mc_0.9z1.8)
integer*4 GOOct1_4z1_8 ! (deg-2) Average counts at 1.4<z<1.8
* for the GOODS-S sample (GOO_aver_1.4z1.8)
integer*4 GOOco1_4z1_8 ! (deg-2) Convolution error on counts at 1.4<z<1.8
* for the GOODS-S sample (GOO_conv_1.4z1.8)
integer*4 GOOpo1_4z1_8 ! (deg-2) Poissonian error on counts at 1.4<z<1.8
* for the GOODS-S sample (GOO_pois_1.4z1.8)
integer*4 GOOmc1_4z1_8 ! (deg-2) Monte Carlo error on counts at 1.4<z<1.8
* for the GOODS-S sample (GOO_mc_1.4z1.8)
integer*4 GOOct0_9z1_8 ! (deg-2) Average counts at 0.9<z<1.8
* for the GOODS-S sample (GOO_aver_0.9z1.8)
integer*4 GOOco0_9z1_8 ! (deg-2) Convolution error on counts at 0.9<z<1.8
* for the GOODS-S sample (GOO_conv_0.9z1.8)
integer*4 GOOpo0_9z1_8 ! (deg-2) Poissonian error on counts at 0.9<z<1.8
* for the GOODS-S sample (GOO_pois_0.9z1.8)
integer*4 GOOmc0_9z1_8 ! (deg-2) Monte Carlo error on counts at 0.9<z<1.8
* for the GOODS-S sample (GOO_mc_0.9z1.8)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! [OII] and [OIII] number counts from COSMOS and
integer*4 nr__3
parameter (nr__3=9) ! Number of records
character*134 ar__3 ! Full-size record
real*4 Flux_1 ! (10-19W/m2) Limiting flux (G1)
integer*4 COSctOII ! (deg-2) Average [OII] counts at 1.4<z<1.8
* for the COSMOS sample (COS_aver_oii)
integer*4 COScoOII ! (deg-2) Convolution error on [OII] counts at
* 1.4<z<1.8 for the COSMOS sample
* (COS_conv_oii)
integer*4 COSpoOII ! (deg-2) Poissonian error on [OII] counts at 1.4<z<1.8
* for the COSMOS sample (COS_pois_oii)
integer*4 COSmcOII ! (deg-2) Monte Carlo error on [OII] counts at
* 1.4<z<1.8 for the COSMOS sample (COS_mc_oii)
integer*4 COSctOIII ! (deg-2) Average [OIII] counts at 1.4<z<1.8
* for the COSMOS sample (COS_aver_oiii)
integer*4 COScoOIII ! (deg-2) Convolution error on [OIII] counts
* at 1.4<z<1.8 for the COSMOS sample
* (COS_conv_oiii)
integer*4 COSpoOIII ! (deg-2) Poissonian error on [OIII] counts at
* 1.4<z<1.8 for the COSMOS sample
* (COS_pois_oiii)
integer*4 COSmcOIII ! (deg-2) Monte Carlo error on [OIII] counts at
* 1.4<z<1.8 for the COSMOS sample
* (COS_mc_oiii)
integer*4 GOOctOII ! (deg-2) Average [OII] counts at 1.4<z<1.8
* for the GOODS-S sample (GOO_aver_oii)
integer*4 GOOcoOII ! (deg-2) Convolution error on [OII] counts at
* 1.4<z<1.8 for the GOODS-S sample
* (GOO_conv_oii)
integer*4 GOOpoOII ! (deg-2) Poissonian error on [OII] counts at 1.4<z<1.8
* for the GOODS-S sample (GOO_pois_oii)
integer*4 GOOmcOII ! (deg-2) Monte Carlo error on [OII] counts at
* 1.4<z<1.8 for the GOODS-S sample
* (GOO_mc_oii)
integer*4 GOOctOIII ! (deg-2) Average [OIII] counts at 1.4<z<1.8
* for the GOODS-S sample (GOO_aver_oiii)
integer*4 GOOcoOIII ! (deg-2) Convolution error on [OIII] counts at
* 1.4<z<1.8 for the GOODS-S sample
* (GOO_conv_oiii)
integer*4 GOOpoOIII ! (deg-2) Poissonian error on [OIII] counts at
* 1.4<z<1.8 for the GOODS-S sample
* (GOO_pois_oiii)
integer*4 GOOmcOIII ! (deg-2) Monte Carlo error on [OIII] counts at
* 1.4<z<1.8 for the GOODS-S sample
* (GOO_mc_oiii)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! H{alpha}+[NII] number counts from COSMOS and
GOODS-S photometric samples
integer*4 nr__4
parameter (nr__4=14) ! Number of records
character*101 ar__4 ! Full-size record
real*4 Flux_2 ! (10-19W/m2) Limiting flux (G1)
integer*4 COSct1_4z1_8_1 ! (deg-2) Average counts at 1.4<z<1.8
* for the COSMOS sample (COS_aver_1.4z1.8)
integer*4 COSco1_4z1_8_1 ! (deg-2) Convolution error on counts at 1.4<z<1.8
* for the COSMOS sample (COSco1.4z1.8)
integer*4 COSpo1_4z1_8_1 ! (deg-2) Poissonian error on counts at 1.4<z<1.8
* for the COSMOS sample (COSpo1.4z1.8)
integer*4 COSct0_9z1_8_1 ! (deg-2) Average counts at 0.9<z<1.8
* for the COSMOS sample (COS_aver_0.9z1.8)
integer*4 COSco0_9z1_8_1 ! (deg-2) Convolution error on counts at 0.9<z<1.8
* for the COSMOS sample (COSco0.9z1.8)
integer*4 COSpo0_9z1_8_1 ! (deg-2) Poissonian error on counts at 0.9<z<1.8
* for the COSMOS sample (COSpo0.9z1.8)
integer*4 GOOct1_4z1_8_1 ! (deg-2) Average counts at 1.4<z<1.8
* for the GOODS-S sample (GOO_aver_1.4z1.8)
integer*4 GOOco1_4z1_8_1 ! (deg-2) Convolution error on counts at 1.4<z<1.8
* for the GOODS-S sample (GOOco1.4z1.8)
integer*4 GOOpo1_4z1_8_1 ! (deg-2) Poissonian error on counts at 1.4<z<1.8
* for the GOODS-S sample (GOOpo1.4z1.8)
integer*4 GOOct0_9z1_8_1 ! (deg-2) Average counts at 0.9<z<1.8
* for the GOODS-S sample (GOO_aver_0.9z1.8)
integer*4 GOOco0_9z1_8_1 ! (deg-2) Convolution error on counts at 0.9<z<1.8
* for the GOODS-S sample (GOOco0.9z1.8)
integer*4 GOOpo0_9z1_8_1 ! (deg-2) Poissonian error on counts at 0.9<z<1.8
* for the GOODS-S sample (GOOpo0.9z1.8)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! H{alpha}, [OII], and [OIII] number counts for
starbursting galaxies from the GOODS-S
photometric samples at 1.4<z<1.8 and 0.9<z<1.8
integer*4 nr__5
parameter (nr__5=14) ! Number of records
character*92 ar__5 ! Full-size record
real*4 Flux_3 ! (10-19W/m2) Limiting flux (G1)
integer*4 SBct1_4z1_8ha ! (deg-2) Average H{alpha} counts at 1.4<z<1.8
* for the GOODS-S sample
* (SB_aver_1.4z1.8_ha)
integer*4 SBco1_4z1_8ha ! (deg-2) Convolution error on H{alpha} counts at
* 1.4<z<1.8 for the GOODS-S sample
* (SB_conv_1.4z1.8_ha)
integer*4 SBpo1_4z1_8ha ! (deg-2) Poissonian error H{alpha} counts at
* 1.4<z<1.8 for the GOODS-S sample
* (SB_pois_1.4z1.8_ha)
integer*4 SBct1_4z1_8OII ! (deg-2) Average [OII] counts at 1.4<z<1.8
* for the GOODS-S sample
* (SB_aver_1.4z1.8_oii)
integer*4 SBco1_4z1_8OII ! (deg-2) Convolution error on [OII] counts at
* 1.4<z<1.8 for the GOODS-S sample
* (SB_conv_1.4z1.8_oii)
integer*4 SBpo1_4z1_8OII ! (deg-2) Poissonian error on [OII] counts at
* 1.4<z<1.8 for the GOODS-S sample
* (SB_pois_1.4z1.8_oii)
integer*4 SBct1_4z1_8OIII ! (deg-2) Average [OIII] counts at 1.4<z<1.8 for
* the GOODS-S sample
* (SB_aver_1.4z1.8_oiii)
integer*4 SBco1_4z1_8OIII ! (deg-2) Convolution error on [OIII] counts at
* 1.4<z<1.8 for the GOODS-S sample
* (SB_conv_1.4z1.8_oiii)
integer*4 SBpo1_4z1_8OIII ! (deg-2) Poissonian error on [OIII] counts at
* 1.4<z<1.8 for the GOODS-S sample
* (SB_pois_1.4z1.8_oiii)
integer*4 SBct0_9z1_8ha ! (deg-2) Average Ha counts at 0.9<z<1.8 for the
* GOODS-S sample (SB_aver_0.9z1.8_ha)
integer*4 SBco0_9z1_8ha ! (deg-2) Convolution error on Ha counts at
* 0.9<z<1.8 for the GOODS-S sample
* (SB_conv_0.9z1.8_ha)
integer*4 SBpo0_9z1_8ha ! (deg-2) Poissonian error on Ha counts at
* 0.9<z<1.8 for the GOODS-S sample
* (SB_pois_0.9z1.8_ha)
C=============================================================================
C Loading file 'cosmos.dat' ! Photometric parameters and line flux predictions
* for 31193 galaxies in COSMOS at 1.4<z<1.8
C Format for file interpretation
1 format(
+ I6,3X,F12.8,5X,F10.8,10X,F5.2,9X,F6.2,11X,F4.2,10X,F5.2,11X,
+ F4.2,11X,F4.2,11X,F4.2,11X,F4.2,11X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'cosmos.dat')
write(6,*) '....Loading file: cosmos.dat'
do i__=1,31193
read(1,'(A171)')ar__
read(ar__,1)
+ ID,RAdeg,DEdeg,logM,logSFR,zph,AV,FHa,F_OII_,F_OIII_,FHb,fabs
c ..............Just test output...........
write(6,1)
+ ID,RAdeg,DEdeg,logM,logSFR,zph,AV,FHa,F_OII_,F_OIII_,FHb,fabs
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'goodss.dat' ! Photometric parameters and line flux predictions
* for 3858 galaxies in GOODS-S at 1.4<z<1.8
C Format for file interpretation
2 format(
+ I4,1X,F11.8,2X,F12.8,5X,F5.2,5X,F5.2,6X,F4.2,6X,F4.2,6X,F4.2,
+ 6X,F4.2,6X,F4.2,6X,F4.2,6X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'goodss.dat')
write(6,*) '....Loading file: goodss.dat'
do i__=1,3858
read(1,'(A120)')ar__1
read(ar__1,2)
+ ID_1,RAdeg_1,DEdeg_1,logM_1,logSFR_1,zphot,AV_1,FHa_1,
+ F_OII__1,F_OIII__1,FHb_1,fabs_1
c ..............Just test output...........
write(6,2)
+ ID_1,RAdeg_1,DEdeg_1,logM_1,logSFR_1,zphot,AV_1,FHa_1,
+ F_OII__1,F_OIII__1,FHb_1,fabs_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table1.dat' ! H{alpha} number counts from COSMOS and GOODS-S
* photometric samples
* GOODS-S photometric samples at 1.4<z<1.8
C Format for file interpretation
3 format(
+ F4.2,4X,I4,4X,I4,4X,I4,4X,I4,4X,I5,4X,I4,4X,I4,4X,I4,4X,I5,4X,
+ I3,4X,I4,4X,I4,4X,I5,4X,I3,4X,I4,4X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,14
read(1,'(A133)')ar__2
read(ar__2,3)
+ Flux,COSct1_4z1_8,COSco1_4z1_8,COSpo1_4z1_8,COSmc1_4z1_8,
+ COSct0_9z1_8,COSco0_9z1_8,COSpo0_9z1_8,COSmc0_9z1_8,
+ GOOct1_4z1_8,GOOco1_4z1_8,GOOpo1_4z1_8,GOOmc1_4z1_8,
+ GOOct0_9z1_8,GOOco0_9z1_8,GOOpo0_9z1_8,GOOmc0_9z1_8
c ..............Just test output...........
write(6,3)
+ Flux,COSct1_4z1_8,COSco1_4z1_8,COSpo1_4z1_8,COSmc1_4z1_8,
+ COSct0_9z1_8,COSco0_9z1_8,COSpo0_9z1_8,COSmc0_9z1_8,
+ GOOct1_4z1_8,GOOco1_4z1_8,GOOpo1_4z1_8,GOOmc1_4z1_8,
+ GOOct0_9z1_8,GOOco0_9z1_8,GOOpo0_9z1_8,GOOmc0_9z1_8
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! [OII] and [OIII] number counts from COSMOS and
C Format for file interpretation
4 format(
+ F4.2,4X,I4,4X,I4,4X,I4,4X,I4,4X,I4,4X,I4,4X,I4,4X,I4,4X,I5,4X,
+ I4,4X,I4,4X,I4,4X,I5,4X,I4,4X,I4,4X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,9
read(1,'(A134)')ar__3
read(ar__3,4)
+ Flux_1,COSctOII,COScoOII,COSpoOII,COSmcOII,COSctOIII,
+ COScoOIII,COSpoOIII,COSmcOIII,GOOctOII,GOOcoOII,GOOpoOII,
+ GOOmcOII,GOOctOIII,GOOcoOIII,GOOpoOIII,GOOmcOIII
c ..............Just test output...........
write(6,4)
+ Flux_1,COSctOII,COScoOII,COSpoOII,COSmcOII,COSctOIII,
+ COScoOIII,COSpoOIII,COSmcOIII,GOOctOII,GOOcoOII,GOOpoOII,
+ GOOmcOII,GOOctOIII,GOOcoOIII,GOOpoOIII,GOOmcOIII
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! H{alpha}+[NII] number counts from COSMOS and
* GOODS-S photometric samples
C Format for file interpretation
5 format(
+ F4.2,4X,I5,4X,I4,4X,I4,4X,I5,4X,I4,4X,I4,4X,I5,4X,I3,4X,I3,4X,
+ I5,4X,I3,4X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,14
read(1,'(A101)')ar__4
read(ar__4,5)
+ Flux_2,COSct1_4z1_8_1,COSco1_4z1_8_1,COSpo1_4z1_8_1,
+ COSct0_9z1_8_1,COSco0_9z1_8_1,COSpo0_9z1_8_1,GOOct1_4z1_8_1,
+ GOOco1_4z1_8_1,GOOpo1_4z1_8_1,GOOct0_9z1_8_1,GOOco0_9z1_8_1,
+ GOOpo0_9z1_8_1
c ..............Just test output...........
write(6,5)
+ Flux_2,COSct1_4z1_8_1,COSco1_4z1_8_1,COSpo1_4z1_8_1,
+ COSct0_9z1_8_1,COSco0_9z1_8_1,COSpo0_9z1_8_1,GOOct1_4z1_8_1,
+ GOOco1_4z1_8_1,GOOpo1_4z1_8_1,GOOct0_9z1_8_1,GOOco0_9z1_8_1,
+ GOOpo0_9z1_8_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! H{alpha}, [OII], and [OIII] number counts for
* starbursting galaxies from the GOODS-S
* photometric samples at 1.4<z<1.8 and 0.9<z<1.8
C Format for file interpretation
6 format(
+ F4.2,4X,I4,4X,I2,4X,I2,4X,I4,4X,I4,4X,I4,4X,I4,4X,I4,4X,I4,4X,
+ I4,4X,I2,4X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,14
read(1,'(A92)')ar__5
read(ar__5,6)
+ Flux_3,SBct1_4z1_8ha,SBco1_4z1_8ha,SBpo1_4z1_8ha,
+ SBct1_4z1_8OII,SBco1_4z1_8OII,SBpo1_4z1_8OII,SBct1_4z1_8OIII,
+ SBco1_4z1_8OIII,SBpo1_4z1_8OIII,SBct0_9z1_8ha,SBco0_9z1_8ha,
+ SBpo0_9z1_8ha
c ..............Just test output...........
write(6,6)
+ Flux_3,SBct1_4z1_8ha,SBco1_4z1_8ha,SBpo1_4z1_8ha,
+ SBct1_4z1_8OII,SBco1_4z1_8OII,SBpo1_4z1_8OII,SBct1_4z1_8OIII,
+ SBco1_4z1_8OIII,SBpo1_4z1_8OIII,SBct0_9z1_8ha,SBco0_9z1_8ha,
+ SBpo0_9z1_8ha
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end