CDS.headStuff2 FORTRAN Generation
()

Conversion of standardized ReadMe file for catalog 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-Mar-28
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/546/694       NICMOS star formation history            (Thompson +, 2001)
*================================================================================
*Star formation history in the NICMOS Northern Hubble Deep Field.
*     Thompson R.I., Weymann R.J., Storrie-Lombardi L.J.
*    <Astrophys. J. 546, 694 (2001)>
*    =2001ApJ...546..694T
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table1.dat'	! Listing of measured quantities

      integer*4 nr__
      parameter (nr__=281)	! Number of records
      character*104 ar__   	! Full-size record

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg      (nr__) ! (deg) Right Ascension J2000
      real*8        DEdeg      (nr__) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      real*8        NICMOS     (nr__) ! NICMOS identification number (1)
      character*8   v_WBD96_   (nr__) ! WFPC identification number from
*                                     Williams et al., 1996, Cat. <J/AJ/112/1335>
      real*4        z          (nr__) ! Redshift
      real*4        E_B_V      (nr__) ! (mag) (B-V) color excess
      real*8        SFR        (nr__) ! (solMass/yr) Star formation rate determined in Section 5
      real*4        Lum        (nr__) ! (solLum) Bolometric luminosity (2)
      real*4        Frac       (nr__) ! Fraction of the luminosity removed by
*                                     extinction and re-emitted in the
*                                     mid & far-IR
      real*4        SCUBA      (nr__) ! (mJy) ? The SCUBA 850 microns flux
      real*4        Temp       (nr__) ! Template number of best fit between
*                                      1.0 (early-cold) and 6.0 (late-hot)
      real*4        Chi2       (nr__) ! Modified Chi^2^ of the fit from Equ. 1 (3)
      real*4        Totmag     (nr__) ! (mag) ? The F160W AB magnitude (4)
      real*4        Apmag      (nr__) ! (mag) ? The 0.6" diameter aperture F160W magnitude
*                                       See note (5)
      integer*4     RAh        (nr__) ! (h) Right Ascension (J2000)
      integer*4     RAm        (nr__) ! (min) Right Ascension (J2000)
      real*4        RAs        (nr__) ! (s) Right Ascension (J2000)
      character*1   DE_        (nr__) ! Declination (J2000)
      integer*4     DEd        (nr__) ! (deg) Declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Declination (J2000)
      real*4        DEs        (nr__) ! (arcsec) Declination (J2000)
*Note (1): Objects with ID numbers less than 1000 are objects that are
*    identified with sources in the original KFOCAS catalog of Thompson et
*    al., 1999, Cat. <J/AJ/117/17>, [TSW99] NNN.NNN in Simbad.
*    Objects with catalog numbers greater than or equal to 1000 are
*    SExtractor objects that did not match in position with an original
*    catalog entry to within 0.3".
*Note (2): The flux of a galaxy is obtained by integrating over the
*    unextincted selected template scaled by the factor A from equation 2.
*    The bolometric luminosity then follows from the redshift and our
*    adopted cosmology.
*Note (3): It should be noted that the distribution of the modified chi^2^
*    values will not rigorously follow a true chi^2^ distribution. The
*    values are provided to give a qualitative indication of the relative
*    goodness of fit for the best fit values of redshift, template type and
*    E(B-V) from object to object.
*Note (4): Derived from the aperture correction method described in
*    Yan et al., 1998ApJ...503L..19Y.
*Note (5): If no magnitude is listed the object had a zero or negative
*    measured F160W flux.

C=============================================================================

C  Loading file 'table1.dat'	! Listing of measured quantities

C  Format for file interpretation

    1 format(
     +  F10.5,1X,A8,1X,F4.2,1X,F4.2,1X,F7.3,1X,E8.2,1X,F4.2,1X,E8.2,
     +  1X,F3.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,
     +  1X,I2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,281
        read(1,'(A104)')ar__
        read(ar__,1)
     +  NICMOS(i__),v_WBD96_(i__),z(i__),E_B_V(i__),SFR(i__),Lum(i__),
     +  Frac(i__),SCUBA(i__),Temp(i__),Chi2(i__),Totmag(i__),
     +  Apmag(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__)
        if(ar__(53:60) .EQ. '') SCUBA(i__) = rNULL__
        if(ar__(71:74) .EQ. '') Totmag(i__) = rNULL__
        if(ar__(76:79) .EQ. '') Apmag(i__) = rNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  NICMOS(i__),v_WBD96_(i__),z(i__),E_B_V(i__),SFR(i__),Lum(i__),
     +  Frac(i__),SCUBA(i__),Temp(i__),Chi2(i__),Totmag(i__),
     +  Apmag(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end