Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/788/119 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/788/119 Properties of the highly ionized gas of quasars (Lehner+, 2014) *================================================================================ *Galactic and circumgalactic O VI and its impact on the cosmological metal *and baryon budgets at 2<z<~3.5. * Lehner N., O'Meara J.M., Fox A.J., Howk J.C., Prochaska J.X., Burns V., * Armstrong A.A. * <Astrophys. J., 788, 119 (2014)> * =2014ApJ...788..119L (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Data Sample Summary integer*4 nr__ parameter (nr__=20) ! Number of records character*77 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*11 ID ! Source identifier integer*4 RAh ! (h) Hour of Right Ascension (J2000) integer*4 RAm ! (min) Minute of Right Ascension (J2000) real*4 RAs ! (s) Second of Right Ascension (J2000) character*1 DE_ ! Sign of the Declination (J2000) integer*4 DEd ! (deg) Degree of Declination (J2000) integer*4 DEm ! (arcmin) Arcminute of Declination (J2000) real*4 DEs ! (arcsec) Arcsecond of Declination (J2000) character*15 PI ! The original PI who acquired the data integer*4 Res ! (10+3) [36/72] Spectral resolution character*1 n_Res ! [d] Note on Res (1) integer*4 S_N ! [6/54] Signal-to-noise ratio (2) real*4 z ! QSO redshift real*8 zabs ! Absorber redshift (3) real*4 deltav ! (10+3km/s) Velocity difference (4) *Note (1): Note as follows: * d = Also, R=36000 in the red (C IV, as well Si IV at high S/N). *Note (2): Signal to noise (S/N) near O VI {lambda}1031; generally, the S/N is * higher at higher wavelengths. *Note (3): The absorber redshift is determined by the centroid relative to the * neutral (e.g., H I or O I) or low ionization (e.g., C II or Si II) species * and will not generally match the exact redshift of the centroid to the * observed O VI. *Note (4): deltav=c(z_em_-z_abs_)/(1+z_abs_). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Results from the Profile Fits integer*4 nr__1 parameter (nr__1=256) ! Number of records character*67 ar__1 ! Full-size record character*11 ID_1 ! Source identifier real*8 z_1 ! Redshift of the absorber character*1 f_z ! [ip] Flag on z (1) character*5 Ion ! Fitted ion integer*4 Num ! [1/16] Component number real*4 Vel ! (km/s) Velocity real*4 e_Vel ! (km/s) Uncertainty in Vel real*4 b ! (km/s) Doppler parameter real*4 e_b ! (km/s) Uncertainty in b real*4 logN ! ([cm-2]) Log column density real*4 e_logN ! ([cm-2]) Uncertainty in logN integer*4 Flag ! [0/3]? Fit flag (2) *Note (1): Flag as follows: * i = intervening; * p = proximate. *Note (2): Flag as follows: * 0 = result is robust; * 1 = uncertain owing to the complexity of the absorption profile or on the * necessity to rely on only one transition of the doublet; * 2 = uncertain owing to the relative low S/N; * 3 = uncertain results because the component is saturated in both transitions * of the fitted doublet. C============================================================================= C Loading file 'table1.dat' ! Data Sample Summary C Format for file interpretation 1 format( + A11,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A15,1X,I2, + A1,1X,I2,1X,F5.3,1X,F7.5,1X,F4.1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,20 read(1,'(A77)')ar__ read(ar__,1) + ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,PI,Res,n_Res,S_N,z,zabs,deltav 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) + ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,PI,Res,n_Res,S_N,z,zabs,deltav write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! Results from the Profile Fits C Format for file interpretation 2 format( + A11,1X,F7.5,1X,A1,1X,A5,1X,I2,1X,F6.1,1X,F4.1,1X,F5.1,1X,F4.1, + 1X,F5.2,1X,F4.2,2X,I1) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,256 read(1,'(A67)')ar__1 read(ar__1,2) + ID_1,z_1,f_z,Ion,Num,Vel,e_Vel,b,e_b,logN,e_logN,Flag if(ar__1(67:67) .EQ. '') Flag = iNULL__ c ..............Just test output........... write(6,2) + ID_1,z_1,f_z,Ion,Num,Vel,e_Vel,b,e_b,logN,e_logN,Flag c .......End.of.Just test output........... end do close(1) C============================================================================= stop end