Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/448/893 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe C============================================================================= C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Mar-29 C============================================================================= * This code was generated from the ReadMe file documenting a catalogue * according to the "Standard for Documentation of Astronomical Catalogues" * currently in use by the Astronomical Data Centers (CDS, ADC, A&A) * (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx) * Please report problems or questions to C============================================================================= implicit none * Unspecified or NULL values, generally corresponding to blank columns, * are assigned one of the following special values: * rNULL__ for unknown or NULL floating-point values * iNULL__ for unknown or NULL integer values real*4 rNULL__ integer*4 iNULL__ parameter (rNULL__=--2147483648.) ! NULL real number parameter (iNULL__=(-2147483647-1)) ! NULL int number integer idig ! testing NULL number C============================================================================= Cat. J/A+A/448/893 Photometry and lines of (z~0.2-1.0) galaxies (Lamareille+, 2006) *================================================================================ *Spectrophotometric properties of galaxies at intermediate redshifts (z~0.2-1.0). *I. Sample description, photometric properties and spectral measurements. * Lamareille F., Contini T., Le Borgne J.-F., Brinchmann J., Charlot S., * Richard J. * <Astron. Astrophys. 448, 893 (2006)> * =2006A&A...448..893L C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'objects.dat' ! List of objects integer*4 nr__ parameter (nr__=141) ! Number of records character*54 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) integer*4 LCL05 (nr__) ! Identification number character*13 OName (nr__) ! Alternative identification if available (1) integer*4 RAh (nr__) ! (h) Right Ascension J2000 (hours) integer*4 RAm (nr__) ! (min) Right Ascension J2000 (minutes) real*4 RAs (nr__) ! (s) Right Ascension J2000 (seconds) character*1 DE_ (nr__) ! Declination J2000 (sign) integer*4 DEd (nr__) ! (deg) Declination J2000 (degrees) integer*4 DEm (nr__) ! (arcmin) Declination J2000 (minutes) real*4 DEs (nr__) ! (arcsec) Declination J2000 (seconds) real*4 z (nr__) ! Redshift *Note (1): Other names use the acronyms: * CFRS = Canada-France Redshift Survey, Cat. <VII/225> * GDDS = Gemini Deep Deep Survey, Cat. <J/AJ/127/2455> * [CBB2001] = Couch et al., Cat. <J/ApJ/549/820> * [CPK2001] = Campusano et al., 2001A&A...378..394C * [LBP2003] = Le Borgne et al., 2003A&A...402..433L * [SKK2001] = Smail et al., 2001MNRAS.323..839S c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Photometric data integer*4 nr__1 parameter (nr__1=141) ! Number of records character*94 ar__1 ! Full-size record integer*4 LCL05_1 (nr__1) ! Identification number real*4 Rmag (nr__1) ! (mag) Observed photometric magnitude integer*4 Filter (nr__1) ! [1-6] Identification of photometric filter (1) real*4 Daper (nr__1) ! (mag) Aperture difference * (photometric-spectroscopic magnitude) real*4 I0mag (nr__1) ! (mag) k-corrected magnitude in I band (2) real*4 AI (nr__1) ! (mag) Foreground dust attenuation in the rest-frame * I-band (2) real*4 v_U_B_0 (nr__1) ! (mag) ?=- U-B k-corrected spectroscopic color index (2) real*4 v_B_V_0 (nr__1) ! (mag) ?=- B-V k-corrected spectroscopic color index (2) real*4 v_V_R_0 (nr__1) ! (mag) V-R k-corrected spectroscopic color index (2) real*4 v_V_I_0 (nr__1) ! (mag) V-I k-corrected spectroscopic color index (2) real*4 v_u_g_0 (nr__1) ! (mag) ?=- u-g k-corrected spectroscopic color index (2) real*4 v_g_r_0 (nr__1) ! (mag) g-r k-corrected spectroscopic color index (2) real*4 v_g_i_0 (nr__1) ! (mag) g-i k-corrected spectroscopic color index (2) real*4 IMAG (nr__1) ! (mag) k-corrected absolute magnitude in I band (2) real*4 e_IMAG (nr__1) ! (mag) Error in the magnitude (2) real*4 Lens (nr__1) ! (mag) Lensing correction when necessary (3) *Note (1): Photometric filter as follows: * 1 = FORS R Gunn * 2 = FORS R Bessel * 3 = TEK2048 R * 4 = BTC R * 5 = CFH12k I * 6 = WFPC F606W *Note (2): All magnitudes are in AB system. * All colors and absolute magnitude are corrected from foreground dust. *Note (3): The lensing correction is already applied to the absolute magnitude. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Spectral classification integer*4 nr__2 parameter (nr__2=141) ! Number of records character*25 ar__2 ! Full-size record integer*4 LCL05_2 (nr__2) ! Identification number integer*4 red (nr__2) ! [0-3] "red" classification (1) integer*4 blue (nr__2) ! [0-4] "blue" classification (1) integer*4 Class (nr__2) ! [0-4] final classification (1) real*4 E_B_V (nr__2) ! (mag) ?=- Extinction coefficient (2) real*4 e_E_B_V (nr__2) ! (mag) ?=- Error on extinction coefficient *Note (1): Classifications as follows: * 0 = not applicable * 1 = star-forming galaxy * 2 = Seyfert 2 galaxy * 3 = Mixed object * 4 = candidate star-forming galaxy *Note (2): Extinction coefficient computed from Halpha/Hbeta emission-line ratio. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Forbidden emission lines measurements integer*4 nr__3 parameter (nr__3=141) ! Number of records character*124 ar__3 ! Full-size record integer*4 LCL05_3 (nr__3) ! Identification number real*4 FO2 (nr__3) ! ?=- [OII]3727 emission line flux (G1) real*4 e_FO2 (nr__3) ! ?=- Error on FO2 (G1) real*4 EWO2 (nr__3) ! (0.1nm) ?=- [OII]3727 emission line * equivalent width (G2) real*4 e_EWO2 (nr__3) ! (0.1nm) ?=- Error on EWO2 real*4 FO3a (nr__3) ! ?=- [OIII]4959 emission line flux (G1) real*4 e_FO3a (nr__3) ! ?=- Error on FO3a (G1) real*4 EWO3a (nr__3) ! (0.1nm) ?=- [OIII]4959 emission line * equivalent width (G2) real*4 e_EWO3a (nr__3) ! (0.1nm) ?=- Error on EWO3a real*4 FO3b (nr__3) ! ?=- [OIII]5007 emission line flux (G1) real*4 e_FO3b (nr__3) ! ?=- Error on FO3b (G1) real*4 EWO3b (nr__3) ! (0.1nm) ?=- [OIII]5007 emission line * equivalent width (G2) real*4 e_EWO3b (nr__3) ! (0.1nm) ?=- Error on EWO3b real*4 FN2 (nr__3) ! ?=- [NII]6584 emission line flux (G1) real*4 e_FN2 (nr__3) ! ?=- Error on FN2 (G1) real*4 EWN2 (nr__3) ! (0.1nm) ?=- [NII]6584 emission line * equivalent width (G2) real*4 e_EWN2 (nr__3) ! (0.1nm) ?=- Error on EWN2 real*4 FS2 (nr__3) ! ?=- [SII]6717+6731 emission line flux (G1) real*4 e_FS2 (nr__3) ! ?=- Error on FS2 (G1) real*4 EWS2 (nr__3) ! (0.1nm) ?=- [SII]6717+6731 emission line * equivalent width (G2) real*4 e_EWS2 (nr__3) ! (0.1nm) ?=- Error on EWS2 c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table9.dat' ! Balmer lines measurements integer*4 nr__4 parameter (nr__4=141) ! Number of records character*125 ar__4 ! Full-size record integer*4 LCL05_4 (nr__4) ! Identification number real*4 FHa (nr__4) ! ?=- Halpha flux of emission line (G1) real*4 e_FHa (nr__4) ! ?=- Error on FHa (G1) real*4 EWHa (nr__4) ! (0.1nm) ?=- Halpha equivalent width of * emission line (G2) real*4 e_EWHa (nr__4) ! (0.1nm) ?=- Error on EWHa (G2) real*4 EWHaabs (nr__4) ! (0.1nm) ?=- Halpha equivalent width of * absorption line (G2) real*4 FHb (nr__4) ! (10-20W/m2) Hbeta flux of emission line real*4 e_FHb (nr__4) ! (10-20W/m2) Error on FHb real*4 EWHb (nr__4) ! (0.1nm) ?=- Hbeta equivalent width of * emission line (G2) real*4 e_EWHb (nr__4) ! (0.1nm) ?=- Error on EWHb (G2) real*4 EWHbabs (nr__4) ! (0.1nm) Hbeta equivalent width of * absorption line (G2) real*4 FHg (nr__4) ! ?=- Hgamma flux of emission line (G1) real*4 e_FHg (nr__4) ! ?=- Error on FHg (G1) real*4 EWHg (nr__4) ! (0.1nm) ?=- Hgamma equivalent width of * emission line (G2) real*4 e_EWHg (nr__4) ! (0.1nm) ?=- Error on EWHg (G2) real*4 EWHgabs (nr__4) ! (0.1nm) ?=- Hgamma equivalent width of * absorption line (G2) real*4 FHd (nr__4) ! ?=- Hdelta flux of emission line (G1) real*4 e_FHd (nr__4) ! ?=- Error on FHd (G1) real*4 EWHd (nr__4) ! (0.1nm) ?=- Hdelta equivalent width of * emission line (G2) real*4 e_EWHd (nr__4) ! (0.1nm) ?=- Error on EWHd (G2) real*4 EWHdabs (nr__4) ! (0.1nm) ?=- Hdelta equivalent width of * absorption line (G2) C============================================================================= C Loading file 'objects.dat' ! List of objects C Format for file interpretation 1 format( + 6X,I3,1X,A13,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X, + F6.4) C Effective file loading open(unit=1,status='old',file= +'objects.dat') write(6,*) '....Loading file: objects.dat' do i__=1,141 read(1,'(A54)')ar__ read(ar__,1) + LCL05(i__),OName(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__), + DEd(i__),DEm(i__),DEs(i__),z(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) + LCL05(i__),OName(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__), + DEd(i__),DEm(i__),DEs(i__),z(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' ! Photometric data C Format for file interpretation 2 format( + 6X,I3,1X,F5.2,1X,I1,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,F6.2,1X,F4.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,141 read(1,'(A94)')ar__1 read(ar__1,2) + LCL05_1(i__),Rmag(i__),Filter(i__),Daper(i__),I0mag(i__), + AI(i__),v_U_B_0(i__),v_B_V_0(i__),v_V_R_0(i__),v_V_I_0(i__), + v_u_g_0(i__),v_g_r_0(i__),v_g_i_0(i__),IMAG(i__),e_IMAG(i__), + Lens(i__) if (idig(ar__1(37:41)).EQ.0) v_U_B_0(i__) = rNULL__ if (idig(ar__1(43:47)).EQ.0) v_B_V_0(i__) = rNULL__ if (idig(ar__1(61:65)).EQ.0) v_u_g_0(i__) = rNULL__ c ..............Just test output........... write(6,2) + LCL05_1(i__),Rmag(i__),Filter(i__),Daper(i__),I0mag(i__), + AI(i__),v_U_B_0(i__),v_B_V_0(i__),v_V_R_0(i__),v_V_I_0(i__), + v_u_g_0(i__),v_g_r_0(i__),v_g_i_0(i__),IMAG(i__),e_IMAG(i__), + Lens(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! Spectral classification C Format for file interpretation 3 format(6X,I3,1X,I1,1X,I1,1X,I1,1X,F4.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,141 read(1,'(A25)')ar__2 read(ar__2,3) + LCL05_2(i__),red(i__),blue(i__),Class(i__),E_B_V(i__), + e_E_B_V(i__) if (idig(ar__2(17:20)).EQ.0) E_B_V(i__) = rNULL__ if (idig(ar__2(22:25)).EQ.0) e_E_B_V(i__) = rNULL__ c ..............Just test output........... write(6,3) + LCL05_2(i__),red(i__),blue(i__),Class(i__),E_B_V(i__), + e_E_B_V(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! Forbidden emission lines measurements C Format for file interpretation 4 format( + 6X,I3,1X,F6.2,1X,F3.1,1X,F6.1,1X,F4.1,1X,F6.2,1X,F3.1,1X,F6.1, + 1X,F4.1,1X,F6.2,1X,F3.1,1X,F6.1,1X,F4.1,1X,F6.2,1X,F3.1,1X, + F6.1,1X,F4.1,1X,F6.2,1X,F3.1,1X,F6.1,1X,F4.1) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,141 read(1,'(A124)')ar__3 read(ar__3,4) + LCL05_3(i__),FO2(i__),e_FO2(i__),EWO2(i__),e_EWO2(i__), + FO3a(i__),e_FO3a(i__),EWO3a(i__),e_EWO3a(i__),FO3b(i__), + e_FO3b(i__),EWO3b(i__),e_EWO3b(i__),FN2(i__),e_FN2(i__), + EWN2(i__),e_EWN2(i__),FS2(i__),e_FS2(i__),EWS2(i__), + e_EWS2(i__) if (idig(ar__3(11:16)).EQ.0) FO2(i__) = rNULL__ if (idig(ar__3(18:20)).EQ.0) e_FO2(i__) = rNULL__ if (idig(ar__3(22:27)).EQ.0) EWO2(i__) = rNULL__ if (idig(ar__3(29:32)).EQ.0) e_EWO2(i__) = rNULL__ if (idig(ar__3(34:39)).EQ.0) FO3a(i__) = rNULL__ if (idig(ar__3(41:43)).EQ.0) e_FO3a(i__) = rNULL__ if (idig(ar__3(45:50)).EQ.0) EWO3a(i__) = rNULL__ if (idig(ar__3(52:55)).EQ.0) e_EWO3a(i__) = rNULL__ if (idig(ar__3(57:62)).EQ.0) FO3b(i__) = rNULL__ if (idig(ar__3(64:66)).EQ.0) e_FO3b(i__) = rNULL__ if (idig(ar__3(68:73)).EQ.0) EWO3b(i__) = rNULL__ if (idig(ar__3(75:78)).EQ.0) e_EWO3b(i__) = rNULL__ if (idig(ar__3(80:85)).EQ.0) FN2(i__) = rNULL__ if (idig(ar__3(87:89)).EQ.0) e_FN2(i__) = rNULL__ if (idig(ar__3(91:96)).EQ.0) EWN2(i__) = rNULL__ if (idig(ar__3(98:101)).EQ.0) e_EWN2(i__) = rNULL__ if (idig(ar__3(103:108)).EQ.0) FS2(i__) = rNULL__ if (idig(ar__3(110:112)).EQ.0) e_FS2(i__) = rNULL__ if (idig(ar__3(114:119)).EQ.0) EWS2(i__) = rNULL__ if (idig(ar__3(121:124)).EQ.0) e_EWS2(i__) = rNULL__ c ..............Just test output........... write(6,4) + LCL05_3(i__),FO2(i__),e_FO2(i__),EWO2(i__),e_EWO2(i__), + FO3a(i__),e_FO3a(i__),EWO3a(i__),e_EWO3a(i__),FO3b(i__), + e_FO3b(i__),EWO3b(i__),e_EWO3b(i__),FN2(i__),e_FN2(i__), + EWN2(i__),e_EWN2(i__),FS2(i__),e_FS2(i__),EWS2(i__), + e_EWS2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table9.dat' ! Balmer lines measurements C Format for file interpretation 5 format( + 6X,I3,1X,F6.2,1X,F3.1,1X,F6.1,1X,F4.1,1X,F5.1,1X,F6.2,1X,F3.1, + 1X,F6.1,1X,F4.1,1X,F5.1,1X,F6.2,1X,F3.1,1X,F6.1,1X,F4.1,1X, + F5.1,1X,F6.2,1X,F3.1,1X,F6.1,1X,F4.1,1X,F5.1) C Effective file loading open(unit=1,status='old',file= +'table9.dat') write(6,*) '....Loading file: table9.dat' do i__=1,141 read(1,'(A125)')ar__4 read(ar__4,5) + LCL05_4(i__),FHa(i__),e_FHa(i__),EWHa(i__),e_EWHa(i__), + EWHaabs(i__),FHb(i__),e_FHb(i__),EWHb(i__),e_EWHb(i__), + EWHbabs(i__),FHg(i__),e_FHg(i__),EWHg(i__),e_EWHg(i__), + EWHgabs(i__),FHd(i__),e_FHd(i__),EWHd(i__),e_EWHd(i__), + EWHdabs(i__) if (idig(ar__4(11:16)).EQ.0) FHa(i__) = rNULL__ if (idig(ar__4(18:20)).EQ.0) e_FHa(i__) = rNULL__ if (idig(ar__4(22:27)).EQ.0) EWHa(i__) = rNULL__ if (idig(ar__4(29:32)).EQ.0) e_EWHa(i__) = rNULL__ if (idig(ar__4(34:38)).EQ.0) EWHaabs(i__) = rNULL__ if (idig(ar__4(51:56)).EQ.0) EWHb(i__) = rNULL__ if (idig(ar__4(58:61)).EQ.0) e_EWHb(i__) = rNULL__ if (idig(ar__4(69:74)).EQ.0) FHg(i__) = rNULL__ if (idig(ar__4(76:78)).EQ.0) e_FHg(i__) = rNULL__ if (idig(ar__4(80:85)).EQ.0) EWHg(i__) = rNULL__ if (idig(ar__4(87:90)).EQ.0) e_EWHg(i__) = rNULL__ if (idig(ar__4(92:96)).EQ.0) EWHgabs(i__) = rNULL__ if (idig(ar__4(98:103)).EQ.0) FHd(i__) = rNULL__ if (idig(ar__4(105:107)).EQ.0) e_FHd(i__) = rNULL__ if (idig(ar__4(109:114)).EQ.0) EWHd(i__) = rNULL__ if (idig(ar__4(116:119)).EQ.0) e_EWHd(i__) = rNULL__ if (idig(ar__4(121:125)).EQ.0) EWHdabs(i__) = rNULL__ c ..............Just test output........... write(6,5) + LCL05_4(i__),FHa(i__),e_FHa(i__),EWHa(i__),e_EWHa(i__), + EWHaabs(i__),FHb(i__),e_FHb(i__),EWHb(i__),e_EWHb(i__), + EWHbabs(i__),FHg(i__),e_FHg(i__),EWHg(i__),e_EWHg(i__), + EWHgabs(i__),FHd(i__),e_FHd(i__),EWHd(i__),e_EWHd(i__), + EWHdabs(i__) 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