Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/530/A93 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/A+A/530/A93 de Houtman, Kepler and Halley star catalogs (Verbunt+ 2011) *================================================================================ *Early star catalogues of the southern sky. * Verbunt, F., van Gent, R.H. * <Astron. Astrophys. 530, A93 (2011)> * =2011A&A...530A..93V C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'houtman.dat' ! Star Catalogue of de Houtman (1603) integer*4 nr__ parameter (nr__=304) ! Number of records character*70 ar__ ! Full-size record C 1601.0 position composed of: RAd RAm DE- DEd DEm real*8 RAdeg (nr__) ! (deg) Right Ascension 1601.0 real*8 DEdeg (nr__) ! (deg) Declination 1601.0 C ---------------------------------- ! (position vector(s) in degrees) integer*4 F (nr__) ! [1/304] Sequence number in de Houtman (1603) * Catalogue integer*4 C (nr__) ! [1/21] Sequence number of constellation character*3 cst (nr__) ! Abbreviation of constellation name integer*4 N (nr__) ! [1/56] Sequence number of star in constellation integer*4 RAd (nr__) ! (deg) [0/360] Degrees of right ascension (1) integer*4 RAm (nr__) ! (arcmin) Minutes of right ascension (1) character*1 DE_ (nr__) ! Sign of declination (1) integer*4 DEd (nr__) ! (deg) Degrees of declination (1) integer*4 DEm (nr__) ! (arcmin) Arcminutes of declination (1) integer*4 Vf (nr__) ! magnitude as given by de Houtman (9=nebulous) integer*4 HIP (nr__) ! Hipparcos number of identification integer*4 I (nr__) ! [1,6] Quality of identification (G3) integer*4 v_K (nr__) ! flag indicating comparison with Knobel (2) real*4 Vmag (nr__) ! Magnitude of Hipparcos identification real*4 D_RA (nr__) ! (arcmin) Difference in RA: Hipparcos-Houtman real*4 D_DE (nr__) ! (arcmin) Difference in DE: Hipparcos-Houtman real*4 Del (nr__) ! (arcmin) Offset between de Houtman star and * identification *Note (1): The right ascension is given by RA=RAd+RAm/60; * the declination, always negative, is given by DE=-(DEd+DEm/60). * The equinox is assumed to be AD 1601.0 *Note (2): Flag for comparison with Knobel (1917MNRAS..77..580K and * 1917MNRAS..77..414K) as follows: * 1 = same identification in Knobel * 2 = identification in Knobel to alternative identification * 3 = identification in Knobel different from our identification c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'classis.dat' ! Secunda Classis and Tertia Classis Kepler (1627) integer*4 nr__1 parameter (nr__1=434) ! Number of records character*82 ar__1 ! Full-size record C 1601 position composed of: Elon.Z Elon.d Elon.m Elat.- Elat.d Elat.m real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension 1601 real*8 DEdeg_1 (nr__1) ! (deg) Declination 1601 C ---------------------------------- ! (position vector(s) in degrees) integer*4 K (nr__1) ! [1005/1438] Sequence number in * Classis Catalogue (1) integer*4 C_1 (nr__1) ! [1/62] Sequence number of constellation (2) character*3 cst_1 (nr__1) ! Abbreviation of constellation name integer*4 N_1 (nr__1) ! [1/63] Sequence number of star in * constellation (3) integer*4 Elon_Z (nr__1) ! (30deg) [1,12] Zodiakal sign of ecliptic longitude, * Equinox AD1601 (G1) integer*4 Elon_d (nr__1) ! (deg) [0,30] Degrees of ecliptic longitude (G1) real*4 Elon_m (nr__1) ! (arcmin) Minutes of ecliptic longitude (G1) integer*4 Elat_d (nr__1) ! (deg) Degrees of ecliptic latitude (G2) integer*4 Elat_m (nr__1) ! (arcmin) Arcminutes of ecliptic latitude (G2) character*1 Elat__ (nr__1) ! [AB] Sign of ecliptic latitude, Australis (-) * or Borealis (+), Equinox AD1601 (G2) integer*4 VK (nr__1) ! [1,9] Magnitude as given by Kepler (4) integer*4 HIP_1 (nr__1) ! Hipparcos number of identification integer*4 I_1 (nr__1) ! [1,6] Quality of identification (G3) real*4 Vmag_1 (nr__1) ! Magnitude of Hipparcos identification real*8 Dlon (nr__1) ! (arcmin) Difference in longitude: Hipparcos-Kepler real*4 Dlat (nr__1) ! (arcmin) Difference in latitude: Hipparcos-Kepler real*4 Del_1 (nr__1) ! (arcmin) Offset between Kepler star and identification integer*4 P (nr__1) ! [0/1028] Corresponding sequence number in * Ptolemaios (5) integer*4 v_P (nr__1) ! Flag indicating comparison with Ptolemaios (5) *Note (1): * * K < 1005 refers to the Kepler (1627) edition of Brahe's catalogue * * 1004 < K < 1303 refers to Secunda Classis, the entries in which * we number starting with 1005 (ending with 1302) * * K > 1302 refers to Tertia Classis *Note (2): * * C < 47 constellation in Kepler (1627) edition of Brahe catalogue, * to which star is added in Secunda Classis * * 46 < C < 51 constellation from catalogue of Ptolemaios, not present * in Brahe Catalogue, for which Secunda Classis lists stars * * C > 50 new constellations from Plancius, in Tertia Classis * It sometimes refers to the Secunda Classis. For detail see paper. *Note (3): Numbering continued for stars in constellation listed in Brahe; * starting with 1 for stars in constellations not listed in Brahe *Note (4): 8 = no magnitude given in catalogue; 9 = nebulous *Note (5): Comparison with our identification for corresponding star in * Toomer's 1998 edition of the star catalogue of Ptolemaios (Princeton * Univ. Press " Ptolemy's Almagest", 2nd ed.) as follows: * 1 = same identification as for star in Classis * 2 = identification to alternative identification * 3 = identification different (including identification for new stars) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'aliter.dat' ! Variants on Secunda Classis (Kepler (1627) integer*4 nr__2 parameter (nr__2=20) ! Number of records character*82 ar__2 ! Full-size record C 1601 position composed of: Elon.Z Elon.d Elon.m Elat.- Elat.d Elat.m real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension 1601 real*8 DEdeg_2 (nr__2) ! (deg) Declination 1601 C ---------------------------------- ! (position vector(s) in degrees) integer*4 K_1 (nr__2) ! [1005/1438] Sequence number in * Classis Catalogue (1) integer*4 C_2 (nr__2) ! [1/62] Sequence number of constellation (2) character*3 cst_2 (nr__2) ! Abbreviation of constellation name integer*4 N_2 (nr__2) ! [1/63] Sequence number of star in * constellation (3) integer*4 Elon_Z_1 (nr__2) ! (30deg) [1,12] Zodiakal sign of ecliptic longitude, * Equinox AD1601 (G1) integer*4 Elon_d_1 (nr__2) ! (deg) [0,30] Degrees of ecliptic longitude (G1) real*4 Elon_m_1 (nr__2) ! (arcmin) Minutes of ecliptic longitude (G1) integer*4 Elat_d_1 (nr__2) ! (deg) Degrees of ecliptic latitude (G2) integer*4 Elat_m_1 (nr__2) ! (arcmin) Arcminutes of ecliptic latitude (G2) character*1 Elat___1 (nr__2) ! [AB] Sign of ecliptic latitude, Australis (-) * or Borealis (+), Equinox AD1601 (G2) integer*4 VK_1 (nr__2) ! [1,9] Magnitude as given by Kepler (4) integer*4 HIP_2 (nr__2) ! Hipparcos number of identification integer*4 I_2 (nr__2) ! [1,6] Quality of identification (G3) real*4 Vmag_2 (nr__2) ! Magnitude of Hipparcos identification real*8 Dlon_1 (nr__2) ! (arcmin) Difference in longitude: Hipparcos-Kepler real*4 Dlat_1 (nr__2) ! (arcmin) Difference in latitude: Hipparcos-Kepler real*4 Del_2 (nr__2) ! (arcmin) Offset between Kepler star and identification integer*4 P_1 (nr__2) ! [0/1028] Corresponding sequence number in * Ptolemaios (5) integer*4 v_P_1 (nr__2) ! Flag indicating comparison with Ptolemaios (5) *Note (1): * * K < 1005 refers to the Kepler (1627) edition of Brahe's catalogue * * 1004 < K < 1303 refers to Secunda Classis, the entries in which * we number starting with 1005 (ending with 1302) * * K > 1302 refers to Tertia Classis *Note (2): * * C < 47 constellation in Kepler (1627) edition of Brahe catalogue, * to which star is added in Secunda Classis * * 46 < C < 51 constellation from catalogue of Ptolemaios, not present * in Brahe Catalogue, for which Secunda Classis lists stars * * C > 50 new constellations from Plancius, in Tertia Classis * It sometimes refers to the Secunda Classis. For detail see paper. *Note (3): Numbering continued for stars in constellation listed in Brahe; * starting with 1 for stars in constellations not listed in Brahe *Note (4): 8 = no magnitude given in catalogue; 9 = nebulous *Note (5): Comparison with our identification for corresponding star in * Toomer's 1998 edition of the star catalogue of Ptolemaios (Princeton * Univ. Press " Ptolemy's Almagest", 2nd ed.) as follows: * 1 = same identification as for star in Classis * 2 = identification to alternative identification * 3 = identification different (including identification for new stars) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'halley.dat' ! Star catalogue of Halley (1679) integer*4 nr__3 parameter (nr__3=341) ! Number of records character*139 ar__3 ! Full-size record C 1678 position composed of: Elon.Z Elon.d Elon.m Elat.- Elat.d Elat.m real*8 RAdeg_3 (nr__3) ! (deg) Right Ascension 1678 real*8 DEdeg_3 (nr__3) ! (deg) Declination 1678 C ---------------------------------- ! (position vector(s) in degrees) integer*4 E (nr__3) ! Sequence number in Halley catalogue integer*4 C_3 (nr__3) ! Sequence number of constellation character*3 cst_3 (nr__3) ! Abbreviation of constellation name integer*4 N_3 (nr__3) ! [1/46] Sequence number of star in constellation integer*4 Elon_Z_2 (nr__3) ! (30deg) [1,12] Zodiakal sign of ecliptic longitude, * Equinox AD 1678 (G1) integer*4 Elon_d_2 (nr__3) ! (deg) [0,30] Degrees of ecliptic longitude (G1) real*4 Elon_m_2 (nr__3) ! (arcmin) Minutes of ecliptic longitude (G1) character*1 Elat___2 (nr__3) ! [-] Sign of ecliptic latitude (G2) integer*4 Elat_d_2 (nr__3) ! (deg) Degrees of ecliptic latitude, Equinox 1678 (G2) real*4 Elat_m_2 (nr__3) ! (arcmin) Arcminutes of ecliptic latitude (G2) real*4 VE (nr__3) ! Magnitude as given by Halley (9=nebulous) integer*4 HIP1 (nr__3) ! ? Hipparcos number of first standard star (2) integer*4 d1_d (nr__3) ! (deg) ? Degrees of angle to first standard star (2) integer*4 d1_m (nr__3) ! (arcmin) ? Arcminutes of angle to first standard star (2) integer*4 d1_s (nr__3) ! (arcsec) ? Arcseconds of angle to first standard star (2) integer*4 HIP2 (nr__3) ! ? Hipparcos number of second standard star (2) integer*4 d2_d (nr__3) ! (deg) ? Degrees of angle to second standard star (2) integer*4 d2_m (nr__3) ! (arcmin) ? Arcminutes of angle to 2nd standard star (2) integer*4 d2_s (nr__3) ! (arcsec) ? Arcseconds of angle to 2nd standard star (2) integer*4 n_HIP2 (nr__3) ! [-3,-1]? only one standard star used (2) integer*4 HIP_3 (nr__3) ! ?=0 Hipparcos number of identification integer*4 I_3 (nr__3) ! [1,6] Quality of identification (G3) real*4 Vmag_3 (nr__3) ! ? Magnitude of Hipparcos identification real*4 Dlon_2 (nr__3) ! (arcmin) ? Difference in longitude: Hipparcos-Kepler real*4 Dlat_2 (nr__3) ! (arcmin) ? Difference in latitude: Hipparcos-Kepler real*4 Del_3 (nr__3) ! (arcmin) ? Offset between Kepler star and identification integer*4 HIP3 (nr__3) ! ? Hipparcos number of third standard star (3) integer*4 d3_d (nr__3) ! (deg) ? Degrees of angle to third standard star (3) integer*4 d3_m (nr__3) ! (arcmin) ? Arcminutes of angle to third standard star (3) integer*4 d3_s (nr__3) ! (arcsec) ? Arcseconds of angle to third standard star (3) integer*4 HIP4 (nr__3) ! ? Hipparcos number of fourth standard star (3) integer*4 d4_d (nr__3) ! (deg) ? Degrees of angle to fourth standard star (3) integer*4 d4_m (nr__3) ! (arcmin) ? Arcminutes of angle to fourth standard star (3) integer*4 d4_s (nr__3) ! (arcsec) ? Arcseconds of angle to fourth standard star (3) *Note (2): Halley give measured angles to two standard stars; for some stars * more than two standard stars are given, sometimes only one standard * star is used, and combined with information from Brahe or Kepler. We * indicates these with following n_HIP2 values: * -1 = longitude from Brahe * -2 = latitude from Brahe * -3 = longitude from Kepler * These values are announced in the paper as "negative value for HIP" *Note (3): These parameters give extra distances to standard stars * for four sequence numbers, announced in the paper as * "at the end of the table, four sequence numbers are repeated". C============================================================================= C Loading file 'houtman.dat' ! Star Catalogue of de Houtman (1603) C Format for file interpretation 1 format( + I3,1X,I2,2X,A3,1X,I2,2X,I3,1X,I2,1X,A1,I2,1X,I2,2X,I1,1X,I6, + 2X,I1,1X,I1,1X,F4.1,1X,F6.1,1X,F6.1,1X,F6.1) C Effective file loading open(unit=1,status='old',file= +'houtman.dat') write(6,*) '....Loading file: houtman.dat' do i__=1,304 read(1,'(A70)')ar__ read(ar__,1) + F(i__),C(i__),cst(i__),N(i__),RAd(i__),RAm(i__),DE_(i__), + DEd(i__),DEm(i__),Vf(i__),HIP(i__),I(i__),v_K(i__),Vmag(i__), + D_RA(i__),D_DE(i__),Del(i__) 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(RAd(i__) .GT. -180) RAdeg(i__)=RAd(i__) if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/60. if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__) if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,1) + F(i__),C(i__),cst(i__),N(i__),RAd(i__),RAm(i__),DE_(i__), + DEd(i__),DEm(i__),Vf(i__),HIP(i__),I(i__),v_K(i__),Vmag(i__), + D_RA(i__),D_DE(i__),Del(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'classis.dat' ! Secunda Classis and Tertia Classis Kepler (1627) C Format for file interpretation 2 format( + I4,2X,I2,2X,A3,1X,I2,1X,I2,1X,I2,1X,F4.1,2X,I2,1X,I2,1X,A1,1X, + I1,1X,I6,1X,I1,1X,F4.1,F7.1,1X,F6.1,1X,F6.1,3X,I4,1X,I1) C Effective file loading open(unit=1,status='old',file= +'classis.dat') write(6,*) '....Loading file: classis.dat' do i__=1,434 read(1,'(A82)')ar__1 read(ar__1,2) + K(i__),C_1(i__),cst_1(i__),N_1(i__),Elon_Z(i__),Elon_d(i__), + Elon_m(i__),Elat_d(i__),Elat_m(i__),Elat__(i__),VK(i__), + HIP_1(i__),I_1(i__),Vmag_1(i__),Dlon(i__),Dlat(i__), + Del_1(i__),P(i__),v_P(i__) RAdeg_1(i__) = rNULL__ DEdeg_1(i__) = rNULL__ c Derive coordinates RAdeg_1 and DEdeg_1 from input data c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown) c ..............Just test output........... write(6,2) + K(i__),C_1(i__),cst_1(i__),N_1(i__),Elon_Z(i__),Elon_d(i__), + Elon_m(i__),Elat_d(i__),Elat_m(i__),Elat__(i__),VK(i__), + HIP_1(i__),I_1(i__),Vmag_1(i__),Dlon(i__),Dlat(i__), + Del_1(i__),P(i__),v_P(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'aliter.dat' ! Variants on Secunda Classis (Kepler (1627) C Format for file interpretation 3 format( + I4,2X,I2,2X,A3,1X,I2,1X,I2,1X,I2,1X,F4.1,2X,I2,1X,I2,1X,A1,1X, + I1,1X,I6,1X,I1,1X,F4.1,F7.1,1X,F6.1,1X,F6.1,3X,I4,1X,I1) C Effective file loading open(unit=1,status='old',file= +'aliter.dat') write(6,*) '....Loading file: aliter.dat' do i__=1,20 read(1,'(A82)')ar__2 read(ar__2,3) + K_1(i__),C_2(i__),cst_2(i__),N_2(i__),Elon_Z_1(i__), + Elon_d_1(i__),Elon_m_1(i__),Elat_d_1(i__),Elat_m_1(i__), + Elat___1(i__),VK_1(i__),HIP_2(i__),I_2(i__),Vmag_2(i__), + Dlon_1(i__),Dlat_1(i__),Del_2(i__),P_1(i__),v_P_1(i__) RAdeg_2(i__) = rNULL__ DEdeg_2(i__) = rNULL__ c Derive coordinates RAdeg_2 and DEdeg_2 from input data c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown) c ..............Just test output........... write(6,3) + K_1(i__),C_2(i__),cst_2(i__),N_2(i__),Elon_Z_1(i__), + Elon_d_1(i__),Elon_m_1(i__),Elat_d_1(i__),Elat_m_1(i__), + Elat___1(i__),VK_1(i__),HIP_2(i__),I_2(i__),Vmag_2(i__), + Dlon_1(i__),Dlat_1(i__),Del_2(i__),P_1(i__),v_P_1(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'halley.dat' ! Star catalogue of Halley (1679) C Format for file interpretation 4 format( + I3,1X,I2,2X,A3,1X,I2,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,F4.1,1X, + F3.1,1X,I6,1X,I2,1X,I2,1X,I2,1X,I6,1X,I2,1X,I2,1X,I2,1X,I2,1X, + I6,1X,I1,1X,F4.1,1X,F6.1,1X,F6.1,1X,F6.1,1X,I6,1X,I2,1X,I2,1X, + I2,1X,I5,1X,I2,1X,I2,1X,I2) C Effective file loading open(unit=1,status='old',file= +'halley.dat') write(6,*) '....Loading file: halley.dat' do i__=1,341 read(1,'(A139)')ar__3 read(ar__3,4) + E(i__),C_3(i__),cst_3(i__),N_3(i__),Elon_Z_2(i__), + Elon_d_2(i__),Elon_m_2(i__),Elat___2(i__),Elat_d_2(i__), + Elat_m_2(i__),VE(i__),HIP1(i__),d1_d(i__),d1_m(i__),d1_s(i__), + HIP2(i__),d2_d(i__),d2_m(i__),d2_s(i__),n_HIP2(i__), + HIP_3(i__),I_3(i__),Vmag_3(i__),Dlon_2(i__),Dlat_2(i__), + Del_3(i__),HIP3(i__),d3_d(i__),d3_m(i__),d3_s(i__),HIP4(i__), + d4_d(i__),d4_m(i__),d4_s(i__) if(ar__3(40:45) .EQ. '') HIP1(i__) = iNULL__ if(ar__3(47:48) .EQ. '') d1_d(i__) = iNULL__ if(ar__3(50:51) .EQ. '') d1_m(i__) = iNULL__ if(ar__3(53:54) .EQ. '') d1_s(i__) = iNULL__ if(ar__3(56:61) .EQ. '') HIP2(i__) = iNULL__ if(ar__3(63:64) .EQ. '') d2_d(i__) = iNULL__ if(ar__3(66:67) .EQ. '') d2_m(i__) = iNULL__ if(ar__3(69:70) .EQ. '') d2_s(i__) = iNULL__ if(ar__3(72:73) .EQ. '') n_HIP2(i__) = iNULL__ if(ar__3(84:87) .EQ. '') Vmag_3(i__) = rNULL__ if(ar__3(89:94) .EQ. '') Dlon_2(i__) = rNULL__ if(ar__3(96:101) .EQ. '') Dlat_2(i__) = rNULL__ if(ar__3(103:108) .EQ. '') Del_3(i__) = rNULL__ if(ar__3(110:115) .EQ. '') HIP3(i__) = iNULL__ if(ar__3(117:118) .EQ. '') d3_d(i__) = iNULL__ if(ar__3(120:121) .EQ. '') d3_m(i__) = iNULL__ if(ar__3(123:124) .EQ. '') d3_s(i__) = iNULL__ if(ar__3(126:130) .EQ. '') HIP4(i__) = iNULL__ if(ar__3(132:133) .EQ. '') d4_d(i__) = iNULL__ if(ar__3(135:136) .EQ. '') d4_m(i__) = iNULL__ if(ar__3(138:139) .EQ. '') d4_s(i__) = iNULL__ RAdeg_3(i__) = rNULL__ DEdeg_3(i__) = rNULL__ c Derive coordinates RAdeg_3 and DEdeg_3 from input data c (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown) c ..............Just test output........... write(6,4) + E(i__),C_3(i__),cst_3(i__),N_3(i__),Elon_Z_2(i__), + Elon_d_2(i__),Elon_m_2(i__),Elat___2(i__),Elat_d_2(i__), + Elat_m_2(i__),VE(i__),HIP1(i__),d1_d(i__),d1_m(i__),d1_s(i__), + HIP2(i__),d2_d(i__),d2_m(i__),d2_s(i__),n_HIP2(i__), + HIP_3(i__),I_3(i__),Vmag_3(i__),Dlon_2(i__),Dlat_2(i__), + Del_3(i__),HIP3(i__),d3_d(i__),d3_m(i__),d3_s(i__),HIP4(i__), + d4_d(i__),d4_m(i__),d4_s(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end