Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/449/951 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 2026-Jun-10
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/449/951 GOODS-MUSIC sample: multicolour catalog (Grazian+, 2006)
*================================================================================
*The GOODS-MUSIC sample:
*a multicolour catalog of near-IR selected galaxies in the GOODS-South field.
* Grazian A., Fontana A., De Santis C., Nonino M., Salimbeni S., Giallongo E.,
* Cristiani S., Gallozzi S., Vanzella E.
* <Astron. Astrophys. 449, 951 (2006)>
* =2006A&A...449..951G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Coordinates, spectroscopy and magnitude limits
integer*4 nr__
parameter (nr__=18296) ! Number of records
character*90 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq (nr__) ! Identification number
real*8 RAdeg (nr__) ! (deg) Right ascension, in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination, in decimal degrees (J2000)
real*4 zspec (nr__) ! ?=-1.0 Spectroscopic Redshift
character*9 SpClass (nr__) ! Spectroscopic class (1)
character*10 r_SpClass (nr__) ! Spectroscopic Catalog (2)
integer*4 q_zspec (nr__) ! ?=99 Quality of Spectroscopic Redshift (3)
real*4 zphot (nr__) ! ?=-1.0 Photometric redshift
integer*4 POS (nr__) ! Position flag (4)
integer*4 star (nr__) ! Star flag (5)
integer*4 AGN (nr__) ! AGN flag (6)
real*4 zlim (nr__) ! z band magnitude limit (7)
real*4 kslim (nr__) ! (mag) Ks band magnitude limit (8)
real*4 S_G (nr__) ! Star-galaxy classifier of SExtractor
* (0 for galaxy, 1 for star)
*Note (1): Spectroscopic classification as follows:
* GALAXY = galaxy
* EARLY = early-type galaxy
* EMISSION = late-type galaxy
* COMPOSITE = early+late type galaxy
* AGN = Active Galactic Nucleus
* STAR = star
* UNKNOWN = spectrum not available
*Note (2): Reference spectroscopic catalog as follows:
* GOODSV1.0 = Vanzella et al., 2005, Cat. <J/A+A/434/53>
* K20 = Mignoli et al., 2005, Cat. <J/A+A/437/883>
* CXO-CDFS = Szokoly et al., 2004, Cat. <J/ApJS/155/271>
* VVDS = Le Fevre et al., 2004A&A...428..1043L
* MASTER = http://www.eso.org/science/goods/spectroscopy/CDFS_Mastercat/
* LCIRS = Doherty et al., 2005MNRAS.361..525D
* COMBO-17 = Wolf et al., 2004A&A...421..913W, Cat. <II/253>
* DADDI = Daddi et al., 2005ApJ...626..680D
* SIC = Cristiani et al., 2000A&A...359..489C
* UNKNOWN = spectrum not available.
*Note (3): Quality of spectroscopic redshift as follows:
* 0 = very good
* 1 = good
* 2 = uncertain
* 3 = bad quality
* 99 = not available
*Note (4): Position flag as follows:
* 1 = inside GOODS-MUSIC area
* 0 = outside GOODS-MUSIC area
*Note (5): Star flag on the basis of spectroscopy, morphology,
* and BzK colours as follows:
* 1 = probable star
* 0 = no star (A galaxy should have star flag=0 and AGN flag=0)
*Note (6): AGN flag, based only on spectroscopy, as follows:
* 1 = probable AGN
* 0 = no AGN (A galaxy should have star flag=0 and AGN flag=0)
*Note (7): Magnitude limit in the z band in 1arcsec^2^ and at S/N=1.
*Note (8): Magnitude limit in the Ks band in 1arcsec^2^ and at S/N=1.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Optical photometry
integer*4 nr__1
parameter (nr__1=18296) ! Number of records
character*110 ar__1 ! Full-size record
integer*4 Seq_1 (nr__1) ! Identification number
character*1 l_U35mag (nr__1) ! Limit flag on U35mag
real*4 U35mag (nr__1) ! (mag) ?=99.000 U35 magnitude, in AB photometric system
character*1 l_U38mag (nr__1) ! Limit flag on U38mag
real*4 U38mag (nr__1) ! (mag) ?=99.000 U38 magnitude, in AB photometric system
character*1 l_UVIMmag (nr__1) ! Limit flag on UVIMmag
real*4 UVIMmag (nr__1) ! (mag) ?=99.000 U-VIMOS magnitude, in AB
* photometric system
character*1 l_Bmag (nr__1) ! Limit flag on Bmag
real*4 Bmag (nr__1) ! (mag) ?=99.000 B magnitude, in AB photometric system
character*1 l_Vmag (nr__1) ! Limit flag on Vmag
real*4 Vmag (nr__1) ! (mag) ?=99.000 V magnitude, in AB photometric system
character*1 l_imag (nr__1) ! Limit flag on imag
real*4 imag (nr__1) ! (mag) ?=99.000 i magnitude, in AB photometric system
character*1 l_zmag (nr__1) ! Limit flag on zmag
real*4 zmag (nr__1) ! (mag) ?=99.000 z magnitude, in AB photometric system
real*4 e_U35mag (nr__1) ! (mag) ?=99.000 U35 magnitude error
real*4 e_U38mag (nr__1) ! (mag) ?=99.000 U38 magnitude error
real*4 e_UVIMmag (nr__1) ! (mag) ?=99.000 U-VIMOS magnitude error
real*4 e_Bmag (nr__1) ! (mag) ?=99.000 B magnitude error
real*4 e_Vmag (nr__1) ! (mag) ?=99.000 V magnitude error
real*4 e_imag (nr__1) ! (mag) ?=99.000 i magnitude error
real*4 e_zmag (nr__1) ! (mag) ?=99.000 z magnitude error
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Infrared photometry
integer*4 nr__2
parameter (nr__2=18296) ! Number of records
character*110 ar__2 ! Full-size record
integer*4 Seq_2 (nr__2) ! Identification number
character*1 l_Jmag (nr__2) ! Limit flag on Jmag
real*4 Jmag (nr__2) ! (mag) ?=99.000 J magnitude, in AB photometric system
character*1 l_Hmag (nr__2) ! Limit flag on Hmag
real*4 Hmag (nr__2) ! (mag) ?=99.000 H magnitude, in AB photometric system
character*1 l_Ksmag (nr__2) ! Limit flag on Ksmag
real*4 Ksmag (nr__2) ! (mag) ?=99.000 Ks magnitude, in AB photometric system
character*1 l_IR36mag (nr__2) ! Limit flag on IR36mag
real*4 IR36mag (nr__2) ! (mag) ?=99.000 IRAC 3.6{mu}m magnitude,
* in AB photometric system
character*1 l_IR45mag (nr__2) ! Limit flag on IR45mag
real*4 IR45mag (nr__2) ! (mag) ?=99.000 IRAC 4.5{mu}m magnitude,
* in AB photometric system
character*1 l_IR58mag (nr__2) ! Limit flag on IR58mag
real*4 IR58mag (nr__2) ! (mag) ?=99.000 IRAC 5.8{mu}m magnitude,
* in AB photometric system
character*1 l_IR80mag (nr__2) ! Limit flag on IR80mag
real*4 IR80mag (nr__2) ! (mag) ?=99.000 IRAC 8.0{mu}m magnitude,
* in AB photometric system
real*4 e_Jmag (nr__2) ! (mag) ?=99.000 J magnitude error
real*4 e_Hmag (nr__2) ! (mag) ?=99.000 H magnitude error
real*4 e_Ksmag (nr__2) ! (mag) ?=99.000 Ks magnitude error
real*4 e_IR36mag (nr__2) ! (mag) ?=99.000 IRAC 3.6 micron magnitude error
real*4 e_IR45mag (nr__2) ! (mag) ?=99.000 IRAC 4.5 micron magnitude error
real*4 e_IR58mag (nr__2) ! (mag) ?=99.000 IRAC 5.8 micron magnitude error
real*4 e_IR80mag (nr__2) ! (mag) ?=99.000 IRAC 8.0 micron magnitude error
C=============================================================================
C Loading file 'table5.dat' ! Coordinates, spectroscopy and magnitude limits
C Format for file interpretation
1 format(
+ I5,1X,F9.6,1X,F10.6,1X,F6.3,1X,A9,1X,A10,1X,I2,1X,F6.3,1X,I1,
+ 1X,I1,1X,I1,1X,F6.3,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,18296
read(1,'(A90)')ar__
read(ar__,1)
+ Seq(i__),RAdeg(i__),DEdeg(i__),zspec(i__),SpClass(i__),
+ r_SpClass(i__),q_zspec(i__),zphot(i__),POS(i__),star(i__),
+ AGN(i__),zlim(i__),kslim(i__),S_G(i__)
c ..............Just test output...........
write(6,1)
+ Seq(i__),RAdeg(i__),DEdeg(i__),zspec(i__),SpClass(i__),
+ r_SpClass(i__),q_zspec(i__),zphot(i__),POS(i__),star(i__),
+ AGN(i__),zlim(i__),kslim(i__),S_G(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Optical photometry
C Format for file interpretation
2 format(
+ I5,1X,A1,F6.3,1X,A1,F6.3,1X,A1,F6.3,1X,A1,F6.3,1X,A1,F6.3,1X,
+ A1,F6.3,1X,A1,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
+ F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,18296
read(1,'(A110)')ar__1
read(ar__1,2)
+ Seq_1(i__),l_U35mag(i__),U35mag(i__),l_U38mag(i__),
+ U38mag(i__),l_UVIMmag(i__),UVIMmag(i__),l_Bmag(i__),Bmag(i__),
+ l_Vmag(i__),Vmag(i__),l_imag(i__),imag(i__),l_zmag(i__),
+ zmag(i__),e_U35mag(i__),e_U38mag(i__),e_UVIMmag(i__),
+ e_Bmag(i__),e_Vmag(i__),e_imag(i__),e_zmag(i__)
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),l_U35mag(i__),U35mag(i__),l_U38mag(i__),
+ U38mag(i__),l_UVIMmag(i__),UVIMmag(i__),l_Bmag(i__),Bmag(i__),
+ l_Vmag(i__),Vmag(i__),l_imag(i__),imag(i__),l_zmag(i__),
+ zmag(i__),e_U35mag(i__),e_U38mag(i__),e_UVIMmag(i__),
+ e_Bmag(i__),e_Vmag(i__),e_imag(i__),e_zmag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Infrared photometry
C Format for file interpretation
3 format(
+ I5,1X,A1,F6.3,1X,A1,F6.3,1X,A1,F6.3,1X,A1,F6.3,1X,A1,F6.3,1X,
+ A1,F6.3,1X,A1,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
+ F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,18296
read(1,'(A110)')ar__2
read(ar__2,3)
+ Seq_2(i__),l_Jmag(i__),Jmag(i__),l_Hmag(i__),Hmag(i__),
+ l_Ksmag(i__),Ksmag(i__),l_IR36mag(i__),IR36mag(i__),
+ l_IR45mag(i__),IR45mag(i__),l_IR58mag(i__),IR58mag(i__),
+ l_IR80mag(i__),IR80mag(i__),e_Jmag(i__),e_Hmag(i__),
+ e_Ksmag(i__),e_IR36mag(i__),e_IR45mag(i__),e_IR58mag(i__),
+ e_IR80mag(i__)
c ..............Just test output...........
write(6,3)
+ Seq_2(i__),l_Jmag(i__),Jmag(i__),l_Hmag(i__),Hmag(i__),
+ l_Ksmag(i__),Ksmag(i__),l_IR36mag(i__),IR36mag(i__),
+ l_IR45mag(i__),IR45mag(i__),l_IR58mag(i__),IR58mag(i__),
+ l_IR80mag(i__),IR80mag(i__),e_Jmag(i__),e_Hmag(i__),
+ e_Ksmag(i__),e_IR36mag(i__),e_IR45mag(i__),e_IR58mag(i__),
+ e_IR80mag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end