FORTRAN Generation
(/./ftp/cats/J/A_A/502/267)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/502/267 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.5, on 2013-May-23
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__=-1.e37)     	! NULL real number
      parameter  (iNULL__=-2147483647)	! NULL int  number

C=============================================================================
Cat. J/A+A/502/267     EWs of 31 giant stars of 10 open clusters   (Smiljanic+, 2009)
*================================================================================
*CNONa and 12C/13C in giant stars of 10 open clusters.
*    Smiljanic R., Gauderon R., North P., Barbuy B., Charbonnel C., Mowlavi N.
*   <Astron. Astrophys. 502, 267 (2009)>
*   =2009A&A...502..267S
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'stars.dat'	! Star list

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

      real*8        RA_        (nr__) ! (deg) Right Ascension J2000.0
      real*8        Dec        (nr__) ! (deg)     Declination J2000.0
      character*8   Cluster    (nr__) ! Cluster name
      integer*4     Star       (nr__) ! Star number within the cluster
      integer*4     RAh        (nr__) ! (h) Simbad Hour of Right Ascension (J2000.0)
      integer*4     RAm        (nr__) ! (min) Simbad Minute of Right Ascension (J2000.0)
      real*4        RAs        (nr__) ! (s) Simbad Second of Right Ascension (J2000.0)
      character*1   DE_        (nr__) ! Simbad Sign of the Declination (J2000.0)
      integer*4     DEd        (nr__) ! (deg) Simbad Degree of Declination (J2000.0)
      integer*4     DEm        (nr__) ! (arcmin) Simbad Arcminute of Declination (J2000.0)
      integer*4     DEs        (nr__) ! (arcsec) Simbad Arcsecond of Declination (J2000.0)
      real*4        Bmag       (nr__) ! (mag) ? Simbad B band magnitude
      real*4        Vmag       (nr__) ! (mag) Simbad V band magnitude
      character*7   SpType     (nr__) ! Simbad MK spectral type

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

C  Declarations for 'tablea1.dat'	! Equivalent widths of 31 giant stars

      integer*4 nr__1
      parameter (nr__1=4929)	! Number of records
      character*46 ar__1  	! Full-size record

      character*8   Cluster_1  (nr__1) ! Cluster name
      integer*4     Star_1     (nr__1) ! Star number within the cluster
      real*8        lambda     (nr__1) ! (0.1nm) Wavelength
      character*5   Ion        (nr__1) ! Element and ionisation
*                                   (I = neutral, II = ionised)
      real*4        EP         (nr__1) ! (eV) Excitation potential
      real*4        loggf      (nr__1) ! Oscillator strength
      real*4        EW         (nr__1) ! (0.1pm) ? Equivalent width

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

C  Loading file 'stars.dat'	! Star list

C  Format for file interpretation

    1 format(
     +  A8,1X,I3,2X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,2X,F6.3,1X,
     +  F6.3,1X,A7)

C  Effective file loading

      open(unit=1,file='stars.dat', status='old')
      write(6,*) '....Loading file: stars.dat'
      do i__=1,31
        read(1,'(A57)')ar__
        read(ar__,1)
     +  Cluster(i__),Star(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),Bmag(i__),Vmag(i__),SpType(i__)
        if(ar__(37:42) .EQ. '') Bmag(i__) = rNULL__
c       Derive coordinates RA_ and Dec from input data
c       (RA_ and Dec are set to rNULL_ when unknown)
        RA_(i__) = RAh(i__)
        if(RAh(i__) .GE. 0) RA_(i__)=RAh(i__)*15.
        if(RAm(i__) .GE. 0) RA_(i__)=RA_(i__)+RAm(i__)/4.
        if(RAs(i__) .GE. 0) RA_(i__)=RA_(i__)+RAs(i__)/240.
        Dec(i__) = DEd(i__)
        if(DEm(i__) .GE. 0) Dec(i__)=Dec(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) Dec(i__)=Dec(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c    ..............Just test output...........
        write(6,1)
     +  Cluster(i__),Star(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),Bmag(i__),Vmag(i__),SpType(i__)
        write(6,'(6H Pos: 2F8.4)') RA_(i__),Dec(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea1.dat'	! Equivalent widths of 31 giant stars

C  Format for file interpretation

    2 format(A8,1X,I3,1X,F7.2,2X,A5,1X,F4.2,2X,F6.3,1X,F5.1)

C  Effective file loading

      open(unit=1,file='tablea1.dat', status='old')
      write(6,*) '....Loading file: tablea1.dat'
      do i__=1,4929
        read(1,'(A46)')ar__1
        read(ar__1,2)
     +  Cluster_1(i__),Star_1(i__),lambda(i__),Ion(i__),EP(i__),
     +  loggf(i__),EW(i__)
        if(ar__1(42:46) .EQ. '') EW(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Cluster_1(i__),Star_1(i__),lambda(i__),Ion(i__),EP(i__),
     +  loggf(i__),EW(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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