FORTRAN Generation
(/./ftp/cats/J/ApJ/507/655)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/507/655 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-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      (nr__) ! (deg) Right Ascension J2000
      real*8        DEdeg      (nr__) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*3   v_FBK98_   (nr__) ! Identification (CNN or VNN) (1)
      real*4        Xpos       (nr__) ! (pix) X position (2)
      real*4        Ypos       (nr__) ! (pix) Y position (2)
      integer*4     RAh        (nr__) ! (h) Right ascension (J2000) (3)
      integer*4     RAm        (nr__) ! (min) Right ascension (J2000) (3)
      real*4        RAs        (nr__) ! (s) Right ascension (J2000) (3)
      character*1   DE_        (nr__) ! Declination sign (J2000)
      integer*4     DEd        (nr__) ! (deg) Declination (J2000) (3)
      integer*4     DEm        (nr__) ! (arcmin) Declination (J2000) (3)
      real*4        DEs        (nr__) ! (arcsec) Declination (J2000) (3)
      character*4   Chip       (nr__) ! 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 (nr__1) ! Identification (CNN or VNN)
      real*8        HJD        (nr__1) ! (d) Heliocentric Julian date
      real*4        mag        (nr__1) ! (mag) V or I magnitude
      real*4        e_mag      (nr__1) ! (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 (nr__2) ! Identification (CNN or VNN)
      real*8        HJD_1      (nr__2) ! (d) Heliocentric Julian date
      real*4        mag_1      (nr__2) ! (mag) V or I magnitude
      real*4        e_mag_1    (nr__2) ! (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 (nr__3) ! Identification (CNN or VNN)
      character*1   l_Per      (nr__3) ! Limit flag on Per
      real*4        Per        (nr__3) ! (d) Period
      character*1   l_Per2     (nr__3) ! Limit flag on Per2
      integer*4     Per2       (nr__3) ! (d) ? Upper limit for period interval
      real*4        F555Wph    (nr__3) ! (mag) Phase-weighted mean F555W magnitude
      real*4        F555Wav    (nr__3) ! (mag) Intensity-averaged mean F555W magnitude
      real*4        F814Wph    (nr__3) ! (mag) ? Phase-weighted mean F814W magnitude
      real*4        F814WphCor (nr__3) ! (mag) ? Correction adopted for F814Wph (1)
      real*4        F814Wav    (nr__3) ! (mag) ? Intensity-average mean F814W magnitude
      real*4        F814WavCor (nr__3) ! (mag) ? Correction adopted for F814Wav (1)
      real*4        Vmagph     (nr__3) ! (mag) ? Phase-weighted mean V magnitude
      real*4        Vmagav     (nr__3) ! (mag) ? Intensity-averaged mean V magnitude
      real*4        Imagph     (nr__3) ! (mag) ? Phase-weighted mean I magnitude
      real*4        Imagav     (nr__3) ! (mag) ? Intensity-averaged mean I magnitude
      real*4        chi2_sigma (nr__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 '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 (nr__4) ! Identification (CNN or VNN)
      character*1   l_Per_1    (nr__4) ! Limit flag on Per
      real*4        Per_1      (nr__4) ! (d) Period
      character*1   l_Per2_1   (nr__4) ! Limit flag on Per2
      integer*4     Per2_1     (nr__4) ! (d) ? Upper limit for period interval
      real*4        F555Wph_1  (nr__4) ! (mag) Phase-weighted mean F555W magnitude
      real*4        F555Wav_1  (nr__4) ! (mag) Intensity-averaged mean F555W magnitude
      real*4        F814Wph_1  (nr__4) ! (mag) ? Phase-weighted mean F814W magnitude
      real*4        F814WphCor_1(nr__4) ! (mag) ? Correction adopted for F814Wph (1)
      real*4        F814Wav_1  (nr__4) ! (mag) ? Intensity-average mean F814W magnitude
      real*4        F814WavCor_1(nr__4) ! (mag) ? Correction adopted for F814Wav (1)
      real*4        Vmagph_1   (nr__4) ! (mag) ? Phase-weighted mean V magnitude
      real*4        Vmagav_1   (nr__4) ! (mag) ? Intensity-averaged mean V magnitude
      real*4        Imagph_1   (nr__4) ! (mag) ? Phase-weighted mean I magnitude
      real*4        Imagav_1   (nr__4) ! (mag) ? Intensity-averaged mean I magnitude
      real*4        chi2_sigma_1(nr__4) ! ? 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    (nr__5) ! (deg) Right Ascension J2000
      real*8        DEdeg_1    (nr__5) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*3   v_FBK98__5 (nr__5) ! Identification (CNN or VNN) (1)
      real*4        Xpos_1     (nr__5) ! (pix) X position (2)
      real*4        Ypos_1     (nr__5) ! (pix) Y position (2)
      integer*4     RAh_1      (nr__5) ! (h) Right ascension (J2000) (3)
      integer*4     RAm_1      (nr__5) ! (min) Right ascension (J2000) (3)
      real*4        RAs_1      (nr__5) ! (s) Right ascension (J2000) (3)
      character*1   DE__1      (nr__5) ! Declination sign (J2000)
      integer*4     DEd_1      (nr__5) ! (deg) Declination (J2000) (3)
      integer*4     DEm_1      (nr__5) ! (arcmin) Declination (J2000) (3)
      real*4        DEs_1      (nr__5) ! (arcsec) Declination (J2000) (3)
      character*4   Chip_1     (nr__5) ! 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 (nr__6) ! Identification (CNN or VNN)
      real*8        HJD_2      (nr__6) ! (d) Heliocentric Julian date
      real*4        mag_2      (nr__6) ! (mag) V or I magnitude
      real*4        e_mag_2    (nr__6) ! (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 (nr__7) ! Identification (CNN or VNN)
      real*8        HJD_3      (nr__7) ! (d) Heliocentric Julian date
      real*4        mag_3      (nr__7) ! (mag) V or I magnitude
      real*4        e_mag_3    (nr__7) ! (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 (nr__8) ! Identification (CNN or VNN)
      character*1   l_Per_2    (nr__8) ! Limit flag on Per
      real*4        Per_2      (nr__8) ! (d) Period
      character*1   l_Per2_2   (nr__8) ! Limit flag on Per2
      integer*4     Per2_2     (nr__8) ! (d) ? Upper limit for period interval
      real*4        F555Wph_2  (nr__8) ! (mag) Phase-weighted mean F555W magnitude
      real*4        F555Wav_2  (nr__8) ! (mag) Intensity-averaged mean F555W magnitude
      real*4        F814Wph_2  (nr__8) ! (mag) ? Phase-weighted mean F814W magnitude
      real*4        F814WphCor_2(nr__8) ! (mag) ? Correction adopted for F814Wph (1)
      real*4        F814Wav_2  (nr__8) ! (mag) ? Intensity-average mean F814W magnitude
      real*4        F814WavCor_2(nr__8) ! (mag) ? Correction adopted for F814Wav (1)
      real*4        Vmagph_2   (nr__8) ! (mag) ? Phase-weighted mean V magnitude
      real*4        Vmagav_2   (nr__8) ! (mag) ? Intensity-averaged mean V magnitude
      real*4        Imagph_2   (nr__8) ! (mag) ? Phase-weighted mean I magnitude
      real*4        Imagav_2   (nr__8) ! (mag) ? Intensity-averaged mean I magnitude
      real*4        chi2_sigma_2(nr__8) ! ? 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 (nr__9) ! Identification (CNN or VNN)
      character*1   l_Per_3    (nr__9) ! Limit flag on Per
      real*4        Per_3      (nr__9) ! (d) Period
      character*1   l_Per2_3   (nr__9) ! Limit flag on Per2
      integer*4     Per2_3     (nr__9) ! (d) ? Upper limit for period interval
      real*4        F555Wph_3  (nr__9) ! (mag) Phase-weighted mean F555W magnitude
      real*4        F555Wav_3  (nr__9) ! (mag) Intensity-averaged mean F555W magnitude
      real*4        F814Wph_3  (nr__9) ! (mag) ? Phase-weighted mean F814W magnitude
      real*4        F814WphCor_3(nr__9) ! (mag) ? Correction adopted for F814Wph (1)
      real*4        F814Wav_3  (nr__9) ! (mag) ? Intensity-average mean F814W magnitude
      real*4        F814WavCor_3(nr__9) ! (mag) ? Correction adopted for F814Wav (1)
      real*4        Vmagph_3   (nr__9) ! (mag) ? Phase-weighted mean V magnitude
      real*4        Vmagav_3   (nr__9) ! (mag) ? Intensity-averaged mean V magnitude
      real*4        Imagph_3   (nr__9) ! (mag) ? Phase-weighted mean I magnitude
      real*4        Imagav_3   (nr__9) ! (mag) ? Intensity-averaged mean I magnitude
      real*4        chi2_sigma_3(nr__9) ! ? 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    (nr__10) ! (deg) Right Ascension J2000
      real*8        DEdeg_2    (nr__10) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      real*4        Xpos_2     (nr__10) ! (pix) X position (1)
      real*4        Ypos_2     (nr__10) ! (pix) Y position (1)
      integer*4     RAh_2      (nr__10) ! (h) Right ascension (J2000) (2)
      integer*4     RAm_2      (nr__10) ! (min) Right ascension (J2000) (2)
      real*4        RAs_2      (nr__10) ! (s) Right ascension (J2000) (2)
      character*1   DE__2      (nr__10) ! Declination sign (J2000) (2)
      integer*4     DEd_2      (nr__10) ! (deg) Declination (J2000) (2)
      integer*4     DEm_2      (nr__10) ! (arcmin) Declination (J2000) (2)
      real*4        DEs_2      (nr__10) ! (arcsec) Declination (J2000) (2)
      real*4        F555W      (nr__10) ! (mag) F555W magnitude(3)
      real*4        e_F555W    (nr__10) ! (mag) rms uncertainty on F555W
      real*4        F814W      (nr__10) ! (mag) F814W magnitude (3)
      real*4        e_F814W    (nr__10) ! (mag) rms uncertainty on F814W
      real*4        Vmag       (nr__10) ! (mag) V magnitude (4)
      real*4        Imag       (nr__10) ! (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    (nr__11) ! (deg) Right Ascension J2000
      real*8        DEdeg_3    (nr__11) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      real*4        Xpos_3     (nr__11) ! (pix) X position (1)
      real*4        Ypos_3     (nr__11) ! (pix) Y position (1)
      integer*4     RAh_3      (nr__11) ! (h) Right ascension (J2000) (2)
      integer*4     RAm_3      (nr__11) ! (min) Right ascension (J2000) (2)
      real*4        RAs_3      (nr__11) ! (s) Right ascension (J2000) (2)
      character*1   DE__3      (nr__11) ! Declination sign (J2000) (2)
      integer*4     DEd_3      (nr__11) ! (deg) Declination (J2000) (2)
      integer*4     DEm_3      (nr__11) ! (arcmin) Declination (J2000) (2)
      real*4        DEs_3      (nr__11) ! (arcsec) Declination (J2000) (2)
      real*4        F555W_1    (nr__11) ! (mag) F555W magnitude(3)
      real*4        e_F555W_1  (nr__11) ! (mag) rms uncertainty on F555W
      real*4        F814W_1    (nr__11) ! (mag) F814W magnitude (3)
      real*4        e_F814W_1  (nr__11) ! (mag) rms uncertainty on F814W
      real*4        Vmag_1     (nr__11) ! (mag) V magnitude (4)
      real*4        Imag_1     (nr__11) ! (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    (nr__12) ! (deg) Right Ascension J2000
      real*8        DEdeg_4    (nr__12) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      real*4        Xpos_4     (nr__12) ! (pix) X position (1)
      real*4        Ypos_4     (nr__12) ! (pix) Y position (1)
      integer*4     RAh_4      (nr__12) ! (h) Right ascension (J2000) (2)
      integer*4     RAm_4      (nr__12) ! (min) Right ascension (J2000) (2)
      real*4        RAs_4      (nr__12) ! (s) Right ascension (J2000) (2)
      character*1   DE__4      (nr__12) ! Declination sign (J2000) (2)
      integer*4     DEd_4      (nr__12) ! (deg) Declination (J2000) (2)
      integer*4     DEm_4      (nr__12) ! (arcmin) Declination (J2000) (2)
      real*4        DEs_4      (nr__12) ! (arcsec) Declination (J2000) (2)
      real*4        F555W_2    (nr__12) ! (mag) F555W magnitude(3)
      real*4        e_F555W_2  (nr__12) ! (mag) rms uncertainty on F555W
      real*4        F814W_2    (nr__12) ! (mag) F814W magnitude (3)
      real*4        e_F814W_2  (nr__12) ! (mag) rms uncertainty on F814W
      real*4        Vmag_2     (nr__12) ! (mag) V magnitude (4)
      real*4        Imag_2     (nr__12) ! (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    (nr__13) ! (deg) Right Ascension J2000
      real*8        DEdeg_5    (nr__13) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      real*4        Xpos_5     (nr__13) ! (pix) X position (1)
      real*4        Ypos_5     (nr__13) ! (pix) Y position (1)
      integer*4     RAh_5      (nr__13) ! (h) Right ascension (J2000) (2)
      integer*4     RAm_5      (nr__13) ! (min) Right ascension (J2000) (2)
      real*4        RAs_5      (nr__13) ! (s) Right ascension (J2000) (2)
      character*1   DE__5      (nr__13) ! Declination sign (J2000) (2)
      integer*4     DEd_5      (nr__13) ! (deg) Declination (J2000) (2)
      integer*4     DEm_5      (nr__13) ! (arcmin) Declination (J2000) (2)
      real*4        DEs_5      (nr__13) ! (arcsec) Declination (J2000) (2)
      real*4        F555W_3    (nr__13) ! (mag) F555W magnitude(3)
      real*4        e_F555W_3  (nr__13) ! (mag) rms uncertainty on F555W
      real*4        F814W_3    (nr__13) ! (mag) F814W magnitude (3)
      real*4        e_F814W_3  (nr__13) ! (mag) rms uncertainty on F814W
      real*4        Vmag_3     (nr__13) ! (mag) V magnitude (4)
      real*4        Imag_3     (nr__13) ! (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(nr__14) ! Identification (CNN or VNN)
      character*76  Note       (nr__14) ! 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(nr__15) ! Identification (CNN or VNN)
      character*76  Note_1     (nr__15) ! 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_(i__),Xpos(i__),Ypos(i__),RAh(i__),RAm(i__),RAs(i__),
     +  DE_(i__),DEd(i__),DEm(i__),DEs(i__),Chip(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,1)
     +  v_FBK98_(i__),Xpos(i__),Ypos(i__),RAh(i__),RAm(i__),RAs(i__),
     +  DE_(i__),DEd(i__),DEm(i__),DEs(i__),Chip(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 '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(i__),HJD(i__),mag(i__),e_mag(i__)
c    ..............Just test output...........
        write(6,2)v_FBK98__1(i__),HJD(i__),mag(i__),e_mag(i__)
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(i__),HJD_1(i__),mag_1(i__),e_mag_1(i__)
c    ..............Just test output...........
        write(6,3)
     +  v_FBK98__2(i__),HJD_1(i__),mag_1(i__),e_mag_1(i__)
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(i__),l_Per(i__),Per(i__),l_Per2(i__),Per2(i__),
     +  F555Wph(i__),F555Wav(i__),F814Wph(i__),F814WphCor(i__),
     +  F814Wav(i__),F814WavCor(i__),Vmagph(i__),Vmagav(i__),
     +  Imagph(i__),Imagav(i__),chi2_sigma(i__)
        if(ar__3(12:13) .EQ. '') Per2(i__) = iNULL__
        if(ar__3(27:31) .EQ. '') F814Wph(i__) = rNULL__
        if(ar__3(33:37) .EQ. '') F814WphCor(i__) = rNULL__
        if(ar__3(39:43) .EQ. '') F814Wav(i__) = rNULL__
        if(ar__3(45:49) .EQ. '') F814WavCor(i__) = rNULL__
        if(ar__3(51:55) .EQ. '') Vmagph(i__) = rNULL__
        if(ar__3(57:61) .EQ. '') Vmagav(i__) = rNULL__
        if(ar__3(63:67) .EQ. '') Imagph(i__) = rNULL__
        if(ar__3(69:73) .EQ. '') Imagav(i__) = rNULL__
        if(ar__3(75:79) .EQ. '') chi2_sigma(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  v_FBK98__3(i__),l_Per(i__),Per(i__),l_Per2(i__),Per2(i__),
     +  F555Wph(i__),F555Wav(i__),F814Wph(i__),F814WphCor(i__),
     +  F814Wav(i__),F814WavCor(i__),Vmagph(i__),Vmagav(i__),
     +  Imagph(i__),Imagav(i__),chi2_sigma(i__)
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(i__),l_Per_1(i__),Per_1(i__),l_Per2_1(i__),
     +  Per2_1(i__),F555Wph_1(i__),F555Wav_1(i__),F814Wph_1(i__),
     +  F814WphCor_1(i__),F814Wav_1(i__),F814WavCor_1(i__),
     +  Vmagph_1(i__),Vmagav_1(i__),Imagph_1(i__),Imagav_1(i__),
     +  chi2_sigma_1(i__)
        if(ar__4(12:13) .EQ. '') Per2_1(i__) = iNULL__
        if(ar__4(27:31) .EQ. '') F814Wph_1(i__) = rNULL__
        if(ar__4(33:37) .EQ. '') F814WphCor_1(i__) = rNULL__
        if(ar__4(39:43) .EQ. '') F814Wav_1(i__) = rNULL__
        if(ar__4(45:49) .EQ. '') F814WavCor_1(i__) = rNULL__
        if(ar__4(51:55) .EQ. '') Vmagph_1(i__) = rNULL__
        if(ar__4(57:61) .EQ. '') Vmagav_1(i__) = rNULL__
        if(ar__4(63:67) .EQ. '') Imagph_1(i__) = rNULL__
        if(ar__4(69:73) .EQ. '') Imagav_1(i__) = rNULL__
        if(ar__4(75:79) .EQ. '') chi2_sigma_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  v_FBK98__4(i__),l_Per_1(i__),Per_1(i__),l_Per2_1(i__),
     +  Per2_1(i__),F555Wph_1(i__),F555Wav_1(i__),F814Wph_1(i__),
     +  F814WphCor_1(i__),F814Wav_1(i__),F814WavCor_1(i__),
     +  Vmagph_1(i__),Vmagav_1(i__),Imagph_1(i__),Imagav_1(i__),
     +  chi2_sigma_1(i__)
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(i__),Xpos_1(i__),Ypos_1(i__),RAh_1(i__),RAm_1(i__),
     +  RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
     +  Chip_1(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)
        if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
        if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
        if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
        if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
        if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
        if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
        if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c    ..............Just test output...........
        write(6,6)
     +  v_FBK98__5(i__),Xpos_1(i__),Ypos_1(i__),RAh_1(i__),RAm_1(i__),
     +  RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
     +  Chip_1(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 '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(i__),HJD_2(i__),mag_2(i__),e_mag_2(i__)
c    ..............Just test output...........
        write(6,7)
     +  v_FBK98__6(i__),HJD_2(i__),mag_2(i__),e_mag_2(i__)
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(i__),HJD_3(i__),mag_3(i__),e_mag_3(i__)
c    ..............Just test output...........
        write(6,8)
     +  v_FBK98__7(i__),HJD_3(i__),mag_3(i__),e_mag_3(i__)
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(i__),l_Per_2(i__),Per_2(i__),l_Per2_2(i__),
     +  Per2_2(i__),F555Wph_2(i__),F555Wav_2(i__),F814Wph_2(i__),
     +  F814WphCor_2(i__),F814Wav_2(i__),F814WavCor_2(i__),
     +  Vmagph_2(i__),Vmagav_2(i__),Imagph_2(i__),Imagav_2(i__),
     +  chi2_sigma_2(i__)
        if(ar__8(12:13) .EQ. '') Per2_2(i__) = iNULL__
        if(ar__8(27:31) .EQ. '') F814Wph_2(i__) = rNULL__
        if(ar__8(33:37) .EQ. '') F814WphCor_2(i__) = rNULL__
        if(ar__8(39:43) .EQ. '') F814Wav_2(i__) = rNULL__
        if(ar__8(45:49) .EQ. '') F814WavCor_2(i__) = rNULL__
        if(ar__8(51:55) .EQ. '') Vmagph_2(i__) = rNULL__
        if(ar__8(57:61) .EQ. '') Vmagav_2(i__) = rNULL__
        if(ar__8(63:67) .EQ. '') Imagph_2(i__) = rNULL__
        if(ar__8(69:73) .EQ. '') Imagav_2(i__) = rNULL__
        if(ar__8(75:79) .EQ. '') chi2_sigma_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  v_FBK98__8(i__),l_Per_2(i__),Per_2(i__),l_Per2_2(i__),
     +  Per2_2(i__),F555Wph_2(i__),F555Wav_2(i__),F814Wph_2(i__),
     +  F814WphCor_2(i__),F814Wav_2(i__),F814WavCor_2(i__),
     +  Vmagph_2(i__),Vmagav_2(i__),Imagph_2(i__),Imagav_2(i__),
     +  chi2_sigma_2(i__)
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(i__),l_Per_3(i__),Per_3(i__),l_Per2_3(i__),
     +  Per2_3(i__),F555Wph_3(i__),F555Wav_3(i__),F814Wph_3(i__),
     +  F814WphCor_3(i__),F814Wav_3(i__),F814WavCor_3(i__),
     +  Vmagph_3(i__),Vmagav_3(i__),Imagph_3(i__),Imagav_3(i__),
     +  chi2_sigma_3(i__)
        if(ar__9(12:13) .EQ. '') Per2_3(i__) = iNULL__
        if(ar__9(27:31) .EQ. '') F814Wph_3(i__) = rNULL__
        if(ar__9(33:37) .EQ. '') F814WphCor_3(i__) = rNULL__
        if(ar__9(39:43) .EQ. '') F814Wav_3(i__) = rNULL__
        if(ar__9(45:49) .EQ. '') F814WavCor_3(i__) = rNULL__
        if(ar__9(51:55) .EQ. '') Vmagph_3(i__) = rNULL__
        if(ar__9(57:61) .EQ. '') Vmagav_3(i__) = rNULL__
        if(ar__9(63:67) .EQ. '') Imagph_3(i__) = rNULL__
        if(ar__9(69:73) .EQ. '') Imagav_3(i__) = rNULL__
        if(ar__9(75:79) .EQ. '') chi2_sigma_3(i__) = rNULL__
c    ..............Just test output...........
        write(6,10)
     +  v_FBK98__9(i__),l_Per_3(i__),Per_3(i__),l_Per2_3(i__),
     +  Per2_3(i__),F555Wph_3(i__),F555Wav_3(i__),F814Wph_3(i__),
     +  F814WphCor_3(i__),F814Wav_3(i__),F814WavCor_3(i__),
     +  Vmagph_3(i__),Vmagav_3(i__),Imagph_3(i__),Imagav_3(i__),
     +  chi2_sigma_3(i__)
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(i__),Ypos_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),
     +  DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),F555W(i__),
     +  e_F555W(i__),F814W(i__),e_F814W(i__),Vmag(i__),Imag(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)
        if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
        if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
        if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
        if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
        if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
        if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
        if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c    ..............Just test output...........
        write(6,11)
     +  Xpos_2(i__),Ypos_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),
     +  DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),F555W(i__),
     +  e_F555W(i__),F814W(i__),e_F814W(i__),Vmag(i__),Imag(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 '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(i__),Ypos_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),
     +  DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),F555W_1(i__),
     +  e_F555W_1(i__),F814W_1(i__),e_F814W_1(i__),Vmag_1(i__),
     +  Imag_1(i__)
        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)
        if(RAh_3(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
        if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
        if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240.
        if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
        if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
        if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600.
        if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c    ..............Just test output...........
        write(6,12)
     +  Xpos_3(i__),Ypos_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),
     +  DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),F555W_1(i__),
     +  e_F555W_1(i__),F814W_1(i__),e_F814W_1(i__),Vmag_1(i__),
     +  Imag_1(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
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(i__),Ypos_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),
     +  DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),F555W_2(i__),
     +  e_F555W_2(i__),F814W_2(i__),e_F814W_2(i__),Vmag_2(i__),
     +  Imag_2(i__)
        RAdeg_4(i__) = rNULL__
        DEdeg_4(i__) = 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(i__) .GT. -180) RAdeg_4(i__)=RAh_4(i__)*15.
        if(RAm_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAm_4(i__)/4.
        if(RAs_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAs_4(i__)/240.
        if(DEd_4(i__) .GE. 0) DEdeg_4(i__)=DEd_4(i__)
        if(DEm_4(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEm_4(i__)/60.
        if(DEs_4(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEs_4(i__)/3600.
        if(DE__4(i__).EQ.'-'.AND.DEdeg_4(i__).GE.0) DEdeg_4(i__)=-DEdeg_4(i__)
c    ..............Just test output...........
        write(6,13)
     +  Xpos_4(i__),Ypos_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),
     +  DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),F555W_2(i__),
     +  e_F555W_2(i__),F814W_2(i__),e_F814W_2(i__),Vmag_2(i__),
     +  Imag_2(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_4(i__),DEdeg_4(i__)
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(i__),Ypos_5(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),
     +  DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__),F555W_3(i__),
     +  e_F555W_3(i__),F814W_3(i__),e_F814W_3(i__),Vmag_3(i__),
     +  Imag_3(i__)
        RAdeg_5(i__) = rNULL__
        DEdeg_5(i__) = 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(i__) .GT. -180) RAdeg_5(i__)=RAh_5(i__)*15.
        if(RAm_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAm_5(i__)/4.
        if(RAs_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAs_5(i__)/240.
        if(DEd_5(i__) .GE. 0) DEdeg_5(i__)=DEd_5(i__)
        if(DEm_5(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEm_5(i__)/60.
        if(DEs_5(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEs_5(i__)/3600.
        if(DE__5(i__).EQ.'-'.AND.DEdeg_5(i__).GE.0) DEdeg_5(i__)=-DEdeg_5(i__)
c    ..............Just test output...........
        write(6,14)
     +  Xpos_5(i__),Ypos_5(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),
     +  DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__),F555W_3(i__),
     +  e_F555W_3(i__),F814W_3(i__),e_F814W_3(i__),Vmag_3(i__),
     +  Imag_3(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_5(i__),DEdeg_5(i__)
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(i__),Note(i__)
c    ..............Just test output...........
        write(6,15)v_FBK98__10(i__),Note(i__)
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(i__),Note_1(i__)
c    ..............Just test output...........
        write(6,16)v_FBK98__11(i__),Note_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end