Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/842/6 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-20
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/842/6 Astro-photometric catalog of the core of NGC 5139 (Bellini+, 2017)
*================================================================================
*The state-of-the-art HST astro-photometric analysis of the core
*of {omega} Centauri.
*I. The catalog.
* Bellini A., Anderson J., Bedin L.R., King I.R., van der Marel R.P.,
* Piotto G., Cool A.
* <Astrophys. J., 842, 6-6 (2017)>
* =2017ApJ...842....6B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! List of HST WFC3 observations of the core
of NGC5139
integer*4 nr__
parameter (nr__=26) ! Number of records
character*131 ar__ ! Full-size record
character*9 Inst (nr__) ! Instrument (WFC3/UVIS or WFC3/IR)
character*6 Filt (nr__) ! Filter name (1)
character*22 Exp (nr__) ! (s) Exposures
integer*4 Time (nr__) ! (s) [908/27600] Total time
integer*4 PID1 (nr__) ! Program ID
integer*4 PID2 (nr__) ! ? Second program ID if any
integer*4 PID3 (nr__) ! ? Third program ID if any
integer*4 PID4 (nr__) ! ? Fourth program ID if any
integer*4 PID5 (nr__) ! ? Fifth program ID if any
integer*4 PID6 (nr__) ! ? Sixth program ID if any
integer*4 PID7 (nr__) ! ? Seventh program ID if any
integer*4 PID8 (nr__) ! ? Eighth program ID if any
character*27 Comm (nr__) ! Description
character*9 Date (nr__) ! Epoch
*Note (1): Photometric-calibration Zero Points (table 2):
* ------------------------------------------------------
* Filter Delmag ZP (VEGA)
* (mag) (mag)
* ------------------------------------------------------
* WFC3/UVIS
* ------------------------------------------------------
* F225W +7.0889 (+7.3856 -0.2967) +22.3808
* F275W +7.0240 (+7.2577 -0.2337) +22.6322
* F336W +6.1648 (+6.3602 -0.1954) +23.4836
* F350LP +6.1925 (+6.3602 -0.1677) +26.7874
* F390M +6.1867 (+6.3602 -0.1735) +23.5377
* F390W +6.2007 (+6.3602 -0.1595) +25.1413
* F438W +6.1961 (+6.3602 -0.1641) +24.9738
* F467M +6.4827 (+6.6330 -0.1503) +23.8362
* F555W +3.8302 (+4.0051 -0.1749) +25.8160
* F606W +3.8517 (+4.0051 -0.1534) +25.9866
* F621M +6.9822 (+7.1251 -0.1429) +24.4539
* F656N +6.6090 (+6.7474 -0.1384) +19.8215
* F658N +6.2200 (+6.3602 -0.1402) +20.6795
* F673N +6.2313 (+6.3602 -0.1289) +22.3297
* F775W +6.4621 (+6.6330 -0.1709) +24.4747
* F814W +3.8445 (+4.0051 -0.1616) +24.6803
* F850LP +4.2503 (+4.4454 -0.1951) +23.3130
* F953N +7.1632 (+7.3235 -0.1603) +19.7549
* ------------------------------------------------------
* WFC3/IR
* ------------------------------------------------------
* F098M -0.0736 +25.1057
* F105W -0.0335 +25.6236
* F110W -0.0743 +26.0628
* F125W -0.0844 +25.3293
* F139M -0.0633 +23.4006
* F140W -0.1262 +25.3761
* F153M -0.0636 +23.2098
* F160W -0.0756 +24.6949
* ------------------------------------------------------
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! *Astrometric data
integer*4 nr__1
parameter (nr__1=478477) ! Number of records
character*214 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 ID (nr__1) ! [1/478477] Identification number
real*8 Xpix (nr__1) ! (pix) X-position on the master frame
real*8 Ypix (nr__1) ! (pix) Y-position on the master frame
integer*4 RAh (nr__1) ! (h) [13] Hour of Right Ascension (J2000)
integer*4 RAm (nr__1) ! (min) [26/27] Minute of Right Ascension (J2000)
real*4 RAs (nr__1) ! (s) Second of Right Ascension (J2000)
character*1 DE_ (nr__1) ! [-] Sign of the Declination (J2000)
integer*4 DEd (nr__1) ! (deg) [47] Degree of Declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) [25/31] Arcminute of Declination (J2000)
real*4 DEs (nr__1) ! (arcsec) Arcsecond of Declination (J2000)
real*8 pmRA (nr__1) ! (mas/yr) [-18/12]?=-999.999 Proper motion along RA
* times cos(DE)
real*8 pmDE (nr__1) ! (mas/yr) [-14/16]?=-999.999 Proper motion along DE
real*8 e_pmRA (nr__1) ! (mas/yr) [0.003/4]?=-999.999 1{sigma} uncertainty
* in pmRA computed using actual residuals
real*8 e_pmDE (nr__1) ! (mas/yr) [0.005/4]?=-999.999 1{sigma} uncertainty
* in pmDE computed using actual residuals
real*8 epmRA (nr__1) ! (mas/yr) [0.004/4]?=-999.999 The 1{sigma} uncertainty
* in pmRA computed using expected errors
real*8 epmDE (nr__1) ! (mas/yr) [0.004/4]?=-999.999 The 1{sigma} uncertainty
* in pmDE computed using expected errors
real*8 QpmRA (nr__1) ! [0.1/67]?=-999.999 Reduced {chi}^2^ of fit
* to pmRA
real*8 QpmDE (nr__1) ! [0.1/87]?=-999.999 Reduced {chi}^2^ of fit
* to pmDE
real*8 sigmaX (nr__1) ! (pix) [0.0005/0.4]?=-999.999 The 1{sigma} uncertainty
* in intercept of PM fit for RA component
* computed using actual residuals
real*8 sigmaY (nr__1) ! (pix) [0.0005/0.4]?=-999.999 The 1{sigma} uncertainty
* in intercept of PM fit for DE component
* computed using actual residuals
real*8 Time_1 (nr__1) ! (yr) [0.7/11]?=-999.999 Time baseline; Julian years
real*8 errX (nr__1) ! (pix) [0.0006/0.3]?=-999.999 The 1{sigma} uncertainty
* in intercept of PM fit for RA component
* using expected errors
real*8 errY (nr__1) ! (pix) [0.0005/0.4]?=-999.999 The 1{sigma} uncertainty
* in intercept of PM fit for DE component
* using expected errors
integer*4 Uref (nr__1) ! [0/1]?=-999 Reference star flag; 1=Yes, 0=No
integer*4 Nfound (nr__1) ! [7/281]?=-999 Total number of data points
* available for PM fits
integer*4 Nused (nr__1) ! [5/281]?=-999 Final number of data points used
* for PM fits
real*8 DelpmRA (nr__1) ! (mas/yr) [-0.5/0.5]?=-999.999 Difference in pmRA (1)
real*8 DelpmDE (nr__1) ! (mas/yr) [-0.5/0.6]?=-999.999 Difference in pmDE (1)
*Note (1): Between locally-corrected and amplifier-based PMs.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'method1.dat' ! *The Method-one photometric file
integer*4 nr__2
parameter (nr__2=12440402) ! Number of records
character*83 ar__2 ! Full-size record
character*6 Filt_1 (nr__2) ! Filter identifier; see Table 1
integer*4 ID_1 (nr__2) ! [1/478477] Identification number as in Table 3;
* column added by CDS
real*8 mag (nr__2) ! (mag) [9.1/32.6]?=-99.999 Apparent Vega magnitude
* in Filt (1)
real*8 e_mag (nr__2) ! (mag) [0/96973]?=-99.999 RMS uncertainty in mag (1)
integer*4 Nf (nr__2) ! [0/71] Number of exposures a star found in
integer*4 Ng (nr__2) ! [0/71] Number of good measurements
real*4 QFIT (nr__2) ! [0/1]? Quality-of-fit parameter
character*1 f_QFIT (nr__2) ! [*] *: QFIT not available (only one measurement
* in method3)
real*8 fL (nr__2) ! [-18.3/20124] Fraction of light in aperture
* from neighbors
real*4 sky (nr__2) ! (ct) [0/999.9]? Local sky-background value
* (only for Method One)
real*4 RMS (nr__2) ! (ct) [0/999.9]? RMS of sky
* (only for Method One)
integer*4 Sat (nr__2) ! [0/1]? Saturation flag, UVIS filters only:
* 1=saturated (only for Method One)
*Note (1): If a star is found in only one exposure, it is not
* possible to compute its photometric rms, and a flag value of -99.999
* is used. If a star is not measured in one particular filter, both its
* magnitude and rms will be flagged at -99.999.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'method2.dat' ! *The Method-two photometric file
integer*4 nr__3
parameter (nr__3=12440402) ! Number of records
character*83 ar__3 ! Full-size record
character*6 Filt_2 (nr__3) ! Filter identifier; see Table 1
integer*4 ID_2 (nr__3) ! [1/478477] Identification number as in Table 3;
* column added by CDS
real*8 mag_1 (nr__3) ! (mag) [9.1/32.6]?=-99.999 Apparent Vega magnitude
* in Filt (1)
real*8 e_mag_1 (nr__3) ! (mag) [0/96973]?=-99.999 RMS uncertainty in mag (1)
integer*4 Nf_1 (nr__3) ! [0/71] Number of exposures a star found in
integer*4 Ng_1 (nr__3) ! [0/71] Number of good measurements
real*4 QFIT_1 (nr__3) ! [0/1]? Quality-of-fit parameter
character*1 f_QFIT_1 (nr__3) ! [*] *: QFIT not available (only one measurement
* in method3)
real*8 fL_1 (nr__3) ! [-18.3/20124] Fraction of light in aperture
* from neighbors
real*4 sky_1 (nr__3) ! (ct) [0/999.9]? Local sky-background value
* (only for Method One)
real*4 RMS_1 (nr__3) ! (ct) [0/999.9]? RMS of sky
* (only for Method One)
integer*4 Sat_1 (nr__3) ! [0/1]? Saturation flag, UVIS filters only:
* 1=saturated (only for Method One)
*Note (1): If a star is found in only one exposure, it is not
* possible to compute its photometric rms, and a flag value of -99.999
* is used. If a star is not measured in one particular filter, both its
* magnitude and rms will be flagged at -99.999.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'method3.dat' ! *The Method-three photometric file
integer*4 nr__4
parameter (nr__4=12440402) ! Number of records
character*83 ar__4 ! Full-size record
character*6 Filt_3 (nr__4) ! Filter identifier; see Table 1
integer*4 ID_3 (nr__4) ! [1/478477] Identification number as in Table 3;
* column added by CDS
real*8 mag_2 (nr__4) ! (mag) [9.1/32.6]?=-99.999 Apparent Vega magnitude
* in Filt (1)
real*8 e_mag_2 (nr__4) ! (mag) [0/96973]?=-99.999 RMS uncertainty in mag (1)
integer*4 Nf_2 (nr__4) ! [0/71] Number of exposures a star found in
integer*4 Ng_2 (nr__4) ! [0/71] Number of good measurements
real*4 QFIT_2 (nr__4) ! [0/1]? Quality-of-fit parameter
character*1 f_QFIT_2 (nr__4) ! [*] *: QFIT not available (only one measurement
* in method3)
real*8 fL_2 (nr__4) ! [-18.3/20124] Fraction of light in aperture
* from neighbors
real*4 sky_2 (nr__4) ! (ct) [0/999.9]? Local sky-background value
* (only for Method One)
real*4 RMS_2 (nr__4) ! (ct) [0/999.9]? RMS of sky
* (only for Method One)
integer*4 Sat_2 (nr__4) ! [0/1]? Saturation flag, UVIS filters only:
* 1=saturated (only for Method One)
*Note (1): If a star is found in only one exposure, it is not
* possible to compute its photometric rms, and a flag value of -99.999
* is used. If a star is not measured in one particular filter, both its
* magnitude and rms will be flagged at -99.999.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Artificial-star-tests file
integer*4 nr__5
parameter (nr__5=550000) ! Number of records
character*107 ar__5 ! Full-size record
real*8 XInpix (nr__5) ! (pix) Input X pixel position
real*8 YInpix (nr__5) ! (pix) Input Y pixel position
real*8 F606WInmag (nr__5) ! (mag) [-14/-3] Input F606W instrumental magnitude
real*8 F438WInmag (nr__5) ! (mag) [-14/-0.08] Input F438W instrumental magnitude
real*8 XRecpix (nr__5) ! (pix) Recovered X pixel position
real*8 YRecpix (nr__5) ! (pix) Recovered Y pixel position
real*8 F606WM1mag (nr__5) ! (mag) [-90.1/-0] Recovered method 1 F606W magnitude
real*8 F438WM1mag (nr__5) ! (mag) [-90.1/-0] Recovered method 1 F438W magnitude
real*8 F606WM2mag (nr__5) ! (mag) [-90.1/-0] Recovered method 2 F606W magnitude
real*8 F438WM2mag (nr__5) ! (mag) [-90.1/-0] Recovered method 2 F438W magnitude
real*8 F606WM3mag (nr__5) ! (mag) [-90.1/-0] Recovered method 3 F606W magnitude
real*8 F438WM3mag (nr__5) ! (mag) [-90.1/-0] Recovered method 3 F438W magnitude
C=============================================================================
C Loading file 'table1.dat' ! List of HST WFC3 observations of the core
* of NGC5139
C Format for file interpretation
1 format(
+ A9,1X,A6,1X,A22,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,
+ I5,1X,I5,1X,A27,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,26
read(1,'(A131)')ar__
read(ar__,1)
+ Inst(i__),Filt(i__),Exp(i__),Time(i__),PID1(i__),PID2(i__),
+ PID3(i__),PID4(i__),PID5(i__),PID6(i__),PID7(i__),PID8(i__),
+ Comm(i__),Date(i__)
if(ar__(53:57) .EQ. '') PID2(i__) = iNULL__
if(ar__(59:63) .EQ. '') PID3(i__) = iNULL__
if(ar__(65:69) .EQ. '') PID4(i__) = iNULL__
if(ar__(71:75) .EQ. '') PID5(i__) = iNULL__
if(ar__(77:81) .EQ. '') PID6(i__) = iNULL__
if(ar__(83:87) .EQ. '') PID7(i__) = iNULL__
if(ar__(89:93) .EQ. '') PID8(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Inst(i__),Filt(i__),Exp(i__),Time(i__),PID1(i__),PID2(i__),
+ PID3(i__),PID4(i__),PID5(i__),PID6(i__),PID7(i__),PID8(i__),
+ Comm(i__),Date(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! *Astrometric data
C Format for file interpretation
2 format(
+ I6,1X,F7.2,1X,F7.2,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,
+ 1X,F9.4,1X,F9.4,1X,F9.4,1X,F9.4,1X,F9.4,1X,F9.4,1X,F9.4,1X,
+ F9.4,1X,F9.4,1X,F9.4,1X,F10.5,1X,F9.4,1X,F9.4,1X,I4,1X,I4,1X,
+ I4,1X,F9.4,1X,F9.4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,478477
read(1,'(A214)')ar__1
read(ar__1,2)
+ ID(i__),Xpix(i__),Ypix(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),pmRA(i__),pmDE(i__),
+ e_pmRA(i__),e_pmDE(i__),epmRA(i__),epmDE(i__),QpmRA(i__),
+ QpmDE(i__),sigmaX(i__),sigmaY(i__),Time_1(i__),errX(i__),
+ errY(i__),Uref(i__),Nfound(i__),Nused(i__),DelpmRA(i__),
+ DelpmDE(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(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,2)
+ ID(i__),Xpix(i__),Ypix(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),pmRA(i__),pmDE(i__),
+ e_pmRA(i__),e_pmDE(i__),epmRA(i__),epmDE(i__),QpmRA(i__),
+ QpmDE(i__),sigmaX(i__),sigmaY(i__),Time_1(i__),errX(i__),
+ errY(i__),Uref(i__),Nfound(i__),Nused(i__),DelpmRA(i__),
+ DelpmDE(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 'method1.dat' ! *The Method-one photometric file
C Format for file interpretation
3 format(
+ A6,1X,I6,2X,F8.4,1X,F10.4,1X,I2,2X,I2,2X,F5.3,A1,1X,F9.2,5X,
+ F5.1,5X,F5.1,3X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'method1.dat')
write(6,*) '....Loading file: method1.dat'
do i__=1,12440402
read(1,'(A83)')ar__2
read(ar__2,3)
+ Filt_1(i__),ID_1(i__),mag(i__),e_mag(i__),Nf(i__),Ng(i__),
+ QFIT(i__),f_QFIT(i__),fL(i__),sky(i__),RMS(i__),Sat(i__)
if(ar__2(44:48) .EQ. '') QFIT(i__) = rNULL__
if(ar__2(65:69) .EQ. '') sky(i__) = rNULL__
if(ar__2(75:79) .EQ. '') RMS(i__) = rNULL__
if(ar__2(83:83) .EQ. '') Sat(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Filt_1(i__),ID_1(i__),mag(i__),e_mag(i__),Nf(i__),Ng(i__),
+ QFIT(i__),f_QFIT(i__),fL(i__),sky(i__),RMS(i__),Sat(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'method2.dat' ! *The Method-two photometric file
C Format for file interpretation
4 format(
+ A6,1X,I6,2X,F8.4,1X,F10.4,1X,I2,2X,I2,2X,F5.3,A1,1X,F9.2,5X,
+ F5.1,5X,F5.1,3X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'method2.dat')
write(6,*) '....Loading file: method2.dat'
do i__=1,12440402
read(1,'(A83)')ar__3
read(ar__3,4)
+ Filt_2(i__),ID_2(i__),mag_1(i__),e_mag_1(i__),Nf_1(i__),
+ Ng_1(i__),QFIT_1(i__),f_QFIT_1(i__),fL_1(i__),sky_1(i__),
+ RMS_1(i__),Sat_1(i__)
if(ar__3(44:48) .EQ. '') QFIT_1(i__) = rNULL__
if(ar__3(65:69) .EQ. '') sky_1(i__) = rNULL__
if(ar__3(75:79) .EQ. '') RMS_1(i__) = rNULL__
if(ar__3(83:83) .EQ. '') Sat_1(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ Filt_2(i__),ID_2(i__),mag_1(i__),e_mag_1(i__),Nf_1(i__),
+ Ng_1(i__),QFIT_1(i__),f_QFIT_1(i__),fL_1(i__),sky_1(i__),
+ RMS_1(i__),Sat_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'method3.dat' ! *The Method-three photometric file
C Format for file interpretation
5 format(
+ A6,1X,I6,2X,F8.4,1X,F10.4,1X,I2,2X,I2,2X,F5.3,A1,1X,F9.2,5X,
+ F5.1,5X,F5.1,3X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'method3.dat')
write(6,*) '....Loading file: method3.dat'
do i__=1,12440402
read(1,'(A83)')ar__4
read(ar__4,5)
+ Filt_3(i__),ID_3(i__),mag_2(i__),e_mag_2(i__),Nf_2(i__),
+ Ng_2(i__),QFIT_2(i__),f_QFIT_2(i__),fL_2(i__),sky_2(i__),
+ RMS_2(i__),Sat_2(i__)
if(ar__4(44:48) .EQ. '') QFIT_2(i__) = rNULL__
if(ar__4(65:69) .EQ. '') sky_2(i__) = rNULL__
if(ar__4(75:79) .EQ. '') RMS_2(i__) = rNULL__
if(ar__4(83:83) .EQ. '') Sat_2(i__) = iNULL__
c ..............Just test output...........
write(6,5)
+ Filt_3(i__),ID_3(i__),mag_2(i__),e_mag_2(i__),Nf_2(i__),
+ Ng_2(i__),QFIT_2(i__),f_QFIT_2(i__),fL_2(i__),sky_2(i__),
+ RMS_2(i__),Sat_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Artificial-star-tests file
C Format for file interpretation
6 format(
+ F9.4,1X,F9.4,1X,F8.4,1X,F8.4,1X,F7.2,1X,F7.2,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,550000
read(1,'(A107)')ar__5
read(ar__5,6)
+ XInpix(i__),YInpix(i__),F606WInmag(i__),F438WInmag(i__),
+ XRecpix(i__),YRecpix(i__),F606WM1mag(i__),F438WM1mag(i__),
+ F606WM2mag(i__),F438WM2mag(i__),F606WM3mag(i__),
+ F438WM3mag(i__)
c ..............Just test output...........
write(6,6)
+ XInpix(i__),YInpix(i__),F606WInmag(i__),F438WInmag(i__),
+ XRecpix(i__),YRecpix(i__),F606WM1mag(i__),F438WM1mag(i__),
+ F606WM2mag(i__),F438WM2mag(i__),F606WM3mag(i__),
+ F438WM3mag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end