Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/507/655 into FORTRAN code for reading data files line by line.
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-19
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/507/655 VI photometry of new Cepheids in NGC 2541 (Ferrarese+, 1998)
*================================================================================
*The Hubble Space Telescope Key Project on the extragalactic distance scale.
*XII. The discovery of cepheids and a new distance to NGC 2541.
* Ferrarese L., Bresolin F., Kennicutt R.C.Jr, Saha A., Stetson P.B.,
* Freedman W.L., Mould J.R., Madore B.F., Sakai S., Ford H.C., Gibson B.K.,
* Graham J.A., Han M., Hoessel J.G., Huchra J., Hughes S.M.,
* Illingworth G.D., Phelps R., Prosser C.F., Silbermann N.A.
* <Astrophys. J. 507, 655 (1998)>
* =1998ApJ...507..655F (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Astrometry for the NGC 2541 Cepheids
integer*4 nr__
parameter (nr__=34) ! Number of records
character*46 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*3 v_FBK98_ ! Identification (CNN or VNN) (1)
real*4 Xpos ! (pix) X position (2)
real*4 Ypos ! (pix) Y position (2)
integer*4 RAh ! (h) Right ascension (J2000) (3)
integer*4 RAm ! (min) Right ascension (J2000) (3)
real*4 RAs ! (s) Right ascension (J2000) (3)
character*1 DE_ ! Declination sign (J2000)
integer*4 DEd ! (deg) Declination (J2000) (3)
integer*4 DEm ! (arcmin) Declination (J2000) (3)
real*4 DEs ! (arcsec) Declination (J2000) (3)
character*4 Chip ! Chip (PC, WF 2, WF 3, WF 4)
*Note (1): CNN for Cepheids stars, VNN for other variable stars
*Note (2): The X- and Y-coordinates are relative to the 1995 October 30 epoch.
* For each 800 x 800 pixels chip, pixel [1, 1] is located at the edge
* of the pyramid.
*Note (3): RA and DE are calculated using the IRAF task
* STSDAS.HST CALIB.WFPC.METRIC, version 1.3.5 (1996 July).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! V (F555W) photometry for NGC 2541 Cepheids
integer*4 nr__1
parameter (nr__1=437) ! Number of records
character*26 ar__1 ! Full-size record
character*3 v_FBK98__1 ! Identification (CNN or VNN)
real*8 HJD ! (d) Heliocentric Julian date
real*4 mag ! (mag) V or I magnitude
real*4 e_mag ! (mag) rms uncertainty on mag
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! I (F814W) photometry for NGC 2541 Cepheids
integer*4 nr__2
parameter (nr__2=159) ! Number of records
character*26 ar__2 ! Full-size record
character*3 v_FBK98__2 ! Identification (CNN or VNN)
real*8 HJD_1 ! (d) Heliocentric Julian date
real*4 mag_1 ! (mag) V or I magnitude
real*4 e_mag_1 ! (mag) rms uncertainty on mag
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! DOPHOT mean photometry of the NGC 2541 Cepheids
integer*4 nr__3
parameter (nr__3=34) ! Number of records
character*79 ar__3 ! Full-size record
character*3 v_FBK98__3 ! Identification (CNN or VNN)
character*1 l_Per ! Limit flag on Per
real*4 Per ! (d) Period
character*1 l_Per2 ! Limit flag on Per2
integer*4 Per2 ! (d) ? Upper limit for period interval
real*4 F555Wph ! (mag) Phase-weighted mean F555W magnitude
real*4 F555Wav ! (mag) Intensity-averaged mean F555W magnitude
real*4 F814Wph ! (mag) ? Phase-weighted mean F814W magnitude
real*4 F814WphCor ! (mag) ? Correction adopted for F814Wph (1)
real*4 F814Wav ! (mag) ? Intensity-average mean F814W magnitude
real*4 F814WavCor ! (mag) ? Correction adopted for F814Wav (1)
real*4 Vmagph ! (mag) ? Phase-weighted mean V magnitude
real*4 Vmagav ! (mag) ? Intensity-averaged mean V magnitude
real*4 Imagph ! (mag) ? Phase-weighted mean I magnitude
real*4 Imagav ! (mag) ? Intensity-averaged mean I magnitude
real*4 chi2_sigma ! ? Reduced {xi}^2^_r_ value for DoPHOT or
* {sigma} for ALLFRAME
*Note (1): Value of the magnitude correction applied to the F814W magnitudes,
* defined as in section 4.4 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Allframe mean photometry of the NGC 2541 Cepheids
integer*4 nr__4
parameter (nr__4=34) ! Number of records
character*79 ar__4 ! Full-size record
character*3 v_FBK98__4 ! Identification (CNN or VNN)
character*1 l_Per_1 ! Limit flag on Per
real*4 Per_1 ! (d) Period
character*1 l_Per2_1 ! Limit flag on Per2
integer*4 Per2_1 ! (d) ? Upper limit for period interval
real*4 F555Wph_1 ! (mag) Phase-weighted mean F555W magnitude
real*4 F555Wav_1 ! (mag) Intensity-averaged mean F555W magnitude
real*4 F814Wph_1 ! (mag) ? Phase-weighted mean F814W magnitude
real*4 F814WphCor_1 ! (mag) ? Correction adopted for F814Wph (1)
real*4 F814Wav_1 ! (mag) ? Intensity-average mean F814W magnitude
real*4 F814WavCor_1 ! (mag) ? Correction adopted for F814Wav (1)
real*4 Vmagph_1 ! (mag) ? Phase-weighted mean V magnitude
real*4 Vmagav_1 ! (mag) ? Intensity-averaged mean V magnitude
real*4 Imagph_1 ! (mag) ? Phase-weighted mean I magnitude
real*4 Imagav_1 ! (mag) ? Intensity-averaged mean I magnitude
real*4 chi2_sigma_1 ! ? Reduced {xi}^2^_r_ value for DoPHOT or
* {sigma} for ALLFRAME
*Note (1): Value of the magnitude correction applied to the F814W magnitudes,
* defined as in section 4.4 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table12.dat' ! Astrometry for other NGC 2541 variable candidates
integer*4 nr__5
parameter (nr__5=23) ! Number of records
character*46 ar__5 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension J2000
real*8 DEdeg_1 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*3 v_FBK98__5 ! Identification (CNN or VNN) (1)
real*4 Xpos_1 ! (pix) X position (2)
real*4 Ypos_1 ! (pix) Y position (2)
integer*4 RAh_1 ! (h) Right ascension (J2000) (3)
integer*4 RAm_1 ! (min) Right ascension (J2000) (3)
real*4 RAs_1 ! (s) Right ascension (J2000) (3)
character*1 DE__1 ! Declination sign (J2000)
integer*4 DEd_1 ! (deg) Declination (J2000) (3)
integer*4 DEm_1 ! (arcmin) Declination (J2000) (3)
real*4 DEs_1 ! (arcsec) Declination (J2000) (3)
character*4 Chip_1 ! Chip (PC, WF 2, WF 3, WF 4)
*Note (1): CNN for Cepheids stars, VNN for other variable stars
*Note (2): The X- and Y-coordinates are relative to the 1995 October 30 epoch.
* For each 800 x 800 pixels chip, pixel [1, 1] is located at the edge
* of the pyramid.
*Note (3): RA and DE are calculated using the IRAF task
* STSDAS.HST CALIB.WFPC.METRIC, version 1.3.5 (1996 July).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table13.dat' ! V photometry for other NGC 2541 variable candidates
integer*4 nr__6
parameter (nr__6=297) ! Number of records
character*26 ar__6 ! Full-size record
character*3 v_FBK98__6 ! Identification (CNN or VNN)
real*8 HJD_2 ! (d) Heliocentric Julian date
real*4 mag_2 ! (mag) V or I magnitude
real*4 e_mag_2 ! (mag) rms uncertainty on mag
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table14.dat' ! I photometry for other NGC 2541 variable candidates
integer*4 nr__7
parameter (nr__7=106) ! Number of records
character*26 ar__7 ! Full-size record
character*3 v_FBK98__7 ! Identification (CNN or VNN)
real*8 HJD_3 ! (d) Heliocentric Julian date
real*4 mag_3 ! (mag) V or I magnitude
real*4 e_mag_3 ! (mag) rms uncertainty on mag
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table15.dat' ! DOPHOT mean photometry of other NGC 2541 variables
integer*4 nr__8
parameter (nr__8=22) ! Number of records
character*79 ar__8 ! Full-size record
character*3 v_FBK98__8 ! Identification (CNN or VNN)
character*1 l_Per_2 ! Limit flag on Per
real*4 Per_2 ! (d) Period
character*1 l_Per2_2 ! Limit flag on Per2
integer*4 Per2_2 ! (d) ? Upper limit for period interval
real*4 F555Wph_2 ! (mag) Phase-weighted mean F555W magnitude
real*4 F555Wav_2 ! (mag) Intensity-averaged mean F555W magnitude
real*4 F814Wph_2 ! (mag) ? Phase-weighted mean F814W magnitude
real*4 F814WphCor_2 ! (mag) ? Correction adopted for F814Wph (1)
real*4 F814Wav_2 ! (mag) ? Intensity-average mean F814W magnitude
real*4 F814WavCor_2 ! (mag) ? Correction adopted for F814Wav (1)
real*4 Vmagph_2 ! (mag) ? Phase-weighted mean V magnitude
real*4 Vmagav_2 ! (mag) ? Intensity-averaged mean V magnitude
real*4 Imagph_2 ! (mag) ? Phase-weighted mean I magnitude
real*4 Imagav_2 ! (mag) ? Intensity-averaged mean I magnitude
real*4 chi2_sigma_2 ! ? Reduced {xi}^2^_r_ value for DoPHOT or
* {sigma} for ALLFRAME
*Note (1): Value of the magnitude correction applied to the F814W magnitudes,
* defined as in section 4.4 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table16.dat' ! Allframe mean photometry of other NGC 2541 variables
integer*4 nr__9
parameter (nr__9=11) ! Number of records
character*79 ar__9 ! Full-size record
character*3 v_FBK98__9 ! Identification (CNN or VNN)
character*1 l_Per_3 ! Limit flag on Per
real*4 Per_3 ! (d) Period
character*1 l_Per2_3 ! Limit flag on Per2
integer*4 Per2_3 ! (d) ? Upper limit for period interval
real*4 F555Wph_3 ! (mag) Phase-weighted mean F555W magnitude
real*4 F555Wav_3 ! (mag) Intensity-averaged mean F555W magnitude
real*4 F814Wph_3 ! (mag) ? Phase-weighted mean F814W magnitude
real*4 F814WphCor_3 ! (mag) ? Correction adopted for F814Wph (1)
real*4 F814Wav_3 ! (mag) ? Intensity-average mean F814W magnitude
real*4 F814WavCor_3 ! (mag) ? Correction adopted for F814Wav (1)
real*4 Vmagph_3 ! (mag) ? Phase-weighted mean V magnitude
real*4 Vmagav_3 ! (mag) ? Intensity-averaged mean V magnitude
real*4 Imagph_3 ! (mag) ? Phase-weighted mean I magnitude
real*4 Imagav_3 ! (mag) ? Intensity-averaged mean I magnitude
real*4 chi2_sigma_3 ! ? Reduced {xi}^2^_r_ value for DoPHOT or
* {sigma} for ALLFRAME
*Note (1): Value of the magnitude correction applied to the F814W magnitudes,
* defined as in section 4.4 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table17.dat' ! Photometry for the PC secondary standards
integer*4 nr__10
parameter (nr__10=24) ! Number of records
character*71 ar__10 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_2 ! (deg) Right Ascension J2000
real*8 DEdeg_2 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
real*4 Xpos_2 ! (pix) X position (1)
real*4 Ypos_2 ! (pix) Y position (1)
integer*4 RAh_2 ! (h) Right ascension (J2000) (2)
integer*4 RAm_2 ! (min) Right ascension (J2000) (2)
real*4 RAs_2 ! (s) Right ascension (J2000) (2)
character*1 DE__2 ! Declination sign (J2000) (2)
integer*4 DEd_2 ! (deg) Declination (J2000) (2)
integer*4 DEm_2 ! (arcmin) Declination (J2000) (2)
real*4 DEs_2 ! (arcsec) Declination (J2000) (2)
real*4 F555W ! (mag) F555W magnitude(3)
real*4 e_F555W ! (mag) rms uncertainty on F555W
real*4 F814W ! (mag) F814W magnitude (3)
real*4 e_F814W ! (mag) rms uncertainty on F814W
real*4 Vmag ! (mag) V magnitude (4)
real*4 Imag ! (mag) I magnitude (4)
*Note (1): The X and Y coordinate are relative to the 1995 October 30 epoch.
* For each 800x800 pixels chip, pixel [1,1] is located at the edge
* of the WFPC2 pyramid
*Note (2): RA and DE are calculated using the IRAF task
* STSDAS.HST_CALIB.WFPC.METRIC, version 1.3.5 (July 1996)
*Note (3): F555W and F814W instrumental magnitudes are intensity averaged
* over 13 and 5 epochs respectively.
*Note (4): V and I magnitudes are calculated from the instrumental
* F555W and F814W magnitudes as described in section 3 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table18.dat' ! Photometry for the WF2 secondary standards
integer*4 nr__11
parameter (nr__11=29) ! Number of records
character*71 ar__11 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_3 ! (deg) Right Ascension J2000
real*8 DEdeg_3 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
real*4 Xpos_3 ! (pix) X position (1)
real*4 Ypos_3 ! (pix) Y position (1)
integer*4 RAh_3 ! (h) Right ascension (J2000) (2)
integer*4 RAm_3 ! (min) Right ascension (J2000) (2)
real*4 RAs_3 ! (s) Right ascension (J2000) (2)
character*1 DE__3 ! Declination sign (J2000) (2)
integer*4 DEd_3 ! (deg) Declination (J2000) (2)
integer*4 DEm_3 ! (arcmin) Declination (J2000) (2)
real*4 DEs_3 ! (arcsec) Declination (J2000) (2)
real*4 F555W_1 ! (mag) F555W magnitude(3)
real*4 e_F555W_1 ! (mag) rms uncertainty on F555W
real*4 F814W_1 ! (mag) F814W magnitude (3)
real*4 e_F814W_1 ! (mag) rms uncertainty on F814W
real*4 Vmag_1 ! (mag) V magnitude (4)
real*4 Imag_1 ! (mag) I magnitude (4)
*Note (1): The X and Y coordinate are relative to the 1995 October 30 epoch.
* For each 800x800 pixels chip, pixel [1,1] is located at the edge
* of the WFPC2 pyramid
*Note (2): RA and DE are calculated using the IRAF task
* STSDAS.HST_CALIB.WFPC.METRIC, version 1.3.5 (July 1996)
*Note (3): F555W and F814W instrumental magnitudes are intensity averaged
* over 13 and 5 epochs respectively.
*Note (4): V and I magnitudes are calculated from the instrumental
* F555W and F814W magnitudes as described in section 3 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table19.dat' ! Photometry for the WF3 secondary standards
integer*4 nr__12
parameter (nr__12=32) ! Number of records
character*71 ar__12 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_4 ! (deg) Right Ascension J2000
real*8 DEdeg_4 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
real*4 Xpos_4 ! (pix) X position (1)
real*4 Ypos_4 ! (pix) Y position (1)
integer*4 RAh_4 ! (h) Right ascension (J2000) (2)
integer*4 RAm_4 ! (min) Right ascension (J2000) (2)
real*4 RAs_4 ! (s) Right ascension (J2000) (2)
character*1 DE__4 ! Declination sign (J2000) (2)
integer*4 DEd_4 ! (deg) Declination (J2000) (2)
integer*4 DEm_4 ! (arcmin) Declination (J2000) (2)
real*4 DEs_4 ! (arcsec) Declination (J2000) (2)
real*4 F555W_2 ! (mag) F555W magnitude(3)
real*4 e_F555W_2 ! (mag) rms uncertainty on F555W
real*4 F814W_2 ! (mag) F814W magnitude (3)
real*4 e_F814W_2 ! (mag) rms uncertainty on F814W
real*4 Vmag_2 ! (mag) V magnitude (4)
real*4 Imag_2 ! (mag) I magnitude (4)
*Note (1): The X and Y coordinate are relative to the 1995 October 30 epoch.
* For each 800x800 pixels chip, pixel [1,1] is located at the edge
* of the WFPC2 pyramid
*Note (2): RA and DE are calculated using the IRAF task
* STSDAS.HST_CALIB.WFPC.METRIC, version 1.3.5 (July 1996)
*Note (3): F555W and F814W instrumental magnitudes are intensity averaged
* over 13 and 5 epochs respectively.
*Note (4): V and I magnitudes are calculated from the instrumental
* F555W and F814W magnitudes as described in section 3 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table20.dat' ! Photometry for the WF4 secondary standards
integer*4 nr__13
parameter (nr__13=29) ! Number of records
character*71 ar__13 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_5 ! (deg) Right Ascension J2000
real*8 DEdeg_5 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
real*4 Xpos_5 ! (pix) X position (1)
real*4 Ypos_5 ! (pix) Y position (1)
integer*4 RAh_5 ! (h) Right ascension (J2000) (2)
integer*4 RAm_5 ! (min) Right ascension (J2000) (2)
real*4 RAs_5 ! (s) Right ascension (J2000) (2)
character*1 DE__5 ! Declination sign (J2000) (2)
integer*4 DEd_5 ! (deg) Declination (J2000) (2)
integer*4 DEm_5 ! (arcmin) Declination (J2000) (2)
real*4 DEs_5 ! (arcsec) Declination (J2000) (2)
real*4 F555W_3 ! (mag) F555W magnitude(3)
real*4 e_F555W_3 ! (mag) rms uncertainty on F555W
real*4 F814W_3 ! (mag) F814W magnitude (3)
real*4 e_F814W_3 ! (mag) rms uncertainty on F814W
real*4 Vmag_3 ! (mag) V magnitude (4)
real*4 Imag_3 ! (mag) I magnitude (4)
*Note (1): The X and Y coordinate are relative to the 1995 October 30 epoch.
* For each 800x800 pixels chip, pixel [1,1] is located at the edge
* of the WFPC2 pyramid
*Note (2): RA and DE are calculated using the IRAF task
* STSDAS.HST_CALIB.WFPC.METRIC, version 1.3.5 (July 1996)
*Note (3): F555W and F814W instrumental magnitudes are intensity averaged
* over 13 and 5 epochs respectively.
*Note (4): V and I magnitudes are calculated from the instrumental
* F555W and F814W magnitudes as described in section 3 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'note5.dat' ! Individual notes on Cepheids
integer*4 nr__14
parameter (nr__14=44) ! Number of records
character*80 ar__14 ! Full-size record
character*3 v_FBK98__10 ! Identification (CNN or VNN)
character*76 Note ! Text of the note (1)
*Note (1): The term "isolated" is use d for stars for which no companions are
* identified within a 3 pixel radius.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'note12.dat' ! Individual notes on other variable stars
integer*4 nr__15
parameter (nr__15=46) ! Number of records
character*80 ar__15 ! Full-size record
character*3 v_FBK98__11 ! Identification (CNN or VNN)
character*76 Note_1 ! Text of the note (1)
*Note (1): The term "isolated" is use d for stars for which no companions are
* identified within a 3 pixel radius.
C=============================================================================
C Loading file 'table5.dat' ! Astrometry for the NGC 2541 Cepheids
C Format for file interpretation
1 format(
+ A3,1X,F6.2,1X,F6.2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,
+ 1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,34
read(1,'(A46)')ar__
read(ar__,1)
+ v_FBK98_,Xpos,Ypos,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Chip
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ v_FBK98_,Xpos,Ypos,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Chip
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! V (F555W) photometry for NGC 2541 Cepheids
C Format for file interpretation
2 format(A3,1X,F11.3,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,437
read(1,'(A26)')ar__1
read(ar__1,2)v_FBK98__1,HJD,mag,e_mag
c ..............Just test output...........
write(6,2)v_FBK98__1,HJD,mag,e_mag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! I (F814W) photometry for NGC 2541 Cepheids
C Format for file interpretation
3 format(A3,1X,F11.3,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,159
read(1,'(A26)')ar__2
read(ar__2,3)v_FBK98__2,HJD_1,mag_1,e_mag_1
c ..............Just test output...........
write(6,3)v_FBK98__2,HJD_1,mag_1,e_mag_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! DOPHOT mean photometry of the NGC 2541 Cepheids
C Format for file interpretation
4 format(
+ A3,1X,A1,F4.1,1X,A1,I2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,34
read(1,'(A79)')ar__3
read(ar__3,4)
+ v_FBK98__3,l_Per,Per,l_Per2,Per2,F555Wph,F555Wav,F814Wph,
+ F814WphCor,F814Wav,F814WavCor,Vmagph,Vmagav,Imagph,Imagav,
+ chi2_sigma
if(ar__3(12:13) .EQ. '') Per2 = iNULL__
if(ar__3(27:31) .EQ. '') F814Wph = rNULL__
if(ar__3(33:37) .EQ. '') F814WphCor = rNULL__
if(ar__3(39:43) .EQ. '') F814Wav = rNULL__
if(ar__3(45:49) .EQ. '') F814WavCor = rNULL__
if(ar__3(51:55) .EQ. '') Vmagph = rNULL__
if(ar__3(57:61) .EQ. '') Vmagav = rNULL__
if(ar__3(63:67) .EQ. '') Imagph = rNULL__
if(ar__3(69:73) .EQ. '') Imagav = rNULL__
if(ar__3(75:79) .EQ. '') chi2_sigma = rNULL__
c ..............Just test output...........
write(6,4)
+ v_FBK98__3,l_Per,Per,l_Per2,Per2,F555Wph,F555Wav,F814Wph,
+ F814WphCor,F814Wav,F814WavCor,Vmagph,Vmagav,Imagph,Imagav,
+ chi2_sigma
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Allframe mean photometry of the NGC 2541 Cepheids
C Format for file interpretation
5 format(
+ A3,1X,A1,F4.1,1X,A1,I2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,34
read(1,'(A79)')ar__4
read(ar__4,5)
+ v_FBK98__4,l_Per_1,Per_1,l_Per2_1,Per2_1,F555Wph_1,F555Wav_1,
+ F814Wph_1,F814WphCor_1,F814Wav_1,F814WavCor_1,Vmagph_1,
+ Vmagav_1,Imagph_1,Imagav_1,chi2_sigma_1
if(ar__4(12:13) .EQ. '') Per2_1 = iNULL__
if(ar__4(27:31) .EQ. '') F814Wph_1 = rNULL__
if(ar__4(33:37) .EQ. '') F814WphCor_1 = rNULL__
if(ar__4(39:43) .EQ. '') F814Wav_1 = rNULL__
if(ar__4(45:49) .EQ. '') F814WavCor_1 = rNULL__
if(ar__4(51:55) .EQ. '') Vmagph_1 = rNULL__
if(ar__4(57:61) .EQ. '') Vmagav_1 = rNULL__
if(ar__4(63:67) .EQ. '') Imagph_1 = rNULL__
if(ar__4(69:73) .EQ. '') Imagav_1 = rNULL__
if(ar__4(75:79) .EQ. '') chi2_sigma_1 = rNULL__
c ..............Just test output...........
write(6,5)
+ v_FBK98__4,l_Per_1,Per_1,l_Per2_1,Per2_1,F555Wph_1,F555Wav_1,
+ F814Wph_1,F814WphCor_1,F814Wav_1,F814WavCor_1,Vmagph_1,
+ Vmagav_1,Imagph_1,Imagav_1,chi2_sigma_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table12.dat' ! Astrometry for other NGC 2541 variable candidates
C Format for file interpretation
6 format(
+ A3,1X,F6.2,1X,F6.2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,
+ 1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table12.dat')
write(6,*) '....Loading file: table12.dat'
do i__=1,23
read(1,'(A46)')ar__5
read(ar__5,6)
+ v_FBK98__5,Xpos_1,Ypos_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,
+ DEs_1,Chip_1
RAdeg_1 = rNULL__
DEdeg_1 = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1 .GT. -180) RAdeg_1=RAh_1*15.
if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
if(DEd_1 .GE. 0) DEdeg_1=DEd_1
if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,6)
+ v_FBK98__5,Xpos_1,Ypos_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,
+ DEs_1,Chip_1
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table13.dat' ! V photometry for other NGC 2541 variable candidates
C Format for file interpretation
7 format(A3,1X,F11.3,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table13.dat')
write(6,*) '....Loading file: table13.dat'
do i__=1,297
read(1,'(A26)')ar__6
read(ar__6,7)v_FBK98__6,HJD_2,mag_2,e_mag_2
c ..............Just test output...........
write(6,7)v_FBK98__6,HJD_2,mag_2,e_mag_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table14.dat' ! I photometry for other NGC 2541 variable candidates
C Format for file interpretation
8 format(A3,1X,F11.3,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table14.dat')
write(6,*) '....Loading file: table14.dat'
do i__=1,106
read(1,'(A26)')ar__7
read(ar__7,8)v_FBK98__7,HJD_3,mag_3,e_mag_3
c ..............Just test output...........
write(6,8)v_FBK98__7,HJD_3,mag_3,e_mag_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table15.dat' ! DOPHOT mean photometry of other NGC 2541 variables
C Format for file interpretation
9 format(
+ A3,1X,A1,F4.1,1X,A1,I2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table15.dat')
write(6,*) '....Loading file: table15.dat'
do i__=1,22
read(1,'(A79)')ar__8
read(ar__8,9)
+ v_FBK98__8,l_Per_2,Per_2,l_Per2_2,Per2_2,F555Wph_2,F555Wav_2,
+ F814Wph_2,F814WphCor_2,F814Wav_2,F814WavCor_2,Vmagph_2,
+ Vmagav_2,Imagph_2,Imagav_2,chi2_sigma_2
if(ar__8(12:13) .EQ. '') Per2_2 = iNULL__
if(ar__8(27:31) .EQ. '') F814Wph_2 = rNULL__
if(ar__8(33:37) .EQ. '') F814WphCor_2 = rNULL__
if(ar__8(39:43) .EQ. '') F814Wav_2 = rNULL__
if(ar__8(45:49) .EQ. '') F814WavCor_2 = rNULL__
if(ar__8(51:55) .EQ. '') Vmagph_2 = rNULL__
if(ar__8(57:61) .EQ. '') Vmagav_2 = rNULL__
if(ar__8(63:67) .EQ. '') Imagph_2 = rNULL__
if(ar__8(69:73) .EQ. '') Imagav_2 = rNULL__
if(ar__8(75:79) .EQ. '') chi2_sigma_2 = rNULL__
c ..............Just test output...........
write(6,9)
+ v_FBK98__8,l_Per_2,Per_2,l_Per2_2,Per2_2,F555Wph_2,F555Wav_2,
+ F814Wph_2,F814WphCor_2,F814Wav_2,F814WavCor_2,Vmagph_2,
+ Vmagav_2,Imagph_2,Imagav_2,chi2_sigma_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table16.dat' ! Allframe mean photometry of other NGC 2541 variables
C Format for file interpretation
10 format(
+ A3,1X,A1,F4.1,1X,A1,I2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table16.dat')
write(6,*) '....Loading file: table16.dat'
do i__=1,11
read(1,'(A79)')ar__9
read(ar__9,10)
+ v_FBK98__9,l_Per_3,Per_3,l_Per2_3,Per2_3,F555Wph_3,F555Wav_3,
+ F814Wph_3,F814WphCor_3,F814Wav_3,F814WavCor_3,Vmagph_3,
+ Vmagav_3,Imagph_3,Imagav_3,chi2_sigma_3
if(ar__9(12:13) .EQ. '') Per2_3 = iNULL__
if(ar__9(27:31) .EQ. '') F814Wph_3 = rNULL__
if(ar__9(33:37) .EQ. '') F814WphCor_3 = rNULL__
if(ar__9(39:43) .EQ. '') F814Wav_3 = rNULL__
if(ar__9(45:49) .EQ. '') F814WavCor_3 = rNULL__
if(ar__9(51:55) .EQ. '') Vmagph_3 = rNULL__
if(ar__9(57:61) .EQ. '') Vmagav_3 = rNULL__
if(ar__9(63:67) .EQ. '') Imagph_3 = rNULL__
if(ar__9(69:73) .EQ. '') Imagav_3 = rNULL__
if(ar__9(75:79) .EQ. '') chi2_sigma_3 = rNULL__
c ..............Just test output...........
write(6,10)
+ v_FBK98__9,l_Per_3,Per_3,l_Per2_3,Per2_3,F555Wph_3,F555Wav_3,
+ F814Wph_3,F814WphCor_3,F814Wav_3,F814WavCor_3,Vmagph_3,
+ Vmagav_3,Imagph_3,Imagav_3,chi2_sigma_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table17.dat' ! Photometry for the PC secondary standards
C Format for file interpretation
11 format(
+ F6.2,1X,F6.2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table17.dat')
write(6,*) '....Loading file: table17.dat'
do i__=1,24
read(1,'(A71)')ar__10
read(ar__10,11)
+ Xpos_2,Ypos_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,F555W,
+ e_F555W,F814W,e_F814W,Vmag,Imag
RAdeg_2 = rNULL__
DEdeg_2 = rNULL__
c Derive coordinates RAdeg_2 and DEdeg_2 from input data
c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
if(RAh_2 .GT. -180) RAdeg_2=RAh_2*15.
if(RAm_2 .GT. -180) RAdeg_2=RAdeg_2+RAm_2/4.
if(RAs_2 .GT. -180) RAdeg_2=RAdeg_2+RAs_2/240.
if(DEd_2 .GE. 0) DEdeg_2=DEd_2
if(DEm_2 .GE. 0) DEdeg_2=DEdeg_2+DEm_2/60.
if(DEs_2 .GE. 0) DEdeg_2=DEdeg_2+DEs_2/3600.
if(DE__2.EQ.'-'.AND.DEdeg_2.GE.0) DEdeg_2=-DEdeg_2
c ..............Just test output...........
write(6,11)
+ Xpos_2,Ypos_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,F555W,
+ e_F555W,F814W,e_F814W,Vmag,Imag
write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table18.dat' ! Photometry for the WF2 secondary standards
C Format for file interpretation
12 format(
+ F6.2,1X,F6.2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table18.dat')
write(6,*) '....Loading file: table18.dat'
do i__=1,29
read(1,'(A71)')ar__11
read(ar__11,12)
+ Xpos_3,Ypos_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,DEs_3,
+ F555W_1,e_F555W_1,F814W_1,e_F814W_1,Vmag_1,Imag_1
RAdeg_3 = rNULL__
DEdeg_3 = rNULL__
c Derive coordinates RAdeg_3 and DEdeg_3 from input data
c (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown)
if(RAh_3 .GT. -180) RAdeg_3=RAh_3*15.
if(RAm_3 .GT. -180) RAdeg_3=RAdeg_3+RAm_3/4.
if(RAs_3 .GT. -180) RAdeg_3=RAdeg_3+RAs_3/240.
if(DEd_3 .GE. 0) DEdeg_3=DEd_3
if(DEm_3 .GE. 0) DEdeg_3=DEdeg_3+DEm_3/60.
if(DEs_3 .GE. 0) DEdeg_3=DEdeg_3+DEs_3/3600.
if(DE__3.EQ.'-'.AND.DEdeg_3.GE.0) DEdeg_3=-DEdeg_3
c ..............Just test output...........
write(6,12)
+ Xpos_3,Ypos_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,DEs_3,
+ F555W_1,e_F555W_1,F814W_1,e_F814W_1,Vmag_1,Imag_1
write(6,'(6H Pos: 2F8.4)') RAdeg_3,DEdeg_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table19.dat' ! Photometry for the WF3 secondary standards
C Format for file interpretation
13 format(
+ F6.2,1X,F6.2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table19.dat')
write(6,*) '....Loading file: table19.dat'
do i__=1,32
read(1,'(A71)')ar__12
read(ar__12,13)
+ Xpos_4,Ypos_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_4,
+ F555W_2,e_F555W_2,F814W_2,e_F814W_2,Vmag_2,Imag_2
RAdeg_4 = rNULL__
DEdeg_4 = rNULL__
c Derive coordinates RAdeg_4 and DEdeg_4 from input data
c (RAdeg_4 and DEdeg_4 are set to rNULL__ when unknown)
if(RAh_4 .GT. -180) RAdeg_4=RAh_4*15.
if(RAm_4 .GT. -180) RAdeg_4=RAdeg_4+RAm_4/4.
if(RAs_4 .GT. -180) RAdeg_4=RAdeg_4+RAs_4/240.
if(DEd_4 .GE. 0) DEdeg_4=DEd_4
if(DEm_4 .GE. 0) DEdeg_4=DEdeg_4+DEm_4/60.
if(DEs_4 .GE. 0) DEdeg_4=DEdeg_4+DEs_4/3600.
if(DE__4.EQ.'-'.AND.DEdeg_4.GE.0) DEdeg_4=-DEdeg_4
c ..............Just test output...........
write(6,13)
+ Xpos_4,Ypos_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_4,
+ F555W_2,e_F555W_2,F814W_2,e_F814W_2,Vmag_2,Imag_2
write(6,'(6H Pos: 2F8.4)') RAdeg_4,DEdeg_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table20.dat' ! Photometry for the WF4 secondary standards
C Format for file interpretation
14 format(
+ F6.2,1X,F6.2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table20.dat')
write(6,*) '....Loading file: table20.dat'
do i__=1,29
read(1,'(A71)')ar__13
read(ar__13,14)
+ Xpos_5,Ypos_5,RAh_5,RAm_5,RAs_5,DE__5,DEd_5,DEm_5,DEs_5,
+ F555W_3,e_F555W_3,F814W_3,e_F814W_3,Vmag_3,Imag_3
RAdeg_5 = rNULL__
DEdeg_5 = rNULL__
c Derive coordinates RAdeg_5 and DEdeg_5 from input data
c (RAdeg_5 and DEdeg_5 are set to rNULL__ when unknown)
if(RAh_5 .GT. -180) RAdeg_5=RAh_5*15.
if(RAm_5 .GT. -180) RAdeg_5=RAdeg_5+RAm_5/4.
if(RAs_5 .GT. -180) RAdeg_5=RAdeg_5+RAs_5/240.
if(DEd_5 .GE. 0) DEdeg_5=DEd_5
if(DEm_5 .GE. 0) DEdeg_5=DEdeg_5+DEm_5/60.
if(DEs_5 .GE. 0) DEdeg_5=DEdeg_5+DEs_5/3600.
if(DE__5.EQ.'-'.AND.DEdeg_5.GE.0) DEdeg_5=-DEdeg_5
c ..............Just test output...........
write(6,14)
+ Xpos_5,Ypos_5,RAh_5,RAm_5,RAs_5,DE__5,DEd_5,DEm_5,DEs_5,
+ F555W_3,e_F555W_3,F814W_3,e_F814W_3,Vmag_3,Imag_3
write(6,'(6H Pos: 2F8.4)') RAdeg_5,DEdeg_5
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'note5.dat' ! Individual notes on Cepheids
C Format for file interpretation
15 format(A3,1X,A76)
C Effective file loading
open(unit=1,status='old',file=
+'note5.dat')
write(6,*) '....Loading file: note5.dat'
do i__=1,44
read(1,'(A80)')ar__14
read(ar__14,15)v_FBK98__10,Note
c ..............Just test output...........
write(6,15)v_FBK98__10,Note
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'note12.dat' ! Individual notes on other variable stars
C Format for file interpretation
16 format(A3,1X,A76)
C Effective file loading
open(unit=1,status='old',file=
+'note12.dat')
write(6,*) '....Loading file: note12.dat'
do i__=1,46
read(1,'(A80)')ar__15
read(ar__15,16)v_FBK98__11,Note_1
c ..............Just test output...........
write(6,16)v_FBK98__11,Note_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end