Conversion of standardized ReadMe file for
file /./ftp/cats/I/346 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-29 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. I/346 FON Astrographic Catalogue Southern Part (FONAC-S) (Yuldoshev+, 2017) *================================================================================ *Catalog of Equatorial Coordinates and B- Magnitudes of Stars of the Kitab Part *of the FON Project. * Yuldoshev Q.X., Ehgamberdiev Sh.A., Muminov M.M., Protsyuk Yu.I., Relke H., * Andruk V.M. * <Kinem. and Physics of Celest. Bodies. 33, 250 (2017)> * =2017KPCB...33..250Y * =2017KFNT...33..250Y * =2018yCat.1346....0Y C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'fs00.dat' ! Star data file (zone from +02.0 to -01.9 degrees) integer*4 nr__ parameter (nr__=2550586) ! Number of records character*164 ar__ ! Full-size record C J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs real*8 RAdeg (nr__) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID (nr__) ! Serial number of the star integer*4 k (nr__) ! Number of determinations (plates) real*8 RAhour (nr__) ! (h) Right Ascension (J2000) real*4 e_RAhour (nr__) ! (arcsec) Error in right ascension, RA*cos(DE) real*8 DEdeg_1 (nr__) ! (deg) Declination (J2000) real*4 e_DEdeg (nr__) ! (arcsec) Error in Declination real*4 Bmag (nr__) ! (mag) Johnson Blue magnitude real*4 e_Bmag (nr__) ! (mag) Error in B magnitude real*8 Epoch (nr__) ! (yr) Epoch of observations real*4 FWHM (nr__) ! (1.46pix) Diameter of star image real*8 cI (nr__) ! Intensity in the center of image * (in max256 unit) integer*4 RAh (nr__) ! (h) Right ascension (J2000) integer*4 RAm (nr__) ! (min) Right ascension (J2000) real*4 RAs (nr__) ! (s) [0/60] Right ascension (J2000) character*1 DE_ (nr__) ! Declination sign (J2000) integer*4 DEd (nr__) ! (deg) Declination (J2000) integer*4 DEm (nr__) ! (arcmin) Declination (J2000) real*4 DEs (nr__) ! (arcsec) [0/60] Declination (J2000) real*4 f_mag (nr__) ! (mag) ?=0 UCAC4 f magnitude (1) real*8 pmRA (nr__) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1) real*8 pmDE (nr__) ! (mas/yr) ?=0 Proper motion in DE (1) real*4 B4mag (nr__) ! (mag) ?=0 UCAC4 B magnitude (1) real*4 Vmag (nr__) ! (mag) ?=0 UCAC4 V magnitude (1) real*4 rmag (nr__) ! (mag) ?=0 UCAC4 r magnitude (1) real*4 imag (nr__) ! (mag) ?=0 UCAC4 i magnitude (1) *Note (1): Proper motions and magnitudes were taken from UCAC4 catalogue. * Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z, * Cat. I/322). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'fs04.dat' ! Star data file (zone from -02.0 to -05.9 degrees) integer*4 nr__1 parameter (nr__1=2233270) ! Number of records character*164 ar__1 ! Full-size record C J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000 real*8 DEdeg_2 (nr__1) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID_1 (nr__1) ! Serial number of the star integer*4 k_1 (nr__1) ! Number of determinations (plates) real*8 RAhour_1 (nr__1) ! (h) Right Ascension (J2000) real*4 e_RAhour_1 (nr__1) ! (arcsec) Error in right ascension, RA*cos(DE) real*8 DEdeg_3 (nr__1) ! (deg) Declination (J2000) real*4 e_DEdeg_1 (nr__1) ! (arcsec) Error in Declination real*4 Bmag_1 (nr__1) ! (mag) Johnson Blue magnitude real*4 e_Bmag_1 (nr__1) ! (mag) Error in B magnitude real*8 Epoch_1 (nr__1) ! (yr) Epoch of observations real*4 FWHM_1 (nr__1) ! (1.46pix) Diameter of star image real*8 cI_1 (nr__1) ! Intensity in the center of image * (in max256 unit) integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000) integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000) real*4 RAs_1 (nr__1) ! (s) [0/60] Right ascension (J2000) character*1 DE__1 (nr__1) ! Declination sign (J2000) integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000) integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000) real*4 DEs_1 (nr__1) ! (arcsec) [0/60] Declination (J2000) real*4 f_mag_1 (nr__1) ! (mag) ?=0 UCAC4 f magnitude (1) real*8 pmRA_1 (nr__1) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1) real*8 pmDE_1 (nr__1) ! (mas/yr) ?=0 Proper motion in DE (1) real*4 B4mag_1 (nr__1) ! (mag) ?=0 UCAC4 B magnitude (1) real*4 Vmag_1 (nr__1) ! (mag) ?=0 UCAC4 V magnitude (1) real*4 rmag_1 (nr__1) ! (mag) ?=0 UCAC4 r magnitude (1) real*4 imag_1 (nr__1) ! (mag) ?=0 UCAC4 i magnitude (1) *Note (1): Proper motions and magnitudes were taken from UCAC4 catalogue. * Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z, * Cat. I/322). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'fs08.dat' ! Star data file (zone from -06.0 to -09.9 degrees) integer*4 nr__2 parameter (nr__2=2486306) ! Number of records character*164 ar__2 ! Full-size record C J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension J2000 real*8 DEdeg_4 (nr__2) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID_2 (nr__2) ! Serial number of the star integer*4 k_2 (nr__2) ! Number of determinations (plates) real*8 RAhour_2 (nr__2) ! (h) Right Ascension (J2000) real*4 e_RAhour_2 (nr__2) ! (arcsec) Error in right ascension, RA*cos(DE) real*8 DEdeg_5 (nr__2) ! (deg) Declination (J2000) real*4 e_DEdeg_2 (nr__2) ! (arcsec) Error in Declination real*4 Bmag_2 (nr__2) ! (mag) Johnson Blue magnitude real*4 e_Bmag_2 (nr__2) ! (mag) Error in B magnitude real*8 Epoch_2 (nr__2) ! (yr) Epoch of observations real*4 FWHM_2 (nr__2) ! (1.46pix) Diameter of star image real*8 cI_2 (nr__2) ! Intensity in the center of image * (in max256 unit) integer*4 RAh_2 (nr__2) ! (h) Right ascension (J2000) integer*4 RAm_2 (nr__2) ! (min) Right ascension (J2000) real*4 RAs_2 (nr__2) ! (s) [0/60] Right ascension (J2000) character*1 DE__2 (nr__2) ! Declination sign (J2000) integer*4 DEd_2 (nr__2) ! (deg) Declination (J2000) integer*4 DEm_2 (nr__2) ! (arcmin) Declination (J2000) real*4 DEs_2 (nr__2) ! (arcsec) [0/60] Declination (J2000) real*4 f_mag_2 (nr__2) ! (mag) ?=0 UCAC4 f magnitude (1) real*8 pmRA_2 (nr__2) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1) real*8 pmDE_2 (nr__2) ! (mas/yr) ?=0 Proper motion in DE (1) real*4 B4mag_2 (nr__2) ! (mag) ?=0 UCAC4 B magnitude (1) real*4 Vmag_2 (nr__2) ! (mag) ?=0 UCAC4 V magnitude (1) real*4 rmag_2 (nr__2) ! (mag) ?=0 UCAC4 r magnitude (1) real*4 imag_2 (nr__2) ! (mag) ?=0 UCAC4 i magnitude (1) *Note (1): Proper motions and magnitudes were taken from UCAC4 catalogue. * Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z, * Cat. I/322). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'fs12.dat' ! Star data file (zone from -10.0 to -13.9 degrees) integer*4 nr__3 parameter (nr__3=2416612) ! Number of records character*164 ar__3 ! Full-size record C J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs real*8 RAdeg_3 (nr__3) ! (deg) Right Ascension J2000 real*8 DEdeg_6 (nr__3) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID_3 (nr__3) ! Serial number of the star integer*4 k_3 (nr__3) ! Number of determinations (plates) real*8 RAhour_3 (nr__3) ! (h) Right Ascension (J2000) real*4 e_RAhour_3 (nr__3) ! (arcsec) Error in right ascension, RA*cos(DE) real*8 DEdeg_7 (nr__3) ! (deg) Declination (J2000) real*4 e_DEdeg_3 (nr__3) ! (arcsec) Error in Declination real*4 Bmag_3 (nr__3) ! (mag) Johnson Blue magnitude real*4 e_Bmag_3 (nr__3) ! (mag) Error in B magnitude real*8 Epoch_3 (nr__3) ! (yr) Epoch of observations real*4 FWHM_3 (nr__3) ! (1.46pix) Diameter of star image real*8 cI_3 (nr__3) ! Intensity in the center of image * (in max256 unit) integer*4 RAh_3 (nr__3) ! (h) Right ascension (J2000) integer*4 RAm_3 (nr__3) ! (min) Right ascension (J2000) real*4 RAs_3 (nr__3) ! (s) [0/60] Right ascension (J2000) character*1 DE__3 (nr__3) ! Declination sign (J2000) integer*4 DEd_3 (nr__3) ! (deg) Declination (J2000) integer*4 DEm_3 (nr__3) ! (arcmin) Declination (J2000) real*4 DEs_3 (nr__3) ! (arcsec) [0/60] Declination (J2000) real*4 f_mag_3 (nr__3) ! (mag) ?=0 UCAC4 f magnitude (1) real*8 pmRA_3 (nr__3) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1) real*8 pmDE_3 (nr__3) ! (mas/yr) ?=0 Proper motion in DE (1) real*4 B4mag_3 (nr__3) ! (mag) ?=0 UCAC4 B magnitude (1) real*4 Vmag_3 (nr__3) ! (mag) ?=0 UCAC4 V magnitude (1) real*4 rmag_3 (nr__3) ! (mag) ?=0 UCAC4 r magnitude (1) real*4 imag_3 (nr__3) ! (mag) ?=0 UCAC4 i magnitude (1) *Note (1): Proper motions and magnitudes were taken from UCAC4 catalogue. * Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z, * Cat. I/322). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'fs16.dat' ! Star data file (zone from -14.0 to -17.9 degrees) integer*4 nr__4 parameter (nr__4=2575447) ! Number of records character*164 ar__4 ! Full-size record C J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs real*8 RAdeg_4 (nr__4) ! (deg) Right Ascension J2000 real*8 DEdeg_8 (nr__4) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID_4 (nr__4) ! Serial number of the star integer*4 k_4 (nr__4) ! Number of determinations (plates) real*8 RAhour_4 (nr__4) ! (h) Right Ascension (J2000) real*4 e_RAhour_4 (nr__4) ! (arcsec) Error in right ascension, RA*cos(DE) real*8 DEdeg_9 (nr__4) ! (deg) Declination (J2000) real*4 e_DEdeg_4 (nr__4) ! (arcsec) Error in Declination real*4 Bmag_4 (nr__4) ! (mag) Johnson Blue magnitude real*4 e_Bmag_4 (nr__4) ! (mag) Error in B magnitude real*8 Epoch_4 (nr__4) ! (yr) Epoch of observations real*4 FWHM_4 (nr__4) ! (1.46pix) Diameter of star image real*8 cI_4 (nr__4) ! Intensity in the center of image * (in max256 unit) integer*4 RAh_4 (nr__4) ! (h) Right ascension (J2000) integer*4 RAm_4 (nr__4) ! (min) Right ascension (J2000) real*4 RAs_4 (nr__4) ! (s) [0/60] Right ascension (J2000) character*1 DE__4 (nr__4) ! Declination sign (J2000) integer*4 DEd_4 (nr__4) ! (deg) Declination (J2000) integer*4 DEm_4 (nr__4) ! (arcmin) Declination (J2000) real*4 DEs_4 (nr__4) ! (arcsec) [0/60] Declination (J2000) real*4 f_mag_4 (nr__4) ! (mag) ?=0 UCAC4 f magnitude (1) real*8 pmRA_4 (nr__4) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1) real*8 pmDE_4 (nr__4) ! (mas/yr) ?=0 Proper motion in DE (1) real*4 B4mag_4 (nr__4) ! (mag) ?=0 UCAC4 B magnitude (1) real*4 Vmag_4 (nr__4) ! (mag) ?=0 UCAC4 V magnitude (1) real*4 rmag_4 (nr__4) ! (mag) ?=0 UCAC4 r magnitude (1) real*4 imag_4 (nr__4) ! (mag) ?=0 UCAC4 i magnitude (1) *Note (1): Proper motions and magnitudes were taken from UCAC4 catalogue. * Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z, * Cat. I/322). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'fs20.dat' ! Star data file (zone from -18.0 to -20.9 degrees) integer*4 nr__5 parameter (nr__5=1151047) ! Number of records character*164 ar__5 ! Full-size record C J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs real*8 RAdeg_5 (nr__5) ! (deg) Right Ascension J2000 real*8 DEdeg_10 (nr__5) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID_5 (nr__5) ! Serial number of the star integer*4 k_5 (nr__5) ! Number of determinations (plates) real*8 RAhour_5 (nr__5) ! (h) Right Ascension (J2000) real*4 e_RAhour_5 (nr__5) ! (arcsec) Error in right ascension, RA*cos(DE) real*8 DEdeg_11 (nr__5) ! (deg) Declination (J2000) real*4 e_DEdeg_5 (nr__5) ! (arcsec) Error in Declination real*4 Bmag_5 (nr__5) ! (mag) Johnson Blue magnitude real*4 e_Bmag_5 (nr__5) ! (mag) Error in B magnitude real*8 Epoch_5 (nr__5) ! (yr) Epoch of observations real*4 FWHM_5 (nr__5) ! (1.46pix) Diameter of star image real*8 cI_5 (nr__5) ! Intensity in the center of image * (in max256 unit) integer*4 RAh_5 (nr__5) ! (h) Right ascension (J2000) integer*4 RAm_5 (nr__5) ! (min) Right ascension (J2000) real*4 RAs_5 (nr__5) ! (s) [0/60] Right ascension (J2000) character*1 DE__5 (nr__5) ! Declination sign (J2000) integer*4 DEd_5 (nr__5) ! (deg) Declination (J2000) integer*4 DEm_5 (nr__5) ! (arcmin) Declination (J2000) real*4 DEs_5 (nr__5) ! (arcsec) [0/60] Declination (J2000) real*4 f_mag_5 (nr__5) ! (mag) ?=0 UCAC4 f magnitude (1) real*8 pmRA_5 (nr__5) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1) real*8 pmDE_5 (nr__5) ! (mas/yr) ?=0 Proper motion in DE (1) real*4 B4mag_5 (nr__5) ! (mag) ?=0 UCAC4 B magnitude (1) real*4 Vmag_5 (nr__5) ! (mag) ?=0 UCAC4 V magnitude (1) real*4 rmag_5 (nr__5) ! (mag) ?=0 UCAC4 r magnitude (1) real*4 imag_5 (nr__5) ! (mag) ?=0 UCAC4 i magnitude (1) *Note (1): Proper motions and magnitudes were taken from UCAC4 catalogue. * Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z, * Cat. I/322). C============================================================================= C Loading file 'fs00.dat' ! Star data file (zone from +02.0 to -01.9 degrees) C Format for file interpretation 1 format( + I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X, + F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3, + 2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'fs00.dat') write(6,*) '....Loading file: fs00.dat' do i__=1,2550586 read(1,'(A164)')ar__ read(ar__,1) + ID(i__),k(i__),RAhour(i__),e_RAhour(i__),DEdeg_1(i__), + e_DEdeg(i__),Bmag(i__),e_Bmag(i__),Epoch(i__),FWHM(i__), + cI(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__), + DEs(i__),f_mag(i__),pmRA(i__),pmDE(i__),B4mag(i__),Vmag(i__), + rmag(i__),imag(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(RAhour(i__) .GT. -180) RAdeg(i__)=RAhour(i__)*15. if(DEdeg_1(i__) .GE. 0) DEdeg(i__)=DEdeg_1(i__) if(DEd(i__) .GE. 0) DEdeg(i__)=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) + ID(i__),k(i__),RAhour(i__),e_RAhour(i__),DEdeg_1(i__), + e_DEdeg(i__),Bmag(i__),e_Bmag(i__),Epoch(i__),FWHM(i__), + cI(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__), + DEs(i__),f_mag(i__),pmRA(i__),pmDE(i__),B4mag(i__),Vmag(i__), + rmag(i__),imag(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 'fs04.dat' ! Star data file (zone from -02.0 to -05.9 degrees) C Format for file interpretation 2 format( + I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X, + F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3, + 2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'fs04.dat') write(6,*) '....Loading file: fs04.dat' do i__=1,2233270 read(1,'(A164)')ar__1 read(ar__1,2) + ID_1(i__),k_1(i__),RAhour_1(i__),e_RAhour_1(i__),DEdeg_3(i__), + e_DEdeg_1(i__),Bmag_1(i__),e_Bmag_1(i__),Epoch_1(i__), + FWHM_1(i__),cI_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__), + DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),f_mag_1(i__), + pmRA_1(i__),pmDE_1(i__),B4mag_1(i__),Vmag_1(i__),rmag_1(i__), + imag_1(i__) RAdeg_1(i__) = rNULL__ DEdeg_2(i__) = rNULL__ c Derive coordinates RAdeg_1 and DEdeg_2 from input data c (RAdeg_1 and DEdeg_2 are set to rNULL__ when unknown) if(RAhour_1(i__) .GT. -180) RAdeg_1(i__)=RAhour_1(i__)*15. if(DEdeg_3(i__) .GE. 0) DEdeg_2(i__)=DEdeg_3(i__) if(DEd_1(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEd_1(i__) if(DEm_1(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_1(i__)/60. if(DEs_1(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_1(i__)/3600. if(DE__1(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__) c ..............Just test output........... write(6,2) + ID_1(i__),k_1(i__),RAhour_1(i__),e_RAhour_1(i__),DEdeg_3(i__), + e_DEdeg_1(i__),Bmag_1(i__),e_Bmag_1(i__),Epoch_1(i__), + FWHM_1(i__),cI_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__), + DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),f_mag_1(i__), + pmRA_1(i__),pmDE_1(i__),B4mag_1(i__),Vmag_1(i__),rmag_1(i__), + imag_1(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'fs08.dat' ! Star data file (zone from -06.0 to -09.9 degrees) C Format for file interpretation 3 format( + I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X, + F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3, + 2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'fs08.dat') write(6,*) '....Loading file: fs08.dat' do i__=1,2486306 read(1,'(A164)')ar__2 read(ar__2,3) + ID_2(i__),k_2(i__),RAhour_2(i__),e_RAhour_2(i__),DEdeg_5(i__), + e_DEdeg_2(i__),Bmag_2(i__),e_Bmag_2(i__),Epoch_2(i__), + FWHM_2(i__),cI_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__), + DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),f_mag_2(i__), + pmRA_2(i__),pmDE_2(i__),B4mag_2(i__),Vmag_2(i__),rmag_2(i__), + imag_2(i__) RAdeg_2(i__) = rNULL__ DEdeg_4(i__) = rNULL__ c Derive coordinates RAdeg_2 and DEdeg_4 from input data c (RAdeg_2 and DEdeg_4 are set to rNULL__ when unknown) if(RAhour_2(i__) .GT. -180) RAdeg_2(i__)=RAhour_2(i__)*15. if(DEdeg_5(i__) .GE. 0) DEdeg_4(i__)=DEdeg_5(i__) if(DEd_2(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEd_2(i__) if(DEm_2(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEm_2(i__)/60. if(DEs_2(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEs_2(i__)/3600. if(DE__2(i__).EQ.'-'.AND.DEdeg_4(i__).GE.0) DEdeg_4(i__)=-DEdeg_4(i__) c ..............Just test output........... write(6,3) + ID_2(i__),k_2(i__),RAhour_2(i__),e_RAhour_2(i__),DEdeg_5(i__), + e_DEdeg_2(i__),Bmag_2(i__),e_Bmag_2(i__),Epoch_2(i__), + FWHM_2(i__),cI_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__), + DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),f_mag_2(i__), + pmRA_2(i__),pmDE_2(i__),B4mag_2(i__),Vmag_2(i__),rmag_2(i__), + imag_2(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'fs12.dat' ! Star data file (zone from -10.0 to -13.9 degrees) C Format for file interpretation 4 format( + I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X, + F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3, + 2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'fs12.dat') write(6,*) '....Loading file: fs12.dat' do i__=1,2416612 read(1,'(A164)')ar__3 read(ar__3,4) + ID_3(i__),k_3(i__),RAhour_3(i__),e_RAhour_3(i__),DEdeg_7(i__), + e_DEdeg_3(i__),Bmag_3(i__),e_Bmag_3(i__),Epoch_3(i__), + FWHM_3(i__),cI_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__), + DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),f_mag_3(i__), + pmRA_3(i__),pmDE_3(i__),B4mag_3(i__),Vmag_3(i__),rmag_3(i__), + imag_3(i__) RAdeg_3(i__) = rNULL__ DEdeg_6(i__) = rNULL__ c Derive coordinates RAdeg_3 and DEdeg_6 from input data c (RAdeg_3 and DEdeg_6 are set to rNULL__ when unknown) if(RAhour_3(i__) .GT. -180) RAdeg_3(i__)=RAhour_3(i__)*15. if(DEdeg_7(i__) .GE. 0) DEdeg_6(i__)=DEdeg_7(i__) if(DEd_3(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEd_3(i__) if(DEm_3(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEm_3(i__)/60. if(DEs_3(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEs_3(i__)/3600. if(DE__3(i__).EQ.'-'.AND.DEdeg_6(i__).GE.0) DEdeg_6(i__)=-DEdeg_6(i__) c ..............Just test output........... write(6,4) + ID_3(i__),k_3(i__),RAhour_3(i__),e_RAhour_3(i__),DEdeg_7(i__), + e_DEdeg_3(i__),Bmag_3(i__),e_Bmag_3(i__),Epoch_3(i__), + FWHM_3(i__),cI_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__), + DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),f_mag_3(i__), + pmRA_3(i__),pmDE_3(i__),B4mag_3(i__),Vmag_3(i__),rmag_3(i__), + imag_3(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_6(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'fs16.dat' ! Star data file (zone from -14.0 to -17.9 degrees) C Format for file interpretation 5 format( + I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X, + F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3, + 2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'fs16.dat') write(6,*) '....Loading file: fs16.dat' do i__=1,2575447 read(1,'(A164)')ar__4 read(ar__4,5) + ID_4(i__),k_4(i__),RAhour_4(i__),e_RAhour_4(i__),DEdeg_9(i__), + e_DEdeg_4(i__),Bmag_4(i__),e_Bmag_4(i__),Epoch_4(i__), + FWHM_4(i__),cI_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__), + DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),f_mag_4(i__), + pmRA_4(i__),pmDE_4(i__),B4mag_4(i__),Vmag_4(i__),rmag_4(i__), + imag_4(i__) RAdeg_4(i__) = rNULL__ DEdeg_8(i__) = rNULL__ c Derive coordinates RAdeg_4 and DEdeg_8 from input data c (RAdeg_4 and DEdeg_8 are set to rNULL__ when unknown) if(RAhour_4(i__) .GT. -180) RAdeg_4(i__)=RAhour_4(i__)*15. if(DEdeg_9(i__) .GE. 0) DEdeg_8(i__)=DEdeg_9(i__) if(DEd_4(i__) .GE. 0) DEdeg_8(i__)=DEdeg_8(i__)+DEd_4(i__) if(DEm_4(i__) .GE. 0) DEdeg_8(i__)=DEdeg_8(i__)+DEm_4(i__)/60. if(DEs_4(i__) .GE. 0) DEdeg_8(i__)=DEdeg_8(i__)+DEs_4(i__)/3600. if(DE__4(i__).EQ.'-'.AND.DEdeg_8(i__).GE.0) DEdeg_8(i__)=-DEdeg_8(i__) c ..............Just test output........... write(6,5) + ID_4(i__),k_4(i__),RAhour_4(i__),e_RAhour_4(i__),DEdeg_9(i__), + e_DEdeg_4(i__),Bmag_4(i__),e_Bmag_4(i__),Epoch_4(i__), + FWHM_4(i__),cI_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__), + DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),f_mag_4(i__), + pmRA_4(i__),pmDE_4(i__),B4mag_4(i__),Vmag_4(i__),rmag_4(i__), + imag_4(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_4(i__),DEdeg_8(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'fs20.dat' ! Star data file (zone from -18.0 to -20.9 degrees) C Format for file interpretation 6 format( + I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X, + F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3, + 2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'fs20.dat') write(6,*) '....Loading file: fs20.dat' do i__=1,1151047 read(1,'(A164)')ar__5 read(ar__5,6) + ID_5(i__),k_5(i__),RAhour_5(i__),e_RAhour_5(i__), + DEdeg_11(i__),e_DEdeg_5(i__),Bmag_5(i__),e_Bmag_5(i__), + Epoch_5(i__),FWHM_5(i__),cI_5(i__),RAh_5(i__),RAm_5(i__), + RAs_5(i__),DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__), + f_mag_5(i__),pmRA_5(i__),pmDE_5(i__),B4mag_5(i__),Vmag_5(i__), + rmag_5(i__),imag_5(i__) RAdeg_5(i__) = rNULL__ DEdeg_10(i__) = rNULL__ c Derive coordinates RAdeg_5 and DEdeg_10 from input data c (RAdeg_5 and DEdeg_10 are set to rNULL__ when unknown) if(RAhour_5(i__) .GT. -180) RAdeg_5(i__)=RAhour_5(i__)*15. if(DEdeg_11(i__) .GE. 0) DEdeg_10(i__)=DEdeg_11(i__) if(DEd_5(i__) .GE. 0) DEdeg_10(i__)=DEdeg_10(i__)+DEd_5(i__) if(DEm_5(i__) .GE. 0) DEdeg_10(i__)=DEdeg_10(i__)+DEm_5(i__)/60. if(DEs_5(i__) .GE. 0) DEdeg_10(i__)=DEdeg_10(i__)+DEs_5(i__)/3600. if(DE__5(i__).EQ.'-'.AND.DEdeg_10(i__).GE.0) DEdeg_10(i__)=-DEdeg_10(i__) c ..............Just test output........... write(6,6) + ID_5(i__),k_5(i__),RAhour_5(i__),e_RAhour_5(i__), + DEdeg_11(i__),e_DEdeg_5(i__),Bmag_5(i__),e_Bmag_5(i__), + Epoch_5(i__),FWHM_5(i__),cI_5(i__),RAh_5(i__),RAm_5(i__), + RAs_5(i__),DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__), + f_mag_5(i__),pmRA_5(i__),pmDE_5(i__),B4mag_5(i__),Vmag_5(i__), + rmag_5(i__),imag_5(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_5(i__),DEdeg_10(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end