Conversion of standardized ReadMe file for
file /./ftp/cats/IV/26 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-Apr-17 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. IV/26 Hercules Messier 13 (M13) photometry (Ludendorff, 1905) *================================================================================ *Der grosse Sternhaufen im Hercules Messier 13. * Ludendorff H. * <Publ. Astrophys. Obs. Potsdam, 50, 1 (1905)> * =1905POPot..50....1L C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'catalog.dat' ! Catalog of stars in Messier 13 integer*4 nr__ parameter (nr__=1136) ! Number of records character*90 ar__ ! Full-size record integer*4 Seq (nr__) ! [1,1136]+ Sequential number of star in cluster character*1 n_Seq (nr__) ! [*] indicates a note in notes.dat file integer*4 ma (nr__) ! Brightness estimate from Plate A (2) character*1 n_ma (nr__) ! [*] when the star is not clearly seen on both * plates integer*4 dm (nr__) ! ? difference (ma-mb+1) between plates (2) real*8 Xpos (nr__) ! X mean position (in 5.sin(1arcmin)) real*8 Ypos (nr__) ! Y mean position (in 5.sin(1arcmin)) integer*4 dXpos (nr__) ! (30mas) ? Position difference on 2 plates (A-B) character*1 u_dXpos (nr__) ! [)] when star not found on the 2 reductions integer*4 dYpos (nr__) ! (30mas) ? Position difference on 2 plates (A-B) character*1 u_dYpos (nr__) ! [)] when star not found on the 2 reductions character*1 oRA_ (nr__) ! Offset in right ascension (sign) (1) integer*4 oRAm (nr__) ! (arcmin) Offset in right ascension (min) (1) real*4 oRAs (nr__) ! (arcsec) Offset in right ascension (s) (1) character*1 oDE_ (nr__) ! Offset in declination (sign) (1) integer*4 oDEm (nr__) ! (arcmin) Offset in declination (arcmin) (1) real*4 oDEs (nr__) ! (arcsec) Offset in declination (arcsec) (1) integer*4 S (nr__) ! ? Number of the star in Scheiner's (1892) integer*4 S2 (nr__) ! ? when 2 Scheiner's numbers represent the star character*1 u_S (nr__) ! [)] when Scheiner's star seen on 1 plate only character*1 n_S (nr__) ! [+] disagreement with Scheiner's, * in notes2.dat file real*4 Gr (nr__) ! (mag) ? Magnitude assigned in Scheiner's catalog real*4 dRA (nr__) ! (arcsec) ? Difference in RA Ludendorff - Scheiner character*1 u_dRA (nr__) ! [)] disagreement with Scheiner's real*4 dDE (nr__) ! (arcsec) ? Difference in DE Ludendorff - Scheiner character*1 u_dDE (nr__) ! [)] disagreement with Scheiner's *Note (1): The reference position is * RA = 16 38 06.14, DE = +36 39 19.06 (1900) *Note (2): the estimate consists in a number in the range 1 to 16; * this "Helligkeitsstuffe" (brightness step) is 1 for a photographic * magnitude of 8.5mag, 2 for 8.7, etc (i.e. mag = 8.3 + 0.2*ma). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'notes.dat' ! Individual notes (in German, english added) integer*4 nr__1 parameter (nr__1=706) ! Number of records character*80 ar__1 ! Full-size record integer*4 Seq_1 (nr__1) ! Sequential number (repeated when several lines) character*75 Note (nr__1) ! Text of the note (English translation added) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'notes2.dat' ! Individual remarks on Scheiner identifications, and Scheiner catalog and charts integer*4 nr__2 parameter (nr__2=234) ! Number of records character*80 ar__2 ! Full-size record integer*4 Seq_2 (nr__2) ! ? Sequential number integer*4 S_1 (nr__2) ! ? Scheiner number character*1 u_S_1 (nr__2) ! [)] when Scheiner's star seen on 1 plate only character*66 Note_1 (nr__2) ! Text of the note C============================================================================= C Loading file 'catalog.dat' ! Catalog of stars in Messier 13 C Format for file interpretation 1 format( + I4,A1,1X,I2,A1,1X,I2,2X,F7.4,1X,F7.4,2X,I3,A1,2X,I3,A1,1X,A1, + I1,1X,F5.2,1X,A1,I1,1X,F5.2,2X,I3,1X,I3,A1,A1,1X,F4.1,2X,F5.2, + A1,1X,F5.2,A1) C Effective file loading open(unit=1,status='old',file= +'catalog.dat') write(6,*) '....Loading file: catalog.dat' do i__=1,1136 read(1,'(A90)')ar__ read(ar__,1) + Seq(i__),n_Seq(i__),ma(i__),n_ma(i__),dm(i__),Xpos(i__), + Ypos(i__),dXpos(i__),u_dXpos(i__),dYpos(i__),u_dYpos(i__), + oRA_(i__),oRAm(i__),oRAs(i__),oDE_(i__),oDEm(i__),oDEs(i__), + S(i__),S2(i__),u_S(i__),n_S(i__),Gr(i__),dRA(i__),u_dRA(i__), + dDE(i__),u_dDE(i__) if(ar__(11:12) .EQ. '') dm(i__) = iNULL__ if(ar__(32:34) .EQ. '') dXpos(i__) = iNULL__ if(ar__(38:40) .EQ. '') dYpos(i__) = iNULL__ if(ar__(62:64) .EQ. '') S(i__) = iNULL__ if(ar__(66:68) .EQ. '') S2(i__) = iNULL__ if(ar__(72:75) .EQ. '') Gr(i__) = rNULL__ if(ar__(78:82) .EQ. '') dRA(i__) = rNULL__ if(ar__(85:89) .EQ. '') dDE(i__) = rNULL__ c ..............Just test output........... write(6,1) + Seq(i__),n_Seq(i__),ma(i__),n_ma(i__),dm(i__),Xpos(i__), + Ypos(i__),dXpos(i__),u_dXpos(i__),dYpos(i__),u_dYpos(i__), + oRA_(i__),oRAm(i__),oRAs(i__),oDE_(i__),oDEm(i__),oDEs(i__), + S(i__),S2(i__),u_S(i__),n_S(i__),Gr(i__),dRA(i__),u_dRA(i__), + dDE(i__),u_dDE(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'notes.dat' ! Individual notes (in German, english added) C Format for file interpretation 2 format(I4,1X,A75) C Effective file loading open(unit=1,status='old',file= +'notes.dat') write(6,*) '....Loading file: notes.dat' do i__=1,706 read(1,'(A80)')ar__1 read(ar__1,2)Seq_1(i__),Note(i__) c ..............Just test output........... write(6,2)Seq_1(i__),Note(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'notes2.dat' ! Individual remarks on Scheiner identifications, * and Scheiner catalog and charts C Format for file interpretation 3 format(I4,1X,I3,4X,A1,1X,A66) C Effective file loading open(unit=1,status='old',file= +'notes2.dat') write(6,*) '....Loading file: notes2.dat' do i__=1,234 read(1,'(A80)')ar__2 read(ar__2,3)Seq_2(i__),S_1(i__),u_S_1(i__),Note_1(i__) if(ar__2(1:4) .EQ. '') Seq_2(i__) = iNULL__ if(ar__2(6:8) .EQ. '') S_1(i__) = iNULL__ c ..............Just test output........... write(6,3)Seq_2(i__),S_1(i__),u_S_1(i__),Note_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end