Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/466/481 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.
++++Only ONE component of position, remove MPOS flag: RAh (the full equatorial position requires RAh RAm RAs DE- DEd DEm DEs) for: table1.dat program load_ReadMe C============================================================================= C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Mar-28 C============================================================================= * This code was generated from the ReadMe file documenting a catalogue * according to the "Standard for Documentation of Astronomical Catalogues" * currently in use by the Astronomical Data Centers (CDS, ADC, A&A) * (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx) * Please report problems or questions to C============================================================================= implicit none * Unspecified or NULL values, generally corresponding to blank columns, * are assigned one of the following special values: * rNULL__ for unknown or NULL floating-point values * iNULL__ for unknown or NULL integer values real*4 rNULL__ integer*4 iNULL__ parameter (rNULL__=--2147483648.) ! NULL real number parameter (iNULL__=(-2147483647-1)) ! NULL int number integer idig ! testing NULL number C============================================================================= Cat. J/A+A/466/481 Near-infrared survey around PKS 1343-601 (Schroeder+, 2007) *================================================================================ *The highly obscured region around PKS 1343-601. *I. Galactic interstellar extinctions using DENIS galaxy colours. * Schroeder A.C., Mamon G.A., Kraan-Korteweg R.C., Woudt P.A. * <Astron. Astrophys., 466, 481-499 (2007)> * =2007A&A...466..481S C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! DENIS strips in the searched area integer*4 nr__ parameter (nr__=43) ! Number of records character*76 ar__ ! Full-size record integer*4 Slot ! Slot number (according to their designation * in the sky) integer*4 Strip ! Strip number (obs. of a given slot with a * unique number according to observing date) (1) character*10 Date ! ("DD/MM/YYYY") Date of observation integer*4 RAh ! (h) Strip central Right Ascension (J2000) (2) integer*4 RAm ! (min) Strip central Right Ascension (J2000) (2) real*4 RAs ! (s) Strip central Right Ascension (J2000) (2) integer*4 Imagef ! Image numbers range for the searched area integer*4 Imagel ! Image numbers range for the searched area character*3 Band ! [IJK -] Observed passbands real*4 PSFI ! ?=- Point-spread function on I image (3) real*4 PSFJ ! ?=- Point-spread function on J image (3) real*4 PSFK ! ?=- Point-spread function on K image (3) integer*4 Qual ! ?=- Seeing quality as estimated from the * photometry of stars (4) character*1 n_Qual ! [)] Note on Qual (5) integer*4 Weath ! ?=- Weather conditions (6) integer*4 q_Strip ! ?=- astrometric quality of the strips (7) integer*4 Prior ! ?=- Order of quality of strips with more than * one observation (8) *Note (1): Strip numbers smaller than 3000 were part of the pre-survey * (these have no I-band counterparts, and the exposure time was * slightly larger); they were used for verification only since they * are not astrometrically and photometrically calibrated. *Note (2): Recall that each image has a width of about 12'. *Note (3): Half-flux radius as calculated by SExtractor multiplied by * 2 (an equivalent to the point-spread function, PSF). The half-flux * radii were determined from all stars in the magnitude ranges * 11.0m<I<15.0m, 9.0m<J<13.0m, and 7.0m<K<10.5m that do not lie * within 50 pixels of the image border. *Note (4): Quality flag as follows: * 1 = good seeing (photometry not affected) * 2 = medium seeing (7"-aperture photometry is affected) * 3 = bad seeing (all magnitudes affected) *Note (5): An ) indicates that the Qual value is an estimate of the * seeing where the astrometry was insufficient for a statistical * comparison of the photometry in the overlap regions. *Note (6): Weather conditions as follows: * 1 = no clouds * 2 = unknown conditions * 3 = possible clouds * 4 = clouds *Note (7): The DENIS standard high accuracy of an rms of 0.2" is * indicated with a 1, lower accuracy with 2. For strips where * astrometric calibration failed we determined the coordinates * from the DSS*-red images (accuracy of about 2-3"). *Note (8): Order of quality of strips with more than one observation * (1 indicates the best strip, higher numbers designate lower priorities). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Galaxies in the search area integer*4 nr__1 parameter (nr__1=126) ! Number of records character*145 ar__1 ! 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*7 DZOA ! DZOA designation (NNNN-NN) (1) integer*4 Nobs ! Total number of sightings of same galaxy (2) character*1 n_Nobs ! [+] + when additional observation from the * pre-survey integer*4 RAh_1 ! (h) Right ascension (J2000) integer*4 RAm_1 ! (min) Right ascension (J2000) real*4 RAs_1 ! (s) Right ascension (J2000) character*1 u_RAs ! Uncertainty flag on Right ascension character*1 DE_ ! Declination sign (J2000) integer*4 DEd ! (deg) Declination (J2000) integer*4 DEm ! (arcmin) Declination (J2000) integer*4 DEs ! (arcsec) Declination (J2000) character*1 u_DEs ! Uncertainty flag on Declination real*4 GLON ! (deg) Galactic longitude real*4 GLAT ! (deg) Galactic latitude real*4 AB ! (mag) B-band extinction character*2 Class ! Classification of candidate (3) integer*4 Bvis ! [0/1] Visibility in B band (G1) integer*4 Ivis ! [0/1] Visibility in I band (G1) integer*4 Jvis ! [0/1] Visibility in J band (G1) integer*4 Kvis ! [0/1] Visibility in K band (G1) character*4 MType ! Morphological type real*4 Imag ! (mag) ?=- Total I magnitude real*4 e_Imag ! (mag) ?=- rms uncertainty on Imag real*4 Jmag ! (mag) ?=- Total J magnitude real*4 e_Jmag ! (mag) ?=- rms uncertainty on Jmag real*4 Kmag ! (mag) ?=- Total K magnitude real*4 e_Kmag ! (mag) ?=- rms uncertainty on Kmag real*4 v_I_J_0 ! (mag) ?=- Dereddened I-J colour index real*4 e__I_J_0 ! (mag) ?=- rms uncertainty on (I-J)0 real*4 v_I_K_0 ! (mag) ?=- Dereddened I-K colour index real*4 e__I_K_0 ! (mag) ?=- rms uncertainty on (I-K)0 real*4 v_J_K_0 ! (mag) ?=- Dereddened J-K colour index real*4 e__J_K_0 ! (mag) ?=- rms uncertainty on (J-K)0 integer*4 q_Imag ! Quality of I-band photometry (4) integer*4 q_Jmag ! Quality of I-band photometry (4) integer*4 q_Kmag ! Quality of I-band photometry (4) integer*4 Idiam ! (arcsec) ?=- Major diameter in I band integer*4 Jdiam ! (arcsec) ?=- Major diameter in J band integer*4 Kdiam ! (arcsec) ?=- Major diameter in K band *Note (1): Identity consisting of the lettering DZOA (for DENIS-ZOA) * followed by the DENIS slot number and a consecutive number for * each initial candidate *Note (2): Total number of sightings of same galaxy * (e.g., on overlap regions or repeat observations) *Note (3): Classification as follows: * DG = DENIS galaxy * UG = uncertain galaxy * BG = galaxy that was identified after consulting the WKK catalogue * NG = B-band galaxy that was found to be non-galaxian with DENIS *Note (4): Quality as follows: * 0 = good photometry * 2 = uncertainty in the 7"-aperture magnitude * 3 = uncertainty in the total magnitude * 4 = uncertainty in both * 5 = means the 7"-aperture magnitude is unreliable * 6 = unreliability in the total magnitude * 7 = both 7"-aperture and total magnitude are unreliable * 9 = the strip was non-photometric (strip 6052) and all photometry * is unreliable (the estimated extinction due to the clouds in * this case is ~0.25-0.3mag. Magnitudes with quality parameters * 5-9 have been excluded in the further analysis. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Galactic objects detected in the search area integer*4 nr__2 parameter (nr__2=27) ! Number of records character*108 ar__2 ! 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*7 DZOA_1 ! DZOA designation (NNNN-NN) character*1 n_DZOA ! [a-d] Individual note (1) integer*4 Nobs_1 ! Total number of sightings of same galaxy character*2 n_Nobs_1 ! [+ ] + when additional observation from the * pre-survey integer*4 RAh_2 ! (h) Right ascension (J2000) integer*4 RAm_2 ! (min) Right ascension (J2000) real*4 RAs_2 ! (s) Right ascension (J2000) character*1 u_RAs_1 ! Uncertainty flag on Right ascension character*1 DE__1 ! Declination sign (J2000) integer*4 DEd_1 ! (deg) Declination (J2000) integer*4 DEm_1 ! (arcmin) Declination (J2000) integer*4 DEs_1 ! (arcsec) Declination (J2000) character*1 u_DEs_1 ! Uncertainty flag on Declination real*4 GLON_1 ! (deg) Galactic longitude real*4 GLAT_1 ! (deg) Galactic latitude real*4 AB_1 ! (mag) B-band extinction character*4 Class_1 ! Classification integer*4 Bvis_1 ! [0/1] Visibility in B band (G1) integer*4 Ivis_1 ! [0/1] Visibility in I band (G1) integer*4 Jvis_1 ! [0/1] Visibility in J band (G1) integer*4 Kvis_1 ! [0/1] Visibility in K band (G1) character*39 ONames ! Other names *Note (1): Individual notes as follows: * a = While this object shows two apparent point sources in the K-band, * the B- and I-bands clearly show an elongated compact extended * emission with the brightest part positioned on one of the point * sources. * b = This object shows diffuse emission, as well as clumps, and * possibly a dust lane. It could be a PN or a YSO seen on-edge. * c = Very faint, it seems to show comparable features to DZOA4652-10 * and DZOA4656-08. * d = An extremely bright point source in the K-band, is also shows * diffuse and clumpy extended emission and possible a dust lane * in all bands. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Cross-identifications with other catalogues integer*4 nr__3 parameter (nr__3=126) ! Number of records character*117 ar__3 ! 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) character*7 DZOA_2 ! DZOA designation (NNNN-NN) integer*4 RAh_3 ! (h) Right ascension (J2000) integer*4 RAm_3 ! (min) Right ascension (J2000) real*4 RAs_3 ! (s) Right ascension (J2000) character*1 u_RAs_2 ! Uncertainty flag on Right ascension character*1 DE__2 ! Declination sign (J2000) integer*4 DEd_2 ! (deg) Declination (J2000) integer*4 DEm_2 ! (arcmin) Declination (J2000) integer*4 DEs_2 ! (arcsec) Declination (J2000) character*1 u_DEs_2 ! Uncertainty flag on Declination real*4 AB_2 ! (mag) B-band extinction character*2 Class_2 ! Classification integer*4 Bvis_2 ! [0/1] Visibility in B band (G1) integer*4 Ivis_2 ! [0/1] Visibility in I band (G1) integer*4 Jvis_2 ! [0/1] Visibility in J band (G1) integer*4 Kvis_2 ! [0/1] Visibility in K band (G1) character*17 v2MASX ! 2MASX designation (JHHMMSSss-DDMMSSs) integer*4 WKK ! ? WKK number (1) character*42 ONames_1 ! Other names *Note (1): WKK number (WKK98 NNN in Simbad) from * -- Fairall et al., 1998, Cat. <J/A+AS/127/463>, * -- Woudt et al., 1998A&A...338....8W; 1999, Cat. <J/A+A/352/39> * and 2001, Cat. <J/A+A/380/441>. C============================================================================= C Loading file 'table1.dat' ! DENIS strips in the searched area C Format for file interpretation 1 format( + I4,1X,I5,1X,A10,1X,I2,1X,I2,1X,F4.1,1X,I7,1X,I7,1X,A3,1X,F4.2, + 1X,F4.2,1X,F4.2,1X,I1,A1,1X,I1,1X,I1,1X,I1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,43 read(1,'(A76)')ar__ read(ar__,1) + Slot,Strip,Date,RAh,RAm,RAs,Imagef,Imagel,Band,PSFI,PSFJ,PSFK, + Qual,n_Qual,Weath,q_Strip,Prior if (idig(ar__(54:57)).EQ.0) PSFI = rNULL__ if (idig(ar__(59:62)).EQ.0) PSFJ = rNULL__ if (idig(ar__(64:67)).EQ.0) PSFK = rNULL__ if (Qual .EQ. 45) Qual = iNULL__ if (Weath .EQ. 45) Weath = iNULL__ if (q_Strip .EQ. 45) q_Strip = iNULL__ if (Prior .EQ. 45) Prior = iNULL__ c ..............Just test output........... write(6,1) + Slot,Strip,Date,RAh,RAm,RAs,Imagef,Imagel,Band,PSFI,PSFJ,PSFK, + Qual,n_Qual,Weath,q_Strip,Prior c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Galaxies in the search area C Format for file interpretation 2 format( + 4X,A7,1X,I2,A1,1X,I2,1X,I2,1X,F4.1,A1,1X,A1,I2,1X,I2,1X,I2,A1, + 1X,F6.2,1X,F5.2,1X,F4.1,1X,A2,1X,I1,I1,I1,I1,1X,A4,1X,F5.2,1X, + F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2, + 1X,F4.2,1X,F4.2,1X,F4.2,1X,I1,I1,I1,1X,I2,1X,I2,1X,I2) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,126 read(1,'(A145)')ar__1 read(ar__1,2) + DZOA,Nobs,n_Nobs,RAh_1,RAm_1,RAs_1,u_RAs,DE_,DEd,DEm,DEs, + u_DEs,GLON,GLAT,AB,Class,Bvis,Ivis,Jvis,Kvis,MType,Imag, + e_Imag,Jmag,e_Jmag,Kmag,e_Kmag,v_I_J_0,e__I_J_0,v_I_K_0, + e__I_K_0,v_J_K_0,e__J_K_0,q_Imag,q_Jmag,q_Kmag,Idiam,Jdiam, + Kdiam if (idig(ar__1(71:75)).EQ.0) Imag = rNULL__ if (idig(ar__1(77:80)).EQ.0) e_Imag = rNULL__ if (idig(ar__1(82:86)).EQ.0) Jmag = rNULL__ if (idig(ar__1(88:91)).EQ.0) e_Jmag = rNULL__ if (idig(ar__1(93:97)).EQ.0) Kmag = rNULL__ if (idig(ar__1(99:102)).EQ.0) e_Kmag = rNULL__ if (idig(ar__1(104:107)).EQ.0) v_I_J_0 = rNULL__ if (idig(ar__1(109:112)).EQ.0) e__I_J_0 = rNULL__ if (idig(ar__1(114:117)).EQ.0) v_I_K_0 = rNULL__ if (idig(ar__1(119:122)).EQ.0) e__I_K_0 = rNULL__ if (idig(ar__1(124:127)).EQ.0) v_J_K_0 = rNULL__ if (idig(ar__1(129:132)).EQ.0) e__J_K_0 = rNULL__ if (Idiam .EQ. 45) Idiam = iNULL__ if (Jdiam .EQ. 45) Jdiam = iNULL__ if (Kdiam .EQ. 45) Kdiam = iNULL__ 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_1 .GT. -180) RAdeg=RAh_1*15. if(RAm_1 .GT. -180) RAdeg=RAdeg+RAm_1/4. if(RAs_1 .GT. -180) RAdeg=RAdeg+RAs_1/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,2) + DZOA,Nobs,n_Nobs,RAh_1,RAm_1,RAs_1,u_RAs,DE_,DEd,DEm,DEs, + u_DEs,GLON,GLAT,AB,Class,Bvis,Ivis,Jvis,Kvis,MType,Imag, + e_Imag,Jmag,e_Jmag,Kmag,e_Kmag,v_I_J_0,e__I_J_0,v_I_K_0, + e__I_K_0,v_J_K_0,e__J_K_0,q_Imag,q_Jmag,q_Kmag,Idiam,Jdiam, + Kdiam write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Galactic objects detected in the search area C Format for file interpretation 3 format( + 4X,A7,A1,1X,I1,A2,1X,I2,1X,I2,1X,F4.1,A1,1X,A1,I2,1X,I2,1X,I2, + A1,1X,F6.2,1X,F5.2,1X,F5.1,1X,A4,1X,I1,I1,I1,I1,1X,A39) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,27 read(1,'(A108)')ar__2 read(ar__2,3) + DZOA_1,n_DZOA,Nobs_1,n_Nobs_1,RAh_2,RAm_2,RAs_2,u_RAs_1,DE__1, + DEd_1,DEm_1,DEs_1,u_DEs_1,GLON_1,GLAT_1,AB_1,Class_1,Bvis_1, + Ivis_1,Jvis_1,Kvis_1,ONames 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_2 .GT. -180) RAdeg_1=RAh_2*15. if(RAm_2 .GT. -180) RAdeg_1=RAdeg_1+RAm_2/4. if(RAs_2 .GT. -180) RAdeg_1=RAdeg_1+RAs_2/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,3) + DZOA_1,n_DZOA,Nobs_1,n_Nobs_1,RAh_2,RAm_2,RAs_2,u_RAs_1,DE__1, + DEd_1,DEm_1,DEs_1,u_DEs_1,GLON_1,GLAT_1,AB_1,Class_1,Bvis_1, + Ivis_1,Jvis_1,Kvis_1,ONames write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Cross-identifications with other catalogues C Format for file interpretation 4 format( + 4X,A7,1X,I2,1X,I2,1X,F4.1,A1,1X,A1,I2,1X,I2,1X,I2,A1,1X,F4.1, + 1X,A2,1X,I1,I1,I1,I1,1X,A17,5X,I4,1X,A42) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,126 read(1,'(A117)')ar__3 read(ar__3,4) + DZOA_2,RAh_3,RAm_3,RAs_3,u_RAs_2,DE__2,DEd_2,DEm_2,DEs_2, + u_DEs_2,AB_2,Class_2,Bvis_2,Ivis_2,Jvis_2,Kvis_2,v2MASX,WKK, + ONames_1 if(ar__3(71:74) .EQ. '') WKK = iNULL__ 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_3 .GT. -180) RAdeg_2=RAh_3*15. if(RAm_3 .GT. -180) RAdeg_2=RAdeg_2+RAm_3/4. if(RAs_3 .GT. -180) RAdeg_2=RAdeg_2+RAs_3/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,4) + DZOA_2,RAh_3,RAm_3,RAs_3,u_RAs_2,DE__2,DEd_2,DEm_2,DEs_2, + u_DEs_2,AB_2,Class_2,Bvis_2,Ivis_2,Jvis_2,Kvis_2,v2MASX,WKK, + ONames_1 write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_2 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end C Locate position of first digit in string; or return 0 integer function idig(c) character*(*) c character*1 c1 integer lc,i lc=len(c) idig=0 do i=1,lc if(c(i:i).ne.' ') go to 1 end do 1 if(i.gt.lc) return c1=c(i:i) if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1 if(i.gt.lc) return c1=c(i:i) if(c1.ge.'0'.and.c1.le.'9') idig=i return end